6 #include "CAFAna/Core/Binning.h" 25 if ( isFHC) polarity =
"fhc";
26 if (!isFHC) polarity =
"rhc";
28 std::cout <<
"\n================================= \n" 29 <<
" make distribution for quantiles" 30 <<
"\n================================= \n" 31 <<
" polarity: " << polarity <<
", period: " <<
period 32 <<
"\n================================= \n" 42 std::string MyDef =
"prod_caf_R19-11-18-prod5reco.f_fd_genie_N1810j0211a_nonswap_"+polarity+
"_nova_v08_"+
period+
"_v1";
53 TH2* hQuantiles = sQuantiles .
ToTH2( sQuantiles .
POT() );
56 std::string Out_OldPresel =
"quantiles_" + polarity +
"_" + period +
"_numu2020.root";
57 TFile *
OutFile = TFile::Open(Out_OldPresel.c_str(),
"RECREATE");
59 hQuantiles ->
Write(
"FDSpec2D" );
64 hQuantiles->Draw(
"colz");
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
TH2 * ToTH2(double exposure, EExposureType expotype=kPOT, EBinType bintype=kBinContent) const
Spectrum must be 2D to obtain TH2.
Cuts and Vars for the 2020 FD DiF Study.
void SetSpillCut(const SpillCut &cut)
const HistAxis kNumuCCOptimisedAxis2020("Reconstructed Neutrino Energy (GeV)", kNumuCCEOptimisedBinning, kNumuE2020)
Representation of a spectrum in any variable, with associated POT.
const HistAxis kHadEFracAxis2020("E_{had.} / E_{#nu}", Binning::Simple(200, 0, 1), kNumuHadEFrac2020)
void make_quantiles_histogram_2020(bool isFHC, std::string period="full")
virtual void Go() override
Load all the registered spectra.
const SystShifts kNoShift
std::vector< double > POT
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const SpillCut kStandardSpillCuts
Apply this unless you're doing something special.