From 14265bf7a8a5675ef0eb923dd87646f5a8637ad4 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Mon, 3 Nov 2025 15:47:01 +0100 Subject: [PATCH] Fix bug in filling deta-dphi at lower and upper limit of kT --- .../Core/FemtoUniverseDetaDphiStar.h | 20 +++++-------------- ...irTaskTrackTrackSpherHarMultKtExtended.cxx | 15 +++++++++++--- 2 files changed, 17 insertions(+), 18 deletions(-) diff --git a/PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h b/PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h index aa57c5b4755..489655c3342 100644 --- a/PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h +++ b/PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h @@ -200,7 +200,7 @@ class FemtoUniverseDetaDphiStar if constexpr (kPartOneType == o2::aod::femtouniverseparticle::ParticleType::kTrack && kPartTwoType == o2::aod::femtouniverseparticle::ParticleType::kTrack) { std::string dirName = static_cast(DirNames[0]); - for (int j = 0; j < static_cast(ktBins.size() - 1); j++) { + for (int j = 1; j < static_cast(ktBins.size() - 1); j++) { std::string histSuffixkT1 = std::to_string(static_cast(ktBins[j] * 100.0)); std::string histSuffixkT2 = std::to_string(static_cast(ktBins[j + 1] * 100.0)); std::string histFolderkT = "kT_" + histSuffixkT1 + "_" + histSuffixkT2 + "/"; @@ -655,15 +655,15 @@ class FemtoUniverseDetaDphiStar return false; } - int ktbinval = -1; - if (ktval >= ktBins[0] && ktval < ktBins[1]) { - ktbinval = 0; - } else if (ktval >= ktBins[1] && ktval < ktBins[2]) { + int ktbinval = 1; + if (ktval >= ktBins[1] && ktval < ktBins[2]) { ktbinval = 1; } else if (ktval >= ktBins[2] && ktval < ktBins[3]) { ktbinval = 2; } else if (ktval >= ktBins[3] && ktval < ktBins[4]) { ktbinval = 3; + } else if (ktval >= ktBins[4] && ktval < ktBins[5]) { + ktbinval = 4; } auto deta = part1.eta() - part2.eta(); @@ -676,16 +676,6 @@ class FemtoUniverseDetaDphiStar LOG(fatal) << "FemtoUniverseDetaDphiStar: passed arguments don't agree with FemtoUniverseDetaDphiStar's type of events! Please provide same or mixed."; } - if (std::pow(dphiAvg, 2) / std::pow(cutDeltaPhiStarMaxVector[ktbinval], 2) + std::pow(deta, 2) / std::pow(cutDeltaEtaMaxVector[ktbinval], 2) > 1.) { - if (ChosenEventType == femto_universe_container::EventType::same) { - histdetadphisameafterkT[ktbinval]->Fill(deta, dphiAvg); - } else if (ChosenEventType == femto_universe_container::EventType::mixed) { - histdetadphimixedafterkT[ktbinval]->Fill(deta, dphiAvg); - } else { - LOG(fatal) << "FemtoUniverseDetaDphiStar: passed arguments don't agree with FemtoUniverseDetaDphiStar's type of events! Please provide same or mixed."; - } - } - if (CircCut && (std::pow(dphiAvg, 2) / std::pow(cutDeltaPhiStarMaxVector[ktbinval], 2) + std::pow(deta, 2) / std::pow(cutDeltaEtaMaxVector[ktbinval], 2) < 1.)) { return true; } else if (!CircCut && (dphiAvg > cutDeltaPhiStarMinVector[ktbinval]) && (dphiAvg < cutDeltaPhiStarMaxVector[ktbinval]) && (deta > cutDeltaEtaMinVector[ktbinval]) && (deta < cutDeltaEtaMaxVector[ktbinval])) { diff --git a/PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx b/PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx index e79c214a8cf..58987f15aca 100644 --- a/PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx +++ b/PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx @@ -526,7 +526,10 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended { } float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2); - if (kT > ConfKtKstarBins.value.back() || kT < ConfKtKstarBins.value.front()) + float lastElement = ConfKtKstarBins.value.back(); + float firstRealElement = ConfKtKstarBins.value[1]; + + if (kT < firstRealElement || kT > lastElement) continue; if (twotracksconfigs.ConfIsCPR.value) { @@ -565,7 +568,10 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended { } float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2); - if (kT > ConfKtKstarBins.value.back() || kT < ConfKtKstarBins.value.front()) + float lastElement = ConfKtKstarBins.value.back(); + float firstRealElement = ConfKtKstarBins.value[1]; + + if (kT < firstRealElement || kT > lastElement) continue; if (twotracksconfigs.ConfIsCPR.value) { @@ -944,7 +950,10 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended { } float kT = FemtoUniverseMath::getkT(p1, mass1, p2, mass2); - if (kT > ConfKtKstarBins.value.back() || kT < ConfKtKstarBins.value.front()) + float lastElement = ConfKtKstarBins.value.back(); + float firstRealElement = ConfKtKstarBins.value[1]; + + if (kT < firstRealElement || kT > lastElement) continue; if (twotracksconfigs.ConfIsCPR.value) {