From 73e37f9ee86da4d7d742f4bf50b84637405d3d72 Mon Sep 17 00:00:00 2001
From: "se.gorbunov" <se.gorbunov@gsi.de>
Date: Sat, 11 Feb 2023 02:19:33 +0000
Subject: [PATCH] L1: remove a hardcoded cut for mcbm/global tracker
---
reco/L1/CbmL1.cxx | 11 ++++++++---
reco/L1/L1Algo/L1CATrackFinder.cxx | 4 ----
2 files changed, 8 insertions(+), 7 deletions(-)
diff --git a/reco/L1/CbmL1.cxx b/reco/L1/CbmL1.cxx
index d03df39f8c..68e7d8df88 100644
--- a/reco/L1/CbmL1.cxx
+++ b/reco/L1/CbmL1.cxx
@@ -732,6 +732,11 @@ InitStatus CbmL1::Init()
trackingIterAllPrimJump.SetExtendTracksFlag(false);
trackingIterAllSec.SetExtendTracksFlag(false);
+ trackingIterFastPrim.SetMaxSlopePV(9.5);
+ trackingIterAllPrim.SetMaxSlopePV(9.5);
+ trackingIterAllPrimJump.SetMaxSlopePV(9.5);
+ trackingIterAllSec.SetMaxSlopePV(9.5);
+
fInitManager.SetCAIterationsNumberCrosscheck(4);
// Initialize CA track finder iterations sequence
fInitManager.PushBackCAIteration(trackingIterFastPrim);
@@ -754,7 +759,7 @@ InitStatus CbmL1::Init()
it.SetPickGather(3.0f);
it.SetTripletLinkChi2(16.0);
it.SetMaxQp(1.0 / 0.1); //(1.0 / 0.5);
- it.SetMaxSlopePV(.1f);
+ it.SetMaxSlopePV(9.5);
it.SetMaxSlope(.5f);
it.SetMaxDZ(0.05);
it.SetTargetPosSigmaXY(1., 1.); //(1, 1);
@@ -774,7 +779,7 @@ InitStatus CbmL1::Init()
trd2dIter2.SetPickGather(3.0f);
trd2dIter2.SetTripletLinkChi2(16.);
trd2dIter2.SetMaxQp(1.0 / 0.05); //(1.0 / 0.5);
- trd2dIter2.SetMaxSlopePV(.5f);
+ trd2dIter2.SetMaxSlopePV(9.5f);
trd2dIter2.SetMaxSlope(.5f);
trd2dIter2.SetMaxDZ(0.05);
trd2dIter2.SetTargetPosSigmaXY(8 * 10, 6 * 10); //(1, 1);
@@ -791,7 +796,7 @@ InitStatus CbmL1::Init()
trd2dIter3.SetPickGather(3.0f);
trd2dIter3.SetTripletLinkChi2(200.);
trd2dIter3.SetMaxQp(1.0 / 0.1); //(1.0 / 0.5);
- trd2dIter3.SetMaxSlopePV(.01);
+ trd2dIter3.SetMaxSlopePV(9.5);
trd2dIter3.SetMaxSlope(.4); //.5f);
trd2dIter3.SetMaxDZ(0.05);
trd2dIter3.SetTargetPosSigmaXY(7 * 10, 6 * 10); //(1, 1);
diff --git a/reco/L1/L1Algo/L1CATrackFinder.cxx b/reco/L1/L1Algo/L1CATrackFinder.cxx
index 3f941d0a96..a35ab55006 100644
--- a/reco/L1/L1Algo/L1CATrackFinder.cxx
+++ b/reco/L1/L1Algo/L1CATrackFinder.cxx
@@ -281,10 +281,6 @@ inline void L1Algo::findSingletsStep1( /// input 1st stage of singlet search
T.C22 = T.C33 = fMaxSlopePV * fMaxSlopePV / fvec(9.);
- if (kGlobal == fTrackingMode || kMcbm == fTrackingMode) T.C22 = T.C33 = fvec(10.);
-
- // TODO: Why 9. and 10.?
-
T.C44 = fMaxInvMom / fvec(3.) * fMaxInvMom / fvec(3.);
T.C55 = timeEr2;
--
GitLab