Skip to content

Commit a3589d5

Browse files
authored
[PWGCF] fix configurable (#17422)
1 parent 712ca39 commit a3589d5

1 file changed

Lines changed: 11 additions & 7 deletions

File tree

PWGCF/Femto/FemtoNuclei/Tasks/piDeFemtoSystematics.cxx

Lines changed: 11 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -144,7 +144,9 @@ struct PiDeFemtoSystematics {
144144
Configurable<float> sharedTPCFractionMax{"sharedTPCFractionMax", 1.f, "Maximum deuteron shared TPC-cluster fraction"};
145145
Configurable<float> tpcInnerParamMin{"tpcInnerParamMin", 0.f, "Minimum deuteron TPC inner parameter"};
146146
Configurable<float> tofMomentumMin{"tofMomentumMin", 1.2f, "TPC inner parameter above which combined TPC+TOF PID is required"};
147-
Configurable<float> nsigmaMax{"nsigmaMax", 2.5f, "Common deuteron TPC, ITS and combined PID threshold"};
147+
Configurable<float> combinedNsigmaMax{"combinedNsigmaMax", 2.5f, "Maximum deuteron combined TPC+TOF n-sigma"};
148+
Configurable<float> tpcNsigmaMax{"tpcNsigmaMax", 2.5f, "Maximum absolute deuteron TPC n-sigma"};
149+
Configurable<float> itsNsigmaMax{"itsNsigmaMax", 2.5f, "Maximum absolute deuteron ITS n-sigma"};
148150
Configurable<bool> requireIndividualNsigma{"requireIndividualNsigma", false, "Also apply individual TPC and TOF cuts in the combined PID branch"};
149151
Configurable<float> dcaXYOffset{"dcaXYOffset", 0.004f, "Deuteron DCAxy cut offset"};
150152
Configurable<float> dcaXYPtCoefficient{"dcaXYPtCoefficient", 0.013f, "Deuteron DCAxy inverse-pT coefficient"};
@@ -358,7 +360,9 @@ struct PiDeFemtoSystematics {
358360
{"deuteronSharedTPCFractionMax", deuteronCuts.sharedTPCFractionMax.value},
359361
{"deuteronTPCInnerParamMin", deuteronCuts.tpcInnerParamMin.value},
360362
{"deuteronTOFMomentumMin", deuteronCuts.tofMomentumMin.value},
361-
{"deuteronNsigmaMax", deuteronCuts.nsigmaMax.value},
363+
{"deuteronCombNsigmaMax", deuteronCuts.combinedNsigmaMax.value},
364+
{"deuteronTPCNsigmaMax", deuteronCuts.tpcNsigmaMax.value},
365+
{"deuteronITSNsigmaMax", deuteronCuts.itsNsigmaMax.value},
362366
{"requireIndividualNsigma", deuteronCuts.requireIndividualNsigma.value},
363367
{"deuteronDCAxyOffset", deuteronCuts.dcaXYOffset.value},
364368
{"deuteronDCAxyPtCoefficient", deuteronCuts.dcaXYPtCoefficient.value},
@@ -536,22 +540,22 @@ struct PiDeFemtoSystematics {
536540
if (track.hasTOF() &&
537541
tpcInnerParam > deuteronCuts.tofMomentumMin.value) {
538542
const float tofNsigma = track.tofNSigmaDe();
539-
if (std::hypot(tpcNsigma, tofNsigma) > deuteronCuts.nsigmaMax.value) {
543+
if (std::hypot(tpcNsigma, tofNsigma) > deuteronCuts.combinedNsigmaMax.value) {
540544
return false;
541545
}
542546
return !deuteronCuts.requireIndividualNsigma.value ||
543-
(std::abs(tpcNsigma) <= deuteronCuts.nsigmaMax.value &&
544-
std::abs(tofNsigma) <= deuteronCuts.nsigmaMax.value);
547+
(std::abs(tpcNsigma) <= deuteronCuts.combinedNsigmaMax.value &&
548+
std::abs(tofNsigma) <= deuteronCuts.combinedNsigmaMax.value);
545549
}
546550
if (tpcInnerParam <= deuteronCuts.tofMomentumMin.value) {
547-
if (std::abs(tpcNsigma) > deuteronCuts.nsigmaMax.value) {
551+
if (std::abs(tpcNsigma) > deuteronCuts.tpcNsigmaMax.value) {
548552
return false;
549553
}
550554
o2::aod::ITSResponse itsResponse;
551555
const float itsNsigma =
552556
itsResponse.nSigmaITS<o2::track::PID::Deuteron>(
553557
track.itsClusterSizes(), track.p(), track.eta());
554-
return std::abs(itsNsigma) <= deuteronCuts.nsigmaMax.value;
558+
return std::abs(itsNsigma) <= deuteronCuts.itsNsigmaMax.value;
555559
}
556560
return false;
557561
}

0 commit comments

Comments
 (0)