21 double dydx = 0., dxdy = 0.;
22 double pt = 0., px = 0., py = 0.;
31 pt = 0.01 * rho * (0.3*
tesla0.
z());
46 if (arg1<0.0 && arg2<0.0) {
49 px = pt/
sqrt(1. + dydx*dydx);
55 }
else if ( arg1>arg2 ) {
60 px = pt/
sqrt(1. + dydx*dydx);
67 py = pt/
sqrt(1. + dxdy*dxdy);
71 if(px*(pMid.
x() - v.
x()) + py*(pMid.
y() - v.
y()) < 0.) {
83 double dzdrphi = -flfit.
n1()/flfit.n2();
84 double pz = pt*dzdrphi;
103 double z_0 = -flfit.c()/flfit.n2();
121 double pt = 0., px = 0., py = 0.;
128 px = pt/
sqrt(1. + dydx*dydx);
131 if (px*(pMid.
x() - v.
x()) + py*(pMid.
y() - v.
y()) < 0.) {
144 double dzdr = -flfit.
n1()/flfit.
n2();
158 double z_0 = -flfit.
c()/flfit.
n2();
FTS stateAtVertex() const
FTS helixStateAtVertex() const
ROOT::Math::SMatrixIdentity AlgebraicMatrixID
ROOT::Math::SMatrix< double, 6, 6, ROOT::Math::MatRepSym< double, 6 > > AlgebraicSymMatrix66
Global3DPoint GlobalPoint
CLHEP::HepSymMatrix AlgebraicSymMatrix
edm::ESHandle< MagneticField > pSetup
Global3DVector GlobalVector
FTS straightLineStateAtVertex() const