Changes in external/TrackCovariance/ObsTrk.cc [a0f5d71:942a705] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
external/TrackCovariance/ObsTrk.cc
ra0f5d71 r942a705 20 20 fGenPar.ResizeTo(5); 21 21 fGenParACTS.ResizeTo(6); 22 fGenParILC.ResizeTo(5);23 22 fObsPar.ResizeTo(5); 24 23 fObsParACTS.ResizeTo(6); 25 fObsParILC.ResizeTo(5);26 24 fCov.ResizeTo(5, 5); 27 25 fCovACTS.ResizeTo(6, 6); 28 fCovILC.ResizeTo(5, 5);29 26 fGenPar = XPtoPar(x,p,Q); 30 27 fGenParACTS = ParToACTS(fGenPar); 31 fGenParILC = ParToILC(fGenPar);32 28 /* 33 29 std::cout << "ObsTrk::ObsTrk: fGenPar"; … … 37 33 fObsPar = GenToObsPar(fGenPar, fGC); 38 34 fObsParACTS = ParToACTS(fObsPar); 39 fObsParILC = ParToILC(fObsPar);40 35 fObsX = ParToX(fObsPar); 41 36 fObsP = ParToP(fObsPar); 42 37 fObsQ = ParToQ(fObsPar); 43 38 fCovACTS = CovToACTS(fCov); 44 fCovILC = CovToILC(fCov);45 39 } 46 40 // … … 212 206 } 213 207 214 // Parameter conversion to ILC format215 TVectorD ObsTrk::ParToILC(TVectorD Par)216 {217 TVectorD pILC(5); // Return vector218 //219 pILC(0) = Par(0)*1.0e3; // d0 in mm220 pILC(1) = Par(1); // phi0 is unchanged221 pILC(2) = -2 * Par(2)*1.0e-3; // w in mm^-1222 pILC(3) = Par(3)*1.0e3; // z0 in mm223 pILC(4) = Par(4); // tan(lambda) = cot(theta)224 //225 return pILC;226 }227 // Covariance conversion to ILC format228 TMatrixDSym ObsTrk::CovToILC(TMatrixDSym Cov)229 {230 TMatrixDSym cILC(5); cILC.Zero();231 //232 // Fill derivative matrix233 TMatrixD A(5, 5); A.Zero();234 //235 A(0, 0) = 1.0e3; // D-d0 in mm236 A(1, 1) = 1.0; // phi0-phi0237 A(2, 2) = -2.0e-3; // w-C238 A(3, 3) = 1.0e3; // z0-z0 conversion to mm239 A(4, 4) = 1.0; // tan(lambda) - cot(theta)240 //241 TMatrixDSym Cv = Cov;242 TMatrixD At(5, 5);243 At.Transpose(A);244 Cv.Similarity(At);245 cILC = Cv;246 //247 return cILC;248 }249 250 208 251 209
Note:
See TracChangeset
for help on using the changeset viewer.