Fork me on GitHub

Changeset d3165fa in git


Ignore:
Timestamp:
Mar 1, 2021, 4:29:02 PM (4 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
master
Children:
a1a25d4, df408d2
Parents:
a6e8029
Message:

comply with ROOT5 TMath definition

File:
1 edited

Legend:

Unmodified
Added
Removed
  • external/TrackCovariance/TrkUtil.cc

    ra6e8029 rd3165fa  
    3030        Double_t r2 = x.Perp2();
    3131        Double_t cross = x(0) * p(1) - x(1) * p(0);
    32         Double_t T = TMath::Sqrt(pt * pt - 2 * a * cross + a * a * r2);
     32        Double_t T = sqrt(pt * pt - 2 * a * cross + a * a * r2);
    3333        Double_t phi0 = TMath::ATan2((p(1) - a * x(0)) / T, (p(0) + a * x(1)) / T);     // Phi0
    3434        Double_t D;                                                     // Impact parameter D
     
    4040        Par(2) = C;             // Store C
    4141        //Longitudinal parameters
    42         Double_t B = C * TMath::Sqrt(TMath::Max(r2 - D * D, 0.0) / (1 + 2 * C * D));
     42        Double_t B = C * sqrt(TMath::Max(r2 - D * D, 0.0) / (1 + 2 * C * D));
    4343        Double_t st = TMath::ASin(B) / C;
    4444        Double_t ct = p(2) / pt;
     
    106106        pACTS(2) = Par(1);                      // Phi0 is unchanged
    107107        pACTS(3) = TMath::ATan2(1.0, Par(4));           // Theta in [0, pi] range
    108         pACTS(4) = Par(2) / (b * TMath::Sqrt(1 + Par(4) * Par(4)));             // q/p in GeV
     108        pACTS(4) = Par(2) / (b * sqrt(1 + Par(4) * Par(4)));            // q/p in GeV
    109109        pACTS(5) = 0.0;                         // Time: currently undefined
    110110        //
     
    123123        A(0, 0) = 1000.;                // D-D  conversion to mm
    124124        A(1, 2) = 1.0;          // phi0-phi0
    125         A(2, 4) = 1.0 / (TMath::Sqrt(1.0 + ct * ct) * b);       // q/p-C
     125        A(2, 4) = 1.0 / (sqrt(1.0 + ct * ct) * b);      // q/p-C
    126126        A(3, 1) = 1000.;                // z0-z0 conversion to mm
    127127        A(4, 3) = -1.0 / (1.0 + ct * ct); // theta - cot(theta)
Note: See TracChangeset for help on using the changeset viewer.