74 double CVNem_CalE = 0.0;
76 for(
unsigned int fuzzkProngIdx = 0; fuzzkProngIdx < srVertexBranch.
elastic.
fuzzyk.
png.size(); ++fuzzkProngIdx){
78 double png_CalE = srVertexBranch.
elastic.
fuzzyk.
png[fuzzkProngIdx].shwlid.calE;
83 if ( emPID < 0 )
continue;
84 if ( emPID >= hadPID ) CVNem_CalE += png_CalE;
87 if( CVNem_CalE == 0.0) CVNem_CalE = srVertexBranch.
elastic.
fuzzyk.
png[0].shwlid.calE;
107 const double wx = (hadclust.
NXCell() > 0) ? hadclust.
W(hadclust.
XCell(0).
get()) : 0;
108 const double wy = (hadclust.
NYCell() > 0) ? hadclust.
W(hadclust.
YCell(0).
get()) : 0;
112 unsigned int cellsFromEdge = -1;
118 unsigned int minCellsFromEdge = 99999999;
121 for(
unsigned int hitIdx = 0; hitIdx < hadclust.
NCell(); ++hitIdx)
128 const int planeNum = chit->
Plane();
131 if(cellsFromEdge < minCellsFromEdge)
133 minCellsFromEdge = cellsFromEdge;
149 if ( (rhit.
PECorr()>100.0) && (rhit.
PECorr()<245.0) ) nMipHit++;
156 if (!(hadclust.
NCell()==0)){
166 TVector3 boxmin, boxmax, meanpos;
168 srhadclust.
boxmin = boxmin;
169 srhadclust.
boxmax = boxmax;
173 if (hadclust.
NXCell()==0){
177 if (hadclust.
NYCell()==0){
190 srhadclust.
calE = -5;
192 srhadclust.
boxmin = TVector3(-5555,-5555,-5555);
193 srhadclust.
boxmax = TVector3(-5555,-5555,-5555);
194 srhadclust.
meanpos = TVector3(-5555,-5555,-5555);
T max(const caf::Proxy< T > &a, T b)
float angleerror
1-sigma error from using quasielastic angle formula for neutrino energy
virtual double W(const rb::CellHit *chit) const
Estimate the unmeasured coordinate of chit.
float neutronid
Likelihood neutron.
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
float rawgevhad
Uncorrected hadronic energy [GeV].
float NDTrkLenCat() const
float pionid
Likelihood pion.
float muonid
Likelihood muon.
float orphCalE
calorimetric energy of hits that don't appear in any FuzzyK prongs
art::Ptr< rb::CellHit > XCell(unsigned int xIdx) const
Get the ith cell in the x-view.
float ndhadcalactE
Near detector – hadronic calorimetric energy NOT on the muon track in active region [GeV]...
float ndtrklenact
Near detector – muon track length in active region [cm].
unsigned int ncellsfromedge
minimum number of cells to edge of detector
float ndtrkcaltranE
Near detector – muon calorimetric energy in transition plane [GeV].
SRHadClust hadclust
Cluster information for hits in slice but not on muon track (track with higest ReMId value) ...
std::vector< SRProng > png2d
Vector of 2D prong objects.
float trknonqeE
Track length non-quasielastic neutrino energy [GeV].
float rawgevslice
Uncorrected slice energy [GeV].
unsigned short Plane() const
std::vector< SRFuzzyKProng > png
Vector of 3D prong objects.
Overarching information for a numu hadronic cluster.
float ndtrkcalcatE
Near detector – muon calorimetric energy in muon catcher [GeV].
float calccE
Calorimetric charged current neutrino energy [GeV].
Vertical planes which measure X.
unsigned int Ncells() const
Number of cells in this plane.
float NDTrkLenAct() const
float truemuoncatcherE
ND value: true energy of muon as it crosses the transition plane into muon catcher [GeV]...
rb::RecoHit MakeRecoHit(rb::CellHit const &cellhit, double w)
unsigned int ncontplanes
number of continuous planes
A collection of associated CellHits.
float NDHadTrkTran() const
void FillHadClustVars(const rb::Cluster &hadclust, caf::SRHadClust &srhadclust)
unsigned int MaxCell(geo::View_t view) const
float ucrecomuonE1trk
Reconstructed muon energy for uncontained single track events at the FarDet from TMVA [GeV]...
SRVector3D meanpos
Mean position of hits in cluster, weighted by charge [cm].
const PlaneGeo * Plane(unsigned int i) const
unsigned int nhit
number of hits
float NDTrkCalAct() const
rb::Cluster HadCluster() const
float RecoTrkCCHadE() const
float ndhadtrkactE
Near detector – hadronic calorimetric energy on the muon track in active region [GeV].
Calibrated quantities relying on position in the orthogonal view. To generate a rb::CellHit from a rb...
unsigned short Cell() const
float photonid
Likelihood photon.
float UCMuonESingle() const
float electronid
Likelihood electron.
float trkccE
Track length cc neutrino energy [GeV].
SRCVNParticleResult cvnpart
CVN prongID information for 4 views.
float ndhadcalcatE
Near detector – hadronic calorimetric energy NOT on the muon track in muon catcher [GeV]...
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
float trkqeE
Track length quasielastic neutrino energy [GeV].
float hadtrkE
Hadronic calorimetric energy on the muon track[GeV].
Vectors of reconstructed vertices found by various algorithms.
art::Ptr< rb::CellHit > YCell(unsigned int yIdx) const
Get the ith cell in the y-view.
void FillNumuLSTMEnergyVars(const LSTME::LSTMEnergy &artE, float &lstmmuon, float &lstmnu)
float ndtrktranx
Near detector – if muon track crosses transition plane, the x location [cm].
Numu energy estimator output.
unsigned int nmiphit
number of minimum ionizing hits
bool IsCalibrated() const
You MUST check here before accessing PECorr, MIP or GeV.
void FillNumuEnergyVars(const numue::NumuE &e, caf::SRNumuEnergy &sre)
float ucrecomuonE
Reconstructed muon energy for uncontained non-single track events at the FarDet from TMVA [GeV]...
float angleE
Quasielastic angle formula neutrino energy [GeV].
SRTrueNumuEnergy mc
True information used for retuning numu energy fits.
float E
Neutrino energy, set to match trkccE [GeV].
float protonid
Likelihood proton.
float UCMuonENonSingle() const
float NDHadCalTran() const
unsigned int NYCell() const
Number of cells in the y-view.
float calE
Calorimetric energy of the cluster [GeV].
unsigned int MinPlane(geo::View_t view=geo::kXorY) const
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
float ndhadcaltranE
Near detector – hadronic calorimetric energy NOT on the muon track in transition plane [GeV]...
float ndtrkcalactE
Near detector – muon calorimetric energy in active region [GeV].
float AngleQEError() const
float rawgevshw
Uncorrected leading shower energy [GeV].
float ndtrklencat
Near detector – muon track length in muon catcher [cm].
float otherid
Likelihood other.
float pizeroid
Likelihood pizero.
unsigned int NXCell() const
Number of cells in the x-view.
unsigned int MinCell(geo::View_t view) const
SRElastic elastic
Single vertex found by Elastic Arms.
int MostContiguousPlanes(geo::View_t view) const
Longest run of adjacent planes with hits.
float NDHadCalAct() const
unsigned int lastcell
last cell
unsigned int ncalhit
number of hits with calibration
float MCTrueMuonE() const
double TotalGeV(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple sum of the estimated GeV of all the hits.
float ndtrktrany
Near detector – if muon track crosses transition plane, the y location [cm].
unsigned int MaxPlane(geo::View_t view=geo::kXorY) const
SRFuzzyK fuzzyk
Primary 3D prong object.
float ndhadtrktranE
Near detector – hadronic calorimetric energy on the muon track in transition plane [GeV]...
float recotrkcchadE
Reconstructed hadronic energy for track cc neutrino energy estimator [GeV].
void MinMaxMeanXYZ(TVector3 &lo, TVector3 &hi, TVector3 &mean, rb::AveragingScheme scheme=kDefaultScheme) const
Gets the min/max/mean all at once, called by the functions above.
This module creates Common Analysis Files.
SRVector3D boxmin
Minimum coordinates box containing all the hits [cm].
T min(const caf::Proxy< T > &a, T b)
bool MCGoodTrueMuon() const
bool truegoodmuon
Have a primary muon that contributes to at least three hits in each view.
float ndhadtrkcatE
Near detector – hadronic calorimetric energy on the muon track in muon catcher [GeV].
float NDTrkCalTran() const
unsigned int firstcell
first cell
float NDHadTrkCat() const
Nue energy estimator output in the standard record.
float MCTrueMuonCatcherE() const
float NDHadCalCat() const
Encapsulate the geometry of one entire detector (near, far, ndos)
float hadcalE
Hadronic calorimetric energy NOT on the muon track[GeV].
unsigned int firstplane
first plane
void FillNueEnergyVars(const rb::Cluster &slice, const SRVertexBranch &srVertexBranch, caf::SRNueEnergy &sre)
float truemuonE
True energy of primary muon [GeV].
unsigned int lastplane
last plane
float recomuonE
Reconstructed muon energy for all neutrino energy estimators [GeV].
SRVector3D boxmax
Maximum coordinates box containing all the hits [cm].
float NDHadTrkAct() const
float NDTrkCalCat() const