readNtuple.h
Go to the documentation of this file.
1 #ifndef readNtuple_h
2 #define readNtuple_h
3 
4 #include <TROOT.h>
5 #include <TChain.h>
6 #include <TFile.h>
7 #include "makeHists.h"
8 #include "iostream"
9 #include "fstream"
10 
11 #include "TMVA/Tools.h"
12 #include "TMVA/Reader.h"
13 #include "TMVA/MethodCuts.h"
14 
15 using namespace std;
16 
17 class readNtuple {
18  public :
19 
20  TChain *fChain; //!poInt_ter to the analyzed TChain or TChain
21  Int_t fCurrent; //!current Tree number in a TChain
22 
23 
24  Int_t sl_runnum;
25  Int_t sl_subrunnum;
26  Int_t sl_evtnum;
27  Int_t sl_goodSpill;
28  Int_t sl_Nslices;
29  Double_t sl_spillPOT;
30  Int_t sliceIndex;
31 
32  Float_t weight;
33 
34  //spill quality
35  Int_t sl_Ndeaddcm;
36 
37  //beam quality
38  Double_t sl_horncurrent;
39  Double_t sl_posX;
40  Double_t sl_posY;
41  Double_t sl_widthX;
42  Double_t sl_widthY;
44 
45  /*** Vertex ***/
46  Float_t Vx;//true vertex X
47  Float_t Vy;//true vertex Y
48  Float_t Vz;//true vertex Z
49  Float_t recoVx;//reco vertex X
50  Float_t recoVy;//reco vertex Y
51  Float_t recoVz;//reco vertex Z
52 
53  /*** true info ***/
54  Float_t truePt;//true Pt
55  Float_t trueE;//true E
56  Float_t trueTheta;//true theta
57  Int_t ccnc;
58  Int_t PDG;
59  Int_t origPDG;
60  Int_t mode;
61  bool isQE;
62  Int_t intType;
63  Int_t hitNucl;
64  Float_t HadX;
65  Float_t HadY;
66  Float_t HadW;
67  Float_t qsqr;
68  Float_t OscWt;
69 
70  Int_t lepPDG;
71  Float_t lepPx;
72  Float_t lepPy;
73  Float_t lepPz;
74  Float_t lepE;
75  Int_t mesonPDG;
76  Float_t mesonPx;
77  Float_t mesonPy;
78  Float_t mesonPz;
79 
80  //for pi0->2gamma from most energetic pi0
81  Int_t HasPi0;
82  Int_t HasP2GG;//has pi0->2gamma
83  Int_t PMPDG;//pdg code for pi0 mother
84  Int_t Ndaus;//# of pi0 daughters
85  Int_t NdausPDG[5];//pdg code for pi0 daughters
86  Float_t Epi0;
87  Float_t gam1E;
88  Float_t gam1Px;
89  Float_t gam1Py;
90  Float_t gam1Pz;
91  Float_t gam2E;
92  Float_t gam2Px;
93  Float_t gam2Py;
94  Float_t gam2Pz;
95 
96  //all pi0, first 20 if possible ?
97  Int_t Npi0;
98  Float_t pi0Px[20];
99  Float_t pi0Py[20];
100  Float_t pi0Pz[20];
101  Float_t pi0E[20];
102  Float_t pi0Vx[20];
103  Float_t pi0Vy[20];
104  Float_t pi0Vz[20];
105 
106  /*** slicer information ***/
107  Int_t Ncells;//# of associate cells
108  Int_t Nplanes;//# of extended planes
109  Float_t recoE;//reco energy
110  Int_t MinPlane;//min plane index
111 
112  /*** longest kalman track ***/
113  Float_t KtrackLength;//length of longest track
114  Int_t KtrackNcells;//# of cells associated to longest track
115  Float_t KtrackEnergy;//energy associated to longest track
116  Float_t KtrackNcells_ratio;//# of cells associated to longest track over Ncells
117  Float_t KtrackEnergy_ratio;//energy associated to longest track over recoE
118  Float_t KtrackNplanes;
119  Float_t KtrackTheta;
120  Float_t KtrackPhi;
121  Float_t KtrackStartX;
122  Float_t KtrackStartY;
123  Float_t KtrackStartZ;
124  Float_t KtrackStopX;
125  Float_t KtrackStopY;
126  Float_t KtrackStopZ;
127  Int_t NKtracks;//# of tracks
128 
130  Float_t trackcellX[500];
131  Float_t trackcellY[500];
132  Float_t trackcellZ[500];
133  Float_t trackcellE[500];
134  Int_t trackcellADC[500];
135  Int_t trackcellPlane[500];
136  Int_t trackcellCell[500];
137  Float_t trackcellPECorr[500];
138  Float_t trackcellW[500];
139 
140  /*** mip information ***/
141  Float_t Nmip;//# of mip cells
142  Float_t Fmip;//# of mip cells over Ncells
143 
144  /***slice associated cells information ***/
145  Int_t slcellNcell;
146  Float_t slcellX[500];
147  Float_t slcellY[500];
148  Float_t slcellZ[500];
149  Float_t slcellE[500];
150  Int_t slcellADC[500];
151  Int_t slcellPlane[500];
152  Int_t slcellCell[500];
153  Float_t slcellPECorr[500];
154  Float_t slcellW[500];
155 
156  /*** prong information ***/
157  Int_t Nprongs;
158  Float_t Ebalance;
159  Float_t Dphi;
160 
161  Float_t prongEnergy[20];
162  Float_t prongPhi[20];
163  Float_t prongTheta[20];
164  Float_t prongStartX[20];
165  Float_t prongStartY[20];
166  Float_t prongStartZ[20];
167  Float_t prongLength[20];
168  Float_t prongStopX[20];
169  Float_t prongStopY[20];
170  Float_t prongStopZ[20];
171  Float_t prongEnergyXView[20];//energy for XView cells
172  Float_t prongEnergyYView[20];//energy for YView cells
173  Int_t prongNCells[20];//# of cells
174  Int_t prongNCellsXView[20];//# of cells for XView
175  Int_t prongNCellsYView[20];//# of cells for YView
176  Int_t prongNplanes[20];//# of extent planes
177 
178  //all prong associated cells
180  Int_t allprongIndex[500];
181  Float_t allprongcellX[500];
182  Float_t allprongcellY[500];
183  Float_t allprongcellZ[500];
184  Float_t allprongcellE[500];
185  Int_t allprongcellADC[500];
186  Int_t allprongcellPlane[500];
187  Int_t allprongcellCell[500];
188  Float_t allprongcellPECorr[500];
189  Float_t allprongcellW[500];
190 
191  //leading prong
193  Float_t prong1Fmip;
194  Float_t prong1cellX[500];
195  Float_t prong1cellY[500];
196  Float_t prong1cellZ[500];
197  Float_t prong1cellE[500];
198  Int_t prong1cellADC[500];
199  Int_t prong1cellPlane[500];
200  Int_t prong1cellCell[500];
201  Float_t prong1cellPECorr[500];
202  Float_t prong1cellW[500];
203 
204  //2D prongs
205  Int_t Nprongs2D;
206  Int_t prong2DInXView[20];//prong in X View or not
207  Float_t prongEnergy2D[20];//energy
208  Int_t prongNCells2D[20];//# of cells
209  Float_t Ebalance2D;
210  Float_t prongLength2D[20];
211  Float_t prongTheta2D[20];
212  Float_t prongStartX2D[20];
213  Float_t prongStartY2D[20];
214  Float_t prongStartZ2D[20];
215  Float_t prongStopX2D[20];
216  Float_t prongStopY2D[20];
217  Float_t prongStopZ2D[20];
218 
220  Int_t allprong2DIndex[500];
221  Float_t allprong2DcellX[500];
222  Float_t allprong2DcellY[500];
223  Float_t allprong2DcellZ[500];
224  Float_t allprong2DcellE[500];
225  Int_t allprong2DcellADC[500];
226  Int_t allprong2DcellPlane[500];
227  Int_t allprong2DcellCell[500];
228  Float_t allprong2DcellPECorr[500];
229 
230  Float_t rvp;
231  Float_t lid;
232  Float_t lide;
233  Float_t lem;
234  Int_t Nshowers;
235  Float_t shwLID[20];
236  Int_t shw3D[20];
237  Int_t shwNCells[20];
238  Int_t shwNCellsXview[20];
239  Int_t shwNCellsYview[20];
240  Float_t shwLength[20];
241  Float_t shwStartX[20];
242  Float_t shwStartY[20];
243  Float_t shwStartZ[20];
244  Float_t shwStopX[20];
245  Float_t shwStopY[20];
246  Float_t shwStopZ[20];
247  Float_t shwTheta[20];
248  Float_t shwEcells[20];
249  Float_t shwEnu[20];
250 
252  Int_t shwIndex[500];
253  Float_t shwcellX[500];
254  Float_t shwcellY[500];
255  Float_t shwcellZ[500];
256  Float_t shwcellE[500];
257  Int_t shwcellADC[500];
258  Int_t shwcellPlane[500];
259  Int_t shwcellCell[500];
260  Float_t shwcellPECorr[500];
261  Float_t shwcellW[500];
262 
263  readNtuple(const char* TypeNameIn, const char* TypeNameOut, TChain *tree=0);
264  virtual ~readNtuple();
265  virtual void Init(TChain *tree);
266  virtual void Loop(const char* TypeNameOut, const char* TypeOpt);
267  makeHists* myhists;
268 };
269 
270 #endif
271 
272 #ifdef readNtuple_cxx
273 readNtuple::readNtuple(const char* TypeNameIn, const char* TypeNameOut, TChain *tree)
274 {
275  //tree = new TChain("NCAna/NCAna");
276  tree = new TChain("Xnue");
277  ifstream fin(TypeNameIn,ios::in);
278  TString fstr;
279  while(fin>>fstr) tree->Add(fstr);
280 
281  Init(tree);
282  myhists = new makeHists();
283  TString histoname = (TString)TypeNameOut + "_results.root";
284  myhists->bookHists(histoname);
285 
286  cout<<"********************************************"<<endl;
287  cout<<"Input File is "<<TypeNameIn<<endl;
288  cout<<"Output File is "<<histoname<<endl;
289  cout<<"********************************************"<<endl;
290 }
291 
293 {
294  if (!fChain) return;
295  delete fChain->GetCurrentFile();
296 
297  if(myhists){
298  myhists->saveHists();
299  delete myhists;
300 
301  }
302 }
303 
304 void readNtuple::Init(TChain *tree)
305 {
306  if (!tree) return;
307  fChain = tree;
308 
309  fChain->SetBranchAddress("sl_runnum",&sl_runnum);
310  fChain->SetBranchAddress("sl_subrunnum",&sl_subrunnum);
311  fChain->SetBranchAddress("sl_evtnum",&sl_evtnum);
312  fChain->SetBranchAddress("sl_Nslices",&sl_Nslices);
313  fChain->SetBranchAddress("sl_goodSpill",&sl_goodSpill);
314  fChain->SetBranchAddress("sl_spillPOT",&sl_spillPOT);
315  fChain->SetBranchAddress("sliceIndex",&sliceIndex);
316 
317  fChain->SetBranchAddress("sl_Ndeaddcm",&sl_Ndeaddcm);
318 
319  fChain->SetBranchAddress("sl_horncurrent",&sl_horncurrent);
320  fChain->SetBranchAddress("sl_posX",&sl_posX);
321  fChain->SetBranchAddress("sl_posY",&sl_posY);
322  fChain->SetBranchAddress("sl_widthX",&sl_widthX);
323  fChain->SetBranchAddress("sl_widthY",&sl_widthY);
324  fChain->SetBranchAddress("sl_deltaspilltimensec",&sl_deltaspilltimensec);
325 
326  fChain->SetBranchAddress("Vx",&Vx);
327  fChain->SetBranchAddress("Vy",&Vy);
328  fChain->SetBranchAddress("Vz",&Vz);
329  fChain->SetBranchAddress("recoVx",&recoVx);
330  fChain->SetBranchAddress("recoVy",&recoVy);
331  fChain->SetBranchAddress("recoVz",&recoVz);
332 
333  fChain->SetBranchAddress("truePt",&truePt);
334  fChain->SetBranchAddress("trueE",&trueE);
335  fChain->SetBranchAddress("trueTheta",&trueTheta);
336  fChain->SetBranchAddress("ccnc",&ccnc);
337  fChain->SetBranchAddress("PDG",&PDG);
338  fChain->SetBranchAddress("origPDG",&origPDG);
339  fChain->SetBranchAddress("mode",&mode);
340  fChain->SetBranchAddress("isQE",&isQE);
341  fChain->SetBranchAddress("intType",&intType);
342  fChain->SetBranchAddress("hitNucl",&hitNucl);
343  fChain->SetBranchAddress("HadX",&HadX);
344  fChain->SetBranchAddress("HadY",&HadY);
345  //fChain->SetBranchAddress("HadW",&HadW);
346  fChain->SetBranchAddress("qsqr",&qsqr);
347  fChain->SetBranchAddress("OscWt",&OscWt);
348 
349  fChain->SetBranchAddress("lepPDG",&lepPDG);
350  fChain->SetBranchAddress("lepE",&lepE);
351  fChain->SetBranchAddress("lepPx",&lepPx);
352  fChain->SetBranchAddress("lepPy",&lepPy);
353  fChain->SetBranchAddress("lepPz",&lepPz);
354  fChain->SetBranchAddress("mesonPDG",&mesonPDG);
355  fChain->SetBranchAddress("mesonPx",&mesonPx);
356  fChain->SetBranchAddress("mesonPy",&mesonPy);
357  fChain->SetBranchAddress("mesonPz",&mesonPz);
358 
359  fChain->SetBranchAddress("Ncells",&Ncells);
360  fChain->SetBranchAddress("Nplanes",&Nplanes);
361  fChain->SetBranchAddress("recoE",&recoE);
362  fChain->SetBranchAddress("MinPlane",&MinPlane);
363 
364  fChain->SetBranchAddress("slcellNcell",&slcellNcell);
365  fChain->SetBranchAddress("slcellX",&slcellX);
366  fChain->SetBranchAddress("slcellY",&slcellY);
367  fChain->SetBranchAddress("slcellZ",&slcellZ);
368  fChain->SetBranchAddress("slcellE",&slcellE);
369  fChain->SetBranchAddress("slcellADC",&slcellADC);
370  fChain->SetBranchAddress("slcellPlane",&slcellPlane);
371  fChain->SetBranchAddress("slcellCell",&slcellCell);
372  fChain->SetBranchAddress("slcellPECorr",&slcellPECorr);
373  fChain->SetBranchAddress("slcellW",&slcellW);
374 
375  fChain->SetBranchAddress("KtrackLength",&KtrackLength);
376  fChain->SetBranchAddress("KtrackEnergy",&KtrackEnergy);
377  fChain->SetBranchAddress("KtrackNcells_ratio",&KtrackNcells_ratio);
378  fChain->SetBranchAddress("KtrackNcells",&KtrackNcells);
379  fChain->SetBranchAddress("KtrackEnergy_ratio",&KtrackEnergy_ratio);
380  fChain->SetBranchAddress("KtrackNplanes",&KtrackNplanes);
381  fChain->SetBranchAddress("KtrackTheta",&KtrackTheta);
382  fChain->SetBranchAddress("KtrackPhi",&KtrackPhi);
383  fChain->SetBranchAddress("KtrackStartX",&KtrackStartX);
384  fChain->SetBranchAddress("KtrackStartY",&KtrackStartY);
385  fChain->SetBranchAddress("KtrackStartZ",&KtrackStartZ);
386  fChain->SetBranchAddress("KtrackStopX",&KtrackStopX);
387  fChain->SetBranchAddress("KtrackStopY",&KtrackStopY);
388  fChain->SetBranchAddress("KtrackStopZ",&KtrackStopZ);
389  fChain->SetBranchAddress("NKtracks",&NKtracks);
390 
391  fChain->SetBranchAddress("trackcellNcell",&trackcellNcell);
392  fChain->SetBranchAddress("trackcellX",&trackcellX);
393  fChain->SetBranchAddress("trackcellY",&trackcellY);
394  fChain->SetBranchAddress("trackcellZ",&trackcellZ);
395  fChain->SetBranchAddress("trackcellE",&trackcellE);
396  fChain->SetBranchAddress("trackcellADC",&trackcellADC);
397  fChain->SetBranchAddress("trackcellPlane",&trackcellPlane);
398  fChain->SetBranchAddress("trackcellCell",&trackcellCell);
399  fChain->SetBranchAddress("trackcellPECorr",&trackcellPECorr);
400  fChain->SetBranchAddress("trackcellW",&trackcellW);
401 
402  fChain->SetBranchAddress("Nmip",&Nmip);
403  fChain->SetBranchAddress("Fmip",&Fmip);
404 
405  fChain->SetBranchAddress("Nprongs",&Nprongs);
406  fChain->SetBranchAddress("Ebalance",&Ebalance);
407  fChain->SetBranchAddress("Dphi",&Dphi);
408 
409  fChain->SetBranchAddress("prongEnergy",&prongEnergy);
410  fChain->SetBranchAddress("prongPhi",&prongPhi);
411  fChain->SetBranchAddress("prongTheta",&prongTheta);
412  fChain->SetBranchAddress("prongLength",&prongLength);
413  fChain->SetBranchAddress("prongStartX",&prongStartX);
414  fChain->SetBranchAddress("prongStartY",&prongStartY);
415  fChain->SetBranchAddress("prongStartZ",&prongStartZ);
416  fChain->SetBranchAddress("prongStopX",&prongStopX);
417  fChain->SetBranchAddress("prongStopY",&prongStopY);
418  fChain->SetBranchAddress("prongStopZ",&prongStopZ);
419  fChain->SetBranchAddress("prongEnergyXView",&prongEnergyXView);
420  fChain->SetBranchAddress("prongEnergyYView",&prongEnergyYView);
421  fChain->SetBranchAddress("prongNCells",&prongNCells);
422  fChain->SetBranchAddress("prongNCellsXView",&prongNCellsXView);
423  fChain->SetBranchAddress("prongNCellsYView",&prongNCellsYView);
424  fChain->SetBranchAddress("prongNplanes",&prongNplanes);
425 
426  fChain->SetBranchAddress("allprongcellNcell",&allprongcellNcell);
427  fChain->SetBranchAddress("allprongIndex",&allprongIndex);
428  fChain->SetBranchAddress("allprongcellX",&allprongcellX);
429  fChain->SetBranchAddress("allprongcellY",&allprongcellY);
430  fChain->SetBranchAddress("allprongcellZ",&allprongcellZ);
431  fChain->SetBranchAddress("allprongcellE",&allprongcellE);
432  fChain->SetBranchAddress("allprongcellADC",&allprongcellADC);
433  fChain->SetBranchAddress("allprongcellPlane",&allprongcellPlane);
434  fChain->SetBranchAddress("allprongcellCell",&allprongcellCell);
435  fChain->SetBranchAddress("allprongcellPECorr",&allprongcellPECorr);
436  fChain->SetBranchAddress("allprongcellW",&allprongcellW);
437 
438  fChain->SetBranchAddress("prong1cellNcell",&prong1cellNcell);
439  fChain->SetBranchAddress("prong1Fmip",&prong1Fmip);
440  fChain->SetBranchAddress("prong1cellX",&prong1cellX);
441  fChain->SetBranchAddress("prong1cellY",&prong1cellY);
442  fChain->SetBranchAddress("prong1cellZ",&prong1cellZ);
443  fChain->SetBranchAddress("prong1cellE",&prong1cellE);
444  fChain->SetBranchAddress("prong1cellADC",&prong1cellADC);
445  fChain->SetBranchAddress("prong1cellPlane",&prong1cellPlane);
446  fChain->SetBranchAddress("prong1cellCell",&prong1cellCell);
447  fChain->SetBranchAddress("prong1cellPECorr",&prong1cellPECorr);
448  fChain->SetBranchAddress("prong1cellW",&prong1cellW);
449 
450  fChain->SetBranchAddress("Nprongs2D",&Nprongs2D);
451  fChain->SetBranchAddress("Ebalance2D",&Ebalance2D);
452  fChain->SetBranchAddress("prong2DInXView",&prong2DInXView);
453  fChain->SetBranchAddress("prongEnergy2D",&prongEnergy2D);
454  fChain->SetBranchAddress("prongNCells2D",&prongNCells2D);
455  fChain->SetBranchAddress("prongLength2D",&prongLength2D);
456  fChain->SetBranchAddress("prongTheta2D",&prongTheta2D);
457  fChain->SetBranchAddress("prongStartX2D",&prongStartX2D);
458  fChain->SetBranchAddress("prongStartY2D",&prongStartY2D);
459  fChain->SetBranchAddress("prongStartZ2D",&prongStartZ2D);
460  fChain->SetBranchAddress("prongStopX2D",&prongStopX2D);
461  fChain->SetBranchAddress("prongStopY2D",&prongStopY2D);
462  fChain->SetBranchAddress("prongStopZ2D",&prongStopZ2D);
463 
464  fChain->SetBranchAddress("allprong2DcellNcell",&allprong2DcellNcell);
465  fChain->SetBranchAddress("allprong2DIndex",&allprong2DIndex);
466  fChain->SetBranchAddress("allprong2DcellX",&allprong2DcellX);
467  fChain->SetBranchAddress("allprong2DcellY",&allprong2DcellY);
468  fChain->SetBranchAddress("allprong2DcellZ",&allprong2DcellZ);
469  fChain->SetBranchAddress("allprong2DcellE",&allprong2DcellE);
470  fChain->SetBranchAddress("allprong2DcellADC",&allprong2DcellADC);
471  fChain->SetBranchAddress("allprong2DcellPlane",&allprong2DcellPlane);
472  fChain->SetBranchAddress("allprong2DcellCell",&allprong2DcellCell);
473  fChain->SetBranchAddress("allprong2DcellPECorr",&allprong2DcellPECorr);
474 
475  fChain->SetBranchAddress("lem",&lem);
476  fChain->SetBranchAddress("lid",&lid);
477  fChain->SetBranchAddress("lide",&lide);
478  fChain->SetBranchAddress("rvp",&rvp);
479 
480  fChain->SetBranchAddress("Nshowers",&Nshowers);
481  fChain->SetBranchAddress("shwLID",&shwLID);
482  fChain->SetBranchAddress("shw3D",&shw3D);
483  fChain->SetBranchAddress("shwNCells",&shwNCells);
484  fChain->SetBranchAddress("shwNCellsXview",&shwNCellsXview);
485  fChain->SetBranchAddress("shwNCellsYview",&shwNCellsYview);
486  fChain->SetBranchAddress("shwLength",&shwLength);
487  fChain->SetBranchAddress("shwStartX",&shwStartX);
488  fChain->SetBranchAddress("shwStartY",&shwStartY);
489  fChain->SetBranchAddress("shwStartZ",&shwStartZ);
490  fChain->SetBranchAddress("shwStopX",&shwStopX);
491  fChain->SetBranchAddress("shwStopY",&shwStopY);
492  fChain->SetBranchAddress("shwStopZ",&shwStopZ);
493  fChain->SetBranchAddress("shwTheta",&shwTheta);
494  fChain->SetBranchAddress("shwEcells",&shwEcells);
495  fChain->SetBranchAddress("shwEnu",&shwEnu);
496 
497  fChain->SetBranchAddress("shwcellNcell",&shwcellNcell);
498  fChain->SetBranchAddress("shwIndex",&shwIndex);
499  fChain->SetBranchAddress("shwcellX",&shwcellX);
500  fChain->SetBranchAddress("shwcellY",&shwcellY);
501  fChain->SetBranchAddress("shwcellZ",&shwcellZ);
502  fChain->SetBranchAddress("shwcellE",&shwcellE);
503  fChain->SetBranchAddress("shwcellADC",&shwcellADC);
504  fChain->SetBranchAddress("shwcellPlane",&shwcellPlane);
505  fChain->SetBranchAddress("shwcellCell",&shwcellCell);
506  fChain->SetBranchAddress("shwcellPECorr",&shwcellPECorr);
507  fChain->SetBranchAddress("shwcellW",&shwcellW);
508 
509  fChain->SetBranchAddress("weight",&weight);
510 }
511 #endif // #ifdef readNtuple_cxx
512 
Double_t sl_posY
Definition: readNtuple.h:40
Float_t KtrackNcells_ratio
Definition: readNtuple.h:116
Double_t sl_spillPOT
Definition: readNtuple.h:29
TString fin
Definition: Style.C:24
Float_t recoVy
Definition: readNtuple.h:50
virtual void Init(TChain *tree)
Float_t Vx
Definition: readNtuple.h:46
Float_t Vz
Definition: readNtuple.h:48
Float_t Vy
Definition: readNtuple.h:47
Float_t KtrackStopX
Definition: readNtuple.h:124
Float_t Epi0
Definition: readNtuple.h:86
Int_t trackcellNcell
Definition: readNtuple.h:129
Float_t HadX
Definition: readNtuple.h:64
Float_t prong1Fmip
Definition: readNtuple.h:193
Float_t lepPx
Definition: readNtuple.h:71
Float_t Fmip
Definition: readNtuple.h:142
const Var weight
Int_t hitNucl
Definition: readNtuple.h:63
Float_t lide
Definition: readNtuple.h:232
Float_t KtrackLength
Definition: readNtuple.h:113
Float_t Ebalance2D
Definition: readNtuple.h:209
Int_t slcellNcell
Definition: readNtuple.h:145
Float_t trueE
Definition: readNtuple.h:55
Int_t allprongcellNcell
Definition: readNtuple.h:179
Float_t gam1Py
Definition: readNtuple.h:89
Float_t gam1Px
Definition: readNtuple.h:88
Float_t Dphi
Definition: readNtuple.h:159
readNtuple(const char *TypeNameIn, const char *TypeNameNu, const char *TypeNameInt, TChain *tree=0)
Float_t lem
Definition: readNtuple.h:233
Int_t Nplanes
Definition: readNtuple.h:108
Float_t recoE
Definition: readNtuple.h:109
Float_t rvp
Definition: readNtuple.h:230
Float_t KtrackEnergy
Definition: readNtuple.h:115
Float_t lid
Definition: readNtuple.h:231
Int_t ccnc
Definition: readNtuple.h:57
Float_t HadW
Definition: readNtuple.h:66
Float_t KtrackTheta
Definition: readNtuple.h:119
Float_t gam2Py
Definition: readNtuple.h:93
Int_t HasP2GG
Definition: readNtuple.h:82
PID
Definition: FillPIDs.h:14
Float_t trueTheta
Definition: readNtuple.h:56
Int_t prong1cellNcell
Definition: readNtuple.h:192
Int_t sliceIndex
Definition: readNtuple.h:30
Float_t KtrackEnergy_ratio
Definition: readNtuple.h:117
Int_t mesonPDG
Definition: readNtuple.h:75
Int_t Ndaus
Definition: readNtuple.h:84
Int_t sl_Nslices
Definition: readNtuple.h:28
Float_t KtrackStartZ
Definition: readNtuple.h:123
Int_t MinPlane
Definition: readNtuple.h:110
Int_t sl_goodSpill
Definition: readNtuple.h:27
Int_t PDG
Definition: readNtuple.h:58
Int_t Nshowers
Definition: readNtuple.h:234
Float_t KtrackStartX
Definition: readNtuple.h:121
Int_t sl_deltaspilltimensec
Definition: readNtuple.h:43
Float_t gam2Pz
Definition: readNtuple.h:94
Float_t lepPz
Definition: readNtuple.h:73
Float_t gam2Px
Definition: readNtuple.h:92
Float_t qsqr
Definition: readNtuple.h:67
Float_t gam1E
Definition: readNtuple.h:87
Int_t sl_subrunnum
Definition: readNtuple.h:25
TChain * fChain
Int_t sl_evtnum
Definition: readNtuple.h:26
Int_t intType
Definition: readNtuple.h:62
Float_t KtrackStopY
Definition: readNtuple.h:125
Int_t Npi0
Definition: readNtuple.h:97
Int_t sl_Ndeaddcm
Definition: readNtuple.h:35
Double_t sl_horncurrent
Definition: readNtuple.h:38
Int_t Nprongs2D
Definition: readNtuple.h:205
void Init(void)
Definition: gXSecComp.cxx:138
Float_t KtrackPhi
Definition: readNtuple.h:120
Float_t lepPy
Definition: readNtuple.h:72
OStream cout
Definition: OStream.cxx:6
virtual ~readNtuple()
Float_t HadY
Definition: readNtuple.h:65
Float_t KtrackStopZ
Definition: readNtuple.h:126
Int_t shwcellNcell
Definition: readNtuple.h:251
Int_t allprong2DcellNcell
Definition: readNtuple.h:219
Float_t OscWt
Definition: readNtuple.h:68
Double_t sl_widthX
Definition: readNtuple.h:41
Int_t lepPDG
Definition: readNtuple.h:70
ifstream in
Definition: comparison.C:7
Int_t NKtracks
Definition: readNtuple.h:127
Float_t Nmip
Definition: readNtuple.h:141
Float_t lepE
Definition: readNtuple.h:74
Float_t mesonPx
Definition: readNtuple.h:76
Int_t Nprongs
Definition: readNtuple.h:157
Double_t sl_posX
Definition: readNtuple.h:39
Float_t gam1Pz
Definition: readNtuple.h:90
Int_t sl_runnum
current Tree number in a TChain
Definition: readNtuple.h:24
Int_t mode
Definition: readNtuple.h:60
bool isQE
Definition: readNtuple.h:61
Int_t Ncells
Definition: readNtuple.h:107
Float_t Ebalance
Definition: readNtuple.h:158
Int_t origPDG
Definition: readNtuple.h:59
Int_t KtrackNcells
Definition: readNtuple.h:114
Float_t recoVx
Definition: readNtuple.h:49
Int_t PMPDG
Definition: readNtuple.h:83
Float_t truePt
Definition: readNtuple.h:54
Double_t sl_widthY
Definition: readNtuple.h:42
Float_t KtrackNplanes
Definition: readNtuple.h:118
Float_t weight
Definition: readNtuple.h:32
Float_t recoVz
Definition: readNtuple.h:51
Float_t KtrackStartY
Definition: readNtuple.h:122
Int_t HasPi0
Definition: readNtuple.h:81
Float_t mesonPz
Definition: readNtuple.h:78
Definition: FillPIDs.h:15
Float_t mesonPy
Definition: readNtuple.h:77
Float_t gam2E
Definition: readNtuple.h:91