RKPropagatorInR Class Reference

#include <RKPropagatorInR.h>

Inheritance diagram for RKPropagatorInR:

Public Member Functions
virtual Propagator *	clone () const
virtual const MagneticField *	magneticField () const
TrajectoryStateOnSurface	myPropagate (const FreeTrajectoryState &, const Plane &) const
TrajectoryStateOnSurface	myPropagate (const FreeTrajectoryState &, const Cylinder &) const
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Plane &) const
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Cylinder &) const
	RKPropagatorInR (const MagVolume &vol, PropagationDirection dir=alongMomentum)
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 Attributes
const MagVolume *	theVolume

Detailed Description

Definition at line 7 of file RKPropagatorInR.h.

Constructor & Destructor Documentation

RKPropagatorInR::RKPropagatorInR ( const MagVolume &  vol, PropagationDirection  dir = alongMomentum  )

inline

Definition at line 10 of file RKPropagatorInR.h.

Referenced by clone().

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

Member Function Documentation

Propagator * RKPropagatorInR::clone ( void  ) const

virtual

Implements Propagator.

Definition at line 73 of file RKPropagatorInR.cc.

References RKPropagatorInR().

{
    return new RKPropagatorInR(*this);
}

virtual const MagneticField* RKPropagatorInR::magneticField ( ) const

inlinevirtual

Implements Propagator.

Definition at line 27 of file RKPropagatorInR.h.

References theVolume.

{return theVolume;}

TrajectoryStateOnSurface RKPropagatorInR::myPropagate	(	const FreeTrajectoryState &	,
		const Plane &
	)		const

Definition at line 56 of file RKPropagatorInR.cc.

{
    return TrajectoryStateOnSurface();
}

TrajectoryStateOnSurface RKPropagatorInR::myPropagate	(	const FreeTrajectoryState &	ts,
		const Cylinder &	cyl
	)		const

Definition at line 15 of file RKPropagatorInR.cc.

References FreeTrajectoryState::charge(), alignCSCRings::e, FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), Cylinder::radius(), dqm_diff::start, relval_parameters_module::step, theVolume, Surface::toGlobal(), and GloballyPositioned< T >::toLocal().

{
  //typedef RK4PreciseSolver<double,5>           Solver;
    typedef RKAdaptiveSolver<double,RKOne4OrderStep, 5>     Solver;
    //typedef RKAdaptiveSolver<Scalar,RKOneCashKarpStep, 5>   Solver;
    typedef double                                          Scalar;
    typedef Solver::Vector                                  RKVector;

    GlobalPoint pos( ts.position());
    GlobalVector mom( ts.momentum());

    LocalPoint startpos = cyl.toLocal(pos);
    LocalVector startmom = cyl.toLocal(mom);

    CylindricalState startState( startpos, startmom, ts.charge());
    RKVector start = startState.parameters();

    RKLocalFieldProvider localField( *theVolume, cyl);
    CylindricalLorentzForce<double,5> deriv(localField);
    RKCylindricalDistance<double,5> dist;
    double eps = 1.e-5;
    Solver solver;
    try {
    Scalar step = cyl.radius() - startState.rho();
    RKVector rkresult = solver( startState.rho(), start, step, deriv, dist, eps);
    CylindricalState endState( cyl.radius(), rkresult, startState.prSign());
    return TrajectoryStateOnSurface( GlobalTrajectoryParameters( cyl.toGlobal( endState.position()), 
                                     cyl.toGlobal( endState.momentum()),
                                     TrackCharge(endState.charge()), 
                                     theVolume),
                     cyl);
    }
    catch (CylindricalLorentzForceException& e) {
        // the propagation failed due to momentum almost parallel to the plane.
        // This does not mean the propagation is impossible, but it should be done
    // in a different parametrization (e.g. s)  
    return TrajectoryStateOnSurface();
    }
}

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInR::propagateWithPath ( const FreeTrajectoryState &  , const Plane &    ) const

virtual

Implements Propagator.

Definition at line 62 of file RKPropagatorInR.cc.

{
    return std::pair< TrajectoryStateOnSurface, double>();
}

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInR::propagateWithPath ( const FreeTrajectoryState &  , const Cylinder &    ) const

virtual

Implements Propagator.

Definition at line 68 of file RKPropagatorInR.cc.

{
    return std::pair< TrajectoryStateOnSurface, double>();
}

Member Data Documentation

const MagVolume* RKPropagatorInR::theVolume

private

Definition at line 31 of file RKPropagatorInR.h.

Referenced by magneticField(), and myPropagate().