SystComponentScale.h

Go to the documentation of this file.

#pragma once

#include "sbnana/CAFAna/Core/Cut.h"
#include "sbnana/CAFAna/Core/ISyst.h"

namespace ana
{
  /// Uncertainty in the scale of a single component of the spectrum
  class SystComponentScale: public ISyst
  {
  public:
    enum ESystType{kLinear, kExponential};

    /// \param cut Select the component to be scaled
    SystComponentScale(const std::string& shortName,
                       const std::string& latexName,
                       const Cut& cut,
                       double oneSigma,
                       ESystType type = kExponential)
      : ISyst(shortName, latexName),
        fCut(cut), fOneSigma(oneSigma), fType(type)
    {
    }

    ~SystComponentScale();

    /// Scaling this component between 1/(1+x) and (1+x) is the 1-sigma range
    double OneSigmaScale() const {return fOneSigma;}

    const Cut& GetCut() const {return fCut;}

    void Shift(double sigma,
               caf::SRSliceProxy* slc,
               double& weight) const override;

    // Some derived classes might have a back-channel allowing them to
    // implement this.
    static std::unique_ptr<SystComponentScale> LoadFrom(TDirectory* dir);
    virtual void SaveTo(TDirectory* dir) const {assert(0 && "unimplemented");}

  protected:
    Cut fCut;
    double fOneSigma;
    ESystType fType;
  };
} // namespace