#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.

References clone(), myPropagate(), and propagateWithPath().

Referenced by clone().

10 :

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

Propagator::Propagator

Propagator(PropagationDirection dir=alongMomentum)

Definition: Propagator.h:46

RKPropagatorInR::theVolume

const MagVolume * theVolume

Definition: RKPropagatorInR.h:31

dir

dbl *** dir

Definition: mlp_gen.cc:35

Member Function Documentation

Propagator * RKPropagatorInR::clone ( void ) const

virtual

Implements Propagator.

Definition at line 73 of file RKPropagatorInR.cc.

References RKPropagatorInR().

Referenced by RKPropagatorInR().

 {
     return new RKPropagatorInR(*this);
 }

virtual const MagneticField* RKPropagatorInR::magneticField ( ) const

inlinevirtual

Implements Propagator.

Definition at line 27 of file RKPropagatorInR.h.

References theVolume.

27 {return theVolume;}

RKPropagatorInR::theVolume

const MagVolume * theVolume

Definition: RKPropagatorInR.h:31

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

Definition at line 56 of file RKPropagatorInR.cc.

Referenced by RKPropagatorInR().

 {
     return TrajectoryStateOnSurface();
 }

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

Definition at line 15 of file RKPropagatorInR.cc.

References FreeTrajectoryState::charge(), MillePedeFileConverter_cfg::e, FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), Cylinder::radius(), command_line::start, 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.

Referenced by RKPropagatorInR().

 {
     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().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation