#include <RKPropagatorInZ.h>

Inheritance diagram for RKPropagatorInZ:

Public Member Functions
Propagator *	clone () const override

const MagneticField *	magneticField () const override

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

TrajectoryStateOnSurface	myPropagate (const FreeTrajectoryState &, const Cylinder &) const

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

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

	RKPropagatorInZ (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 RKPropagatorInZ.h.

Constructor & Destructor Documentation

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

inline

Definition at line 10 of file RKPropagatorInZ.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

RKPropagatorInZ::theVolume

const MagVolume * theVolume

Definition: RKPropagatorInZ.h:31

dir

dbl *** dir

Definition: mlp_gen.cc:35

Member Function Documentation

Propagator * RKPropagatorInZ::clone ( void ) const

overridevirtual

Implements Propagator.

Definition at line 79 of file RKPropagatorInZ.cc.

References RKPropagatorInZ().

Referenced by RKPropagatorInZ().

 {
     return new RKPropagatorInZ(*this);
 }

const MagneticField* RKPropagatorInZ::magneticField ( ) const

inlineoverridevirtual

Implements Propagator.

Definition at line 27 of file RKPropagatorInZ.h.

References theVolume.

27 {return theVolume;}

RKPropagatorInZ::theVolume

const MagVolume * theVolume

Definition: RKPropagatorInZ.h:31

TrajectoryStateOnSurface RKPropagatorInZ::myPropagate	(	const FreeTrajectoryState &	ts,
		const Plane &	plane
	)		const

Definition at line 13 of file RKPropagatorInZ.cc.

References FreeTrajectoryState::charge(), MillePedeFileConverter_cfg::e, PV3DBase< T, PVType, FrameType >::mag(), FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), command_line::start, theVolume, GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by RKPropagatorInZ().

 {
   //typedef RK4PreciseSolver<double,5>           Solver;
   //typedef RKAdaptiveSolver<double,RKOne4OrderStep, 5>   Solver;
     typedef RKAdaptiveSolver<double,RKOneCashKarpStep, 5>   Solver;
     typedef Solver::Vector                       RKVector;
 
     GlobalPoint pos( ts.position());
     GlobalVector mom( ts.momentum());
 
     // cout << "RKPropagatorInZ: starting from FTS " << ts << endl;
 
     LocalPoint startpos = plane.toLocal(pos);
     LocalVector startmom = plane.toLocal(mom);
     double pzSign = startmom.z() > 0 ? 1.0 : -1.0;
 
     // cout << "In local plane coordinates: " << startpos << ", momentum " << startmom << endl;
 
     RKVector start;
     start(0) = startpos.x();
     start(1) = startpos.y();
     start(2) = startmom.x()/startmom.z();
     start(3) = startmom.y()/startmom.z();
     start(4) = pzSign * ts.charge() / startmom.mag();
 
     // cout << "RKPropagatorInZ: Solving with par " <<  startpos.z() << " and state " << start << endl;
 
     RKLocalFieldProvider localField( *theVolume, plane);
 
     CurvilinearLorentzForce<double,5> deriv(localField);
     RKCurvilinearDistance<double,5> dist;
     double eps = 1.e-5;
     Solver solver;
     try {
     RKVector rkresult = solver( startpos.z(), start, -startpos.z(), deriv, dist, eps);
 
     return TrajectoryStateOnSurface( LocalTrajectoryParameters( rkresult(4), rkresult(2), rkresult(3),
                                     rkresult(0), rkresult(1), pzSign),
                      plane, theVolume);
     }
     catch (CurvilinearLorentzForceException& 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();
     }
 }

TrajectoryStateOnSurface RKPropagatorInZ::myPropagate	(	const FreeTrajectoryState &	,
		const Cylinder &
	)		const

Definition at line 62 of file RKPropagatorInZ.cc.

 {
     return TrajectoryStateOnSurface();
 }

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

overridevirtual

Implements Propagator.

Definition at line 68 of file RKPropagatorInZ.cc.

Referenced by RKPropagatorInZ().

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

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

overridevirtual

Implements Propagator.

Definition at line 74 of file RKPropagatorInZ.cc.

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

Member Data Documentation

const MagVolume* RKPropagatorInZ::theVolume

private

Definition at line 31 of file RKPropagatorInZ.h.

Referenced by magneticField(), and myPropagate().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation