Skip to content

Commit bcf9975

Browse files
dsekihatdsekihat
authored
PWGEM/Dilepton: add a task to select EoI (AliceO2Group#7927)
1 parent 5d2f00d commit bcf9975

16 files changed

+435
-103
lines changed

‎PWGEM/Dilepton/Core/Dilepton.h‎

Lines changed: 51 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -101,7 +101,7 @@ struct Dilepton {
101101
Configurable<int> cfgAnalysisType{"cfgAnalysisType", static_cast<int>(o2::aod::pwgem::dilepton::utils::pairutil::DileptonAnalysisType::kQC), "kQC:0, kUPC:1, kFlowV2:2, kFlowV3:3, kPolarization:4, kVM:5, kHFll:6"};
102102
Configurable<int> cfgEP2Estimator_for_Mix{"cfgEP2Estimator_for_Mix", 3, "FT0M:0, FT0A:1, FT0C:2, BTot:3, BPos:4, BNeg:5"};
103103
Configurable<int> cfgQvecEstimator{"cfgQvecEstimator", 0, "FT0M:0, FT0A:1, FT0C:2, BTot:3, BPos:4, BNeg:5"};
104-
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2, NTPV:3"};
104+
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2"};
105105
Configurable<float> cfgCentMin{"cfgCentMin", 0, "min. centrality"};
106106
Configurable<float> cfgCentMax{"cfgCentMax", 999.f, "max. centrality"};
107107
Configurable<bool> cfgDoMix{"cfgDoMix", true, "flag for event mixing"};
@@ -170,8 +170,8 @@ struct Dilepton {
170170
Configurable<int> cfg_min_ncrossedrows{"cfg_min_ncrossedrows", 100, "min ncrossed rows"};
171171
Configurable<float> cfg_max_chi2tpc{"cfg_max_chi2tpc", 4.0, "max chi2/NclsTPC"};
172172
Configurable<float> cfg_max_chi2its{"cfg_max_chi2its", 5.0, "max chi2/NclsITS"};
173-
Configurable<float> cfg_max_dcaxy{"cfg_max_dcaxy", 1.0, "max dca XY for single track in cm"};
174-
Configurable<float> cfg_max_dcaz{"cfg_max_dcaz", 1.0, "max dca Z for single track in cm"};
173+
Configurable<float> cfg_max_dcaxy{"cfg_max_dcaxy", 0.2, "max dca XY for single track in cm"};
174+
Configurable<float> cfg_max_dcaz{"cfg_max_dcaz", 0.2, "max dca Z for single track in cm"};
175175
Configurable<bool> cfg_require_itsib_any{"cfg_require_itsib_any", false, "flag to require ITS ib any hits"};
176176
Configurable<bool> cfg_require_itsib_1st{"cfg_require_itsib_1st", true, "flag to require ITS ib 1st hit"};
177177

@@ -369,6 +369,7 @@ struct Dilepton {
369369
// fwdfitter.setTGeoMat(false);
370370
}
371371

372+
fRegistry.addClone("Event/before/hCollisionCounter", "Event/norm/hCollisionCounter");
372373
fRegistry.add("Pair/mix/hDiffBC", "diff. global BC in mixed event;|BC_{current} - BC_{mixed}|", kTH1D, {{1001, -0.5, 1000.5}}, true);
373374
if (doprocessTriggerAnalysis) {
374375
fRegistry.add("Event/hNInspectedTVX", "N inspected TVX;run number;N_{TVX}", kTProfile, {{80000, 520000.5, 600000.5}}, true);
@@ -1052,7 +1053,7 @@ struct Dilepton {
10521053
{
10531054
for (auto& collision : collisions) {
10541055
initCCDB<isTriggerAnalysis>(collision);
1055-
const float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
1056+
const float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
10561057
float centrality = centralities[cfgCentEstimator];
10571058
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
10581059
continue;
@@ -1409,7 +1410,7 @@ struct Dilepton {
14091410

14101411
for (auto& collision : collisions) {
14111412
initCCDB<isTriggerAnalysis>(collision);
1412-
const float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
1413+
const float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
14131414
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
14141415
continue;
14151416
}
@@ -1543,6 +1544,51 @@ struct Dilepton {
15431544
}
15441545
PROCESS_SWITCH(Dilepton, processTriggerAnalysis, "run dilepton analysis on triggered data", false);
15451546

1547+
void processNorm(aod::EMEventNormInfos const& collisions)
1548+
{
1549+
for (auto& collision : collisions) {
1550+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 1.0);
1551+
if (collision.selection_bit(o2::aod::evsel::kIsTriggerTVX)) {
1552+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 2.0);
1553+
}
1554+
if (collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
1555+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 3.0);
1556+
}
1557+
if (collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
1558+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 4.0);
1559+
}
1560+
if (collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
1561+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 5.0);
1562+
}
1563+
if (collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
1564+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 6.0);
1565+
}
1566+
if (collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC)) {
1567+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 7.0);
1568+
}
1569+
if (collision.selection_bit(o2::aod::evsel::kIsVertexTRDmatched)) {
1570+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 8.0);
1571+
}
1572+
if (collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched)) {
1573+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 9.0);
1574+
}
1575+
if (collision.sel8()) {
1576+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 10.0);
1577+
}
1578+
if (abs(collision.posZ()) < 10.0) {
1579+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 11.0);
1580+
}
1581+
if (collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
1582+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 12.0);
1583+
}
1584+
if (!fEMEventCut.IsSelected(collision)) {
1585+
continue;
1586+
}
1587+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), o2::aod::pwgem::dilepton::utils::eventhistogram::nbin_ev); // accepted
1588+
} // end of collision loop
1589+
}
1590+
PROCESS_SWITCH(Dilepton, processNorm, "process normalization info", true);
1591+
15461592
void processDummy(MyCollisions const&) {}
15471593
PROCESS_SWITCH(Dilepton, processDummy, "Dummy function", false);
15481594
};

‎PWGEM/Dilepton/Core/DileptonMC.h‎

Lines changed: 53 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -91,7 +91,7 @@ struct DileptonMC {
9191
Configurable<float> d_bz_input{"d_bz_input", -999, "bz field in kG, -999 is automatic"};
9292

9393
Configurable<int> cfgEventGeneratorType{"cfgEventGeneratorType", -1, "if positive, select event generator type. i.e. gap or signal"};
94-
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2, NTPV:3"};
94+
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2"};
9595
Configurable<float> cfgCentMin{"cfgCentMin", 0, "min. centrality"};
9696
Configurable<float> cfgCentMax{"cfgCentMax", 999.f, "max. centrality"};
9797
Configurable<int> cfgNtracksPV08Min{"cfgNtracksPV08Min", -1, "min. multNTracksPV"};
@@ -151,8 +151,8 @@ struct DileptonMC {
151151
Configurable<int> cfg_min_ncrossedrows{"cfg_min_ncrossedrows", 100, "min ncrossed rows"};
152152
Configurable<float> cfg_max_chi2tpc{"cfg_max_chi2tpc", 4.0, "max chi2/NclsTPC"};
153153
Configurable<float> cfg_max_chi2its{"cfg_max_chi2its", 5.0, "max chi2/NclsITS"};
154-
Configurable<float> cfg_max_dcaxy{"cfg_max_dcaxy", 1.0, "max dca XY for single track in cm"};
155-
Configurable<float> cfg_max_dcaz{"cfg_max_dcaz", 1.0, "max dca Z for single track in cm"};
154+
Configurable<float> cfg_max_dcaxy{"cfg_max_dcaxy", 0.2, "max dca XY for single track in cm"};
155+
Configurable<float> cfg_max_dcaz{"cfg_max_dcaz", 0.2, "max dca Z for single track in cm"};
156156
Configurable<bool> cfg_require_itsib_any{"cfg_require_itsib_any", false, "flag to require ITS ib any hits"};
157157
Configurable<bool> cfg_require_itsib_1st{"cfg_require_itsib_1st", true, "flag to require ITS ib 1st hit"};
158158

@@ -401,6 +401,7 @@ struct DileptonMC {
401401
// fwdfitter.setUseAbsDCA(true);
402402
// fwdfitter.setTGeoMat(false);
403403
}
404+
fRegistry.addClone("Event/before/hCollisionCounter", "Event/norm/hCollisionCounter");
404405
}
405406

406407
template <typename TCollision>
@@ -911,7 +912,7 @@ struct DileptonMC {
911912
{
912913
for (auto& collision : collisions) {
913914
initCCDB(collision);
914-
float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
915+
float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
915916
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
916917
continue;
917918
}
@@ -983,7 +984,7 @@ struct DileptonMC {
983984
// all MC tracks which belong to the MC event corresponding to the current reconstructed event
984985

985986
for (auto& collision : collisions) {
986-
float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
987+
float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
987988
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
988989
continue;
989990
}
@@ -1521,7 +1522,7 @@ struct DileptonMC {
15211522

15221523
for (auto& collision : collisions) {
15231524
initCCDB(collision);
1524-
const float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
1525+
const float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
15251526
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
15261527
continue;
15271528
}
@@ -1682,7 +1683,7 @@ struct DileptonMC {
16821683
{
16831684
// for oemga, phi efficiency
16841685
for (auto& collision : collisions) {
1685-
float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
1686+
float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
16861687
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
16871688
continue;
16881689
}
@@ -1728,6 +1729,51 @@ struct DileptonMC {
17281729
}
17291730
PROCESS_SWITCH(DileptonMC, processGen_VM, "process generated info for vector mesons", false);
17301731

1732+
void processNorm(aod::EMEventNormInfos const& collisions)
1733+
{
1734+
for (auto& collision : collisions) {
1735+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 1.0);
1736+
if (collision.selection_bit(o2::aod::evsel::kIsTriggerTVX)) {
1737+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 2.0);
1738+
}
1739+
if (collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
1740+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 3.0);
1741+
}
1742+
if (collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
1743+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 4.0);
1744+
}
1745+
if (collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
1746+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 5.0);
1747+
}
1748+
if (collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
1749+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 6.0);
1750+
}
1751+
if (collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC)) {
1752+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 7.0);
1753+
}
1754+
if (collision.selection_bit(o2::aod::evsel::kIsVertexTRDmatched)) {
1755+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 8.0);
1756+
}
1757+
if (collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched)) {
1758+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 9.0);
1759+
}
1760+
if (collision.sel8()) {
1761+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 10.0);
1762+
}
1763+
if (abs(collision.posZ()) < 10.0) {
1764+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 11.0);
1765+
}
1766+
if (collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
1767+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 12.0);
1768+
}
1769+
if (!fEMEventCut.IsSelected(collision)) {
1770+
continue;
1771+
}
1772+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), o2::aod::pwgem::dilepton::utils::eventhistogram::nbin_ev); // accepted
1773+
} // end of collision loop
1774+
}
1775+
PROCESS_SWITCH(DileptonMC, processNorm, "process normalization info", true);
1776+
17311777
void processDummy(MyCollisions const&) {}
17321778
PROCESS_SWITCH(DileptonMC, processDummy, "Dummy function", false);
17331779
};

‎PWGEM/Dilepton/Core/PhotonHBT.h‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -96,7 +96,7 @@ struct PhotonHBT {
9696

9797
Configurable<bool> cfgDo3D{"cfgDo3D", false, "enable 3D analysis"};
9898
Configurable<int> cfgEP2Estimator_for_Mix{"cfgEP2Estimator_for_Mix", 3, "FT0M:0, FT0A:1, FT0C:2, BTot:3, BPos:4, BNeg:5"};
99-
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2, NTPV:3"};
99+
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2"};
100100
Configurable<float> cfgCentMin{"cfgCentMin", 0, "min. centrality"};
101101
Configurable<float> cfgCentMax{"cfgCentMax", 999, "max. centrality"};
102102
// Configurable<float> cfgSpherocityMin{"cfgSpherocityMin", -999.f, "min. spherocity"};
@@ -643,7 +643,7 @@ struct PhotonHBT {
643643
for (auto& collision : collisions) {
644644
initCCDB<isTriggerAnalysis>(collision);
645645
int ndiphoton = 0;
646-
const float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
646+
const float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
647647
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
648648
continue;
649649
}
@@ -1228,7 +1228,7 @@ struct PhotonHBT {
12281228

12291229
for (auto& collision : collisions) {
12301230
initCCDB<isTriggerAnalysis>(collision);
1231-
const float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
1231+
const float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
12321232
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
12331233
continue;
12341234
}

‎PWGEM/Dilepton/Core/SingleTrackQC.h‎

Lines changed: 52 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -72,7 +72,7 @@ struct SingleTrackQC {
7272
// Configurables
7373
Configurable<std::string> ccdburl{"ccdb-url", "http://alice-ccdb.cern.ch", "url of the ccdb repository"};
7474

75-
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2, NTPV:3"};
75+
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2"};
7676
Configurable<float> cfgCentMin{"cfgCentMin", 0, "min. centrality"};
7777
Configurable<float> cfgCentMax{"cfgCentMax", 999.f, "max. centrality"};
7878
Configurable<std::string> cfg_swt_name{"cfg_swt_name", "fHighTrackMult", "desired software trigger name"}; // 1 trigger name per 1 task. fHighTrackMult, fHighFt0Mult
@@ -208,6 +208,7 @@ struct SingleTrackQC {
208208
fRegistry.add("Track/positive/hChi2TPC", "chi2/number of TPC clusters", kTH1F, {{100, 0, 10}}, false);
209209
fRegistry.add("Track/positive/hTPCNcr2Nf", "TPC Ncr/Nfindable", kTH1F, {{200, 0, 2}}, false);
210210
fRegistry.add("Track/positive/hTPCNcls2Nf", "TPC Ncls/Nfindable", kTH1F, {{200, 0, 2}}, false);
211+
fRegistry.add("Track/positive/hTPCNclsShared", "TPC Ncls shared/Ncls;p_{T} (GeV/c);N_{cls}^{shared}/N_{cls} in TPC", kTH2F, {{1000, 0, 10}, {100, 0, 1}}, false);
211212
fRegistry.add("Track/positive/hNclsITS", "number of ITS clusters", kTH1F, {{8, -0.5, 7.5}}, false);
212213
fRegistry.add("Track/positive/hChi2ITS", "chi2/number of ITS clusters", kTH1F, {{100, 0, 10}}, false);
213214
fRegistry.add("Track/positive/hITSClusterMap", "ITS cluster map", kTH1F, {{128, -0.5, 127.5}}, false);
@@ -264,6 +265,7 @@ struct SingleTrackQC {
264265
addhistograms();
265266
mRunNumber = 0;
266267

268+
fRegistry.addClone("Event/before/hCollisionCounter", "Event/norm/hCollisionCounter");
267269
if (doprocessQC_TriggeredData) {
268270
fRegistry.add("Event/hNInspectedTVX", "N inspected TVX;run number;N_{TVX}", kTProfile, {{80000, 520000.5, 600000.5}}, true);
269271
}
@@ -398,6 +400,7 @@ struct SingleTrackQC {
398400
fRegistry.fill(HIST("Track/positive/hNcrTPC"), track.tpcNClsCrossedRows());
399401
fRegistry.fill(HIST("Track/positive/hTPCNcr2Nf"), track.tpcCrossedRowsOverFindableCls());
400402
fRegistry.fill(HIST("Track/positive/hTPCNcls2Nf"), track.tpcFoundOverFindableCls());
403+
fRegistry.fill(HIST("Track/positive/hTPCNclsShared"), track.pt(), track.tpcFractionSharedCls());
401404
fRegistry.fill(HIST("Track/positive/hChi2TPC"), track.tpcChi2NCl());
402405
fRegistry.fill(HIST("Track/positive/hChi2ITS"), track.itsChi2NCl());
403406
fRegistry.fill(HIST("Track/positive/hITSClusterMap"), track.itsClusterMap());
@@ -427,6 +430,7 @@ struct SingleTrackQC {
427430
fRegistry.fill(HIST("Track/negative/hNcrTPC"), track.tpcNClsCrossedRows());
428431
fRegistry.fill(HIST("Track/negative/hTPCNcr2Nf"), track.tpcCrossedRowsOverFindableCls());
429432
fRegistry.fill(HIST("Track/negative/hTPCNcls2Nf"), track.tpcFoundOverFindableCls());
433+
fRegistry.fill(HIST("Track/negative/hTPCNclsShared"), track.pt(), track.tpcFractionSharedCls());
430434
fRegistry.fill(HIST("Track/negative/hChi2TPC"), track.tpcChi2NCl());
431435
fRegistry.fill(HIST("Track/negative/hChi2ITS"), track.itsChi2NCl());
432436
fRegistry.fill(HIST("Track/negative/hITSClusterMap"), track.itsClusterMap());
@@ -493,7 +497,7 @@ struct SingleTrackQC {
493497
{
494498
for (auto& collision : collisions) {
495499
initCCDB<isTriggerAnalysis>(collision);
496-
float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
500+
float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
497501
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
498502
continue;
499503
}
@@ -545,7 +549,7 @@ struct SingleTrackQC {
545549
passed_trackIds.reserve(tracks.size());
546550
for (auto& collision : collisions) {
547551
initCCDB<isTriggerAnalysis>(collision);
548-
float centralities[4] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C(), collision.centNTPV()};
552+
float centralities[3] = {collision.centFT0M(), collision.centFT0A(), collision.centFT0C()};
549553
if (centralities[cfgCentEstimator] < cfgCentMin || cfgCentMax < centralities[cfgCentEstimator]) {
550554
continue;
551555
}
@@ -670,6 +674,51 @@ struct SingleTrackQC {
670674
}
671675
PROCESS_SWITCH(SingleTrackQC, processQC_TriggeredData, "run single track QC on triggered data", false);
672676

677+
void processNorm(aod::EMEventNormInfos const& collisions)
678+
{
679+
for (auto& collision : collisions) {
680+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 1.0);
681+
if (collision.selection_bit(o2::aod::evsel::kIsTriggerTVX)) {
682+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 2.0);
683+
}
684+
if (collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
685+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 3.0);
686+
}
687+
if (collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
688+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 4.0);
689+
}
690+
if (collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
691+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 5.0);
692+
}
693+
if (collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
694+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 6.0);
695+
}
696+
if (collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC)) {
697+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 7.0);
698+
}
699+
if (collision.selection_bit(o2::aod::evsel::kIsVertexTRDmatched)) {
700+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 8.0);
701+
}
702+
if (collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched)) {
703+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 9.0);
704+
}
705+
if (collision.sel8()) {
706+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 10.0);
707+
}
708+
if (abs(collision.posZ()) < 10.0) {
709+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 11.0);
710+
}
711+
if (collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
712+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), 12.0);
713+
}
714+
if (!fEMEventCut.IsSelected(collision)) {
715+
continue;
716+
}
717+
fRegistry.fill(HIST("Event/norm/hCollisionCounter"), o2::aod::pwgem::dilepton::utils::eventhistogram::nbin_ev); // accepted
718+
} // end of collision loop
719+
}
720+
PROCESS_SWITCH(SingleTrackQC, processNorm, "process normalization info", true);
721+
673722
void processDummy(MyCollisions const&) {}
674723
PROCESS_SWITCH(SingleTrackQC, processDummy, "Dummy function", false);
675724
};

0 commit comments

Comments
 (0)