3 #include "CAFAna/Core/Cut.h" const NuTruthCut kTrueNueCC_Nu([](const caf::SRNeutrinoProxy *truth){return( truth->iscc && abs(truth->pdg)==12 );})
Cuts and Vars for the 2020 FD DiF Study.
const NuTruthCut kOverlaidData_Nu([](const caf::SRNeutrinoProxy *nu){return true;})
const Cut kAnaTrueECut_Numu
const NuTruthCut kOverlaidMC_Nu([](const caf::SRNeutrinoProxy *nu){return(nu->mode > 9999);})
const NuTruthCut kAnaTrueECut_Nu_Nue([](const caf::SRNeutrinoProxy *nu){return(nu->E >=0.&&nu->E< 6.0);})
const NuTruthCut kNoCut_NT([](const caf::SRNeutrinoProxy *nu){return true;})
const NuTruthCut kTrueNumuCC_Nu([](const caf::SRNeutrinoProxy *truth){return( truth->iscc && abs(truth->pdg)==14 );})
const Cut kTrueOFiducialLoose
const NuTruthCut kTrueOFiducial_Nu([](const caf::SRNeutrinoProxy *nu){return(nu->vtx.X()< ovtxmax.X()&& nu->vtx.X() > ovtxmin.X()&& nu->vtx.Y() > ovtxmin.Y()&& nu->vtx.Y()< ovtxmax.Y()&& nu->vtx.Z() > ovtxmin.Z()&& nu->vtx.Z()< ovtxmax.Z());})
const Cut kOverlaidNoEffCut([](const caf::SRProxy *sr){if(sr->mc.nallnus==0) return false;for(int j=0;j< sr->mc.nallnus;j++) if((sr->mc.allnus[j].mode > 9999 ||!sr->spill.ismc)) return true;return false;})
const Cut kAnaTrueECut_Nue
const SpillCut kNoCut_Sp([](const caf::SRSpillProxy *sp){return true;})
const Cut kOverlaidEff51([](const caf::SRProxy *sr){if(sr->mc.nallnus==0) return false;for(int j=0;j< sr->mc.nallnus;j++) if((sr->mc.allnus[j].mode > 9999
||!sr->spill.ismc)
&&sr->mc.allnus[j].eff > 0.51) return true;return false;})
const NuTruthCut kTrueE250MeV_Nu([](const caf::SRNeutrinoProxy *nu){return nu->E > 0.25;})
const NuTruthCut kTrueOFiducialLoose_Nu([](const caf::SRNeutrinoProxy *nu){return(nu->vtx.X()< ovtxmaxLoose.X()&& nu->vtx.X() > ovtxminLoose.X()&& nu->vtx.Y() > ovtxminLoose.Y()&& nu->vtx.Y()< ovtxmaxLoose.Y()&& nu->vtx.Z() > ovtxminLoose.Z()&& nu->vtx.Z()< ovtxmaxLoose.Z());})
Template for Cut and SpillCut.
const Cut kNotOverlaidData([](const caf::SRProxy *sr){if(sr->mc.nallnus==0) return true;return false;})
const NuTruthCut kAnaTrueECut_Nu_Numu([](const caf::SRNeutrinoProxy *nu){return(nu->E >=0.&&nu->E< 7.);})
1.8.11
