From 604f776db6aea57944fcc37f9c012f5d0c95420d Mon Sep 17 00:00:00 2001
From: Felix Schlepper <felix.schlepper@cern.ch>
Date: Wed, 19 Aug 2026 16:13:26 +0200
Subject: [PATCH] Change division to use float for massInv

---
 .../src/TrackParametrizationWithError.cxx        | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/DataFormats/Reconstruction/src/TrackParametrizationWithError.cxx b/DataFormats/Reconstruction/src/TrackParametrizationWithError.cxx
index eb8071ec0073d..e029f76c98543 100644
--- a/DataFormats/Reconstruction/src/TrackParametrizationWithError.cxx
+++ b/DataFormats/Reconstruction/src/TrackParametrizationWithError.cxx
@@ -1364,12 +1364,12 @@ GPUd() bool TrackParametrizationWithError<value_T>::correctForMaterial(value_t x
   }
   auto m = this->getPID().getMass();
   int charge2 = this->getAbsCharge() * this->getAbsCharge();
-  value_t p = this->getP(), p0 = p, p02 = p * p, e2 = p02 + this->getPID().getMass2(), massInv = 1. / m, bg = p * massInv, dETot = 0.;
+  value_t p = this->getP(), p0 = p, p02 = p * p, e2 = p02 + this->getPID().getMass2(), massInv = 1.f / m, bg = p * massInv, dETot = 0.f;
   value_t e = gpu::CAMath::Sqrt(e2), e0 = e;
   if (m > 0 && xrho != 0.f) {
     value_t ekin = e - m, dedx = this->getdEdxBBOpt(bg);
 #ifdef _BB_NONCONST_CORR_
-    value_t dedxDer = 0., dedx1 = dedx;
+    value_t dedxDer = 0.f, dedx1 = dedx;
 #endif
     if (charge2 != 1) {
       dedx *= charge2;
@@ -1391,11 +1391,11 @@ GPUd() bool TrackParametrizationWithError<value_T>::correctForMaterial(value_t x
     }
     while (na--) {
 #ifdef _BB_NONCONST_CORR_
-      if (dedxDer != 0.) { // correction for non-constantness of dedx vs beta*gamma (in linear approximation): for a single step dE -> dE * [(exp(dedxDer) - 1)/dedxDer]
+      if (dedxDer != 0.f) { // correction for non-constantness of dedx vs beta*gamma (in linear approximation): for a single step dE -> dE * [(exp(dedxDer) - 1)/dedxDer]
         if (xrho < 0) {
           dedxDer = -dedxDer; // E.loss ( -> positive derivative)
         }
-        auto corrC = (gpu::CAMath::Exp(dedxDer) - 1.) / dedxDer;
+        auto corrC = (gpu::CAMath::Exp(dedxDer) - 1.f) / dedxDer;
         dE *= corrC;
       }
 #endif
@@ -1499,12 +1499,12 @@ GPUd() bool TrackParametrizationWithError<value_T>::correctForMaterial(TrackPara
   auto pid = linRef.getPID();
   auto m = pid.getMass();
   int charge2 = linRef.getAbsCharge() * linRef.getAbsCharge();
-  value_t p = linRef.getP(), p0 = p, p02 = p * p, e2 = p02 + pid.getMass2(), massInv = 1. / m, bg = p * massInv, dETot = 0.;
+  value_t p = linRef.getP(), p0 = p, p02 = p * p, e2 = p02 + pid.getMass2(), massInv = 1.f / m, bg = p * massInv, dETot = 0.f;
   value_t e = gpu::CAMath::Sqrt(e2), e0 = e;
   if (m > 0 && xrho != 0.f) {
     value_t ekin = e - m, dedx = this->getdEdxBBOpt(bg);
 #ifdef _BB_NONCONST_CORR_
-    value_t dedxDer = 0., dedx1 = dedx;
+    value_t dedxDer = 0.f, dedx1 = dedx;
 #endif
     if (charge2 != 1) {
       dedx *= charge2;
@@ -1526,11 +1526,11 @@ GPUd() bool TrackParametrizationWithError<value_T>::correctForMaterial(TrackPara
     }
     while (na--) {
 #ifdef _BB_NONCONST_CORR_
-      if (dedxDer != 0.) { // correction for non-constantness of dedx vs beta*gamma (in linear approximation): for a single step dE -> dE * [(exp(dedxDer) - 1)/dedxDer]
+      if (dedxDer != 0.f) { // correction for non-constantness of dedx vs beta*gamma (in linear approximation): for a single step dE -> dE * [(exp(dedxDer) - 1)/dedxDer]
         if (xrho < 0) {
           dedxDer = -dedxDer; // E.loss ( -> positive derivative)
         }
-        auto corrC = (gpu::CAMath::Exp(dedxDer) - 1.) / dedxDer;
+        auto corrC = (gpu::CAMath::Exp(dedxDer) - 1.f) / dedxDer;
         dE *= corrC;
       }
 #endif