Changeset 92b8d11 in git

Timestamp:

Mar 2, 2021, 3:34:21 PM (4 years ago)

Author:

Franco BEDESCHI <bed@…>

Branches:

master

Children:

537d24f, a1a25d4

Parents:

df408d2

Message:

Fix Sqrt for ROOT5

File:

• external/TrackCovariance/TrkUtil.cc (modified) (4 diffs)

external/TrackCovariance/TrkUtil.cc

@@ -31 +31 @@
         Double_t r2 = x.Perp2();
         Double_t cross = x(0) * p(1) - x(1) * p(0);
-        Double_t T = TMath::Sqrt(pt * pt - 2 * a * cross + a * a * r2);
+        Double_t T = sqrt(pt * pt - 2 * a * cross + a * a * r2);
         Double_t phi0 = atan2((p(1) - a * x(0)) / T, (p(0) + a * x(1)) / T);    // Phi0
         Double_t D;                                                     // Impact parameter D
@@ -41 +41 @@
         Par(2) = C;             // Store C
         //Longitudinal parameters
-        Double_t B = C * TMath::Sqrt(TMath::Max(r2 - D * D, 0.0) / (1 + 2 * C * D));
+        Double_t B = C * sqrt(TMath::Max(r2 - D * D, 0.0) / (1 + 2 * C * D));
         Double_t st = asin(B) / C;
         Double_t ct = p(2) / pt;
@@ -115 +115 @@
         pACTS(2) = Par(1);                      // Phi0 is unchanged
         pACTS(3) = atan2(1.0, Par(4));          // Theta in [0, pi] range
-        pACTS(4) = Par(2) / (b * TMath::Sqrt(1 + Par(4) * Par(4)));             // q/p in GeV
+        pACTS(4) = Par(2) / (b * sqrt(1 + Par(4) * Par(4)));            // q/p in GeV
         pACTS(5) = 0.0;                         // Time: currently undefined
         //
@@ -132 +132 @@
         A(0, 0) = 1000.;                // D-D  conversion to mm
         A(1, 2) = 1.0;          // phi0-phi0
-        A(2, 4) = 1.0 / (TMath::Sqrt(1.0 + ct * ct) * b);       // q/p-C
+        A(2, 4) = 1.0 / (sqrt(1.0 + ct * ct) * b);      // q/p-C
         A(3, 1) = 1000.;                // z0-z0 conversion to mm
         A(4, 3) = -1.0 / (1.0 + ct * ct); // theta - cot(theta)