includes_starplots.h
Go to the documentation of this file.
2 
3 std::vector<const ISyst *> MyGenieRwgtSysts = {
4  GetGenieKnobSyst(rwgt::fReweightMaCCRES),
5  GetGenieKnobSyst(rwgt::fReweightMaNCRES),
6  GetGenieKnobSyst(rwgt::fReweightMvCCRES),
7  GetGenieKnobSyst(rwgt::fReweightCCQEPauliSupViaKF),
8  GetGenieKnobSyst(rwgt::fReweightFormZone),
9 };
10 
11 
12 std::vector<const ISyst *> MyGenieRwgtFrSysts = {
13  GetGenieKnobSyst(rwgt::fReweightFrElas_N),
14  GetGenieKnobSyst(rwgt::fReweightFrCEx_N),
15  GetGenieKnobSyst(rwgt::fReweightFrAbs_N),
16  GetGenieKnobSyst(rwgt::fReweightFrInel_pi),
17 };
18 
19 
20 std::vector<const ISyst *> MyMECSysts = {
27 };
28 
29 std::vector<const ISyst *> MyXSecSysts = {
35 };
36 
37 
38 std::vector<const ISyst *> MyGeniePCASysts = {
44 };
45 
46 
47 std::vector<const ISyst *> MyFluxSysts = {
53 };
54 
55 
56 std::vector<const ISyst *> MyCalibSysts = {
62 };
63 
64 
65 std::vector<const ISyst *> MyOtherSysts = {
71 };
72 
73 
74 // These two clamp th23 to the uo / lo
75 class FitSinSqTheta23UO: public IConstrainedFitVar
76 {
77  public:
78  virtual double GetValue(const osc::IOscCalcAdjustable* osc) const
79  {
80  return util::sqr(sin(osc->GetTh23()));
81  };
82  virtual void SetValue(osc::IOscCalcAdjustable* osc, double val) const
83  {
84  osc->SetTh23(asin(sqrt(Clamp(val))));
85  };
86  virtual std::string ShortName() const {return "sinsqth23_UO";}
87  virtual std::string LatexName() const {return "sin^{2}#theta_{23} Constrained to Upper Octant";}
88 
89  // And now set max/min values for upper octant boundaries
90  virtual double LowLimit() const {return 0.514;}
91  virtual double HighLimit() const {return 1.;}
92 };
94 
95 class FitSinSqTheta23LO: public IConstrainedFitVar
96 {
97  public:
98  virtual double GetValue(const osc::IOscCalcAdjustable* osc) const
99  {
100  return util::sqr(sin(osc->GetTh23()));
101  };
102  virtual void SetValue(osc::IOscCalcAdjustable* osc, double val) const
103  {
104  osc->SetTh23(asin(sqrt(Clamp(val))));
105  };
106  virtual std::string ShortName() const {return "sinsqth23_LO";}
107  virtual std::string LatexName() const {return "sin^{2}#theta_{23} Constraint to Lower Octant";}
108 
109  // And now set max/min values for lower octant boundaries
110  virtual double LowLimit() const {return 0.;}
111  //virtual double HighLimit() const {return 0.514;}
112  virtual double HighLimit() const {return 0.49;}
113 };
115 
116 
const NOvARwgtSyst kMECShapeSyst2018AntiNu("MECShape2018AntiNu","MEC 2018 (q_{0}, |#vec{q}|) response, antineutrinos", novarwgt::kMECQ0Q3RespSystNubar2018)
Definition: MECSysts.h:26
std::vector< const ISyst * > MyMECSysts
std::vector< const ISyst * > MyFluxSysts
std::vector< const ISyst * > MyGeniePCASysts
virtual void SetValue(osc::IOscCalcAdjustable *osc, double val) const
const NOvARwgtSyst kRPACCQEEnhSyst2018("RPAShapeenh2018","RPA shape: higher-Q^{2} enhancement (2018)", novarwgt::kRPACCQEEnhSyst2018)
Definition: RPASysts.h:13
const NOvARwgtSyst kMAQEGenieReducedSyst2018(genie::rew::GSyst::AsString(genie::rew::EGSyst(rwgt::fReweightMaCCQE))++"_reduced_2018", genie::rew::GSyst::AsString(genie::rew::EGSyst(rwgt::fReweightMaCCQE))++"_reduced_2018", novarwgt::kMAQEGenieReducedSyst2018)
2018 &#39;reduced&#39; M_A^QE shift. See documentation in NOvARwgt
Definition: XSecSysts.h:51
T sqrt(T number)
Definition: d0nt_math.hpp:156
const NuTruthSystComponentScale kTauScaleSyst("NuTauScale","#nu_{#tau} Scale", kIsTau_NT &&!kIsNC_NT, 0.6, NuTruthSystComponentScale::kLinear)
100% uncertainty scale on taus
Definition: Systs.h:176
const NOvARwgtSyst kMECEnuShapeSyst2018AntiNu("MECEnuShapeAntiNu","MEC E_{#nu} shape, antineutrinos", novarwgt::kMECEnuShapeSyst2017_NubarOnly)
Definition: MECSysts.h:32
const NOvARwgtSyst kMECShapeSyst2018Nu("MECShape2018Nu","MEC 2018 (q_{0}, |#vec{q}|) response, neutrinos", novarwgt::kMECQ0Q3RespSystNu2018)
Definition: MECSysts.h:25
virtual std::string LatexName() const
virtual std::string ShortName() const
const NOvARwgtSyst kMECInitStateNPFracSyst2018AntiNu("MECInitStateNPFracAntiNu","MEC initial state np fraction, antineutrinos", novarwgt::kMECInitStateNPFracSyst2017_NubarOnly)
Definition: MECSysts.h:35
std::vector< const ISyst * > MyCalibSysts
virtual T GetTh23() const
Definition: IOscCalc.h:94
const NeutronVisEScaleSyst2018 kNeutronVisEScalePrimariesSyst2018(false)
const DummyAnaSyst kAnaRelativeCalibSyst("RelativeCalib","RelCalib")
const NOvARwgtSyst kMECEnuShapeSyst2018Nu("MECEnuShapeNu","MEC E_{#nu} shape, neutrinos", novarwgt::kMECEnuShapeSyst2017_NuOnly)
Definition: MECSysts.h:31
T sqr(T x)
More efficient square function than pow(x,2)
Definition: MathUtil.h:23
BeamSyst * GetFluxPrincipals2018(int PCIdx)
Definition: BeamSysts.cxx:89
virtual std::string LatexName() const
const MuEScaleSyst2017 kMuEScaleSyst2017(0.0074, 0.0012)
const FitSinSqTheta23LO kFitSinSqTheta23LO
std::vector< const ISyst * > MyOtherSysts
const DummyAnaSyst kAnaLightlevelSyst("Lightlevel","Lightlevel")
const DummyAnaSyst kAnaCalibrationSyst("Calibration","AbsCalib")
const DummyAnaSyst kAnaCherenkovSyst("Cherenkov","Cherenkov")
virtual double LowLimit() const
std::vector< const ISyst * > MyGenieRwgtSysts
const NOvARwgtSyst kMECInitStateNPFracSyst2018Nu("MECInitStateNPFracNu","MEC initial state np fraction, neutrinos", novarwgt::kMECInitStateNPFracSyst2017_NuOnly)
Definition: MECSysts.h:34
std::vector< const ISyst * > MyGenieRwgtFrSysts
T Clamp(T val) const
Definition: IFitVar.h:60
virtual double GetValue(const osc::IOscCalcAdjustable *osc) const
const NOvARwgtSyst kCOHCCScaleSyst2018("COHCCScale2018","Coherent CC Scale", novarwgt::kCOHNormCCSystKnob)
20% uncertainty scale on coherent CC events. See documentation in NOvARwgt
Definition: XSecSysts.h:40
Oscillation probability calculators.
Definition: Calcs.h:5
GeniePCASyst * GetGeniePrincipals2018Small(int PCIdx)
virtual double HighLimit() const
virtual double LowLimit() const
const NOvARwgtSyst kCOHNCScaleSyst2018("COHNCScale2018","Coherent NC Scale", novarwgt::kCOHNormNCSystKnob)
20% uncertainty scale on coherent NC events. See documentation in NOvARwgt
Definition: XSecSysts.h:43
T sin(T number)
Definition: d0nt_math.hpp:132
const DummyAnaSyst kAna2018NormRHC("NormRHC2018","RHC. Norm.")
const FitSinSqTheta23UO kFitSinSqTheta23UO
virtual std::string ShortName() const
virtual void SetTh23(const T &th23)=0
std::vector< const ISyst * > MyXSecSysts
virtual double GetValue(const osc::IOscCalcAdjustable *osc) const
const NOvARwgtSyst * GetGenieKnobSyst(rwgt::ReweightLabel_t knobIdx, std::string altName, std::string altLabel)
Convenience function to get a GENIE knob syst. (Allows using the GENIE knob name & description as the...
Definition: XSecSysts.cxx:119
virtual void SetValue(osc::IOscCalcAdjustable *osc, double val) const
virtual double HighLimit() const
const DummyAnaSyst kAna2018NormFHC("NormFHC2018","FHC. Norm.")
const DummyAnaSyst kAnaCalibShapeSyst("CalibShape","CalibShape")
T asin(T number)
Definition: d0nt_math.hpp:60
const NOvARwgtSyst kRPARESSyst2018("RPAShapeRES2018","RPA shape: resonance events (2018)", novarwgt::kRPARESSyst2018)
Definition: RPASysts.h:21
enum BeamMode string