28 std::string polarity =
"fhc";
if (!isFHC) polarity =
"rhc";
32 std::cout <<
"\n================================= \n" 33 <<
"\n\t isFHC --> " << isFHC <<
" --> " << polarity
34 <<
"\n\t Period --> " <<
period 35 <<
"\n\t Far detector isFD --> " << detector
36 <<
"\n================================= \n" 40 std::vector<std::pair<Cut, std::string> > MyCuts =
BasicPIDPlots_Cuts( isFHC,
true,
false );
41 const size_t NCuts = MyCuts.size();
47 std::string MyNonSwapDef =
"prod_caf_R19-11-18-prod5reco.f_fd_genie_N1810j0211a_nonswap_" +polarity+
"_nova_v08_"+
period+
"_v1";
48 std::string MyFluxSwapDef =
"prod_caf_R19-11-18-prod5reco.f_fd_genie_N1810j0211a_fluxswap_"+polarity+
"_nova_v08_"+period+
"_v1";
49 std::string MyTauSwapDef =
"prod_caf_R19-11-18-prod5reco.f_fd_genie_N1810j0211a_tau_" +polarity+
"_nova_v08_"+period+
"_v1";
68 Cut kThisCut = MyCuts[
cc].first;
92 std::string OutName =
"BasicPIDPlots_2020_"+detector+
"_"+data+
"_"+polarity+
"_"+period+
".root";
93 TFile *
OutFile = TFile::Open(OutName.c_str(),
"RECREATE");
97 sNuRecoE [
cc] ->
SaveTo( OutFile, TString(
"RecoNuE" )+TString(ThisCutNa) ) ;
109 NoCutRecoNu ->
SaveTo( OutFile, TString(
"NoCutRecoNu" ) ) ;
Far Detector at Ash River.
Cuts and Vars for the 2020 FD DiF Study.
void SaveTo(const osc::IOscCalc &x, TDirectory *dir, const std::string &name)
Collection of SpectrumLoaders for many configurations.
std::vector< std::pair< Cut, std::string > > BasicPIDPlots_Cuts(bool isFHC, bool isFD, bool isData)
void SetSpillCut(const SpillCut &cut)
Representation of a spectrum in any variable, with associated POT.
void Go()
Call Go() on all the loaders.
const XML_Char const XML_Char * data
void SetSpillCut(const SpillCut &cut)
void NuMu2020_BasicPIDPlots_FD(bool isFHC, std::string period="full")
std::vector< float > Spectrum
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.
void SetLoaderPath(const std::string &path, caf::Det_t det, DataMC datamc, DataSource src=kBeam, SwappingConfig swap=kNonSwap)
Configure loader via wildcard path.
const Cut kNoCut
The simplest possible cut: pass everything, used as a default.