14 #include "CAFAna/Core/SpectrumLoaderMockData.h" 27 #include "TFitResultPtr.h" 32 #include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Bins.h" 33 #include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Cuts.h" 34 #include "/nova/app/users/acedeno/tag_releaseS17-09-24/NDAna/ncpi0_Systematics/ncpi0Vars.h" 49 const std::string MC =
"dataset_def_name_newest_snapshot prod_caf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_v1";
94 TFile *
out =
new TFile(
"xyzcontainmentstudy1DSpectrumsApril2018.root",
"RECREATE");
98 ndMC.
SaveTo(out->mkdir(
"ndMC"));
99 ndMC3P.
SaveTo(out->mkdir(
"ndMC3P"));
100 ndNCPi0.
SaveTo(out->mkdir(
"ndNCPi0"));
101 ndCCPi0.
SaveTo(out->mkdir(
"ndCCPi0"));
102 ndNCNonPi0.
SaveTo(out->mkdir(
"ndNCNonPi0"));
103 ndNCPi0L.
SaveTo(out->mkdir(
"ndNCPi0L"));
104 ndCCNonPi0.
SaveTo(out->mkdir(
"ndCCNonPi0"));
106 ndMCy.
SaveTo(out->mkdir(
"ndMCy"));
107 ndMC3Py.
SaveTo(out->mkdir(
"ndMC3Py"));
108 ndNCPi0y.
SaveTo(out->mkdir(
"ndNCPi0y"));
109 ndCCPi0y.
SaveTo(out->mkdir(
"ndCCPi0y"));
110 ndNCNonPi0y.
SaveTo(out->mkdir(
"ndNCNonPi0y"));
111 ndNCPi0Ly.
SaveTo(out->mkdir(
"ndNCPi0Ly"));
112 ndCCNonPi0y.
SaveTo(out->mkdir(
"ndCCNonPi0y"));
115 ndMCz.
SaveTo(out->mkdir(
"ndMCz"));
116 ndMC3Pz.
SaveTo(out->mkdir(
"ndMC3Pz"));
117 ndNCPi0z.
SaveTo(out->mkdir(
"ndNCPi0z"));
118 ndCCPi0z.
SaveTo(out->mkdir(
"ndCCPi0z"));
119 ndNCNonPi0z.
SaveTo(out->mkdir(
"ndNCNonPi0z"));
120 ndNCPi0Lz.
SaveTo(out->mkdir(
"ndNCPi0Lz"));
121 ndCCNonPi0z.
SaveTo(out->mkdir(
"ndCCNonPi0z"));
Cuts and Vars for the 2020 FD DiF Study.
const Var ContStopZ([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return-10.0f;if(sr->vtx.nelastic==0) return-10.0f;return float(sr->vtx.elastic[0].fuzzyk.png[0].shwlid.stop.z);})
const Var ContStopX([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return-10.0f;if(sr->vtx.nelastic==0) return-10.0f;return float(sr->vtx.elastic[0].fuzzyk.png[0].shwlid.stop.x);})
Representation of a spectrum in any variable, with associated POT.
const Var ContStopY([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return-10.0f;if(sr->vtx.nelastic==0) return-10.0f;return float(sr->vtx.elastic[0].fuzzyk.png[0].shwlid.stop.y);})
virtual void Go() override
Load all the registered spectra.
void SaveTo(TDirectory *dir, const std::string &name) const
const Cut k2ProngPreSel2CVN
const Binning bins_for_vtxXY
const Cut kPreSelect
3 prong analysis cuts////
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const Binning bins_for_vtxZ