From 372f09bd2e26d412887436b539f7d2350c238bcc Mon Sep 17 00:00:00 2001
From: Alvaro Tolosa Delgado <alvaro.tolosa.delgado@cern.ch>
Date: Mon, 6 May 2024 13:27:54 +0200
Subject: [PATCH] fix missing unit conversion factors
---
DDRec/include/DDRec/DCH_info.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/DDRec/include/DDRec/DCH_info.h b/DDRec/include/DDRec/DCH_info.h
index 17bb22198..6dbec27b8 100644
--- a/DDRec/include/DDRec/DCH_info.h
+++ b/DDRec/include/DDRec/DCH_info.h
@@ -122,18 +122,18 @@ public:
/// tan(stereoangle) = R(z=0) / (L/2) * tan( twist_angle/2)
DCH_angle_t stereoangle_z0(DCH_length_t r_z0) const {
- return atan( r_z0/Lhalf*tan(twist_angle/2));
+ return atan( r_z0/Lhalf*tan(twist_angle/2/dd4hep::rad));
}
/// tan(stereoangle) = R(z=L/2) / (L/2) * sin( twist_angle/2)
DCH_angle_t stereoangle_zLhalf(DCH_length_t r_zLhalf) const {
- return atan( r_zLhalf/Lhalf*sin(twist_angle/2));
+ return atan( r_zLhalf/Lhalf*sin(twist_angle/2/dd4hep::rad));
}
/// WireLength = 2*dch_Lhalf/cos(atan(Pitch_z0(r_z0)/(2*dch_Lhalf)))/cos(stereoangle_z0(r_z0))
DCH_length_t WireLength(int nlayer, DCH_length_t r_z0) const {
auto Pitch_z0 = database.at(nlayer).Pitch_z0(r_z0);
- return 2*Lhalf/cos(atan(Pitch_z0/(2*Lhalf)))/cos(stereoangle_z0(r_z0)) ;
+ return 2*Lhalf/cos(atan(Pitch_z0/(2*Lhalf)))/cos(stereoangle_z0(r_z0)/dd4hep::rad) ;
};
/// Internal helper struct for defining the layer layout
--
GitLab