FillPIDs.cxx
Go to the documentation of this file.
1 #include "CAFMaker/FillPIDs.h"
2 
3 #include "CAFMaker/FillTruth.h" // for FillParticleTruth
4 
6 
7 #include "RecoBase/Shower.h"
8 
9 #include "CosRej/CosRejObj.h"
10 #include "CosRej/NueCosRej.h"
11 #include "NCID/NCCosRej.h"
13 #include "RecoJMShower/EID.h"
14 #include "RecoJMShower/JMShower.h"
15 #include "RecVarPID/RVP.h"
16 #include "XnuePID/Xnue.h"
17 #include "ReMId/ReMId.h"
18 #include "MuonID/MuonID.h"
19 #include "LEM/PIDDetails.h"
20 #include "ShowerLID/EventLID.h"
21 #include "ShowerLID/ShowerLID.h"
22 #include "ShowerLID/ShowerPID.h"
25 #include "Preselection/PreselObj.h"
26 #include "Preselection/Veto.h"
28 #include "CVN/func/Result.h"
29 #include "CVN/func/CVNFeatures.h"
30 #include "CVN/func/PixelMap.h"
31 #include "CVN/func/TrainingData.h"
33 
34 namespace caf
35 {
36  //......................................................................
38  {
39  srTrk.remavededx = remid.AvgDedx();
40  srTrk.rempid = remid.Value();
41  srTrk.scatllh = remid.ScatSeparation();
42  srTrk.dedxllh = remid.DedxSeparation();
43  srTrk.measfrac = remid.MeasurementFraction();
44  srTrk.remcont = remid.Contained();
45  }
46 
47  //......................................................................
48  void FillMuonIDVars(const muonid::MuonID& muid, SRKalmanTrack& srTrk)
49  {
50  srTrk.muonid = muid.Value();
51  }
52 
53  //......................................................................
54  void FillJMEIDVars(const jmshower::EID& eid,
55  caf::SRJMEid& sreid)
56  {
57  sreid.ann = eid.ANN();
58  sreid.annE = eid.ANNE();
59  sreid.annNoCos = eid.ANNNoCos();
60  sreid.annENoCos = eid.ANNENoCos();
61  sreid.annEL = eid.ANNEL();
62  sreid.annEPi0 = eid.ANNEPi0();
63  sreid.annEPi0EL = eid.ANNEPi0EL();
64  sreid.ismuon = eid.IsMuon();
65  sreid.dedx0 = eid.Dedx0();
66  sreid.dedx1 = eid.Dedx1();
67  sreid.dedx2 = eid.Dedx2();
68  sreid.dedx3 = eid.Dedx3();
69  sreid.dedx4 = eid.Dedx4();
70  sreid.dedx5 = eid.Dedx5();
71  sreid.eglll = eid.EGLLL();
72  sreid.egllt = eid.EGLLT();
73  sreid.emulll = eid.EMuLLL();
74  sreid.emullt = eid.EMuLLT();
75  sreid.epi0lll = eid.EPi0LLL();
76  sreid.epi0llt = eid.EPi0LLT();
77  sreid.eplll = eid.EPLLL();
78  sreid.epllt = eid.EPLLT();
79  sreid.enlll = eid.ENLLL();
80  sreid.enllt = eid.ENLLT();
81  sreid.epilll = eid.EPiLLL();
82  sreid.epillt = eid.EPiLLT();
83  sreid.vtxgev = eid.VtxGeV();
84  sreid.pi0mass = eid.Pi0Mass();
85  sreid.shwEFrac = eid.ShE();
86  sreid.gap = eid.Gap();
87  sreid.elll = eid.ELLL();
88  sreid.ellt = eid.ELLT();
89  sreid.eellll = eid.EelLLL();
90  sreid.eelllt = eid.EelLLT();
91  sreid.mulll = eid.MuLLL();
92  sreid.mullt = eid.MuLLT();
93  }
94  //......................................................................
95  void FillCVNNeutronDaughterResultVars(const std::vector<art::Ptr<rb::PID> >& cvnneutrons,
97  {
98  for (unsigned int i=0; i<cvnneutrons.size(); ++i){
99  float val = cvnneutrons[i]->Value();
100  int id = cvnneutrons[i]->Pdg();
101  if(id==2212)
102  cvnneutron.protonid = val;
103  else if(id==22)
104  cvnneutron.photonid = val;
105  else
106  cvnneutron.otherid = val;
107  }
108 
109  }
110 
111  //......................................................................
112  void FillCVNParticleResultVars(const std::vector<art::Ptr<rb::PID> >& cvnparts,
113  caf::SRCVNParticleResult& cvnpart)
114  {
115 
116  float largestv = 0;
117  int largestid = 0;
118  for (unsigned int i=0; i<cvnparts.size(); ++i){
119  float val = cvnparts[i]->Value();
120  int id = cvnparts[i]->Pdg();
121  if (val > largestv){
122  largestv = val;
123  largestid = id;
124  }
125  switch (id) {
126  case 11:
127  cvnpart.electronid = val;
128  break;
129  case 13:
130  cvnpart.muonid = val;
131  break;
132  case 2212:
133  cvnpart.protonid = val;
134  break;
135  case 2112:
136  cvnpart.neutronid = val;
137  break;
138  case 211:
139  cvnpart.pionid = val;
140  break;
141  case 111:
142  cvnpart.pizeroid = val;
143  break;
144  case 22:
145  cvnpart.photonid = val;
146  break;
147  case 22 * 11:
148  cvnpart.emid = val;
149  break;
150  case 2212 * 211:
151  cvnpart.hadronid = val;
152  break;
153  default:
154  cvnpart.otherid = val;
155  break;
156  }
157  }
158  cvnpart.maxval = largestv;
159  cvnpart.pdgmax = largestid;
160 
161  if(cvnpart.emid < 0)
162  cvnpart.emid = cvnpart.electronid + cvnpart.photonid;
163  if(cvnpart.hadronid < 0)
164  cvnpart.hadronid = cvnpart.pionid + cvnpart.protonid;
165 
166  }
167 
168  //......................................................................
169  void FillShowerVars(const rb::Shower& shw,
170  caf::SRFuzzyKProng& srshw,
171  // const art::PtrVector<rb::Cluster>& sliceList,
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,
177  const int& sliceIdx)
178  {
179  srshw.shwlid.nhit = shw.NCell();
180  srshw.shwlid.nhitx = shw.NXCell();
181  srshw.shwlid.nhity = shw.NYCell();
182  srshw.shwlid.nplane = shw.ExtentPlane();
183  srshw.shwlid.nplanex = shw.ExtentPlane(geo::kX);
184  srshw.shwlid.nplaney = shw.ExtentPlane(geo::kY);
187  srshw.shwlid.len = shw.TotalLength();
188  srshw.shwlid.calE = shw.CalorimetricEnergy();
189  srshw.shwlid.start = shw.Start();
190  srshw.shwlid.dir = shw.Dir();
191  srshw.shwlid.stop = shw.Stop();
192 
193  srshw.truth = FillParticleTruth(sliceList, shw.AllCells(), allTracks, sliceTracks, allTracksBirks, sliceTracksBirks, sliceIdx);
194  srshw.truthXView = FillParticleTruth(sliceList, shw.XCells(),allTracks, sliceTracks, allTracksBirks, sliceTracksBirks, sliceIdx);
195  srshw.truthYView = FillParticleTruth(sliceList, shw.YCells(), allTracks, sliceTracks, allTracksBirks, sliceTracksBirks, sliceIdx);
196  if (shw.NXCell() == 0) srshw.shwlid.view = kY;
197  else if (shw.NYCell() == 0) srshw.shwlid.view = kX;
198  else srshw.shwlid.view = kXorY;
199  }
200 
201  //......................................................................
203  {
204  shwlid.width = slid.Radius();
205  shwlid.gap = slid.Gap();
206  shwlid.lid.ann = slid.Value();
207  shwlid.lid.anne = slid.ValueWithE();
208  shwlid.lid.annepi0 = slid.ValueEPi0();
209  shwlid.lid.annecos = slid.ValueECos();
210  shwlid.lid.ismuon = slid.IsMuon();
211  shwlid.lid.eglll = slid.EGLLL();
212  shwlid.lid.egllt = slid.EGLLT();
213  shwlid.lid.emulll = slid.EMuLLL();
214  shwlid.lid.emullt = slid.EMuLLT();
215  shwlid.lid.epi0lll = slid.EPi0LLL();
216  shwlid.lid.epi0llt = slid.EPi0LLT();
217  shwlid.lid.eplll = slid.EPLLL();
218  shwlid.lid.epllt = slid.EPLLT();
219  shwlid.lid.enlll = slid.ENLLL();
220  shwlid.lid.enllt = slid.ENLLT();
221  shwlid.lid.epilll = slid.EPiLLL();
222  shwlid.lid.epillt = slid.EPiLLT();
223  shwlid.lid.vtxgev = slid.VertexEnergy();
224  shwlid.lid.pi0mass = slid.Pi0mass();
225  shwlid.lid.shwEFrac = slid.ShowerEFrac();
226  shwlid.lid.gap = slid.Gap();
227  shwlid.lid.costheta = slid.CosTheta();
228 
229  std::map<int, float> tPartLongLL = slid.fPartLongLL;
230  std::map<int, float> tPartTransLL = slid.fPartTransLL;
231  std::map<int, float>::iterator itr;
232 
233  itr = tPartLongLL.find(slid::DedxParticleType::kELECTRON);
234  if (itr != tPartLongLL.end()) shwlid.lid.elll = itr->second;
235  else shwlid.lid.elll = -5;
236 
237  itr = tPartTransLL.find(slid::DedxParticleType::kELECTRON);
238  if (itr != tPartTransLL.end()) shwlid.lid.ellt = itr->second;
239  else shwlid.lid.ellt = -5;
240 
241  itr = tPartLongLL.find(slid::DedxParticleType::kMUON);
242  if (itr != tPartLongLL.end()) shwlid.lid.mulll = itr->second;
243  else shwlid.lid.mulll = -5;
244 
245  itr = tPartTransLL.find(slid::DedxParticleType::kMUON);
246  if (itr != tPartTransLL.end()) shwlid.lid.mullt = itr->second;
247  else shwlid.lid.mullt = -5;
248  shwlid.lid.invglll = slid.InverseGLLL();
249 
250  shwlid.lid.dedx0 = slid.Dedx0();
251  shwlid.lid.dedx1 = slid.Dedx1();
252  shwlid.lid.dedx2 = slid.Dedx2();
253  shwlid.lid.dedx3 = slid.Dedx3();
254  shwlid.lid.dedx4 = slid.Dedx4();
255  shwlid.lid.dedx5 = slid.Dedx5();
256  shwlid.lid.dedxp0c0 = slid.Dedxp0c0();
257  shwlid.lid.dedxp0c1 = slid.Dedxp0c1();
258  shwlid.lid.dedxp1c0 = slid.Dedxp1c0();
259  shwlid.lid.dedxp1c1 = slid.Dedxp1c1();
260  shwlid.lid.dedxp2c0 = slid.Dedxp2c0();
261  shwlid.lid.dedxp2c1 = slid.Dedxp2c1();
262  shwlid.lid.dedxp3c0 = slid.Dedxp3c0();
263  shwlid.lid.dedxp3c1 = slid.Dedxp3c1();
264  shwlid.lid.dedxp4c0 = slid.Dedxp4c0();
265  shwlid.lid.dedxp4c1 = slid.Dedxp4c1();
266  shwlid.lid.dedxp5c0 = slid.Dedxp5c0();
267  shwlid.lid.dedxp5c1 = slid.Dedxp5c1();
268 
269  shwlid.lidE.E = slid.NueEnergy();
270  shwlid.lidE.depE = slid.ShowerDepEnergy();
271  // shwlid.lidE = slid.ShowerEnergy();
272  // shwlid.lidE = slid.HadronicEnergy();
273  }
274 
275  //......................................................................
276  void FillSpidVars(const slid::ShowerPID& spid,
277  caf::SRSPid& srspid)
278  {
279  srspid.annepi0 = spid.ValueEPi0();
280  srspid.annepi0el = spid.ValueEPi0EL();
281 
282  srspid.elll = spid.ELLL();
283  srspid.ellt = spid.ELLT();
284  srspid.glll = spid.GLLL();
285  srspid.gllt = spid.GLLT();
286  srspid.mulll = spid.MuLLL();
287  srspid.mullt = spid.MuLLT();
288  srspid.pi0lll = spid.Pi0LLL();
289  srspid.pi0llt = spid.Pi0LLT();
290  srspid.plll = spid.PLLL();
291  srspid.pllt = spid.PLLT();
292  srspid.nlll = spid.NLLL();
293  srspid.nllt = spid.NLLT();
294  srspid.pilll = spid.PiLLL();
295  srspid.pillt = spid.PiLLT();
296  srspid.dedx0 = spid.Dedx0();
297  srspid.dedx1 = spid.Dedx1();
298  srspid.dedx2 = spid.Dedx2();
299  srspid.dedx3 = spid.Dedx3();
300  srspid.dedx4 = spid.Dedx4();
301  srspid.dedx5 = spid.Dedx5();
302 
303  }
304 
305  //......................................................................
307  caf::SRFuzzyKProng& png)
308  {
309  png.shwlid.lidE.E = slid.NueEnergy();
310  png.shwlid.lidE.depE = slid.ShowerDepEnergy();
311  png.shwlid.lidE.shwE = slid.ShowerEnergy();
312  png.shwlid.lidE.hadE = slid.HadronicEnergy();
313  }
314 
315  //......................................................................
317  caf::SRLem& srlem)
318  {
319  srlem.pid = lem.Value();
320  srlem.pidfit = lem.vars.fPIDFit;
321  srlem.pidexp = lem.vars.fPIDExp;
322  srlem.pidexpfit = lem.vars.fPIDExpFit;
323  srlem.meany = lem.vars.fMeanY;
324  srlem.meanyfit = lem.vars.fMeanYFit;
325  srlem.meanyexp = lem.vars.fMeanYExp;
326  srlem.meanysig = lem.vars.fMeanYSig;
327  srlem.meanysigexp = lem.vars.fMeanYSigExp;
328  srlem.meanqfrac = lem.vars.fMeanQFrac;
329  srlem.meanqfracfit = lem.vars.fMeanQFracFit;
330  srlem.meanqfracexp = lem.vars.fMeanQFracExp;
331  srlem.meanqfracsig = lem.vars.fMeanQFracSig;
333  srlem.meanqfracbkg = lem.vars.fMeanQFracBkg;
335  srlem.energydiff = lem.vars.fEnergyDiff;
336  srlem.energydiffexp = lem.vars.fEnergyDiffExp;
337  srlem.enrichfrac = lem.vars.fEnrichFrac;
338  srlem.enrichfracexp = lem.vars.fEnrichFracExp;
339  srlem.fannid = lem.vars.fFannId;
340  srlem.fannidenrich = lem.vars.fFannIdEnrich;
341  srlem.chisig = lem.vars.fChiSig;
342  srlem.chibkg = lem.vars.fChiBkg;
343  srlem.chisigexp = lem.vars.fChiSigExp;
344  srlem.chibkgexp = lem.vars.fChiBkgExp;
345  srlem.avginvE = lem.vars.fAvgInvE;
346  srlem.avgexpE = lem.vars.fAvgExpE;
347  srlem.avgsigE = lem.vars.fAvgSigE;
348  srlem.tmvabdtg = lem.vars.fTMVABDTG;
349  srlem.tmvabdt = lem.vars.fTMVABDT;
350  srlem.tmvabdtd = lem.vars.fTMVABDTD;
351  srlem.tmvabdtmf = lem.vars.fTMVABDTMF;
352  srlem.tmvamlp = lem.vars.fTMVAMLP;
353  srlem.dectree = lem.vars.fDecTreePID;
354  }
355 
356  //......................................................................
357  void FillRVPVars(const rvp::RVP& rvp, caf::SRRvp& srrvp)
358  {
359  srrvp.pid = rvp.Value();
360  srrvp.ncell = rvp.NCells();
361  srrvp.recoE = rvp.RecoE();
362  srrvp.longtr = rvp.LongTr();
363  srrvp.longtrfrac = rvp.LongTrFrac();
364  srrvp.mipfrac = rvp.MipFrac();
365  srrvp.miphits = rvp.MipHits();
366  srrvp.epl20frac = rvp.EPl20Frac();
367  srrvp.efrac2plwin = rvp.EFrac2PlWin();
368  srrvp.efrac6plwin = rvp.EFrac6PlWin();
369  srrvp.efrac2sigrd = rvp.EFrac2SigRd();
370  srrvp.prongs3D = rvp.Prongs3D();
371  srrvp.prongEbal3D = rvp.ProngEBal3D();
372  srrvp.prongs2D = rvp.Prongs2D();
373  srrvp.prongEbal2D = rvp.ProngEBal2D();
374  srrvp.eiso3sig = rvp.EIso3Sig();
375  srrvp.rvp12 = rvp.RVP12();
376  }
377 
378  //......................................................................
379  void FillXnueVars(const xnue::Xnue& xnue, caf::SRXnue& srxnue)
380  {
381  srxnue.bdt = xnue.Value();
382  srxnue.p1energy = xnue.P1Energy();
383  srxnue.p1Fmip = xnue.P1Fmip();
384  srxnue.p2Fmip = xnue.P2Fmip();
385  srxnue.emaxfrac6p = xnue.Emaxfrac6p();
386  srxnue.efrac10p = xnue.Efrac10p();
387  srxnue.efracp2 = xnue.Efracp2();
388  srxnue.efracp3 = xnue.Efracp3();
389  srxnue.efracp4 = xnue.Efracp4();
390  srxnue.efrac2sig = xnue.Efrac2sig();
391  }
392 
393  //......................................................................
395  caf::SRCosRej& srcosrej,
396  caf::SRContain& srcontain)
397  {
398  srcosrej.anglekal = cosrej.AngleKal();
399  srcosrej.anglecos = cosrej.AngleCos();
400 
401  srcosrej.numucontpid2020 = cosrej.ConCosPID();
402  srcosrej.numucontpid2019 = cosrej.OldCosPID();
403  srcosrej.numuunconttunedpid = cosrej.UnconTunedCosPID();
404 
405  srcosrej.nkal3d = cosrej.NTracks3D();
406  srcosrej.mindist = cosrej.MinDist();
407  srcosrej.mincell = cosrej.MinCell();
408  srcosrej.eratio = cosrej.ERatio();
409  srcosrej.scatt = cosrej.Scatt();
410  srcosrej.hadE = cosrej.StartAct();
411  srcosrej.costhetatrue = cosrej.CosThetaTrue();
412  srcosrej.kalthetatrue = cosrej.KalThetaTrue();
413  srcosrej.pdgbest = cosrej.PDGBest();
414  srcosrej.fscattmax = cosrej.FScattMax();
415  srcosrej.fscattsum = cosrej.FScattSum();
416  srcosrej.fscattext = cosrej.FScattExt();
417  srcosrej.fscattsig = cosrej.FScattSig();
418  srcosrej.kfitspeed = cosrej.KFitSpeed();
419  srcosrej.cfitspeed = cosrej.CFitSpeed();
420  srcosrej.kdirscore = cosrej.KDirScore();
421  srcosrej.cdirscore = cosrej.CDirScore();
422  srcosrej.kscorediff = cosrej.KScoreDiff();
423  srcosrej.cscorediff = cosrej.CScoreDiff();
424  srcosrej.kalchisq = cosrej.KalChisq();
425  srcosrej.kalslope = cosrej.KalSlope();
426  srcosrej.kalchidiff = cosrej.KalChiDiff();
427  srcosrej.coschisq = cosrej.CosChisq();
428  srcosrej.cosslope = cosrej.CosSlope();
429  srcosrej.coschidiff = cosrej.CosChiDiff();
430 
431  srcontain.vtxdist = cosrej.VtxDist();
432  srcontain.enddist = cosrej.EndDist();
433  srcontain.cosfwddist = cosrej.CosFwdDist();
434  srcontain.cosbakdist = cosrej.CosBakDist();
435  srcontain.cosfwdcell = cosrej.CosFwdCell();
436  srcontain.cosfwdcellnd = cosrej.CosFwdCellND();
437  srcontain.cosbakcell = cosrej.CosBakCell();
438  srcontain.cosbakcellnd = cosrej.CosBakCellND();
439  srcontain.cosfwdair = cosrej.CosFwdAir();
440  srcontain.cosfwdsteel = cosrej.CosFwdSteel();
441  srcontain.cosbakair = cosrej.CosBakAir();
442  srcontain.cosbaksteel = cosrej.CosBakSteel();
443  srcontain.kalfwdair = cosrej.KalFwdAir();
444  srcontain.kalfwdsteel = cosrej.KalFwdSteel();
445  srcontain.kalbakair = cosrej.KalBakAir();
446  srcontain.kalbaksteel = cosrej.KalBakSteel();
447  srcontain.kalfwddist = cosrej.KalFwdDist();
448  srcontain.kalbakdist = cosrej.KalBakDist();
449  srcontain.kalfwdcell = cosrej.KalFwdCell();
450  srcontain.kalfwdcellnd = cosrej.KalFwdCellND();
451  srcontain.kalbakcell = cosrej.KalBakCell();
452  srcontain.kalbakcellnd = cosrej.KalBakCellND();
453  srcontain.cosyposattrans = cosrej.CosYPosAtTrans();
454  srcontain.kalyposattrans = cosrej.KalYPosAtTrans();
455  }
456 
457  //......................................................................
458  void FillNueCosRejVars(const cosrej::NueCosRej& nuecosrej,
459  caf::SRNueCosRej& srnuecosrej)
460  {
461  srnuecosrej.hitsperplane = nuecosrej.HitsPerPlane();
462  srnuecosrej.pngptp = nuecosrej.ProngTransMom();
463  srnuecosrej.partptp = nuecosrej.ParticleShowerTransMom();
464  srnuecosrej.photptp = nuecosrej.PhotonShowerTransMom();
465  srnuecosrej.photpxp = nuecosrej.PhotonShowerMomX();
466  srnuecosrej.photpyp = nuecosrej.PhotonShowerMomY();
467  srnuecosrej.starttop = nuecosrej.StartDistToTop();
468  srnuecosrej.startbottom = nuecosrej.StartDistToBottom();
469  srnuecosrej.startfront = nuecosrej.StartDistToFront();
470  srnuecosrej.startback = nuecosrej.StartDistToBack();
471  srnuecosrej.startwest = nuecosrej.StartDistToWest();
472  srnuecosrej.starteast = nuecosrej.StartDistToEast();
473  srnuecosrej.stoptop = nuecosrej.StopDistToTop();
474  srnuecosrej.stopbottom = nuecosrej.StopDistToBottom();
475  srnuecosrej.stopfront = nuecosrej.StopDistToFront();
476  srnuecosrej.stopback = nuecosrej.StopDistToBack();
477  srnuecosrej.stopwest = nuecosrej.StopDistToWest();
478  srnuecosrej.stopeast = nuecosrej.StopDistToEast();
479  srnuecosrej.cospidcontain= nuecosrej.CosPIDContain();
480  srnuecosrej.cospidcontainxy= nuecosrej.CosPIDContainXY();
481  srnuecosrej.cospidlight = nuecosrej.CosPIDLight();
482  srnuecosrej.cospidperibdt = nuecosrej.PeriBDT();
483  srnuecosrej.cospidcorebdt = nuecosrej.CoreBDT();
484  srnuecosrej.cosdang = nuecosrej.CosAngleToNextProng();
485  srnuecosrej.vtxdoca = nuecosrej.ProngDistToVtx();
486  srnuecosrej.prongmaxx = nuecosrej.ProngMaxX();
487  srnuecosrej.prongmaxy = nuecosrej.ProngMaxY();
488  srnuecosrej.prongmaxz = nuecosrej.ProngMaxZ();
489  srnuecosrej.prongminx = nuecosrej.ProngMinX();
490  srnuecosrej.prongminy = nuecosrej.ProngMinY();
491  srnuecosrej.prongminz = nuecosrej.ProngMinZ();
492  srnuecosrej.distallpngtop = std::min(nuecosrej.StartMinDistToTop(),
493  nuecosrej.StopMinDistToTop());
494  srnuecosrej.distallpngbottom = std::min(nuecosrej.StartMinDistToBottom(),
495  nuecosrej.StopMinDistToBottom());
496  srnuecosrej.distallpngfront = std::min(nuecosrej.StartMinDistToFront(),
497  nuecosrej.StopMinDistToFront());
498  srnuecosrej.distallpngback = std::min(nuecosrej.StartMinDistToBack(),
499  nuecosrej.StopMinDistToBack());
500  srnuecosrej.distallpngwest = std::min(nuecosrej.StartMinDistToWest(),
501  nuecosrej.StopMinDistToWest());
502  srnuecosrej.distallpngeast = std::min(nuecosrej.StartMinDistToEast(),
503  nuecosrej.StopMinDistToEast());
504 
505  srnuecosrej.sparsenessasymm = nuecosrej.SparsenessAsymm();
506  srnuecosrej.hitsperplaneasymm = nuecosrej.HitsPerPlaneAsymm();
507  srnuecosrej.sparsenessasymmslice = nuecosrej.SparsenessAsymmSlice();
508  srnuecosrej.hitsperplaneasymmslice = nuecosrej.HitsPerPlaneAsymmSlice();
509  srnuecosrej.musliceidxbydist = nuecosrej.MuSliceIdxByDist();
510  srnuecosrej.muanglediffbydist = nuecosrej.MuAngleDiffByDist();
511  srnuecosrej.mutimediffbydist = nuecosrej.MuTimeDiffByDist();
512  srnuecosrej.muclosestapproachbydist = nuecosrej.MuClosestApproachByDist();
513  srnuecosrej.musliceidxbytime = nuecosrej.MuSliceIdxByTime();
514  srnuecosrej.muanglediffbytime = nuecosrej.MuAngleDiffByTime();
515  srnuecosrej.mutimediffbytime = nuecosrej.MuTimeDiffByTime();
516  srnuecosrej.muclosestapproachbytime = nuecosrej.MuClosestApproachByTime();
517 
518  }
519 
520  //......................................................................
521  void FillNCCosRejVars(const ncid::NCCosRej& nccosrej,
522  caf::SRNCCosRej& srnccosrej)
523  {
524  srnccosrej.pngptp = nccosrej.ProngTransMom();
525  srnccosrej.cospidbdt = nccosrej.CosPIDDT();
526  srnccosrej.cospidbdtg = nccosrej.CosPIDDTG();
527  }
528 
529  //......................................................................
530  void FillNCPi0BkgRejVars(const ncpi0::NCPi0BkgRej& ncpi0bkgrej,
531  caf::SRNCPi0BkgRej& srncpi0bkgrej)
532  {
533  srncpi0bkgrej.ncpidbdtg = ncpi0bkgrej.NCPIDBDTG();
534  srncpi0bkgrej.ncpidbdtglt = ncpi0bkgrej.NCPIDBDTGLT();
535  }
536 
537 
538  //......................................................................
540  caf::SRPresel& srnuepre)
541  {
542  srnuepre.passpre = nuepre.PassPresel();
543  }
544 
545 
546  //......................................................................
548  caf::SRPresel& srrockpre)
549  {
550  srrockpre.passpre = rockpre.PassPresel();
551  }
552 
553  //......................................................................
554  void FillVetoVars(const presel::Veto& veto,
555  const presel::Veto& nueveto,
556  caf::SRVeto& srveto)
557  {
558  srveto.ncell = veto.NCell();
559  srveto.fwddist = veto.FwdDist();
560  srveto.bakdist = veto.BakDist();
561  srveto.ratio = veto.Ratio();
562  srveto.diry = veto.DirY();
563  srveto.angle = veto.Angle();
564  srveto.anglevar = veto.AngleVar();
565  srveto.anglevarold = veto.AngleVarOld();
566  srveto.passthru = veto.PassThru();
567  srveto.passthruold = veto.PassThruOld();
568  srveto.passangle = veto.PassAngle();
569  srveto.passanglefirstana = veto.PassAngleFirstAna();
570  srveto.passnumicut = veto.PassNuMICut();
571  srveto.keep = veto.Keep();
572  srveto.keepnue = nueveto.Keep();
573  }
574 
575  //......................................................................
577  caf::SRNuonEResult& srnuone)
578  {
579  srnuone.setDefault();
580 
581  srnuone.output = result.fOutput;
582 
583  srnuone.noutput = result.fOutput.size();
584 
585  if(srnuone.noutput > 0) srnuone.rawgevslice = result.fOutput[0];
586  if(srnuone.noutput > 1) srnuone.nuoneid = result.fOutput[1];
587  if(srnuone.noutput > 2) srnuone.pi0id = result.fOutput[2];
588  if(srnuone.noutput > 3) srnuone.nueccid = result.fOutput[3];
589  if(srnuone.noutput > 4) srnuone.otherid = result.fOutput[4];
590  if(srnuone.noutput > 5) srnuone.epi0nuoneid = result.fOutput[5];
591  if(srnuone.noutput > 6) srnuone.epi0pi0id = result.fOutput[6];
592  }
593 
594  //......................................................................
596  caf::SRCVNResult& srcvn,
597  bool NoCosmics)
598  {
599  srcvn.output = result.fOutput;
600 
601  std::vector<float>::const_iterator maxIt =
602  std::max_element(result.fOutput.cbegin(), result.fOutput.cend());
603 
604  srcvn.maxval = *maxIt;
605  srcvn.argmax = std::distance(result.fOutput.cbegin(), maxIt);
606  srcvn.noutput = result.fOutput.size();
607 
608  srcvn.numuid = result.fOutput[0];
609  srcvn.nueid = result.fOutput[1];
610  srcvn.nutauid = -5;
611  srcvn.ncid = result.fOutput[2];
612  srcvn.cosmicid = result.fOutput[3];
613  if (NoCosmics) srcvn.cosmicid = -5;
614  }
615 
616  //......................................................................
618  caf::SRPixelObjMap& srcvnmap, bool useGeV)
619  {
620  srcvnmap.nchan = 2; //two views to pixel map
621  srcvnmap.nplanes = map.fNPlane/2;
622  srcvnmap.ncells = map.fNCell;
623  srcvnmap.hitfracx = map.fFracX;
624  srcvnmap.hitfracy = map.fFracY;
625  srcvnmap.firstplane = map.fBound.FirstPlane();
626  srcvnmap.lastplane = map.fBound.LastPlane();
627  srcvnmap.firstcellx = map.fBound.FirstCell(0);
628  srcvnmap.lastcellx = map.fBound.LastCell(0);
629  srcvnmap.firstcelly = map.fBound.FirstCell(1);
630  srcvnmap.lastcelly = map.fBound.LastCell(1);
631 
632  std::vector<unsigned char> tval = map.PixelMapToVector(useGeV);
633  std::copy(tval.begin(), tval.end(), srcvnmap.cvnmap);
634  std::vector<unsigned char> tlab = map.PixelMapToObjVector(true);
635  std::copy(tlab.begin(), tlab.end(), srcvnmap.cvnlabmap);
636  std::vector<unsigned char> tobj = map.PixelMapToObjVector(false);
637  std::copy(tobj.begin(), tobj.end(), srcvnmap.cvnobjmap);
638  }
639 
640  //......................................................................
642  caf::SRPixelMap& srcvnmap, bool useGeV)
643  {
644  srcvnmap.nchan = 2; //two views to pixel map
645  srcvnmap.nplanes = map.fNPlane/2;
646  srcvnmap.ncells = map.fNCell;
647  srcvnmap.hitfracx = map.fFracX;
648  srcvnmap.hitfracy = map.fFracY;
649  srcvnmap.firstplane = map.fBound.FirstPlane();
650  srcvnmap.lastplane = map.fBound.LastPlane();
651  srcvnmap.firstcellx = map.fBound.FirstCell(0);
652  srcvnmap.lastcellx = map.fBound.LastCell(0);
653  srcvnmap.firstcelly = map.fBound.FirstCell(1);
654  srcvnmap.lastcelly = map.fBound.LastCell(1);
655 
656  std::vector<unsigned char> tval = map.PixelMapToVector(useGeV);
657  std::copy(tval.begin(), tval.end(), srcvnmap.cvnmap);
658  }
659 
660  //......................................................................
662  caf::SRTrainingData& srtdata)
663  {
664  srtdata.parent = cvntdata.fPPt;
665  srtdata.interaction = cvntdata.fInt;
666  srtdata.finalstate = cvntdata.fFSt;
667  srtdata.finalstateprong = cvntdata.fFSp;
668  srtdata.particles = cvntdata.fFSlong;
669  srtdata.nuenergy = cvntdata.fNuEnergy;
670  srtdata.lepenergy = cvntdata.fLepEnergy;
671  srtdata.vtxx = cvntdata.fVtxX;
672  srtdata.vtxy = cvntdata.fVtxY;
673  srtdata.vtxz = cvntdata.fVtxZ;
674  }
675 
676  //......................................................................
678  caf::SRProngTrainingData& srpdata)
679  {
680  srpdata.label3d = cvnpdata.fLabel3D;
681  srpdata.labelx = cvnpdata.fLabelX;
682  srpdata.labely = cvnpdata.fLabelY;
683  srpdata.isprimary = cvnpdata.fIsPrimary;
684  srpdata.purity3d = cvnpdata.fPurity3D;
685  srpdata.purityx = cvnpdata.fPurityX;
686  srpdata.purityy = cvnpdata.fPurityY;
687  srpdata.rece = cvnpdata.fRecE;
688  srpdata.ncellx = cvnpdata.fNCellX;
689  srpdata.ncelly = cvnpdata.fNCellY;
690  srpdata.vertx = cvnpdata.fVertX;
691  srpdata.verty = cvnpdata.fVertY;
692  srpdata.vertz = cvnpdata.fVertZ;
693  srpdata.plength = cvnpdata.fPLength;
694  srpdata.pangle = cvnpdata.fPAngle;
695  srpdata.pgap = cvnpdata.fPGap;
696  }
697 
698  //......................................................................
699  void FillCVNFeaturesVars(const cvn::Features& features,
700  caf::SRCVNFeatures& srcvnfeatures,
701  int maxcomponents)
702  {
703  srcvnfeatures.output = features.fOutput;
704  srcvnfeatures.noutput = features.fOutput.size();
705  srcvnfeatures.components = features.fPrincipalComponents;
706  if(int(srcvnfeatures.components.size()) > maxcomponents)
707  srcvnfeatures.components.resize(maxcomponents);
708  srcvnfeatures.ncomponents = srcvnfeatures.components.size();
709  }
710 
711  //......................................................................
713  {
714  remid.pid = srTrk.rempid;
715  remid.scatllh = srTrk.scatllh;
716  remid.dedxllh = srTrk.dedxllh;
717  remid.len = srTrk.len;
718  remid.measfrac = srTrk.measfrac;
719  remid.cont = srTrk.remcont;
720  }
721 
722  //......................................................................
723  void CopyMuonIDVars(const SRKalmanTrack& srTrk, caf::SRMuonID& muid)
724  {
725  muid.pid = srTrk.muonid;
726  }
727 
729  const SliceLID::Prediction& artSliceLID, caf::SRSliceLID& cafSliceLID
730  )
731  {
732  cafSliceLID.ncid = artSliceLID.nc;
733  cafSliceLID.nueid = artSliceLID.nue;
734  cafSliceLID.numuid = artSliceLID.numu;
735  cafSliceLID.nutauid = artSliceLID.nutau;
736  cafSliceLID.cosmicid = artSliceLID.cosmic;
737  }
738 
739 
740 } // end namespace caf
741 ////////////////////////////////////////////////////////////////////////
int cosfwdcell
cosmictrack projected # cells from end point forwards to det edge
Definition: SRContain.h:30
NuonE output.
Definition: SRNuonEResult.h:13
float ELLT() const
Definition: ShowerPID.h:68
float ncid
Likelihood Neutral Current.
Definition: SRCVNResult.h:23
::xsd::cxx::tree::id< char, ncname > id
Definition: Database.h:165
double ANNEL() const
Definition: EID.cxx:34
float efracp4
Definition: SRXnue.h:24
unsigned char cvnlabmap[16000]
store 1D array with hit labels
Definition: SRPixelObjMap.h:33
int ncell
number of cells in slice
Definition: SRRvp.h:20
int ismuon
Boolean to distinguish muons from electrons.
Definition: SRSLid.h:26
int PDGBest() const
Definition: CosRejObj.cxx:736
SRShowerLID shwlid
Shower information.
Definition: SRFuzzyKProng.h:19
double CosAngleToNextProng() const
Cosine of Angle between the leading and the next leading prong.
Definition: NueCosRej.h:164
float KalThetaTrue() const
Definition: CosRejObj.cxx:730
float dedxp1c1
Definition: SRSLid.h:79
float pllt
Proton ll for ltransverse shower.
Definition: SRSPid.h:33
int FirstCell(const unsigned int &view) const
Definition: Boundary.h:44
double ANNNoCos() const
Definition: EID.cxx:61
float annENoCos
ann output, with E without cosTheta
Definition: SRJMEid.h:27
float EIso3Sig() const
Definition: RVP.h:74
double ELLL() const
Definition: EID.cxx:798
double MuClosestApproachByDist() const
Distance of closest approach between this slice and the possible muon parent (selected by closest app...
Definition: NueCosRej.h:224
float gap
gap from vertex to start of shower
Definition: SRJMEid.h:57
float KalChisq() const
Definition: CosRejObj.cxx:538
double EPiLLT() const
Definition: EID.cxx:748
SRParticleTruth truthXView
Truth information for the prong.
Definition: SRProng.h:47
float pilll
Pion ll for longitudinal shower.
Definition: SRSPid.h:36
float dedx4
Plane Dedx.
Definition: SRJMEid.h:40
float dedxp4c1
Definition: SRSLid.h:85
Contains the SPID pid information for a shower (slid::ShowerLID objects).
Definition: SRSPid.h:15
float MuLLL() const
Definition: ShowerPID.h:71
double bdt
Definition: SRXnue.h:26
unsigned int noutput
Number of entries in output vector.
Definition: SRNuonEResult.h:27
double NCPIDBDTGLT() const
Definition: NCPi0BkgRej.h:22
double ProngMaxY() const
The maximum Y between the start and stop points of the leading prong in the slice.
Definition: NueCosRej.h:176
float ValueEPi0() const
Definition: ShowerPID.h:58
float fscattsum
sum of scattering angles (Fernanda)
Definition: SRCosRej.h:43
float scatllh
Log-likelihood value from scattering angle.
Definition: SRRemid.h:26
unsigned int interaction
Class of the event.
unsigned int nplanes
number of planes in pixel map
Definition: SRPixelObjMap.h:17
float ann
ann output
Definition: SRJMEid.h:21
float EFrac2PlWin() const
Definition: RVP.h:67
unsigned short maxplanecont
maximum number of contiguous planes in prong
Definition: SRShower.h:29
float enllt
Electron - neutron LL for transverse shower.
Definition: SRSLid.h:45
void FillCVNNeutronDaughterResultVars(const std::vector< art::Ptr< rb::PID > > &cvnneutrons, caf::SRCVNNeutronDaughterResult &cvnneutron)
Definition: FillPIDs.cxx:95
float neutronid
Likelihood neutron.
int labelx
Class of the X view prong.
float ValueWithE() const
Definition: ShowerLID.h:103
float Efracp2() const
Definition: Xnue.h:52
float shwE
Energy of shower [GeV].
Definition: SRSLidEnergy.h:21
unsigned int NCells() const
Definition: RVP.h:60
Variables describing Michel E&#39;s found around the end of a track.
Definition: SRPixelMap.h:13
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
Definition: Cluster.cxx:134
float CosSlope() const
Definition: CosRejObj.cxx:544
unsigned int firstcellx
global position of first x cell in map
Definition: SRPixelMap.h:26
double MuLLT() const
Definition: EID.cxx:850
float ellt
Electron ll for ltransverse shower.
Definition: SRSPid.h:25
float anne
ann with energy output
Definition: SRSLid.h:23
InteractionType fInt
Class of the event.
Definition: TrainingData.h:60
float KalSlope() const
Definition: CosRejObj.cxx:532
float kalchisq
chisq value of slope fit (kalmantrack)
Definition: SRCosRej.h:62
float pionid
Likelihood pion.
double Dedx4() const
Definition: EID.cxx:622
float epi0llt
Electron - Pi0 ll for ltransverse shower.
Definition: SRJMEid.h:47
float glll
Gamma ll for longitudinal shower.
Definition: SRSPid.h:26
float Dedx0() const
Definition: ShowerLID.h:126
float PLLL() const
Definition: ShowerPID.h:75
double MuClosestApproachByTime() const
Distance of closest approach between this slice and the possible muon parent (selected by difference ...
Definition: NueCosRej.h:241
float P1Fmip() const
Definition: Xnue.h:48
float efrac6plwin
largest fraction of energy in a 6 plane window out of total slice
Definition: SRRvp.h:28
const art::PtrVector< rb::CellHit > & XCells() const
Get all cells from the x-view.
Definition: Cluster.h:124
float muonid
Likelihood muon.
double StopMinDistToBottom() const
Minimum perpendicular distance of all prongs stop points to the bottom edge of the detector...
Definition: NueCosRej.h:124
bool PassPresel() const
Definition: PreselObj.cxx:30
SRVector3D dir
Shower direction at start point [unit vector recommended].
Definition: SRShower.h:33
float dedx0
Plane dedx for e/photon/pi0 identification.
Definition: SRSPid.h:38
float elll
electron ll for longitudinal shower
Definition: SRJMEid.h:61
float nueid
Likelihood of Charge Current NuE.
Definition: SRSliceLID.h:14
float Pi0LLL() const
Definition: ShowerPID.h:73
std::vector< unsigned char > PixelMapToObjVector(bool useLabels) const
Definition: PixelMap.cxx:417
double ANNEPi0EL() const
Definition: EID.cxx:52
unsigned int ncells
number of cells in pixel map
Definition: SRPixelObjMap.h:21
float shwEFrac
fraction of energy of leading shower out of total energy of slice
Definition: SRSLid.h:50
float ENLLT() const
Definition: ShowerLID.h:114
double EGLLL() const
Definition: EID.cxx:638
float distallpngbottom
Definition: SRNueCosRej.h:95
double EPiLLL() const
Definition: EID.cxx:738
float Dedxp2c1() const
Definition: ShowerLID.h:138
float emullt
Electron - muon ll for ltransverse shower.
Definition: SRJMEid.h:45
std::vector< float > components
Principal components of output.
Definition: SRCVNFeatures.h:18
float prongEbal3D
the energy balance between the two highest energy 3D prongs in slice
Definition: SRRvp.h:31
This is a helper class for ParticleIDAlg that provides a tidy structure in which to hold the dE/dx hi...
float CosChisq() const
Definition: CosRejObj.cxx:550
float CosFwdDist() const
Definition: CosRejObj.cxx:556
float KalYPosAtTrans() const
Definition: CosRejObj.cxx:682
float KalBakDist() const
Definition: CosRejObj.cxx:610
float Dedx3() const
Definition: ShowerLID.h:129
float meanyexp
Hadronic y of matches, "exp".
Definition: SRLem.h:32
float costhetatrue
cosine of angle between cosmic trk dir and true dir of most contributing particle ...
Definition: SRCosRej.h:26
float ann
ann output
Definition: SRSLid.h:22
void FillVetoVars(const presel::Veto &veto, const presel::Veto &nueveto, caf::SRVeto &srveto)
Definition: FillPIDs.cxx:554
float EGLLL() const
Definition: ShowerLID.h:105
float cfitspeed
fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for cosmic track ...
Definition: SRCosRej.h:56
float mutimediffbydist
Definition: SRNueCosRej.h:193
float Dedx2() const
Definition: ShowerPID.h:62
float eglll
Plane Dedx.
Definition: SRJMEid.h:42
float Dedx0() const
Definition: ShowerPID.h:60
float muanglediffbytime
Definition: SRNueCosRej.h:206
unsigned int lastcelly
global position of last y cell in map
Definition: SRPixelMap.h:29
float vtxgev
Energy of slice in vertex region.
Definition: SRJMEid.h:54
float Dedx2() const
Definition: ShowerLID.h:128
bool IsMuon() const
Definition: EID.cxx:547
double MuAngleDiffByDist() const
Minimum angle difference between this slice and the possible muon parent (selected by closest approac...
Definition: NueCosRej.h:220
float dedx1
Definition: SRSLid.h:67
float EPiLLT() const
Definition: ShowerLID.h:116
X or Y views.
Definition: PlaneGeo.h:30
bool PassThru() const
Definition: Veto.cxx:100
float KalBakSteel() const
Definition: CosRejObj.cxx:622
float cospidbdtg
NC cosrej PID. A Keras using 15 variables.
Definition: SRNCCosRej.h:23
float E
Energy [GeV].
Definition: SRSLidEnergy.h:18
float mullt
Muon ll for ltransverse shower.
Definition: SRSPid.h:29
float len
Track length of identified track.
Definition: SRRemid.h:28
float p1Fmip
Definition: SRXnue.h:18
float UnconTunedCosPID() const
Definition: CosRejObj.cxx:478
float epillt
Electron - pion ll for ltransverse shower.
Definition: SRJMEid.h:53
double MuAngleDiffByTime() const
Minimum angle difference between this slice and the possible muon parent (selected by difference in t...
Definition: NueCosRej.h:237
float KalFwdSteel() const
Definition: CosRejObj.cxx:604
float Dedxp4c0() const
Definition: ShowerLID.h:141
float cosyposattrans
Y position of cosmic track at transition (ND only, use to check if went through air gap) ...
Definition: SRContain.h:49
double StopDistToTop() const
Perpendicular distance of stop point of the leading prong to the top edge of the detector.
Definition: NueCosRej.h:72
void FillMuonIDVars(const muonid::MuonID &muid, SRKalmanTrack &srTrk)
Definition: FillPIDs.cxx:48
unsigned short nplaney
number of planes spanned in y view
Definition: SRShower.h:39
double CoreBDT() const
Definition: NueCosRej.h:158
float nllt
Neutron ll for ltransverse shower.
Definition: SRSPid.h:35
float eplll
Electron - proton ll for longitudinal shower.
Definition: SRJMEid.h:48
double EelLLT() const
Definition: EID.cxx:828
float hitfracy
fraction of y hits in slice in y view pixel map
Definition: SRPixelMap.h:31
double ScatSeparation() const
Return the scattering separation variable used as an input to the kNN that determines the pid value...
Definition: ReMId.cxx:206
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...
Definition: NueCosRej.h:136
float maxval
Maximum value among net outputs.
Definition: SRCVNResult.h:25
float StartAct() const
Definition: CosRejObj.cxx:712
void CopyMuonIDVars(const SRKalmanTrack &srTrk, caf::SRMuonID &muid)
Definition: FillPIDs.cxx:723
double StartMinDistToTop() const
Minimum perpendicular distance of all prongs start points.
Definition: NueCosRej.h:96
void FillCVNFeaturesVars(const cvn::Features &features, caf::SRCVNFeatures &srcvnfeatures, int maxcomponents)
Definition: FillPIDs.cxx:699
float gap
gap from shower start to ElasticArms vertex [cm]
Definition: SRShower.h:40
float pid
The default PID value - normally use this one.
Definition: SRLem.h:24
float p1energy
Definition: SRXnue.h:17
unsigned short nhit
number of hits
Definition: SRShower.h:25
bool passanglefirstana
Would this slice pass the old angle cut used for first analysis.
Definition: SRVeto.h:29
std::map< int, float > fPartLongLL
Map of the longitudinal ll by paricle type.
Definition: ShowerLID.h:238
Vertical planes which measure X.
Definition: PlaneGeo.h:28
double ENLLL() const
Definition: EID.cxx:718
float tmvabdtmf
TMVA boosted decision tree, Fisher discriminant.
Definition: SRLem.h:65
float cosslope
slope of track timing fit (cosmictrack)
Definition: SRCosRej.h:64
float muclosestapproachbydist
Definition: SRNueCosRej.h:197
float kalyposattrans
Y position of Kalman track and transition (ND only, use to check if went through air gap) ...
Definition: SRContain.h:50
float Dedxp1c0() const
Definition: ShowerLID.h:135
unsigned int fNPlane
Number of planes, length of pixel map.
Definition: PixelMap.h:117
float tmvabdt
TMVA boosted decision tree, adaptive boost.
Definition: SRLem.h:63
float epllt
Electron - proton ll for ltransverse shower.
Definition: SRJMEid.h:49
int CosBakCellND() const
Definition: CosRejObj.cxx:646
float dedxp0c1
Definition: SRSLid.h:77
unsigned int ncells
number of cells in pixel map
Definition: SRPixelMap.h:22
float remavededx
average dE/dx value as computed by remid
Definition: SRKalmanTrack.h:15
float Efracp4() const
Definition: Xnue.h:54
int pdgbest
pdg code of most contributing true particle to best kalman trk
Definition: SRCosRej.h:28
float LongTrFrac() const
Definition: RVP.h:63
ParticlesType fFSlong
Class of the event by primary n prongs.
Definition: TrainingData.h:63
void FillJMEIDVars(const jmshower::EID &eid, caf::SRJMEid &sreid)
Definition: FillPIDs.cxx:54
double Dedx3() const
Definition: EID.cxx:615
float epilll
Electron - pion ll for longitudinal shower.
Definition: SRJMEid.h:52
double CosPIDDTG() const
Gradient BDT trained, based on simple containment variables.
Definition: NCCosRej.h:26
float eplll
Electron - proton LL for longitudinal shower.
Definition: SRSLid.h:42
float dedxp3c0
Definition: SRSLid.h:82
float Scatt() const
Definition: CosRejObj.cxx:706
float FScattMax() const
Definition: CosRejObj.cxx:767
float fscattsig
sigma of scattering variable distribution (Fernanda)
Definition: SRCosRej.h:45
Details of processing cuts made by the veto modules.
Definition: SRVeto.h:12
float dedx2
Plane dedx for e/photon/pi0 identification.
Definition: SRSPid.h:40
unsigned int parent
Class of the event by parent particle.
float energydiffexp
Pot. diff between sig and bkg matches, "exp".
Definition: SRLem.h:45
unsigned short nhitx
number of hits in x-view
Definition: SRShower.h:26
float dedx5
Plane dedx for e/photon/pi0 identification.
Definition: SRSPid.h:43
float Emaxfrac6p() const
Definition: Xnue.h:50
float EPl20Frac() const
Definition: RVP.h:66
double PeriBDT() const
Definition: NueCosRej.h:157
ProngType fLabelX
Class of the X view prong.
float ellt
Electron LL for transverse shower.
Definition: SRSLid.h:59
float meanyfit
Hadronic y of matches, "fit".
Definition: SRLem.h:31
float mipfrac
fraction of mip hits out of total number of hits
Definition: SRRvp.h:24
int KalBakCellND() const
Definition: CosRejObj.cxx:670
float cospidbdt
NC Cosrej PID. A real adaptive BDT using 13 variables.
Definition: SRNCCosRej.h:20
Output from Cosmic Rejection (CosRej) module.
Definition: SRCosRej.h:8
float hadE
Hadronic energy [GeV].
Definition: SRSLidEnergy.h:22
float meanqfracbkg
Mean fraction matched in background matches.
Definition: SRLem.h:41
float KalFwdDist() const
Definition: CosRejObj.cxx:592
double fMeanQFracBkgExp
Definition: PIDExtraVars.h:47
std::vector< float > fOutput
Vector of outputs from neural net.
Definition: CVNFeatures.h:23
float enrichfrac
Fraction of matches that are "enriched".
Definition: SRLem.h:47
preselection information
Definition: SRPresel.h:12
double StopDistToWest() const
Perpendicular distance of stop point of the leading prong to the west edge of the detector (west is p...
Definition: NueCosRej.h:88
double NCPIDBDTG() const
Gradient BDT trained, based on prong/shower variables.
Definition: NCPi0BkgRej.h:19
float annepi0
ann for e/pi0 separation
Definition: SRSPid.h:21
float dedxp5c1
Definition: SRSLid.h:87
float calE
energy based on summed calibrated deposited charge [GeV]
Definition: SRShower.h:31
float vtxdist
shortest distance to wall from primary Kalman track start position (proxy for vertex) ...
Definition: SRContain.h:26
unsigned int particles
Class of the event by primary n prongs.
float EPLLL() const
Definition: ShowerLID.h:111
float AngleCos() const
Definition: CosRejObj.cxx:459
double StopDistToFront() const
Perpendicular distance of stop point of the leading prong to the front edge of the detector...
Definition: NueCosRej.h:80
float dedxp0c0
Definition: SRSLid.h:76
bool Keep() const
Definition: Veto.cxx:125
double HitsPerPlaneAsymmSlice() const
Asymmetry in the hitsperplane at the start and end of the slice.
Definition: NueCosRej.h:207
float annepi0
e/pi0 PID for oscillation analysis
Definition: SRSLid.h:24
int cosbakcell
cosmictrack projected # cells from start point backwards to det edge
Definition: SRContain.h:35
float NLLL() const
Definition: ShowerPID.h:77
float Dedxp5c0() const
Definition: ShowerLID.h:143
double PhotonShowerTransMom() const
Transverse component of the event momentum. Particle masses in the momentum calculation are assumed t...
Definition: NueCosRej.h:36
float muanglediffbydist
Definition: SRNueCosRej.h:189
int LastPlane() const
Definition: Boundary.h:43
bool passpre
Bool saying if event passed Preselection.
Definition: SRPresel.h:18
virtual TVector3 Start() const
Definition: Prong.h:73
float kalfwdair
for Kalmantrack projected distance forwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:42
Features, basic output of CVN neural net.
Definition: CVNFeatures.h:15
PixelMap for CVN.
bool PassAngleFirstAna() const
Definition: Veto.cxx:115
float CosBakSteel() const
Definition: CosRejObj.cxx:586
unsigned distance(const T &t1, const T &t2)
float costheta
cosine of track with respect to beam direction
Definition: SRSLid.h:52
float scatt
sum of all all angular deviation (looping over kaltrk traj pts) / trk len
Definition: SRCosRej.h:41
float chibkgexp
chisq of a fit assuming best match is bkg, "exp"
Definition: SRLem.h:56
float CosTheta() const
Definition: ShowerLID.h:121
float cscorediff
chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit modu...
Definition: SRCosRej.h:60
int KalFwdCell() const
Definition: CosRejObj.cxx:652
float Dedxp3c1() const
Definition: ShowerLID.h:140
double CosPIDLight() const
A BDT trained based on simple containment variables and the nue output of CVN. Includes pX/p and pY/p...
Definition: NueCosRej.h:155
int kalbakcellnd
Kalmantrack projected # cells from start point backwards to det edge, including muon catcher...
Definition: SRContain.h:46
bool passthru
Would this slice pass the (fixed) through-going check?
Definition: SRVeto.h:26
double EMuLLL() const
Definition: EID.cxx:658
float nutauid
Likelihood Charge Current NuTau.
Definition: SRCVNResult.h:22
std::vector< float > fOutput
Vector of outputs from neural net.
Definition: Result.h:30
float BakDist() const
Definition: Veto.cxx:70
float KalBakAir() const
Definition: CosRejObj.cxx:616
double AvgDedx() const
Definition: ReMId.cxx:68
bool passangle
Would this slice pass the most recent angle cut?
Definition: SRVeto.h:28
bool passthruold
Would this slice pass the (not fixed, no ratio check) through-going check?
Definition: SRVeto.h:27
double Value() const
Definition: PID.h:22
Horizontal planes which measure Y.
Definition: PlaneGeo.h:29
void FillSlidVars(const slid::ShowerLID &slid, SRShowerLID &shwlid)
Definition: FillPIDs.cxx:202
double Dedx1() const
Definition: EID.cxx:601
float emulll
Electron - muon LL for longitudinal shower.
Definition: SRSLid.h:38
float Efrac10p() const
Definition: Xnue.h:51
float dedx3
Plane Dedx.
Definition: SRJMEid.h:39
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
Definition: Cluster.cxx:180
double Gap() const
Definition: EID.cxx:788
float nueid
Likelihood Charge Current NuE.
Definition: SRCVNResult.h:21
float tmvabdtd
TMVA boosted decision tree, decorrelation.
Definition: SRLem.h:64
float pngptp
Event transverse momentum fraction, based on prongs. See ncid::NCCosRej::ProngTransMom.
Definition: SRNCCosRej.h:17
void FillSliceLID(const SliceLID::Prediction &artSliceLID, caf::SRSliceLID &cafSliceLID)
Definition: FillPIDs.cxx:728
float invglll
photon LL for longitudinal shower, assuming shower started at opposite end
Definition: SRSLid.h:57
float kalslope
slope of track timing fit (kalmantrack)
Definition: SRCosRej.h:61
float LongTr() const
Definition: RVP.h:62
float annNoCos
ann output, without cosTheta
Definition: SRJMEid.h:26
float EPi0LLT() const
Definition: ShowerLID.h:110
float emullt
Electron - muon LL for transverse shower.
Definition: SRSLid.h:39
double Dedx2() const
Definition: EID.cxx:608
float kalfwdsteel
for Kalmantrack projected distance forwards, how much is through steel (ND only, currently is just di...
Definition: SRContain.h:43
float dedxp5c0
Definition: SRSLid.h:86
float epilll
Electron - pion LL for longitudinal shower.
Definition: SRSLid.h:46
float sparsenessasymmslice
Definition: SRNueCosRej.h:175
float meanqfracbkgexp
Fraction matched in background matches, "exp".
Definition: SRLem.h:42
float cospidcontain
Nue cosrej PID. A BDT using 5 variables relating to containment.
Definition: SRNueCosRej.h:115
float efracp2
Definition: SRXnue.h:22
unsigned int noutput
Number of entries in output vector.
Definition: SRCVNResult.h:27
virtual double TotalLength() const
Length (cm) of a shower.
Definition: Shower.cxx:43
float kalfwddist
Kalmantrack projected distance (cm) from end point forwards to det edge.
Definition: SRContain.h:39
float gllt
Gamma ll for ltransverse shower.
Definition: SRSPid.h:27
float annecos
e/cos PID for oscillation analysis
Definition: SRSLid.h:25
void CopyRemidVars(const SRKalmanTrack &srTrk, caf::SRRemid &remid)
Definition: FillPIDs.cxx:712
PID
Definition: FillPIDs.h:14
double ShE() const
Definition: EID.cxx:778
int LastCell(const unsigned int &view) const
Definition: Boundary.h:45
Store BDT variables for the short-baseline oscillation study.
Definition: SRXnue.h:12
float ellt
electron ll for transverse shower
Definition: SRJMEid.h:62
float dedx3
Plane dedx for e/photon/pi0 identification.
Definition: SRSPid.h:41
float angle
Cosine of CosmicTrack with respect to the beam direction.
Definition: SRVeto.h:23
double ANN() const
Definition: EID.cxx:16
float hadE
GeV sum of non (Kalman best ReMId) track E.
Definition: SRCosRej.h:51
void FillNCPi0BkgRejVars(const ncpi0::NCPi0BkgRej &ncpi0bkgrej, caf::SRNCPi0BkgRej &srncpi0bkgrej)
Definition: FillPIDs.cxx:530
float Dedxp4c1() const
Definition: ShowerLID.h:142
float cosmicid
Likelihood Cosmic Interaction.
Definition: SRCVNResult.h:24
double EMuLLT() const
Definition: EID.cxx:668
double fVtxY
True vertex of Y position.
Definition: TrainingData.h:67
float mindist
minimum projected distance to edge of any Kalman track with > 15 hits
Definition: SRCosRej.h:17
float epl20frac
fraction of energy in first 20 planes out of total slice energy
Definition: SRRvp.h:26
A class detailing the cuts made on a particular slice.
Definition: Veto.h:12
unsigned int lastplane
global position of last plane in map
Definition: SRPixelObjMap.h:24
float pid
The output value from the BDT.
Definition: SRRvp.h:19
int cont
True if track is contained as defined by ReMId.
Definition: SRRemid.h:30
float elll
Electron LL for longitudinal shower.
Definition: SRSLid.h:58
double CosPIDDT() const
Real Adaptive BDT trained, based on simple containment variables.
Definition: NCCosRej.h:23
float photonid
Likelihood photon.
float cosbakdist
cosmictrack projected distance (cm) from start point backwards to det edge
Definition: SRContain.h:34
void FillCVNProngTrainingData(const cvn::ProngTrainingData &cvnpdata, caf::SRProngTrainingData &srpdata)
Definition: FillPIDs.cxx:677
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.
Definition: FillTruth.cxx:25
std::vector< float > output
List of net output nodes.
Definition: SRCVNFeatures.h:15
float PLLT() const
Definition: ShowerPID.h:76
int MostMissingPlanes(geo::View_t view) const
Longest run of adjacent planes with no hits.
Definition: Cluster.cxx:668
double StopMinDistToTop() const
Minimum perpendicular distance of all prongs stop points to the top edge of the detector.
Definition: NueCosRej.h:120
float cospidperibdt
Nue cosrej PID for the peripheral sample for 2020+.
Definition: SRNueCosRej.h:126
double vtxx
True vertex of X position.
float electronid
Likelihood electron.
double ANNENoCos() const
Definition: EID.cxx:70
float ProngEBal3D() const
Definition: RVP.h:71
int Prongs2D() const
Definition: RVP.h:72
float pid
PID value output by kNN.
Definition: SRRemid.h:25
bool keep
Definition: SRVeto.h:34
double AngleVarOld() const
Definition: Veto.cxx:95
int kalfwdcellnd
Kalmantrack projected # cells from end point forwards to det edge, including muon catcher...
Definition: SRContain.h:41
float p2Fmip
Definition: SRXnue.h:19
float pi0mass
best pi0 mass hypothesis coming from combinations of JMShowers
Definition: SRSLid.h:49
bool keepnue
Does the event past nue version of veto? (see doc-db 14654)
Definition: SRVeto.h:36
float dedxp3c1
Definition: SRSLid.h:83
void FillLEMVars(const lem::PIDDetails &lem, caf::SRLem &srlem)
Definition: FillPIDs.cxx:316
void FillNuePreselVars(const presel::PreselObj &nuepre, caf::SRPresel &srnuepre)
Definition: FillPIDs.cxx:539
void FillCosRejVars(const cosrej::CosRejObj &cosrej, caf::SRCosRej &srcosrej, caf::SRContain &srcontain)
Definition: FillPIDs.cxx:394
float P2Fmip() const
Definition: Xnue.h:49
float fLepEnergy
True energy of outgoing lepton.
Definition: TrainingData.h:65
void FillXnueVars(const xnue::Xnue &xnue, caf::SRXnue &srxnue)
Definition: FillPIDs.cxx:379
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...
Definition: NueCosRej.h:116
double StopMinDistToEast() const
Minimum perpendicular distance of all prongs stop points to the east edge of the detector (east is ne...
Definition: NueCosRej.h:140
float CScoreDiff() const
Definition: CosRejObj.cxx:520
Result for CVN.
float CosBakAir() const
Definition: CosRejObj.cxx:580
double CalorimetricEnergy(EEnergyCalcScheme escheme=kRecomputeEnergy) const
Simple estimate of neutrino energy.
Definition: Cluster.cxx:439
float fNuEnergy
True energy of neutrino event.
Definition: TrainingData.h:64
if(dump)
float meanqfracfit
Fraction of charge matched, "fit".
Definition: SRLem.h:37
float distallpngfront
Definition: SRNueCosRej.h:99
double StartMinDistToBottom() const
Minimum perpendicular distance of all prongs start points to the bottom edge of the detector...
Definition: NueCosRej.h:100
float anglekal
cos of angle of best ReMId Kalman track
Definition: SRCosRej.h:23
float CosChiDiff() const
Definition: CosRejObj.cxx:749
unsigned int nplanes
number of planes in pixel map
Definition: SRPixelMap.h:18
float efrac10p
Definition: SRXnue.h:21
float dedxp4c0
Definition: SRSLid.h:84
float epi0lll
Electron - Pi0 LL for longitudinal shower.
Definition: SRSLid.h:40
int mincell
minimum projected cells to edge of any Kalman track with > 15 hits
Definition: SRCosRej.h:18
float CosBakDist() const
Definition: CosRejObj.cxx:574
int ismuon
Boolean to distinguish muons from electrons.
Definition: SRJMEid.h:28
std::vector< unsigned char > PixelMapToVector(bool useGeV) const
Definition: PixelMap.cxx:449
float FScattSig() const
Definition: CosRejObj.cxx:785
float dedxp2c1
Definition: SRSLid.h:81
void FillNueCosRejVars(const cosrej::NueCosRej &nuecosrej, caf::SRNueCosRej &srnuecosrej)
Definition: FillPIDs.cxx:458
float MipFrac() const
Definition: RVP.h:64
double StartMinDistToBack() const
Minimum perpendicular distance of all prongs start points to the back edge of the detector...
Definition: NueCosRej.h:108
float dedx5
Definition: SRSLid.h:71
float annEPi0EL
ann output, EPi0 Elastic scattering
Definition: SRJMEid.h:25
float meanysig
Mean hadronic y of matches that are signal.
Definition: SRLem.h:33
unsigned char cvnobjmap[16000]
store 1D array with hit object id
Definition: SRPixelObjMap.h:34
float bakdist
Projected live distance to the edge projected backwards from track start.
Definition: SRVeto.h:20
float PiLLT() const
Definition: ShowerPID.h:80
void FillRVPVars(const rvp::RVP &rvp, caf::SRRvp &srrvp)
Definition: FillPIDs.cxx:357
float kdirscore
chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for best ...
Definition: SRCosRej.h:57
unsigned int firstplane
global position of first plane in map
Definition: SRPixelMap.h:24
bool PassAngle() const
Definition: Veto.cxx:110
bool IsMuon() const
Definition: ShowerLID.h:122
const art::PtrVector< rb::CellHit > & YCells() const
Get all cells from the x-view.
Definition: Cluster.h:126
float fannid
ANN output, "fann" package.
Definition: SRLem.h:50
float VtxDist() const
Definition: CosRejObj.cxx:755
float meanqfracsig
Fraction matched in signal matches.
Definition: SRLem.h:39
unsigned int pdgmax
pdg of largest likelihood
float RVP12() const
Definition: RVP.h:75
double ParticleShowerTransMom() const
Transverse component of the event momentum. Particle masses in the momentum calculation are inferred ...
Definition: NueCosRej.h:32
double DirY() const
Definition: Veto.cxx:80
float ConCosPID() const
Definition: CosRejObj.cxx:465
float FScattSum() const
Definition: CosRejObj.cxx:773
float cosmicid
Likelihood of Cosmic.
Definition: SRSliceLID.h:16
virtual TVector3 Dir() const
Unit vector describing prong direction.
Definition: Prong.h:77
float KDirScore() const
Definition: CosRejObj.cxx:502
TVector3 Stop() const
Endpoint of the shower.
Definition: Shower.cxx:55
float mulll
Muon LL for longitudinal shower.
Definition: SRSLid.h:60
float Dedx1() const
Definition: ShowerPID.h:61
float CosFwdAir() const
Definition: CosRejObj.cxx:562
double fMeanQFracSigExp
Definition: PIDExtraVars.h:44
unsigned int ncomponents
Size of components vector.
Definition: SRCVNFeatures.h:19
float Dedxp2c0() const
Definition: ShowerLID.h:137
int KalBakCell() const
Definition: CosRejObj.cxx:664
double CosPIDContain() const
A BDT trained based on simple containment variables Used in the peripheral sample of the 2017 analysi...
Definition: NueCosRej.h:145
float cospidcorebdt
Nue cosrej PID for the core sample for 2020+.
Definition: SRNueCosRej.h:128
float dedxp1c0
Definition: SRSLid.h:78
float annEL
ann output,
Definition: SRJMEid.h:23
float elll
Electron ll for longitudinal shower.
Definition: SRSPid.h:24
float dedxp2c0
Definition: SRSLid.h:80
float Radius() const
Definition: ShowerLID.h:148
float emaxfrac6p
Definition: SRXnue.h:20
std::vector< float > output
List of net output nodes.
Definition: SRCVNResult.h:19
float epllt
Electron - proton LL for transverse shower.
Definition: SRSLid.h:43
float cosfwddist
cosmictrack projected distance (cm) from end point forwards to det edge
Definition: SRContain.h:29
float measfrac
Fraction of planes used to measure dE/dx.
Definition: SRRemid.h:29
float dedx0
Definition: SRSLid.h:66
int cosfwdcellnd
cosmictrack projected # cells from end point forwards to det edge, including muon catcher...
Definition: SRContain.h:31
CVN PID output for a single particle.
float HadronicEnergy() const
Definition: ShowerLID.h:159
float depE
Total energy deposited in shower [GeV].
Definition: SRSLidEnergy.h:20
float pidexp
Fraction of matches that are signal, "exp".
Definition: SRLem.h:27
float Dedxp5c1() const
Definition: ShowerLID.h:144
float EFrac6PlWin() const
Definition: RVP.h:68
float enlll
Electron - neutron LL for longitudinal shower.
Definition: SRSLid.h:44
float CosThetaTrue() const
Definition: CosRejObj.cxx:724
Definition: RVP.h:16
float EMuLLL() const
Definition: ShowerLID.h:107
double HitsPerPlaneAsymm() const
Asymmetry in the hitsperplane at the start and end of the leading shower.
Definition: NueCosRej.h:199
float energydiff
Potential diff between sig and bkg matches.
Definition: SRLem.h:44
int CosBakCell() const
Definition: CosRejObj.cxx:640
float meanqfracsigexp
Fraction matched in signal matches "exp".
Definition: SRLem.h:40
bool Contained() const
Definition: ReMId.cxx:243
float Dedx5() const
Definition: ShowerPID.h:65
void FillCVNParticleResultVars(const std::vector< art::Ptr< rb::PID > > &cvnparts, caf::SRCVNParticleResult &cvnpart)
Definition: FillPIDs.cxx:112
float epillt
Electron - pion LL for transverse shower.
Definition: SRSLid.h:47
int CosFwdCell() const
Definition: CosRejObj.cxx:628
unsigned int lastcellx
global position of last x cell in map
Definition: SRPixelMap.h:27
float cosfwdair
for cosmictrack projected distance forwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:32
float fFracY
fraction of y view hits contained in map
Definition: PixelMap.h:120
double ENLLT() const
Definition: EID.cxx:728
unsigned int nchan
number of channels in pixel map
Definition: SRPixelMap.h:23
Contains the reco muon PID (ReMId) output.
Definition: SRRemid.h:19
float EndDist() const
Definition: CosRejObj.cxx:761
std::map< int, float > fPartTransLL
Map of the transverse ll by particle type.
Definition: ShowerLID.h:241
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.
Definition: NueCosRej.h:187
unsigned short nplanex
number of planes spanned in x view
Definition: SRShower.h:38
double CosPIDContainXY() const
A BDT trained based on simple containment variables. Includes pX/p and pY/p separately instead of pT/...
Definition: NueCosRej.h:150
float plll
Proton ll for longitudinal shower.
Definition: SRSPid.h:32
float cosfwdsteel
for cosmictrack projected distance forwards, how much is through steel (ND only, currently is just di...
Definition: SRContain.h:33
This class contains the LEM PID output.
Definition: SRLem.h:18
float annE
ann output, with E
Definition: SRJMEid.h:22
float avgsigE
Fraction of signal matches, sigmoid weighting.
Definition: SRLem.h:60
float dedx1
Plane Dedx.
Definition: SRJMEid.h:37
float rawgevslice
event energy in GeV
Definition: SRNuonEResult.h:20
double EPLLL() const
Definition: EID.cxx:698
double StartMinDistToWest() const
Minimum perpendicular distance of all prongs start points to the west edge of the detector (west is p...
Definition: NueCosRej.h:112
float chisig
chi-square of a fit assuming best match is signal
Definition: SRLem.h:53
double StartDistToBottom() const
Perpendicular distance of start point of the leading prong to the bottom edge of the detector...
Definition: NueCosRej.h:52
CVN PID output.
Definition: SRCVNResult.h:13
float shwEFrac
fraction of energy of leading shower out of total energy of slice
Definition: SRJMEid.h:56
double StopMinDistToFront() const
Minimum perpendicular distance of all prongs stop points to the front edge of the detector...
Definition: NueCosRej.h:128
float MuLLT() const
Definition: ShowerPID.h:72
float efrac2plwin
largest fraction of energy in a 2 plane window out of total slice
Definition: SRRvp.h:27
SRVector3D start
Shower start point in detector coordinates. [cm].
Definition: SRShower.h:32
float CDirScore() const
Definition: CosRejObj.cxx:508
float EPi0LLL() const
Definition: ShowerLID.h:109
unsigned int fNCell
Number of cells, width of pixel map.
Definition: PixelMap.h:118
double StopDistToBack() const
Perpendicular distance of stop point of the leading prong to the back edge of the detector...
Definition: NueCosRej.h:84
float fwddist
Projected live distance to the edge projected forwards from track end.
Definition: SRVeto.h:19
float meanysigexp
Hadronic y of matches that are signal, "exp".
Definition: SRLem.h:34
double ProngDistToVtx() const
Perpendicular distance of closes approach between the leading prong and the event vertex...
Definition: NueCosRej.h:168
float nlll
Neutron ll for longitudinal shower.
Definition: SRSPid.h:34
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
Definition: SRPixelObjMap.h:32
float protonid
Likelihood proton.
float dedx4
Plane dedx for e/photon/pi0 identification.
Definition: SRSPid.h:42
double AngleVar() const
Definition: Veto.cxx:90
float kalbaksteel
for Kalmantrack projected distance backwards, how much is through steel (ND only, currently is just d...
Definition: SRContain.h:48
Output from Cosmic Rejection (Nuecosrej) module.
Definition: SRNueCosRej.h:10
double PhotonShowerMomY() const
The cosine between the reconstructed momentum and the y-axis. Particle masses in the momentum calcula...
Definition: NueCosRej.h:44
float hitfracx
faction of x hits in slice in x view pixel map
Definition: SRPixelMap.h:30
float ValueECos() const
Definition: ShowerLID.h:146
double ProngMinZ() const
The minimium Z between the start and stop points of the leading prong in the slice.
Definition: NueCosRej.h:191
float enrichfracexp
Fraction of matches that are "enriched", "exp".
Definition: SRLem.h:48
double StartMinDistToFront() const
Minimum perpendicular distance of all prongs start points to the front edge of the detector...
Definition: NueCosRej.h:104
CVN PID output for a single particle.
int MinCell() const
Definition: CosRejObj.cxx:688
float emulll
Electron - muon ll for longitudinal shower.
Definition: SRJMEid.h:44
float CosFwdSteel() const
Definition: CosRejObj.cxx:568
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)
Definition: SRCosRej.h:63
float pi0lll
Pi0 ll for longitudinal shower.
Definition: SRSPid.h:30
float KScoreDiff() const
Definition: CosRejObj.cxx:514
float hitsperplaneasymmslice
Definition: SRNueCosRej.h:180
void FillNCCosRejVars(const ncid::NCCosRej &nccosrej, caf::SRNCCosRej &srnccosrej)
Definition: FillPIDs.cxx:521
float RecoE() const
Definition: RVP.h:61
float avgexpE
Fraction of signal matches, "exp" weighting.
Definition: SRLem.h:59
float FScattExt() const
Definition: CosRejObj.cxx:779
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)
Definition: FillPIDs.cxx:169
Contains the RVP PID output.
Definition: SRRvp.h:13
int ncell
Number of hits in the slice.
Definition: SRVeto.h:18
bool passnumicut
Does this event pass the NuMI timing cuts?
Definition: SRVeto.h:32
float miphits
number of mip hits in slice
Definition: SRRvp.h:25
float pi0id
Likelihood primary prong is a pi0.
Definition: SRNuonEResult.h:22
float Efracp3() const
Definition: Xnue.h:53
float coschidiff
chisq difference of track timing fits (cosmictrack)
Definition: SRCosRej.h:66
unsigned short nplane
number of planes spanned
Definition: SRShower.h:28
unsigned int NYCell() const
Number of cells in the y-view.
Definition: Cluster.h:108
float coschisq
chisq value of slope fit (cosmictrack)
Definition: SRCosRej.h:65
float EPiLLL() const
Definition: ShowerLID.h:115
float Dedx5() const
Definition: ShowerLID.h:131
float longtrfrac
fraction of cells in longest track out of total cells in slice
Definition: SRRvp.h:23
float emid
Likelihood electron or photon.
float Dedx3() const
Definition: ShowerPID.h:63
double vtxy
True vertex of Y position.
void FillMuIdVars(const remid::ReMId &remid, SRKalmanTrack &srTrk)
Definition: FillPIDs.cxx:37
X or Y views.
Definition: SREnums.h:24
unsigned short nhity
number of hits in y-view
Definition: SRShower.h:27
float diry
Cosine of CosmicTrack with respect to the y-direction.
Definition: SRVeto.h:22
float enllt
Electron - neutron ll for ltransverse shower.
Definition: SRJMEid.h:51
unsigned int lastplane
global position of last plane in map
Definition: SRPixelMap.h:25
double Angle() const
Definition: Veto.cxx:85
double fVtxX
True vertex of X position.
Definition: TrainingData.h:66
float efrac2sig
Definition: SRXnue.h:25
SRParticleTruth truth
Truth information for the prong.
Definition: SRProng.h:46
void FillLIDEnergyVars(const slid::ShowerLID &slid, caf::SRFuzzyKProng &png)
Definition: FillPIDs.cxx:306
PIDExtraVars vars
Definition: PIDDetails.h:32
double ANNE() const
Definition: EID.cxx:25
double StartDistToTop() const
Perpendicular distance of start point of the leading prong to the top edge of the detector...
Definition: NueCosRej.h:48
float Pi0mass() const
Definition: ShowerLID.h:118
Cosmic Rejection PIDs for Numu analysis.
Definition: FillParentInfo.h:9
float NueEnergy() const
Definition: ShowerLID.h:158
float dedxllh
Log-likelihood value from dE/dx.
Definition: SRRemid.h:27
float meanqfracexp
Fraction of charge matched, "exp".
Definition: SRLem.h:38
float P1Energy() const
Definition: Xnue.h:47
virtual void setDefault()
unsigned int firstcelly
global position of first y cell in map
Definition: SRPixelMap.h:28
float epi0nuoneid
Likelihood nu-on-e from epi0 classifier.
Definition: SRNuonEResult.h:25
float len
track length [cm]
Definition: SRTrack.h:41
Output of the jmshower::NueSel module.
Definition: SRJMEid.h:15
float fannidenrich
ANN output, including enrichfracexp.
Definition: SRLem.h:51
float maxval
Maximum likelihood among net outputs.
float eglll
Electron - gamma LL for longitudinal shower.
Definition: SRSLid.h:36
ProngType fLabel3D
Class of the 3D prong.
float muclosestapproachbytime
Definition: SRNueCosRej.h:214
float eiso3sig
fraction of slice energy more then 3sigma from mean energy weighted slice position ...
Definition: SRRvp.h:34
float cosbaksteel
for cosmictrack projected distance backwards, how much is through steel (ND only, currently is just d...
Definition: SRContain.h:38
Result, basic output of CVN neural net.
Definition: Result.h:15
float enlll
Electron - neutron ll for longitudinal shower.
Definition: SRJMEid.h:50
float prongEbal2D
the energy balance between highest energy 2D prong and 3D prong in the opposite view ...
Definition: SRRvp.h:33
float mutimediffbytime
Definition: SRNueCosRej.h:210
float Dedxp3c0() const
Definition: ShowerLID.h:139
float ShowerDepEnergy() const
Definition: ShowerLID.h:156
int kalfwdcell
Kalmantrack projected # cells from end point forwards to det edge.
Definition: SRContain.h:40
float dedx2
Definition: SRSLid.h:68
int label3d
Class of the 3D prong.
float VertexEnergy() const
Definition: ShowerLID.h:155
float anglevar
Defined as: abs(angle)*(diry+1) - new angle variable for new angle cut.
Definition: SRVeto.h:24
float pidexpfit
Fraction of matches that are signal, "expfit".
Definition: SRLem.h:28
float egllt
Electron - gamma LL for transverse shower.
Definition: SRSLid.h:37
Containment variables.
Definition: SRContain.h:12
float fFracX
fraction of x view hits contained in map
Definition: PixelMap.h:119
FinalStateProngsType fFSp
Class of the event by primary n prongs.
Definition: TrainingData.h:62
float enddist
shortest distance to wall from primary Kalman track end position
Definition: SRContain.h:27
float ShowerEnergy() const
Definition: ShowerLID.h:157
float muonid
MuonID classifier value.
Definition: SRKalmanTrack.h:19
double MuLLL() const
Definition: EID.cxx:840
double ProngMaxX() const
The maximum X between the start and stop points of the leading prong in the slice.
Definition: NueCosRej.h:172
float hitsperplaneasymm
Definition: SRNueCosRej.h:170
double ProngTransMom() const
Transverse component of the energy-weighted average of prong directions.
Definition: NueCosRej.h:26
float pi0mass
best pi0 mass hypothesis coming from combinations of JMShowers
Definition: SRJMEid.h:55
Horizontal planes which measure Y.
Definition: SREnums.h:23
float meany
Mean hadronic y of matches.
Definition: SRLem.h:30
float annepi0el
particle PID.
Definition: SRSPid.h:22
float pid
MuonID value output by ND-trained BDT.
Definition: SRMuonID.h:18
float Dedx4() const
Definition: ShowerPID.h:64
float epi0pi0id
Likelihood pi0 from epi0 classifier.
Definition: SRNuonEResult.h:26
float otherid
Likelihood other.
float mulll
Muon ll for longitudinal shower.
Definition: SRSPid.h:28
float pi0llt
Pi0 ll for ltransverse shower.
Definition: SRSPid.h:31
float pizeroid
Likelihood pizero.
float dedx1
Plane dedx for e/photon/pi0 identification.
Definition: SRSPid.h:39
int NTracks3D() const
Definition: CosRejObj.cxx:484
double EPi0LLT() const
Definition: EID.cxx:688
float CosYPosAtTrans() const
Definition: CosRejObj.cxx:676
int nkal3d
number of 3D tracks as determined by Kalman tracker
Definition: SRCosRej.h:25
float CFitSpeed() const
Definition: CosRejObj.cxx:496
double SparsenessAsymm() const
Asymmetry in the zero-hit planes at the start and end of the leading shower.
Definition: NueCosRej.h:195
float GLLT() const
Definition: ShowerPID.h:70
SRParticleTruth truthYView
Truth information for the prong.
Definition: SRProng.h:48
The TrainingData objects contains a PixelMap and the output class type, and any other bit that goes i...
Definition: TrainingData.h:23
A rb::Prong with a length.
Definition: Shower.h:18
unsigned int nchan
number of channels in pixel map
Definition: SRPixelObjMap.h:22
double Ratio() const
Definition: Veto.cxx:75
unsigned int NXCell() const
Number of cells in the x-view.
Definition: Cluster.h:106
float eratio
ratio of best kalman track GeV / slice GeV
Definition: SRCosRej.h:50
double EelLLL() const
Definition: EID.cxx:818
PixelMap, basic input to CVN neural net.
Definition: PixelMap.h:23
float recoE
reco energy of slice
Definition: SRRvp.h:21
int MostContiguousPlanes(geo::View_t view) const
Longest run of adjacent planes with hits.
Definition: Cluster.cxx:635
float ValueEPi0() const
Definition: ShowerLID.h:124
float MipHits() const
Definition: RVP.h:65
float rvp12
PID value for 12 variable RVP for comparison purposes.
Definition: SRRvp.h:35
float PiLLL() const
Definition: ShowerPID.h:79
std::vector< float > output
List of net output nodes.
Definition: SRNuonEResult.h:19
double StartDistToEast() const
Perpendicular distance of start point of the leading prong to the east edge of the detector (east is ...
Definition: NueCosRej.h:68
float numuid
Likelihood Charge Current NuMu.
Definition: SRCVNResult.h:20
Contains the reco muon PID (ReMId) output.
Definition: SRMuonID.h:12
float ValueEPi0EL() const
Definition: ShowerPID.h:59
Variables describing Michel E&#39;s found around the end of a track.
Definition: SRPixelObjMap.h:12
float Dedxp0c0() const
Definition: ShowerLID.h:133
float ProngEBal2D() const
Definition: RVP.h:73
unsigned int ExtentPlane(geo::View_t view=geo::kXorY) const
Definition: Cluster.h:250
float Efrac2sig() const
Definition: Xnue.h:55
double PhotonShowerMomX() const
The cosine between the reconstructed momentum and the x-axis. Particle masses in the momentum calcula...
Definition: NueCosRej.h:40
Definition: FillPIDs.h:16
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
Definition: SRPixelMap.h:33
float MinDist() const
Definition: CosRejObj.cxx:694
float hadronid
Likelihood Proton or Pion.
Vertical planes which measure X.
Definition: SREnums.h:22
bool PassThruOld() const
Definition: Veto.cxx:105
unsigned int NCell() const
Definition: Veto.cxx:60
float InverseGLLL() const
Definition: ShowerLID.h:150
float epi0llt
Electron - Pi0 LL for transverse shower.
Definition: SRSLid.h:41
double StopDistToBottom() const
Perpendicular distance of stop point of the leading prong to the bottom edge of the detector...
Definition: NueCosRej.h:76
float nuenergy
True energy of neutrino event.
int kalbakcell
Kalmantrack projected # cells from start point backwards to det edge.
Definition: SRContain.h:45
float EMuLLT() const
Definition: ShowerLID.h:108
double Dedx0() const
Definition: EID.cxx:594
float ENLLL() const
Definition: ShowerLID.h:113
float kalthetatrue
cosine of angle between kalman trk dir and true dir of most contributing particle ...
Definition: SRCosRej.h:27
float nueccid
Likelihood nuecc events.
Definition: SRNuonEResult.h:23
float GLLL() const
Definition: ShowerPID.h:69
CVN features.
Definition: SRCVNFeatures.h:9
double ProngTransMom() const
Transverse component of the energy-weighted average of prong directions.
Definition: NCCosRej.h:20
float EFrac2SigRd() const
Definition: RVP.h:69
unsigned int firstcellx
global position of first x cell in map
Definition: SRPixelObjMap.h:25
int KalFwdCellND() const
Definition: CosRejObj.cxx:658
float KFitSpeed() const
Definition: CosRejObj.cxx:490
double fVtxZ
True vertex of Z position.
Definition: TrainingData.h:68
float KalFwdAir() const
Definition: CosRejObj.cxx:598
unsigned int lastcelly
global position of last y cell in map
Definition: SRPixelObjMap.h:28
int prongs3D
number of 3D FuzzyK prongs associated with slice
Definition: SRRvp.h:30
Attach LEM-specific info to the base PID object.
Definition: PIDDetails.h:20
unsigned int lastcellx
global position of last x cell in map
Definition: SRPixelObjMap.h:26
float numuunconttunedpid
tuned cosmic rejection PID for uncontained events - Jose&#39;s;
Definition: SRCosRej.h:36
double MeasurementFraction() const
Return the measurement fraction variable used as an input to the kNN that determines the pid value...
Definition: ReMId.cxx:213
float tmvamlp
TMVA ANN ("multi layer perceptron")
Definition: SRLem.h:66
double HitsPerPlane() const
Number of hits in slice per number of planes in leading prong. This is a handle on rejecting FEB flas...
Definition: NueCosRej.h:22
double ELLT() const
Definition: EID.cxx:808
float gap
gap from vertex to start of shower
Definition: SRSLid.h:51
int MuSliceIdxByTime() const
The index of the slice that is most likely to be the parent muon according to difference in time...
Definition: NueCosRej.h:229
double StartDistToBack() const
Perpendicular distance of start point of the leading prong to the back edge of the detector...
Definition: NueCosRej.h:60
float Dedx4() const
Definition: ShowerLID.h:130
void FillRockPreselVars(const presel::PreselObj &rockpre, caf::SRPresel &srrockpre)
Definition: FillPIDs.cxx:547
void FillSpidVars(const slid::ShowerPID &spid, caf::SRSPid &srspid)
Definition: FillPIDs.cxx:276
Build slid::LID objects to store electron ID, if asked for, otherwise, calculate LID info and make av...
Definition: FillPIDs.h:13
SRVector3D stop
shower stop point
Definition: SRShower.h:41
float mullt
Muon LL for transverse shower.
Definition: SRSLid.h:61
float dedx3
Definition: SRSLid.h:69
float FwdDist() const
Definition: Veto.cxx:65
double ProngMaxZ() const
The maximum Z between the start and stop points of the leading prong in the slice.
Definition: NueCosRej.h:180
int CosFwdCellND() const
Definition: CosRejObj.cxx:634
float Dedxp0c1() const
Definition: ShowerLID.h:134
double EPLLT() const
Definition: EID.cxx:708
void FillCVNPixelMaps(const cvn::PixelMap &map, caf::SRPixelObjMap &srcvnmap, bool useGeV)
Definition: FillPIDs.cxx:617
float epi0lll
Electron - Pi0 ll for longitudinal shower.
Definition: SRJMEid.h:46
unsigned int finalstate
Class of the event by final state.
float eelllt
electron el ll for transverse shower
Definition: SRJMEid.h:64
void FillCVNTrainingData(const cvn::TrainingData &cvntdata, caf::SRTrainingData &srtdata)
Definition: FillPIDs.cxx:661
Output from Cosmic Rejection (Nuecosrej) module.
Definition: SRNCCosRej.h:10
float ERatio() const
Definition: CosRejObj.cxx:700
A PID for muons.
Definition: FillPIDs.h:10
float longtr
longest discrete track merge track
Definition: SRRvp.h:22
float dedx0
Definition: SRJMEid.h:36
int MuSliceIdxByDist() const
The index of the slice that is most likely to be the parent muon according to distance of closest app...
Definition: NueCosRej.h:212
This module creates Common Analysis Files.
Definition: FileReducer.h:10
double StartDistToFront() const
Perpendicular distance of start point of the leading prong to the front edge of the detector...
Definition: NueCosRej.h:56
ParentParticleType fPPt
Class of the event by parent particle.
Definition: TrainingData.h:59
float nutauid
Likelihood of Charge Current NuTau.
Definition: SRSliceLID.h:15
float efrac2sigrd
fraction of slice energy within 2sigma of mean energy weighted slice position
Definition: SRRvp.h:29
SRSLidEnergy lidE
Definition: SRShowerLID.h:34
double EPi0LLL() const
Definition: EID.cxx:678
float ncid
Likelihood of Neutral Current.
Definition: SRSliceLID.h:12
float numucontpid2019
cosmic rejection PID for contained events; 2019 Analysis
Definition: SRCosRej.h:35
void FillCVNResultVars(const cvn::Result &result, caf::SRCVNResult &srcvn, bool NoCosmics)
Definition: FillPIDs.cxx:595
double StopMinDistToBack() const
Minimum perpendicular distance of all prongs stop points to the back edge of the detector.
Definition: NueCosRej.h:132
Boundary fBound
Definition: PixelMap.h:153
T min(const caf::Proxy< T > &a, T b)
float kalbakdist
Kalmantrack projected distance (cm) from start point backwards to det edge.
Definition: SRContain.h:44
float numuid
Likelihood of Charge Current NuMu.
Definition: SRSliceLID.h:13
void FillNuonEResultVars(const cvn::Result &result, caf::SRNuonEResult &srnuone)
Definition: FillPIDs.cxx:576
float meanqfrac
Mean fraction of charge matched.
Definition: SRLem.h:36
float cosbakair
for cosmictrack projected distance backwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:37
int Prongs3D() const
Definition: RVP.h:70
float NLLT() const
Definition: ShowerPID.h:78
double MuTimeDiffByDist() const
Minimum time difference between this slice and the possible muon parent (selected by closest approach...
Definition: NueCosRej.h:216
double MuTimeDiffByTime() const
Minimum time difference between this slice and the possible muon parent (selected by difference in ti...
Definition: NueCosRej.h:233
Object collecting Preselection variables.
float numucontpid2020
cosmic rejection PID for contained events; 2020 Analysis
Definition: SRCosRej.h:34
float kfitspeed
fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for best ReMId Kalman track ...
Definition: SRCosRej.h:55
float dedx2
Plane Dedx.
Definition: SRJMEid.h:38
float Gap() const
Definition: ShowerLID.h:117
float EGLLT() const
Definition: ShowerLID.h:106
double StopDistToEast() const
Perpendicular distance of stop point of the leading prong to the east edge of the detector (east is n...
Definition: NueCosRej.h:92
float ratio
Number of hits in the track relative to the slice.
Definition: SRVeto.h:21
double Dedx5() const
Definition: EID.cxx:628
float anglecos
cos of angle of Cosmic Track
Definition: SRCosRej.h:24
unsigned int argmax
Index of maximum value.
Definition: SRCVNResult.h:26
bool PassNuMICut() const
Definition: Veto.cxx:120
float avginvE
Fraction of signal matches, 1/E weighting.
Definition: SRLem.h:58
float len
track length [cm]
Definition: SRShower.h:34
float ELLL() const
Definition: ShowerPID.h:67
float Dedx1() const
Definition: ShowerLID.h:127
unsigned short maxplanegap
maximum number of gapped planes in prong
Definition: SRShower.h:30
float chisigexp
chisq of a fit assuming best match is signal, "exp"
Definition: SRLem.h:55
float pidfit
Fraction of matches that are signal, "fit".
Definition: SRLem.h:26
int cosbakcellnd
cosmictrack projected # cells from start point backwards to det edge, including muon catcher...
Definition: SRContain.h:36
float chibkg
chisq of a fit assuming best match is background
Definition: SRLem.h:54
float Pi0LLT() const
Definition: ShowerPID.h:74
Definition: fwd.h:28
int prongs2D
number of 2D FuzzyK prongs associated with slice
Definition: SRRvp.h:32
double SparsenessAsymmSlice() const
Asymmetry in the zero-hit planes at the start and end of the slice.
Definition: NueCosRej.h:203
float egllt
Electron - gamma ll for ltransverse shower.
Definition: SRJMEid.h:43
unsigned int noutput
Number of entries in output vector.
Definition: SRCVNFeatures.h:16
float width
Shower width [cm].
Definition: SRShower.h:37
View_t view
Prong view caf::kX = 0, caf::kY = 1 or caf::kXorY = 2.
Definition: SRShower.h:35
unsigned int firstcelly
global position of first y cell in map
Definition: SRPixelObjMap.h:27
float nuoneid
Likelihood nu-on-e.
Definition: SRNuonEResult.h:21
float KalChiDiff() const
Definition: CosRejObj.cxx:743
Definition: FillPIDs.h:15
float mullt
muon ll for transverse shower
Definition: SRJMEid.h:66
float kscorediff
chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit modu...
Definition: SRCosRej.h:59
float fscattext
scattering variable (Fernanda)
Definition: SRCosRej.h:44
double DedxSeparation() const
Return the dE/dx separation variable used as an input to the kNN that determines the pid value...
Definition: ReMId.cxx:199
float OldCosPID() const
Definition: CosRejObj.cxx:471
float pillt
Pion ll for ltransverse shower.
Definition: SRSPid.h:37
float dedx4
Definition: SRSLid.h:70
float dectree
Decision tree PID.
Definition: SRLem.h:68
float dedx5
Plane Dedx.
Definition: SRJMEid.h:41
double Pi0Mass() const
Definition: EID.cxx:768
float annEPi0
ann output, EPi0
Definition: SRJMEid.h:24
double ProngMinX() const
The minimium X between the start and stop points of the leading prong in the slice.
Definition: NueCosRej.h:183
double VtxGeV() const
Definition: EID.cxx:758
float cdirscore
chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for cosmi...
Definition: SRCosRej.h:58
float vtxgev
Energy of slice in vertex region.
Definition: SRSLid.h:48
float hitfracy
fraction of y hits in slice in y view pixel map
Definition: SRPixelObjMap.h:30
int FirstPlane() const
Definition: Boundary.h:42
FinalStateType fFSt
Class of the event by final state.
Definition: TrainingData.h:61
float fscattmax
maximum scattering angle (Fernanda)
Definition: SRCosRej.h:42
float efracp3
Definition: SRXnue.h:23
float otherid
Likelihood other.
Definition: SRNuonEResult.h:24
unsigned int finalstateprong
float ShowerEFrac() const
Definition: ShowerLID.h:120
float EPLLT() const
Definition: ShowerLID.h:112
double ANNEPi0() const
Definition: EID.cxx:43
float Dedxp1c1() const
Definition: ShowerLID.h:136
std::vector< float > fPrincipalComponents
Linear combinations of fOutput.
Definition: CVNFeatures.h:25
float anglevarold
Difference of squares between beam angle and Y-direction (Also labeled CosmicVar AnglePID and AngleVa...
Definition: SRVeto.h:25
float mulll
muon ll for longitudinal shower
Definition: SRJMEid.h:65
double EGLLT() const
Definition: EID.cxx:648
double StartDistToWest() const
Perpendicular distance of start point of the leading prong to the west edge of the detector (west is ...
Definition: NueCosRej.h:64
float eellll
electron el ll for longitudinal shower
Definition: SRJMEid.h:63
float lepenergy
True energy of outgoing lepton.
float hitfracx
faction of x hits in slice in x view pixel map
Definition: SRPixelObjMap.h:29
float tmvabdtg
TMVA boosted decision tree, gradient boost.
Definition: SRLem.h:62
float AngleKal() const
Definition: CosRejObj.cxx:453
unsigned int firstplane
global position of first plane in map
Definition: SRPixelObjMap.h:23
float kalbakair
for Kalmantrack projected distance backwards how much is through air (ND only, NYI) ...
Definition: SRContain.h:47