diff --git a/PWGCF/Femto/FemtoNuclei/Tasks/piDeFemtoSystematics.cxx b/PWGCF/Femto/FemtoNuclei/Tasks/piDeFemtoSystematics.cxx index c12af02d7d2..515dbde5ba7 100644 --- a/PWGCF/Femto/FemtoNuclei/Tasks/piDeFemtoSystematics.cxx +++ b/PWGCF/Femto/FemtoNuclei/Tasks/piDeFemtoSystematics.cxx @@ -144,7 +144,9 @@ struct PiDeFemtoSystematics { Configurable sharedTPCFractionMax{"sharedTPCFractionMax", 1.f, "Maximum deuteron shared TPC-cluster fraction"}; Configurable tpcInnerParamMin{"tpcInnerParamMin", 0.f, "Minimum deuteron TPC inner parameter"}; Configurable tofMomentumMin{"tofMomentumMin", 1.2f, "TPC inner parameter above which combined TPC+TOF PID is required"}; - Configurable nsigmaMax{"nsigmaMax", 2.5f, "Common deuteron TPC, ITS and combined PID threshold"}; + Configurable combinedNsigmaMax{"combinedNsigmaMax", 2.5f, "Maximum deuteron combined TPC+TOF n-sigma"}; + Configurable tpcNsigmaMax{"tpcNsigmaMax", 2.5f, "Maximum absolute deuteron TPC n-sigma"}; + Configurable itsNsigmaMax{"itsNsigmaMax", 2.5f, "Maximum absolute deuteron ITS n-sigma"}; Configurable requireIndividualNsigma{"requireIndividualNsigma", false, "Also apply individual TPC and TOF cuts in the combined PID branch"}; Configurable dcaXYOffset{"dcaXYOffset", 0.004f, "Deuteron DCAxy cut offset"}; Configurable dcaXYPtCoefficient{"dcaXYPtCoefficient", 0.013f, "Deuteron DCAxy inverse-pT coefficient"}; @@ -358,7 +360,9 @@ struct PiDeFemtoSystematics { {"deuteronSharedTPCFractionMax", deuteronCuts.sharedTPCFractionMax.value}, {"deuteronTPCInnerParamMin", deuteronCuts.tpcInnerParamMin.value}, {"deuteronTOFMomentumMin", deuteronCuts.tofMomentumMin.value}, - {"deuteronNsigmaMax", deuteronCuts.nsigmaMax.value}, + {"deuteronCombNsigmaMax", deuteronCuts.combinedNsigmaMax.value}, + {"deuteronTPCNsigmaMax", deuteronCuts.tpcNsigmaMax.value}, + {"deuteronITSNsigmaMax", deuteronCuts.itsNsigmaMax.value}, {"requireIndividualNsigma", deuteronCuts.requireIndividualNsigma.value}, {"deuteronDCAxyOffset", deuteronCuts.dcaXYOffset.value}, {"deuteronDCAxyPtCoefficient", deuteronCuts.dcaXYPtCoefficient.value}, @@ -536,22 +540,22 @@ struct PiDeFemtoSystematics { if (track.hasTOF() && tpcInnerParam > deuteronCuts.tofMomentumMin.value) { const float tofNsigma = track.tofNSigmaDe(); - if (std::hypot(tpcNsigma, tofNsigma) > deuteronCuts.nsigmaMax.value) { + if (std::hypot(tpcNsigma, tofNsigma) > deuteronCuts.combinedNsigmaMax.value) { return false; } return !deuteronCuts.requireIndividualNsigma.value || - (std::abs(tpcNsigma) <= deuteronCuts.nsigmaMax.value && - std::abs(tofNsigma) <= deuteronCuts.nsigmaMax.value); + (std::abs(tpcNsigma) <= deuteronCuts.combinedNsigmaMax.value && + std::abs(tofNsigma) <= deuteronCuts.combinedNsigmaMax.value); } if (tpcInnerParam <= deuteronCuts.tofMomentumMin.value) { - if (std::abs(tpcNsigma) > deuteronCuts.nsigmaMax.value) { + if (std::abs(tpcNsigma) > deuteronCuts.tpcNsigmaMax.value) { return false; } o2::aod::ITSResponse itsResponse; const float itsNsigma = itsResponse.nSigmaITS( track.itsClusterSizes(), track.p(), track.eta()); - return std::abs(itsNsigma) <= deuteronCuts.nsigmaMax.value; + return std::abs(itsNsigma) <= deuteronCuts.itsNsigmaMax.value; } return false; }