1 #ifndef ALIANALYSISTASKEMCALTRIGGERQA_H
2 #define ALIANALYSISTASKEMCALTRIGGERQA_H
11 class AliEMCALGeometry;
15 class AliEMCALRecoUtils;
16 #include "AliEMCALGeoParams.h"
17 #include "AliAnalysisTaskSE.h"
19 class AliAnalysisTaskEMCALTriggerQA : public AliAnalysisTaskSE
22 AliAnalysisTaskEMCALTriggerQA(); // default constructor
24 AliAnalysisTaskEMCALTriggerQA(const char *name); // named constructor
26 virtual ~AliAnalysisTaskEMCALTriggerQA() { ; } // destructor
30 void LocalInit() { Init() ; }
32 void UserCreateOutputObjects();
34 void UserExec(Option_t *option);
36 AliEMCALRecoUtils* GetRecoUtils() { if(!fRecoUtils) fRecoUtils = new AliEMCALRecoUtils ;
37 else return fRecoUtils ; }
39 // OADB and geometry settings
43 void SetGeometryName(TString name) { fGeoName = name ; }
47 void SwitchOnEMCALOADB() { fAccessOADB = kTRUE ; }
48 void SwitchOffEMCALOADB() { fAccessOADB = kFALSE ; }
50 void SetOADBFilePath(TString path) { fOADBFilePath = path ; }
54 void SetTRUTotalSignalHistogramsRange(Int_t nbins, Float_t max) { fNBinsTRUSignal = nbins; fMaxTRUSignal = max ; }
55 void SetSTUTotalSignalHistogramsRange(Int_t nbins, Float_t max) { fNBinsSTUSignal = nbins; fMaxSTUSignal = max ; }
56 void SetV0TotalSignalHistogramsRange (Int_t nbins, Float_t max) { fNBinsV0Signal = nbins; fMaxV0Signal = max ; }
57 void SetSTUFEERatioHistogramsRange (Int_t nbins, Float_t max) { fNBinsSTUFEERatio = nbins; fMaxSTUFEERatio = max ; }
58 void SetSTUTRURatioHistogramsRange (Int_t nbins, Float_t max) { fNBinsSTUTRURatio = nbins; fMaxSTUFEERatio = max ; }
59 void SetClusterEHistogramsRange (Int_t nbins, Float_t max) { fNBinsClusterE = nbins; fMaxClusterE = max ; }
62 TList *fOutputList; //! Output list
64 AliEMCALRecoUtils *fRecoUtils; // RecoUtils
66 Bool_t fGeoSet ; // Geometry already set
67 AliEMCALGeometry *fGeometry; // Access to EMCAL geometry utils
68 TString fGeoName; // Name of geometry used
70 Bool_t fOADBSet ; // AODB parameters already set
71 Bool_t fAccessOADB ; // Get calibration from OADB for EMCAL
72 TString fOADBFilePath ; // Default path $ALICE_ROOT/OADB/EMCAL, if needed change
74 TH1F *fhNEvents; //! Number of selected events
75 TH2F *fhFORAmp; //! FEE cells deposited energy, grouped like FastOR 2x2 per Row and Column
76 TH2F *fhFORAmpL1G; //! FEE cells deposited energy, grouped like FastOR 2x2 per Row and Column, with L1 Gamma trigger event
77 TH2F *fhFORAmpL1J; //! FEE cells deposited energy, grouped like FastOR 2x2 per Row and Column, with L1 Jet trigger event
78 TH2F *fhL0Amp; //! FALTRO signal per Row and Column for FOR involves L0 patch
79 TH2F *fhL0AmpL1G; //! FALTRO signal per Row and Column for FOR involves L0 patch, with L1G trigger event
80 TH2F *fhL0AmpL1J; //! FALTRO signal per Row and Column for FOR involves L0 patch, with L1J trigger event
81 TH2F *fhL1Amp; //! STU signal per Row and Column for FOR involves L0 patch
82 TH2F *fhL1GAmp; //! STU signal per Row and Column for FOR position of L1 Gamma patch (top-left)
83 TH2F *fhL1JAmp; //! STU signal per Row and Column for FOR position of L1 Jet patch (top-left)
84 TH2F *fhL0Patch; //! FOR with L0 patch associated
85 TH2F *fhL1GPatch; //! FOR with L1 Gamma patch associated
86 TH2F *fhL1JPatch; //! FOR with L1 Jet patch associated
87 TH2F *fhFEESTU; //! Correlation FEE vs STU
88 TH2F *fhTRUSTU; //! Correlation TRU vs STU
89 TH2I *fhV0STU; //! Total signal STU vs V0C+V0S
91 TH2F *fhGPMaxVV0TT; //! V0 signal vs maximum gamma L1 patch
92 TH2F *fhJPMaxVV0TT; //! V0 signal vs maximum jet L1 patch
93 TProfile2D *fhFORMeanAmp; //! Mean FastOR(FEE) signal per Row and Column
94 TProfile2D *fhL0MeanAmp; //! Mean FastOR(TRU) signal per Row and Column
95 TProfile2D *fhL1MeanAmp; //! Mean FastOR(STU) signal per Row and Column
96 TH1F *fhV0MB; //! V0 distribution for MB triggered event
97 TH1F *fhV0L1G; //! V0 distribution for L1G triggered event
98 TH1F *fhV0L1J; //! V0 distribution for L1J triggered event
99 TH2F *fhL1GPatchMax; //! FOR of max. amplitude patch with L1 Gamma patch associated
100 TH2F *fhL1JPatchMax; //! FOR of max. amplitude patch with L1 Jet patch associated
102 // Cluster vs trigger histograms
104 TH1F *fhClusMB; //! Clusters distribution for MB trigger
105 TH1F *fhClusMBPure; //! Clusters distribution for MB trigger
106 TH1F *fhClusL0; //! Clusters distribution for L0 trigger
107 TH1F *fhClusL1G; //! Clusters distribution for L1G trigger
108 TH1F *fhClusL1J; //! Clusters distribution for L1J trigger
109 TH1F *fhClusL1GOnly; //! Clusters distribution for L1G trigger and not L1J
110 TH1F *fhClusL1JOnly; //! Clusters distribution for L1J trigger and not L1G
111 TH1F *fhClusMaxMB; //! Maximum E Cluster per event distribution for MB trigger
112 TH1F *fhClusMaxMBPure; //! Maximum E Cluster per event distribution for MB trigger
113 TH1F *fhClusMaxL0; //! Maximum E Cluster per event distribution for L0 trigger
114 TH1F *fhClusMaxL1G; //! Maximum E Cluster per event distribution for L1G trigger
115 TH1F *fhClusMaxL1J; //! Maximum E Cluster per event distribution for L1J trigger
116 TH1F *fhClusMaxL1GOnly; //! Maximum E Cluster per event distribution for L1G trigger and not L1J
117 TH1F *fhClusMaxL1JOnly; //! Maximum E Cluster per event distribution for L1J trigger and not L1G
119 TH2F *fhClusCenMB; //! Clusters Centrality vs E distribution for MB trigger
120 TH2F *fhClusCenL0; //! Clusters Centrality vs E distribution for L0 trigger
121 TH2F *fhClusCenL1G; //! Clusters Centrality vs E distribution for L1G trigger
122 TH2F *fhClusCenL1J; //! Clusters Centrality vs E distribution for L1J trigger
123 TH2F *fhClusCenL1GOnly; //! Clusters Centrality vs E distribution for L1G trigger and not L1J
124 TH2F *fhClusCenL1JOnly; //! Clusters Centrality vs E distribution for L1J trigger and not L1G
125 TH2F *fhClusCenMaxMB; //! Maximum E Cluster vs Centrality per event distribution for MB trigger
126 TH2F *fhClusCenMaxL0; //! Maximum E Cluster vs Centrality per event distribution for L0 trigger
127 TH2F *fhClusCenMaxL1G; //! Maximum E Cluster vs Centrality per event distribution for L1G trigger
128 TH2F *fhClusCenMaxL1J; //! Maximum E Cluster vs Centrality per event distribution for L1J trigger
129 TH2F *fhClusCenMaxL1GOnly; //! Maximum E Cluster vs Centrality per event distribution for L1G trigger and not L1J
130 TH2F *fhClusCenMaxL1JOnly; //! Maximum E Cluster vs Centrality per event distribution for L1J trigger and not L1G
132 TH2F *fhClusV0MB; //! Clusters Centrality vs E distribution for MB trigger
133 TH2F *fhClusV0L0; //! Clusters Centrality vs E distribution for L0 trigger
134 TH2F *fhClusV0L1G; //! Clusters Centrality vs E distribution for L1G trigger
135 TH2F *fhClusV0L1J; //! Clusters Centrality vs E distribution for L1J trigger
136 TH2F *fhClusV0L1GOnly; //! Clusters Centrality vs E distribution for L1G trigger and not L1J
137 TH2F *fhClusV0L1JOnly; //! Clusters Centrality vs E distribution for L1J trigger and not L1G
138 TH2F *fhClusV0MaxMB; //! Maximum E Cluster vs Centrality per event distribution for MB trigger
139 TH2F *fhClusV0MaxL0; //! Maximum E Cluster vs Centrality per event distribution for L0 trigger
140 TH2F *fhClusV0MaxL1G; //! Maximum E Cluster vs Centrality per event distribution for L1G trigger
141 TH2F *fhClusV0MaxL1J; //! Maximum E Cluster vs Centrality per event distribution for L1J trigger
142 TH2F *fhClusV0MaxL1GOnly; //! Maximum E Cluster vs Centrality per event distribution for L1G trigger and not L1J
143 TH2F *fhClusV0MaxL1JOnly; //! Maximum E Cluster vs Centrality per event distribution for L1J trigger and not L1G
145 TH2F *fhClusEtaMB; //! Clusters eta vs E distribution for MB trigger
146 TH2F *fhClusEtaL0; //! Clusters eta vs E distribution for L0 trigger
147 TH2F *fhClusEtaL1G; //! Clusters eta vs E distribution for L1G trigger
148 TH2F *fhClusEtaL1J; //! Clusters eta vs E distribution for L1J trigger
149 TH2F *fhClusEtaL1GOnly; //! Clusters eta vs E distribution for L1G trigger and not L1J
150 TH2F *fhClusEtaL1JOnly; //! Clusters eta vs E distribution for L1J trigger and not L1G
151 TH2F *fhClusEtaMaxMB; //! Maximum E Cluster vs Eta per event distribution for MB trigger
152 TH2F *fhClusEtaMaxL0; //! Maximum E Cluster vs Eta per event distribution for L0 trigger
153 TH2F *fhClusEtaMaxL1G; //! Maximum E Cluster vs Eta per event distribution for L1G trigger
154 TH2F *fhClusEtaMaxL1J; //! Maximum E Cluster vs Eta per event distribution for L1J trigger
155 TH2F *fhClusEtaMaxL1GOnly; //! Maximum E Cluster vs Eta per event distribution for L1G trigger and not L1J
156 TH2F *fhClusEtaMaxL1JOnly; //! Maximum E Cluster vs Eta per event distribution for L1J trigger and not L1G
158 TH2F *fhClusPhiMB; //! Clusters Phi vs E distribution for MB trigger
159 TH2F *fhClusPhiL0; //! Clusters Phi vs E distribution for L0 trigger
160 TH2F *fhClusPhiL1G; //! Clusters Phi vs E distribution for L1G trigger
161 TH2F *fhClusPhiL1J; //! Clusters Phi vs E distribution for L1J trigger
162 TH2F *fhClusPhiL1GOnly; //! Clusters Phi vs E distribution for L1G trigger and not L1J
163 TH2F *fhClusPhiL1JOnly; //! Clusters Phi vs E distribution for L1J trigger and not L1G
164 TH2F *fhClusPhiMaxMB; //! Maximum E Cluster vs Phi per event distribution for MB trigger
165 TH2F *fhClusPhiMaxL0; //! Maximum E Cluster vs Phi per event distribution for L0 trigger
166 TH2F *fhClusPhiMaxL1G; //! Maximum E Cluster vs Phi per event distribution for L1G trigger
167 TH2F *fhClusPhiMaxL1J; //! Maximum E Cluster vs Phi per event distribution for L1J trigger
168 TH2F *fhClusPhiMaxL1GOnly; //! Maximum E Cluster vs Phi per event distribution for L1G trigger and not L1J
169 TH2F *fhClusPhiMaxL1JOnly; //! Maximum E Cluster vs Phi per event distribution for L1J trigger and not L1G
171 TH2F *fhClusEtaPhiHighMB; //! Clusters eta vs phi distribution for MB trigger, energy above 10 GeV
172 TH2F *fhClusEtaPhiHighL0; //! Clusters eta vs phi distribution for L0 trigger, energy above 10 GeV
173 TH2F *fhClusEtaPhiHighL1G; //! Clusters eta vs phi distribution for L1G trigger, energy above 10 GeV
174 TH2F *fhClusEtaPhiHighL1J; //! Clusters eta vs phi distribution for L1J trigger, energy above 10 GeV
175 TH2F *fhClusEtaPhiHighL1GOnly; //! Clusters eta vs phi distribution for L1G trigger and not L1J, energy above 10 GeV
176 TH2F *fhClusEtaPhiHighL1JOnly; //! Clusters eta vs phi distribution for L1J trigger and not L1G, energy above 10 GeV
177 TH2F *fhClusEtaPhiHighCluMaxMB; //! Maximum E Cluster, Phi vs Eta per event distribution for MB trigger, energy above 10 GeV
178 TH2F *fhClusEtaPhiHighCluMaxL0; //! Maximum E Cluster, Phi vs Eta per event distribution for L0 trigger, energy above 10 GeV
179 TH2F *fhClusEtaPhiHighCluMaxL1G; //! Maximum E Cluster, Phi vs Eta per event distribution for L1G trigger, energy above 10 GeV
180 TH2F *fhClusEtaPhiHighCluMaxL1J; //! Maximum E Cluster, Phi vs Eta per event distribution for L1J trigger, energy above 10 GeV
181 TH2F *fhClusEtaPhiHighCluMaxL1GOnly; //! Maximum E Cluster, Phi vs Eta per event distribution for L1G trigger and not L1J, energy above 10 GeV
182 TH2F *fhClusEtaPhiHighCluMaxL1JOnly; //! Maximum E Cluster, Phi vs Eta per event distribution for L1J trigger and not L1G, energy above 10 GeV
184 TH2F *fhClusEtaPhiHighCellMaxMB; //! Clusters maximum energy cell index eta vs phi distribution for MB trigger, energy above 10 GeV
185 TH2F *fhClusEtaPhiHighCellMaxL0; //! Clusters maximum energy cell index eta vs phi distribution for L0 trigger, energy above 10 GeV
186 TH2F *fhClusEtaPhiHighCellMaxL1G; //! Clusters maximum energy cell index eta vs phi distribution for L1G trigger, energy above 10 GeV
187 TH2F *fhClusEtaPhiHighCellMaxL1J; //! Clusters maximum energy cell index eta vs phi distribution for L1J trigger, energy above 10 GeV
188 TH2F *fhClusEtaPhiHighCellMaxL1GOnly; //! Clusters maximum energy cell index eta vs phi distribution for L1G trigger and not L1J, energy above 10 GeV
189 TH2F *fhClusEtaPhiHighCellMaxL1JOnly; //! Clusters maximum energy cell index eta vs phi distribution for L1J trigger and not L1G, energy above 10 GeV
190 TH2F *fhClusEtaPhiHighCellMaxCluMaxMB; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for MB trigger, energy above 10 GeV
191 TH2F *fhClusEtaPhiHighCellMaxCluMaxL0; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L0 trigger, energy above 10 GeV
192 TH2F *fhClusEtaPhiHighCellMaxCluMaxL1G; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1G trigger, energy above 10 GeV
193 TH2F *fhClusEtaPhiHighCellMaxCluMaxL1J; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1J trigger, energy above 10 GeV
194 TH2F *fhClusEtaPhiHighCellMaxCluMaxL1GOnly; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1G trigger and not L1J, energy above 10 GeV
195 TH2F *fhClusEtaPhiHighCellMaxCluMaxL1JOnly; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1J trigger and not L1G, energy above 10 GeV
197 TH2F *fhClusEtaPhiLowMB; //! Clusters eta vs phi distribution for MB trigger, energy below 10 GeV
198 TH2F *fhClusEtaPhiLowL0; //! Clusters eta vs phi distribution for L0 trigger, energy below 10 GeV
199 TH2F *fhClusEtaPhiLowL1G; //! Clusters eta vs phi distribution for L1G trigger, energy below 10 GeV
200 TH2F *fhClusEtaPhiLowL1J; //! Clusters eta vs phi distribution for L1J trigger, energy below 10 GeV
201 TH2F *fhClusEtaPhiLowL1GOnly; //! Clusters eta vs phi distribution for L1G trigger and not L1J, energy below 10 GeV
202 TH2F *fhClusEtaPhiLowL1JOnly; //! Clusters eta vs phi distribution for L1J trigger and not L1G, energy below 10 GeV
203 TH2F *fhClusEtaPhiLowCluMaxMB; //! Maximum E Cluster, Phi vs Eta per event distribution for MB trigger, energy below 10 GeV
204 TH2F *fhClusEtaPhiLowCluMaxL0; //! Maximum E Cluster, Phi vs Eta per event distribution for L0 trigger, energy below 10 GeV
205 TH2F *fhClusEtaPhiLowCluMaxL1G; //! Maximum E Cluster, Phi vs Eta per event distribution for L1G trigger, energy below 10 GeV
206 TH2F *fhClusEtaPhiLowCluMaxL1J; //! Maximum E Cluster, Phi vs Eta per event distribution for L1J trigger, energy below 10 GeV
207 TH2F *fhClusEtaPhiLowCluMaxL1GOnly; //! Maximum E Cluster, Phi vs Eta per event distribution for L1G trigger and not L1J, energy below 10 GeV
208 TH2F *fhClusEtaPhiLowCluMaxL1JOnly; //! Maximum E Cluster, Phi vs Eta per event distribution for L1J trigger and not L1G, energy below 10 GeV
210 TH2F *fhClusEtaPhiLowCellMaxMB; //! Clusters maximum energy cell index eta vs phi distribution for MB trigger, energy below 10 GeV
211 TH2F *fhClusEtaPhiLowCellMaxL0; //! Clusters maximum energy cell index eta vs phi distribution for L0 trigger, energy below 10 GeV
212 TH2F *fhClusEtaPhiLowCellMaxL1G; //! Clusters maximum energy cell index eta vs phi distribution for L1G trigger, energy below 10 GeV
213 TH2F *fhClusEtaPhiLowCellMaxL1J; //! Clusters maximum energy cell index eta vs phi distribution for L1J trigger, energy below 10 GeV
214 TH2F *fhClusEtaPhiLowCellMaxL1GOnly; //! Clusters maximum energy cell index eta vs phi distribution for L1G trigger and not L1J, energy below 10 GeV
215 TH2F *fhClusEtaPhiLowCellMaxL1JOnly; //! Clusters maximum energy cell index eta vs phi distribution for L1J trigger and not L1G, energy below 10 GeV
216 TH2F *fhClusEtaPhiLowCellMaxCluMaxMB; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for MB trigger, energy below 10 GeV
217 TH2F *fhClusEtaPhiLowCellMaxCluMaxL0; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L0 trigger, energy below 10 GeV
218 TH2F *fhClusEtaPhiLowCellMaxCluMaxL1G; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1G trigger, energy below 10 GeV
219 TH2F *fhClusEtaPhiLowCellMaxCluMaxL1J; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1J trigger, energy below 10 GeV
220 TH2F *fhClusEtaPhiLowCellMaxCluMaxL1GOnly; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1G trigger and not L1J, energy below 10 GeV
221 TH2F *fhClusEtaPhiLowCellMaxCluMaxL1JOnly; //! Maximum E Cluster, maximum energy cell index Phi vs Eta per event distribution for L1J trigger and not L1G, energy below 10 GeV
227 Int_t fNBinsSTUSignal ; // Number of bins for STU total signal histograms
228 Float_t fMaxSTUSignal ; // Maximum value for TRU total signal histograms
229 Int_t fNBinsTRUSignal ; // Number of bins for TRU total signal histograms
230 Float_t fMaxTRUSignal ; // Maximum value for TRU total signal histograms
231 Int_t fNBinsV0Signal ; // Number of bins for V0 total signal histograms
232 Float_t fMaxV0Signal ; // Maximum value for V0 total signal histograms
233 Int_t fNBinsSTUFEERatio ; // Number of bins for STU/FEE ratios histograms
234 Float_t fMaxSTUFEERatio ; // Maximum value for STU/FEE ratios histograms
235 Int_t fNBinsSTUTRURatio ; // Number of bins for STU/TRU ratios histograms
236 Float_t fMaxSTUTRURatio ; // Maximum value for STU/TRU ratios histograms
237 Int_t fNBinsClusterE ; // Number of bins for E cluster histograms
238 Float_t fMaxClusterE ; // Maximum value for E cluster histograms
240 //Constants needed by the class: EMCAL
241 static const int fgkFALTRORows = AliEMCALGeoParams::fgkEMCALRows*(AliEMCALGeoParams::fgkEMCALModules-7)/2; // total number
242 // of fake altro rows in EMCAL
243 // (ALTRO channels in one SM times 5 SM divided by 2 per FALTRO)
245 static const int fgkFALTROCols = AliEMCALGeoParams::fgkEMCALCols; // total number of fake altro columns in EMCAL
246 // (ALTRO channels in one SM times 2 SM divided by 2 per FALTRO)
249 AliAnalysisTaskEMCALTriggerQA (const AliAnalysisTaskEMCALTriggerQA&); // not implemented
251 AliAnalysisTaskEMCALTriggerQA& operator=(const AliAnalysisTaskEMCALTriggerQA&); // not implemented
253 ClassDef(AliAnalysisTaskEMCALTriggerQA, 10);