Skip to content

Commit cddbd8c

Browse files
authored
[PWGCF] FemtoUniverse: Fix bug in accessing deltaeta and deltaphi ranges (#13737)
1 parent 207dafe commit cddbd8c

File tree

1 file changed

+7
-2
lines changed

1 file changed

+7
-2
lines changed

PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h

Lines changed: 7 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -668,6 +668,11 @@ class FemtoUniverseDetaDphiStar
668668

669669
auto deta = part1.eta() - part2.eta();
670670
auto dphiAvg = averagePhiStar(part1, part2, 0);
671+
auto DeltaPhiStarMax = static_cast<float>(cutDeltaPhiStarMaxVector[ktbinval]);
672+
auto DeltaPhiStarMin = static_cast<float>(cutDeltaPhiStarMinVector[ktbinval]);
673+
auto DeltaEtaMax = static_cast<float>(cutDeltaEtaMaxVector[ktbinval]);
674+
auto DeltaEtaMin = static_cast<float>(cutDeltaEtaMinVector[ktbinval]);
675+
671676
if (ChosenEventType == femto_universe_container::EventType::same) {
672677
histdetadphisamebeforekT[ktbinval]->Fill(deta, dphiAvg);
673678
} else if (ChosenEventType == femto_universe_container::EventType::mixed) {
@@ -676,9 +681,9 @@ class FemtoUniverseDetaDphiStar
676681
LOG(fatal) << "FemtoUniverseDetaDphiStar: passed arguments don't agree with FemtoUniverseDetaDphiStar's type of events! Please provide same or mixed.";
677682
}
678683

679-
if (CircCut && (std::pow(dphiAvg, 2) / std::pow(cutDeltaPhiStarMaxVector[ktbinval], 2) + std::pow(deta, 2) / std::pow(cutDeltaEtaMaxVector[ktbinval], 2) < 1.)) {
684+
if (CircCut && (std::pow(dphiAvg, 2) / std::pow(DeltaPhiStarMax, 2) + std::pow(deta, 2) / std::pow(DeltaEtaMax, 2) < 1.)) {
680685
return true;
681-
} else if (!CircCut && (dphiAvg > cutDeltaPhiStarMinVector[ktbinval]) && (dphiAvg < cutDeltaPhiStarMaxVector[ktbinval]) && (deta > cutDeltaEtaMinVector[ktbinval]) && (deta < cutDeltaEtaMaxVector[ktbinval])) {
686+
} else if (!CircCut && (dphiAvg > DeltaPhiStarMin) && (dphiAvg < DeltaPhiStarMax) && (deta > DeltaEtaMin) && (deta < DeltaEtaMax)) {
682687
return true;
683688
} else {
684689
if (ChosenEventType == femto_universe_container::EventType::same) {

0 commit comments

Comments
 (0)