From 3f94c699d7e5101242f168732bb73e68fad29128 Mon Sep 17 00:00:00 2001
From: yuanzhe <yuanzhe.wang@cern.ch>
Date: Tue, 10 Jun 2025 18:04:22 +0200
Subject: [PATCH] "Add option to disable SV radius cut for kink"

---
 PWGLF/TableProducer/Common/kinkBuilder.cxx | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/PWGLF/TableProducer/Common/kinkBuilder.cxx b/PWGLF/TableProducer/Common/kinkBuilder.cxx
index d2b25be26a0..51be861741a 100644
--- a/PWGLF/TableProducer/Common/kinkBuilder.cxx
+++ b/PWGLF/TableProducer/Common/kinkBuilder.cxx
@@ -106,6 +106,7 @@ struct kinkBuilder {
   Configurable<float> etaMax{"etaMax", 1., "eta daughter"};
   Configurable<float> nTPCClusMinDaug{"nTPCClusMinDaug", 80, "daug NTPC clusters cut"};
   Configurable<bool> askTOFforDaug{"askTOFforDaug", false, "If true, ask for TOF signal"};
+  Configurable<bool> doSVRadiusCut{"doSVRadiusCut", true, "If true, apply the cut on the radius of the secondary vertex and tracksIU"};
 
   o2::vertexing::DCAFitterN<2> fitter;
   o2::base::MatLayerCylSet* lut = nullptr;
@@ -329,7 +330,7 @@ struct kinkBuilder {
 
       // cut on decay radius to 17 cm
       float decRad2 = kinkCand.decVtx[0] * kinkCand.decVtx[0] + kinkCand.decVtx[1] * kinkCand.decVtx[1];
-      if (decRad2 < LayerRadii[3] * LayerRadii[3]) {
+      if (doSVRadiusCut && decRad2 < LayerRadii[3] * LayerRadii[3]) {
         continue;
       }
 
@@ -348,11 +349,11 @@ struct kinkBuilder {
         }
       }
 
-      if (lastLayerMoth >= firstLayerDaug) {
+      if (doSVRadiusCut && lastLayerMoth >= firstLayerDaug) {
         continue;
       }
 
-      if (decRad2 < LayerRadii[lastLayerMoth] * LayerRadii[lastLayerMoth]) {
+      if (doSVRadiusCut && decRad2 < LayerRadii[lastLayerMoth] * LayerRadii[lastLayerMoth]) {
         continue;
       }