includes_starplots.h

Go to the documentation of this file.

#include "3FlavorAna/Systs/GeniePCASyst.h"

std::vector<const ISyst *> MyGenieRwgtSysts = {
  GetGenieKnobSyst(rwgt::fReweightMaCCRES),
  GetGenieKnobSyst(rwgt::fReweightMaNCRES),
  GetGenieKnobSyst(rwgt::fReweightMvCCRES),
  GetGenieKnobSyst(rwgt::fReweightCCQEPauliSupViaKF),
  GetGenieKnobSyst(rwgt::fReweightFormZone),
};


std::vector<const ISyst *> MyGenieRwgtFrSysts = {
  GetGenieKnobSyst(rwgt::fReweightFrElas_N),
  GetGenieKnobSyst(rwgt::fReweightFrCEx_N),
  GetGenieKnobSyst(rwgt::fReweightFrAbs_N),
  GetGenieKnobSyst(rwgt::fReweightFrInel_pi),
};


std::vector<const ISyst *> MyMECSysts = {
  &kMECEnuShapeSyst2018Nu,
  &kMECEnuShapeSyst2018AntiNu,
  &kMECShapeSyst2018Nu,
  &kMECShapeSyst2018AntiNu,
  &kMECInitStateNPFracSyst2018Nu,
  &kMECInitStateNPFracSyst2018AntiNu
};

std::vector<const ISyst *> MyXSecSysts = {
  &kMAQEGenieReducedSyst2018,
  &kRPARESSyst2018,
  &kRPACCQEEnhSyst2018,
  &kCOHCCScaleSyst2018,
  &kCOHNCScaleSyst2018
};


std::vector<const ISyst *> MyGeniePCASysts = {
  GetGeniePrincipals2018Small(0),
  GetGeniePrincipals2018Small(1),
  GetGeniePrincipals2018Small(2),
  GetGeniePrincipals2018Small(3),
  GetGeniePrincipals2018Small(4),
};


std::vector<const ISyst *> MyFluxSysts = {
  GetFluxPrincipals2018(0),
  GetFluxPrincipals2018(1),
  GetFluxPrincipals2018(2),
  GetFluxPrincipals2018(3),
  GetFluxPrincipals2018(4)
};


std::vector<const ISyst *> MyCalibSysts = {
  &kAnaCalibrationSyst,
  &kAnaRelativeCalibSyst,  
  &kAnaCalibShapeSyst,
  &kAnaLightlevelSyst,
  &kAnaCherenkovSyst 
};


std::vector<const ISyst *> MyOtherSysts = {
  &kNeutronVisEScalePrimariesSyst2018,
  &kMuEScaleSyst2017,
  &kAna2018NormFHC,
  &kAna2018NormRHC,
  &kTauScaleSyst
};


// These two clamp th23 to the uo / lo
class FitSinSqTheta23UO: public IConstrainedFitVar
{
 public:
  virtual double GetValue(const osc::IOscCalcAdjustable* osc) const
  {
    return util::sqr(sin(osc->GetTh23()));
  };
  virtual void SetValue(osc::IOscCalcAdjustable* osc, double val) const
  {
    osc->SetTh23(asin(sqrt(Clamp(val))));
  };
  virtual std::string ShortName() const {return "sinsqth23_UO";}
  virtual std::string LatexName() const {return "sin^{2}#theta_{23} Constrained to Upper Octant";}

  // And now set max/min values for upper octant boundaries
  virtual double LowLimit() const {return 0.514;}
  virtual double HighLimit() const {return 1.;}
};
const FitSinSqTheta23UO kFitSinSqTheta23UO = FitSinSqTheta23UO();

class FitSinSqTheta23LO: public IConstrainedFitVar
{
 public:
  virtual double GetValue(const osc::IOscCalcAdjustable* osc) const
  {
    return util::sqr(sin(osc->GetTh23()));
  };
  virtual void SetValue(osc::IOscCalcAdjustable* osc, double val) const
  {
    osc->SetTh23(asin(sqrt(Clamp(val))));
  };
  virtual std::string ShortName() const {return "sinsqth23_LO";}
  virtual std::string LatexName() const {return "sin^{2}#theta_{23} Constraint to Lower Octant";}

  // And now set max/min values for lower octant boundaries
  virtual double LowLimit() const {return 0.;}
  //virtual double HighLimit() const {return 0.514;}
  virtual double HighLimit() const {return 0.49;}
};
const FitSinSqTheta23LO kFitSinSqTheta23LO = FitSinSqTheta23LO();