#include <SmartPropagator.h>

Inheritance diagram for SmartPropagator:

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

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

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 BoundCylinder &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
	return the magneticField More...

TrajectoryStateOnSurface	propagate (const FreeTrajectoryState &fts, const Surface &surface) const

TrajectoryStateOnSurface	propagate (const TrajectoryStateOnSurface &tsos, const Surface &surface) const

TrajectoryStateOnSurface	propagate (const FreeTrajectoryState &fts, const Plane &plane) const

TrajectoryStateOnSurface	propagate (const TrajectoryStateOnSurface &tsos, const Plane &plane) const

TrajectoryStateOnSurface	propagate (const FreeTrajectoryState &fts, const Cylinder &cylinder) const

TrajectoryStateOnSurface	propagate (const TrajectoryStateOnSurface &tsos, const Cylinder &cylinder) const

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

std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Surface &surface) const

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

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

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

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

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

	SmartPropagator (Propagator aTkProp, 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
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &, const reco::BeamSpot &) const

virtual std::pair < FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &, const GlobalPoint &, const GlobalPoint &) const

virtual PropagationDirection	propagationDirection () const

	Propagator (PropagationDirection dir=alongMomentum)

virtual bool	setMaxDirectionChange (float phiMax)

virtual	~Propagator ()

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

static ReferenceCountingPointer < BoundCylinder > &	theTkVolume ()

Private Attributes
const MagneticField *	theField

Propagator *	theGenProp

Propagator *	theTkProp

Detailed Description

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

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

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 37 of file SmartPropagator.h.

Constructor & Destructor Documentation

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

Defines which propagator is used inside Tk and which outside.

Definition at line 44 of file SmartPropagator.cc.

References initTkVolume(), and theTkVolume().

Referenced by clone().

                                                                           :
   Propagator(dir), theTkProp(aTkProp->clone()), theGenProp(aGenProp->clone()), theField(field) { 
 
   if (theTkVolume()==0) 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 53 of file SmartPropagator.cc.

References initTkVolume(), and theTkVolume().

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

SmartPropagator::SmartPropagator ( const SmartPropagator & aProp )

Copy constructor.

Definition at line 62 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
     //if (theTkVolume==0) initTkVolume(epsilon);
 
   }

SmartPropagator::~SmartPropagator ( )

virtual

virtual destructor

Definition at line 76 of file SmartPropagator.cc.

References theGenProp, and theTkProp.

                                   {
 
   delete theTkProp;
   delete theGenProp;
 
 }

Member Function Documentation

virtual SmartPropagator* SmartPropagator::clone ( void ) const

inlinevirtual

Virtual constructor (using copy c'tor)

Implements Propagator.

Definition at line 57 of file SmartPropagator.h.

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

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

Propagator * SmartPropagator::getGenPropagator ( ) const

return the propagator used outside tracker

Definition at line 205 of file SmartPropagator.cc.

References theGenProp.

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

                                                     {
 
   return theGenProp;
 
 }

Propagator * SmartPropagator::getTkPropagator ( ) const

return the propagator used inside tracker

Definition at line 198 of file SmartPropagator.cc.

References theTkProp.

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

                                                    {
 
   return theTkProp;
 
 }

void SmartPropagator::initTkVolume ( float epsilon )

staticprivate

build the tracker volume

Definition at line 85 of file SmartPropagator.cc.

References BoundCylinder::build(), epsilon, TrackerBounds::halfLength(), pos, TrackerBounds::radius(), CosmicsPD_Skims::radius, 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() = BoundCylinder::build(pos, rot, radius, 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 154 of file SmartPropagator.cc.

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

Referenced by propagate(), and 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 164 of file SmartPropagator.cc.

References TrackerBounds::halfLength(), 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 BoundCylinder & cylin ) const

true if a cylinder is inside tracker volume

Definition at line 176 of file SmartPropagator.cc.

References BoundSurface::bounds(), TrackerBounds::halfLength(), Bounds::length(), PV3DBase< T, PVType, FrameType >::perp(), TrackerBounds::radius(), Cylinder::radius(), and PV3DBase< T, PVType, FrameType >::z().

                                                                     {
 
   GlobalPoint gp(cylin.radius(),0.,(cylin.bounds().length())/2.);
 //  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 Plane & plane ) const

true if a plane is inside tracker volume

Definition at line 187 of file SmartPropagator.cc.

References 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

inlinevirtual

return the magneticField

Implements Propagator.

Definition at line 141 of file SmartPropagator.h.

References theField.

Referenced by clone().

141 {return theField;}

SmartPropagator::theField

const MagneticField * theField

Definition: SmartPropagator.h:149

TrajectoryStateOnSurface SmartPropagator::propagate	(	const FreeTrajectoryState &	state,
		const Surface &	sur
	)		const

virtual

Propagate from a free state (e.g. position and momentum in in global cartesian coordinates) to a surface.Only use the generic method if the surface type (plane or cylinder) is not known at the calling point.

Reimplemented from Propagator.

Definition at line 104 of file SmartPropagator.cc.

References Propagator::propagate().

                                                                                   {
   return Propagator::propagate( fts, surface);
 }

TrajectoryStateOnSurface SmartPropagator::propagate	(	const TrajectoryStateOnSurface &	state,
		const Surface &	sur
	)		const

inlinevirtual

The following three methods are equivalent to the corresponding methods above, but if the starting state is a TrajectoryStateOnSurface, it's better to use it as such rather than use just the FreeTrajectoryState part. It may help some concrete propagators.Only use the generic method if the surface type (plane or cylinder) is not known at the calling point.

Reimplemented from Propagator.

Definition at line 74 of file SmartPropagator.h.

References Propagator::propagate().

                                                                      {
       return Propagator::propagate(tsos,surface);
     }

TrajectoryStateOnSurface SmartPropagator::propagate	(	const FreeTrajectoryState &	fts,
		const Plane &	plane
	)		const

virtual

Implements Propagator.

Definition at line 110 of file SmartPropagator.cc.

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

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

TrajectoryStateOnSurface SmartPropagator::propagate	(	const TrajectoryStateOnSurface &	tsos,
		const Plane &	plane
	)		const

inlinevirtual

Reimplemented from Propagator.

Definition at line 82 of file SmartPropagator.h.

References Propagator::propagate().

                                                                  {
       return Propagator::propagate(tsos, plane);
     }

TrajectoryStateOnSurface SmartPropagator::propagate	(	const FreeTrajectoryState &	fts,
		const Cylinder &	cylinder
	)		const

virtual

Implements Propagator.

Definition at line 122 of file SmartPropagator.cc.

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

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

TrajectoryStateOnSurface SmartPropagator::propagate	(	const TrajectoryStateOnSurface &	tsos,
		const Cylinder &	cylinder
	)		const

inlinevirtual

Reimplemented from Propagator.

Definition at line 90 of file SmartPropagator.h.

References Propagator::propagate().

                                                                        {
       return Propagator::propagate(tsos, cylinder);
     }

std::pair<TrajectoryStateOnSurface,double> SmartPropagator::propagateWithPath	(	const FreeTrajectoryState &	state,
		const Surface &	sur
	)		const

inlinevirtual

The methods propagateWithPath() are identical to the corresponding methods propagate() in what concerns the resulting TrajectoryStateOnSurface, but they provide in addition the exact path length along the trajectory.Only use the generic method if the surface type (plane or cylinder) is not known at the calling point.

Reimplemented from Propagator.

Definition at line 96 of file SmartPropagator.h.

References Propagator::propagateWithPath().

                                                       {
         return Propagator::propagateWithPath(fts,surface);
       }

std::pair<TrajectoryStateOnSurface,double> SmartPropagator::propagateWithPath	(	const TrajectoryStateOnSurface &	state,
		const Surface &	sur
	)		const

inlinevirtual

The following three methods are equivalent to the corresponding methods above, but if the starting state is a TrajectoryStateOnSurface, it's better to use it as such rather than use just the FreeTrajectoryState part. It may help some concrete propagators.Only use the generic method if the surface type (plane or cylinder) is not known at the calling point.

Reimplemented from Propagator.

Definition at line 102 of file SmartPropagator.h.

References Propagator::propagateWithPath().

                                                       {
         return Propagator::propagateWithPath(tsos,surface);
       }

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

virtual

Implements Propagator.

Definition at line 133 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 TrajectoryStateOnSurface &	tsos,
		const Plane &	plane
	)		const

inlinevirtual

Reimplemented from Propagator.

Definition at line 112 of file SmartPropagator.h.

References Propagator::propagateWithPath().

                                                   {
         return Propagator::propagateWithPath(tsos, plane);
       }

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

virtual

Implements Propagator.

Definition at line 144 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 Cylinder &	cylinder
	)		const

inlinevirtual

Reimplemented from Propagator.

Definition at line 122 of file SmartPropagator.h.

References Propagator::propagateWithPath().

                                                         {
         return Propagator::propagateWithPath(tsos, cylinder);
       }

void SmartPropagator::setPropagationDirection ( PropagationDirection dir ) const

inlinevirtual

setting the direction fo both components

Reimplemented from Propagator.

Definition at line 62 of file SmartPropagator.h.

References getGenPropagator(), getTkPropagator(), and Propagator::setPropagationDirection().

     {
       Propagator::setPropagationDirection (dir);
       getTkPropagator()->setPropagationDirection(dir);
       getGenPropagator()->setPropagationDirection(dir);
     }

ReferenceCountingPointer< BoundCylinder > & SmartPropagator::theTkVolume ( )

staticprivate

Definition at line 36 of file SmartPropagator.cc.

Referenced by initTkVolume(), and SmartPropagator().

                                                                        {
   static ReferenceCountingPointer<BoundCylinder> local=0;
   return local;
 }