Fork me on GitHub

Changeset f9517a5 in git for external/TrackCovariance


Ignore:
Timestamp:
Mar 4, 2021, 11:31:58 AM (4 years ago)
Author:
Pavel Demin <pavel.demin@…>
Branches:
master
Children:
dc883b4
Parents:
4586f8f
Message:

include TMath.h to TrkUtil.cc

File:
1 edited

Legend:

Unmodified
Added
Removed
  • external/TrackCovariance/TrkUtil.cc

    r4586f8f rf9517a5  
    11#include "TrkUtil.h"
     2#include <TMath.h>
    23#include <iostream>
    34
     
    3132        Double_t r2 = x.Perp2();
    3233        Double_t cross = x(0) * p(1) - x(1) * p(0);
    33         Double_t T = sqrt(pt * pt - 2 * a * cross + a * a * r2);
    34         Double_t phi0 = atan2((p(1) - a * x(0)) / T, (p(0) + a * x(1)) / T);    // Phi0
     34        Double_t T = TMath::Sqrt(pt * pt - 2 * a * cross + a * a * r2);
     35        Double_t phi0 = TMath::ATan2((p(1) - a * x(0)) / T, (p(0) + a * x(1)) / T);     // Phi0
    3536        Double_t D;                                                     // Impact parameter D
    3637        if (pt < 10.0) D = (T - pt) / a;
     
    4142        Par(2) = C;             // Store C
    4243        //Longitudinal parameters
    43         Double_t B = C * sqrt(TMath::Max(r2 - D * D, 0.0) / (1 + 2 * C * D));
    44         Double_t st = asin(B) / C;
     44        Double_t B = C * TMath::Sqrt(TMath::Max(r2 - D * D, 0.0) / (1 + 2 * C * D));
     45        Double_t st = TMath::ASin(B) / C;
    4546        Double_t ct = p(2) / pt;
    4647        Double_t z0 = x(2) - ct * st;
     
    6970        //
    7071        TVector3 Xval;
    71         Xval(0) = -D * sin(phi0);
    72         Xval(1) = D * cos(phi0);
     72        Xval(0) = -D * TMath::Sin(phi0);
     73        Xval(1) = D * TMath::Cos(phi0);
    7374        Xval(2) = z0;
    7475        //
     
    9293        TVector3 Pval;
    9394        Double_t pt = Bz * cSpeed() / TMath::Abs(2 * C);
    94         Pval(0) = pt * cos(phi0);
    95         Pval(1) = pt * sin(phi0);
     95        Pval(0) = pt * TMath::Cos(phi0);
     96        Pval(1) = pt * TMath::Sin(phi0);
    9697        Pval(2) = pt * ct;
    9798        //
     
    114115        pACTS(1) = 1000 * Par(3);               // z0 from m to mm
    115116        pACTS(2) = Par(1);                      // Phi0 is unchanged
    116         pACTS(3) = atan2(1.0, Par(4));          // Theta in [0, pi] range
    117         pACTS(4) = Par(2) / (b * sqrt(1 + Par(4) * Par(4)));            // q/p in GeV
     117        pACTS(3) = TMath::ATan2(1.0, Par(4));           // Theta in [0, pi] range
     118        pACTS(4) = Par(2) / (b * TMath::Sqrt(1 + Par(4) * Par(4)));             // q/p in GeV
    118119        pACTS(5) = 0.0;                         // Time: currently undefined
    119120        //
     
    132133        A(0, 0) = 1000.;                // D-D  conversion to mm
    133134        A(1, 2) = 1.0;          // phi0-phi0
    134         A(2, 4) = 1.0 / (sqrt(1.0 + ct * ct) * b);      // q/p-C
     135        A(2, 4) = 1.0 / (TMath::Sqrt(1.0 + ct * ct) * b);       // q/p-C
    135136        A(3, 1) = 1000.;                // z0-z0 conversion to mm
    136137        A(4, 3) = -1.0 / (1.0 + ct * ct); // theta - cot(theta)
    137         A(4, 4) = -C * ct / (b * pow(1.0 + ct * ct, 3.0 / 2.0)); // q/p-cot(theta)
     138        A(4, 4) = -C * ct / (b * TMath::Power(1.0 + ct * ct, 3.0 / 2.0)); // q/p-cot(theta)
    138139        //
    139140        TMatrixDSym Cv = Cov;
Note: See TracChangeset for help on using the changeset viewer.