16 std::string production =
"sumdecaf_R17-11-14-prod4reco";
17 if(polarity ==
"fhc") production +=
".d";
18 if(polarity ==
"rhc") production +=
".e";
30 std::string fd_genie_nonswap =
"prod_" + production +
"_fd_genie_nonswap_" + polarity +
"_nova_v08_" +
period +
"_v1_numu2018";
33 std::string OutFileName = OutDir +
"quantiles__" + polarity +
"_" + period +
"__numu2018.root";
35 TFile *
f = TFile::Open(OutFileName.c_str(),
"RECREATE");
36 hQuantiles->Write(
"FDSpec2D");
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
Cuts and Vars for the 2020 FD DiF Study.
const HistAxis kHadEFracAxis("E_{had.} / E_{#nu}", Binning::Simple(200, 0, 1), kHadEFrac)
HistAxis that implements Hadronic Energy fraction binning used by L.Vinton to derive Hadronic Energy ...
const Var kXSecCVWgt2018Fix
void make_ehade_histogram(std::string polarity="rhc", std::string period="full")
const HistAxis kNumuCCOptimisedAxis("Reconstructed Neutrino Energy (GeV)", kNumuCCEOptimisedBinning, kCCE)
HistAxis that implements optimised numuCCE from L. Vinton. See docdb 16332. This was close to 'custC'...
TH2 * MakeQuantileHistogram(SpectrumLoader &loader, const HistAxis &independentAxis, const HistAxis &quantileAxis, const Cut &cut, const SpillCut &spillCut, const SystShifts &shift, const Var &wei)
Returns a TH2D with xAxis as independentAxis and yAxis as quantileAxis Quantiles in yAxis can be calc...
const SystShifts kNoShift
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const SpillCut kStandardSpillCuts
Apply this unless you're doing something special.