Skip to content

Commit 66e41b2

Browse files
coppediscoppedis
committed
following Anton's suggestion
1 parent baccf0c commit 66e41b2

1 file changed

Lines changed: 40 additions & 38 deletions

File tree

‎PWGMM/UE/Tasks/ue-zdc-analysys.cxx‎

Lines changed: 40 additions & 38 deletions
Original file line numberDiff line numberDiff line change
@@ -196,46 +196,47 @@ struct ZDCAnalysis {
196196
for (auto& collision : cols) {
197197
const auto& foundBC = collision.foundBC_as<BCsRun3>();
198198
if (foundBC.has_zdc()) {
199+
const auto& zdcread = foundBC.zdc();
199200
if (TDCcut) {
200-
if ((foundBC.zdc().timeZNA() >= tdcZNmincut) && (foundBC.zdc().timeZNA() <= tdcZNmaxcut))
201-
registry.get<TH1>(HIST("ZNAcoll"))->Fill(foundBC.zdc().amplitudeZNA());
202-
if ((foundBC.zdc().timeZNC() >= tdcZNmincut) && (foundBC.zdc().timeZNC() <= tdcZNmaxcut))
203-
registry.get<TH1>(HIST("ZNCcoll"))->Fill(foundBC.zdc().amplitudeZNC());
204-
if ((foundBC.zdc().timeZPA() >= tdcZPmincut) && (foundBC.zdc().timeZPA() <= tdcZPmaxcut))
205-
registry.get<TH1>(HIST("ZPAcoll"))->Fill(foundBC.zdc().amplitudeZPA());
206-
if ((foundBC.zdc().timeZPC() >= tdcZPmincut) && (foundBC.zdc().timeZPC() <= tdcZPmaxcut))
207-
registry.get<TH1>(HIST("ZPCcoll"))->Fill(foundBC.zdc().amplitudeZPC());
208-
if (((foundBC.zdc().timeZNC() >= tdcZNmincut) && (foundBC.zdc().timeZNC() <= tdcZNmaxcut)) && ((foundBC.zdc().timeZNA() >= tdcZNmincut) && (foundBC.zdc().timeZNA() <= tdcZNmaxcut)))
209-
registry.get<TH2>(HIST("ZNvsZEMcoll"))->Fill(foundBC.zdc().amplitudeZEM1() + foundBC.zdc().amplitudeZEM2(), foundBC.zdc().amplitudeZNA() + foundBC.zdc().amplitudeZNC());
210-
if (((foundBC.zdc().timeZNC() >= tdcZNmincut) && (foundBC.zdc().timeZNC() <= tdcZNmaxcut)) && ((foundBC.zdc().timeZNA() >= tdcZNmincut) && (foundBC.zdc().timeZNA() <= tdcZNmaxcut)))
211-
registry.get<TH2>(HIST("ZNAvsZNCcoll"))->Fill(foundBC.zdc().amplitudeZNC(), foundBC.zdc().amplitudeZNA());
212-
if (((foundBC.zdc().timeZPC() >= tdcZPmincut) && (foundBC.zdc().timeZPC() <= tdcZPmaxcut)) && ((foundBC.zdc().timeZPA() >= tdcZPmincut) && (foundBC.zdc().timeZPA() <= tdcZPmaxcut)))
213-
registry.get<TH2>(HIST("ZPAvsZPCcoll"))->Fill(foundBC.zdc().amplitudeZPC(), foundBC.zdc().amplitudeZPA());
214-
if ((foundBC.zdc().timeZNA() >= tdcZNmincut) && (foundBC.zdc().timeZNA() <= tdcZNmaxcut))
215-
registry.get<TH2>(HIST("ZNAvsZPAcoll"))->Fill(foundBC.zdc().amplitudeZPA(), foundBC.zdc().amplitudeZNA());
216-
if ((foundBC.zdc().timeZNC() >= tdcZNmincut) && (foundBC.zdc().timeZNC() <= tdcZNmaxcut))
217-
registry.get<TH2>(HIST("ZNCvsZPCcoll"))->Fill(foundBC.zdc().amplitudeZPC(), foundBC.zdc().amplitudeZNC());
201+
if ((zdcread.timeZNA() >= tdcZNmincut) && (zdcread.timeZNA() <= tdcZNmaxcut))
202+
registry.get<TH1>(HIST("ZNAcoll"))->Fill(zdcread.amplitudeZNA());
203+
if ((zdcread.timeZNC() >= tdcZNmincut) && (zdcread.timeZNC() <= tdcZNmaxcut))
204+
registry.get<TH1>(HIST("ZNCcoll"))->Fill(zdcread.amplitudeZNC());
205+
if ((zdcread.timeZPA() >= tdcZPmincut) && (zdcread.timeZPA() <= tdcZPmaxcut))
206+
registry.get<TH1>(HIST("ZPAcoll"))->Fill(zdcread.amplitudeZPA());
207+
if ((zdcread.timeZPC() >= tdcZPmincut) && (zdcread.timeZPC() <= tdcZPmaxcut))
208+
registry.get<TH1>(HIST("ZPCcoll"))->Fill(zdcread.amplitudeZPC());
209+
if (((zdcread.timeZNC() >= tdcZNmincut) && (zdcread.timeZNC() <= tdcZNmaxcut)) && ((zdcread.timeZNA() >= tdcZNmincut) && (zdcread.timeZNA() <= tdcZNmaxcut)))
210+
registry.get<TH2>(HIST("ZNvsZEMcoll"))->Fill(zdcread.amplitudeZEM1() + zdcread.amplitudeZEM2(), zdcread.amplitudeZNA() + zdcread.amplitudeZNC());
211+
if (((zdcread.timeZNC() >= tdcZNmincut) && (zdcread.timeZNC() <= tdcZNmaxcut)) && ((zdcread.timeZNA() >= tdcZNmincut) && (zdcread.timeZNA() <= tdcZNmaxcut)))
212+
registry.get<TH2>(HIST("ZNAvsZNCcoll"))->Fill(zdcread.amplitudeZNC(), zdcread.amplitudeZNA());
213+
if (((zdcread.timeZPC() >= tdcZPmincut) && (zdcread.timeZPC() <= tdcZPmaxcut)) && ((zdcread.timeZPA() >= tdcZPmincut) && (zdcread.timeZPA() <= tdcZPmaxcut)))
214+
registry.get<TH2>(HIST("ZPAvsZPCcoll"))->Fill(zdcread.amplitudeZPC(), zdcread.amplitudeZPA());
215+
if ((zdcread.timeZNA() >= tdcZNmincut) && (zdcread.timeZNA() <= tdcZNmaxcut))
216+
registry.get<TH2>(HIST("ZNAvsZPAcoll"))->Fill(zdcread.amplitudeZPA(), zdcread.amplitudeZNA());
217+
if ((zdcread.timeZNC() >= tdcZNmincut) && (zdcread.timeZNC() <= tdcZNmaxcut))
218+
registry.get<TH2>(HIST("ZNCvsZPCcoll"))->Fill(zdcread.amplitudeZPC(), zdcread.amplitudeZNC());
218219
//
219220
} else {
220-
registry.get<TH1>(HIST("ZNAcoll"))->Fill(foundBC.zdc().amplitudeZNA());
221-
registry.get<TH1>(HIST("ZNCcoll"))->Fill(foundBC.zdc().amplitudeZNC());
222-
registry.get<TH1>(HIST("ZPAcoll"))->Fill(foundBC.zdc().amplitudeZPA());
223-
registry.get<TH1>(HIST("ZPCcoll"))->Fill(foundBC.zdc().amplitudeZPC());
224-
registry.get<TH2>(HIST("ZNvsZEMcoll"))->Fill(foundBC.zdc().amplitudeZEM1() + foundBC.zdc().amplitudeZEM2(), foundBC.zdc().amplitudeZNA() + foundBC.zdc().amplitudeZNC());
225-
registry.get<TH2>(HIST("ZNAvsZNCcoll"))->Fill(foundBC.zdc().amplitudeZNC(), foundBC.zdc().amplitudeZNA());
226-
registry.get<TH2>(HIST("ZPAvsZPCcoll"))->Fill(foundBC.zdc().amplitudeZPC(), foundBC.zdc().amplitudeZPA());
227-
registry.get<TH2>(HIST("ZNAvsZPAcoll"))->Fill(foundBC.zdc().amplitudeZPA(), foundBC.zdc().amplitudeZNA());
228-
registry.get<TH2>(HIST("ZNCvsZPCcoll"))->Fill(foundBC.zdc().amplitudeZPC(), foundBC.zdc().amplitudeZNC());
221+
registry.get<TH1>(HIST("ZNAcoll"))->Fill(zdcread.amplitudeZNA());
222+
registry.get<TH1>(HIST("ZNCcoll"))->Fill(zdcread.amplitudeZNC());
223+
registry.get<TH1>(HIST("ZPAcoll"))->Fill(zdcread.amplitudeZPA());
224+
registry.get<TH1>(HIST("ZPCcoll"))->Fill(zdcread.amplitudeZPC());
225+
registry.get<TH2>(HIST("ZNvsZEMcoll"))->Fill(zdcread.amplitudeZEM1() + zdcread.amplitudeZEM2(), zdcread.amplitudeZNA() + zdcread.amplitudeZNC());
226+
registry.get<TH2>(HIST("ZNAvsZNCcoll"))->Fill(zdcread.amplitudeZNC(), zdcread.amplitudeZNA());
227+
registry.get<TH2>(HIST("ZPAvsZPCcoll"))->Fill(zdcread.amplitudeZPC(), zdcread.amplitudeZPA());
228+
registry.get<TH2>(HIST("ZNAvsZPAcoll"))->Fill(zdcread.amplitudeZPA(), zdcread.amplitudeZNA());
229+
registry.get<TH2>(HIST("ZNCvsZPCcoll"))->Fill(zdcread.amplitudeZPC(), zdcread.amplitudeZNC());
229230
}
230-
registry.get<TH1>(HIST("ZEM1coll"))->Fill(foundBC.zdc().amplitudeZEM1());
231-
registry.get<TH1>(HIST("ZEM2coll"))->Fill(foundBC.zdc().amplitudeZEM2());
231+
registry.get<TH1>(HIST("ZEM1coll"))->Fill(zdcread.amplitudeZEM1());
232+
registry.get<TH1>(HIST("ZEM2coll"))->Fill(zdcread.amplitudeZEM2());
232233
//
233-
registry.get<TH2>(HIST("ZNCvstdccoll"))->Fill(foundBC.zdc().timeZNC(), foundBC.zdc().amplitudeZNC());
234-
registry.get<TH2>(HIST("ZNAvstdccoll"))->Fill(foundBC.zdc().timeZNA(), foundBC.zdc().amplitudeZNA());
235-
registry.get<TH2>(HIST("ZPCvstdccoll"))->Fill(foundBC.zdc().timeZPC(), foundBC.zdc().amplitudeZPC());
236-
registry.get<TH2>(HIST("ZPAvstdccoll"))->Fill(foundBC.zdc().timeZPA(), foundBC.zdc().amplitudeZPA());
237-
registry.get<TH2>(HIST("ZEM1vstdccoll"))->Fill(foundBC.zdc().timeZEM1(), foundBC.zdc().amplitudeZEM1());
238-
registry.get<TH2>(HIST("ZEM2vstdccoll"))->Fill(foundBC.zdc().timeZEM2(), foundBC.zdc().amplitudeZEM2());
234+
registry.get<TH2>(HIST("ZNCvstdccoll"))->Fill(zdcread.timeZNC(), zdcread.amplitudeZNC());
235+
registry.get<TH2>(HIST("ZNAvstdccoll"))->Fill(zdcread.timeZNA(), zdcread.amplitudeZNA());
236+
registry.get<TH2>(HIST("ZPCvstdccoll"))->Fill(zdcread.timeZPC(), zdcread.amplitudeZPC());
237+
registry.get<TH2>(HIST("ZPAvstdccoll"))->Fill(zdcread.timeZPA(), zdcread.amplitudeZPA());
238+
registry.get<TH2>(HIST("ZEM1vstdccoll"))->Fill(zdcread.timeZEM1(), zdcread.amplitudeZEM1());
239+
registry.get<TH2>(HIST("ZEM2vstdccoll"))->Fill(zdcread.timeZEM2(), zdcread.amplitudeZEM2());
239240
}
240241
}
241242
}
@@ -288,9 +289,10 @@ struct ZDCAnalysis {
288289
}
289290

290291
if (foundBC.has_zdc()) {
291-
registry.get<TH2>(HIST("ZNvsFV0Acorrel"))->Fill(multV0A / 100., foundBC.zdc().amplitudeZNA() + foundBC.zdc().amplitudeZNC());
292-
registry.get<TH2>(HIST("ZNvsFT0correl"))->Fill((multT0A + multT0C) / 100., foundBC.zdc().amplitudeZNC() + foundBC.zdc().amplitudeZNA());
293-
registry.get<TH2>(HIST("ZNvsFDDcorrel"))->Fill(multFDC + multFDA, foundBC.zdc().amplitudeZNC() + foundBC.zdc().amplitudeZNA());
292+
const auto& zdcread = foundBC.zdc();
293+
registry.get<TH2>(HIST("ZNvsFV0Acorrel"))->Fill(multV0A / 100., zdcread.amplitudeZNA() + zdcread.amplitudeZNC());
294+
registry.get<TH2>(HIST("ZNvsFT0correl"))->Fill((multT0A + multT0C) / 100., zdcread.amplitudeZNC() + zdcread.amplitudeZNA());
295+
registry.get<TH2>(HIST("ZNvsFDDcorrel"))->Fill(multFDC + multFDA, zdcread.amplitudeZNC() + zdcread.amplitudeZNA());
294296
}
295297
}
296298
PROCESS_SWITCH(ZDCAnalysis, processZdcCorrela, "Processing ZDC vs. mult. w. collision association", true);

0 commit comments

Comments
 (0)