XSecSysts.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"

#include "CAFAna/Systs/SystComponentScale.h"

// Names for the GENIE systs
#include "nugen/NuReweight/ReweightLabels.h"

#include "NOvARwgt/rwgt/ISystKnob.h"

namespace ana
{
  //----------------------------------------------------------------------------------

  /// Base class for systs actually implemented in NOvARwgt.
  class NOvARwgtSyst : public ISyst
  {
    public:
      NOvARwgtSyst(std::string shortName, std::string latexName, const novarwgt::ISystKnob * knob)
        : ISyst(shortName, latexName),  // todo: should this be pulling the shortName from the knob's Name()?
          fNOvARwgtKnob(knob)
      {}

      void TruthShift(double sigma, caf::SRNeutrinoProxy* nu, double& weight) const override;

      const novarwgt::ISystKnob * GetNOvARwgtKnob() const { return fNOvARwgtKnob; }

    private:
      const novarwgt::ISystKnob * fNOvARwgtKnob;
  };

  //----------------------------------------------------------------------

  /// Convenience function to get a GENIE knob syst.  (Allows using the GENIE knob name & description as the default.)
  const NOvARwgtSyst * GetGenieKnobSyst(rwgt::ReweightLabel_t knobIdx, std::string altShortName="", std::string altLatexName="");

  //----------------------------------------------------------------------
  /// \brief 20% uncertainty scale on coherent CC events.  See documentation in NOvARwgt
  extern const NOvARwgtSyst kCOHCCScaleSyst2018;

  /// \brief 20% uncertainty scale on coherent NC events.  See documentation in NOvARwgt
  extern const NOvARwgtSyst kCOHNCScaleSyst2018;

  //----------------------------------------------------------------------
  /// 2017 'reduced' M_A^QE shift.  See documentation in NOvARwgt.
  extern const NOvARwgtSyst kMAQEGenieReducedSyst2017;

  //----------------------------------------------------------------------
  /// 2018 'reduced' M_A^QE shift.  See documentation in NOvARwgt
  extern const NOvARwgtSyst kMAQEGenieReducedSyst2018;

  //----------------------------------------------------------------------
  ///@{
  /// 2020 z-expansion knobs with proper correlations.  See documentation in NOvARwgt
  extern const NOvARwgtSyst kZExpEV1Syst2020;
  extern const NOvARwgtSyst kZExpEV2Syst2020;
  extern const NOvARwgtSyst kZExpEV3Syst2020;
  extern const NOvARwgtSyst kZExpEV4Syst2020;
  ///@}

  //----------------------------------------------------------------------
  /// Radiative corrections syst (nues).  See documentation in NOvARwgt
  extern const NOvARwgtSyst    kRadCorrNue;

  /// Radiative corrections syst (nuebars).  See documentation in NOvARwgt
  extern const NOvARwgtSyst kRadCorrNuebar;

  //----------------------------------------------------------------------
  /// Second-class current syst.  See documentation in NOvARwgt
  extern const NOvARwgtSyst  k2ndClassCurrs;


  //----------------------------------------------------------------------
  /// Nuclear Optical Potential. See documentation in NOvARwgt
  extern const NOvARwgtSyst kNucOptPtlSyst;
} // namespace