#include <SmartPropagator.h>

Inheritance diagram for SmartPropagator:

Public Member Functions
virtual 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...

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

virtual	~SmartPropagator ()

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

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

virtual 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 35 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 38 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 47 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 56 of file SmartPropagator.cc.

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

                                                              :
   Propagator(aProp.propagationDirection()), theTkProp(0), theGenProp(0) { 
     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 ( )

virtual

virtual destructor

Definition at line 70 of file SmartPropagator.cc.

References theGenProp, and theTkProp.

                                   {
 
   delete theTkProp;
   delete theGenProp;
 
 }

Member Function Documentation

virtual SmartPropagator* SmartPropagator::clone ( void ) const

inlineoverridevirtual

Virtual constructor (using copy c'tor)

Implements Propagator.

Definition at line 55 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 195 of file SmartPropagator.cc.

References theGenProp.

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

                                                           {
 
   return theGenProp;
 
 }

const Propagator * SmartPropagator::getTkPropagator ( ) const

return the propagator used inside tracker

Definition at line 188 of file SmartPropagator.cc.

References theTkProp.

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

                                                          {
 
   return theTkProp;
 
 }

void SmartPropagator::initTkVolume ( float epsilon )

private

build the tracker volume

Definition at line 79 of file SmartPropagator.cc.

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

Referenced by magneticField(), and 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 144 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 insideTkVol(), propagateWithPath(), and setPropagationDirection().

                                                                       {
 
   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 154 of file SmartPropagator.cc.

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

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

                                                             {
 
   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.) );
 
 
 }

virtual const MagneticField* SmartPropagator::magneticField ( ) const

inlineoverridevirtual

return the magneticField

Implements Propagator.

Definition at line 105 of file SmartPropagator.h.

References geometryDiff::epsilon, initTkVolume(), and theField.

Referenced by clone().

105 {return theField;}

SmartPropagator::theField

const MagneticField * theField

Definition: SmartPropagator.h:113

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

overrideprivatevirtual

Implements Propagator.

Definition at line 100 of file SmartPropagator.cc.

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

Referenced by setPropagationDirection().

 {
   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 111 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 123 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 134 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 60 of file SmartPropagator.h.

References getGenPropagator(), getTkPropagator(), insideTkVol(), Propagator::propagate(), Propagator::propagateWithPath(), propagateWithPath(), Propagator::setPropagationDirection(), theGenProp, and theTkProp.

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