Skip to content

Commit fc6bcf0

Browse files
authored
Refresh the event selection module
1 parent 9c08acb commit fc6bcf0

File tree

1 file changed

+26
-34
lines changed

1 file changed

+26
-34
lines changed

PWGLF/Tasks/GlobalEventProperties/heavyionMultiplicity.cxx

Lines changed: 26 additions & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -127,7 +127,7 @@ static constexpr TrackSelectionFlags::flagtype TrackSelectionDca =
127127
static constexpr TrackSelectionFlags::flagtype TrackSelectionDcaxyOnly =
128128
TrackSelectionFlags::kDCAxy;
129129

130-
AxisSpec axisEvent{10, 0.5, 10.5, "#Event", "EventAxis"};
130+
AxisSpec axisEvent{15, 0.5, 15.5, "#Event", "EventAxis"};
131131
AxisSpec axisVtxZ{40, -20, 20, "Vertex Z", "VzAxis"};
132132
AxisSpec axisEta{40, -2, 2, "#eta", "EtaAxis"};
133133
AxisSpec axisEtaExtended{100, -5, 5, "#eta", "EtaAxisExtended"};
@@ -189,15 +189,14 @@ struct HeavyionMultiplicity {
189189
ConfigurableAxis binsMult{"binsMult", {500, 0.0f, +500.0f}, ""};
190190
ConfigurableAxis binsDCA{"binsDCA", {500, -10.0f, 10.0f}, ""};
191191

192+
Configurable<bool> isApplyTFcut{"isApplyTFcut", true, "Enable TimeFrameBorder cut"};
193+
Configurable<bool> isApplyITSROcut{"isApplyITSROcut", true, "Enable ITS ReadOutFrameBorder cut"};
192194
Configurable<bool> isApplySameBunchPileup{"isApplySameBunchPileup", true, "Enable SameBunchPileup cut"};
193195
Configurable<bool> isApplyGoodZvtxFT0vsPV{"isApplyGoodZvtxFT0vsPV", true, "Enable GoodZvtxFT0vsPV cut"};
194-
Configurable<bool> isApplyExtraCorrCut{"isApplyExtraCorrCut", false, "Enable extra NPVtracks vs FTOC correlation cut"};
195-
Configurable<bool> isApplyExtraPhiCut{"isApplyExtraPhiCut", false, "Enable extra phi cut"};
196-
Configurable<float> npvTracksCut{"npvTracksCut", 1.0f, "Apply extra NPVtracks cut"};
197-
Configurable<float> ft0cCut{"ft0cCut", 1.0f, "Apply extra FT0C cut"};
198196
Configurable<bool> isApplyNoCollInTimeRangeStandard{"isApplyNoCollInTimeRangeStandard", true, "Enable NoCollInTimeRangeStandard cut"};
199197
Configurable<bool> isApplyNoCollInRofStandard{"isApplyNoCollInRofStandard", false, "Enable NoCollInRofStandard cut"};
200198
Configurable<bool> isApplyNoHighMultCollInPrevRof{"isApplyNoHighMultCollInPrevRof", false, "Enable NoHighMultCollInPrevRof cut"};
199+
Configurable<bool> isApplyInelgt0{"isApplyInelgt0", false, "Enable INEL > 0 condition"};
201200
Configurable<bool> isApplyFT0CbasedOccupancy{"isApplyFT0CbasedOccupancy", false, "Enable FT0CbasedOccupancy cut"};
202201
Configurable<bool> isApplyCentFT0C{"isApplyCentFT0C", true, "Centrality based on FT0C"};
203202
Configurable<bool> isApplyCentFV0A{"isApplyCentFV0A", false, "Centrality based on FV0A"};
@@ -206,8 +205,8 @@ struct HeavyionMultiplicity {
206205
Configurable<bool> isApplyCentFT0M{"isApplyCentFT0M", false, "Centrality based on FT0A + FT0C"};
207206
Configurable<bool> isApplyCentNGlobal{"isApplyCentNGlobal", false, "Centrality based on global tracks"};
208207
Configurable<bool> isApplyCentMFT{"isApplyCentMFT", false, "Centrality based on MFT tracks"};
209-
Configurable<bool> isApplyInelgt0{"isApplyInelgt0", false, "Enable INEL > 0 condition"};
210208
Configurable<bool> isApplyTVX{"isApplyTVX", false, "Enable TVX trigger sel"};
209+
Configurable<bool> isApplyExtraPhiCut{"isApplyExtraPhiCut", false, "Enable extra phi cut"};
211210

212211
void init(InitContext const&)
213212
{
@@ -230,14 +229,15 @@ struct HeavyionMultiplicity {
230229
auto hstat = histos.get<TH1>(HIST("EventHist"));
231230
auto* x = hstat->GetXaxis();
232231
x->SetBinLabel(1, "All events");
233-
x->SetBinLabel(2, "sel8");
234-
x->SetBinLabel(3, "kNoSameBunchPileup"); // reject collisions in case of pileup with another collision in the same foundBC
235-
x->SetBinLabel(4, "kIsGoodZvtxFT0vsPV"); // small difference between z-vertex from PV and from FT0
236-
x->SetBinLabel(5, "ApplyExtraCorrCut");
237-
x->SetBinLabel(6, "ApplyNoCollInTimeRangeStandard");
238-
x->SetBinLabel(7, "ApplyNoCollInRofStandard");
239-
x->SetBinLabel(8, "ApplyNoHighMultCollInPrevRof");
240-
x->SetBinLabel(9, "INEL > 0");
232+
x->SetBinLabel(2, "kIsTriggerTVX");
233+
x->SetBinLabel(3, "kNoTimeFrameBorder");
234+
x->SetBinLabel(4, "kNoITSROFrameBorder");
235+
x->SetBinLabel(5, "kNoSameBunchPileup"); // reject collisions in case of pileup with another collision in the same foundBC
236+
x->SetBinLabel(6, "kIsGoodZvtxFT0vsPV"); // small difference between z-vertex from PV and from FT0
237+
x->SetBinLabel(7, "ApplyNoCollInTimeRangeStandard");
238+
x->SetBinLabel(8, "ApplyNoCollInRofStandard");
239+
x->SetBinLabel(9, "ApplyNoHighMultCollInPrevRof");
240+
x->SetBinLabel(10, "INEL > 0");
241241

242242
if (doprocessData) {
243243
histos.add("hdcaxy", "dca to pv in the xy plane", kTH1D, {dcaAxis}, false);
@@ -363,50 +363,42 @@ struct HeavyionMultiplicity {
363363
bool isEventSelected(CheckCol const& col)
364364
{
365365
histos.fill(HIST("EventHist"), 1);
366-
367-
if (!col.sel8()) {
366+
if (!col.selection_bit(o2::aod::evsel::kIsTriggerTVX)) {
368367
return false;
369368
}
370369
histos.fill(HIST("EventHist"), 2);
371-
372-
if (isApplyTVX && !col.selection_bit(o2::aod::evsel::kIsTriggerTVX)) {
373-
return false;
374-
}
375-
376-
if (isApplySameBunchPileup && !col.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
370+
if (isApplyTFcut && !col.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
377371
return false;
378372
}
379373
histos.fill(HIST("EventHist"), 3);
380-
381-
if (isApplyGoodZvtxFT0vsPV && !col.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
374+
if (isApplyITSROcut && !col.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
382375
return false;
383376
}
384377
histos.fill(HIST("EventHist"), 4);
385-
386-
if (isApplyExtraCorrCut && col.multNTracksPV() > npvTracksCut && col.multFT0C() < (10 * col.multNTracksPV() - ft0cCut)) {
378+
if (isApplySameBunchPileup && !col.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
387379
return false;
388380
}
389381
histos.fill(HIST("EventHist"), 5);
390-
391-
if (isApplyNoCollInTimeRangeStandard && !col.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
382+
if (isApplyGoodZvtxFT0vsPV && !col.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
392383
return false;
393384
}
394385
histos.fill(HIST("EventHist"), 6);
395-
396-
if (isApplyNoCollInRofStandard && !col.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
386+
if (isApplyNoCollInTimeRangeStandard && !col.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
397387
return false;
398388
}
399389
histos.fill(HIST("EventHist"), 7);
400-
401-
if (isApplyNoHighMultCollInPrevRof && !col.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof)) {
390+
if (isApplyNoCollInRofStandard && !col.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
402391
return false;
403392
}
404393
histos.fill(HIST("EventHist"), 8);
405-
406-
if (isApplyInelgt0 && !col.isInelGt0()) {
394+
if (isApplyNoHighMultCollInPrevRof && !col.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof)) {
407395
return false;
408396
}
409397
histos.fill(HIST("EventHist"), 9);
398+
if (isApplyInelgt0 && !col.isInelGt0()) {
399+
return false;
400+
}
401+
histos.fill(HIST("EventHist"), 10);
410402
return true;
411403
}
412404

0 commit comments

Comments
 (0)