Update HFFilterHelpers.h

Author: zhangbiao-phy

EventFiltering/PWGHF/HFFilterHelpers.h

@@ -285,7 +285,7 @@ constexpr float cutsNsigma[4][7] = {
   {3., 3., 3., 5., 3., 3., 5.},            // TPC proton from Lc, pi/K from D0, K from 3-prong, femto selected proton, pi/K from Xic/Omegac, K from Xic*->SigmaC-Kaon, femto selected deuteron
   {3., 3., 3., 2.5, 3., 3., 5.},           // TOF proton from Lc, pi/K from D0, K from 3-prong, femto selected proton, pi/K from Xic/Omegac, K from Xic*->SigmaC-Kaon, femto selected deuteron
   {999., 999., 999., 2.5, 999., 999., 5.}, // Sum in quadrature of TPC and TOF (used only for femto selected proton and deuteron for pT < 4 GeV/c)
-  {999., 999., 999., 999., 999., 999., 4.} // used only for femto selected deuteron for less than pt threshold
+  {999., 999., 999., 999., 999., 999., -3.} // ITS used only for femto selected deuteron for less than pt threshold
 };
 static const std::vector<std::string> labelsColumnsNsigma = {"PrFromLc", "PiKaFromDZero", "KaFrom3Prong", "PrForFemto", "PiKaFromCharmBaryon", "SoftKaonFromXicResoToSigmaC", "DeForFemto"};
 static const std::vector<std::string> labelsRowsNsigma = {"TPC", "TOF", "Comb", "ITS"};
@@ -636,8 +636,8 @@ class HfFilterHelper
   float mPtMinSigmaC2520Zero{0.f};                                                // pt min SigmaC(2520)0 candidate
   float mPtMinSigmaCPlusPlus{0.f};                                                // pt min SigmaC++ candidate
   float mPtMinSigmaC2520PlusPlus{0.f};                                            // pt min SigmaC(2520)++ candidate
-  std::array<float, 4> mNSigmaPrCutsForFemto{3., 3., 3., 4.};                     // cut values for Nsigma TPC, TOF, combined, ITS for femto protons
-  std::array<float, 4> mNSigmaDeCutsForFemto{3., 3., 3., 4.};                     // cut values for Nsigma TPC, TOF, combined, ITS for femto deuterons
+  std::array<float, 4> mNSigmaPrCutsForFemto{3., 3., 3., -3.};                    // cut values for Nsigma TPC, TOF, combined, ITS for femto protons
+  std::array<float, 4> mNSigmaDeCutsForFemto{3., 3., 3., -3.};                    // cut values for Nsigma TPC, TOF, combined, ITS for femto deuterons
   float mNSigmaTpcPrCutForCharmBaryons{3.};                                       // maximum Nsigma TPC for protons in Lc and Xic decays
   float mNSigmaTofPrCutForCharmBaryons{3.};                                       // maximum Nsigma TOF for protons in Lc and Xic decays
   float mNSigmaTpcKaCutFor3Prongs{3.};                                            // maximum Nsigma TPC for kaons in 3-prong decays
@@ -876,7 +876,7 @@ inline bool HfFilterHelper::isSelectedTrack4Femto(const T1& track, const T2& tra
   if (trackSpecies == kDeuteronForFemto) {
     // Apply different PID strategy in different pt range
     if (pt <= ptThresholdPidStrategy) {
-      if (std::fabs(NSigmaTPC) > nSigmaCuts[0] && std::fabs(NSigmaITS) > nSigmaCuts[3]) { // Use TPC and ITS below the threshold
+      if (std::fabs(NSigmaTPC) > nSigmaCuts[0] && NSigmaITS > nSigmaCuts[3]) { // Use TPC and ITS below the threshold, NSigmaITS for deuteron with a lower limit
         return false;
       }
     } else {