Skip to content

Commit 8da160d

Browse files
authored
[PWGLF] Fix TOF PID selection and decay mode check for hypertriton kink analysis (#12500)
1 parent 0f3b1ea commit 8da160d

File tree

1 file changed

+5
-4
lines changed

1 file changed

+5
-4
lines changed

PWGLF/TableProducer/Nuspex/hyperkinkRecoTask.cxx

Lines changed: 5 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -142,9 +142,9 @@ struct H3LDecay {
142142

143143
if ((haveTriton && havePion0) || (haveAntiTriton && havePion0)) {
144144
return H3LDecay::k2bodyNeutral;
145-
} else if ((haveHelium3 && havePion0) || (haveAntiHelium3 && havePion0)) {
145+
} else if ((haveHelium3 && havePionMinus) || (haveAntiHelium3 && havePionPlus)) {
146146
return H3LDecay::k2bodyCharged;
147-
} else if ((haveDeuteron && haveProton && havePionPlus) || (haveAntiDeuteron && haveAntiProton && havePionMinus)) {
147+
} else if ((haveDeuteron && haveProton && havePionMinus) || (haveAntiDeuteron && haveAntiProton && havePionPlus)) {
148148
return H3LDecay::k3bodyCharged;
149149
} else {
150150
return kNChannel;
@@ -408,6 +408,7 @@ struct HyperkinkRecoTask {
408408
Configurable<float> maxZVertex{"maxZVertex", 10.0f, "Accepted z-vertex range (cm)"};
409409
Configurable<float> cutTPCNSigmaDaug{"cutTPCNSigmaDaug", 5, "TPC NSigma cut for daughter tracks"};
410410
Configurable<float> cutTOFNSigmaDaug{"cutTOFNSigmaDaug", 1000, "TOF NSigma cut for daughter tracks"};
411+
Configurable<bool> askTOFForDaug{"askTOFForDaug", false, "If true, ask for TOF signal"};
411412

412413
// CCDB options
413414
Configurable<double> inputBz{"inputBz", -999, "bz field, -999 is automatic"};
@@ -671,7 +672,7 @@ struct HyperkinkRecoTask {
671672
auto originalDaugCol = daugTrack.collision_as<CollisionsFull>();
672673
nSigmaTOF = getTOFNSigma<false>(mRespParamsV3, daugTrack, collision, originalDaugCol);
673674
}
674-
if (std::abs(nSigmaTOF) > cutTOFNSigmaDaug) {
675+
if ((daugTrack.hasTOF() || askTOFForDaug) && std::abs(nSigmaTOF) > cutTOFNSigmaDaug) {
675676
continue;
676677
}
677678

@@ -793,7 +794,7 @@ struct HyperkinkRecoTask {
793794
registry.fill(HIST("hDaugNewTOFNSigma_WrongCol"), nSigmaTOF);
794795
}
795796
}
796-
if (std::abs(nSigmaTOF) > cutTOFNSigmaDaug) {
797+
if ((daugTrack.hasTOF() || askTOFForDaug) && std::abs(nSigmaTOF) > cutTOFNSigmaDaug) {
797798
continue;
798799
}
799800
registry.fill(HIST("hCandidateCounter"), 3);

0 commit comments

Comments
 (0)