#include <SmartPropagator.h>

Inheritance diagram for SmartPropagator:

Public Member Functions
SmartPropagator *	clone () const override
	Virtual constructor (using copy c'tor) More...

const Propagator *	getGenPropagator () const
	return the propagator used outside tracker More...

const Propagator *	getTkPropagator () const
	return the propagator used inside tracker More...

bool	insideTkVol (const FreeTrajectoryState &fts) const
	true if a fts is inside tracker volume More...

bool	insideTkVol (const Surface &surface) const
	true if a surface is inside tracker volume More...

bool	insideTkVol (const Cylinder &cylin) const
	true if a cylinder is inside tracker volume More...

bool	insideTkVol (const Plane &plane) const
	true if a plane is inside tracker volume More...

const MagneticField *	magneticField () const override
	return the magneticField More...

void	setPropagationDirection (PropagationDirection dir) override
	setting the direction fo both components More...

	SmartPropagator (const Propagator aTkProp, const Propagator aGenProp, const MagneticField *field, PropagationDirection dir=alongMomentum, float epsilon=5)
	Defines which propagator is used inside Tk and which outside. More...

	SmartPropagator (const Propagator &aTkProp, const Propagator &aGenProp, const MagneticField *field, PropagationDirection dir=alongMomentum, float epsilon=5)
	Defines which propagator is used inside Tk and which outside. More...

	SmartPropagator (const SmartPropagator &)
	Copy constructor. More...

	~SmartPropagator () 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 < 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	~Propagator ()

Private Member Functions
void	initTkVolume (float epsilon)
	build the tracker volume More...

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

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

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

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

Private Attributes
const MagneticField *	theField

Propagator *	theGenProp

Propagator *	theTkProp

ReferenceCountingPointer < Cylinder >	theTkVolume

Detailed Description

A propagator which use different algorithm to propagate inside or outside tracker

Author: Stefano Lacaprara - INFN Padova author Chang Liu - Purdue University

Modification: 26-Jun-2002 SL: theTkVolume is now a static ReferenceCountingPointer<BoundCylinder> 28-Aug-2002 SL: added methods to unhide Propagator methods 29-Oct-2002 SL: fixed clone and copy constructor, and BoundCylinder are build with CylinderBuilder to enforce the referencePointer

Definition at line 34 of file SmartPropagator.h.

Constructor & Destructor Documentation

SmartPropagator::SmartPropagator	(	const Propagator *	aTkProp,
		const Propagator *	aGenProp,
		const MagneticField *	field,
		PropagationDirection	dir = `alongMomentum`,
		float	epsilon = `5`
	)

Defines which propagator is used inside Tk and which outside.

Definition at line 36 of file SmartPropagator.cc.

References initTkVolume().

Referenced by clone().

     : Propagator(dir), theTkProp(aTkProp->clone()), theGenProp(aGenProp->clone()), theField(field) {
   initTkVolume(epsilon);
 }

SmartPropagator::SmartPropagator	(	const Propagator &	aTkProp,
		const Propagator &	aGenProp,
		const MagneticField *	field,
		PropagationDirection	dir = `alongMomentum`,
		float	epsilon = `5`
	)

Defines which propagator is used inside Tk and which outside.

Definition at line 45 of file SmartPropagator.cc.

References initTkVolume().

     : Propagator(dir), theTkProp(aTkProp.clone()), theGenProp(aGenProp.clone()), theField(field) {
   initTkVolume(epsilon);
 }

SmartPropagator::SmartPropagator ( const SmartPropagator & aProp )

Copy constructor.

Definition at line 54 of file SmartPropagator.cc.

References Propagator::clone(), getGenPropagator(), getTkPropagator(), theGenProp, and theTkProp.

     : Propagator(aProp.propagationDirection()), theTkProp(nullptr), theGenProp(nullptr) {
   if (aProp.theTkProp)
     theTkProp = aProp.getTkPropagator()->clone();
   if (aProp.theGenProp)
     theTkProp = aProp.getGenPropagator()->clone();
 
   //SL since it's a copy constructor, then the TkVolume has been already
   //initialized
   // initTkVolume(epsilon);
 }

SmartPropagator::~SmartPropagator ( )

override

virtual destructor

Definition at line 67 of file SmartPropagator.cc.

References theGenProp, and theTkProp.

                                   {
   delete theTkProp;
   delete theGenProp;
 }

Member Function Documentation

SmartPropagator* SmartPropagator::clone ( void ) const

inlineoverridevirtual

Virtual constructor (using copy c'tor)

Implements Propagator.

Definition at line 58 of file SmartPropagator.h.

References getGenPropagator(), getTkPropagator(), magneticField(), and SmartPropagator().

                                           {
     return new SmartPropagator(getTkPropagator(), getGenPropagator(), magneticField());
   }

const Propagator * SmartPropagator::getGenPropagator ( ) const

return the propagator used outside tracker

Definition at line 157 of file SmartPropagator.cc.

References theGenProp.

Referenced by clone(), propagateWithPath(), and SmartPropagator().

157 { return theGenProp; }

SmartPropagator::theGenProp

Propagator * theGenProp

Definition: SmartPropagator.h:107

const Propagator * SmartPropagator::getTkPropagator ( ) const

return the propagator used inside tracker

Definition at line 155 of file SmartPropagator.cc.

References theTkProp.

Referenced by clone(), propagateWithPath(), and SmartPropagator().

155 { return theTkProp; }

SmartPropagator::theTkProp

Propagator * theTkProp

Definition: SmartPropagator.h:106

void SmartPropagator::initTkVolume ( float epsilon )

private

build the tracker volume

Definition at line 73 of file SmartPropagator.cc.

References Cylinder::build(), geometryDiff::epsilon, TrackerBounds::halfLength(), TrackerBounds::radius(), CosmicsPD_Skims::radius, makeMuonMisalignmentScenario::rot, and theTkVolume.

Referenced by SmartPropagator().

                                                 {
   //
   // fill tracker dimensions
   //
   float radius = TrackerBounds::radius();
   float r_out = radius + epsilon / 2;
   float r_in = r_out - epsilon;
   float z_max = TrackerBounds::halfLength();
   float z_min = -z_max;
 
   Surface::PositionType pos(0, 0, 0);  // centered at the global origin
   Surface::RotationType rot;           // unit matrix - barrel cylinder orientation
 
   theTkVolume = Cylinder::build(radius, pos, rot, new SimpleCylinderBounds(r_in, r_out, z_min, z_max));
 }

bool SmartPropagator::insideTkVol ( const FreeTrajectoryState & fts ) const

true if a fts is inside tracker volume

Definition at line 125 of file SmartPropagator.cc.

References runTauDisplay::gp, TrackerBounds::halfLength(), PV3DBase< T, PVType, FrameType >::perp(), FreeTrajectoryState::position(), TrackerBounds::radius(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by propagateWithPath().

                                                                       {
   GlobalPoint gp = fts.position();
   //  LocalPoint lp = theTkVolume()->toLocal(gp);
   //  return theTkVolume()->bounds().inside(lp);
 
   return ((gp.perp() <= TrackerBounds::radius() + 10.) && (fabs(gp.z()) <= TrackerBounds::halfLength() + 10.));
 }

bool SmartPropagator::insideTkVol ( const Surface & surface ) const

true if a surface is inside tracker volume

Definition at line 133 of file SmartPropagator.cc.

References runTauDisplay::gp, TrackerBounds::halfLength(), PV3DBase< T, PVType, FrameType >::perp(), GloballyPositioned< T >::position(), TrackerBounds::radius(), and PV3DBase< T, PVType, FrameType >::z().

                                                               {
   const GlobalPoint& gp = surface.position();
   // LocalPoint lp = theTkVolume()->toLocal(gp);
 
   // return theTkVolume()->bounds().inside(lp);
   return ((gp.perp() <= TrackerBounds::radius() + 10.) && (fabs(gp.z()) <= TrackerBounds::halfLength() + 10.));
 }

bool SmartPropagator::insideTkVol ( const Cylinder & cylin ) const

true if a cylinder is inside tracker volume

bool SmartPropagator::insideTkVol ( const Plane & plane ) const

true if a plane is inside tracker volume

Definition at line 148 of file SmartPropagator.cc.

References runTauDisplay::gp, TrackerBounds::halfLength(), PV3DBase< T, PVType, FrameType >::perp(), GloballyPositioned< T >::position(), TrackerBounds::radius(), and PV3DBase< T, PVType, FrameType >::z().

                                                           {
   const GlobalPoint& gp = plane.position();
   //  LocalPoint lp = theTkVolume()->toLocal(gp);
   //  return theTkVolume()->bounds().inside(lp);
   return ((gp.perp() <= TrackerBounds::radius() + 10.) && (fabs(gp.z()) <= TrackerBounds::halfLength() + 10.));
 }

const MagneticField* SmartPropagator::magneticField ( ) const

inlineoverridevirtual

return the magneticField

Implements Propagator.

Definition at line 100 of file SmartPropagator.h.

References theField.

Referenced by clone().

100 { return theField; }

SmartPropagator::theField

const MagneticField * theField

Definition: SmartPropagator.h:108

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

overrideprivatevirtual

Implements Propagator.

Definition at line 89 of file SmartPropagator.cc.

References getGenPropagator(), getTkPropagator(), insideTkVol(), and Propagator::propagateWithPath().

                                                                                                          {
   if (insideTkVol(fts) && insideTkVol(plane)) {
     return getTkPropagator()->propagateWithPath(fts, plane);
   } else {
     return getGenPropagator()->propagateWithPath(fts, plane);
   }
 }

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

overrideprivatevirtual

Implements Propagator.

Definition at line 98 of file SmartPropagator.cc.

References getGenPropagator(), getTkPropagator(), insideTkVol(), and Propagator::propagateWithPath().

                                                                                                                {
   if (insideTkVol(fts) && insideTkVol(cylinder)) {
     return getTkPropagator()->propagateWithPath(fts, cylinder);
   } else {
     return getGenPropagator()->propagateWithPath(fts, cylinder);
   }
 }

std::pair< TrajectoryStateOnSurface, double > SmartPropagator::propagateWithPath	(	const TrajectoryStateOnSurface &	tsos,
		const Plane &	sur
	)		const

overrideprivatevirtual

Reimplemented from Propagator.

Definition at line 107 of file SmartPropagator.cc.

References TrajectoryStateOnSurface::freeState(), getGenPropagator(), getTkPropagator(), insideTkVol(), and Propagator::propagateWithPath().

                                                                                                          {
   if (insideTkVol(*fts.freeState()) && insideTkVol(plane)) {
     return getTkPropagator()->propagateWithPath(fts, plane);
   } else {
     return getGenPropagator()->propagateWithPath(fts, plane);
   }
 }

std::pair< TrajectoryStateOnSurface, double > SmartPropagator::propagateWithPath	(	const TrajectoryStateOnSurface &	tsos,
		const Cylinder &	sur
	)		const

overrideprivatevirtual

Reimplemented from Propagator.

Definition at line 116 of file SmartPropagator.cc.

References TrajectoryStateOnSurface::freeState(), getGenPropagator(), getTkPropagator(), insideTkVol(), and Propagator::propagateWithPath().

                                                                                                                {
   if (insideTkVol(*fts.freeState()) && insideTkVol(cylinder)) {
     return getTkPropagator()->propagateWithPath(fts, cylinder);
   } else {
     return getGenPropagator()->propagateWithPath(fts, cylinder);
   }
 }

void SmartPropagator::setPropagationDirection ( PropagationDirection dir )

inlineoverridevirtual

setting the direction fo both components

Reimplemented from Propagator.

Definition at line 63 of file SmartPropagator.h.

References Propagator::setPropagationDirection(), theGenProp, and theTkProp.

                                                                   {
     Propagator::setPropagationDirection(dir);
     theTkProp->setPropagationDirection(dir);
     theGenProp->setPropagationDirection(dir);
   }