Skip to content

Commit c049190

Browse files
authored
[PWGEM-10] createEMReducedEvent add bin to check for sel8 and MinNGamma (#3502)
- added an extra bin to the event histogram to check if all three min number of gamma candiadates and the sel8 condidations are met at the same time. Right now we can only see if either sel8 is selected or if the minimum nuber of photons for each of the methods are there. But there could be the case where we have no overlap between the events statisfying each of the conditions.
1 parent 835a938 commit c049190

1 file changed

Lines changed: 5 additions & 1 deletion

File tree

PWGEM/PhotonMeson/TableProducer/createEMReducedEvent.cxx

Lines changed: 5 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -55,13 +55,14 @@ struct createEMReducedEvent {
5555

5656
void init(o2::framework::InitContext&)
5757
{
58-
auto hEventCounter = registry.add<TH1>("hEventCounter", "hEventCounter", kTH1I, {{6, 0.5f, 6.5f}});
58+
auto hEventCounter = registry.add<TH1>("hEventCounter", "hEventCounter", kTH1I, {{7, 0.5f, 7.5f}});
5959
hEventCounter->GetXaxis()->SetBinLabel(1, "all");
6060
hEventCounter->GetXaxis()->SetBinLabel(2, "has > minN_PCM");
6161
hEventCounter->GetXaxis()->SetBinLabel(3, "has > minN_PHOS");
6262
hEventCounter->GetXaxis()->SetBinLabel(4, "has > minN_EMC");
6363
hEventCounter->GetXaxis()->SetBinLabel(5, "has > minN_any");
6464
hEventCounter->GetXaxis()->SetBinLabel(6, "sel8");
65+
hEventCounter->GetXaxis()->SetBinLabel(7, "sel8 & minN");
6566
registry.add<TH1>("hNGammas_PCM", ";#it{N}_{#gamma,PCM};#it{count}", kTH1I, {{21, -0.5, 20.5}});
6667
registry.add<TH1>("hNGammas_PHOS", ";#it{N}_{#gamma,PHOS};#it{count}", kTH1I, {{21, -0.5, 20.5}});
6768
registry.add<TH1>("hNGammas_EMC", ";#it{N}_{#gamma,EMC};#it{count}", kTH1I, {{21, -0.5, 20.5}});
@@ -126,6 +127,9 @@ struct createEMReducedEvent {
126127
if (collision.sel8()) {
127128
registry.fill(HIST("hEventCounter"), 6);
128129
}
130+
if (collision.sel8() && ng_emc >= minN_EMC && ng_phos >= minN_PHOS && ng_pcm >= minN_PCM) {
131+
registry.fill(HIST("hEventCounter"), 7);
132+
}
129133

130134
// store event selection decisions
131135
uint64_t tag = collision.selection_raw();

0 commit comments

Comments
 (0)