Add configurable for TOFPID cut in hyperkinkRecoTask

Author: wang-yuanzhe

PWGLF/TableProducer/Nuspex/hyperkinkRecoTask.cxx

@@ -406,7 +406,8 @@ struct HyperkinkRecoTask {
   // Configurable for event selection
   Configurable<bool> doEventCut{"doEventCut", true, "Apply event selection"};
   Configurable<float> maxZVertex{"maxZVertex", 10.0f, "Accepted z-vertex range (cm)"};
-  Configurable<float> cutNSigmaDaug{"cutNSigmaDaug", 5, "TPC NSigmaTPC cut for daughter tracks"};
+  Configurable<float> cutTPCNSigmaDaug{"cutTPCNSigmaDaug", 5, "TPC NSigma cut for daughter tracks"};
+  Configurable<float> cutTOFNSigmaDaug{"cutTOFNSigmaDaug", 1000, "TOF NSigma cut for daughter tracks"};
 
   // CCDB options
   Configurable<double> inputBz{"inputBz", -999, "bz field, -999 is automatic"};
@@ -651,7 +652,7 @@ struct HyperkinkRecoTask {
 
       auto daugTrack = kinkCand.trackDaug_as<FullTracksExtIU>();
       float tpcNSigmaDaug = getTPCNSigma(daugTrack, pidTypeDaug);
-      if (std::abs(tpcNSigmaDaug) > cutNSigmaDaug) {
+      if (std::abs(tpcNSigmaDaug) > cutTPCNSigmaDaug) {
         continue;
       }
       float invMass = RecoDecay::m(std::array{std::array{kinkCand.pxDaug(), kinkCand.pyDaug(), kinkCand.pzDaug()}, std::array{kinkCand.pxDaugNeut(), kinkCand.pyDaugNeut(), kinkCand.pzDaugNeut()}}, std::array{massChargedDaug, massNeutralDaug});
@@ -669,6 +670,9 @@ struct HyperkinkRecoTask {
         auto originalDaugCol = daugTrack.collision_as<CollisionsFull>();
         nSigmaTOF = getTOFNSigma<false>(mRespParamsV3, daugTrack, collision, originalDaugCol);
       }
+      if (std::abs(nSigmaTOF) > cutTOFNSigmaDaug) {
+        continue;
+      }
       hypkinkCand.nSigmaTOFDaug = nSigmaTOF;
 
       o2::dataformats::VertexBase primaryVtx = {{collision.posX(), collision.posY(), collision.posZ()}, {collision.covXX(), collision.covXY(), collision.covYY(), collision.covXZ(), collision.covYZ(), collision.covZZ()}};
@@ -756,7 +760,7 @@ struct HyperkinkRecoTask {
       }
 
       float tpcNSigmaDaug = getTPCNSigma(daugTrack, pidTypeDaug);
-      if (std::abs(tpcNSigmaDaug) > cutNSigmaDaug) {
+      if (std::abs(tpcNSigmaDaug) > cutTPCNSigmaDaug) {
         continue;
       }
       float invMass = RecoDecay::m(std::array{std::array{kinkCand.pxDaug(), kinkCand.pyDaug(), kinkCand.pzDaug()}, std::array{kinkCand.pxDaugNeut(), kinkCand.pyDaugNeut(), kinkCand.pzDaugNeut()}}, std::array{massChargedDaug, massNeutralDaug});
@@ -784,6 +788,9 @@ struct HyperkinkRecoTask {
           registry.fill(HIST("hDaugNewTOFNSigma_WrongCol"), nSigmaTOF);
         }
       }
+      if (std::abs(nSigmaTOF) > cutTOFNSigmaDaug) {
+        continue;
+      }
       hypkinkCand.nSigmaTOFDaug = nSigmaTOF;
 
       std::array<float, 3> posDecVtx = {kinkCand.xDecVtx() + collision.posX(), kinkCand.yDecVtx() + collision.posY(), kinkCand.zDecVtx() + collision.posZ()};