AliPhysics  v5-06-40-01 (42bb456)
 All Classes Namespaces Files Functions Variables Enumerations Enumerator Macros
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  // Setters and getters
106 
107  Float_t GetEMCALCellAmpMin() const { return fEMCALCellAmpMin ; }
108  void SetEMCALCellAmpMin(Float_t amp) { fEMCALCellAmpMin = amp ; }
109 
110  Float_t GetPHOSCellAmpMin() const { return fPHOSCellAmpMin ; }
111  void SetPHOSCellAmpMin (Float_t amp) { fPHOSCellAmpMin = amp ; }
112 
115 
118 
119  Float_t GetInvMassMinECut() const { return fInvMassMinECut ; }
120  void SetInvMassMinECut(Float_t cut) { fInvMassMinECut = cut ; }
121 
122  Float_t GetInvMassMaxECut() const { return fInvMassMaxECut ; }
123  void SetInvMassMaxECut(Float_t cut) { fInvMassMaxECut = cut ; }
124 
125  Float_t GetInvMassMinM02Cut() const { return fInvMassMinM02Cut ; }
126  void SetInvMassMinM02Cut(Float_t cut) { fInvMassMinM02Cut = cut ; }
127 
128  Float_t GetInvMassMaxM02Cut() const { return fInvMassMaxM02Cut ; }
129  void SetInvMassMaxM02Cut(Float_t cut) { fInvMassMaxM02Cut = cut ; }
130 
131  Float_t GetInvMassMaxOpenAngle() const { return fInvMassMaxOpenAngle; }
133 
136 
137  Double_t GetTimeCutMin() const { return fTimeCutMin ; }
138  Double_t GetTimeCutMax() const { return fTimeCutMax ; }
139  void SetTimeCut(Double_t min, Double_t max) {
140  fTimeCutMin = min ; fTimeCutMax = max ; }
141 
142  // Histogram switchs
143 
146 
149 
152 
155 
158 
161 
164 
167 
170 
171  void SwitchOnCorrelation() { fCorrelate = kTRUE ; }
172  void SwitchOffCorrelation() { fCorrelate = kFALSE ; }
173 
176 
177  // Analysis not to be used in QA
178 
181 
182  void SwitchOnStudyWeight() { fStudyWeight = kTRUE ; }
183  void SwitchOffStudyWeight() { fStudyWeight = kFALSE ; }
184 
185  void SwitchOnStudyExotic() { fStudyExotic = kTRUE ; }
186  void SwitchOffStudyExotic() { fStudyExotic = kFALSE ; }
187 
188  void SetNECrossCuts(Int_t n) { fExoNECrossCuts = n ; }
189  void SetNDTimeCuts (Int_t n) { fExoNDTimeCuts = n ; }
190 
191  void SetExoECrossCuts (Int_t i, Float_t c) { if (i<fExoNECrossCuts) fExoECrossCuts[i] = c ; }
192  void SetExoDTimeCuts (Int_t i, Float_t c) { if (i<fExoNDTimeCuts ) fExoDTimeCuts [i] = c ; }
193 
194  private:
195 
196  // Switches
197 
201  Bool_t fFillAllTH3 ;
202  Bool_t fFillAllTMHisto ;
203 
208 
209  Bool_t fCorrelate ;
211 
212  // Analysis not to be used in QA
213 
215  Bool_t fStudyExotic;
216  Bool_t fStudyWeight;
217 
218  // Parameters
219 
220  Int_t fNModules ;
221  Int_t fNRCU ;
222  Int_t fNMaxCols ;
223  Int_t fNMaxRows ;
224 
225  // Cuts
226 
227  Double_t fTimeCutMin ;
228  Double_t fTimeCutMax ;
229  Float_t fCellAmpMin;
231  Float_t fPHOSCellAmpMin ;
234 
235  // Invariant mass analysis
236 
237  Float_t fInvMassMinECut;
238  Float_t fInvMassMaxECut;
243  // Exotic studies
244 
245  Float_t fExoNECrossCuts ;
246  Float_t fExoECrossCuts[10];
247  Float_t fExoNDTimeCuts ;
248  Float_t fExoDTimeCuts[5] ;
249 
250  TLorentzVector fClusterMomentum;
251  TLorentzVector fClusterMomentum2;
252  TLorentzVector fPrimaryMomentum;
253 
254  // Calorimeter Clusters
255 
256  TH1F * fhE ;
257  TH1F * fhPt ;
258  TH1F * fhPhi;
259  TH1F * fhEta;
260  TH2F * fhEtaPhi;
261  TH3F * fhEtaPhiE ;
262  TH1F * fhECharged ;
263  TH1F * fhPtCharged ;
264  TH1F * fhPhiCharged;
265  TH1F * fhEtaCharged;
268 
269  TH2F * fhIM;
270  TH2F * fhIMSame;
271  TH2F * fhIMDiff;
272  TH2F * fhIMDCAL;
273  TH2F * fhIMDCALSame;
274  TH2F * fhIMDCALDiff;
275  TH2F * fhIMDCALPHOS;
277  TH2F * fhIMEMCALPHOS;
279 
280  TH2F * fhAsym;
281 
282  TH2F* fhOpAngle;
284 
287 
290 
293 
294  TH1F * fhNClusters;
295 
301 
306 
307 // TH2F * fhClusterMaxCellDiffAverageTime; //!<! Difference between cluster average time and time of cell with more energy
308 // TH2F * fhClusterMaxCellDiffWeightedTime; //!<! Difference between cluster weighted time and time of cell with more energy
310 
311  TH2F * fhLambda0;
312  TH2F * fhLambda1;
313 // TH2F * fhDispersion; //!<! Cluster Dispersion vs Energy
314 
315  // Bad clusters histograms
316 
322 
326 
330 
333 
334  // Cluster cell size
335 
339  TH2F * fhDeltaIA[2];
340  TH2F * fhDeltaIAL0[2];
341  TH2F * fhDeltaIAL1[2];
342  TH2F * fhDeltaIANCells[2] ;
343  TH2F * fhDeltaIAMC[4];
348 
349  // Cluster/cell Position
350 
351  TH2F * fhRNCells ;
352  TH2F * fhXNCells ;
353  TH2F * fhYNCells ;
354  TH2F * fhZNCells ;
355 
356  TH2F * fhRE ;
357  TH2F * fhXE ;
358  TH2F * fhYE ;
359  TH2F * fhZE ;
360  TH3F * fhXYZ;
361 
362  TH2F * fhRCellE ;
363  TH2F * fhXCellE ;
364  TH2F * fhYCellE ;
365  TH2F * fhZCellE ;
366  TH3F * fhXYZCell;
367 
372 
377 
378  // Calorimeter cells
379 
380  TH1F * fhNCells;
382  TH1F * fhAmplitude;
383  TH2F * fhAmpId;
385  TH2F * fhEtaPhiCell;
386 
387  TH1F * fhTime;
388  //TH2F * fhTimeVz; //!<! Time measured in towers/crystals vs vertex z component, for E > 0.5
389  TH2F * fhTimeId;
390  TH2F * fhTimeAmp;
391 
392  TH2F * fhAmpIdLowGain;
395 
396  TH2F * fhCellECross;
397 
398  // Calorimeters Correlation
399 
404 
409 
414 
415  // V0 Correlation
416 
425 
426  // Track Correlation
427 
432 
433  // Centrality
434 
439 
440  // Event plane
441 
446 
447  // Module histograms
448 
449  TH2F * fhEMod ;
450  TH2F * fhAmpMod ;
451  TH2F * fhEWeirdMod ;
452  TH2F * fhAmpWeirdMod ;
453  TH2F * fhTimeMod ;
454  TH2F * fhNClustersMod ;
455  TH2F * fhNCellsMod ;
464  TH2F * fhGridCells ;
465  TH2F * fhGridCellsE ;
466  TH2F * fhGridCellsTime ;
470  TH2F ** fhTimeAmpPerRCU;
471  TH2F ** fhIMMod;
472 
473  // Weight studies
474 
479 
480  TH2F * fhLambda0ForW0[14];
481 //TH2F * fhLambda1ForW0[7]; //!<! L1 for 7 defined w0= 3, 3.5 ... 6
482 
483  TH2F * fhLambda0ForW0MC[14][5];
484 //TH2F * fhLambda1ForW0MC[7][5]; //!<! L1 for 7 defined w0= 3, 3.5 ... 6, depending on the particle of origin
485 
490 
491  // Exotic studies
492 
493  TH2F * fhExoNCell [10][5] ;
494  TH2F * fhExoL0 [10][5] ;
495  TH2F * fhExoL1 [10][5] ;
496  TH2F * fhExoECross [10][5] ;
497  TH2F * fhExoTime [10][5] ;
498  TH2F * fhExoDTime [10] ;
499  TH2F * fhExoL0NCell[10][5] ;
500  TH2F * fhExoL0ECross ;
501  TH2F * fhExoL1NCell[10][5] ;
502  TH2F * fhExoL1ECross ;
503 
504  // Pure MC histograms
505 
507  enum mcTypes { kmcPhoton = 0, kmcPi0 = 1, kmcEta = 2,
510 
511  TH2F * fhRecoMCE[7][2] ;
512  TH2F * fhRecoMCPhi[7][2] ;
513  TH2F * fhRecoMCEta[7][2] ;
514  TH2F * fhRecoMCDeltaE[7][2] ;
515  TH2F * fhRecoMCRatioE[7][2] ;
516  TH2F * fhRecoMCDeltaPhi[7][2];
517  TH2F * fhRecoMCDeltaEta[7][2];
518 
519  TH1F * fhGenMCE [4] ;
520  TH1F * fhGenMCPt[4] ;
521  TH2F * fhGenMCEtaPhi[4] ;
522  TH1F * fhGenMCAccE [4] ;
523  TH1F * fhGenMCAccPt[4] ;
524  TH2F * fhGenMCAccEtaPhi[4] ;
525 
526  TH2F * fhEMVxyz ;
527  TH2F * fhEMR ;
528  TH2F * fhHaVxyz ;
529  TH2F * fhHaR ;
530 
531  // Histograms for MC track-matching
532 
533  TH2F * fh1EOverP;
534  TH2F * fh2dR;
535  TH2F * fh2EledEdx;
536  TH2F * fh2MatchdEdx;
537  TH2F * fh1EOverPR02;
538  TH2F * fh1EleEOverP;
539 
540  TH2F * fh1EOverPMod;
541  TH2F * fh2dRMod;
542  TH2F * fh2EledEdxMod;
546 
547  TH2F * fhMCEle1EOverP;
548  TH1F * fhMCEle1dR;
550 
552  TH1F * fhMCChHad1dR;
554 
556  TH1F * fhMCNeutral1dR;
558 
562 
566 
570 
574 
579 
582 
585 
587  ClassDef(AliAnaCalorimeterQA,32) ;
589 
590 } ;
591 
592 #endif //ALIANACALORIMETERQA_H
593 
594 
595 
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 * 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 * fhLambda0ForW0MC[14][5]
! L0 for 7 defined w0= 3, 3.5 ... 6, depending on the particle of origin
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()
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)
Calculate cluster weights.
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.
TH2F * fhLambda0ForW0[14]
! L0 for 7 defined w0= 3, 3.5 ... 6
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)
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 * 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