98 for (
unsigned int i=0;
i<cvnneutrons.size(); ++
i){
99 float val = cvnneutrons[
i]->Value();
100 int id = cvnneutrons[
i]->Pdg();
118 for (
unsigned int i=0;
i<cvnparts.size(); ++
i){
119 float val = cvnparts[
i]->Value();
120 int id = cvnparts[
i]->Pdg();
158 cvnpart.
maxval = largestv;
159 cvnpart.
pdgmax = largestid;
172 const std::vector<rb::Cluster>& sliceList,
173 std::vector<cheat::TrackIDE>& allTracks,
174 std::vector<cheat::TrackIDE>& sliceTracks,
175 std::vector<cheat::TrackIDE>& allTracksBirks,
176 std::vector<cheat::TrackIDE>& sliceTracksBirks,
229 std::map<int, float> tPartLongLL = slid.
fPartLongLL;
231 std::map<int, float>::iterator itr;
234 if (itr != tPartLongLL.end()) shwlid.
lid.
elll = itr->second;
238 if (itr != tPartTransLL.end()) shwlid.
lid.
ellt = itr->second;
242 if (itr != tPartLongLL.end()) shwlid.
lid.
mulll = itr->second;
246 if (itr != tPartTransLL.end()) shwlid.
lid.
mullt = itr->second;
601 std::vector<float>::const_iterator maxIt =
633 std::copy(tval.begin(), tval.end(), srcvnmap.
cvnmap);
635 std::copy(tlab.begin(), tlab.end(), srcvnmap.
cvnlabmap);
637 std::copy(tobj.begin(), tobj.end(), srcvnmap.
cvnobjmap);
657 std::copy(tval.begin(), tval.end(), srcvnmap.
cvnmap);
706 if(
int(srcvnfeatures.
components.size()) > maxcomponents)
707 srcvnfeatures.
components.resize(maxcomponents);
732 cafSliceLID.
ncid = artSliceLID.
nc;
733 cafSliceLID.
nueid = artSliceLID.
nue;
int cosfwdcell
cosmictrack projected # cells from end point forwards to det edge
float ncid
Likelihood Neutral Current.
::xsd::cxx::tree::id< char, ncname > id
unsigned char cvnlabmap[16000]
store 1D array with hit labels
int ncell
number of cells in slice
int ismuon
Boolean to distinguish muons from electrons.
SRShowerLID shwlid
Shower information.
double CosAngleToNextProng() const
Cosine of Angle between the leading and the next leading prong.
float KalThetaTrue() const
float pllt
Proton ll for ltransverse shower.
int FirstCell(const unsigned int &view) const
float annENoCos
ann output, with E without cosTheta
double MuClosestApproachByDist() const
Distance of closest approach between this slice and the possible muon parent (selected by closest app...
float gap
gap from vertex to start of shower
SRParticleTruth truthXView
Truth information for the prong.
float pilll
Pion ll for longitudinal shower.
Contains the SPID pid information for a shower (slid::ShowerLID objects).
unsigned int noutput
Number of entries in output vector.
double NCPIDBDTGLT() const
double ProngMaxY() const
The maximum Y between the start and stop points of the leading prong in the slice.
float fscattsum
sum of scattering angles (Fernanda)
float scatllh
Log-likelihood value from scattering angle.
unsigned int interaction
Class of the event.
unsigned int nplanes
number of planes in pixel map
float EFrac2PlWin() const
unsigned short maxplanecont
maximum number of contiguous planes in prong
float enllt
Electron - neutron LL for transverse shower.
void FillCVNNeutronDaughterResultVars(const std::vector< art::Ptr< rb::PID > > &cvnneutrons, caf::SRCVNNeutronDaughterResult &cvnneutron)
float neutronid
Likelihood neutron.
int labelx
Class of the X view prong.
float shwE
Energy of shower [GeV].
unsigned int NCells() const
Variables describing Michel E's found around the end of a track.
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
unsigned int firstcellx
global position of first x cell in map
float ellt
Electron ll for ltransverse shower.
float anne
ann with energy output
InteractionType fInt
Class of the event.
float kalchisq
chisq value of slope fit (kalmantrack)
float pionid
Likelihood pion.
float epi0llt
Electron - Pi0 ll for ltransverse shower.
float glll
Gamma ll for longitudinal shower.
double MuClosestApproachByTime() const
Distance of closest approach between this slice and the possible muon parent (selected by difference ...
float efrac6plwin
largest fraction of energy in a 6 plane window out of total slice
const art::PtrVector< rb::CellHit > & XCells() const
Get all cells from the x-view.
float muonid
Likelihood muon.
double StopMinDistToBottom() const
Minimum perpendicular distance of all prongs stop points to the bottom edge of the detector...
SRVector3D dir
Shower direction at start point [unit vector recommended].
float dedx0
Plane dedx for e/photon/pi0 identification.
float elll
electron ll for longitudinal shower
float nueid
Likelihood of Charge Current NuE.
std::vector< unsigned char > PixelMapToObjVector(bool useLabels) const
unsigned int ncells
number of cells in pixel map
float shwEFrac
fraction of energy of leading shower out of total energy of slice
float emullt
Electron - muon ll for ltransverse shower.
std::vector< float > components
Principal components of output.
float prongEbal3D
the energy balance between the two highest energy 3D prongs in slice
This is a helper class for ParticleIDAlg that provides a tidy structure in which to hold the dE/dx hi...
float KalYPosAtTrans() const
float meanyexp
Hadronic y of matches, "exp".
float costhetatrue
cosine of angle between cosmic trk dir and true dir of most contributing particle ...
void FillVetoVars(const presel::Veto &veto, const presel::Veto &nueveto, caf::SRVeto &srveto)
float cfitspeed
fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for cosmic track ...
unsigned int lastcelly
global position of last y cell in map
float vtxgev
Energy of slice in vertex region.
double MuAngleDiffByDist() const
Minimum angle difference between this slice and the possible muon parent (selected by closest approac...
float KalBakSteel() const
float cospidbdtg
NC cosrej PID. A Keras using 15 variables.
float mullt
Muon ll for ltransverse shower.
float len
Track length of identified track.
float UnconTunedCosPID() const
float epillt
Electron - pion ll for ltransverse shower.
double MuAngleDiffByTime() const
Minimum angle difference between this slice and the possible muon parent (selected by difference in t...
float KalFwdSteel() const
float cosyposattrans
Y position of cosmic track at transition (ND only, use to check if went through air gap) ...
double StopDistToTop() const
Perpendicular distance of stop point of the leading prong to the top edge of the detector.
void FillMuonIDVars(const muonid::MuonID &muid, SRKalmanTrack &srTrk)
unsigned short nplaney
number of planes spanned in y view
float nllt
Neutron ll for ltransverse shower.
float eplll
Electron - proton ll for longitudinal shower.
float hitfracy
fraction of y hits in slice in y view pixel map
double ScatSeparation() const
Return the scattering separation variable used as an input to the kNN that determines the pid value...
ProngType fLabelY
Class of the Y view prong.
double StopMinDistToWest() const
Minimum perpendicular distance of all prongs stop points to the west edge of the detector (west is po...
float maxval
Maximum value among net outputs.
void CopyMuonIDVars(const SRKalmanTrack &srTrk, caf::SRMuonID &muid)
double StartMinDistToTop() const
Minimum perpendicular distance of all prongs start points.
void FillCVNFeaturesVars(const cvn::Features &features, caf::SRCVNFeatures &srcvnfeatures, int maxcomponents)
float gap
gap from shower start to ElasticArms vertex [cm]
float pid
The default PID value - normally use this one.
unsigned short nhit
number of hits
bool passanglefirstana
Would this slice pass the old angle cut used for first analysis.
std::map< int, float > fPartLongLL
Map of the longitudinal ll by paricle type.
Vertical planes which measure X.
float tmvabdtmf
TMVA boosted decision tree, Fisher discriminant.
float cosslope
slope of track timing fit (cosmictrack)
float muclosestapproachbydist
float kalyposattrans
Y position of Kalman track and transition (ND only, use to check if went through air gap) ...
unsigned int fNPlane
Number of planes, length of pixel map.
float tmvabdt
TMVA boosted decision tree, adaptive boost.
float epllt
Electron - proton ll for ltransverse shower.
unsigned int ncells
number of cells in pixel map
float remavededx
average dE/dx value as computed by remid
int pdgbest
pdg code of most contributing true particle to best kalman trk
ParticlesType fFSlong
Class of the event by primary n prongs.
void FillJMEIDVars(const jmshower::EID &eid, caf::SRJMEid &sreid)
float epilll
Electron - pion ll for longitudinal shower.
double CosPIDDTG() const
Gradient BDT trained, based on simple containment variables.
float eplll
Electron - proton LL for longitudinal shower.
float fscattsig
sigma of scattering variable distribution (Fernanda)
Details of processing cuts made by the veto modules.
float dedx2
Plane dedx for e/photon/pi0 identification.
unsigned int parent
Class of the event by parent particle.
float energydiffexp
Pot. diff between sig and bkg matches, "exp".
unsigned short nhitx
number of hits in x-view
float dedx5
Plane dedx for e/photon/pi0 identification.
ProngType fLabelX
Class of the X view prong.
float ellt
Electron LL for transverse shower.
float meanyfit
Hadronic y of matches, "fit".
float mipfrac
fraction of mip hits out of total number of hits
float cospidbdt
NC Cosrej PID. A real adaptive BDT using 13 variables.
Output from Cosmic Rejection (CosRej) module.
float hadE
Hadronic energy [GeV].
float meanqfracbkg
Mean fraction matched in background matches.
std::vector< float > fOutput
Vector of outputs from neural net.
float enrichfrac
Fraction of matches that are "enriched".
double StopDistToWest() const
Perpendicular distance of stop point of the leading prong to the west edge of the detector (west is p...
double NCPIDBDTG() const
Gradient BDT trained, based on prong/shower variables.
float annepi0
ann for e/pi0 separation
float calE
energy based on summed calibrated deposited charge [GeV]
float vtxdist
shortest distance to wall from primary Kalman track start position (proxy for vertex) ...
unsigned int particles
Class of the event by primary n prongs.
double StopDistToFront() const
Perpendicular distance of stop point of the leading prong to the front edge of the detector...
double HitsPerPlaneAsymmSlice() const
Asymmetry in the hitsperplane at the start and end of the slice.
float annepi0
e/pi0 PID for oscillation analysis
int cosbakcell
cosmictrack projected # cells from start point backwards to det edge
double PhotonShowerTransMom() const
Transverse component of the event momentum. Particle masses in the momentum calculation are assumed t...
bool passpre
Bool saying if event passed Preselection.
virtual TVector3 Start() const
float kalfwdair
for Kalmantrack projected distance forwards how much is through air (ND only, NYI) ...
Features, basic output of CVN neural net.
bool PassAngleFirstAna() const
float CosBakSteel() const
unsigned distance(const T &t1, const T &t2)
float costheta
cosine of track with respect to beam direction
float scatt
sum of all all angular deviation (looping over kaltrk traj pts) / trk len
float chibkgexp
chisq of a fit assuming best match is bkg, "exp"
float otherid
Likelihood other.
float cscorediff
chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit modu...
double CosPIDLight() const
A BDT trained based on simple containment variables and the nue output of CVN. Includes pX/p and pY/p...
int kalbakcellnd
Kalmantrack projected # cells from start point backwards to det edge, including muon catcher...
bool passthru
Would this slice pass the (fixed) through-going check?
float nutauid
Likelihood Charge Current NuTau.
std::vector< float > fOutput
Vector of outputs from neural net.
bool passangle
Would this slice pass the most recent angle cut?
bool passthruold
Would this slice pass the (not fixed, no ratio check) through-going check?
Horizontal planes which measure Y.
void FillSlidVars(const slid::ShowerLID &slid, SRShowerLID &shwlid)
float emulll
Electron - muon LL for longitudinal shower.
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
float nueid
Likelihood Charge Current NuE.
float tmvabdtd
TMVA boosted decision tree, decorrelation.
float pngptp
Event transverse momentum fraction, based on prongs. See ncid::NCCosRej::ProngTransMom.
void FillSliceLID(const SliceLID::Prediction &artSliceLID, caf::SRSliceLID &cafSliceLID)
float invglll
photon LL for longitudinal shower, assuming shower started at opposite end
float kalslope
slope of track timing fit (kalmantrack)
float annNoCos
ann output, without cosTheta
float emullt
Electron - muon LL for transverse shower.
float kalfwdsteel
for Kalmantrack projected distance forwards, how much is through steel (ND only, currently is just di...
float epilll
Electron - pion LL for longitudinal shower.
float sparsenessasymmslice
float meanqfracbkgexp
Fraction matched in background matches, "exp".
float cospidcontain
Nue cosrej PID. A BDT using 5 variables relating to containment.
unsigned int noutput
Number of entries in output vector.
virtual double TotalLength() const
Length (cm) of a shower.
float kalfwddist
Kalmantrack projected distance (cm) from end point forwards to det edge.
float gllt
Gamma ll for ltransverse shower.
float annecos
e/cos PID for oscillation analysis
void CopyRemidVars(const SRKalmanTrack &srTrk, caf::SRRemid &remid)
int LastCell(const unsigned int &view) const
Store BDT variables for the short-baseline oscillation study.
float ellt
electron ll for transverse shower
float dedx3
Plane dedx for e/photon/pi0 identification.
float angle
Cosine of CosmicTrack with respect to the beam direction.
float hadE
GeV sum of non (Kalman best ReMId) track E.
void FillNCPi0BkgRejVars(const ncpi0::NCPi0BkgRej &ncpi0bkgrej, caf::SRNCPi0BkgRej &srncpi0bkgrej)
float cosmicid
Likelihood Cosmic Interaction.
double fVtxY
True vertex of Y position.
float mindist
minimum projected distance to edge of any Kalman track with > 15 hits
float epl20frac
fraction of energy in first 20 planes out of total slice energy
A class detailing the cuts made on a particular slice.
unsigned int lastplane
global position of last plane in map
float pid
The output value from the BDT.
int cont
True if track is contained as defined by ReMId.
float elll
Electron LL for longitudinal shower.
double CosPIDDT() const
Real Adaptive BDT trained, based on simple containment variables.
float photonid
Likelihood photon.
float cosbakdist
cosmictrack projected distance (cm) from start point backwards to det edge
void FillCVNProngTrainingData(const cvn::ProngTrainingData &cvnpdata, caf::SRProngTrainingData &srpdata)
SRParticleTruth FillParticleTruth(const std::vector< rb::Cluster > &sliceList, const art::PtrVector< rb::CellHit > &hits, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, int sliceIdx)
Function to fill the particle truth for a set of hits.
std::vector< float > output
List of net output nodes.
int MostMissingPlanes(geo::View_t view) const
Longest run of adjacent planes with no hits.
double StopMinDistToTop() const
Minimum perpendicular distance of all prongs stop points to the top edge of the detector.
float cospidperibdt
Nue cosrej PID for the peripheral sample for 2020+.
double vtxx
True vertex of X position.
float electronid
Likelihood electron.
float ProngEBal3D() const
float pid
PID value output by kNN.
double AngleVarOld() const
int kalfwdcellnd
Kalmantrack projected # cells from end point forwards to det edge, including muon catcher...
float pi0mass
best pi0 mass hypothesis coming from combinations of JMShowers
bool keepnue
Does the event past nue version of veto? (see doc-db 14654)
void FillLEMVars(const lem::PIDDetails &lem, caf::SRLem &srlem)
void FillNuePreselVars(const presel::PreselObj &nuepre, caf::SRPresel &srnuepre)
void FillCosRejVars(const cosrej::CosRejObj &cosrej, caf::SRCosRej &srcosrej, caf::SRContain &srcontain)
float fLepEnergy
True energy of outgoing lepton.
void FillXnueVars(const xnue::Xnue &xnue, caf::SRXnue &srxnue)
int labely
Class of the Y view prong.
double StartMinDistToEast() const
Minimum perpendicular distance of all prongs start points to the east edge of the detector (east is n...
double StopMinDistToEast() const
Minimum perpendicular distance of all prongs stop points to the east edge of the detector (east is ne...
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
float fNuEnergy
True energy of neutrino event.
float meanqfracfit
Fraction of charge matched, "fit".
double StartMinDistToBottom() const
Minimum perpendicular distance of all prongs start points to the bottom edge of the detector...
float anglekal
cos of angle of best ReMId Kalman track
unsigned int nplanes
number of planes in pixel map
float epi0lll
Electron - Pi0 LL for longitudinal shower.
int mincell
minimum projected cells to edge of any Kalman track with > 15 hits
int ismuon
Boolean to distinguish muons from electrons.
std::vector< unsigned char > PixelMapToVector(bool useGeV) const
void FillNueCosRejVars(const cosrej::NueCosRej &nuecosrej, caf::SRNueCosRej &srnuecosrej)
double StartMinDistToBack() const
Minimum perpendicular distance of all prongs start points to the back edge of the detector...
float annEPi0EL
ann output, EPi0 Elastic scattering
float meanysig
Mean hadronic y of matches that are signal.
unsigned char cvnobjmap[16000]
store 1D array with hit object id
float bakdist
Projected live distance to the edge projected backwards from track start.
float protonid
Likelihood proton.
void FillRVPVars(const rvp::RVP &rvp, caf::SRRvp &srrvp)
float kdirscore
chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for best ...
unsigned int firstplane
global position of first plane in map
const art::PtrVector< rb::CellHit > & YCells() const
Get all cells from the x-view.
float fannid
ANN output, "fann" package.
float meanqfracsig
Fraction matched in signal matches.
unsigned int pdgmax
pdg of largest likelihood
double ParticleShowerTransMom() const
Transverse component of the event momentum. Particle masses in the momentum calculation are inferred ...
float cosmicid
Likelihood of Cosmic.
virtual TVector3 Dir() const
Unit vector describing prong direction.
TVector3 Stop() const
Endpoint of the shower.
float mulll
Muon LL for longitudinal shower.
unsigned int ncomponents
Size of components vector.
double CosPIDContain() const
A BDT trained based on simple containment variables Used in the peripheral sample of the 2017 analysi...
float cospidcorebdt
Nue cosrej PID for the core sample for 2020+.
float elll
Electron ll for longitudinal shower.
std::vector< float > output
List of net output nodes.
float epllt
Electron - proton LL for transverse shower.
float cosfwddist
cosmictrack projected distance (cm) from end point forwards to det edge
float measfrac
Fraction of planes used to measure dE/dx.
int cosfwdcellnd
cosmictrack projected # cells from end point forwards to det edge, including muon catcher...
CVN PID output for a single particle.
float HadronicEnergy() const
float depE
Total energy deposited in shower [GeV].
float pidexp
Fraction of matches that are signal, "exp".
float EFrac6PlWin() const
float enlll
Electron - neutron LL for longitudinal shower.
float CosThetaTrue() const
double HitsPerPlaneAsymm() const
Asymmetry in the hitsperplane at the start and end of the leading shower.
float energydiff
Potential diff between sig and bkg matches.
float meanqfracsigexp
Fraction matched in signal matches "exp".
void FillCVNParticleResultVars(const std::vector< art::Ptr< rb::PID > > &cvnparts, caf::SRCVNParticleResult &cvnpart)
float epillt
Electron - pion LL for transverse shower.
unsigned int lastcellx
global position of last x cell in map
float cosfwdair
for cosmictrack projected distance forwards how much is through air (ND only, NYI) ...
float fFracY
fraction of y view hits contained in map
unsigned int nchan
number of channels in pixel map
Contains the reco muon PID (ReMId) output.
std::map< int, float > fPartTransLL
Map of the transverse ll by particle type.
double vtxz
True vertex of Z position.
double ProngMinY() const
The minimium Y between the start and stop points of the leading prong in the slice.
unsigned short nplanex
number of planes spanned in x view
double CosPIDContainXY() const
A BDT trained based on simple containment variables. Includes pX/p and pY/p separately instead of pT/...
float plll
Proton ll for longitudinal shower.
float cosfwdsteel
for cosmictrack projected distance forwards, how much is through steel (ND only, currently is just di...
This class contains the LEM PID output.
float annE
ann output, with E
float avgsigE
Fraction of signal matches, sigmoid weighting.
float rawgevslice
event energy in GeV
double StartMinDistToWest() const
Minimum perpendicular distance of all prongs start points to the west edge of the detector (west is p...
float chisig
chi-square of a fit assuming best match is signal
double StartDistToBottom() const
Perpendicular distance of start point of the leading prong to the bottom edge of the detector...
float shwEFrac
fraction of energy of leading shower out of total energy of slice
double StopMinDistToFront() const
Minimum perpendicular distance of all prongs stop points to the front edge of the detector...
float efrac2plwin
largest fraction of energy in a 2 plane window out of total slice
SRVector3D start
Shower start point in detector coordinates. [cm].
unsigned int fNCell
Number of cells, width of pixel map.
double StopDistToBack() const
Perpendicular distance of stop point of the leading prong to the back edge of the detector...
float fwddist
Projected live distance to the edge projected forwards from track end.
float meanysigexp
Hadronic y of matches that are signal, "exp".
double ProngDistToVtx() const
Perpendicular distance of closes approach between the leading prong and the event vertex...
float nlll
Neutron ll for longitudinal shower.
unsigned char cvnmap[16000]
store 1D array the size of the pixel map, currently 80 cells x 100 planes x 2 views. Although sparse, store this way for easier conversion and compression in hdf5. First element is channel 0, plane 0, cell 0, then channel 0, plane 0, cell 1, etc. Channel 0 is x view, 1 is y view
float protonid
Likelihood proton.
float dedx4
Plane dedx for e/photon/pi0 identification.
float kalbaksteel
for Kalmantrack projected distance backwards, how much is through steel (ND only, currently is just d...
Output from Cosmic Rejection (Nuecosrej) module.
double PhotonShowerMomY() const
The cosine between the reconstructed momentum and the y-axis. Particle masses in the momentum calcula...
float hitfracx
faction of x hits in slice in x view pixel map
double ProngMinZ() const
The minimium Z between the start and stop points of the leading prong in the slice.
float enrichfracexp
Fraction of matches that are "enriched", "exp".
double StartMinDistToFront() const
Minimum perpendicular distance of all prongs start points to the front edge of the detector...
CVN PID output for a single particle.
float emulll
Electron - muon ll for longitudinal shower.
float CosFwdSteel() const
The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes i...
float kalchidiff
chisq difference of track timing fits (kalmantrack)
float pi0lll
Pi0 ll for longitudinal shower.
float hitsperplaneasymmslice
void FillNCCosRejVars(const ncid::NCCosRej &nccosrej, caf::SRNCCosRej &srnccosrej)
float avgexpE
Fraction of signal matches, "exp" weighting.
void FillShowerVars(const rb::Shower &shw, caf::SRFuzzyKProng &srshw, const std::vector< rb::Cluster > &sliceList, std::vector< cheat::TrackIDE > &allTracks, std::vector< cheat::TrackIDE > &sliceTracks, std::vector< cheat::TrackIDE > &allTracksBirks, std::vector< cheat::TrackIDE > &sliceTracksBirks, const int &sliceIdx)
Contains the RVP PID output.
int ncell
Number of hits in the slice.
bool passnumicut
Does this event pass the NuMI timing cuts?
float miphits
number of mip hits in slice
float pi0id
Likelihood primary prong is a pi0.
float coschidiff
chisq difference of track timing fits (cosmictrack)
unsigned short nplane
number of planes spanned
unsigned int NYCell() const
Number of cells in the y-view.
float coschisq
chisq value of slope fit (cosmictrack)
float longtrfrac
fraction of cells in longest track out of total cells in slice
float emid
Likelihood electron or photon.
double vtxy
True vertex of Y position.
void FillMuIdVars(const remid::ReMId &remid, SRKalmanTrack &srTrk)
unsigned short nhity
number of hits in y-view
float diry
Cosine of CosmicTrack with respect to the y-direction.
float enllt
Electron - neutron ll for ltransverse shower.
unsigned int lastplane
global position of last plane in map
double fVtxX
True vertex of X position.
SRParticleTruth truth
Truth information for the prong.
void FillLIDEnergyVars(const slid::ShowerLID &slid, caf::SRFuzzyKProng &png)
double StartDistToTop() const
Perpendicular distance of start point of the leading prong to the top edge of the detector...
Cosmic Rejection PIDs for Numu analysis.
float dedxllh
Log-likelihood value from dE/dx.
float meanqfracexp
Fraction of charge matched, "exp".
virtual void setDefault()
unsigned int firstcelly
global position of first y cell in map
float epi0nuoneid
Likelihood nu-on-e from epi0 classifier.
float len
track length [cm]
Output of the jmshower::NueSel module.
float fannidenrich
ANN output, including enrichfracexp.
float maxval
Maximum likelihood among net outputs.
float eglll
Electron - gamma LL for longitudinal shower.
ProngType fLabel3D
Class of the 3D prong.
float muclosestapproachbytime
float eiso3sig
fraction of slice energy more then 3sigma from mean energy weighted slice position ...
float cosbaksteel
for cosmictrack projected distance backwards, how much is through steel (ND only, currently is just d...
Result, basic output of CVN neural net.
float enlll
Electron - neutron ll for longitudinal shower.
float prongEbal2D
the energy balance between highest energy 2D prong and 3D prong in the opposite view ...
float ShowerDepEnergy() const
int kalfwdcell
Kalmantrack projected # cells from end point forwards to det edge.
int label3d
Class of the 3D prong.
float VertexEnergy() const
float anglevar
Defined as: abs(angle)*(diry+1) - new angle variable for new angle cut.
float pidexpfit
Fraction of matches that are signal, "expfit".
float egllt
Electron - gamma LL for transverse shower.
float fFracX
fraction of x view hits contained in map
FinalStateProngsType fFSp
Class of the event by primary n prongs.
float enddist
shortest distance to wall from primary Kalman track end position
float ShowerEnergy() const
float muonid
MuonID classifier value.
double ProngMaxX() const
The maximum X between the start and stop points of the leading prong in the slice.
double ProngTransMom() const
Transverse component of the energy-weighted average of prong directions.
float pi0mass
best pi0 mass hypothesis coming from combinations of JMShowers
Horizontal planes which measure Y.
float meany
Mean hadronic y of matches.
float annepi0el
particle PID.
float pid
MuonID value output by ND-trained BDT.
float epi0pi0id
Likelihood pi0 from epi0 classifier.
float otherid
Likelihood other.
float mulll
Muon ll for longitudinal shower.
float pi0llt
Pi0 ll for ltransverse shower.
float pizeroid
Likelihood pizero.
float dedx1
Plane dedx for e/photon/pi0 identification.
float CosYPosAtTrans() const
int nkal3d
number of 3D tracks as determined by Kalman tracker
double SparsenessAsymm() const
Asymmetry in the zero-hit planes at the start and end of the leading shower.
SRParticleTruth truthYView
Truth information for the prong.
The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes i...
A rb::Prong with a length.
unsigned int nchan
number of channels in pixel map
unsigned int NXCell() const
Number of cells in the x-view.
float eratio
ratio of best kalman track GeV / slice GeV
PixelMap, basic input to CVN neural net.
float recoE
reco energy of slice
int MostContiguousPlanes(geo::View_t view) const
Longest run of adjacent planes with hits.
float rvp12
PID value for 12 variable RVP for comparison purposes.
std::vector< float > output
List of net output nodes.
double StartDistToEast() const
Perpendicular distance of start point of the leading prong to the east edge of the detector (east is ...
float numuid
Likelihood Charge Current NuMu.
Contains the reco muon PID (ReMId) output.
float ValueEPi0EL() const
Variables describing Michel E's found around the end of a track.
float ProngEBal2D() const
unsigned int ExtentPlane(geo::View_t view=geo::kXorY) const
double PhotonShowerMomX() const
The cosine between the reconstructed momentum and the x-axis. Particle masses in the momentum calcula...
unsigned char cvnmap[16000]
store 1D array the size of the pixel map, currently 80 cells x 100 planes x 2 views. Although sparse, store this way for easier conversion and compression in hdf5. First element is channel 0, plane 0, cell 0, then channel 0, plane 0, cell 1, etc. Channel 0 is x view, 1 is y view
float hadronid
Likelihood Proton or Pion.
Vertical planes which measure X.
unsigned int NCell() const
float InverseGLLL() const
float epi0llt
Electron - Pi0 LL for transverse shower.
double StopDistToBottom() const
Perpendicular distance of stop point of the leading prong to the bottom edge of the detector...
float nuenergy
True energy of neutrino event.
int kalbakcell
Kalmantrack projected # cells from start point backwards to det edge.
float kalthetatrue
cosine of angle between kalman trk dir and true dir of most contributing particle ...
float nueccid
Likelihood nuecc events.
double ProngTransMom() const
Transverse component of the energy-weighted average of prong directions.
float EFrac2SigRd() const
unsigned int firstcellx
global position of first x cell in map
double fVtxZ
True vertex of Z position.
unsigned int lastcelly
global position of last y cell in map
int prongs3D
number of 3D FuzzyK prongs associated with slice
Attach LEM-specific info to the base PID object.
unsigned int lastcellx
global position of last x cell in map
float numuunconttunedpid
tuned cosmic rejection PID for uncontained events - Jose's;
double MeasurementFraction() const
Return the measurement fraction variable used as an input to the kNN that determines the pid value...
float tmvamlp
TMVA ANN ("multi layer perceptron")
double HitsPerPlane() const
Number of hits in slice per number of planes in leading prong. This is a handle on rejecting FEB flas...
float gap
gap from vertex to start of shower
int MuSliceIdxByTime() const
The index of the slice that is most likely to be the parent muon according to difference in time...
double StartDistToBack() const
Perpendicular distance of start point of the leading prong to the back edge of the detector...
void FillRockPreselVars(const presel::PreselObj &rockpre, caf::SRPresel &srrockpre)
void FillSpidVars(const slid::ShowerPID &spid, caf::SRSPid &srspid)
Build slid::LID objects to store electron ID, if asked for, otherwise, calculate LID info and make av...
SRVector3D stop
shower stop point
float mullt
Muon LL for transverse shower.
double ProngMaxZ() const
The maximum Z between the start and stop points of the leading prong in the slice.
void FillCVNPixelMaps(const cvn::PixelMap &map, caf::SRPixelObjMap &srcvnmap, bool useGeV)
float epi0lll
Electron - Pi0 ll for longitudinal shower.
unsigned int finalstate
Class of the event by final state.
float eelllt
electron el ll for transverse shower
void FillCVNTrainingData(const cvn::TrainingData &cvntdata, caf::SRTrainingData &srtdata)
Output from Cosmic Rejection (Nuecosrej) module.
float longtr
longest discrete track merge track
int MuSliceIdxByDist() const
The index of the slice that is most likely to be the parent muon according to distance of closest app...
This module creates Common Analysis Files.
double StartDistToFront() const
Perpendicular distance of start point of the leading prong to the front edge of the detector...
ParentParticleType fPPt
Class of the event by parent particle.
float nutauid
Likelihood of Charge Current NuTau.
float efrac2sigrd
fraction of slice energy within 2sigma of mean energy weighted slice position
float ncid
Likelihood of Neutral Current.
float numucontpid2019
cosmic rejection PID for contained events; 2019 Analysis
void FillCVNResultVars(const cvn::Result &result, caf::SRCVNResult &srcvn, bool NoCosmics)
double StopMinDistToBack() const
Minimum perpendicular distance of all prongs stop points to the back edge of the detector.
T min(const caf::Proxy< T > &a, T b)
float kalbakdist
Kalmantrack projected distance (cm) from start point backwards to det edge.
float numuid
Likelihood of Charge Current NuMu.
void FillNuonEResultVars(const cvn::Result &result, caf::SRNuonEResult &srnuone)
float meanqfrac
Mean fraction of charge matched.
float cosbakair
for cosmictrack projected distance backwards how much is through air (ND only, NYI) ...
double MuTimeDiffByDist() const
Minimum time difference between this slice and the possible muon parent (selected by closest approach...
double MuTimeDiffByTime() const
Minimum time difference between this slice and the possible muon parent (selected by difference in ti...
Object collecting Preselection variables.
float numucontpid2020
cosmic rejection PID for contained events; 2020 Analysis
float kfitspeed
fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for best ReMId Kalman track ...
double StopDistToEast() const
Perpendicular distance of stop point of the leading prong to the east edge of the detector (east is n...
float ratio
Number of hits in the track relative to the slice.
float anglecos
cos of angle of Cosmic Track
unsigned int argmax
Index of maximum value.
float photonid
Likelihood photon.
float avginvE
Fraction of signal matches, 1/E weighting.
float len
track length [cm]
unsigned short maxplanegap
maximum number of gapped planes in prong
float chisigexp
chisq of a fit assuming best match is signal, "exp"
float pidfit
Fraction of matches that are signal, "fit".
int cosbakcellnd
cosmictrack projected # cells from start point backwards to det edge, including muon catcher...
float chibkg
chisq of a fit assuming best match is background
int prongs2D
number of 2D FuzzyK prongs associated with slice
double SparsenessAsymmSlice() const
Asymmetry in the zero-hit planes at the start and end of the slice.
float egllt
Electron - gamma ll for ltransverse shower.
unsigned int noutput
Number of entries in output vector.
float width
Shower width [cm].
View_t view
Prong view caf::kX = 0, caf::kY = 1 or caf::kXorY = 2.
unsigned int firstcelly
global position of first y cell in map
float nuoneid
Likelihood nu-on-e.
float mullt
muon ll for transverse shower
float kscorediff
chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit modu...
float fscattext
scattering variable (Fernanda)
double DedxSeparation() const
Return the dE/dx separation variable used as an input to the kNN that determines the pid value...
float pillt
Pion ll for ltransverse shower.
float dectree
Decision tree PID.
float annEPi0
ann output, EPi0
double ProngMinX() const
The minimium X between the start and stop points of the leading prong in the slice.
float cdirscore
chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for cosmi...
float vtxgev
Energy of slice in vertex region.
float hitfracy
fraction of y hits in slice in y view pixel map
FinalStateType fFSt
Class of the event by final state.
float fscattmax
maximum scattering angle (Fernanda)
float otherid
Likelihood other.
unsigned int finalstateprong
float ShowerEFrac() const
std::vector< float > fPrincipalComponents
Linear combinations of fOutput.
float anglevarold
Difference of squares between beam angle and Y-direction (Also labeled CosmicVar AnglePID and AngleVa...
float mulll
muon ll for longitudinal shower
double StartDistToWest() const
Perpendicular distance of start point of the leading prong to the west edge of the detector (west is ...
float eellll
electron el ll for longitudinal shower
float lepenergy
True energy of outgoing lepton.
float hitfracx
faction of x hits in slice in x view pixel map
float tmvabdtg
TMVA boosted decision tree, gradient boost.
unsigned int firstplane
global position of first plane in map
float kalbakair
for Kalmantrack projected distance backwards how much is through air (ND only, NYI) ...