Skip to content

[PWGLF] Apply a different event selection for PbPb in the feeddown part of hStrangenessCorrelation.cxx #17951

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
lhusova merged 3 commits into from
Sep 17, 2026
+21 −12
Merged

[PWGLF] Apply a different event selection for PbPb in the feeddown part of hStrangenessCorrelation.cxx #17951

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
535a59e
Add files via upload
gengjiabin12 Sep 17, 2026
7a5d3b3
Add files via upload
gengjiabin12 Sep 17, 2026
4ce1f01
Add files via upload
gengjiabin12 Sep 17, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
33 changes: 21 additions & 12 deletions PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2882,8 +2882,8 @@ struct HStrangeCorrelation {
}
if (!masterConfigurations.doPPAnalysis) {
// event selections in Pb-Pb
histos.add("hEventSelection", "hEventSelection", kTH1F, {{10, 0, 10}});
std::array<TString, 10> eventSelLabel = {"all", "sel8", "kIsTriggerTVX", "PV_{z}", "kIsGoodITSLayersAll", "kIsGoodZvtxFT0vsPV", "OccupCut", "kNoTimeFrameBorder", "kNoITSROFrameBorder", "kNoSameBunchPileup "};
histos.add("hEventSelection", "hEventSelection", kTH1F, {{11, 0, 11}});
std::array<TString, 11> eventSelLabel = {"all", "sel8", "kIsTriggerTVX", "Centrality", "PV_{z}", "kIsGoodITSLayersAll", "kIsGoodZvtxFT0vsPV", "OccupCut", "kNoTimeFrameBorder", "kNoITSROFrameBorder", "kNoSameBunchPileup "};
for (int i = 1; i <= histos.get<TH1>(HIST("hEventSelection"))->GetNbinsX(); i++) {
histos.get<TH1>(HIST("hEventSelection"))->GetXaxis()->SetBinLabel(i, eventSelLabel[i - 1].Data());
}
Expand Down Expand Up @@ -3336,19 +3336,26 @@ struct HStrangeCorrelation {
histos.fill(HIST("hEventSelection"), 2.5 /* FT0 vertex (acceptable FT0C-FT0A time difference) collisions */);
}

if (collision.centFT0C() > axisRanges[5][1] || collision.centFT0C() < axisRanges[5][0]) {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 3.5 /* collisions after sel centrality sel*/);
}

if (std::abs(collision.posZ()) > masterConfigurations.zVertexCut) {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 3.5 /* collisions after sel pvz sel*/);
histos.fill(HIST("hEventSelection"), 4.5 /* collisions after sel pvz sel*/);
}

if (!collision.selection_bit(aod::evsel::kIsGoodITSLayersAll) && masterConfigurations.requireAllGoodITSLayers) {
// cut time intervals with dead ITS staves
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 4.5 /* collisions after cut time intervals with dead ITS staves*/);
histos.fill(HIST("hEventSelection"), 5.5 /* collisions after cut time intervals with dead ITS staves*/);
}

if (!collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV) && masterConfigurations.requireGoodZvtxFT0vsPV) {
Expand All @@ -3357,15 +3364,15 @@ struct HStrangeCorrelation {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 5.5 /* removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference*/);
histos.fill(HIST("hEventSelection"), 6.5 /* removes collisions with large differences between z of PV by tracks and z of PV from FT0 A-C time difference*/);
}

auto occupancy = collision.trackOccupancyInTimeRange();
if (occupancy < cfgCutOccupancyLow || occupancy > cfgCutOccupancyHigh) {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 6.5 /* Below min occupancy and Above max occupancy*/);
histos.fill(HIST("hEventSelection"), 7.5 /* Below min occupancy and Above max occupancy*/);
}

/*
Expand All @@ -3381,7 +3388,7 @@ struct HStrangeCorrelation {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 7.5 /* reject collisions close to Time Frame borders*/);
histos.fill(HIST("hEventSelection"), 8.5 /* reject collisions close to Time Frame borders*/);
}

if (!collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) {
Expand All @@ -3390,7 +3397,7 @@ struct HStrangeCorrelation {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 8.5 /* reject events affected by the ITS ROF border*/);
histos.fill(HIST("hEventSelection"), 9.5 /* reject events affected by the ITS ROF border*/);
}

if (!collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) {
Expand All @@ -3399,7 +3406,7 @@ struct HStrangeCorrelation {
return false;
}
if (fillHists) {
histos.fill(HIST("hEventSelection"), 9.5 /* rejects collisions which are associated with the same "found-by-T0" bunch crossing*/);
histos.fill(HIST("hEventSelection"), 10.5 /* rejects collisions which are associated with the same "found-by-T0" bunch crossing*/);
}
return true;
}
Expand Down Expand Up @@ -6402,12 +6409,12 @@ struct HStrangeCorrelation {
}
}

void processFeedDown(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>::iterator const& collision, aod::AssocV0s const& associatedV0s, aod::McParticles const&, V0DatasWithoutTrackXMC const&, TracksComplete const&, aod::BCsWithTimestamps const&)
void processFeedDown(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::CentFT0Cs, aod::PVMults>::iterator const& collision, aod::AssocV0s const& associatedV0s, aod::McParticles const&, V0DatasWithoutTrackXMC const&, TracksComplete const&, aod::BCsWithTimestamps const&)
{

// ________________________________________________
// Perform basic event selection
if (!isCollisionSelected(collision)) {
if (!(masterConfigurations.doPPAnalysis ? isCollisionSelected(collision) : isCollisionSelectedPbPb(collision, false))) {
return;
}

Expand All @@ -6420,7 +6427,9 @@ struct HStrangeCorrelation {
if (postrack.tpcNClsCrossedRows() < trackSelection.minTPCNCrossedRowsAssociated || negtrack.tpcNClsCrossedRows() < trackSelection.minTPCNCrossedRowsAssociated) {
continue;
}

if (trackSelection.requireClusterInITS && (postrack.itsNCls() < trackSelection.minITSClustersForDaughterTracks || negtrack.itsNCls() < trackSelection.minITSClustersForDaughterTracks)) {
continue;
}
//---] syst cuts [---
if (v0Data.v0radius() < v0Selection.v0RadiusMin || v0Data.v0radius() > v0Selection.v0RadiusMax ||
std::abs(v0Data.dcapostopv()) < v0Selection.dcapostopv || std::abs(v0Data.dcanegtopv()) < v0Selection.dcanegtopv ||
Expand Down
Loading