Changes in modules/TrackCovariance.cc [3051ea17:46c8df8] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TrackCovariance.cc
r3051ea17 r46c8df8 96 96 Candidate *candidate, *mother; 97 97 Double_t mass, p, pt, q, ct; 98 Double_t dd0, ddz, dphi, dct, dp, dpt , dC;98 Double_t dd0, ddz, dphi, dct, dp, dpt; 99 99 100 100 … … 114 114 candidate->Momentum.SetVectM(track.GetObsP(), mass); 115 115 candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T()); 116 117 // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta))118 candidate->TrackCovariance = track.GetCov();119 116 120 117 pt = candidate->Momentum.Pt(); … … 128 125 129 126 candidate->D0 = track.GetObsPar()[0]; 127 candidate->DZ = track.GetObsPar()[3]; 128 candidate->P = track.GetObsP().Mag(); 129 candidate->CtgTheta = track.GetObsPar()[4]; 130 130 candidate->Phi = track.GetObsPar()[1]; 131 candidate->C = track.GetObsPar()[2]; 132 candidate->DZ = track.GetObsPar()[3]; 133 candidate->CtgTheta = track.GetObsPar()[4]; 134 candidate->P = track.GetObsP().Mag(); 131 135 132 candidate->PT = pt; 136 133 candidate->Charge = q; … … 142 139 dpt = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz); 143 140 dp = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct)); 144 dC = TMath::Sqrt(track.GetCov()(2, 2));145 141 146 142 candidate->ErrorD0 = dd0; 147 143 candidate->ErrorDZ = ddz; 148 144 candidate->ErrorP = dp; 149 candidate->ErrorC = dC;150 145 candidate->ErrorCtgTheta = dct; 151 146 candidate->ErrorPhi = dphi; … … 153 148 //candidate->TrackResolution = dpt / pt; 154 149 candidate->TrackResolution = dp / p; 155 156 150 157 151 candidate->AddCandidate(mother);
Note:
See TracChangeset
for help on using the changeset viewer.