Skip to content

Commit 3c9cd5d

Browse files
authored
Update FemtoUniverseDetaDphiStar.h
1 parent 1144180 commit 3c9cd5d

File tree

1 file changed

+7
-6
lines changed

1 file changed

+7
-6
lines changed

PWGCF/FemtoUniverse/Core/FemtoUniverseDetaDphiStar.h

Lines changed: 7 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -378,13 +378,12 @@ class FemtoUniverseDetaDphiStar
378378
// End: Get the charge from cutcontainer using masks
379379
float pt = part.pt();
380380
for (size_t i = 0; i < 9; i++) {
381-
auto arg = 0.3 * charge * magfield * tmpRadiiTPC[i] * 0.01 / (2. * pt);
382-
// for very low pT particles, this value goes outside of range -1 to 1 at
383-
// at large tpc radius; asin fails
381+
double arg = 0.3 * charge * magfield * tmpRadiiTPC[i] * 0.01 / (2. * pt);
384382
if (abs(arg) < 1.0) {
385383
tmpVec.push_back(phi0 - std::asin(arg));
386-
}
384+
}
387385
}
386+
}
388387

389388
/// Calculate average phi
390389
template <typename T1, typename T2>
@@ -435,10 +434,12 @@ class FemtoUniverseDetaDphiStar
435434
// double deltaphiconstAF = 0.15;
436435
double afsi0b = deltaphiconstFD * magfield * charge1 * ChosenRadii / part1.pt();
437436
double afsi1b = deltaphiconstFD * magfield * charge2 * ChosenRadii / part2.pt();
437+
double dphis = 0.0;
438438

439-
if ((abs(afsi0b) < 1.0) && (abs(afsi1b) < 1.0)) {
440-
double dphis = part2.phi() - part1.phi() + TMath::ASin(afsi1b) - TMath::ASin(afsi0b);
439+
if (abs(afsi0b) < 1.0 && abs(afsi0b) < 1.0) {
440+
dphis = part2.phi() - part1.phi() + TMath::ASin(afsi1b) - TMath::ASin(afsi0b);
441441
}
442+
return dphis;
442443
}
443444
};
444445

0 commit comments

Comments
 (0)