30 double cosl = tn.
perp();
if (cosl < 1.
e-30) cosl = 1.e-30;
31 double cosl1 = 1./cosl;
35 double uj = un.
dot(dj);
36 double uk = un.
dot(dk);
37 double vj = vn.dot(dj);
38 double vk = vn.dot(dk);
41 double t1r = 1./tnl.
z();
57 double sinz =-un.
dot(hdir);
58 double cosz = vn.dot(hdir);
59 double ui = un.
dot(di)*(q*t3r);
60 double vi = vn.dot(di)*(q*t3r);
Basic3DVector< T > z() const
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
CLHEP::HepMatrix asHepMatrix(const ROOT::Math::SMatrix< double, N1, N2, typename ROOT::Math::MatRepStd< double, N1, N2 > > &rm)
Basic3DVector< T > y() const
LocalPoint position() const
Local x and y position coordinates.
virtual GlobalVector inTesla(const GlobalPoint &gp) const =0
Field value ad specified global point, in Tesla.
PreciseFloatType< T, U >::Type dot(const Vector3DBase< U, FrameTag > &v) const
JacobianCurvilinearToLocal(const Surface &surface, const LocalTrajectoryParameters &localParameters, const MagneticField &magField)
AlgebraicMatrix55 theJacobian
Basic3DVector< T > x() const
CLHEP::HepMatrix AlgebraicMatrix
const AlgebraicMatrix55 & jacobian() const
LocalVector momentum() const
Momentum vector in the local frame.
Vector3DBase unit() const
const RotationType & rotation() const
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
ROOT::Math::SMatrix< double, 5, 5, ROOT::Math::MatRepStd< double, 5, 5 > > AlgebraicMatrix55
double signedInverseMomentum() const
Signed inverse momentum q/p (zero for neutrals).
const AlgebraicMatrix jacobian_old() const