#include "sbnana/SBNAna/Cuts/NueCuts.h"
#include "sbnana/SBNAna/Vars/NueVars.h"
#include "sbnanaobj/StandardRecord/Proxy/SRProxy.h"
#include "sbnana/SBNAna/Cuts/VolumeDefinitions.h"

Go to the source code of this file.

Namespaces
	ana
	Oscillation analysis framework, runs over CAF files outside of ART.

Variables
const Cut	ana::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);})

const Cut	ana::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);})

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

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

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

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

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

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

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

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

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

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

const Cut	ana::kShowerRazzleElectronCut = kRecoShowerRazzlePID == 11

const Cut	ana::kShowerRazzleElectronScoreCut = kRecoShowerRazzleElectronScore > 0.8

const Cut	ana::kNueLongestTrackDazzleMuonCut = kLongestTrackDazzlePID != 13

const Cut	ana::kNueLongestTrackDazzleMuonScoreCut = kLongestTrackDazzleMuonScore < 0.8

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

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

const Cut	ana::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);})

const Cut	ana::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);})

Namespaces

Variables