From 3ad0004bf61d5b4d3142c9d9fe67d2d8bc0e13f8 Mon Sep 17 00:00:00 2001 From: ddobrigk Date: Mon, 8 Jul 2024 10:38:17 +0200 Subject: [PATCH 1/4] PWGLF: add INEL>0 selection to hStr corr analysis --- .../Tasks/Strangeness/hStrangeCorrelation.cxx | 18 +++++++++++++++++- 1 file changed, 17 insertions(+), 1 deletion(-) diff --git a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx index 9408bfe2a07..acefe706430 100644 --- a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx +++ b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx @@ -59,6 +59,7 @@ struct correlateStrangeness { Configurable doCorrelationOmegaMinus{"doCorrelationOmegaMinus", false, "do OmegaMinus correlation"}; Configurable doCorrelationOmegaPlus{"doCorrelationOmegaPlus", false, "do OmegaPlus correlation"}; Configurable doCorrelationPion{"doCorrelationPion", false, "do Pion correlation"}; + Configurable selectINELgtZERO{"selectINELgtZERO", true, "select INEL>0 events"}; Configurable zVertexCut{"zVertexCut", 10, "Cut on PV position"}; Configurable skipUnderOverflowInTHn{"skipUnderOverflowInTHn", false, "skip under/overflow in THns"}; Configurable mixingParameter{"mixingParameter", 10, "how many events are mixed"}; @@ -650,6 +651,9 @@ struct correlateStrangeness { if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) { return; } + if (!collision.isInelGt0()&&selectINELgtZERO){ + return; + } // ________________________________________________ if (!doprocessSameEventHCascades) { histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()})); @@ -709,6 +713,9 @@ struct correlateStrangeness { if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) { return; } + if (!collision.isInelGt0()&&selectINELgtZERO){ + return; + } // ________________________________________________ histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()})); histos.fill(HIST("EventQA/hMult"), collision.centFT0M()); @@ -749,7 +756,7 @@ struct correlateStrangeness { // Do hadron - cascade correlations fillCorrelationsCascade(triggerTracks, associatedCascades, false, collision.posZ(), collision.centFT0M()); } - void processSameEventHPions(soa::Join::iterator const& collision, + void processSameEventHPions(soa::Join::iterator const& collision, aod::AssocPions const& associatedPions, aod::TriggerTracks const& triggerTracks, TracksComplete const&, aod::BCsWithTimestamps const&) { @@ -764,6 +771,9 @@ struct correlateStrangeness { if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) { return; } + if (!collision.isInelGt0()&&selectINELgtZERO){ + return; + } // ________________________________________________ if (!doprocessSameEventHCascades && !doprocessSameEventHV0s) { histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()})); @@ -803,6 +813,8 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; + if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){ + continue; if (!doprocessMixedEventHCascades) { if (collision1.globalIndex() == collision2.globalIndex()) { @@ -838,6 +850,8 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; + if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){ + continue; if (collision1.globalIndex() == collision2.globalIndex()) { histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting @@ -870,6 +884,8 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; + if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){ + continue; if (collision1.globalIndex() == collision2.globalIndex()) { histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting From d2461c60cf9257b60d08557f0807a8fcbce273a8 Mon Sep 17 00:00:00 2001 From: ALICE Builder Date: Mon, 8 Jul 2024 10:40:25 +0200 Subject: [PATCH 2/4] Please consider the following formatting changes (#309) --- .../Tasks/Strangeness/hStrangeCorrelation.cxx | 96 +++++++++---------- 1 file changed, 48 insertions(+), 48 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx index acefe706430..619af5a9e0e 100644 --- a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx +++ b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx @@ -651,7 +651,7 @@ struct correlateStrangeness { if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) { return; } - if (!collision.isInelGt0()&&selectINELgtZERO){ + if (!collision.isInelGt0() && selectINELgtZERO) { return; } // ________________________________________________ @@ -713,7 +713,7 @@ struct correlateStrangeness { if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) { return; } - if (!collision.isInelGt0()&&selectINELgtZERO){ + if (!collision.isInelGt0() && selectINELgtZERO) { return; } // ________________________________________________ @@ -771,7 +771,7 @@ struct correlateStrangeness { if (collision.centFT0M() > axisRanges[5][1] || collision.centFT0M() < axisRanges[5][0]) { return; } - if (!collision.isInelGt0()&&selectINELgtZERO){ + if (!collision.isInelGt0() && selectINELgtZERO) { return; } // ________________________________________________ @@ -813,25 +813,25 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; - if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){ + if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) { continue; - if (!doprocessMixedEventHCascades) { - if (collision1.globalIndex() == collision2.globalIndex()) { - histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting + if (!doprocessMixedEventHCascades) { + if (collision1.globalIndex() == collision2.globalIndex()) { + histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting + } + histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); + histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); + histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); } - histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); - histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); - histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); + // ________________________________________________ + // Do slicing + auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); + auto slicedAssocV0s = associatedV0s.sliceBy(collisionSliceV0s, collision2.globalIndex()); + // ________________________________________________ + // Do hadron - V0 correlations + fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M()); } - // ________________________________________________ - // Do slicing - auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); - auto slicedAssocV0s = associatedV0s.sliceBy(collisionSliceV0s, collision2.globalIndex()); - // ________________________________________________ - // Do hadron - V0 correlations - fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M()); - } } void processMixedEventHCascades(soa::Join const& collisions, aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks, @@ -850,24 +850,24 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; - if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){ + if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) { continue; - if (collision1.globalIndex() == collision2.globalIndex()) { - histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting - } + if (collision1.globalIndex() == collision2.globalIndex()) { + histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting + } - histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); - histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); - histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); - // ________________________________________________ - // Do slicing - auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); - auto slicedAssocCascades = associatedCascades.sliceBy(collisionSliceCascades, collision2.globalIndex()); - // ________________________________________________ - // Do hadron - cascade correlations - fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M()); - } + histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); + histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); + histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); + // ________________________________________________ + // Do slicing + auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); + auto slicedAssocCascades = associatedCascades.sliceBy(collisionSliceCascades, collision2.globalIndex()); + // ________________________________________________ + // Do hadron - cascade correlations + fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M()); + } } void processMixedEventHPions(soa::Join const& collisions, aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks, @@ -884,24 +884,24 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; - if ((!collision1.isInelGt0()||!collision2.isInelGt0())&&selectINELgtZERO){ + if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) { continue; - if (collision1.globalIndex() == collision2.globalIndex()) { - histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting - } + if (collision1.globalIndex() == collision2.globalIndex()) { + histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting + } - histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); - histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); - histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); - // ________________________________________________ - // Do slicing - auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); - auto slicedAssocPions = assocPions.sliceBy(collisionSlicePions, collision2.globalIndex()); - // ________________________________________________ - // Do hadron - cascade correlations - fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M()); - } + histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); + histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); + histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); + // ________________________________________________ + // Do slicing + auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); + auto slicedAssocPions = assocPions.sliceBy(collisionSlicePions, collision2.globalIndex()); + // ________________________________________________ + // Do hadron - cascade correlations + fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M()); + } } void processMCGenerated(aod::McCollision const&, soa::SmallGroups> const& collisions, aod::McParticles const& mcParticles) { From 2d13491c7761e43d900b9efe5b787bf4460bcffa Mon Sep 17 00:00:00 2001 From: ddobrigk Date: Tue, 9 Jul 2024 11:41:13 +0200 Subject: [PATCH 3/4] Bugfixes --- .../Tasks/Strangeness/hStrangeCorrelation.cxx | 22 +++++++++++-------- 1 file changed, 13 insertions(+), 9 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx index 619af5a9e0e..aad596935ce 100644 --- a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx +++ b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx @@ -636,7 +636,7 @@ struct correlateStrangeness { } } - void processSameEventHV0s(soa::Join::iterator const& collision, + void processSameEventHV0s(soa::Join::iterator const& collision, aod::AssocV0s const& associatedV0s, aod::TriggerTracks const& triggerTracks, V0DatasWithoutTrackX const&, aod::V0sLinked const&, TracksComplete const&, aod::BCsWithTimestamps const&) { @@ -698,7 +698,7 @@ struct correlateStrangeness { fillCorrelationsV0(triggerTracks, associatedV0s, false, collision.posZ(), collision.centFT0M()); } - void processSameEventHCascades(soa::Join::iterator const& collision, + void processSameEventHCascades(soa::Join::iterator const& collision, aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks, V0DatasWithoutTrackX const&, aod::V0sLinked const&, aod::CascDatas const&, TracksComplete const&, aod::BCsWithTimestamps const&) { @@ -796,7 +796,7 @@ struct correlateStrangeness { // Do hadron - Pion correlations fillCorrelationsPion(triggerTracks, associatedPions, false, collision.posZ(), collision.centFT0M()); } - void processMixedEventHV0s(soa::Join const& collisions, + void processMixedEventHV0s(soa::Join const& collisions, aod::AssocV0s const& associatedV0s, aod::TriggerTracks const& triggerTracks, V0DatasWithoutTrackX const&, aod::V0sLinked const&, TracksComplete const&, aod::BCsWithTimestamps const&) { @@ -813,7 +813,7 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; - if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) { + if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) continue; if (!doprocessMixedEventHCascades) { @@ -833,7 +833,7 @@ struct correlateStrangeness { fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M()); } } - void processMixedEventHCascades(soa::Join const& collisions, + void processMixedEventHCascades(soa::Join const& collisions, aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks, V0DatasWithoutTrackX const&, aod::V0sLinked const&, aod::CascDatas const&, TracksComplete const&, aod::BCsWithTimestamps const&) { @@ -850,7 +850,7 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; - if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) { + if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) continue; if (collision1.globalIndex() == collision2.globalIndex()) { @@ -869,7 +869,7 @@ struct correlateStrangeness { fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M()); } } - void processMixedEventHPions(soa::Join const& collisions, + void processMixedEventHPions(soa::Join const& collisions, aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks, TracksComplete const&) { @@ -884,7 +884,7 @@ struct correlateStrangeness { continue; if (collision2.centFT0M() > axisRanges[5][1] || collision2.centFT0M() < axisRanges[5][0]) continue; - if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) { + if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) continue; if (collision1.globalIndex() == collision2.globalIndex()) { @@ -903,7 +903,7 @@ struct correlateStrangeness { fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M()); } } - void processMCGenerated(aod::McCollision const&, soa::SmallGroups> const& collisions, aod::McParticles const& mcParticles) + void processMCGenerated(aod::McCollision const&, soa::SmallGroups> const& collisions, aod::McParticles const& mcParticles) { for (auto const& mcParticle : mcParticles) { if (!mcParticle.isPhysicalPrimary()) @@ -934,12 +934,14 @@ struct correlateStrangeness { int bestCollisionFT0Mpercentile = -1; float bestCollisionVtxZ = 0.0f; bool bestCollisionSel8 = false; + bool bestCollisionINELgtZERO = false; for (auto& collision : collisions) { if (biggestNContribs < collision.numContrib()) { biggestNContribs = collision.numContrib(); bestCollisionFT0Mpercentile = collision.centFT0M(); bestCollisionSel8 = collision.sel8(); bestCollisionVtxZ = collision.posZ(); + bestCollisionINELgtZERO = collision.isInelGt0(); } } @@ -975,6 +977,8 @@ struct correlateStrangeness { return; if (std::abs(bestCollisionVtxZ) > 10.0f) return; + if (!bestCollisionINELgtZERO) + return; for (auto const& mcParticle : mcParticles) { if (!mcParticle.isPhysicalPrimary()) From 90ce4e10f10ef4ad88e11c3715ab506363b0d357 Mon Sep 17 00:00:00 2001 From: ALICE Builder Date: Tue, 9 Jul 2024 21:42:30 +0200 Subject: [PATCH 4/4] Please consider the following formatting changes (#310) --- .../Tasks/Strangeness/hStrangeCorrelation.cxx | 84 +++++++++---------- 1 file changed, 42 insertions(+), 42 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx index 3a289572d7d..6d07430c300 100644 --- a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx +++ b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx @@ -856,22 +856,22 @@ struct correlateStrangeness { if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) continue; - if (!doprocessMixedEventHCascades) { - if (collision1.globalIndex() == collision2.globalIndex()) { - histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting - } - histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); - histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); - histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); + if (!doprocessMixedEventHCascades) { + if (collision1.globalIndex() == collision2.globalIndex()) { + histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting } - // ________________________________________________ - // Do slicing - auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); - auto slicedAssocV0s = associatedV0s.sliceBy(collisionSliceV0s, collision2.globalIndex()); - // ________________________________________________ - // Do hadron - V0 correlations - fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M()); + histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); + histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); + histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); } + // ________________________________________________ + // Do slicing + auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); + auto slicedAssocV0s = associatedV0s.sliceBy(collisionSliceV0s, collision2.globalIndex()); + // ________________________________________________ + // Do hadron - V0 correlations + fillCorrelationsV0(slicedTriggerTracks, slicedAssocV0s, true, collision1.posZ(), collision1.centFT0M()); + } } void processMixedEventHCascades(soa::Join const& collisions, aod::AssocV0s const&, aod::AssocCascades const& associatedCascades, aod::TriggerTracks const& triggerTracks, @@ -895,21 +895,21 @@ struct correlateStrangeness { if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) continue; - if (collision1.globalIndex() == collision2.globalIndex()) { - histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting - } - - histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); - histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); - histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); - // ________________________________________________ - // Do slicing - auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); - auto slicedAssocCascades = associatedCascades.sliceBy(collisionSliceCascades, collision2.globalIndex()); - // ________________________________________________ - // Do hadron - cascade correlations - fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M()); + if (collision1.globalIndex() == collision2.globalIndex()) { + histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting } + + histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); + histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); + histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); + // ________________________________________________ + // Do slicing + auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); + auto slicedAssocCascades = associatedCascades.sliceBy(collisionSliceCascades, collision2.globalIndex()); + // ________________________________________________ + // Do hadron - cascade correlations + fillCorrelationsCascade(slicedTriggerTracks, slicedAssocCascades, true, collision1.posZ(), collision1.centFT0M()); + } } void processMixedEventHPions(soa::Join const& collisions, aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks, @@ -929,21 +929,21 @@ struct correlateStrangeness { if ((!collision1.isInelGt0() || !collision2.isInelGt0()) && selectINELgtZERO) continue; - if (collision1.globalIndex() == collision2.globalIndex()) { - histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting - } - - histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); - histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); - histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); - // ________________________________________________ - // Do slicing - auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); - auto slicedAssocPions = assocPions.sliceBy(collisionSlicePions, collision2.globalIndex()); - // ________________________________________________ - // Do hadron - cascade correlations - fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M()); + if (collision1.globalIndex() == collision2.globalIndex()) { + histos.fill(HIST("MixingQA/hMixingQA"), 0.0f); // same-collision pair counting } + + histos.fill(HIST("MixingQA/hMEpvz1"), collision1.posZ()); + histos.fill(HIST("MixingQA/hMEpvz2"), collision2.posZ()); + histos.fill(HIST("MixingQA/hMECollisionBins"), colBinning.getBin({collision1.posZ(), collision1.centFT0M()})); + // ________________________________________________ + // Do slicing + auto slicedTriggerTracks = triggerTracks.sliceBy(collisionSliceTracks, collision1.globalIndex()); + auto slicedAssocPions = assocPions.sliceBy(collisionSlicePions, collision2.globalIndex()); + // ________________________________________________ + // Do hadron - cascade correlations + fillCorrelationsPion(slicedTriggerTracks, slicedAssocPions, true, collision1.posZ(), collision1.centFT0M()); + } } void processMCGenerated(aod::McCollision const& mccollision, soa::SmallGroups> const& collisions, aod::McParticles const& mcParticles)