From ee7ffcac05ff3b98015e3ef01e404057c2d57da2 Mon Sep 17 00:00:00 2001
From: "se.gorbunov" <se.gorbunov@gsi.de>
Date: Wed, 5 Apr 2023 14:57:33 +0000
Subject: [PATCH] L1: make tracking independent of the MVD time error
---
reco/L1/CbmL1ReadEvent.cxx | 4 ++++
reco/L1/L1Algo/L1CaTrackFinder.cxx | 7 ++++++-
reco/L1/L1Algo/L1CaTrackFinderSlice.cxx | 2 +-
3 files changed, 11 insertions(+), 2 deletions(-)
diff --git a/reco/L1/CbmL1ReadEvent.cxx b/reco/L1/CbmL1ReadEvent.cxx
index fdd5f17244..04ed4ea15d 100644
--- a/reco/L1/CbmL1ReadEvent.cxx
+++ b/reco/L1/CbmL1ReadEvent.cxx
@@ -694,6 +694,10 @@ void CbmL1::ReadEvent(CbmEvent* event)
const L1Station& st = fpAlgo->GetParameters()->GetStation(th.iStation);
th.u = th.x * st.frontInfo.cos_phi[0] + th.y * st.frontInfo.sin_phi[0];
th.v = th.x * st.backInfo.cos_phi[0] + th.y * st.backInfo.sin_phi[0];
+
+ // Get time
+ th.time = h->GetTime(); // currently ignored by the tracking
+ th.dt = h->GetTimeError(); // currently ignored by the tracking
}
th.Det = 0;
th.iMC = fPerformance ? MatchHitWithMc<L1DetectorID::kMvd>(hitIndex) : -1;
diff --git a/reco/L1/L1Algo/L1CaTrackFinder.cxx b/reco/L1/L1Algo/L1CaTrackFinder.cxx
index eba4866229..113a0cc7e4 100644
--- a/reco/L1/L1Algo/L1CaTrackFinder.cxx
+++ b/reco/L1/L1Algo/L1CaTrackFinder.cxx
@@ -66,7 +66,7 @@ void L1Algo::CaTrackFinder()
for (int iStream = 0; iStream < nDataStreams; ++iStream) {
- fscal maxTimeBeforeHit = std::numeric_limits<fscal>::min();
+ fscal maxTimeBeforeHit = std::numeric_limits<fscal>::lowest();
int nStreamHits = fInputData.GetStreamNhits(iStream);
@@ -92,6 +92,11 @@ void L1Algo::CaTrackFinder()
info.fEventTimeMin = h.t - dt - timeOfFlightMax;
info.fEventTimeMax = h.t + dt - timeOfFlightMin;
+ if (!st.timeInfo) {
+ info.fEventTimeMin = std::numeric_limits<fscal>::lowest();
+ info.fEventTimeMax = std::numeric_limits<fscal>::max();
+ }
+
if (maxTimeBeforeHit < info.fEventTimeMax) { maxTimeBeforeHit = info.fEventTimeMax; }
info.fMaxTimeBeforeHit = maxTimeBeforeHit;
diff --git a/reco/L1/L1Algo/L1CaTrackFinderSlice.cxx b/reco/L1/L1Algo/L1CaTrackFinderSlice.cxx
index 360cb31b1f..9f211669b3 100644
--- a/reco/L1/L1Algo/L1CaTrackFinderSlice.cxx
+++ b/reco/L1/L1Algo/L1CaTrackFinderSlice.cxx
@@ -271,7 +271,7 @@ inline void L1Algo::findSingletsStep1( /// input 1st stage of singlet search
fvec txErr2 = fMaxSlopePV * fMaxSlopePV / fvec(9.);
fvec qpErr2 = fMaxInvMom * fMaxInvMom / fvec(9.);
- T.ResetErrors(1., 1., txErr2, txErr2, qpErr2, timeEr2, 1.e2);
+ T.ResetErrors(1., 1., txErr2, txErr2, qpErr2, (stal.timeInfo ? timeEr2 : 1.e6), 1.e2);
T.chi2 = 0.;
T.NDF = (fpCurrentIteration->GetPrimaryFlag()) ? fvec(2.) : fvec(0.);
--
GitLab