13 #include "CAFAna/Cuts/NueCutsSecondAna.h" 14 #include "CAFAna/Cuts/NumuCuts.h" 20 #include "CAFAna/Systs/NumuSysts.h" 23 #include "CAFAna/Systs/NueExtrapSysts.h" 39 TString
fileName =
"pred_tutorial_nue_numu_rhc-test.root";
47 std::string ldrFDNonSwap =
"prod_decaf_R17-03-01-prod3reco.l_fd_genie_nonswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
49 std::string ldrFDFlxSwap =
"prod_decaf_R17-03-01-prod3reco.l_fd_genie_fluxswap_rhc_nova_v08_full_nue_or_numu_or_nus_contain_v1";
51 ldrFDNonSwap =
"def_snapshot " + ldrFDNonSwap +
" with stride 100";
52 ldrFDFlxSwap =
"defname: "+ldrFDFlxSwap +
" with stride 100";
68 const HistAxis nueAxis = kNueAxisSecondAnaCVN;
70 const Cut nueCutFD = kNueSAFDCVNSsb;
87 auto numuPrediction = numuGenerator.
Generate(loaders);
88 auto nuePrediction = nueGenerator.
Generate(loaders);
104 TFile *
file =
new TFile (fileName,
"recreate");
106 numuPrediction->SaveTo(file->mkdir(
"prediction_numu"));
107 nuePrediction->SaveTo(file->mkdir(
"prediction_nue"));
108 numuPredSysts->SaveTo(file->mkdir(
"prediction_systs_numu"));
109 nuePredSysts->
SaveTo(file->mkdir(
"prediction_systs_nue"));
Implements systematic errors by interpolation between shifted templates.
Far Detector at Ash River.
Cuts and Vars for the 2020 FD DiF Study.
const FDHadEnergyScaleSyst kFDHadEnergyScaleSyst(.05)
void make_prediction_rhc()
Collection of SpectrumLoaders for many configurations.
const NOvARwgtSyst kMECScaleSystSA("mecScale","MEC Scale", novarwgt::kMECScaleSystSA)
osc::IOscCalcAdjustable * DefaultOscCalc()
Create a new calculator with default assumptions for all parameters.
void Go()
Call Go() on all the loaders.
void SetSpillCut(const SpillCut &cut)
virtual void SaveTo(TDirectory *dir, const std::string &name) const override
const HistAxis kNumuCCAxis("Reconstructed Neutrino Energy (GeV)", kNumuEnergyBinning, kCCE)
const BeamSyst kBeamAllTransport((FindCAFAnaDir()+"/data/beam/TABeamSyst_2018Dec11.root"),"totErr","beamTransportComb","Combined Beam Transport Systematics")
All Beam Transport systematics combined in quadratures.
static const BeamSyst & kBeamAll
const SpillCut kStandardSpillCuts
Apply this unless you're doing something special.
void SetLoaderPath(const std::string &path, caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Configure loader via wildcard path.
virtual IOscCalc * Copy() const override