diff --git a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx index d7bc6ed2a1c..1e67da5fac5 100644 --- a/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx +++ b/PWGLF/Tasks/Strangeness/hStrangeCorrelation.cxx @@ -3204,6 +3204,13 @@ struct HStrangeCorrelation { validCollisions[i].reserve(masterConfigurations.mixingParameter); } } + if (masterConfigurations.mixingEventForClosureTest) { + genLevelCollisions.resize(histos.get(HIST("axes/hMultAxis"))->GetNbinsX() * + histos.get(HIST("axes/hVertexZAxis"))->GetNbinsX()); + for (size_t i = 0; i < genLevelCollisions.size(); ++i) { + genLevelCollisions[i].reserve(masterConfigurations.mixingParameter); + } + } if (!masterConfigurations.doPPAnalysis) { // event selections in Pb-Pb histos.add("hEventSelection", "hEventSelection", kTH1F, {{11, 0, 11}}); @@ -3494,7 +3501,7 @@ struct HStrangeCorrelation { if (TESTBIT(doCorrelation, i)) { histos.add(fmt::format("ClosureTest/sameEvent/{}", Particlenames[i]).c_str(), "", kTHnF, {axisDeltaPhiNDim, axisDeltaEtaNDim, axisPtAssocNDim, axisPtTriggerNDim, axisVtxZNDim, axisMultNDim}); if (masterConfigurations.mixingEventForClosureTest) { - histos.addClone("ClosureTest/sameEvent/", "ClosureTest/mixedEvent/"); + histos.add(fmt::format("ClosureTest/mixedEvent/{}", Particlenames[i]).c_str(), "", kTHnF, {axisDeltaPhiNDim, axisDeltaEtaNDim, axisPtAssocNDim, axisPtTriggerNDim, axisVtxZNDim, axisMultNDim}); } if (masterConfigurations.doClosureTestTriggerWithRecoTrackMatch) { histos.add(fmt::format("ClosureTest/TriggerWithRecoTrackMatch/sameEvent/{}", Particlenames[i]).c_str(), "truth pairs whose trigger has at least one reconstructed track with a matching MC label", kTHnF, {axisDeltaPhiNDim, axisDeltaEtaNDim, axisPtAssocNDim, axisPtTriggerNDim, axisVtxZNDim, axisMultNDim}); @@ -7026,7 +7033,8 @@ struct HStrangeCorrelation { TriggerInGen.pt = gpt; TriggerInGen.eta = geta; TriggerInGen.phi = gphi; - genLevelCollisions[binnumb].back().emplace_back(TriggerInGen); + if (gpt <= axisRanges[3][1] && gpt >= axisRanges[3][0]) + genLevelCollisions[binnumb].back().emplace_back(TriggerInGen); } if (masterConfigurations.doClosureTestTriggerWithRecoTrackMatch && mcParticleIdsWithRecoTrackMatch.contains(mcParticle.globalIndex())) { histos.fill(HIST("ClosureTest/TriggerWithRecoTrackMatch/hTrigger"), gpt, geta, bestCollisionCentpercentile); @@ -7155,7 +7163,7 @@ struct HStrangeCorrelation { if (!masterConfigurations.mixingEventForClosureTest) { return; } - for (std::size_t iCollision = 0; iCollision < genLevelCollisions[binnumb].size() - 1; iCollision++) { + for (std::size_t iCollision = 0; iCollision < genLevelCollisions[binnumb].size(); iCollision++) { for (std::size_t iTrigger = 0; iTrigger < genLevelCollisions[binnumb][iCollision].size(); iTrigger++) { auto triggerParticle = genLevelCollisions[binnumb][iCollision][iTrigger]; double getatrigger = triggerParticle.eta;