39 std::vector<std::string> NuMu_CutNames; std::vector<Cut> NuMu_TieredCuts;
40 std::vector<std::string> Nue_CutNames; std::vector<Cut> Nue_TieredCuts;
41 MakeCutFlowVecs( NuMu_CutNames, NuMu_TieredCuts, Nue_CutNames, Nue_TieredCuts );
43 size_t NuMuCuts = NuMu_TieredCuts.size();
44 size_t NueCuts = Nue_TieredCuts .size();
47 Spectrum* NuMu_ReconEnergy[NuMuCuts];
48 Spectrum* Nue_ReconEnergy [NueCuts ];
54 for(
size_t cut = 0;
cut < NuMuCuts; ++
cut) {
59 for(
size_t cut = 0;
cut < NueCuts; ++
cut) {
69 std::string OutName =
"CosmicCutFlow_"+sFHC+
".root";
70 TFile *
OutFile =
new TFile(OutName.c_str(),
"RECREATE");
72 for(
size_t cut = 0;
cut < NuMuCuts; ++
cut){
73 NuMu_ReconEnergy[
cut] ->
SaveTo( OutFile, TString(
"NuMuCCEnergy_")+TString(NuMu_CutNames[
cut]) ) ;
75 for(
size_t cut = 0;
cut < NueCuts; ++
cut) {
76 Nue_ReconEnergy [
cut] ->
SaveTo( OutFile, TString(
"NueCCEnergy_" )+TString(Nue_CutNames [
cut]) ) ;
Represent the binning of a Spectrum's x-axis.
Cuts and Vars for the 2020 FD DiF Study.
void SaveTo(const osc::IOscCalc &x, TDirectory *dir, const std::string &name)
void MakeCutFlowVecs(std::vector< std::string > &NuMu_CutNames, std::vector< Cut > &NuMu_TieredCuts, std::vector< std::string > &Nue_CutNames, std::vector< Cut > &Nue_TieredCuts)
void SetSpillCut(const SpillCut &cut)
void CutFlow_Cosmic(bool isFHC)
Representation of a spectrum in any variable, with associated POT.
const Cut kInCosmicTimingWindow
Is the event far from the start and ends of the spill ? For FD cosmic selection.
virtual void Go() override
Load all the registered spectra.
std::vector< float > Spectrum
const SystShifts kNoShift
Collaborates with Spectrum and OscillatableSpectrum to fill spectra from CAF files.
const Binning kNumuCCEOptimisedBinning
Optimised binning for numuCCE from L. Vinton. See docdb 16332. This was close to 'custC' in that talk...
const SpillCut kStandardSpillCuts
Apply this unless you're doing something special.
static Binning Simple(int n, double lo, double hi, const std::vector< std::string > &labels={})
const Var kUnweighted
The simplest possible Var, always 1. Used as a default weight.