Changeset d3165fa in git for external/TrackCovariance
- Timestamp:
- Mar 1, 2021, 4:29:02 PM (4 years ago)
- Branches:
- master
- Children:
- a1a25d4, df408d2
- Parents:
- a6e8029
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
external/TrackCovariance/TrkUtil.cc
ra6e8029 rd3165fa 30 30 Double_t r2 = x.Perp2(); 31 31 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); 33 33 Double_t phi0 = TMath::ATan2((p(1) - a * x(0)) / T, (p(0) + a * x(1)) / T); // Phi0 34 34 Double_t D; // Impact parameter D … … 40 40 Par(2) = C; // Store C 41 41 //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)); 43 43 Double_t st = TMath::ASin(B) / C; 44 44 Double_t ct = p(2) / pt; … … 106 106 pACTS(2) = Par(1); // Phi0 is unchanged 107 107 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 GeV108 pACTS(4) = Par(2) / (b * sqrt(1 + Par(4) * Par(4))); // q/p in GeV 109 109 pACTS(5) = 0.0; // Time: currently undefined 110 110 // … … 123 123 A(0, 0) = 1000.; // D-D conversion to mm 124 124 A(1, 2) = 1.0; // phi0-phi0 125 A(2, 4) = 1.0 / ( TMath::Sqrt(1.0 + ct * ct) * b); // q/p-C125 A(2, 4) = 1.0 / (sqrt(1.0 + ct * ct) * b); // q/p-C 126 126 A(3, 1) = 1000.; // z0-z0 conversion to mm 127 127 A(4, 3) = -1.0 / (1.0 + ct * ct); // theta - cot(theta)
Note:
See TracChangeset
for help on using the changeset viewer.