Skip to content

Commit dcb1598

Browse files
authored
[PWGLF] Fix DCA determination for tracked cascades (#13934)
1 parent bcb48a4 commit dcb1598

File tree

1 file changed

+14
-8
lines changed

1 file changed

+14
-8
lines changed

PWGLF/Tasks/Strangeness/strangeCascTrack.cxx

Lines changed: 14 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -283,18 +283,24 @@ struct StrangeCascTrack {
283283
return false;
284284
return true;
285285
}
286-
// checks general selection criteria for cascades
286+
// checks general selection criteria for cascades - from std casc
287287
template <typename TEvent, typename TCascade>
288-
bool isValidCasc(TEvent collision, TCascade cascade, TString particle)
288+
bool isValidCascGen(TEvent collision, TCascade cascade)
289289
{
290-
if (cascade.dcaXYCascToPV() > selCuts.cutDCAtoPVxy)
291-
return false;
292-
if (cascade.dcaZCascToPV() > selCuts.cutDCAtoPVz)
293-
return false;
294290
if (cascade.v0cosPA(collision.posX(), collision.posY(), collision.posZ()) < selCuts.cutV0CosPA)
295291
return false;
296292
if (cascade.bachBaryonCosPA() < selCuts.cutBachCosPA)
297293
return false;
294+
return true;
295+
}
296+
// checks general selection criteria for cascades - from std or tracked casc
297+
template <typename TCascade>
298+
bool isValidCascSpec(TCascade cascade, TString particle)
299+
{
300+
if (cascade.dcaXYCascToPV() > selCuts.cutDCAtoPVxy)
301+
return false;
302+
if (cascade.dcaZCascToPV() > selCuts.cutDCAtoPVz)
303+
return false;
298304
ROOT::Math::PxPyPzMVector momentum;
299305
if (particle == "xi")
300306
momentum.SetCoordinates(cascade.px(), cascade.py(), cascade.pz(), o2::constants::physics::MassXiMinus);
@@ -603,11 +609,11 @@ struct StrangeCascTrack {
603609
}
604610
}
605611
if (doApplyGenCutsXi) {
606-
if (!isValidCasc(collision, stdCasc, "xi"))
612+
if (!isValidCascGen(collision, stdCasc) || !isValidCascSpec(cascade, "xi"))
607613
passedAllSelsXi = false;
608614
}
609615
if (doApplyGenCutsOmega) {
610-
if (!isValidCasc(collision, stdCasc, "omega"))
616+
if (!isValidCascGen(collision, stdCasc) || !isValidCascSpec(cascade, "omega"))
611617
passedAllSelsOmega = false;
612618
}
613619
// apply pt cuts

0 commit comments

Comments
 (0)