AliPhysics  e59a9ba (e59a9ba)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnaCalorimeterQA.h
Go to the documentation of this file.
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 
6 //_________________________________________________________________________
30 
31 // --- Root system ---
32 class TH3F;
33 class TH2F;
34 class TH1F;
35 class TObjString;
36 class TObjArray;
37 
38 // --- Analysis system ---
39 class AliVCaloCells;
40 class AliVCaloCluster;
41 class AliVTrack;
42 
44 
46 
47 public:
48 
50 
52  virtual ~AliAnaCalorimeterQA() { ; }
53 
54  // General methods
55 
56  TObjString * GetAnalysisCuts();
57 
58  TList * GetCreateOutputObjects();
59 
60  void Init();
61 
62  void InitParameters();
63 
65 
66  void Print(const Option_t * opt) const;
67 
68  // Main methods
69 
70  void BadClusterHistograms(AliVCluster* clus, const TObjArray *caloClusters, AliVCaloCells * cells,
71  Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Double_t tmax);
72 
73  void CalculateAverageTime(AliVCluster *clus, AliVCaloCells *cells, Double_t timeAverages[2]);
74 
75  void CellHistograms(AliVCaloCells * cells);
76 
77  void CellInClusterPositionHistograms(AliVCluster* cluster);
78 
79  void ClusterAsymmetryHistograms(AliVCluster* clus, Int_t absIdMax, Bool_t goodCluster );
80 
81  void ClusterHistograms(AliVCluster* cluster, const TObjArray *caloClusters, AliVCaloCells * cells,
82  Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Double_t tmax);
83 
84  void ClusterLoopHistograms(const TObjArray * clusters, AliVCaloCells * cells);
85 
86  Bool_t ClusterMCHistograms(Bool_t matched, const Int_t * labels, Int_t nLabels, Int_t & pdg );
87 
88  void ClusterMatchedWithTrackHistograms(AliVCluster* clus, Bool_t mcOK, Int_t pdg);
89 
90  void Correlate();
91 
92  void ExoticHistograms(Int_t absIdMax, Float_t ampMax,
93  AliVCluster *clus, AliVCaloCells* cells);
94 
95  Float_t GetECross(Int_t absId, AliVCaloCells* cells,Float_t dtcut = 10000);
96 
97  void InvariantMassHistograms(Int_t iclus, Int_t nModule, const TObjArray* caloClusters, AliVCaloCells * cells);
98 
99  Bool_t IsGoodCluster(Int_t absIdMax, Float_t m02, Int_t nCellsPerCluster, AliVCaloCells *cells);
100 
101  void MCHistograms();
102 
103  void WeightHistograms(AliVCluster *clus, AliVCaloCells* cells);
104 
105 
106  // Setters and getters
107 
108  Float_t GetEMCALCellAmpMin() const { return fEMCALCellAmpMin ; }
109  void SetEMCALCellAmpMin(Float_t amp) { fEMCALCellAmpMin = amp ; }
110 
111  Float_t GetPHOSCellAmpMin() const { return fPHOSCellAmpMin ; }
112  void SetPHOSCellAmpMin (Float_t amp) { fPHOSCellAmpMin = amp ; }
113 
116 
119 
120  Float_t GetInvMassMinECut() const { return fInvMassMinECut ; }
121  void SetInvMassMinECut(Float_t cut) { fInvMassMinECut = cut ; }
122 
123  Float_t GetInvMassMaxECut() const { return fInvMassMaxECut ; }
124  void SetInvMassMaxECut(Float_t cut) { fInvMassMaxECut = cut ; }
125 
126  Float_t GetInvMassMinM02Cut() const { return fInvMassMinM02Cut ; }
127  void SetInvMassMinM02Cut(Float_t cut) { fInvMassMinM02Cut = cut ; }
128 
129  Float_t GetInvMassMaxM02Cut() const { return fInvMassMaxM02Cut ; }
130  void SetInvMassMaxM02Cut(Float_t cut) { fInvMassMaxM02Cut = cut ; }
131 
132  Float_t GetInvMassMaxOpenAngle() const { return fInvMassMaxOpenAngle; }
134 
137 
138  Double_t GetTimeCutMin() const { return fTimeCutMin ; }
139  Double_t GetTimeCutMax() const { return fTimeCutMax ; }
140  void SetTimeCut(Double_t min, Double_t max) {
141  fTimeCutMin = min ; fTimeCutMax = max ; }
142 
143  // Histogram switchs
144 
147 
150 
153 
156 
159 
162 
165 
168 
171 
172  void SwitchOnCorrelation() { fCorrelate = kTRUE ; }
173  void SwitchOffCorrelation() { fCorrelate = kFALSE ; }
174 
177 
178  // Analysis not to be used in QA
179 
182 
183  void SwitchOnStudyWeight() { fStudyWeight = kTRUE ; }
184  void SwitchOffStudyWeight() { fStudyWeight = kFALSE ; }
185 
186  void SwitchOnStudyExotic() { fStudyExotic = kTRUE ; }
187  void SwitchOffStudyExotic() { fStudyExotic = kFALSE ; }
188 
189  void SetNECrossCuts(Int_t n) { fExoNECrossCuts = n ; }
190  void SetNDTimeCuts (Int_t n) { fExoNDTimeCuts = n ; }
191 
192  void SetExoECrossCuts (Int_t i, Float_t c) { if (i<fExoNECrossCuts) fExoECrossCuts[i] = c ; }
193  void SetExoDTimeCuts (Int_t i, Float_t c) { if (i<fExoNDTimeCuts ) fExoDTimeCuts [i] = c ; }
194 
195  void SetConstantTimeShift(Float_t shift) { fConstantTimeShift = shift ; }
196 
197 
198  private:
199 
200  // Switches
201 
205  Bool_t fFillAllTH3 ;
206  Bool_t fFillAllTMHisto ;
207 
212 
213  Bool_t fCorrelate ;
215 
216  // Analysis not to be used in QA
217 
219  Bool_t fStudyExotic;
220  Bool_t fStudyWeight;
221 
222  // Parameters
223 
224  Int_t fNModules ;
225  Int_t fNRCU ;
226  Int_t fNMaxCols ;
227  Int_t fNMaxRows ;
228 
229  // Cuts
230 
231  Double_t fTimeCutMin ;
232  Double_t fTimeCutMax ;
233  Float_t fCellAmpMin;
235  Float_t fPHOSCellAmpMin ;
238 
239  // Invariant mass analysis
240 
241  Float_t fInvMassMinECut;
242  Float_t fInvMassMaxECut;
247  // Exotic studies
248 
249  Float_t fExoNECrossCuts ;
250  Float_t fExoECrossCuts[10];
251  Float_t fExoNDTimeCuts ;
252  Float_t fExoDTimeCuts[5] ;
253 
254  TLorentzVector fClusterMomentum;
255  TLorentzVector fClusterMomentum2;
256  TLorentzVector fPrimaryMomentum;
257 
259 
260 
261  // Calorimeter Clusters
262 
263  TH1F * fhE ;
264  TH1F * fhPt ;
265  TH1F * fhPhi;
266  TH1F * fhEta;
267  TH2F * fhEtaPhi;
268  TH3F * fhEtaPhiE ;
269  TH1F * fhECharged ;
270  TH1F * fhPtCharged ;
271  TH1F * fhPhiCharged;
272  TH1F * fhEtaCharged;
275 
276  TH2F * fhIM;
277  TH2F * fhIMSame;
278  TH2F * fhIMDiff;
279  TH2F * fhIMDCAL;
280  TH2F * fhIMDCALSame;
281  TH2F * fhIMDCALDiff;
282  TH2F * fhIMDCALPHOS;
284  TH2F * fhIMEMCALPHOS;
286 
287  TH2F * fhAsym;
288 
289  TH2F* fhOpAngle;
291 
294 
297 
300 
301  TH1F * fhNClusters;
302 
308 
313 
314 // TH2F * fhClusterMaxCellDiffAverageTime; //!<! Difference between cluster average time and time of cell with more energy
315 // TH2F * fhClusterMaxCellDiffWeightedTime; //!<! Difference between cluster weighted time and time of cell with more energy
317 
318  TH2F * fhLambda0;
319  TH2F * fhLambda1;
320 // TH2F * fhDispersion; //!<! Cluster Dispersion vs Energy
321 
322  // Bad clusters histograms
323 
329 
333 
337 
340 
341  // Cluster cell size
342 
346  TH2F * fhDeltaIA[2];
347  TH2F * fhDeltaIAL0[2];
348  TH2F * fhDeltaIAL1[2];
349  TH2F * fhDeltaIANCells[2] ;
350  TH2F * fhDeltaIAMC[4];
355 
356  // Cluster/cell Position
357 
358  TH2F * fhRNCells ;
359  TH2F * fhXNCells ;
360  TH2F * fhYNCells ;
361  TH2F * fhZNCells ;
362 
363  TH2F * fhRE ;
364  TH2F * fhXE ;
365  TH2F * fhYE ;
366  TH2F * fhZE ;
367  TH3F * fhXYZ;
368 
369  TH2F * fhRCellE ;
370  TH2F * fhXCellE ;
371  TH2F * fhYCellE ;
372  TH2F * fhZCellE ;
373  TH3F * fhXYZCell;
374 
379 
384 
385  // Calorimeter cells
386 
387  TH1F * fhNCells;
389  TH1F * fhAmplitude;
390  TH2F * fhAmpId;
392  TH2F * fhEtaPhiCell;
393 
394  TH1F * fhTime;
395  //TH2F * fhTimeVz; //!<! Time measured in towers/crystals vs vertex z component, for E > 0.5
396  TH2F * fhTimeId;
397  TH2F * fhTimeAmp;
398 
399  TH2F * fhAmpIdLowGain;
402 
403  TH2F * fhCellECross;
404 
405  // Calorimeters Correlation
406 
411 
416 
421 
422  // V0 Correlation
423 
432 
433  // Track Correlation
434 
439 
440  // Centrality
441 
446 
447  // Event plane
448 
453 
454  // Module histograms
455 
456  TH2F * fhEMod ;
457  TH2F * fhAmpMod ;
458  TH2F * fhEWeirdMod ;
459  TH2F * fhAmpWeirdMod ;
460  TH2F * fhTimeMod ;
461  TH2F * fhNClustersMod ;
462  TH2F * fhNCellsMod ;
471  TH2F * fhGridCells ;
472  TH2F * fhGridCellsE ;
473  TH2F * fhGridCellsTime ;
477  TH2F ** fhTimeAmpPerRCU;
478  TH2F ** fhIMMod;
479 
480  // Weight studies
481 
486 
487  TH2F * fhLambda0ForW0AndCellCuts [12][4][3];
488  TH2F * fhLambda0ForW0AndCellCutsEta0[12][4][3];
489 //TH2F * fhLambda1ForW0AndCellCuts [12][4][3]; //!<! L1 for different w0 and cell cuts
490 
491  TH2F * fhLambda0ForW0MC[12][5];
492 //TH2F * fhLambda1ForW0MC[12][5]; //!<! L1 for different w0, depending on the particle of origin
493 
498 
499  // Exotic studies
500 
501  TH2F * fhExoNCell [10][5] ;
502  TH2F * fhExoL0 [10][5] ;
503  TH2F * fhExoL1 [10][5] ;
504  TH2F * fhExoECross [10][5] ;
505  TH2F * fhExoTime [10][5] ;
506  TH2F * fhExoDTime [10] ;
507  TH2F * fhExoL0NCell[10][5] ;
508  TH2F * fhExoL0ECross ;
509  TH2F * fhExoL1NCell[10][5] ;
510  TH2F * fhExoL1ECross ;
511 
512  // Pure MC histograms
513 
515  enum mcTypes { kmcPhoton = 0, kmcPi0 = 1, kmcEta = 2,
518 
519  TH2F * fhRecoMCE[7][2] ;
520  TH2F * fhRecoMCPhi[7][2] ;
521  TH2F * fhRecoMCEta[7][2] ;
522  TH2F * fhRecoMCDeltaE[7][2] ;
523  TH2F * fhRecoMCRatioE[7][2] ;
524  TH2F * fhRecoMCDeltaPhi[7][2];
525  TH2F * fhRecoMCDeltaEta[7][2];
526 
527  TH1F * fhGenMCE [4] ;
528  TH1F * fhGenMCPt[4] ;
529  TH2F * fhGenMCEtaPhi[4] ;
530  TH1F * fhGenMCAccE [4] ;
531  TH1F * fhGenMCAccPt[4] ;
532  TH2F * fhGenMCAccEtaPhi[4] ;
533 
534  TH2F * fhEMVxyz ;
535  TH2F * fhEMR ;
536  TH2F * fhHaVxyz ;
537  TH2F * fhHaR ;
538 
539  // Histograms for MC track-matching
540 
541  TH2F * fh1EOverP;
542  TH2F * fh2dR;
543  TH2F * fh2EledEdx;
544  TH2F * fh2MatchdEdx;
545  TH2F * fh1EOverPR02;
546  TH2F * fh1EleEOverP;
547 
548  TH2F * fh1EOverPMod;
549  TH2F * fh2dRMod;
550  TH2F * fh2EledEdxMod;
554 
555  TH2F * fhMCEle1EOverP;
556  TH1F * fhMCEle1dR;
558 
560  TH1F * fhMCChHad1dR;
562 
564  TH1F * fhMCNeutral1dR;
566 
570 
574 
578 
582 
587 
590 
593 
595  ClassDef(AliAnaCalorimeterQA,33) ;
597 
598 } ;
599 
600 #endif //ALIANACALORIMETERQA_H
601 
602 
603 
void SetEMCALNCellsPerClusterMin(Int_t n)
TH2F * fhCaloV0MCorrEClusters
! Calo vs V0 multiplicity, total measured cluster energy
TH2F * fhTimeIdLowGain
! Time vs Absolute cell Id, low gain
TH2F * fhDeltaCellClusterRE
! R cluster - R cell distribution (cm) vs cluster energy
Bool_t IsGoodCluster(Int_t absIdMax, Float_t m02, Int_t nCellsPerCluster, AliVCaloCells *cells)
Identify cluster as exotic or not.
TH2F * fhDCALPHOSCorrNCells
! DCAL vs PHOS, number of cells
Int_t pdg
TH2F * fhNCellsPerClusterNoCut
! N cells per cluster vs cluster energy vs eta of cluster
TH2F * fhLambda0ForW0AndCellCutsEta0[12][4][3]
! L0 for different w0 and cell cuts, |eta| < 0.15
TH2F * fhMCNeutral1EOverPR02
! p/E for track-cluster matches, dR < 0.2, MC neutral
TH2F * fhNCellsMod
! Number of towers/crystals with signal for different module, Reco
TH3F * fhXYZ
! cluster X vs Y vs Z (cm)
TH2F * fh1EOverPMod
! p/E for track-cluster matches, per SM
TH2F * fh1EOverPR02Mod
! p/E for track-cluster matches, dR < 0.2, per SM
TH2F * fhDeltaCellClusterYNCells
! Y cluster - Y cell distribution (cm) vs N cells in cluster
TH1F * fhPt
! pT distribution, Reco
void SetPHOSNCellsPerClusterMin(Int_t n)
void SetTimeCut(Double_t min, Double_t max)
TH2F * fhClusterMaxCellECross
! 1 - Energy in cross around max energy cell / max energy cell vs cluster energy, good clusters ...
Bool_t fFillPi0PairDiffTime
Fill time difference histograms of cluster pairs in pi0 mass window, only if fFillAllPi0Histo=kTRUE.
Float_t GetInvMassMaxOpenAngle() const
TH2F * fhGridCellsE
! Cells ordered in column/row for different module, weighted with energy, Reco
TH1F * fhPtCharged
! pT distribution, Reco, matched with track
TH2F * fhClusterMaxCellDiff
! Difference between cluster energy and energy of cell with more energy, good clusters only ...
TH2F * fhBadClusterDeltaIEtaDeltaIPhiE2
! Difference between max cell index and farthest cell, eta vs phi, 2 < E < 6 GeV, with and without ma...
TH2F * fh2dR
! distance between projected track and cluster (eta-phi units)
Float_t fInvMassMinECut
Minimum energy cut value for clusters entering the invariant mass calculation.
Bool_t fStudyClustersAsymmetry
Study asymmetry of clusters, not QA related.
TH2F * fh2dRMod
! distance between projected track and cluster (eta-phi units), per SM
TH2F * fhExoL1[10][5]
! Short shower shape axis for exotic
void SwitchOnFillAllTrackMatchingHistogram()
TH2F * fhExoNCell[10][5]
! Number of cells per cluster for different cuts
TH2F * fhBadClusterEtaPhi
! Time Max cell of bad cluster
Float_t fEMCALCellAmpMin
Amplitude Threshold on EMCal cells.
TH2F * fhIMDCALPHOS
! Cluster pairs invariant mass vs pair pT, for DCal-PHOS pairs
TH1F * fhEta
! eta distribution, Reco
TH2F * fhTimeMod
! Cell time distribution for different module, Reco
void SetNECrossCuts(Int_t n)
TH2F * fhMCChHad1EleEOverP
! p/E for track-cluster matches, dR < 0.2, 60 < dEdx < 100, MC charged hadrons
Bool_t fStudyWeight
Study the energy weight used in different cluster calculations, not QA related.
virtual ~AliAnaCalorimeterQA()
Virtual destructor. Not implemented.
TH2F * fh2MatchdEdxMod
! dE/dx for all matches, per SM
TH2F * fh2EledEdx
! dE/dx vs. momentum for electron candidates
TH2F * fhTrackMatchedDPhiNeg
! Phi distance between track and cluster vs cluster E, after and before photon cuts ...
TH2F * fhSumCellsAmpMod
! Sum of towers/crystals signal for different module, Reco
Bool_t fFillInvMassOpenAngle
Fill opening angle histograms of cluster pairs, only if fFillAllPi0Histo=kTRUE.
TH2F * fhExoL0NCell[10][5]
! Lambda0 vs n cells in cluster for several E cross cuts and cluster with E > 5
TH2F * fhCaloTrackMCorrNClusters
! Calo vs Track Multiplicity, number of clusters
TH3F * fhEtaPhiE
! eta vs phi vs E, Reco
Float_t fPHOSCellAmpMin
Amplitude Threshold on PHOS cells.
TH1F * fhBadClusterEnergy
! Energy of bad cluster
TH2F * fhBadClusterMaxCellDiffWeightedTime
! Difference between cluster weighted time and time of cell with more energy
Float_t GetInvMassMaxECut() const
TH2F * fhRecoMCRatioE[7][2]
! Reco/Gen E generated particle vs reconstructed E
TH2F * fhBadClusterMaxCellDiffAverageTime
! Difference between cluster average time and time of cell with more energy
TH2F * fhCaloCenECells
! Calo vs centrality, total measured cell energy
TH2F * fhIMDCALDiff
! Cluster pairs invariant mass vs pair pT, for DCal pairs
TH2F ** fhNClustersSumEnergyPerMod
! N clusters vs sum of energies in different module, Reco
TH2F * fhEMCALDCALCorrEClusters
! EMCAL vs DCAL, total measured cluster energy
TH1F * fhGenMCAccE[4]
! pt of primary particle, in acceptance
TH2F * fhXNCells
! X (cm) cluster distribution vs N cells in cluster
TH2F * fhRCellE
! R=sqrt(x^2+y^2) (cm) cell distribution vs cell energy
TH1F * fhGenMCPt[4]
! pt of primary particle
TH2F * fhClusterPairDiffTimeESameMod
! Pair of clusters time difference vs E, in same Mod
void ClusterMatchedWithTrackHistograms(AliVCluster *clus, Bool_t mcOK, Int_t pdg)
TH2F * fhEtaPhi
! eta-phi distribution, Reco
TH2F * fhIMSame
! Cluster pairs invariant mass vs pair pT, for EMCAL or PHOS pairs
TH1F * fhGenMCAccPt[4]
! pt of primary particle, in acceptance
Double_t fTimeCutMax
Remove clusters/cells with time larger than this value, in ns.
Float_t fCellAmpMin
Amplitude Threshold on calorimeter cells, set at execution time.
TH2F * fhIMEMCALPHOSSame
! Cluster pairs invariant mass vs pair pT, for EMCAL(DCal eta acceptance)-EMCAL (PHOS eta acceptance)...
Int_t fNModules
Number of EMCAL/PHOS modules.
TH2F * fhCaloEvPECells
! Calo vs event plane angle, total measured cell energy
void SetInvMassMinM02Cut(Float_t cut)
TH1F * fhNCellsCutAmpMin
! Number of towers/crystals with signal, with min amplitude
TH2F ** fhNCellsPerClusterModNoCut
! N cells per clusters different module, Reco, No cut
TH2F * fhClusterMaxCellCloseCellRatio
! Ratio between max cell energy and cell energy of the same cluster
void InvariantMassHistograms(Int_t iclus, Int_t nModule, const TObjArray *caloClusters, AliVCaloCells *cells)
TH2F * fhRecoMCDeltaE[7][2]
! Gen-Reco E generated particle vs reconstructed E
TH2F ** fhTimeAmpPerRCU
! Time vs Amplitude measured in towers/crystals different RCU
Class for the Calorimeter QA analysis.
TH2F * fhExoDTime[10]
! Difference in time between cell with max energy and rest of cells for exotic
void SwitchOffFillInvMassOpAngleHistogram()
TH2F * fhLambda0ForW0AndCellCuts[12][4][3]
! L0 for different w0 and cell cuts
void SetInvMassMaxM02Cut(Float_t cut)
Float_t fInvMassMaxECut
Maximum energy cut value for clusters entering the invariant mass calculation.
void WeightHistograms(AliVCluster *clus, AliVCaloCells *cells)
Bool_t fFillAllPosHisto
Fill all the position related histograms.
TH2F ** fhECellTotalRatioMod
! e cell / e total vs e total, per SM
TH2F * fhEMCALPHOSCorrEClusters
! EMCAL vs PHOS, total measured cluster energy
mcTypes
Enumerator with indeces for MC histograms array indicating the particle type generating the cluster...
TH2F * fhNCellsPerClusterWeirdModNoCut
! N cells per clusters different module, Reco, No cut, ridiculously large energy
Bool_t fFillAllTH3
Fill TH3 histograms.
TH2F * fhClusterTimeEnergy
! Cluster Time vs Energy
TH2F * fhDeltaCellClusterRNCells
! R cluster - R cell distribution (cm) vs N cells in cluster
TH1F * fhMCChHad1dR
! distance between projected track and cluster, MC charged hadrons
TH1F * fhTime
! Time measured in towers/crystals
TH2F * fhCaloCenNClusters
! Calo vs centrality, number of clusters
TH2F * fhBadClusterMaxCellDiff
! Difference between cluster energy and energy of cell with more energy
TH2F * fhDeltaIEtaDeltaIPhiE6[2]
! Difference between max cell index and farthest cell, eta vs phi, E > 6 GeV, with and without matchi...
TH2F * fhBadClusterMaxCellECross
! 1 - Energy in cross around max energy cell / max energy cell vs cluster energy, bad clusters ...
TH2F * fhTrackMatchedDEtaDPhiNeg
! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before ...
TH1F * fhPhiCharged
! phi distribution, Reco, matched with track
TH2F * fh1EOverPR02
! p/E for track-cluster matches, dR < 0.2
TH2F * fhMCEle1EleEOverP
! p/E for track-cluster matches, dR < 0.2, 60 < dEdx < 100, MC electrons
TH2F * fhBadClusterMaxCellCloseCellDiff
! Difference between max cell energy and cell energy of the same cluster for bad clusters ...
void MakeAnalysisFillHistograms()
Main task method, call all the methods filling QA histograms.
TH2F * fhCellECross
! 1 - Energy in cross around cell / cell energy
TH2F * fhClusterPairDiffTimeE
! Pair of clusters time difference vs E
TH2F * fhBadClusterDeltaIEtaDeltaIPhiE0
! Difference between max cell index and farthest cell, eta vs phi, E < 2 GeV, with and without matchi...
TH2F * fhClusterPairDiffTimeEPi0Mass
! EMCal/PHOS Cluster time TOF difference, for pairs in 0.1 < mass < 0.18
Double_t GetTimeCutMin() const
Int_t fNRCU
Number of EMCAL/PHOS RCU.
TH2F * fhBadClusterMaxCellCloseCellRatio
! Ratio between max cell energy and cell energy of the same cluster for bad clusters ...
TH1F * fhCellIdCellLargeTimeSpread
! Cells with large time respect to max (diff > 100 ns)
void SwitchOffFillInvMassReducedEMCALHistogram()
TH2F * fhEMCALPHOSCorrNClusters
! EMCAL vs PHOS, number of clusters
void SwitchOnFillPi0PairDiffTimeHistogram()
TH2F * fhIMDCALSame
! Cluster pairs invariant mass vs pair pT, for DCal pairs
TH2F * fhCaloTrackMCorrEClusters
! Calo vs Track Multiplicity, total measured cluster energy
TH2F * fhZNCells
! Z (cm) cluster distribution vs N cells in cluster
TH2F * fhExoL0[10][5]
! Long shower shape axis for exotic
Float_t GetEMCALCellAmpMin() const
TH1F * fhE
! E distribution, Reco
TH2F * fhTrackMatchedDPhiPosMod
! Phi distance between positive track and cluster vs module for E > 0.5 GeV
TH2F * fhDeltaCellClusterYE
! Y cluster - Y cell distribution (cm) vs cluster energy
TH2F * fhZE
! Z (cm) cluster distribution vs cluster energy
TH2F * fhEMCALDCALCorrNClusters
! EMCAL vs DCAL, number of clusters
TH1F * fhEtaCharged
! eta-phi distribution, Reco, matched with track
Bool_t fFillAllPi0Histo
Fill invariant mass histograms.
AliAnaCalorimeterQA()
Default Constructor. Initialize parameters.
TH1F * fhPhi
! phi distribution, Reco
Base class for CaloTrackCorr analysis algorithms.
Int_t GetEMCALNCellsPerClusterMin() const
Bool_t fFillInvMassInEMCALWithPHOSDCalAcc
Fill invariant mass histograms of EMCal clusters in DCal and PHOS eta acceptance each, only if fFillAllPi0Histo=kTRUE.
TH2F * fhBadClusterPairDiffTimeE
! Pair of clusters time difference vs E, bad cluster
TLorentzVector fClusterMomentum
! Cluster momentum, temporary container
TH2F * fhEtaPhiCharged
! eta distribution, Reco, matched with track
TH2F * fhRecoMCDeltaPhi[7][2]
! Gen-Reco phi generated particle vs reconstructed E
TLorentzVector fClusterMomentum2
! Cluster momentum, temporary container
TH2F * fhECellTotalLogRatio
! log (e cell / e total) vs e total
TH2F * fhClusterMaxCellDiffNoCut
! Difference between cluster energy and energy of cell with more energy, no bad cluster rejection ...
Float_t fExoECrossCuts[10]
List of ecross cuts.
TH2F * fhMCNeutral2MatchdEdx
! dE/dx vs. momentum for all matches, MC neutral
TH2F * fhCaloCenEClusters
! Calo vs centrality, total measured cluster energy
void ExoticHistograms(Int_t absIdMax, Float_t ampMax, AliVCluster *clus, AliVCaloCells *cells)
Fill histograms with exoticity parameters.
TH3F * fhEtaPhiECharged
! eta vs phi vs E, Reco, matched with track
TH2F * fhEMaxCellClusterLogRatio
! log (e max cell / e cluster) vs e cluster
TH2F * fhBadCellTimeSpreadRespectToCellMax
! Difference of the time of cell with maximum dep energy and the rest of cells for bad clusters ...
TH2F * fhDeltaCellClusterXE
! X cluster - X cell distribution (cm) vs cluster energy
TH2F * fhEMCALPHOSCorrNCells
! EMCAL vs PHOS, number of cells
TH2F * fhDeltaIANCells[2]
! Cluster "asymmetry" in cell units vs number of cells in cluster for E > 0.5, with and without match...
Int_t fNMaxCols
Number of EMCAL/PHOS rows.
TH2F * fhEtaPhiCell
! eta vs phi, cells
void SetExoECrossCuts(Int_t i, Float_t c)
Float_t fInvMassMinM02Cut
Minimum M02 shower shape cut value for clusters entering the invariant mass calculation.
Float_t GetInvMassMinECut() const
Double_t fTimeCutMin
Remove clusters/cells with time smaller than this value, in ns.
Bool_t fStudyExotic
Study the exotic cluster for different cuts, not QA related.
TH2F * fhTrackMatchedDPhiNegMod
! Phi distance between negative track and cluster vs module for E > 0.5 GeV
void InitParameters()
Initialize the parameters of the analysis.
TH2F * fhNCellsPerCluster
! N cells per cluster vs cluster energy vs eta of cluster
TH2F * fhNCellsPerClusterWeirdMod
! N cells per clusters different module, Reco, ridiculously large energy
TH2F * fhGenMCAccEtaPhi[4]
! eta vs phi of primary particle, in acceptance
TH2F * fhAmpMod
! Cell amplitude distribution for different module, Reco
TH2F * fhCaloV0SCorrNCells
! Calo vs V0 signal, number of cells
Float_t GetECross(Int_t absId, AliVCaloCells *cells, Float_t dtcut=10000)
TH2F ** fhIMMod
! cluster pairs invariant mass, different module,
TH2F ** fhNCellsPerClusterMod
! N cells per clusters different module, Reco
TH2F * fhTrackMatchedDEtaPos
! Eta distance between track and cluster vs cluster E, after and before photon cuts ...
TH2F * fhECellClusterLogRatio
! log (e cell / e cluster) vs e cluster
void SwitchOffFillAllTrackMatchingHistogram()
TH3F * fhEtaPhiAmpCell
! eta vs phi vs amplitude, cells
TH2F * fh1EleEOverPMod
! p/E for track-cluster matches, dR < 0.2, 60 < dEdx < 100, per SM
TH2F * fhCaloEvPNCells
! Calo vs event plane angle, number of cells
TH2F * fhBadClusterLambda0
! Cluster Lambda0 vs Energy, clusters declared bad
TH2F * fhClusterPairDiffTimeEPi0MassDCalSame
! DCal Cluster time TOF difference, for pairs in 0.1 < mass < 0.18, pairs in same Module ...
TH2F * fhRecoMCPhi[7][2]
! phi generated particle vs reconstructed phi
TH2F * fhExoL1NCell[10][5]
! Lambda1 vs n cells in cluster for several E cross cuts and cluster with E > 5
TH2F * fhTimeId
! Time vs Absolute cell Id
TH2F * fh1EOverP
! p/E for track-cluster matches
TH2F * fhIMvsOpAngle
! Cluster pairs opening angle vs mass
TH2F * fhRecoMCEta[7][2]
! eta generated particle vs reconstructed Eta
TH2F * fh2EledEdxMod
! dE/dx for electron candidates, per SM
TH2F * fhBadClusterLambda1
! Cluster Lambda1 vs Energy, clusters declared bad
TH2F * fhMCEle1EOverPR02
! p/E for track-cluster matches, dR < 0.2, MC electrons
TH2F * fhECellTotalRatio
! e cell / e total vs e total
Int_t fEMCALClusterNCellMin
Minimum number of cells on EMCal clusters.
TH1F * fhNClusters
! Number of clusters
TH2F * fhMCEle1EOverP
! p/E for track-cluster matches, MC electrons
void SetPHOSCellAmpMin(Float_t amp)
TH2F * fhExoECross[10][5]
! E cross for max cell in cluster, for different cuts
TH1F * fhMCEle1dR
! distance between projected track and cluster, MC electrons
TH2F * fhOpAngle
! Cluster pairs opening angle vs pair pT
TH1F * fhGenMCE[4]
! pt of primary particle
TH2F * fhIMDiff
! Cluster pairs invariant mass vs pair pT, for EMCAL or PHOS pairs
Float_t fExoDTimeCuts[5]
List of time cuts.
TH2F * fhEMR
! Electromagnetic distance to vertex vs rec energy
TH2F * fhClusterPairDiffTimeEPi0MassDCal
! DCal Cluster time TOF difference, for pairs in 0.1 < mass < 0.18
TH2F * fhExoL0ECross
! Lambda0 vs E cross fraction for clusters with E > 5 GeV
TH2F * fhCaloTrackMCorrNCells
! Calo vs V0 Track Multiplicity, number of cells
void ClusterAsymmetryHistograms(AliVCluster *clus, Int_t absIdMax, Bool_t goodCluster)
void SetConstantTimeShift(Float_t shift)
TH2F * fhIMDCALPHOSSame
! Cluster pairs invariant mass vs pair pT, for DCal-PHOS pairs
TH2F * fhDeltaCellClusterXNCells
! X cluster - X cell distribution (cm) vs N cells in cluster
TH2F * fhTrackMatchedDEtaDPhiPos
! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV, after and before ...
TH1F * fhNCells
! Number of towers/crystals with signal
TH2F * fhIMDCAL
! Cluster pairs invariant mass vs pair pT, for DCal pairs
TH2F * fhTrackMatchedDPhiPos
! Phi distance between track and cluster vs cluster E, after and before photon cuts ...
TH2F * fhCaloEvPNClusters
! Calo vs event plane angle, number of clusters
TH2F * fhDCALPHOSCorrNClusters
! DCAL vs PHOS, number of clusters
TH2F * fhXCellE
! X (cm) cell distribution vs cell energy
TH2F * fhBadClusterDeltaIA
! Cluster "asymmetry" in cell terms vs E, with and without matching; bad clusters.
TH2F * fhIMEMCALPHOS
! Cluster pairs invariant mass vs pair pT, for EMCAL(DCal eta acceptance)-EMCAL (PHOS eta acceptance)...
TH2F * fhHaR
! Hadron distance to vertex vs rec energy
Float_t GetInvMassMinM02Cut() const
TH2F * fhAmpWeirdMod
! Cell amplitude distribution for different module, very large Amp, Reco
TH2F * fhRE
! R=sqrt(x^2+y^2) (cm) cluster distribution vs cluster energy
void SetInvMassMaxTimeDifference(Float_t c)
TH2F * fhIM
! Cluster pairs invariant mass vs pair pT, for EMCAL or PHOS pairs
TH2F * fhExoTime[10][5]
! Time of exotic cluster, for different cuts
TH2F * fhTimeAmp
! Time vs Amplitude
void SwitchOnFillInvMassOpAngleHistogram()
void CalculateAverageTime(AliVCluster *clus, AliVCaloCells *cells, Double_t timeAverages[2])
TH2F * fhMCEle2MatchdEdx
! dE/dx vs. momentum for all matches, MC electrons
AliAnaCalorimeterQA & operator=(const AliAnaCalorimeterQA &qa)
Copy constructor not implemented.
TH1F * fhECharged
! E distribution, Reco, matched with track
TH2F * fhNClustersMod
! Number of clusters for different module, Reco
TH2F * fhHaVxyz
! Hadron production vertex
Float_t GetInvMassMaxTimeDifference() const
TH2F * fhAmpId
! Amplitude measured in towers/crystals vs id of tower.
TH2F * fhCaloV0SCorrECells
! Calo vs V0 signal, total measured cell energy
TH2F * fhDeltaIAL1[2]
! Cluster "asymmetry" in cell units vs Lambda1 for E > 0.5 GeV, n cells in cluster > 3...
void SetEMCALCellAmpMin(Float_t amp)
TH1F * fhAmplitude
! Amplitude measured in towers/crystals
TH2F * fhEMod
! Cluster E distribution for different module, Reco
TH2F * fhRNCells
! R=sqrt(x^2+y^2) (cm) cluster distribution vs N cells in cluster
Int_t fPHOSClusterNCellMin
Minimum number of cells on PHOS clusters.
TH2F * fhLambda0ForW0MC[12][5]
! L0 for different w0, depending on the particle of origin
TH2F * fhTrackMatchedDEtaNeg
! Eta distance between track and cluster vs cluster E, after and before photon cuts ...
void SetInvMassMaxECut(Float_t cut)
TH2F * fhBadClusterDeltaIEtaDeltaIPhiE6
! Difference between max cell index and farthest cell, eta vs phi, E > 6 GeV, with and without matchi...
TH2F * fhZCellE
! Z (cm) cell distribution vs cell energy
Float_t GetPHOSCellAmpMin() const
TH2F * fhGridCellsTime
! Cells ordered in column/row for different module, weighted with time, Reco
Float_t fExoNECrossCuts
Number of ecross cuts.
Double_t GetTimeCutMax() const
void SetExoDTimeCuts(Int_t i, Float_t c)
TH3F * fhXYZCell
! cell X vs Y vs Z (cm)
Float_t fInvMassMaxOpenAngle
Combine clusters within with a maximum opening angle between them. In radians.
TLorentzVector fPrimaryMomentum
! Primary MC momentum, temporary container
TH2F * fhCaloEvPEClusters
! Calo vs event plane angle, total measured cluster energy
Bool_t fCorrelate
Correlate PHOS/EMCAL cells/clusters, also with V0 and track multiplicity.
TH2F * fhDeltaIEtaDeltaIPhiE2[2]
! Difference between max cell index and farthest cell, eta vs phi, 2 < E < 6 GeV, with and without ma...
Bool_t fFillAllCellTimeHisto
Fill all cell time histo.
TH2F ** fhECellTotalLogRatioMod
! log (e cell / e total) vs e total, per SM
TH2F * fhTrackMatchedDEtaPosMod
! Eta distance between positive track and cluster vs module for E > 0.5 GeV
TH2F * fhDeltaCellClusterZNCells
! Z cluster - Z cell distribution (cm) vs N cells in cluster
void ClusterHistograms(AliVCluster *cluster, const TObjArray *caloClusters, AliVCaloCells *cells, Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Double_t tmax)
Bool_t fFillAllPosHisto2
Fill all the position related histograms 2.
TObjString * GetAnalysisCuts()
Save parameters used for analysis in a string.
TH2F * fhYCellE
! Y (cm) cell distribution vs cell energy
TH2F * fhCaloV0MCorrNClusters
! Calo vs V0 multiplicity , number of clusters
TH2F * fh1EleEOverP
! p/E for track-cluster matches, dR < 0.2, 60 < dEdx < 100
void SetInvMassMaxOpenAngle(Float_t c)
TH2F ** fhNCellsSumAmpPerMod
! N cells vs sum of amplitude in different modules, Reco
TH2F * fhDeltaCellClusterZE
! Z cluster - Z cell distribution (cm) vs cluster energy
TH2F * fhAmpIdLowGain
! Amplitude measured in towers/crystals vs id of tower, low gain towers
void CellInClusterPositionHistograms(AliVCluster *cluster)
Fill histograms releated to cluster cell position.
Bool_t ClusterMCHistograms(Bool_t matched, const Int_t *labels, Int_t nLabels, Int_t &pdg)
TH2F * fhGridCellsELowGain
! Cells ordered in column/row for different module, weighted with energy, Reco, low gain ...
TH2F * fh2MatchdEdx
! dE/dx vs. momentum for all matches
void BadClusterHistograms(AliVCluster *clus, const TObjArray *caloClusters, AliVCaloCells *cells, Int_t absIdMax, Double_t maxCellFraction, Float_t eCrossFrac, Double_t tmax)
TH2F * fhMCChHad1EOverP
! p/E for track-cluster matches, MC charged hadrons
TH2F * fhRecoMCDeltaEta[7][2]
! Gen-Reco eta generated particle vs reconstructed E
TH2F * fhMCChHad2MatchdEdx
! dE/dx vs. momentum for all matches, MC charged
void SwitchOffFillPi0PairDiffTimeHistogram()
TH2F * fhClusterMaxCellCloseCellDiff
! Difference between max cell energy and cell energy of the same cluster
TH2F * fhXE
! X (cm) cluster distribution vs cluster energy
TH2F * fhGridCells
! Cells ordered in column/row for different module, Reco
TH2F * fhGridCellsTimeLowGain
! Cells ordered in column/row for different module, weighted with time, Reco, low gain ...
TH2F * fhCaloTrackMCorrECells
! Calo vs V0 Track Multipliticy, total measured cell energy
Float_t fInvMassMaxM02Cut
Maximum M02 shower shape cut value for clusters entering the invariant mass calculation.
void SetInvMassMinECut(Float_t cut)
TH2F * fhEMVxyz
! Electromagnetic particle production vertex
TH2F * fhLambda0
! Cluster Lambda0 vs Energy
void CellHistograms(AliVCaloCells *cells)
Fill histograms related to cells only.
TH2F * fhCaloV0MCorrECells
! Calo vs V0 multiplicity, total measured cell energy
TH2F * fhDeltaIA[2]
! Cluster "asymmetry" in cell terms vs E, with and without matching
void SwitchOnFillInvMassReducedEMCALHistogram()
TH2F * fhEMCALDCALCorrNCells
! EMCAL vs DCAL, number of cells
TH2F * fhMCNeutral1EleEOverP
! p/E for track-cluster matches, dR < 0.2, 60 < dEdx < 100, MC neutral
TH2F * fhClusterPairDiffTimeEPi0MassSame
! EMCal/PHOS Cluster time TOF difference, for pairs in 0.1 < mass < 0.18, pairs in same Module ...
Int_t fNMaxRows
Number of EMCAL/PHOS columns.
TH2F * fhDeltaIAMC[4]
! Cluster "asymmetry" in cell terms vs E, from MC photon, electron, conversion or hadron...
TH2F * fhEMaxCellClusterRatio
! e max cell / e cluster vs e cluster
Float_t fExoNDTimeCuts
Number of time cuts.
TH2F * fhYE
! Y (cm) cluster distribution vs cluster energy
TH1F * fhMCNeutral1dR
! Distance between projected track and cluster, MC neutral
TH2F * fhBadClusterTimeEnergy
! Time Max cell of bad cluster
TH2F * fhCaloV0SCorrEClusters
! Calo vs V0 signal, total measured cluster energy
TH2F * fhDCALPHOSCorrEClusters
! DCAL vs PHOS, total measured cluster energy
Bool_t fStudyBadClusters
Study bad clusters not passing selection criteria (exotic, shower shape, n cells).
TH2F * fhCaloCenNCells
! Calo vs centrality, number of cells
TH2F * fhYNCells
! Y (cm) cluster distribution vs N cells in cluster
TH2F * fhCaloV0MCorrNCells
! Calo vs V0 multiplicity, number of cells
TH2F * fhMCNeutral1EOverP
! p/E for track-cluster matches, MC neutral
TH2F * fhMCChHad1EOverPR02
! p/E for track-cluster matches, dR < 0.2, MC charged hadrons
Float_t GetInvMassMaxM02Cut() const
TH2F * fhGridCellsLowGain
! Cells ordered in column/row for different module, Reco, low gain
TH2F * fhAsym
! Cluster pairs invariant mass vs pair pT
TH2F * fhRecoMCE[7][2]
! E generated particle vs reconstructed E
Int_t GetPHOSNCellsPerClusterMin() const
TH2F * fhTrackMatchedDEtaNegMod
! Eta distance between negative track and cluster vs module for E > 0.5 GeV
TH2F * fhEMCALPHOSCorrECells
! EMCAL vs PHOS, total measured cell energy
TH2F * fhDeltaIEtaDeltaIPhiE0[2]
! Difference between max cell index and farthest cell, eta vs phi, E < 2 GeV, with and without matchi...
void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
TH2F * fhTimeAmpLowGain
! Time vs Amplitude, low gain
void ClusterLoopHistograms(const TObjArray *clusters, AliVCaloCells *cells)
TH2F * fhDCALPHOSCorrECells
! DCAL vs PHOS, total measured cell energy
TH2F * fhEMCALDCALCorrECells
! EMCAL vs DCAL, total measured cell energy
Float_t fInvMassMaxTimeDifference
Maximum difference between the time of the 2 clusters to be considered in invariant mass...
TH2F * fhEWeirdMod
! Cluster E distribution for different module, very large E, Reco
TH2F * fhCellTimeSpreadRespectToCellMax
! Difference of the time of cell with maximum dep energy and the rest of cells
Bool_t fFillAllTMHisto
Fill track matching histograms.
TH2F * fhSumClustersEnergyMod
! Sum of clusters energy for different module, Reco
TH2F * fhECellClusterRatio
! e cell / e cluster vs e cluster
TH2F * fhExoL1ECross
! Lambda1 vs E cross fraction for clusters with E > 5 GeV
TH2F * fhDeltaIAL0[2]
! Cluster "asymmetry" in cell units vs Lambda0 for E > 0.5 GeV, n cells in cluster > 3...
TH2F * fhCaloV0SCorrNClusters
! Calo vs V0 signal , number of clusters
TH2F * fhGenMCEtaPhi[4]
! eta vs phi of primary particle
TH2F * fhLambda1
! Cluster Lambda1 vs Energy
Float_t fConstantTimeShift
Apply a 600 ns time shift in case of simulation, shift in ns.