@@ -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 > CombnsigmaMax{" CombnsigmaMax" , 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.CombnsigmaMax .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.CombnsigmaMax .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.CombnsigmaMax .value &&
548+ std::abs (tofNsigma) <= deuteronCuts.CombnsigmaMax .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