Select clusters or cells in a given time window in order to reject noisy channels...
[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaCalorimeterQA.h
1 #ifndef ALIANACALORIMETERQA_H
2 #define ALIANACALORIMETERQA_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice     */
5 /* $Id: $ */
6
7 //_________________________________________________________________________
8 // Class to check results from simulations or reconstructed real data. 
9 // Fill few histograms and do some checking plots
10 //
11 //-- Author: Gustavo Conesa (INFN-LNF)a
12
13 // --- Root system ---
14 class TH3F;
15 class TH2F;
16 class TH1F;
17
18 // --- Analysis system --- 
19 class AliESDCaloCluster;
20 class AliAODCaloCluster;
21
22 #include "AliAnaPartCorrBaseClass.h"
23  
24 class AliAnaCalorimeterQA : public AliAnaPartCorrBaseClass {
25   
26  public: 
27   
28   AliAnaCalorimeterQA() ; // default ctor
29   AliAnaCalorimeterQA(const AliAnaCalorimeterQA & g) ; // cpy ctor
30         
31 private:
32   AliAnaCalorimeterQA & operator = (const AliAnaCalorimeterQA & g) ;//cpy assignment
33         
34 public:
35   virtual ~AliAnaCalorimeterQA() {;} //virtual dtor
36   
37   void ClusterHistograms(const TLorentzVector mom, const Double_t tof, Float_t *pos, Float_t * showerShape, 
38                                                  const Int_t nCaloCellsPerCluster, const Int_t nModule,
39                                                  const Int_t nTracksMatched, const TObject* track, 
40                                                  const Int_t * labels, const Int_t nLabels);
41         
42   TList * GetCreateOutputObjects();
43
44   void Init();
45   void InitParameters();
46   
47   void Print(const Option_t * opt) const;
48     
49   void MakeAnalysisFillHistograms() ; 
50   
51   void MCHistograms(const TLorentzVector mom, const Int_t pdg);
52
53   TString GetCalorimeter() const {return fCalorimeter ;}
54   void SetCalorimeter( TString calo ) {fCalorimeter = calo; }
55   TString GetStyleMacro() const {return fStyleMacro ;}
56   void SetStyleMacro( TString macro ) {fStyleMacro = macro; }
57   
58   void SwitchOnPlotsMaking()  {fMakePlots = kTRUE;}
59   void SwitchOffPlotsMaking() {fMakePlots = kFALSE;}
60         
61   void SwitchOnCalorimetersCorrelation()  {fCorrelateCalos = kTRUE;}
62   void SwitchOffCalorimetersCorrelation() {fCorrelateCalos = kFALSE;}
63   void CorrelateCalorimeters(TRefArray* caloClusters);
64         
65   void Terminate(TList * outputList);
66   void ReadHistograms(TList * outputList); //Fill histograms with histograms in ouput list, needed in Terminate.
67
68   Int_t GetModuleNumber(AliESDCaloCluster * cluster);
69   Int_t GetModuleNumber(AliAODCaloCluster * cluster);
70   Int_t GetModuleNumberCellIndexes(const Int_t absId, Int_t & icol, Int_t & irow, Int_t &iRCU);
71         
72   void SetNumberOfModules(Int_t nmod) {fNModules = nmod;}
73
74   void SetTimeCut(Double_t min, Double_t max) {fTimeCutMin = min; fTimeCutMax = max;}
75   Double_t GetTimeCutMin() const {return fTimeCutMin;}
76   Double_t GetTimeCutMax() const {return fTimeCutMax;}
77
78   //Histogram binning setters
79
80   Int_t GetNewRebinForRePlotting(TH1D*histo, const Float_t newXmin, const Float_t newXmax,
81                                  const Int_t newNbins) const;
82
83   virtual void SetHistoPOverERangeAndNBins(Float_t min, Float_t max, Int_t n) {
84         fHistoPOverEBins  = n ;
85         fHistoPOverEMax   = max ;
86         fHistoPOverEMin   = min ;
87   }
88         
89   Int_t   GetHistoPOverEBins()  const { return fHistoPOverEBins ; }
90   Float_t GetHistoPOverEMin()   const { return fHistoPOverEMin ; }
91   Float_t GetHistoPOverEMax()   const { return fHistoPOverEMax ; }      
92         
93         virtual void SetHistodEdxRangeAndNBins(Float_t min, Float_t max, Int_t n) {
94                 fHistodEdxBins  = n ;
95                 fHistodEdxMax   = max ;
96                 fHistodEdxMin   = min ;
97         }
98         
99         Int_t   GetHistodEdxBins()  const { return fHistodEdxBins ; }
100         Float_t GetHistodEdxMin()   const { return fHistodEdxMin ; }
101         Float_t GetHistodEdxMax()   const { return fHistodEdxMax ; }    
102
103         virtual void SetHistodRRangeAndNBins(Float_t min, Float_t max, Int_t n) {
104                 fHistodRBins  = n ;
105                 fHistodRMax   = max ;
106                 fHistodRMin   = min ;
107         }
108         
109         Int_t   GetHistodRBins()  const { return fHistodRBins ; }
110         Float_t GetHistodRMin()   const { return fHistodRMin ; }
111         Float_t GetHistodRMax()   const { return fHistodRMax ; }        
112
113         virtual void SetHistoTimeRangeAndNBins(Float_t min, Float_t max, Int_t n) {
114                 fHistoTimeBins  = n ;
115                 fHistoTimeMax   = max ;
116                 fHistoTimeMin   = min ;
117         }       
118         
119         Int_t   GetHistoTimeBins()  const { return fHistoTimeBins ; }
120         Float_t GetHistoTimeMin()   const { return fHistoTimeMin ; }
121         Float_t GetHistoTimeMax()   const { return fHistoTimeMax ; }    
122
123         virtual void SetHistoNClusterCellRangeAndNBins(Int_t min, Int_t max, Int_t n) {
124                 fHistoNBins  = n ;
125                 fHistoNMax   = max ;
126                 fHistoNMin   = min ;
127         }
128         
129         Int_t GetHistoNClusterCellBins()  const { return fHistoNBins ; }
130         Int_t GetHistoNClusterCellMin()   const { return fHistoNMin ; }
131         Int_t GetHistoNClusterCellMax()   const { return fHistoNMax ; } 
132
133         virtual void SetHistoRatioRangeAndNBins(Float_t min, Float_t max, Int_t n) {
134                 fHistoRatioBins  = n ;
135                 fHistoRatioMax   = max ;
136                 fHistoRatioMin   = min ;
137         }
138         
139         Int_t   GetHistoRatioBins()  const { return fHistoRatioBins ; }
140         Float_t GetHistoRatioMin()   const { return fHistoRatioMin ; }
141         Float_t GetHistoRatioMax()   const { return fHistoRatioMax ; }  
142
143         virtual void SetHistoVertexDistRangeAndNBins(Float_t min, Float_t max, Int_t n) {
144                 fHistoVertexDistBins = n ;
145                 fHistoVertexDistMax   = max ;
146                 fHistoVertexDistMin   = min ;
147         }
148         
149         Int_t   GetHistoVertexDistBins()  const { return fHistoVertexDistBins ; }
150         Float_t GetHistoVertexDistMin()   const { return fHistoVertexDistMin ; }
151         Float_t GetHistoVertexDistMax()   const { return fHistoVertexDistMax ; }        
152         
153         
154         virtual void SetHistoXRangeAndNBins(Float_t min, Float_t max, Int_t n) {
155                 fHistoXBins  = n ;
156                 fHistoXMax   = max ;
157                 fHistoXMin   = min ;
158         }
159         
160         Int_t   GetHistoXBins()  const { return fHistoXBins ; }
161         Float_t GetHistoXMin()   const { return fHistoXMin ; }
162         Float_t GetHistoXMax()   const { return fHistoXMax ; }  
163
164         
165         virtual void SetHistoYRangeAndNBins(Float_t min, Float_t max, Int_t n) {
166                 fHistoYBins  = n ;
167                 fHistoYMax   = max ;
168                 fHistoYMin   = min ;
169         }
170         
171         Int_t   GetHistoYBins()  const { return fHistoYBins ; }
172         Float_t GetHistoYMin()   const { return fHistoYMin ; }
173         Float_t GetHistoYMax()   const { return fHistoYMax ; }  
174
175         
176         virtual void SetHistoZRangeAndNBins(Float_t min, Float_t max, Int_t n) {
177                 fHistoZBins  = n ;
178                 fHistoZMax   = max ;
179                 fHistoZMin   = min ;
180         }
181         
182         Int_t   GetHistoZBins()  const { return fHistoZBins ; }
183         Float_t GetHistoZMin()   const { return fHistoZMin ; }
184         Float_t GetHistoZMax()   const { return fHistoZMax ; }  
185         
186         virtual void SetHistoRRangeAndNBins(Float_t min, Float_t max, Int_t n) {
187                 fHistoRBins  = n ;
188                 fHistoRMax   = max ;
189                 fHistoRMin   = min ;
190         }
191         
192         Int_t   GetHistoRBins()  const { return fHistoRBins ; }
193         Float_t GetHistoRMin()   const { return fHistoRMin ; }
194         Float_t GetHistoRMax()   const { return fHistoRMax ; }  
195         
196         virtual void SetHistoShowerShapeRangeAndNBins(Float_t min, Float_t max, Int_t n) {
197                 fHistoSSBins  = n ;
198                 fHistoSSMax   = max ;
199                 fHistoSSMin   = min ;
200         }
201         
202         Int_t   GetHistoShowerShapeBins()  const { return fHistoSSBins ; }
203         Float_t GetHistoShowerShapeMin()   const { return fHistoSSMin ; }
204         Float_t GetHistoShowerShapeMax()   const { return fHistoSSMax ; }       
205         
206         
207  private:
208   
209   TString  fCalorimeter ;    // Calorimeter selection
210   TString  fStyleMacro  ;    // Location of macro for plots style
211   Bool_t   fMakePlots   ;    // Print plots
212   Bool_t   fCorrelateCalos ; // Correlate PHOS/EMCAL clusters
213   Int_t    fNModules    ;    // Number of EMCAL/PHOS modules, set as many histogras as modules 
214   Int_t    fNRCU        ;    // Number of EMCAL/PHOS RCU, set as many histogras as RCU 
215   Double_t fTimeCutMin  ;    // Remove clusters/cells with time smaller than this value, in ns
216   Double_t fTimeCutMax  ;    // Remove clusters/cells with time larger than this value, in ns
217         
218   //Histograms
219   //Histogram Bins
220   Int_t   fHistoPOverEBins;        // p/E histogram number of bins
221   Float_t fHistoPOverEMax;         // p/E maximum value
222   Float_t fHistoPOverEMin;         // p/E minimum value
223   Int_t   fHistodEdxBins;          // dEdx histogram number of bins
224   Float_t fHistodEdxMax;           // dEdx maximum value
225   Float_t fHistodEdxMin;           // dEdx minimum value
226   Int_t   fHistodRBins;            // dR histogram number of bins
227   Float_t fHistodRMax;             // dR maximum value
228   Float_t fHistodRMin;             // dR minimum value
229   Int_t   fHistoTimeBins;          // cell time histogram number of bins
230   Float_t fHistoTimeMax;           // cell time maximum value
231   Float_t fHistoTimeMin;           // cell time minimum value
232   Int_t   fHistoNBins;             // number of clusters/cells histogram number of bins
233   Int_t   fHistoNMax;              // number maximum value
234   Int_t   fHistoNMin;              // number minimum value
235   Int_t   fHistoRatioBins;         // ratio histogram number of bins
236   Float_t fHistoRatioMax;          // ratio maximum value
237   Float_t fHistoRatioMin;          // ratio minimum value
238   Int_t   fHistoVertexDistBins;    // vertex distance histogram number of bins
239   Float_t fHistoVertexDistMax;     // vertex distance maximum value
240   Float_t fHistoVertexDistMin;     // vertex distance minimum value     
241   Int_t   fHistoRBins;             // r =sqrt(x^2+y^2+z^2) (cm) position histogram number of bins
242   Float_t fHistoRMax;              // r =sqrt(x^2+y^2+z^2) (cm)  maximum value
243   Float_t fHistoRMin;              // r =sqrt(x^2+y^2+z^2) (cm)  minimum value  
244   Int_t   fHistoXBins;             // x (cm) position histogram number of bins
245   Float_t fHistoXMax;              // x (cm) position maximum value
246   Float_t fHistoXMin;              // x (cm) position minimum value
247   Int_t   fHistoYBins;             // y (cm) position histogram number of bins
248   Float_t fHistoYMax;              // y (cm) position maximum value
249   Float_t fHistoYMin;              // y (cm) position minimum value
250   Int_t   fHistoZBins;             // z (cm) position histogram number of bins
251   Float_t fHistoZMax;              // z (cm) position maximum value
252   Float_t fHistoZMin;              // z (cm) position minimum value
253   Int_t   fHistoSSBins;            // Shower Shape parameter histogram number of bins
254   Float_t fHistoSSMax;             // Shower Shape parameter position maximum value
255   Float_t fHistoSSMin;             // Shower Shape parameter position minimum value
256         
257   //CaloClusters 
258   TH1F * fhE  ; //! E distribution, Reco
259   TH1F * fhPt ; //! pT distribution, Reco
260   TH1F * fhPhi; //! phi distribution, Reco 
261   TH1F * fhEta; //! eta distribution, Reco 
262   TH3F * fhEtaPhiE  ; //! eta vs phi vs E, Reco
263   TH1F * fhECharged  ; //! E distribution, Reco, matched with track
264   TH1F * fhPtCharged ; //! pT distribution, Reco, matched with track
265   TH1F * fhPhiCharged; //! phi distribution, Reco, matched with track 
266   TH1F * fhEtaCharged; //! eta distribution, Reco, matched with track 
267   TH3F * fhEtaPhiECharged  ; //! eta vs phi vs E, Reco, matched with track 
268   TH1F * fhEChargedNoOut  ; //! E distribution, Reco, matched with track, no outer param
269   TH1F * fhPtChargedNoOut ; //! pT distribution, Reco, matched with track, no outer param
270   TH1F * fhPhiChargedNoOut; //! phi distribution, Reco, matched with track, no outer param
271   TH1F * fhEtaChargedNoOut; //! eta distribution, Reco, matched with track, no outer param
272   TH2F * fhEtaPhiChargedNoOut  ; //! eta vs phi, Reco, matched with track, no outer param
273   TH1F * fhDeltaE  ; //! MC-Reco E distribution 
274   TH1F * fhDeltaPt ; //! MC-Reco pT distribution
275   TH1F * fhDeltaPhi; //! MC-Reco phi distribution
276   TH1F * fhDeltaEta; //! MC-Reco eta distribution
277   TH1F * fhRatioE  ; //! Reco/MC E distribution 
278   TH1F * fhRatioPt ; //! Reco/MC pT distribution
279   TH1F * fhRatioPhi; //! Reco/MC phi distribution
280   TH1F * fhRatioEta; //! Reco/MC eta distribution
281   TH2F * fh2E  ; //! E distribution, Reco vs MC
282   TH2F * fh2Pt ; //! pT distribution, Reco vs MC
283   TH2F * fh2Phi; //! phi distribution, Reco vs MC
284   TH2F * fh2Eta; //! eta distribution, Reco vs MC
285   
286   TH3F * fhLambda  ;     //! Shower ellipse axis Lambda 0 vs vs Lambda 1 vs E
287   TH2F * fhDispersion  ; //! Shower dispersion vs E
288         
289   TH2F * fhIM; //! cluster pairs invariant mass
290   TH2F * fhIMCellCut; //! cluster pairs invariant mass, n cells > 1 per cluster
291   TH2F * fhAsym; //! cluster pairs invariant mass       
292   
293   TH3F * fhNCellsPerCluster;           //! N cells per cluster vs cluster energy vs eta of cluster      
294   TH3F * fhNCellsPerClusterMIP;        //! N cells per cluster vs cluster energy vs eta of cluster, finer fixed pT bin for MIP search.
295   TH3F * fhNCellsPerClusterMIPCharged; //! N cells per cluster vs cluster energy vs eta of cluster, finer fixed pT bin for MIP search, cluster matched with track.      
296         
297   TH1F * fhNClusters; //! Number of clusters
298         
299   TH2F * fhClusterTimeEnergy;   //! Cluster Time vs Energy 
300   TH1F * fhCellTimeSpreadRespectToCellMax; //! Difference of the time of cell with maximum dep energy and the rest of cells
301   TH1F * fhCellIdCellLargeTimeSpread;      //! Cells with large time respect to max (diff > 100 ns)
302         
303   TH2F * fhRNCells ; //! R=sqrt(x^2+y^2+z^2) (cm) cluster distribution vs N cells in cluster
304   TH2F * fhXNCells ; //! X (cm) cluster distribution vs N cells in cluster
305   TH2F * fhYNCells ; //! Y (cm) cluster distribution vs N cells in cluster
306   TH2F * fhZNCells ; //! Z (cm) cluster distribution vs N cells in cluster
307         
308   TH2F * fhRE ; //! R=sqrt(x^2+y^2+z^2) (cm) cluster distribution vs cluster energy
309   TH2F * fhXE ; //! X (cm) cluster distribution vs cluster energy
310   TH2F * fhYE ; //! Y (cm) cluster distribution vs cluster energy
311   TH2F * fhZE ; //! Z (cm) cluster distribution vs cluster energy
312   TH3F * fhXYZ; //! cluster X vs Y vs Z (cm)
313         
314   TH2F * fhRCellE ; //! R=sqrt(x^2+y^2+z^2) (cm) cell distribution vs cell energy
315   TH2F * fhXCellE ; //! X (cm) cell distribution vs cell energy
316   TH2F * fhYCellE ; //! Y (cm) cell distribution vs cell energy
317   TH2F * fhZCellE ; //! Z (cm) cell distribution vs cell energy
318   TH3F * fhXYZCell; //! cell X vs Y vs Z (cm)
319                 
320   TH2F * fhDeltaCellClusterRNCells ; //! R cluster - R cell distribution (cm) vs N cells in cluster
321   TH2F * fhDeltaCellClusterXNCells ; //! X cluster - X cell distribution (cm) vs N cells in cluster
322   TH2F * fhDeltaCellClusterYNCells ; //! Y cluster - Y cell distribution (cm) vs N cells in cluster
323   TH2F * fhDeltaCellClusterZNCells ; //! Z cluster - Z cell distribution (cm) vs N cells in cluster
324         
325   TH2F * fhDeltaCellClusterRE ; //! R cluster - R cell distribution (cm) vs cluster energy
326   TH2F * fhDeltaCellClusterXE ; //! X cluster - X cell distribution (cm) vs cluster energy
327   TH2F * fhDeltaCellClusterYE ; //! Y cluster - Y cell distribution (cm) vs cluster energy
328   TH2F * fhDeltaCellClusterZE ; //! Z cluster - Z cell distribution (cm) vs cluster energy
329         
330         
331   //Calo Cells
332   TH1F * fhNCells;    //! Number of towers/crystals with signal
333   TH1F * fhAmplitude; //! Amplitude measured in towers/crystals
334   TH2F * fhAmpId;     //! Amplitude measured in towers/crystals vs id of tower. 
335   TH3F * fhEtaPhiAmp; //! eta vs phi vs amplitude, cells
336
337   TH1F * fhTime;      //! Time measured in towers/crystals
338   TH2F * fhTimeId;    //! Time vs Absolute cell Id
339   TH2F * fhTimeAmp;   //! Time vs Amplitude 
340 //  TH1F * fhT0Time;    //! T0 - EMCAL Time measured in towers/crystals
341 //  TH2F * fhT0TimeId;  //! T0 - EMCAL Time vs Absolute cell Id
342 //  TH2F * fhT0TimeAmp; //! T0 - EMCAL Time vs Amplitude 
343
344
345   //Calorimeters Correlation
346   TH2F * fhCaloCorrNClusters; // EMCAL vs PHOS, number of clusters      
347   TH2F * fhCaloCorrEClusters; // EMCAL vs PHOS, total measured cluster energy
348   TH2F * fhCaloCorrNCells; // EMCAL vs PHOS, number of cells
349   TH2F * fhCaloCorrECells; // EMCAL vs PHOS,  total measured cell energy
350         
351   //Module histograms
352   TH1F ** fhEMod  ;               //! E distribution for different module, Reco
353   TH1F ** fhNClustersMod ;        //! Number of clusters for different module, Reco
354   TH2F ** fhNCellsPerClusterMod ; //! N cells per clusters different module, Reco
355   TH1F ** fhNCellsMod ;           //! Number of towers/crystals with signal different module, Reco
356   TH2F ** fhGridCellsMod ;        //! Cells ordered in column/row for different module, Reco
357   TH2F ** fhGridCellsEMod ;       //! Cells ordered in column/row for different module, weighted with energy, Reco
358   TH2F ** fhGridCellsTimeMod ;    //! Cells ordered in column/row for different module, weighted with time, Reco
359   TH1F ** fhAmplitudeMod ;        //! Amplitude measured in towers/crystals different module, Reco
360   TH1F ** fhAmplitudeModFraction; //! Amplitude measured in towers/crystals different fractions of module, Reco
361   TH2F ** fhTimeAmpPerRCU;        //! Time vs Amplitude measured in towers/crystals different RCU
362   //TH2F ** fhT0TimeAmpPerRCU;      //! T0 - EMCAL Time vs Amplitude measured in towers/crystals different RCU
363   //TH2F ** fhTimeCorrRCU;          //! Correlate time entries in the different RCU, E > 0.3
364   TH2F ** fhIMMod;                //! cluster pairs invariant mass, different module,
365   TH2F ** fhIMCellCutMod;         //! cluster pairs invariant mass, n cells > 1 per cluster, different module
366         
367   //MC  
368   TH1F *fhGenGamPt  ; // pt of primary gamma
369   TH1F *fhGenGamEta ; // eta of primart gamma
370   TH1F *fhGenGamPhi ; // phi of primary gamma   
371   TH1F *fhGenPi0Pt  ; // pt of primary pi0
372   TH1F *fhGenPi0Eta ; // eta of primart pi0
373   TH1F *fhGenPi0Phi ; // phi of primary pi0     
374   TH1F *fhGenEtaPt  ; // pt of primary eta
375   TH1F *fhGenEtaEta ; // eta of primart eta
376   TH1F *fhGenEtaPhi ; // phi of primary eta
377   TH1F *fhGenOmegaPt  ; // pt of primary omega
378   TH1F *fhGenOmegaEta ; // eta of primart omega
379   TH1F *fhGenOmegaPhi ; // phi of primary omega
380   TH1F *fhGenElePt  ; // pt of primary electron
381   TH1F *fhGenEleEta ; // eta of primart electron
382   TH1F *fhGenElePhi ; // phi of primary electron
383         
384   //TH3F * fhEMVxyz    ; // Electromagnetic particle production vertex
385   TH2F * fhEMVxyz    ; // Electromagnetic particle production vertex
386   TH2F * fhEMR       ; // Electromagnetic distance to vertex vs rec energy  
387   //TH3F * fhHaVxyz    ; // Hadron production vertex
388   TH2F * fhHaVxyz    ; // Hadron production vertex
389   TH2F * fhHaR       ; // Hadron distance to vertex vs rec energy  
390
391   TH2F * fhGamE  ; //! E distribution of generated photons, Reco
392   TH2F * fhGamPt ; //! pT distribution of generated photons, Reco
393   TH2F * fhGamPhi; //! phi distribution of generated photon, Reco 
394   TH2F * fhGamEta; //! eta distribution of generated photons, Reco 
395   TH1F * fhGamDeltaE  ; //! MC-Reco E distribution of generated photons 
396   TH1F * fhGamDeltaPt ; //! MC-Reco pT distribution of generated photons
397   TH1F * fhGamDeltaPhi; //! MC-Reco phi distribution of generated photons
398   TH1F * fhGamDeltaEta; //! MC-Reco eta distribution of generated photons
399   TH1F * fhGamRatioE  ; //! Reco/MC E distribution of generated photons 
400   TH1F * fhGamRatioPt ; //! Reco/MC pT distribution of generated photons
401   TH1F * fhGamRatioPhi; //! Reco/MC phi distribution of generated photons
402   TH1F * fhGamRatioEta; //! Reco/MC eta distribution of generated photons
403   TH2F * fhEleE  ; //! E distribution of generated electrons, Reco
404   TH2F * fhElePt ; //! pT distribution of generated electrons, Reco
405   TH2F * fhElePhi; //! phi distribution of generated electron, Reco 
406   TH2F * fhEleEta; //! eta distribution of generated electrons, Reco            
407   TH2F * fhPi0E  ; //! E distribution of generated pi0, Reco, gamma decay overlapped
408   TH2F * fhPi0Pt ; //! pT distribution of generated pi0, Reco, gamma decay overlapped
409   TH2F * fhPi0Phi; //! phi distribution of generated pi0, Reco, gamma decay overlapped
410   TH2F * fhPi0Eta; //! eta distribution of generated pi0, Reco, gamma decay overlapped
411   TH2F * fhNeHadE  ; //! E distribution of generated neutral hadron, Reco
412   TH2F * fhNeHadPt ; //! pT distribution of generated neutral hadron, Reco
413   TH2F * fhNeHadPhi; //! phi distribution of generated neutral hadron, Reco 
414   TH2F * fhNeHadEta; //! eta distribution of generated neutral hadron, Reco     
415   TH2F * fhChHadE  ; //! E distribution of generated charged hadron, Reco
416   TH2F * fhChHadPt ; //! pT distribution of generated charged hadron, Reco
417   TH2F * fhChHadPhi; //! phi distribution of generated charged hadron, Reco 
418   TH2F * fhChHadEta; //! eta distribution of generated charged hadron, Reco 
419
420   TH2F * fhGamECharged  ; //! E distribution of generated photons, Reco, track matched cluster
421   TH2F * fhGamPtCharged ; //! pT distribution of generated photons, Reco, track matched cluster
422   TH2F * fhGamPhiCharged; //! phi distribution of generated photon, Reco, track matched cluster 
423   TH2F * fhGamEtaCharged; //! eta distribution of generated photons, Reco, track matched cluster 
424   TH2F * fhEleECharged  ; //! E distribution of generated electrons, Reco, track matched cluster
425   TH2F * fhElePtCharged ; //! pT distribution of generated electrons, Reco, track matched cluster
426   TH2F * fhElePhiCharged; //! phi distribution of generated electron, Reco, track matched cluster 
427   TH2F * fhEleEtaCharged; //! eta distribution of generated electrons, Reco, track matched cluster              
428   TH2F * fhPi0ECharged  ; //! E distribution of generated pi0, Reco, gamma decay overlapped, track matched cluster
429   TH2F * fhPi0PtCharged ; //! pT distribution of generated pi0, Reco, gamma decay overlapped, track matched cluster
430   TH2F * fhPi0PhiCharged; //! phi distribution of generated pi0, Reco, gamma decay overlapped, track matched cluster
431   TH2F * fhPi0EtaCharged; //! eta distribution of generated pi0, Reco, gamma decay overlapped, track matched cluster
432   TH2F * fhNeHadECharged  ; //! E distribution of generated neutral hadron, Reco, track matched cluster
433   TH2F * fhNeHadPtCharged ; //! pT distribution of generated neutral hadron, Reco, track matched cluster
434   TH2F * fhNeHadPhiCharged; //! phi distribution of generated neutral hadron, Reco , track matched cluster
435   TH2F * fhNeHadEtaCharged; //! eta distribution of generated neutral hadron, Reco, track matched cluster       
436   TH2F * fhChHadECharged  ; //! E distribution of generated charged hadron, Reco, track matched cluster
437   TH2F * fhChHadPtCharged ; //! pT distribution of generated charged hadron, Reco, track matched cluster
438   TH2F * fhChHadPhiCharged; //! phi distribution of generated charged hadron, Reco, track matched cluster 
439   TH2F * fhChHadEtaCharged; //! eta distribution of generated charged hadron, Reco, track matched cluster       
440         
441   TH1F *fhGenGamAccE   ; // E of primary gamma
442   TH1F *fhGenGamAccPt  ; // pt of primary gamma
443   TH1F *fhGenGamAccEta ; // eta of primart gamma
444   TH1F *fhGenGamAccPhi ; // phi of primary gamma        
445   TH1F *fhGenPi0AccE   ; // E of primary pi0
446   TH1F *fhGenPi0AccPt  ; // pt of primary pi0
447   TH1F *fhGenPi0AccEta ; // eta of primart pi0
448   TH1F *fhGenPi0AccPhi ; // phi of primary pi0          
449         
450   //Histograms for track-matching
451   TH2F *fh1pOverE;     //! p/E for track-cluster matches
452   TH1F *fh1dR;         //! distance between projected track and cluster
453   TH2F *fh2EledEdx;    //! dE/dx vs. momentum for electron candidates
454   TH2F *fh2MatchdEdx;  //! dE/dx vs. momentum for all matches
455         
456   TH2F *fhMCEle1pOverE;     //! p/E for track-cluster matches, MC electrons
457   TH1F *fhMCEle1dR;         //! distance between projected track and cluster, MC electrons
458   TH2F *fhMCEle2MatchdEdx;  //! dE/dx vs. momentum for all matches, MC electrons        
459         
460   TH2F *fhMCChHad1pOverE;     //! p/E for track-cluster matches, MC charged hadrons
461   TH1F *fhMCChHad1dR;         //! distance between projected track and cluster, MC charged hadrons
462   TH2F *fhMCChHad2MatchdEdx;  //! dE/dx vs. momentum for all matches, MC charged
463         
464   TH2F *fhMCNeutral1pOverE;     //! p/E for track-cluster matches, MC neutral
465   TH1F *fhMCNeutral1dR;         //! distance between projected track and cluster, MC neutral
466   TH2F *fhMCNeutral2MatchdEdx;  //! dE/dx vs. momentum for all matches, MC neutral      
467         
468   TH2F *fh1pOverER02;           //! p/E for track-cluster matches, dR > 0.2     
469   TH2F *fhMCEle1pOverER02;      //! p/E for track-cluster matches, dR > 0.2, MC electrons
470   TH2F *fhMCChHad1pOverER02;    //! p/E for track-cluster matches, dR > 0.2, MC charged hadrons
471   TH2F *fhMCNeutral1pOverER02;  //! p/E for track-cluster matches, dR > 0.2, MC neutral
472         
473         ClassDef(AliAnaCalorimeterQA,9)
474 } ;
475
476
477 #endif //ALIANACALORIMETERQA_H
478
479
480