207{
208 HepVector tdsHelix = kalTrack->
getFHelix(pid);
209 HepSymMatrix tdsMatrix = kalTrack->
getFError(pid);
210
212
213 KalFitTrack fitTrack(IP, tdsHelix, tdsMatrix, pid, 0, 0);
214
215
216 const double rp = point.perp();
217
218
219 const double radius = m_BesKalmanExtWalls[0].radius();
220 const double dphi = fitTrack.intersect_cylinder(std::max(rp, radius));
221 const HepPoint3D lastPivot = fitTrack.x(dphi);
222 fitTrack.pivot(lastPivot);
223 if (rp <= radius) extToAnyPoint(fitTrack, point);
224
225
226
227 fitTrack.pivot(IP);
228 setHelixVector(fitTrack.a());
229 setErrorMatrix(fitTrack.Ea());
230}
HepGeom::Point3D< double > HepPoint3D
const HepSymMatrix & getFError(const int pid) const
const HepVector & getFHelix(const int pid) const