EnergySysts2017.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"

namespace ana
{
  //----------------------------------------------------------------------
  /// Absolute fit muon energy scale systematic
  class MuEScaleSyst2017: public ISyst
  {
  public:
    MuEScaleSyst2017(const std::string& shortName,
                     const std::string& latexName,
                     double scale, double scalemc)
      : ISyst(shortName, latexName), fScale(scale), fScaleMC(scalemc)
    {
    }

    MuEScaleSyst2017(double scale, double scalemc)
      : ISyst("AbsMuEScale2017", "Absolute Muon Energy Scale 2017"),
        fScale(scale), fScaleMC(scalemc)
    {
    }

    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;

  private:
    const double fScale;
    const double fScaleMC;
  };

  extern const MuEScaleSyst2017 kMuEScaleSyst2017;

  //----------------------------------------------------------------------
  // Relative fit muon energy scale systematic.  While this is named 2017, I
  // have modified it from the 2017 version to use a shift in cm for the muon
  // catcher rather than a fractional scaling.
  class RelMuEScaleSyst2017: public ISyst
  {
  public:
    RelMuEScaleSyst2017(const std::string& shortName,
                        const std::string& latexName,
                        double scale, double shiftmucat)
      : ISyst(shortName, latexName), fScale(scale), fShiftMuCat(shiftmucat) {}

    RelMuEScaleSyst2017(double scale, double shiftmucat)
      : ISyst("RelMuEScale2017", "Rel Muon Energy Scale 2017"),
        fScale(scale), fShiftMuCat(shiftmucat) {}

    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;

  private:
    const double fScale;
    const double fShiftMuCat;
  };

  extern const RelMuEScaleSyst2017 kRelMuEScaleSyst2017;

  //----------------------------------------------------------------------
  /// Absolute fit hadronic energy scale systematic
  class DirectHadEScaleSyst2017: public ISyst
  {
  public:
    DirectHadEScaleSyst2017(const std::string& shortName,
                            const std::string& latexName,
                            double scale)
      : ISyst(shortName, latexName), fScale(scale)
    {
    }

    DirectHadEScaleSyst2017(double scale)
      : ISyst("AbsHadEScale2017", "Absolute Hadronic Energy Scale 2017"),
        fScale(scale)
    {
    }

    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;

  private:
    const double fScale;
  };

  extern const DirectHadEScaleSyst2017 kDirectHadEScaleSyst2017;

  //----------------------------------------------------------------------
  /// Relative fit hadronic energy scale systematic
  class DirectRelHadEScaleSyst2017: public ISyst
  {
  public:
    DirectRelHadEScaleSyst2017(const std::string& shortName,
                               const std::string& latexName,
                               double scale)
      : ISyst(shortName, latexName), fScale(scale)
    {
    }

  DirectRelHadEScaleSyst2017(const double scale)
      : ISyst("RelHadEScale2017", "Relative Hadronic Energy Scale 2017"),
        fScale(scale)
    {
    }

    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;

  private:
    const double fScale;
  };

  extern const DirectRelHadEScaleSyst2017 kDirectRelHadEScaleSyst2017;
}