|
13 | 13 | /// \brief Task for ZDC |
14 | 14 | /// \author chiara.oppedisano@cern.ch |
15 | 15 |
|
16 | | -#include "Framework/AnalysisDataModel.h" |
17 | | -#include "Framework/AnalysisTask.h" |
18 | | -#include "Framework/HistogramRegistry.h" |
19 | | -#include "Framework/runDataProcessing.h" |
| 16 | +#include "PWGMM/DataModel/ZDCdmOxygen.h" |
20 | 17 |
|
21 | 18 | #include "Common/CCDB/EventSelectionParams.h" |
22 | 19 | #include "Common/CCDB/TriggerAliases.h" |
23 | | -#include "Common/DataModel/Centrality.h" |
24 | | -#include "Common/DataModel/Multiplicity.h" |
25 | | -#include "Common/DataModel/EventSelection.h" |
26 | 20 | #include "Common/Core/TrackSelection.h" |
27 | 21 | #include "Common/Core/trackUtilities.h" |
| 22 | +#include "Common/DataModel/Centrality.h" |
| 23 | +#include "Common/DataModel/EventSelection.h" |
| 24 | +#include "Common/DataModel/Multiplicity.h" |
28 | 25 |
|
29 | | -#include "PWGMM/DataModel/ZDCdmOxygen.h" |
| 26 | +#include "Framework/AnalysisDataModel.h" |
| 27 | +#include "Framework/AnalysisTask.h" |
| 28 | +#include "Framework/HistogramRegistry.h" |
| 29 | +#include "Framework/runDataProcessing.h" |
30 | 30 |
|
31 | 31 | #include "TH1F.h" |
32 | 32 | #include "TH2F.h" |
@@ -80,7 +80,7 @@ struct ZdcTaskOxygen { |
80 | 80 | { |
81 | 81 | registry.add("debunchHist", "ZN sum vs. diff; ZNA-ZNC (ns); ZNA+ZNC (ns)", {HistType::kTH2D, {{nBinsTiming, -20., 20.}, {nBinsTiming, -20., 20.}}}); |
82 | 82 |
|
83 | | - if(doprocessALICEcoll) { |
| 83 | + if (doprocessALICEcoll) { |
84 | 84 | registry.add("hEventCount", "Number of Event; Cut; #Events Passed Cut", {HistType::kTH1D, {{nEventSelections, 0, nEventSelections}}}); |
85 | 85 | registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(evSel_allEvents + 1, "All events"); |
86 | 86 | registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(evSel_zvtx + 1, "vtxZ"); |
@@ -156,63 +156,63 @@ struct ZdcTaskOxygen { |
156 | 156 |
|
157 | 157 | void processZDCautoTrig(aod::Zdc const& zdc) |
158 | 158 | { |
159 | | - // auto-triggered events for ZDC |
160 | | - |
161 | | - auto tdcZNA = zdc.timeZNA(); |
162 | | - auto tdcZNC = zdc.timeZNC(); |
163 | | - auto tdcZPA = zdc.timeZPA(); |
164 | | - auto tdcZPC = zdc.timeZPC(); |
165 | | - auto tdcZEM1 = zdc.timeZEM1(); |
166 | | - auto tdcZEM2 = zdc.timeZEM2(); |
167 | | - // |
168 | | - double zna = zdc.amplitudeZNA(); |
169 | | - double znc = zdc.amplitudeZNC(); |
170 | | - double zpa = zdc.amplitudeZPA(); |
171 | | - double zpc = zdc.amplitudeZPC(); |
172 | | - double zem1 = zdc.amplitudeZEM1(); |
173 | | - double zem2 = zdc.amplitudeZEM2(); |
174 | | - // |
175 | | - double pmcZNA = zdc.energyCommonZNA(); |
176 | | - double pmcZNC = zdc.energyCommonZNC(); |
177 | | - double pmqZNC[4] = { |
178 | | - 0, |
179 | | - 0, |
180 | | - 0, |
181 | | - 0, |
182 | | - }; |
183 | | - double pmqZNA[4] = { |
184 | | - 0, |
185 | | - 0, |
186 | | - 0, |
187 | | - 0, |
188 | | - }; |
189 | | - const int noofZNsectors = 4; |
190 | | - for (int itow = 0; itow < noofZNsectors; itow++) { |
191 | | - pmqZNA[itow] = (zdc.energySectorZNA())[itow]; |
192 | | - pmqZNC[itow] = (zdc.energySectorZNC())[itow]; |
193 | | - } |
| 159 | + // auto-triggered events for ZDC |
| 160 | + |
| 161 | + auto tdcZNA = zdc.timeZNA(); |
| 162 | + auto tdcZNC = zdc.timeZNC(); |
| 163 | + auto tdcZPA = zdc.timeZPA(); |
| 164 | + auto tdcZPC = zdc.timeZPC(); |
| 165 | + auto tdcZEM1 = zdc.timeZEM1(); |
| 166 | + auto tdcZEM2 = zdc.timeZEM2(); |
| 167 | + // |
| 168 | + double zna = zdc.amplitudeZNA(); |
| 169 | + double znc = zdc.amplitudeZNC(); |
| 170 | + double zpa = zdc.amplitudeZPA(); |
| 171 | + double zpc = zdc.amplitudeZPC(); |
| 172 | + double zem1 = zdc.amplitudeZEM1(); |
| 173 | + double zem2 = zdc.amplitudeZEM2(); |
| 174 | + // |
| 175 | + double pmcZNA = zdc.energyCommonZNA(); |
| 176 | + double pmcZNC = zdc.energyCommonZNC(); |
| 177 | + double pmqZNC[4] = { |
| 178 | + 0, |
| 179 | + 0, |
| 180 | + 0, |
| 181 | + 0, |
| 182 | + }; |
| 183 | + double pmqZNA[4] = { |
| 184 | + 0, |
| 185 | + 0, |
| 186 | + 0, |
| 187 | + 0, |
| 188 | + }; |
| 189 | + const int noofZNsectors = 4; |
| 190 | + for (int itow = 0; itow < noofZNsectors; itow++) { |
| 191 | + pmqZNA[itow] = (zdc.energySectorZNA())[itow]; |
| 192 | + pmqZNC[itow] = (zdc.energySectorZNC())[itow]; |
| 193 | + } |
194 | 194 |
|
195 | | - bool isZNChit = false, isZNAhit = false; |
196 | | - if (tdcCut) { // a narrow TDC window is set |
197 | | - if ((tdcZNC >= tdcZNmincut) && (tdcZNC <= tdcZNmaxcut)) { |
198 | | - isZNChit = true; |
199 | | - } |
200 | | - if ((tdcZNA >= tdcZNmincut) && (tdcZNA <= tdcZNmaxcut)) { |
201 | | - isZNAhit = true; |
202 | | - } |
203 | | - } else { // if no window on TDC is set |
204 | | - if (pmcZNC > 0.) { |
205 | | - isZNChit = true; |
206 | | - } |
207 | | - if (pmcZNA > 0.) { |
208 | | - isZNAhit = true; |
209 | | - } |
| 195 | + bool isZNChit = false, isZNAhit = false; |
| 196 | + if (tdcCut) { // a narrow TDC window is set |
| 197 | + if ((tdcZNC >= tdcZNmincut) && (tdcZNC <= tdcZNmaxcut)) { |
| 198 | + isZNChit = true; |
| 199 | + } |
| 200 | + if ((tdcZNA >= tdcZNmincut) && (tdcZNA <= tdcZNmaxcut)) { |
| 201 | + isZNAhit = true; |
| 202 | + } |
| 203 | + } else { // if no window on TDC is set |
| 204 | + if (pmcZNC > 0.) { |
| 205 | + isZNChit = true; |
210 | 206 | } |
211 | | - if (isZNChit && isZNAhit) { |
212 | | - registry.get<TH1>(HIST("debunchHist"))->Fill(zna-znc, zna+znc); |
| 207 | + if (pmcZNA > 0.) { |
| 208 | + isZNAhit = true; |
213 | 209 | } |
| 210 | + } |
| 211 | + if (isZNChit && isZNAhit) { |
| 212 | + registry.get<TH1>(HIST("debunchHist"))->Fill(zna - znc, zna + znc); |
| 213 | + } |
214 | 214 |
|
215 | | - zdcTableoo(tdcZNA, zna, pmcZNA, pmqZNA[0], pmqZNA[1], pmqZNA[2], pmqZNA[3], |
| 215 | + zdcTableoo(tdcZNA, zna, pmcZNA, pmqZNA[0], pmqZNA[1], pmqZNA[2], pmqZNA[3], |
216 | 216 | tdcZNC, znc, pmcZNC, pmqZNC[0], pmqZNC[1], pmqZNC[2], pmqZNC[3], |
217 | 217 | tdcZPA, zpa, tdcZPC, zpc, tdcZEM1, zem1, tdcZEM2, zem2, |
218 | 218 | -1, -1, -1, |
@@ -253,7 +253,7 @@ struct ZdcTaskOxygen { |
253 | 253 | for (auto const& amplitude : foundBC.fv0a().amplitude()) { |
254 | 254 | multV0A += amplitude; |
255 | 255 | } |
256 | | - } |
| 256 | + } |
257 | 257 |
|
258 | 258 | if (foundBC.has_zdc()) { |
259 | 259 | const auto& zdc = foundBC.zdc(); |
@@ -309,16 +309,16 @@ struct ZdcTaskOxygen { |
309 | 309 | } |
310 | 310 | } |
311 | 311 | if (isZNChit && isZNAhit) { |
312 | | - registry.get<TH1>(HIST("debunchHist"))->Fill(zna-znc, zna+znc); |
| 312 | + registry.get<TH1>(HIST("debunchHist"))->Fill(zna - znc, zna + znc); |
313 | 313 | } |
314 | 314 |
|
315 | 315 | zdcTableoo(tdcZNA, zna, pmcZNA, pmqZNA[0], pmqZNA[1], pmqZNA[2], pmqZNA[3], |
316 | | - tdcZNC, znc, pmcZNC, pmqZNC[0], pmqZNC[1], pmqZNC[2], pmqZNC[3], |
317 | | - tdcZPA, zpa, tdcZPC, zpc, tdcZEM1, zem1, tdcZEM2, zem2, |
318 | | - multFT0A, multFT0C, multV0A, |
319 | | - zv, |
320 | | - centralityFT0C, centralityFT0A, centralityFT0M, |
321 | | - evSelection); |
| 316 | + tdcZNC, znc, pmcZNC, pmqZNC[0], pmqZNC[1], pmqZNC[2], pmqZNC[3], |
| 317 | + tdcZPA, zpa, tdcZPC, zpc, tdcZEM1, zem1, tdcZEM2, zem2, |
| 318 | + multFT0A, multFT0C, multV0A, |
| 319 | + zv, |
| 320 | + centralityFT0C, centralityFT0A, centralityFT0M, |
| 321 | + evSelection); |
322 | 322 | } |
323 | 323 | } |
324 | 324 | } |
|
0 commit comments