Skip to content

Commit 3f94c69

Browse files
committed
"Add option to disable SV radius cut for kink"
1 parent 4c030d6 commit 3f94c69

File tree

1 file changed

+4
-3
lines changed

1 file changed

+4
-3
lines changed

PWGLF/TableProducer/Common/kinkBuilder.cxx

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -106,6 +106,7 @@ struct kinkBuilder {
106106
Configurable<float> etaMax{"etaMax", 1., "eta daughter"};
107107
Configurable<float> nTPCClusMinDaug{"nTPCClusMinDaug", 80, "daug NTPC clusters cut"};
108108
Configurable<bool> askTOFforDaug{"askTOFforDaug", false, "If true, ask for TOF signal"};
109+
Configurable<bool> doSVRadiusCut{"doSVRadiusCut", true, "If true, apply the cut on the radius of the secondary vertex and tracksIU"};
109110

110111
o2::vertexing::DCAFitterN<2> fitter;
111112
o2::base::MatLayerCylSet* lut = nullptr;
@@ -329,7 +330,7 @@ struct kinkBuilder {
329330

330331
// cut on decay radius to 17 cm
331332
float decRad2 = kinkCand.decVtx[0] * kinkCand.decVtx[0] + kinkCand.decVtx[1] * kinkCand.decVtx[1];
332-
if (decRad2 < LayerRadii[3] * LayerRadii[3]) {
333+
if (doSVRadiusCut && decRad2 < LayerRadii[3] * LayerRadii[3]) {
333334
continue;
334335
}
335336

@@ -348,11 +349,11 @@ struct kinkBuilder {
348349
}
349350
}
350351

351-
if (lastLayerMoth >= firstLayerDaug) {
352+
if (doSVRadiusCut && lastLayerMoth >= firstLayerDaug) {
352353
continue;
353354
}
354355

355-
if (decRad2 < LayerRadii[lastLayerMoth] * LayerRadii[lastLayerMoth]) {
356+
if (doSVRadiusCut && decRad2 < LayerRadii[lastLayerMoth] * LayerRadii[lastLayerMoth]) {
356357
continue;
357358
}
358359

0 commit comments

Comments
 (0)