diff --git a/PWGCF/TwoParticleCorrelations/Tasks/nucleibalance.cxx b/PWGCF/TwoParticleCorrelations/Tasks/nucleibalance.cxx index ded525fd309..ee466081ddb 100644 --- a/PWGCF/TwoParticleCorrelations/Tasks/nucleibalance.cxx +++ b/PWGCF/TwoParticleCorrelations/Tasks/nucleibalance.cxx @@ -1710,6 +1710,7 @@ struct Lambdastarproxy { Configurable lstarCutNsigmaTOFKaon{"lstarCutNsigmaTOFKaon", float{NsigmaTOFDefault}, "|nSigma^{TOF}_{K}| cut"}; Configurable lstarCutNsigmaTPCDe{"lstarCutNsigmaTPCDe", float{NsigmaTPCDefault}, "|nSigma^{TPC}_{d}| cut"}; Configurable lstarCutNsigmaTOFDe{"lstarCutNsigmaTOFDe", float{NsigmaTOFDefault}, "|nSigma^{TOF}_{d}| cut"}; + Configurable lstarEnableTOFNsigmaCutDe{"lstarEnableTOFNsigmaCutDe", 0, "Enable deuteron-only TOF nSigma cut in PID strategy 2"}; // Optional deuteron-only TOF auxiliary selections. // Defaults are OFF, so strategy 2 remains TPC nSigma only for deuterons. Configurable lstarEnableBetaCutDe{"lstarEnableBetaCutDe", 0, "Enable deuteron-only TOF beta cut using beta() > lstarBetaCutDe"}; @@ -2453,7 +2454,16 @@ struct Lambdastarproxy { // For deuterons, follow the default idea of the official nuclei task: // use TPC nσ as the main hard PID selection. if (isDeuteron) { - return std::abs(nsTPC) < tpcCut; + if (std::abs(nsTPC) >= tpcCut) { + return false; + } + if (lstarEnableTOFNsigmaCutDe.value != 0) { + if (!hasTof) { + return false; + } + return std::abs(nsTOF) < tofCut; + } + return true; } // For kaons/protons, use the Lambda(1520)-like pT-ref circular TPC+TOF logic.