Skip to content

Commit e20c38a

Browse files
authored
PWGEM/PhotonMeson: Adapt event histogram in EMC QC task (#5107)
- EMC QC task was broken by changes to event histogram from PR 395efdb - This commit defines a new Event histogram to be filled by the EMC QC task ONLY
1 parent 788c7b7 commit e20c38a

1 file changed

Lines changed: 18 additions & 12 deletions

File tree

PWGEM/PhotonMeson/Tasks/emcalQC.cxx

Lines changed: 18 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -76,6 +76,14 @@ struct emcalQC {
7676
THashList* list_ev = reinterpret_cast<THashList*>(fMainList->FindObject("Event"));
7777
o2::aod::pwgem::photon::histogram::DefineHistograms(list_ev, "Event");
7878

79+
list_ev->Add(new TH1F("hNEvents", "hNEvents", 6, 0.5f, 6.5f));
80+
reinterpret_cast<TH2F*>(list_ev->FindObject("hNEvents"))->GetXaxis()->SetBinLabel(1, "all cols");
81+
reinterpret_cast<TH2F*>(list_ev->FindObject("hNEvents"))->GetXaxis()->SetBinLabel(2, "sel8");
82+
reinterpret_cast<TH2F*>(list_ev->FindObject("hNEvents"))->GetXaxis()->SetBinLabel(3, "emc readout");
83+
reinterpret_cast<TH2F*>(list_ev->FindObject("hNEvents"))->GetXaxis()->SetBinLabel(5, "1+ Contributor");
84+
reinterpret_cast<TH2F*>(list_ev->FindObject("hNEvents"))->GetXaxis()->SetBinLabel(6, "z<10cm");
85+
reinterpret_cast<TH2F*>(list_ev->FindObject("hNEvents"))->GetXaxis()->SetBinLabel(7, "unique col");
86+
7987
o2::aod::pwgem::photon::histogram::AddHistClass(fMainList, "Cluster");
8088
THashList* list_cluster = reinterpret_cast<THashList*>(fMainList->FindObject("Cluster"));
8189

@@ -139,12 +147,6 @@ struct emcalQC {
139147
DefineCuts();
140148
addhistograms(); // please call this after DefinCuts();
141149

142-
reinterpret_cast<TH2F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(1, "all cols");
143-
reinterpret_cast<TH2F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(2, "sel8 + readout");
144-
reinterpret_cast<TH2F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(3, "1+ Contributor");
145-
reinterpret_cast<TH2F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(4, "z<10cm");
146-
reinterpret_cast<TH2F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(5, "unique col");
147-
148150
fOutputEvent.setObject(reinterpret_cast<THashList*>(fMainList->FindObject("Event")));
149151
fOutputCluster.setObject(reinterpret_cast<THashList*>(fMainList->FindObject("Cluster")));
150152
}
@@ -157,27 +159,31 @@ struct emcalQC {
157159
THashList* list_cluster = static_cast<THashList*>(fMainList->FindObject("Cluster"));
158160

159161
for (auto& collision : collisions) {
160-
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(1.0);
161-
if (!collision.sel8() || collision.isEMCreadout() == 0) { // Check sel8 and whether EMC was read out
162+
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hNEvents"))->Fill(1.0);
163+
if (!collision.sel8()) {
164+
continue;
165+
}
166+
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hNEvents"))->Fill(2.0);
167+
if (collision.isEMCreadout() == 0) {
162168
continue;
163169
}
164-
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(2.0);
170+
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hNEvents"))->Fill(3.0);
165171

166172
if (collision.numContrib() < 0.5) {
167173
continue;
168174
}
169-
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(3.0);
175+
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hNEvents"))->Fill(4.0);
170176

171177
if (abs(collision.posZ()) > 10.0) {
172178
continue;
173179
}
174-
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(4.0);
180+
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hNEvents"))->Fill(5.0);
175181
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hZvtx"))->Fill(collision.posZ());
176182

177183
if (collision.ncollsPerBC() != 1) { // Check that the collision is unique (the only one in the bc)
178184
continue;
179185
}
180-
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(5.0);
186+
reinterpret_cast<TH1F*>(fMainList->FindObject("Event")->FindObject("hNEvents"))->Fill(6.0);
181187

182188
o2::aod::pwgem::photon::histogram::FillHistClass<EMHistType::kEvent>(list_ev, "", collision); // Fill event histograms
183189

0 commit comments

Comments
 (0)