Skip to content

Commit f81b945

Browse files
fix rapidity
1 parent 41a38b5 commit f81b945

File tree

1 file changed

+4
-3
lines changed

1 file changed

+4
-3
lines changed

PWGJE/Tasks/nucleiInJets.cxx

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1643,7 +1643,7 @@ struct nucleiInJets {
16431643

16441644
auto rapidityData = [&](float m2z) {
16451645
const float rap = trk.rapidity(m2z);
1646-
return rap > std::abs(cfgtrkMaxRap);
1646+
return rap < std::abs(cfgtrkMaxRap);
16471647
};
16481648

16491649
auto prRapidityWithinRange = rapidityData(o2::constants::physics::MassProton);
@@ -2314,8 +2314,9 @@ struct nucleiInJets {
23142314
if (std::fabs(track.dcaXY()) > dcaXYPtDepCut(track.pt()) && useDcaxyPtDepCut)
23152315
continue;
23162316

2317-
auto mass = TDatabasePDG::Instance()->GetParticle(abs(mcTrack.pdgCode()))->Mass();
2318-
auto rapidity = RecoDecay::y(std::array{track.px(), track.py(), track.pz()}, mass);
2317+
// auto mass = TDatabasePDG::Instance()->GetParticle(abs(mcTrack.pdgCode()))->Mass();
2318+
// auto rapidity = RecoDecay::y(std::array{track.px(), track.py(), track.pz()}, mass);
2319+
auto rapidity = mcTrack.y();
23192320

23202321
if (std::abs(rapidity) > cfgtrkMaxRap)
23212322
continue;

0 commit comments

Comments
 (0)