Skip to content

Commit 9690a59

Browse files
LuzhiyonggLuzhiyongg
authored
[PWGCF] add kNoITSROFrameBorder and kNoTimeFrameBorder EvSel (#11943)
1 parent 9a1890a commit 9690a59

File tree

2 files changed

+85
-51
lines changed

2 files changed

+85
-51
lines changed

‎PWGCF/Flow/Tasks/flowTask.cxx‎

Lines changed: 54 additions & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -75,6 +75,8 @@ struct FlowTask {
7575
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalEventCut, bool, false, "Use additional event cut on mult correlations")
7676
O2_DEFINE_CONFIGURABLE(cfgUseTentativeEventCounter, bool, false, "After sel8(), count events regardless of real event selection")
7777
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoSameBunchPileup, bool, false, "rejects collisions which are associated with the same found-by-T0 bunch crossing")
78+
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoITSROFrameBorder, bool, false, "reject events at ITS ROF border")
79+
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoTimeFrameBorder, bool, false, "reject events at TF border")
7880
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodZvtxFT0vsPV, bool, false, "removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference, use this cut at low multiplicities with caution")
7981
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoCollInTimeRangeStandard, bool, false, "no collisions in specified time range")
8082
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodITSLayersAll, bool, true, "cut time intervals with dead ITS staves")
@@ -207,29 +209,33 @@ struct FlowTask {
207209
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(3, "after supicious Runs removal");
208210
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(4, "after additional event cut");
209211
registry.get<TH1>(HIST("hEventCount"))->GetXaxis()->SetBinLabel(5, "after correction loads");
210-
registry.add("hEventCountSpecific", "Number of Event;; Count", {HistType::kTH1D, {{10, 0, 10}}});
212+
registry.add("hEventCountSpecific", "Number of Event;; Count", {HistType::kTH1D, {{12, 0, 12}}});
211213
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(1, "after sel8");
212214
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(2, "kNoSameBunchPileup");
213-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(3, "kIsGoodZvtxFT0vsPV");
214-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(4, "kNoCollInTimeRangeStandard");
215-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(5, "kIsGoodITSLayersAll");
216-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(6, "kNoCollInRofStandard");
217-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(7, "kNoHighMultCollInPrevRof");
218-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(8, "occupancy");
219-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(9, "MultCorrelation");
220-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(10, "cfgEvSelV0AT0ACut");
215+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(3, "kNoITSROFrameBorder");
216+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(4, "kNoTimeFrameBorder");
217+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(5, "kIsGoodZvtxFT0vsPV");
218+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(6, "kNoCollInTimeRangeStandard");
219+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(7, "kIsGoodITSLayersAll");
220+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(8, "kNoCollInRofStandard");
221+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(9, "kNoHighMultCollInPrevRof");
222+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(10, "occupancy");
223+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(11, "MultCorrelation");
224+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(12, "cfgEvSelV0AT0ACut");
221225
if (cfgUseTentativeEventCounter) {
222-
registry.add("hEventCountTentative", "Number of Event;; Count", {HistType::kTH1D, {{10, 0, 10}}});
226+
registry.add("hEventCountTentative", "Number of Event;; Count", {HistType::kTH1D, {{12, 0, 12}}});
223227
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(1, "after sel8");
224228
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(2, "kNoSameBunchPileup");
225-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(3, "kIsGoodZvtxFT0vsPV");
226-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(4, "kNoCollInTimeRangeStandard");
227-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(5, "kIsGoodITSLayersAll");
228-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(6, "kNoCollInRofStandard");
229-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(7, "kNoHighMultCollInPrevRof");
230-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(8, "occupancy");
231-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(9, "MultCorrelation");
232-
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(10, "cfgEvSelV0AT0ACut");
229+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(3, "kNoITSROFrameBorder");
230+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(4, "kNoTimeFrameBorder");
231+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(5, "kIsGoodZvtxFT0vsPV");
232+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(6, "kNoCollInTimeRangeStandard");
233+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(7, "kIsGoodITSLayersAll");
234+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(8, "kNoCollInRofStandard");
235+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(9, "kNoHighMultCollInPrevRof");
236+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(10, "occupancy");
237+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(11, "MultCorrelation");
238+
registry.get<TH1>(HIST("hEventCountTentative"))->GetXaxis()->SetBinLabel(12, "cfgEvSelV0AT0ACut");
233239
}
234240
registry.add("hVtxZ", "Vexter Z distribution", {HistType::kTH1D, {axisVertex}});
235241
registry.add("hMult", "Multiplicity distribution", {HistType::kTH1D, {{3000, 0.5, 3000.5}}});
@@ -589,44 +595,54 @@ struct FlowTask {
589595
}
590596
if (cfgEvSelkNoSameBunchPileup)
591597
registry.fill(HIST("hEventCountSpecific"), 1.5);
598+
if (cfgEvSelkNoITSROFrameBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
599+
return 0;
600+
}
601+
if (cfgEvSelkNoITSROFrameBorder)
602+
registry.fill(HIST("hEventCountSpecific"), 2.5);
603+
if (cfgEvSelkNoTimeFrameBorder && !collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
604+
return 0;
605+
}
606+
if (cfgEvSelkNoTimeFrameBorder)
607+
registry.fill(HIST("hEventCountSpecific"), 3.5);
592608
if (cfgEvSelkIsGoodZvtxFT0vsPV && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
593609
// removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference
594610
// use this cut at low multiplicities with caution
595611
return 0;
596612
}
597613
if (cfgEvSelkIsGoodZvtxFT0vsPV)
598-
registry.fill(HIST("hEventCountSpecific"), 2.5);
614+
registry.fill(HIST("hEventCountSpecific"), 4.5);
599615
if (cfgEvSelkNoCollInTimeRangeStandard && !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
600616
// no collisions in specified time range
601617
return 0;
602618
}
603619
if (cfgEvSelkNoCollInTimeRangeStandard)
604-
registry.fill(HIST("hEventCountSpecific"), 3.5);
620+
registry.fill(HIST("hEventCountSpecific"), 5.5);
605621
if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
606622
// from Jan 9 2025 AOT meeting
607623
// cut time intervals with dead ITS staves
608624
return 0;
609625
}
610626
if (cfgEvSelkIsGoodITSLayersAll)
611-
registry.fill(HIST("hEventCountSpecific"), 4.5);
627+
registry.fill(HIST("hEventCountSpecific"), 6.5);
612628
if (cfgEvSelkNoCollInRofStandard && !collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
613629
// no other collisions in this Readout Frame with per-collision multiplicity above threshold
614630
return 0;
615631
}
616632
if (cfgEvSelkNoCollInRofStandard)
617-
registry.fill(HIST("hEventCountSpecific"), 5.5);
633+
registry.fill(HIST("hEventCountSpecific"), 7.5);
618634
if (cfgEvSelkNoHighMultCollInPrevRof && !collision.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof)) {
619635
// veto an event if FT0C amplitude in previous ITS ROF is above threshold
620636
return 0;
621637
}
622638
if (cfgEvSelkNoHighMultCollInPrevRof)
623-
registry.fill(HIST("hEventCountSpecific"), 6.5);
639+
registry.fill(HIST("hEventCountSpecific"), 8.5);
624640
auto multNTracksPV = collision.multNTracksPV();
625641
auto occupancy = collision.trackOccupancyInTimeRange();
626642
if (cfgEvSelOccupancy && (occupancy < cfgCutOccupancyLow || occupancy > cfgCutOccupancyHigh))
627643
return 0;
628644
if (cfgEvSelOccupancy)
629-
registry.fill(HIST("hEventCountSpecific"), 7.5);
645+
registry.fill(HIST("hEventCountSpecific"), 9.5);
630646

631647
if (cfgEvSelMultCorrelation) {
632648
if (multNTracksPV < fMultPVCutLow->Eval(centrality))
@@ -639,14 +655,14 @@ struct FlowTask {
639655
return 0;
640656
}
641657
if (cfgEvSelMultCorrelation)
642-
registry.fill(HIST("hEventCountSpecific"), 8.5);
658+
registry.fill(HIST("hEventCountSpecific"), 10.5);
643659

644660
// V0A T0A 5 sigma cut
645661
float sigma = 5.0;
646662
if (cfgEvSelV0AT0ACut && (std::fabs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > sigma * fT0AV0ASigma->Eval(collision.multFT0A())))
647663
return 0;
648664
if (cfgEvSelV0AT0ACut)
649-
registry.fill(HIST("hEventCountSpecific"), 9.5);
665+
registry.fill(HIST("hEventCountSpecific"), 11.5);
650666

651667
return 1;
652668
}
@@ -658,25 +674,29 @@ struct FlowTask {
658674
// Regradless of the event selection, fill the event counter histograms
659675
if (collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup))
660676
registry.fill(HIST("hEventCountTentative"), 1.5);
661-
if (collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV))
677+
if (collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))
662678
registry.fill(HIST("hEventCountTentative"), 2.5);
663-
if (collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard))
679+
if (collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder))
664680
registry.fill(HIST("hEventCountTentative"), 3.5);
665-
if (collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll))
681+
if (collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV))
666682
registry.fill(HIST("hEventCountTentative"), 4.5);
667-
if (collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard))
683+
if (collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard))
668684
registry.fill(HIST("hEventCountTentative"), 5.5);
669-
if (collision.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof))
685+
if (collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll))
670686
registry.fill(HIST("hEventCountTentative"), 6.5);
687+
if (collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard))
688+
registry.fill(HIST("hEventCountTentative"), 7.5);
689+
if (collision.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof))
690+
registry.fill(HIST("hEventCountTentative"), 8.5);
671691
auto multNTracksPV = collision.multNTracksPV();
672692
auto occupancy = collision.trackOccupancyInTimeRange();
673693
if (!(occupancy < cfgCutOccupancyLow || occupancy > cfgCutOccupancyHigh))
674-
registry.fill(HIST("hEventCountTentative"), 7.5);
694+
registry.fill(HIST("hEventCountTentative"), 9.5);
675695
if (!((multNTracksPV < fMultPVCutLow->Eval(centrality)) || (multNTracksPV > fMultPVCutHigh->Eval(centrality)) || (multTrk < fMultCutLow->Eval(centrality)) || (multTrk > fMultCutHigh->Eval(centrality))))
676-
registry.fill(HIST("hEventCountTentative"), 8.5);
696+
registry.fill(HIST("hEventCountTentative"), 10.5);
677697
float sigma = 5.0;
678698
if (!(std::fabs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > sigma * fT0AV0ASigma->Eval(collision.multFT0A())))
679-
registry.fill(HIST("hEventCountTentative"), 9.5);
699+
registry.fill(HIST("hEventCountTentative"), 11.5);
680700
}
681701

682702
template <typename TTrack>

‎PWGCF/TwoParticleCorrelations/Tasks/diHadronCor.cxx‎

Lines changed: 31 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -80,6 +80,8 @@ struct DiHadronCor {
8080
O2_DEFINE_CONFIGURABLE(cfgCentTableUnavailable, bool, false, "if a dataset does not provide centrality information")
8181
O2_DEFINE_CONFIGURABLE(cfgUseAdditionalEventCut, bool, false, "Use additional event cut on mult correlations")
8282
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoSameBunchPileup, bool, false, "rejects collisions which are associated with the same found-by-T0 bunch crossing")
83+
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoITSROFrameBorder, bool, false, "reject events at ITS ROF border")
84+
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoTimeFrameBorder, bool, false, "reject events at TF border")
8385
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodZvtxFT0vsPV, bool, false, "removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference, use this cut at low multiplicities with caution")
8486
O2_DEFINE_CONFIGURABLE(cfgEvSelkNoCollInTimeRangeStandard, bool, false, "no collisions in specified time range")
8587
O2_DEFINE_CONFIGURABLE(cfgEvSelkIsGoodITSLayersAll, bool, true, "cut time intervals with dead ITS staves")
@@ -188,17 +190,19 @@ struct DiHadronCor {
188190

189191
// Event Counter
190192
if (doprocessSame && cfgUseAdditionalEventCut) {
191-
registry.add("hEventCountSpecific", "Number of Event;; Count", {HistType::kTH1D, {{10, 0, 10}}});
193+
registry.add("hEventCountSpecific", "Number of Event;; Count", {HistType::kTH1D, {{12, 0, 12}}});
192194
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(1, "after sel8");
193195
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(2, "kNoSameBunchPileup");
194-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(3, "kIsGoodZvtxFT0vsPV");
195-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(4, "kNoCollInTimeRangeStandard");
196-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(5, "kIsGoodITSLayersAll");
197-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(6, "kNoCollInRofStandard");
198-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(7, "kNoHighMultCollInPrevRof");
199-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(8, "occupancy");
200-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(9, "MultCorrelation");
201-
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(10, "cfgEvSelV0AT0ACut");
196+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(3, "kNoITSROFrameBorder");
197+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(4, "kNoTimeFrameBorder");
198+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(5, "kIsGoodZvtxFT0vsPV");
199+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(6, "kNoCollInTimeRangeStandard");
200+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(7, "kIsGoodITSLayersAll");
201+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(8, "kNoCollInRofStandard");
202+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(9, "kNoHighMultCollInPrevRof");
203+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(10, "occupancy");
204+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(11, "MultCorrelation");
205+
registry.get<TH1>(HIST("hEventCountSpecific"))->GetXaxis()->SetBinLabel(12, "cfgEvSelV0AT0ACut");
202206
}
203207

204208
if (cfgUseAdditionalEventCut) {
@@ -578,43 +582,53 @@ struct DiHadronCor {
578582
}
579583
if (fillCounter && cfgEvSelkNoSameBunchPileup)
580584
registry.fill(HIST("hEventCountSpecific"), 1.5);
585+
if (cfgEvSelkNoITSROFrameBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
586+
return 0;
587+
}
588+
if (fillCounter && cfgEvSelkNoITSROFrameBorder)
589+
registry.fill(HIST("hEventCountSpecific"), 2.5);
590+
if (cfgEvSelkNoTimeFrameBorder && !collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) {
591+
return 0;
592+
}
593+
if (fillCounter && cfgEvSelkNoTimeFrameBorder)
594+
registry.fill(HIST("hEventCountSpecific"), 3.5);
581595
if (cfgEvSelkIsGoodZvtxFT0vsPV && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) {
582596
// removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference
583597
// use this cut at low multiplicities with caution
584598
return 0;
585599
}
586600
if (fillCounter && cfgEvSelkIsGoodZvtxFT0vsPV)
587-
registry.fill(HIST("hEventCountSpecific"), 2.5);
601+
registry.fill(HIST("hEventCountSpecific"), 4.5);
588602
if (cfgEvSelkNoCollInTimeRangeStandard && !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
589603
// no collisions in specified time range
590604
return 0;
591605
}
592606
if (fillCounter && cfgEvSelkNoCollInTimeRangeStandard)
593-
registry.fill(HIST("hEventCountSpecific"), 3.5);
607+
registry.fill(HIST("hEventCountSpecific"), 5.5);
594608
if (cfgEvSelkIsGoodITSLayersAll && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
595609
// from Jan 9 2025 AOT meeting
596610
// cut time intervals with dead ITS staves
597611
return 0;
598612
}
599613
if (fillCounter && cfgEvSelkIsGoodITSLayersAll)
600-
registry.fill(HIST("hEventCountSpecific"), 4.5);
614+
registry.fill(HIST("hEventCountSpecific"), 6.5);
601615
if (cfgEvSelkNoCollInRofStandard && !collision.selection_bit(o2::aod::evsel::kNoCollInRofStandard)) {
602616
// no other collisions in this Readout Frame with per-collision multiplicity above threshold
603617
return 0;
604618
}
605619
if (fillCounter && cfgEvSelkNoCollInRofStandard)
606-
registry.fill(HIST("hEventCountSpecific"), 5.5);
620+
registry.fill(HIST("hEventCountSpecific"), 7.5);
607621
if (cfgEvSelkNoHighMultCollInPrevRof && !collision.selection_bit(o2::aod::evsel::kNoHighMultCollInPrevRof)) {
608622
// veto an event if FT0C amplitude in previous ITS ROF is above threshold
609623
return 0;
610624
}
611625
if (fillCounter && cfgEvSelkNoHighMultCollInPrevRof)
612-
registry.fill(HIST("hEventCountSpecific"), 6.5);
626+
registry.fill(HIST("hEventCountSpecific"), 8.5);
613627
auto occupancy = collision.trackOccupancyInTimeRange();
614628
if (cfgEvSelOccupancy && (occupancy < cfgCutOccupancyLow || occupancy > cfgCutOccupancyHigh))
615629
return 0;
616630
if (fillCounter && cfgEvSelOccupancy)
617-
registry.fill(HIST("hEventCountSpecific"), 7.5);
631+
registry.fill(HIST("hEventCountSpecific"), 9.5);
618632

619633
auto multNTracksPV = collision.multNTracksPV();
620634
if (cfgEvSelMultCorrelation) {
@@ -628,14 +642,14 @@ struct DiHadronCor {
628642
return 0;
629643
}
630644
if (fillCounter && cfgEvSelMultCorrelation)
631-
registry.fill(HIST("hEventCountSpecific"), 8.5);
645+
registry.fill(HIST("hEventCountSpecific"), 10.5);
632646

633647
// V0A T0A 5 sigma cut
634648
float sigma = 5.0;
635649
if (cfgEvSelV0AT0ACut && (std::fabs(collision.multFV0A() - fT0AV0AMean->Eval(collision.multFT0A())) > sigma * fT0AV0ASigma->Eval(collision.multFT0A())))
636650
return 0;
637651
if (fillCounter && cfgEvSelV0AT0ACut)
638-
registry.fill(HIST("hEventCountSpecific"), 9.5);
652+
registry.fill(HIST("hEventCountSpecific"), 11.5);
639653

640654
return 1;
641655
}

0 commit comments

Comments
 (0)