RKPropagatorInS Class Reference

#include <RKPropagatorInS.h>

Inheritance diagram for RKPropagatorInS:

Public Member Functions
Propagator *	clone () const override
const MagneticField *	magneticField () const override
	RKPropagatorInS (const MagVolume &vol, PropagationDirection dir=alongMomentum, double tolerance=5.e-5)
	~RKPropagatorInS () override
Public Member Functions inherited from Propagator
template<typename STA , typename SUR >
TrajectoryStateOnSurface	propagate (STA const &state, SUR const &surface) const
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest) const final
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest1, const GlobalPoint &pDest2) const final
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &ftsStart, const reco::BeamSpot &beamSpot) const final
virtual std::pair< TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Surface &) const final
virtual std::pair< TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Surface &sur) const final
virtual std::pair< TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Plane &sur) const
virtual std::pair< TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Cylinder &sur) const
virtual std::pair< FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest) const
virtual std::pair< FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest1, const GlobalPoint &pDest2) const
	Propagate to PCA to a line (given by 2 points) given a starting point. More...
virtual std::pair< FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &ftsStart, const reco::BeamSpot &beamSpot) const
	Propagate to PCA to a line (given by beamSpot position and slope) given a starting point. More...
virtual PropagationDirection	propagationDirection () const final
	Propagator (PropagationDirection dir=alongMomentum)
virtual bool	setMaxDirectionChange (float phiMax)
virtual void	setPropagationDirection (PropagationDirection dir)
virtual	~Propagator ()

Private Types
typedef std::pair< TrajectoryStateOnSurface, double >	TsosWP

Private Member Functions
RKLocalFieldProvider	fieldProvider () const
RKLocalFieldProvider	fieldProvider (const Cylinder &cyl) const
GlobalVector	globalMomentum (const Basic3DVector< float > &mom) const
GlobalPoint	globalPosition (const Basic3DVector< float > &pos) const
GlobalTrajectoryParameters	gtpFromLocal (const Basic3DVector< float > &lpos, const Basic3DVector< float > &lmom, TrackCharge ch, const Surface &surf) const
GlobalTrajectoryParameters	gtpFromVolumeLocal (const CartesianStateAdaptor &state, TrackCharge charge) const
PropagationDirection	invertDirection (PropagationDirection dir) const
GlobalParametersWithPath	propagateParametersOnCylinder (const FreeTrajectoryState &ts, const Cylinder &cyl) const
GlobalParametersWithPath	propagateParametersOnPlane (const FreeTrajectoryState &ts, const Plane &plane) const
std::pair< TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Plane &) const override
std::pair< TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Cylinder &) const override
Basic3DVector< double >	rkMomentum (const GlobalVector &mom) const
Basic3DVector< double >	rkPosition (const GlobalPoint &pos) const

Private Attributes
double	theTolerance
const MagVolume *	theVolume

Detailed Description

Definition at line 20 of file RKPropagatorInS.h.

Member Typedef Documentation

typedef std::pair<TrajectoryStateOnSurface,double> RKPropagatorInS::TsosWP

private

Definition at line 52 of file RKPropagatorInS.h.

Constructor & Destructor Documentation

RKPropagatorInS::RKPropagatorInS ( const MagVolume &  vol, PropagationDirection  dir = alongMomentum, double  tolerance = 5.e-5  )

inlineexplicit

Definition at line 27 of file RKPropagatorInS.h.

Referenced by clone().

                                          : 
    Propagator(dir), theVolume( &vol), theTolerance( tolerance) {}

RKPropagatorInS::~RKPropagatorInS ( )

inlineoverride

Definition at line 31 of file RKPropagatorInS.h.

References clone(), Propagator::propagate(), propagateWithPath(), and Propagator::propagateWithPath().

{}

Member Function Documentation

Propagator * RKPropagatorInS::clone ( void  ) const

overridevirtual

Implements Propagator.

Definition at line 332 of file RKPropagatorInS.cc.

References RKPropagatorInS().

Referenced by ~RKPropagatorInS().

{
    return new RKPropagatorInS(*this);
}

RKLocalFieldProvider RKPropagatorInS::fieldProvider ( ) const

private

Definition at line 345 of file RKPropagatorInS.cc.

References theVolume.

Referenced by propagateParametersOnCylinder(), and propagateParametersOnPlane().

{
  return RKLocalFieldProvider( *theVolume);
}

RKLocalFieldProvider RKPropagatorInS::fieldProvider ( const Cylinder &  cyl ) const

private

Definition at line 350 of file RKPropagatorInS.cc.

References theVolume.

{
  return RKLocalFieldProvider( *theVolume, cyl);
}

GlobalVector RKPropagatorInS::globalMomentum ( const Basic3DVector< float > &  mom ) const

private

Definition at line 379 of file RKPropagatorInS.cc.

References theVolume, and GloballyPositioned< T >::toGlobal().

Referenced by gtpFromVolumeLocal().

{
  if (theVolume != nullptr) return theVolume->toGlobal( LocalVector(mom));
  else return GlobalVector(mom);
}

GlobalPoint RKPropagatorInS::globalPosition ( const Basic3DVector< float > &  pos ) const

private

Definition at line 373 of file RKPropagatorInS.cc.

References theVolume, and GloballyPositioned< T >::toGlobal().

Referenced by gtpFromVolumeLocal(), and propagateParametersOnPlane().

{
  if (theVolume != nullptr) return theVolume->toGlobal( LocalPoint(pos));
  else return GlobalPoint(pos);
}

GlobalTrajectoryParameters RKPropagatorInS::gtpFromLocal ( const Basic3DVector< float > &  lpos, const Basic3DVector< float > &  lmom, TrackCharge  ch, const Surface &  surf  ) const

private

Definition at line 337 of file RKPropagatorInS.cc.

References theVolume, and Surface::toGlobal().

Referenced by propagateParametersOnCylinder().

{
    return GlobalTrajectoryParameters( surf.toGlobal( LocalPoint( lpos)),
                       surf.toGlobal( LocalVector( lmom)), ch, theVolume);
}

GlobalTrajectoryParameters RKPropagatorInS::gtpFromVolumeLocal ( const CartesianStateAdaptor &  state, TrackCharge  charge  ) const

private

Definition at line 387 of file RKPropagatorInS.cc.

References globalMomentum(), globalPosition(), CartesianStateAdaptor::momentum(), CartesianStateAdaptor::position(), and theVolume.

Referenced by propagateParametersOnPlane().

{
  return GlobalTrajectoryParameters( globalPosition(state.position()), 
                     globalMomentum(state.momentum()), 
                     charge, theVolume);
}

PropagationDirection RKPropagatorInS::invertDirection ( PropagationDirection  dir ) const

private

Definition at line 355 of file RKPropagatorInS.cc.

References alongMomentum, anyDirection, dir, and oppositeToMomentum.

Referenced by propagateParametersOnCylinder(), and propagateParametersOnPlane().

{
  if (dir == anyDirection) return dir;
  return ( dir == alongMomentum ? oppositeToMomentum : alongMomentum);
}

const MagneticField* RKPropagatorInS::magneticField ( ) const

inlineoverridevirtual

Implements Propagator.

Definition at line 48 of file RKPropagatorInS.h.

References theVolume.

{return theVolume;}

GlobalParametersWithPath RKPropagatorInS::propagateParametersOnCylinder ( const FreeTrajectoryState &  ts, const Cylinder &  cyl  ) const

private

Definition at line 198 of file RKPropagatorInS.cc.

References funct::abs(), FreeTrajectoryState::charge(), MillePedeFileConverter_cfg::e, fieldProvider(), gtpFromLocal(), invertDirection(), likely, LogDebug, CartesianStateAdaptor::momentum(), FreeTrajectoryState::momentum(), callgraph::path, StraightLineBarrelCylinderCrossing::pathLength(), StraightLineCylinderCrossing::pathLength(), Basic3DVector< T >::perp(), perp(), CartesianStateAdaptor::position(), StraightLineBarrelCylinderCrossing::position(), GloballyPositioned< T >::position(), FreeTrajectoryState::position(), Propagator::propagationDirection(), Cylinder::radius(), rho, CartesianStateAdaptor::rkstate(), alignCSCRings::s, command_line::start, theTolerance, theVolume, GloballyPositioned< T >::toLocal(), FreeTrajectoryState::transverseCurvature(), unlikely, x, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and Basic3DVector< T >::z().

Referenced by propagateWithPath().

                                                   {

  typedef RKAdaptiveSolver<double,RKOneCashKarpStep, 6>   Solver;
  typedef Solver::Vector                                  RKVector;
  
  
  const GlobalPoint& sp = cyl.position();
  if unlikely(sp.x()!=0. || sp.y()!=0.) {
      throw PropagationException("Cannot propagate to an arbitrary cylinder");
    }
  
  GlobalPoint gpos( ts.position());
  GlobalVector gmom( ts.momentum());
  LocalPoint pos(cyl.toLocal(gpos));
  LocalVector mom(cyl.toLocal(gmom));
  double startR = cyl.radius() - pos.perp();
  
  // LogDebug("RKPropagatorInS")  << "RKPropagatorInS: starting from FTS " << ts ;
  
  
  // (transverse) curvature
  double rho = ts.transverseCurvature();
  //
  // Straight line approximation? |rho|<1.e-10 equivalent to ~ 1um 
  // difference in transversal position at 10m.
  //
  if  unlikely( fabs(rho)<1.e-10 ) {
      //
      // Instantiate auxiliary object for finding intersection.
      // Frame-independant point and vector are created explicitely to 
      // avoid confusing gcc (refuses to compile with temporary objects
      // in the constructor).
      //
      
      StraightLineBarrelCylinderCrossing cylCrossing(gpos,gmom,propagationDirection());
      
      //
      // get solution
      //
      std::pair<bool,double> propResult = cylCrossing.pathLength(cyl);
      if  likely( propResult.first && theVolume !=nullptr) {
      double s = propResult.second;
      // point (reconverted to GlobalPoint)
      GlobalPoint x (cylCrossing.position(s));
      GlobalTrajectoryParameters res( x, gmom, ts.charge(), theVolume);
      LogDebug("RKPropagatorInS")  << "Straight line propagation to cylinder succeeded !!";
      return GlobalParametersWithPath( res, s);
    } 
      
      //do someting 
      edm::LogError("RKPropagatorInS")  <<"Straight line propagation to cylinder failed !!";
      return GlobalParametersWithPath( );
      
    }
  
  
  RKLocalFieldProvider field( fieldProvider(cyl));
  // StraightLineCylinderCrossing pathLength( pos, mom, propagationDirection());
  CartesianLorentzForce deriv(field, ts.charge());
  
  RKCartesianDistance dist;
  double eps = theTolerance;
  Solver solver;
  double stot = 0;
  PropagationDirection currentDirection = propagationDirection();
  
  RKVector start( CartesianStateAdaptor::rkstate( pos.basicVector(), mom.basicVector()));
  int safeGuard = 0;
  while (safeGuard++<100) {
    CartesianStateAdaptor startState(start);
    StraightLineCylinderCrossing pathLength( LocalPoint(startState.position()), 
                         LocalVector(startState.momentum()), 
                         currentDirection, eps);
    
    std::pair<bool,double> path = pathLength.pathLength( cyl);
    if  unlikely(!path.first) { 
    LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path length calculation to cylinder failed!" 
                     << "Radius " << cyl.radius() << " pos.perp() " << LocalPoint(startState.position()).perp() ;
    return GlobalParametersWithPath();
      }
    
    LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path lenght to cylinder is " << path.second 
                 << " from point (R,z) " << startState.position().perp() 
                 << ", " << startState.position().z()
                 << " to R " << cyl.radius() 
      ;
    
    
    double sstep = path.second;
    if  unlikely( std::abs(sstep) < eps) {
    LogDebug("RKPropagatorInS")  << "accuracy not reached, but pathLength calculation says we are there! "
                     << path.second ;
    
    GlobalTrajectoryParameters res( gtpFromLocal( startState.position(), 
                              startState.momentum(), 
                              ts.charge(), cyl));
    return GlobalParametersWithPath( res, stot);
      }
    
    LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Solving for " << sstep 
                 << " current distance to cylinder is " << startR ;
    
    RKVector rkresult = solver( 0, start, sstep, deriv, dist, eps);
    stot += sstep;
    CartesianStateAdaptor cur( rkresult);
    double remainingR = cyl.radius() - cur.position().perp();

    if ( fabs(remainingR) < eps) {
      LogDebug("RKPropagatorInS")  << "Accuracy reached! " << remainingR ;
      GlobalTrajectoryParameters res( gtpFromLocal( cur.position(), 
                            cur.momentum(), 
                            ts.charge(), cyl));
      return GlobalParametersWithPath( res, stot);
    }
    
    start = rkresult;
    if (remainingR * startR > 0) {
      LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in same direction again " 
                   << remainingR ;
    }
    else {
      LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in opposite direction " 
                   << remainingR ;
      currentDirection = invertDirection( currentDirection);
    }
    startR = remainingR;
  }
  
  edm::LogError("FailedPropagation") << " too many iterations trying to reach cylinder ";
  return GlobalParametersWithPath();
}

GlobalParametersWithPath RKPropagatorInS::propagateParametersOnPlane ( const FreeTrajectoryState &  ts, const Plane &  plane  ) const

private

Definition at line 51 of file RKPropagatorInS.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::basicVector(), FreeTrajectoryState::charge(), cross(), dir, MillePedeFileConverter_cfg::e, fieldProvider(), RKLocalFieldProvider::frame(), globalPosition(), gtpFromVolumeLocal(), RKLocalFieldProvider::inTesla(), invertDirection(), likely, Plane::localZ(), LogDebug, CartesianStateAdaptor::momentum(), FreeTrajectoryState::momentum(), callgraph::path, StraightLinePlaneCrossing::pathLength(), CartesianStateAdaptor::position(), GloballyPositioned< T >::position(), StraightLinePlaneCrossing::position(), FreeTrajectoryState::position(), Propagator::propagationDirection(), rho, rkMomentum(), rkPosition(), CartesianStateAdaptor::rkstate(), alignCSCRings::s, command_line::start, theTolerance, theVolume, GloballyPositioned< T >::toGlobal(), GloballyPositioned< T >::toLocal(), FrameChanger::transformPlane(), FreeTrajectoryState::transverseCurvature(), unlikely, and x.

Referenced by propagateWithPath().

{

  GlobalPoint gpos( ts.position());
  GlobalVector gmom( ts.momentum());
  double startZ = plane.localZ(gpos);
  // (transverse) curvature
  double rho = ts.transverseCurvature();
  //
  // Straight line approximation? |rho|<1.e-10 equivalent to ~ 1um 
  // difference in transversal position at 10m.
  //
  if  unlikely( fabs(rho)<1.e-10 ) {
    //
    // Instantiate auxiliary object for finding intersection.
    // Frame-independant point and vector are created explicitely to 
    // avoid confusing gcc (refuses to compile with temporary objects
    // in the constructor).
    //
    LogDebug("RKPropagatorInS")<<" startZ = "<<startZ;

    if unlikely(fabs(startZ) < 1e-5){  
      LogDebug("RKPropagatorInS")<< "Propagation is not performed: state is already on final surface.";
      GlobalTrajectoryParameters res( gpos, gmom, ts.charge(), 
                      theVolume);
      return GlobalParametersWithPath( res, 0.0);
    }

    StraightLinePlaneCrossing::PositionType pos(gpos);
    StraightLinePlaneCrossing::DirectionType dir(gmom);
    StraightLinePlaneCrossing planeCrossing(pos,dir, propagationDirection());
    //
    // get solution
    //
    std::pair<bool,double> propResult = planeCrossing.pathLength(plane);
    if likely( propResult.first && theVolume !=nullptr) {
      double s = propResult.second;
      // point (reconverted to GlobalPoint)
      GlobalPoint x (planeCrossing.position(s));
      GlobalTrajectoryParameters res( x, gmom, ts.charge(), theVolume);
      return GlobalParametersWithPath( res, s);
      }
    //do someting
    LogDebug("RKPropagatorInS")<< "Straight line propgation to plane failed !!"; 
    return GlobalParametersWithPath( );

  }


#ifdef EDM_ML_DEBUG
  if (theVolume != 0) {
    LogDebug("RKPropagatorInS")  << "RKPropagatorInS: starting prop to plane in volume with pos " << theVolume->position()
          << " Z axis " << theVolume->toGlobal( LocalVector(0,0,1)) ;

    LogDebug("RKPropagatorInS")  << "The starting position is " << ts.position() << " (global) "
          << theVolume->toLocal(ts.position()) << " (local) " ;
  
    FrameChanger changer;
    auto localPlane = changer.transformPlane( plane, *theVolume);
    LogDebug("RKPropagatorInS")  << "The plane position is " << plane.position() << " (global) "
          << localPlane.position() << " (local) " ;

    LogDebug("RKPropagatorInS")  << "The initial distance to plane is " << plane.localZ( ts.position()) ;

    StraightLinePlaneCrossing cross( ts.position().basicVector(), ts.momentum().basicVector());
    std::pair<bool,double> res3 = cross.pathLength(plane);
    LogDebug("RKPropagatorInS")  << "straight line distance " << res3.first << " " << res3.second ;
  }
#endif

  typedef RKAdaptiveSolver<double,RKOneCashKarpStep, 6>   Solver;
  typedef Solver::Vector                                  RKVector;
  
  RKLocalFieldProvider field( fieldProvider());
  PathToPlane2Order pathLength( field, &field.frame());
  CartesianLorentzForce deriv(field, ts.charge());

  RKCartesianDistance dist;
  double eps = theTolerance;
  Solver solver;
  double stot = 0;
  PropagationDirection currentDirection = propagationDirection();

  // in magVolume frame
  RKVector start( CartesianStateAdaptor::rkstate( rkPosition(gpos), rkMomentum(gmom)));
  int safeGuard = 0;
  while (safeGuard++<100) {
    CartesianStateAdaptor startState(start);

    std::pair<bool,double> path = pathLength( plane, startState.position(), 
                          startState.momentum(), 
                          (double) ts.charge(), currentDirection);
    if  unlikely(!path.first) { 
      LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path length calculation to plane failed!" 
                   << "...distance to plane " << plane.localZ( globalPosition(startState.position()))
                   << "...Local starting position in volume " << startState.position() 
                   << "...Magnetic field " << field.inTesla( startState.position()) ;


      return GlobalParametersWithPath();
    }

    LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path lenght to plane is " << path.second ;
   

    double sstep = path.second;
    if  unlikely( std::abs(sstep) < eps) {
      LogDebug("RKPropagatorInS")  << "On-surface accuracy not reached, but pathLength calculation says we are there! "
         << "path " << path.second << " distance to plane is " << startZ ;
      GlobalTrajectoryParameters res( gtpFromVolumeLocal( startState, ts.charge()));
      return GlobalParametersWithPath( res, stot);
    }

    LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Solving for " << sstep 
          << " current distance to plane is " << startZ ;

    RKVector rkresult = solver( 0, start, sstep, deriv, dist, eps);
    stot += sstep;
    CartesianStateAdaptor cur( rkresult);
    double remainingZ = plane.localZ( globalPosition(cur.position()));

    if ( fabs(remainingZ) < eps) {
      LogDebug("RKPropagatorInS")  << "On-surface accuracy reached! " << remainingZ ;
      GlobalTrajectoryParameters res( gtpFromVolumeLocal( cur, ts.charge()));
      return GlobalParametersWithPath( res, stot);
    }

    start = rkresult;

    if (remainingZ * startZ > 0) {
      LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in same direction again " 
        << remainingZ ;
    }
    else {
      LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in opposite direction " 
        << remainingZ ;
      currentDirection = invertDirection( currentDirection);
    }
    startZ = remainingZ;
  }

  edm::LogError("FailedPropagation") << " too many iterations trying to reach plane ";
  return GlobalParametersWithPath();
}

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInS::propagateWithPath ( const FreeTrajectoryState &  fts, const Plane &  plane  ) const

overrideprivatevirtual

Implements Propagator.

Definition at line 27 of file RKPropagatorInS.cc.

References SurfaceSideDefinition::afterSurface, alongMomentum, analyticalErrorPropagation(), SurfaceSideDefinition::beforeSurface, runTauDisplay::gp, GlobalParametersWithPath::parameters(), propagateParametersOnPlane(), Propagator::propagationDirection(), GlobalParametersWithPath::s(), and unlikely.

Referenced by ~RKPropagatorInS().

{
  GlobalParametersWithPath gp =  propagateParametersOnPlane( fts, plane);
  if unlikely(!gp) return TsosWP(TrajectoryStateOnSurface(),0.);

  SurfaceSideDefinition::SurfaceSide side = PropagationDirectionFromPath()(gp.s(),propagationDirection())==alongMomentum 
    ? SurfaceSideDefinition::beforeSurface : SurfaceSideDefinition::afterSurface;
  return  analyticalErrorPropagation( fts, plane, side, gp.parameters(),gp.s());
}

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInS::propagateWithPath ( const FreeTrajectoryState &  fts, const Cylinder &  cyl  ) const

overrideprivatevirtual

Implements Propagator.

Definition at line 39 of file RKPropagatorInS.cc.

References SurfaceSideDefinition::afterSurface, alongMomentum, analyticalErrorPropagation(), SurfaceSideDefinition::beforeSurface, runTauDisplay::gp, GlobalParametersWithPath::parameters(), propagateParametersOnCylinder(), Propagator::propagationDirection(), GlobalParametersWithPath::s(), and unlikely.

{
  GlobalParametersWithPath gp =  propagateParametersOnCylinder( fts, cyl);
  if  unlikely(!gp) return TsosWP(TrajectoryStateOnSurface(),0.);

  SurfaceSideDefinition::SurfaceSide side = PropagationDirectionFromPath()(gp.s(),propagationDirection())==alongMomentum 
    ? SurfaceSideDefinition::beforeSurface : SurfaceSideDefinition::afterSurface;
  return analyticalErrorPropagation( fts, cyl, side, gp.parameters(),gp.s());
  
}

Basic3DVector< double > RKPropagatorInS::rkMomentum ( const GlobalVector &  mom ) const

private

Definition at line 367 of file RKPropagatorInS.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), theVolume, and GloballyPositioned< T >::toLocal().

Referenced by propagateParametersOnPlane().

{
  if (theVolume != nullptr) return theVolume->toLocal( mom).basicVector();
  else return mom.basicVector();
}

Basic3DVector< double > RKPropagatorInS::rkPosition ( const GlobalPoint &  pos ) const

private

Definition at line 361 of file RKPropagatorInS.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), theVolume, and GloballyPositioned< T >::toLocal().

Referenced by propagateParametersOnPlane().

{
  if (theVolume != nullptr) return theVolume->toLocal( pos).basicVector();
  else return pos.basicVector();
}

Member Data Documentation

double RKPropagatorInS::theTolerance

private

Definition at line 55 of file RKPropagatorInS.h.

Referenced by propagateParametersOnCylinder(), and propagateParametersOnPlane().

const MagVolume* RKPropagatorInS::theVolume

private

Definition at line 54 of file RKPropagatorInS.h.

Referenced by fieldProvider(), globalMomentum(), globalPosition(), gtpFromLocal(), gtpFromVolumeLocal(), magneticField(), propagateParametersOnCylinder(), propagateParametersOnPlane(), rkMomentum(), and rkPosition().