de/dbd/NueCuts_8h_source.html

 #pragma once


 #include "sbnana/CAFAna/Core/Cut.h"

 #include "sbnana/SBNAna/Cuts/Cuts.h"

 #include "sbnana/SBNAna/Vars/NueVars.h"


 namespace ana {

 extern const Cut kRecoShower;

 extern const Cut kNueBasicCut;


 extern const Cut kNueHasTrackCut;

 extern const Cut kNueTrackContainmentCut;

 extern const Cut kNueTrackLenCut;

 extern const Cut kNueMuonCutOLD;

 extern const Cut kNueMuonCutNEW;

 extern const Cut kNueMuonCut;


 extern const Cut kNueNumShowersCut;


 extern const Cut kShowerEnergyCut;

 extern const Cut kShowerdEdxCut;

 extern const Cut kShowerConvGapCut;

 extern const Cut kShowerDensityCut;

 extern const Cut kShowerOpenAngleCut;


 extern const Cut kShowerRazzleElectronCut;

 extern const Cut kShowerRazzleElectronScoreCut;

 extern const Cut kNueLongestTrackDazzleMuonCut;

 extern const Cut kNueLongestTrackDazzleMuonScoreCut;

 extern const Cut kNueAllTrackDazzleMuonCut;

 extern const Cut kNueAllTrackDazzleMuonScoreCut;


 extern const Cut kNueContainedND;

 extern const Cut kNueContainedFD;


 const Cut kPreNueSelND = kFiducialVolumeND && kSlcNuScoreCut && kSlcFlashMatchCut;

 const Cut kRecoNueSel = kRecoShower && kShowerEnergyCut;

 const Cut kFullNueSel = kNueTrackLenCut && kShowerConvGapCut && kShowerdEdxCut && kShowerDensityCut;

 const Cut kRazzleDazzleNueSel = kNueAllTrackDazzleMuonCut && kShowerRazzleElectronScoreCut;


 // Specific FD cuts although they're currenly basically the same as for the ND

 const Cut kNueFlashScoreFDCut   = kSlcFlashMatchCut;

 const Cut kNuePandoraScoreFDCut = kSlcNuScoreCut;

 const Cut kRecoShowerFD = kRecoShower && kNueNumShowersCut && kShowerdEdxCut && kShowerConvGapCut && kNueTrackLenCut && kShowerDensityCut && kShowerEnergyCut;

 const Cut kNueFD        = kContainedFD && kNueFlashScoreFDCut && kNuePandoraScoreFDCut && kRecoShowerFD;


 }

Cut.h

ana::kShowerConvGapCut
const Cut kShowerConvGapCut([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].conversion_gap< 3.25f);})
Definition: NueCuts.h:22

ana::kShowerDensityCut
const Cut kShowerDensityCut([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].density > 4.5);})
Definition: NueCuts.h:23

ana::kNueContainedFD
const Cut kNueContainedFD([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;double this_endx=slc->reco.shw[largestShwIdx].start.x+(slc->reco.shw[largestShwIdx].dir.x *slc->reco.shw[largestShwIdx].len);double this_endy=slc->reco.shw[largestShwIdx].start.y+(slc->reco.shw[largestShwIdx].dir.y *slc->reco.shw[largestShwIdx].len);double this_endz=slc->reco.shw[largestShwIdx].start.z+(slc->reco.shw[largestShwIdx].dir.z *slc->reco.shw[largestShwIdx].len);bool startx=(fvfd_cryo1.xmin< slc->reco.shw[largestShwIdx].start.x)&&(slc->reco.shw[largestShwIdx].start.x< fvfd_cryo1.xmax);bool endx=(fvfd_cryo1.xmin< this_endx)&&(this_endx< fvfd_cryo1.xmax);bool starty=(fvfd_cryo1.ymin< slc->reco.shw[largestShwIdx].start.y)&&(slc->reco.shw[largestShwIdx].start.y< fvfd_cryo1.ymax);bool endy=(fvfd_cryo1.ymin< this_endy)&&(this_endy< fvfd_cryo1.ymax);bool startz=(fvfd_cryo1.zmin< slc->reco.shw[largestShwIdx].start.z)&&(slc->reco.shw[largestShwIdx].start.z< fvfd_cryo1.zmax);bool endz=(fvfd_cryo1.zmin< this_endz)&&(this_endz< fvfd_cryo1.zmax);return(startx &&endx &&starty &&endy &&startz &&endz);})
Definition: NueCuts.h:34

ana::kNueMuonCutOLD
const Cut kNueMuonCutOLD

ana::kPreNueSelND
const Cut kPreNueSelND
Definition: NueCuts.h:36

ana::kSlcFlashMatchCut
const Cut kSlcFlashMatchCut([](const caf::SRSliceProxy *slc){return(kSlcHasFlashMatch(slc)&&slc->fmatch.score >0 &&slc->fmatch.score< 6);})
Definition: SBNAna/Cuts/Cuts.h:28

NueVars.h

ana
process_name opflashCryoW ana
Definition: stage0_icarus_light_only.fcl:60

ana::kShowerEnergyCut
const Cut kShowerEnergyCut([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_energy > 0.2f);})
Definition: NueCuts.h:20

ana::kShowerRazzleElectronScoreCut
const Cut kShowerRazzleElectronScoreCut
Definition: NueCuts.cxx:129

ana::kContainedFD
const Cut kContainedFD
Definition: SBNAna/Cuts/Cuts.h:30

ana::kFullNueSel
const Cut kFullNueSel
Definition: NueCuts.h:38

ana::kNueAllTrackDazzleMuonCut
const Cut kNueAllTrackDazzleMuonCut([](const caf::SRSliceProxy *slc){for(auto const &trk:slc->reco.trk){return false;}return true;})
Definition: NueCuts.h:30

ana::kFiducialVolumeND
const Cut kFiducialVolumeND([](const caf::SRSliceProxy *slc){return PtInVolAbsX(slc->vertex, fvndAbs);})
Definition: SBNAna/Cuts/Cuts.h:17

ana::kNueFlashScoreFDCut
const Cut kNueFlashScoreFDCut
Definition: NueCuts.h:42

Cuts.h

ana::kNueContainedND
const Cut kNueContainedND([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;double this_endx=slc->reco.shw[largestShwIdx].start.x+(slc->reco.shw[largestShwIdx].dir.x *slc->reco.shw[largestShwIdx].len);double this_endy=slc->reco.shw[largestShwIdx].start.y+(slc->reco.shw[largestShwIdx].dir.y *slc->reco.shw[largestShwIdx].len);double this_endz=slc->reco.shw[largestShwIdx].start.z+(slc->reco.shw[largestShwIdx].dir.z *slc->reco.shw[largestShwIdx].len);bool startx=(fvndAbs.xmin< std::abs(slc->reco.shw[largestShwIdx].start.x))&&(std::abs(slc->reco.shw[largestShwIdx].start.x)< fvndAbs.xmax);bool endx=(fvndAbs.xmin< std::abs(this_endx))&&(std::abs(this_endx)< fvndAbs.xmax);bool starty=(fvndAbs.ymin< slc->reco.shw[largestShwIdx].start.y)&&(slc->reco.shw[largestShwIdx].start.y< fvndAbs.ymax);bool endy=(fvndAbs.ymin< this_endy)&&(this_endy< fvndAbs.ymax);bool startz=(fvndAbs.zmin< slc->reco.shw[largestShwIdx].start.z)&&(slc->reco.shw[largestShwIdx].start.z< fvndAbs.zmax);bool endz=(fvndAbs.zmin< this_endz)&&(this_endz< fvndAbs.zmax);return(startx &&endx &&starty &&endy &&startz &&endz);})
Definition: NueCuts.h:33

ana::kNueMuonCut
const Cut kNueMuonCut([](const caf::SRSliceProxy *slc){return(kLongestTrackLength(slc)< 80.f||kLongestTrackChi2Muon(slc) > 30.f||kLongestTrackChi2Proton(slc)< 60.f);})
Definition: NueCuts.h:16

ana::kNueHasTrackCut
const Cut kNueHasTrackCut([](const caf::SRSliceProxy *slc){return slc->reco.ntrk > 0;})
Definition: NueCuts.h:11

ana::kSlcNuScoreCut
const Cut kSlcNuScoreCut([](const caf::SRSliceProxy *slc){return(kSlcIsRecoNu(slc)&&slc->nu_score >0.4);})
Definition: SBNAna/Cuts/Cuts.h:25

ana::kNueTrackLenCut
const Cut kNueTrackLenCut([](const caf::SRSliceProxy *slc){return kLongestTrackLength(slc)< 110.f;})
Definition: NueCuts.h:13

ana::kNueBasicCut
const Cut kNueBasicCut([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_energy > 0.2f &&slc->reco.shw[largestShwIdx].bestplane_dEdx< 3 &&slc->reco.shw[largestShwIdx].conversion_gap< 3);})
Definition: NueCuts.h:9

ana::kNuePandoraScoreFDCut
const Cut kNuePandoraScoreFDCut
Definition: NueCuts.h:43

ana::kRecoNueSel
const Cut kRecoNueSel
Definition: NueCuts.h:37

ana::kNueMuonCutNEW
const Cut kNueMuonCutNEW

ana::kRazzleDazzleNueSel
const Cut kRazzleDazzleNueSel
Definition: NueCuts.h:39

ana::kShowerdEdxCut
const Cut kShowerdEdxCut([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_dEdx< 3.625f);})
Definition: NueCuts.h:21

ana::kNueLongestTrackDazzleMuonScoreCut
const Cut kNueLongestTrackDazzleMuonScoreCut
Definition: NueCuts.cxx:132

ana::kRecoShower
const Cut kRecoShower([](const caf::SRSliceProxy *slc){const int largestShwIdx(kLargestRecoShowerIdx(slc));return false;return(slc->reco.shw[largestShwIdx].bestplane_energy > 0.f &&slc->reco.shw[largestShwIdx].bestplane_dEdx > 0.f &&slc->reco.shw[largestShwIdx].conversion_gap > 0.f);})
Definition: NueCuts.h:8

ana::kRecoShowerFD
const Cut kRecoShowerFD
Definition: NueCuts.h:44

ana::kShowerOpenAngleCut
const Cut kShowerOpenAngleCut([](const caf::SRSliceProxy *slc){return kRecoShower_OpenAngle(slc)< 12.f;})
Definition: NueCuts.h:24

ana::kShowerRazzleElectronCut
const Cut kShowerRazzleElectronCut
Definition: NueCuts.cxx:128

ana::_Cut
Template for Cut and SpillCut.
Definition: Cut.h:16

ana::kNueFD
const Cut kNueFD
Definition: NueCuts.h:45

ana::kNueAllTrackDazzleMuonScoreCut
const Cut kNueAllTrackDazzleMuonScoreCut([](const caf::SRSliceProxy *slc){for(auto const &trk:slc->reco.trk){return false;}return true;})
Definition: NueCuts.h:31

ana::kNueTrackContainmentCut
const Cut kNueTrackContainmentCut([](const caf::SRSliceProxy *slc){const int longestTrackIdx(kLongestTrackIdx(slc));return true;return PtInVol(slc->reco.trk[longestTrackIdx].end, fvndExit);})
Definition: NueCuts.h:12

ana::kNueNumShowersCut
const Cut kNueNumShowersCut([](const caf::SRSliceProxy *slc){return((unsigned int) kRecoShowers_EnergyCut(slc)==1);})
Definition: NueCuts.h:18

ana::kNueLongestTrackDazzleMuonCut
const Cut kNueLongestTrackDazzleMuonCut
Definition: NueCuts.cxx:131