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