From 399eff9d540aa0e602ea337082eadd45271e5a46 Mon Sep 17 00:00:00 2001
From: "se.gorbunov" <se.gorbunov@gsi.de>
Date: Wed, 31 Aug 2022 15:58:33 +0000
Subject: [PATCH] L1: bugfix in time extrapolation with Runge-Kutta method
---
reco/L1/L1Algo/L1TrackParFit.cxx | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/reco/L1/L1Algo/L1TrackParFit.cxx b/reco/L1/L1Algo/L1TrackParFit.cxx
index 7aaacce283..958d197e26 100644
--- a/reco/L1/L1Algo/L1TrackParFit.cxx
+++ b/reco/L1/L1Algo/L1TrackParFit.cxx
@@ -410,9 +410,9 @@ void L1TrackParFit::Extrapolate // extrapolates track parameters and returns ja
fvec p2 = 1.f / (qp0 * qp0);
fvec m2 = fMass2;
fvec v = 29.9792458f * sqrt(p2 / (m2 + p2));
- At[step] = sqrt(1.f + tx * tx + ty * ty) / v;
- At_tx[step] = tx / sqrt(1.f + tx * tx + ty * ty) / v;
- At_ty[step] = ty / sqrt(1.f + tx * tx + ty * ty) / v;
+ At[step] = h * sqrt(1.f + tx * tx + ty * ty) / v;
+ At_tx[step] = h * tx / sqrt(1.f + tx * tx + ty * ty) / v;
+ At_ty[step] = h * ty / sqrt(1.f + tx * tx + ty * ty) / v;
// cout<<Ax[step]<<" Ax[step] "<<Ay[step]<<" ay "<<At[step]<<" At[step] "<<qp0<<" qp0 "<<h<<" h"<<endl;
@@ -421,7 +421,7 @@ void L1TrackParFit::Extrapolate // extrapolates track parameters and returns ja
k[step4 + 1] = ty * h;
k[step4 + 2] = Ax[step] * qp0;
k[step4 + 3] = Ay[step] * qp0;
- k[step4 + 4] = At[step] * h;
+ k[step4 + 4] = At[step];
} // end of Runge-Kutta steps
fvec initialised = ZERO;
--
GitLab