db/db5/NumuVarsIcarus202106_8h_source.html

 #pragma once


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


 namespace ana

 {

   extern const Var kTruthIndex;


   extern const Var kPrimaryEnergy;


   extern const Var kMuMaxTrack;


   extern const Var kPTrackInd;


   extern const Var kRecoMuonP;


   extern const Var kTrueMuonP;

 }

ana::kTrueMuonP
const Var kTrueMuonP([](const caf::SRSliceProxy *slc) -> float{float p(-5.f);{auto const &trk=slc->reco.trk.at(kPTrackInd(slc));p=std::hypot(trk.truth.p.genp.x, trk.truth.p.genp.y, trk.truth.p.genp.z);}return p;})
Definition: NumuVarsIcarus202106.h:17

ana::kRecoMuonP
const Var kRecoMuonP([](const caf::SRSliceProxy *slc) -> float{float p(-5.f);{auto const &trk=slc->reco.trk.at(kPTrackInd(slc));const bool Contained=(!isnan(trk.end.x)&&(trk.end.x< -71.1-25 &&trk.end.x >-369.33+25)&&!isnan(trk.end.y)&&(trk.end.y >-181.7+25 &&trk.end.y< 134.8-25)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+30 &&trk.end.z< 895.95-50));if(Contained) p=trk.rangeP.p_muon;else p=trk.mcsP.fwdP_muon;}return p;})
Definition: NumuVarsIcarus202106.h:15

ana::kMuMaxTrack
const Var kMuMaxTrack([](const caf::SRSliceProxy *slc) -> float{float len(-5.f);for(auto const &trk:slc->reco.trk){if((trk.truth.p.pdg==13||trk.truth.p.pdg==-13)&&trk.truth.p.length > len) len=trk.truth.p.length;}return len;})
Definition: NumuVarsIcarus202106.h:11

ana
process_name opflashCryoW ana
Definition: stage0_icarus_light_only.fcl:60

ana::Var
_Var< caf::SRSliceProxy > Var
Representation of a variable to be retrieved from a caf::StandardRecord object.
Definition: Var.h:73

ana::kPrimaryEnergy
const Var kPrimaryEnergy
Definition: NumuVarsIcarus202106.cxx:9

ana::kTruthIndex
const Var kTruthIndex
Definition: NumuVarsIcarus202106.cxx:7

Var.h

ana::_Var< caf::SRSliceProxy >

ana::kPTrackInd
const Var kPTrackInd([](const caf::SRSliceProxy *slc) -> int{float Longest(0);int PTrackInd(-1);for(std::size_t i(0);i< slc->reco.trk.size();++i){auto const &trk=slc->reco.trk.at(i);if(trk.bestplane==-1) continue;const float Atslc=std::hypot(slc->vertex.x-trk.start.x, slc->vertex.y-trk.start.y, slc->vertex.z-trk.start.z);const bool AtSlice=(Atslc< 10.0 &&trk.pfp.parent_is_primary);const float Chi2Proton=trk.chi2pid[trk.bestplane].chi2_proton;const float Chi2Muon=trk.chi2pid[trk.bestplane].chi2_muon;const bool Contained=(!isnan(trk.end.x)&&(trk.end.x< -71.1-25 &&trk.end.x >-369.33+25)&&!isnan(trk.end.y)&&(trk.end.y >-181.7+25 &&trk.end.y< 134.8-25)&&!isnan(trk.end.z)&&(trk.end.z >-895.95+30 &&trk.end.z< 895.95-50));const bool MaybeMuonExiting=(!Contained &&trk.len > 100);const bool MaybeMuonContained=(Contained &&Chi2Proton > 60 &&Chi2Muon< 30 &&trk.len > 50);{Longest=trk.len;PTrackInd=i;}}return PTrackInd;})
Definition: NumuVarsIcarus202106.h:13