25 std::pair<TrajectoryStateOnSurface,double>
27 const Plane& plane)
const
35 return errorprop( fts, plane, side, gp.
parameters(),gp.
s());
39 std::pair< TrajectoryStateOnSurface, double>
48 return errorprop( fts, cyl, side, gp.
parameters(),gp.
s());
55 const Plane& plane)
const
60 double startZ = plane.
localZ(gpos);
67 if( fabs(rho)<1.
e-10 ) {
74 LogDebug(
"RKPropagatorInS")<<
" startZ = "<<startZ;
76 if (fabs(startZ) < 1
e-5){
77 LogDebug(
"RKPropagatorInS")<<
"Propagation is not performed: state is already on final surface.";
89 std::pair<bool,double> propResult = planeCrossing.
pathLength(plane);
91 double s = propResult.second;
98 LogDebug(
"RKPropagatorInS")<<
"Straight line propgation to plane failed !!";
109 LogDebug(
"RKPropagatorInS") <<
"The starting position is " << ts.
position() <<
" (global) "
114 LogDebug(
"RKPropagatorInS") <<
"The plane position is " << plane.
position() <<
" (global) "
115 << localPlane->position() <<
" (local) " ;
120 std::pair<bool,double> res3 =
cross.pathLength(plane);
121 LogDebug(
"RKPropagatorInS") <<
"straight line distance " << res3.first <<
" " << res3.second ;
141 while (safeGuard++<100) {
144 std::pair<bool,double>
path = pathLength( plane, startState.
position(),
146 (double) ts.
charge(), currentDirection);
148 LogDebug(
"RKPropagatorInS") <<
"RKPropagatorInS: Path length calculation to plane failed!"
150 <<
"...Local starting position in volume " << startState.
position()
157 LogDebug(
"RKPropagatorInS") <<
"RKPropagatorInS: Path lenght to plane is " << path.second ;
160 double sstep = path.second;
162 LogDebug(
"RKPropagatorInS") <<
"On-surface accuracy not reached, but pathLength calculation says we are there! "
163 <<
"path " << path.second <<
" distance to plane is " << startZ ;
168 LogDebug(
"RKPropagatorInS") <<
"RKPropagatorInS: Solving for " << sstep
169 <<
" current distance to plane is " << startZ ;
171 RKVector rkresult = solver( 0, start, sstep, deriv, dist, eps);
176 if ( fabs(remainingZ) < eps) {
177 LogDebug(
"RKPropagatorInS") <<
"On-surface accuracy reached! " << remainingZ ;
184 if (remainingZ * startZ > 0) {
185 LogDebug(
"RKPropagatorInS") <<
"Accuracy not reached yet, trying in same direction again "
189 LogDebug(
"RKPropagatorInS") <<
"Accuracy not reached yet, trying in opposite direction "
195 edm::LogError(
"FailedPropagation") <<
" too many iterations trying to reach plane ";
208 if (sp.
x()!=0. || sp.
y()!=0.) {
216 double startR = cyl.
radius() -
pos.perp();
227 if( fabs(rho)<1.
e-10 ) {
240 std::pair<bool,double> propResult = cylCrossing.
pathLength(cyl);
241 if ( propResult.first &&
theVolume !=0) {
242 double s = propResult.second;
246 LogDebug(
"RKPropagatorInS") <<
"Straight line propagation to cylinder succeeded !!";
250 edm::LogError(
"RKPropagatorInS") <<
"Straight line propagation to cylinder failed !!";
268 while (safeGuard++<100) {
272 currentDirection, eps);
276 LogDebug(
"RKPropagatorInS") <<
"RKPropagatorInS: Path length calculation to cylinder failed!"
281 LogDebug(
"RKPropagatorInS") <<
"RKPropagatorInS: Path lenght to cylinder is " << path.second
284 <<
" to R " << cyl.
radius()
288 double sstep = path.second;
290 LogDebug(
"RKPropagatorInS") <<
"accuracy not reached, but pathLength calculation says we are there! "
299 LogDebug(
"RKPropagatorInS") <<
"RKPropagatorInS: Solving for " << sstep
300 <<
" current distance to cylinder is " << startR ;
302 RKVector rkresult = solver( 0, start, sstep, deriv, dist, eps);
306 if ( fabs(remainingR) < eps) {
307 LogDebug(
"RKPropagatorInS") <<
"Accuracy reached! " << remainingR ;
315 if (remainingR * startR > 0) {
316 LogDebug(
"RKPropagatorInS") <<
"Accuracy not reached yet, trying in same direction again "
320 LogDebug(
"RKPropagatorInS") <<
"Accuracy not reached yet, trying in opposite direction "
326 edm::LogError(
"FailedPropagation") <<
" too many iterations trying to reach cylinder ";
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
virtual Propagator * clone() const
Local3DVector LocalVector
GlobalTrajectoryParameters gtpFromVolumeLocal(const CartesianStateAdaptor &state, TrackCharge charge) const
virtual PropagationDirection propagationDirection() const
ROOT::Math::Plane3D::Vector Vector
GlobalTrajectoryParameters gtpFromLocal(const Basic3DVector< double > &lpos, const Basic3DVector< double > &lmom, TrackCharge ch, const Surface &surf) const
GlobalParametersWithPath propagateParametersOnCylinder(const FreeTrajectoryState &ts, const Cylinder &cyl) const
T perp() const
Magnitude of transverse component.
RKLocalFieldProvider fieldProvider() const
std::pair< bool, double > pathLength(const Cylinder &cyl) const
Global3DPoint GlobalPoint
Basic3DVector cross(const Basic3DVector &v) const
Vector product, or "cross" product, with a vector of same type.
float localZ(const GlobalPoint &gp) const
Fast access to distance from plane for a point.
const GlobalTrajectoryParameters & parameters() const
TrackCharge charge() const
Basic3DVector< double > rkMomentum(const GlobalVector &mom) const
std::pair< TrajectoryStateOnSurface, double > TsosWP
Estimator of the distance between two state vectors, e.g. for convergence test.
Scalar radius() const
Radius of the cylinder.
T z() const
Cartesian z coordinate.
PropagationDirection invertDirection(PropagationDirection dir) const
LocalPoint toLocal(const GlobalPoint &gp) const
virtual std::pair< TrajectoryStateOnSurface, double > propagateWithPath(const FreeTrajectoryState &, const Plane &) const
Derivative calculation for the 6D cartesian case.
static RKSmallVector< double, 6 > rkstate(const Vector3D &pos, const Vector3D &mom)
const Vector3D & momentum() const
virtual TrajectoryStateOnSurface propagate(const FreeTrajectoryState &, const Plane &) const
const MagVolume * theVolume
T perp() const
Magnitude of transverse component.
GlobalVector momentum() const
std::pair< bool, double > pathLength(const Plane &plane) const
PositionType position(double s) const
GlobalPoint toGlobal(const LocalPoint &lp) const
GlobalPoint position() const
const Frame & frame() const
The reference frame in which the field is defined.
static Plane transformPlane(const Plane &plane, const GloballyPositioned< T > &frame)
GlobalPoint position(const double s) const
double transverseCurvature() const
GlobalPoint globalPosition(const Basic3DVector< double > &pos) const
GlobalParametersWithPath propagateParametersOnPlane(const FreeTrajectoryState &ts, const Plane &plane) const
RKPropagatorInS(const MagVolume &vol, PropagationDirection dir=alongMomentum, double tolerance=5.e-5)
const Vector3D & position() const
std::pair< bool, double > pathLength(const Cylinder &cyl) const
const BasicVectorType & basicVector() const
const PositionType & position() const
GlobalVector globalMomentum(const Basic3DVector< double > &mom) const
Basic3DVector< double > rkPosition(const GlobalPoint &pos) const
Global3DVector GlobalVector
Vector inTesla(const LocalPoint &lp) const
the argument lp is in the local frame specified in the constructor