#include <StraightLinePropagator.h>

Inheritance diagram for StraightLinePropagator:

Public Member Functions
virtual StraightLinePropagator *	clone () const override

virtual const MagneticField *	magneticField () const override

std::pair< TSOS, double >	propagateWithPath (const FreeTrajectoryState &fts, const Plane &surface) const override

std::pair< TSOS, double >	propagateWithPath (const FreeTrajectoryState &fts, const Cylinder &surface) const override

	StraightLinePropagator (const MagneticField *field, PropagationDirection aDir=alongMomentum)

	~StraightLinePropagator ()

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 FreeTrajectoryState	FTS

typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
AlgebraicMatrix55	jacobian (double &s) const

TrajectoryStateOnSurface	propagatedState (const FreeTrajectoryState &fts, const Surface &surface, const AlgebraicMatrix55 &jacobian, const GlobalPoint &x, const GlobalVector &p) const

TrajectoryStateOnSurface	propagatedState (const FreeTrajectoryState &fts, const Surface &surface, const AlgebraicMatrix55 &jacobian, const LocalPoint &x, const LocalVector &p) const

bool	propagateParametersOnCylinder (const FreeTrajectoryState &fts, const Cylinder &cylinder, GlobalPoint &x, GlobalVector &p, double &s) const

bool	propagateParametersOnPlane (const FreeTrajectoryState &fts, const Plane &plane, LocalPoint &x, LocalVector &p, double &s) const

Private Attributes
const MagneticField *	theField

Detailed Description

As the name indicates, propagates track parameters according to a straight line model. Intended for test beams without magnetic field and similar cases.

Warning: The errors are NOT propagated.

Definition at line 17 of file StraightLinePropagator.h.

Member Typedef Documentation

typedef FreeTrajectoryState StraightLinePropagator::FTS

private

Definition at line 21 of file StraightLinePropagator.h.

typedef TrajectoryStateOnSurface StraightLinePropagator::TSOS

private

Definition at line 22 of file StraightLinePropagator.h.

Constructor & Destructor Documentation

StraightLinePropagator::StraightLinePropagator	(	const MagneticField *	field,
		PropagationDirection	aDir = `alongMomentum`
	)

inline

Definition at line 26 of file StraightLinePropagator.h.

Referenced by clone().

27 :

28 Propagator(aDir), theField(field) {}

Propagator::Propagator

Propagator(PropagationDirection dir=alongMomentum)

Definition: Propagator.h:46

StraightLinePropagator::theField

const MagneticField * theField

Definition: StraightLinePropagator.h:51

StraightLinePropagator::~StraightLinePropagator ( )

inline

Definition at line 30 of file StraightLinePropagator.h.

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

30 {}

Member Function Documentation

virtual StraightLinePropagator* StraightLinePropagator::clone ( void ) const

inlineoverridevirtual

Implements Propagator.

Definition at line 42 of file StraightLinePropagator.h.

References StraightLinePropagator().

                                                           {
     return new StraightLinePropagator(*this);
   }

AlgebraicMatrix55 StraightLinePropagator::jacobian ( double & s ) const

private

Definition at line 93 of file StraightLinePropagator.cc.

References alongMomentum, dir, Propagator::propagationDirection(), and alignCSCRings::s.

Referenced by propagateWithPath().

                                                                   {
   //Jacobian for 5*5 local error matrix
   AlgebraicMatrix55 j = AlgebraicMatrixID(); //Jacobian
   
   double dir = (propagationDirection() == alongMomentum) ? 1. : -1.;
   if (s*dir < 0.) return j;
 
   j(3,1) = s; 
   j(4,2) = s; 
 
   return j;
 }

virtual const MagneticField* StraightLinePropagator::magneticField ( ) const

inlineoverridevirtual

Implements Propagator.

Definition at line 47 of file StraightLinePropagator.h.

References theField.

47 {return theField;}

StraightLinePropagator::theField

const MagneticField * theField

Definition: StraightLinePropagator.h:51

TrajectoryStateOnSurface StraightLinePropagator::propagatedState	(	const FreeTrajectoryState &	fts,
		const Surface &	surface,
		const AlgebraicMatrix55 &	jacobian,
		const GlobalPoint &	x,
		const GlobalVector &	p
	)		const

private

Definition at line 71 of file StraightLinePropagator.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::hasError(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LocalTrajectoryError::matrix(), theField, and tmp.

Referenced by propagateWithPath().

                                                  {
 
   if(fts.hasError()) {
     // propagate error
     TSOS tmp(fts, surface);
     const AlgebraicSymMatrix55 & eLocal =tmp.localError().matrix();
     AlgebraicSymMatrix55 lte = ROOT::Math::Similarity(jacobian,eLocal);
     LocalTrajectoryError eloc(lte);
 
     TSOS tmp2(tmp.localParameters(), eloc, surface, theField);
     GlobalTrajectoryParameters gtp(x, p, fts.charge(), theField);
     return TSOS(gtp, tmp2.cartesianError(), surface);
   } else {
     // return state without errors
     return TSOS(GlobalTrajectoryParameters(x, p, fts.charge(), theField), surface);
   }
 }

TrajectoryStateOnSurface StraightLinePropagator::propagatedState	(	const FreeTrajectoryState &	fts,
		const Surface &	surface,
		const AlgebraicMatrix55 &	jacobian,
		const LocalPoint &	x,
		const LocalVector &	p
	)		const

private

Definition at line 51 of file StraightLinePropagator.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::hasError(), TrajectoryStateOnSurface::localError(), LocalTrajectoryError::matrix(), theField, and tmp.

                                                 {
   if(fts.hasError()) {
     // propagate error
     TSOS tmp( fts, surface);
     const AlgebraicSymMatrix55 & eLocal =tmp.localError().matrix();
     AlgebraicSymMatrix55 lte = ROOT::Math::Similarity(jacobian,eLocal);
     LocalTrajectoryError eloc(lte);
     LocalTrajectoryParameters ltp(x, p, fts.charge());
     return TSOS(ltp, eloc, surface, theField);
   } else {
     // return state without errors
     return TSOS(LocalTrajectoryParameters(x, p, fts.charge()), surface, theField);
   }
 }

bool StraightLinePropagator::propagateParametersOnCylinder	(	const FreeTrajectoryState &	fts,
		const Cylinder &	cylinder,
		GlobalPoint &	x,
		GlobalVector &	p,
		double &	s
	)		const

private

Definition at line 106 of file StraightLinePropagator.cc.

References alongMomentum, dir, dt, FreeTrajectoryState::momentum(), PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), Propagator::propagationDirection(), Cylinder::radius(), Surface::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by propagateWithPath().

                                                 {
   GlobalPoint sp = cylinder.toGlobal(LocalPoint(0., 0.));
   if (sp.x()!=0. || sp.y()!=0.) {
     throw PropagationException("Cannot propagate to an arbitrary cylinder");
   }
 
   x = fts.position();
   p = fts.momentum();
   s = cylinder.radius() - x.perp();
 
   double dir = (propagationDirection() == alongMomentum) ? 1. : -1.;
   if(s*dir < 0.) return false;
 
   double dt = s/p.perp();
   x = GlobalPoint(x.x() + p.x()*dt, 
                   x.y() + p.y()*dt, 
                   x.z() + p.z()*dt);
 
   return true;
 }

bool StraightLinePropagator::propagateParametersOnPlane	(	const FreeTrajectoryState &	fts,
		const Plane &	plane,
		LocalPoint &	x,
		LocalVector &	p,
		double &	s
	)		const

private

Definition at line 131 of file StraightLinePropagator.cc.

References FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), alignCSCRings::s, GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by propagateWithPath().

                                              {
   
   //Do extrapolation in local frame of plane
   //  LocalPoint sp = plane.toLocal(plane.position());
   x = plane.toLocal(fts.position());
   p = plane.toLocal(fts.momentum());
   s = -x.z(); // sp.z() - x.z(); local z of plane always 0
   
   //double dir = (propagationDirection() == alongMomentum) ? 1. : -1.;
   //if(s*dir < 0.) return false;
   if ((p.x() != 0 || p.y() != 0) && p.z() == 0 && s!= 0) return false;
 
   x = LocalPoint( x.x() + (p.x()/p.z())*s,
                   x.y() + (p.y()/p.z())*s,
                   x.z() + s);    
 
   return true;
 }

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

overridevirtual

Implements Propagator.

Definition at line 9 of file StraightLinePropagator.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::hasError(), jacobian(), AlCaHLTBitMon_ParallelJobs::p, propagatedState(), propagateParametersOnPlane(), alignCSCRings::s, theField, and x.

Referenced by ~StraightLinePropagator().

 {
   // propagate parameters
   LocalPoint x;
   LocalVector p;
   double s = 0;
   bool parametersOK = propagateParametersOnPlane(fts, plane, x, p, s);
   if(!parametersOK) return std::make_pair(TrajectoryStateOnSurface(), 0.);
 
   // compute propagated state
   if (fts.hasError()) {
     return std::make_pair( propagatedState(fts, plane, jacobian(s), x,  p), s);
   } else {
     // return state without errors
     return std::make_pair(TSOS(LocalTrajectoryParameters(x, p, fts.charge()), 
               plane, theField), s);
   }
 }

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

overridevirtual

Implements Propagator.

Definition at line 30 of file StraightLinePropagator.cc.

References FreeTrajectoryState::charge(), FreeTrajectoryState::hasError(), jacobian(), AlCaHLTBitMon_ParallelJobs::p, propagatedState(), propagateParametersOnCylinder(), alignCSCRings::s, theField, and x.

 {
   // propagate parameters
   GlobalPoint x;
   GlobalVector p;
   double s = 0;
   bool parametersOK = propagateParametersOnCylinder(fts, cylinder,  x, p, s);
   if(!parametersOK) return std::make_pair(TrajectoryStateOnSurface(), 0.);
 
   // compute propagated state
   if (fts.hasError()) {
     return std::make_pair( propagatedState(fts, cylinder, jacobian(s), x,  p), s);
   } else {
     // return state without errors
     return std::make_pair(TSOS(GlobalTrajectoryParameters(x, p, fts.charge(),theField),
               cylinder), s);
   }
 }

Member Data Documentation

const MagneticField* StraightLinePropagator::theField

private

Definition at line 51 of file StraightLinePropagator.h.

Referenced by magneticField(), propagatedState(), and propagateWithPath().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation