CMSSW Reference Manual

typedef std::pair<TrajectoryStateOnSurface,double> RKPropagatorInS::TsosWP [private]

Definition at line 60 of file RKPropagatorInS.h.

RKPropagatorInS::RKPropagatorInS	(	const MagVolume &	vol,
		PropagationDirection	dir = `alongMomentum`,
		double	tolerance = `5.e-5`
	)			`[inline, explicit]`

Definition at line 26 of file RKPropagatorInS.h.

Referenced by clone().

00027                                                       : 
00028     Propagator(dir), theVolume( &vol), theTolerance( tolerance) {}

RKPropagatorInS::~RKPropagatorInS ( ) [inline]

Definition at line 30 of file RKPropagatorInS.h.

00030 {}

Propagator * RKPropagatorInS::clone ( void ) const [virtual]

Implements Propagator.

Definition at line 330 of file RKPropagatorInS.cc.

References RKPropagatorInS().

00331 {
00332     return new RKPropagatorInS(*this);
00333 }

RKLocalFieldProvider RKPropagatorInS::fieldProvider ( const Cylinder & cyl ) const [private]

Definition at line 348 of file RKPropagatorInS.cc.

References theVolume.

00349 {
00350   return RKLocalFieldProvider( *theVolume, cyl);
00351 }

RKLocalFieldProvider RKPropagatorInS::fieldProvider ( ) const [private]

Definition at line 343 of file RKPropagatorInS.cc.

References theVolume.

Referenced by propagateParametersOnCylinder(), and propagateParametersOnPlane().

00344 {
00345   return RKLocalFieldProvider( *theVolume);
00346 }

GlobalVector RKPropagatorInS::globalMomentum ( const Basic3DVector< double > & mom ) const [private]

Definition at line 377 of file RKPropagatorInS.cc.

References theVolume, and GloballyPositioned< T >::toGlobal().

Referenced by gtpFromVolumeLocal().

00379 {
00380   if (theVolume != 0) return theVolume->toGlobal( LocalVector(mom));
00381   else return GlobalVector(mom);
00382 }

GlobalPoint RKPropagatorInS::globalPosition ( const Basic3DVector< double > & pos ) const [private]

Definition at line 371 of file RKPropagatorInS.cc.

References theVolume, and GloballyPositioned< T >::toGlobal().

Referenced by gtpFromVolumeLocal(), and propagateParametersOnPlane().

00372 {
00373   if (theVolume != 0) return theVolume->toGlobal( LocalPoint(pos));
00374   else return GlobalPoint(pos);
00375 }

GlobalTrajectoryParameters RKPropagatorInS::gtpFromLocal	(	const Basic3DVector< double > &	lpos,
		const Basic3DVector< double > &	lmom,
		TrackCharge	ch,
		const Surface &	surf
	)			const `[private]`

Definition at line 335 of file RKPropagatorInS.cc.

References theVolume, and Surface::toGlobal().

Referenced by propagateParametersOnCylinder().

00338 {
00339     return GlobalTrajectoryParameters( surf.toGlobal( LocalPoint( lpos)),
00340                                        surf.toGlobal( LocalVector( lmom)), ch, theVolume);
00341 }

GlobalTrajectoryParameters RKPropagatorInS::gtpFromVolumeLocal	(	const CartesianStateAdaptor &	state,
		TrackCharge	charge
	)			const `[private]`

Definition at line 385 of file RKPropagatorInS.cc.

References globalMomentum(), globalPosition(), CartesianStateAdaptor::momentum(), CartesianStateAdaptor::position(), and theVolume.

Referenced by propagateParametersOnPlane().

00387 {
00388   return GlobalTrajectoryParameters( globalPosition(state.position()), 
00389                                      globalMomentum(state.momentum()), 
00390                                      charge, theVolume);
00391 }

PropagationDirection RKPropagatorInS::invertDirection ( PropagationDirection dir ) const [private]

Definition at line 353 of file RKPropagatorInS.cc.

References alongMomentum, anyDirection, and oppositeToMomentum.

Referenced by propagateParametersOnCylinder(), and propagateParametersOnPlane().

00354 {
00355   if (dir == anyDirection) return dir;
00356   return ( dir == alongMomentum ? oppositeToMomentum : alongMomentum);
00357 }

virtual const MagneticField* RKPropagatorInS::magneticField ( ) const [inline, virtual]

Implements Propagator.

Definition at line 56 of file RKPropagatorInS.h.

References theVolume.

Referenced by NavPropagator::propagateInVolume().

00056 {return theVolume;}

TrajectoryStateOnSurface RKPropagatorInS::propagate	(	const TrajectoryStateOnSurface &	ts,
		const Plane &	plane
	)			const `[inline, virtual]`

Reimplemented from Propagator.

Definition at line 47 of file RKPropagatorInS.h.

References first, TrajectoryStateOnSurface::freeState(), and propagateWithPath().

00048                                                                {
00049     return propagateWithPath( *ts.freeState(),plane).first;
00050   }

TrajectoryStateOnSurface RKPropagatorInS::propagate	(	const FreeTrajectoryState &	fts,
		const Cylinder &	cyl
	)			const `[virtual]`

Implements Propagator.

Definition at line 325 of file RKPropagatorInS.cc.

References propagateWithPath().

00326 {
00327   return propagateWithPath( fts, cyl).first;
00328 }

TrajectoryStateOnSurface RKPropagatorInS::propagate	(	const FreeTrajectoryState &	fts,
		const Plane &	plane
	)			const `[virtual]`

Implements Propagator.

Definition at line 319 of file RKPropagatorInS.cc.

References propagateWithPath().

00320 {
00321   return propagateWithPath( fts, plane).first;
00322 }

GlobalParametersWithPath RKPropagatorInS::propagateParametersOnCylinder	(	const FreeTrajectoryState &	ts,
		const Cylinder &	cyl
	)			const `[private]`

Definition at line 191 of file RKPropagatorInS.cc.

References funct::abs(), FreeTrajectoryState::charge(), dist(), e, fieldProvider(), gtpFromLocal(), invertDirection(), LogDebug, FreeTrajectoryState::momentum(), CartesianStateAdaptor::momentum(), path(), StraightLineBarrelCylinderCrossing::pathLength(), StraightLineCylinderCrossing::pathLength(), muonGeometry::perp(), StraightLineBarrelCylinderCrossing::position(), CartesianStateAdaptor::position(), FreeTrajectoryState::position(), Propagator::propagationDirection(), Cylinder::radius(), res, rho, CartesianStateAdaptor::rkstate(), s, theTolerance, theVolume, Surface::toGlobal(), GloballyPositioned< T >::toLocal(), FreeTrajectoryState::transverseCurvature(), PV3DBase< T, PVType, FrameType >::x(), x, and PV3DBase< T, PVType, FrameType >::y().

Referenced by propagateWithPath().

00193 {
00194     typedef RKAdaptiveSolver<double,RKOneCashKarpStep, 6>   Solver;
00195     typedef Solver::Vector                                  RKVector;
00196 
00197     
00198     GlobalPoint sp = cyl.toGlobal(LocalPoint(0.,0.));
00199     if (sp.x()!=0. || sp.y()!=0.) {
00200       throw PropagationException("Cannot propagate to an arbitrary cylinder");
00201     }
00202 
00203     GlobalPoint gpos( ts.position());
00204     GlobalVector gmom( ts.momentum());
00205     LocalPoint pos(cyl.toLocal(gpos));
00206     LocalVector mom(cyl.toLocal(gmom));
00207     double startR = cyl.radius() - pos.perp();
00208 
00209     // LogDebug("RKPropagatorInS")  << "RKPropagatorInS: starting from FTS " << ts ;
00210 
00211 
00212     // (transverse) curvature
00213     double rho = ts.transverseCurvature();
00214     //
00215     // Straight line approximation? |rho|<1.e-10 equivalent to ~ 1um 
00216     // difference in transversal position at 10m.
00217     //
00218     if( fabs(rho)<1.e-10 ) {
00219       //
00220       // Instantiate auxiliary object for finding intersection.
00221       // Frame-independant point and vector are created explicitely to 
00222       // avoid confusing gcc (refuses to compile with temporary objects
00223       // in the constructor).
00224       //
00225       
00226       StraightLineBarrelCylinderCrossing cylCrossing(gpos,gmom,propagationDirection());
00227 
00228       //
00229       // get solution
00230       //
00231       std::pair<bool,double> propResult = cylCrossing.pathLength(cyl);
00232       if ( propResult.first && theVolume !=0) {
00233         double s = propResult.second;
00234         // point (reconverted to GlobalPoint)
00235         GlobalPoint x (cylCrossing.position(s));
00236         GlobalTrajectoryParameters res( x, gmom, ts.charge(), theVolume);
00237         LogDebug("RKPropagatorInS")  << "Straight line propagation to cylinder succeeded !!";
00238         return GlobalParametersWithPath( res, s);
00239       } else {
00240         //do someting 
00241         edm::LogError("RKPropagatorInS")  <<"Straight line propagation to cylinder failed !!";
00242         return GlobalParametersWithPath( );
00243       }
00244     }
00245     
00246 
00247     RKLocalFieldProvider field( fieldProvider(cyl));
00248     // StraightLineCylinderCrossing pathLength( pos, mom, propagationDirection());
00249     CartesianLorentzForce deriv(field, ts.charge());
00250 
00251     RKCartesianDistance dist;
00252     double eps = theTolerance;
00253     Solver solver;
00254     double stot = 0;
00255     PropagationDirection currentDirection = propagationDirection();
00256 
00257     RKVector start( CartesianStateAdaptor::rkstate( pos.basicVector(), mom.basicVector()));
00258     while (true) {
00259       CartesianStateAdaptor startState(start);
00260       StraightLineCylinderCrossing pathLength( LocalPoint(startState.position()), 
00261                                                LocalVector(startState.momentum()), 
00262                                                currentDirection, eps);
00263 
00264       std::pair<bool,double> path = pathLength.pathLength( cyl);
00265       if (!path.first) { 
00266         LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path length calculation to cylinder failed!" 
00267                                      << "Radius " << cyl.radius() << " pos.perp() " << LocalPoint(startState.position()).perp() ;
00268         return GlobalParametersWithPath();
00269       }
00270       else {
00271         LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path lenght to cylinder is " << path.second 
00272                   << " from point (R,z) " << startState.position().perp() 
00273                   << ", " << startState.position().z()
00274                   << " to R " << cyl.radius() 
00275                   ;
00276       }
00277 
00278       double sstep = path.second;
00279       if ( std::abs(sstep) < eps) {
00280         LogDebug("RKPropagatorInS")  << "accuracy not reached, but pathLength calculation says we are there! "
00281              << path.second ;
00282 
00283         GlobalTrajectoryParameters res( gtpFromLocal( startState.position(), 
00284                                                       startState.momentum(), 
00285                                                       ts.charge(), cyl));
00286         return GlobalParametersWithPath( res, stot);
00287       }
00288 
00289       LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Solving for " << sstep 
00290            << " current distance to cylinder is " << startR ;
00291 
00292       RKVector rkresult = solver( 0, start, sstep, deriv, dist, eps);
00293       stot += sstep;
00294       CartesianStateAdaptor cur( rkresult);
00295       double remainingR = cyl.radius() - cur.position().perp();
00296       if ( fabs(remainingR) < eps) {
00297         LogDebug("RKPropagatorInS")  << "Accuracy reached! " << remainingR ;
00298         GlobalTrajectoryParameters res( gtpFromLocal( cur.position(), 
00299                                                       cur.momentum(), 
00300                                                       ts.charge(), cyl));
00301         return GlobalParametersWithPath( res, stot);
00302       }
00303 
00304       start = rkresult;
00305       if (remainingR * startR > 0) {
00306         LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in same direction again " 
00307              << remainingR ;
00308       }
00309       else {
00310         LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in opposite direction " 
00311              << remainingR ;
00312         currentDirection = invertDirection( currentDirection);
00313       }
00314       startR = remainingR;
00315     }
00316 }

GlobalParametersWithPath RKPropagatorInS::propagateParametersOnPlane	(	const FreeTrajectoryState &	ts,
		const Plane &	plane
	)			const `[private]`

Definition at line 54 of file RKPropagatorInS.cc.

References funct::abs(), PV3DBase< T, PVType, FrameType >::basicVector(), FreeTrajectoryState::charge(), dir, dist(), e, fieldProvider(), RKLocalFieldProvider::frame(), globalPosition(), gtpFromVolumeLocal(), RKLocalFieldProvider::inTesla(), invertDirection(), Plane::localZ(), LogDebug, FreeTrajectoryState::momentum(), CartesianStateAdaptor::momentum(), path(), StraightLinePlaneCrossing::pathLength(), CartesianStateAdaptor::position(), StraightLinePlaneCrossing::position(), GloballyPositioned< T >::position(), FreeTrajectoryState::position(), Propagator::propagationDirection(), res, rho, rkMomentum(), rkPosition(), CartesianStateAdaptor::rkstate(), s, theTolerance, theVolume, GloballyPositioned< T >::toGlobal(), GloballyPositioned< T >::toLocal(), FrameChanger::transformPlane(), FreeTrajectoryState::transverseCurvature(), and x.

Referenced by propagateWithPath().

00056 {
00057 
00058   GlobalPoint gpos( ts.position());
00059   GlobalVector gmom( ts.momentum());
00060   double startZ = plane.localZ(gpos);
00061   // (transverse) curvature
00062   double rho = ts.transverseCurvature();
00063   //
00064   // Straight line approximation? |rho|<1.e-10 equivalent to ~ 1um 
00065   // difference in transversal position at 10m.
00066   //
00067   if( fabs(rho)<1.e-10 ) {
00068     //
00069     // Instantiate auxiliary object for finding intersection.
00070     // Frame-independant point and vector are created explicitely to 
00071     // avoid confusing gcc (refuses to compile with temporary objects
00072     // in the constructor).
00073     //
00074     LogDebug("RKPropagatorInS")<<" startZ = "<<startZ;
00075 
00076     if (fabs(startZ) < 1e-5){  
00077       LogDebug("RKPropagatorInS")<< "Propagation is not performed: state is already on final surface.";
00078       GlobalTrajectoryParameters res( gpos, gmom, ts.charge(), 
00079                                       theVolume);
00080       return GlobalParametersWithPath( res, 0.0);
00081     }
00082 
00083     StraightLinePlaneCrossing::PositionType pos(gpos);
00084     StraightLinePlaneCrossing::DirectionType dir(gmom);
00085     StraightLinePlaneCrossing planeCrossing(pos,dir, propagationDirection());
00086     //
00087     // get solution
00088     //
00089     std::pair<bool,double> propResult = planeCrossing.pathLength(plane);
00090     if ( propResult.first && theVolume !=0) {
00091       double s = propResult.second;
00092       // point (reconverted to GlobalPoint)
00093       GlobalPoint x (planeCrossing.position(s));
00094       GlobalTrajectoryParameters res( x, gmom, ts.charge(), theVolume);
00095       return GlobalParametersWithPath( res, s);
00096     } else {
00097       //do someting
00098       LogDebug("RKPropagatorInS")<< "Straight line propgation to plane failed !!"; 
00099       return GlobalParametersWithPath( );
00100     }
00101   }
00102 
00103   if (theVolume != 0) {
00104     LogDebug("RKPropagatorInS")  << "RKPropagatorInS: starting prop to plane in volume with pos " << theVolume->position()
00105               << " Z axis " << theVolume->toGlobal( LocalVector(0,0,1)) ;
00106 
00107     LogDebug("RKPropagatorInS")  << "The starting position is " << ts.position() << " (global) "
00108               << theVolume->toLocal(ts.position()) << " (local) " ;
00109     FrameChanger changer;
00110     FrameChanger::PlanePtr localPlane = changer.transformPlane( plane, *theVolume);
00111     LogDebug("RKPropagatorInS")  << "The plane position is " << plane.position() << " (global) "
00112               << localPlane->position() << " (local) " ;
00113 
00114     LogDebug("RKPropagatorInS")  << "The initial distance to plane is " << plane.localZ( ts.position()) ;
00115 
00116     StraightLinePlaneCrossing cross( ts.position().basicVector(), ts.momentum().basicVector());
00117     std::pair<bool,double> res3 = cross.pathLength(plane);
00118     LogDebug("RKPropagatorInS")  << "straight line distance " << res3.first << " " << res3.second ;
00119   }
00120 
00121   typedef RKAdaptiveSolver<double,RKOneCashKarpStep, 6>   Solver;
00122   typedef Solver::Vector                                  RKVector;
00123   
00124   RKLocalFieldProvider field( fieldProvider());
00125   PathToPlane2Order pathLength( field, &field.frame());
00126   CartesianLorentzForce deriv(field, ts.charge());
00127 
00128   RKCartesianDistance dist;
00129   double eps = theTolerance;
00130   Solver solver;
00131   double stot = 0;
00132   PropagationDirection currentDirection = propagationDirection();
00133 
00134   // in magVolume frame
00135   RKVector start( CartesianStateAdaptor::rkstate( rkPosition(gpos), rkMomentum(gmom)));
00136   while (true) {
00137     CartesianStateAdaptor startState(start);
00138 
00139     std::pair<bool,double> path = pathLength( plane, startState.position(), 
00140                                               startState.momentum(), 
00141                                               (double) ts.charge(), currentDirection);
00142     if (!path.first) { 
00143       LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path length calculation to plane failed!" 
00144                                    << "...distance to plane " << plane.localZ( globalPosition(startState.position()))
00145                                    << "...Local starting position in volume " << startState.position() 
00146                                    << "...Magnetic field " << field.inTesla( startState.position()) ;
00147 
00148 
00149       return GlobalParametersWithPath();
00150     }
00151     else {
00152       LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Path lenght to plane is " << path.second ;
00153     }
00154 
00155     double sstep = path.second;
00156     if ( std::abs(sstep) < eps) {
00157       LogDebug("RKPropagatorInS")  << "On-surface accuracy not reached, but pathLength calculation says we are there! "
00158                  << "path " << path.second << " distance to plane is " << startZ ;
00159       GlobalTrajectoryParameters res( gtpFromVolumeLocal( startState, ts.charge()));
00160       return GlobalParametersWithPath( res, stot);
00161     }
00162 
00163     LogDebug("RKPropagatorInS")  << "RKPropagatorInS: Solving for " << sstep 
00164               << " current distance to plane is " << startZ ;
00165 
00166     RKVector rkresult = solver( 0, start, sstep, deriv, dist, eps);
00167     stot += sstep;
00168     CartesianStateAdaptor cur( rkresult);
00169     double remainingZ = plane.localZ( globalPosition(cur.position()));
00170     if ( fabs(remainingZ) < eps) {
00171       LogDebug("RKPropagatorInS")  << "On-surface accuracy reached! " << remainingZ ;
00172       GlobalTrajectoryParameters res( gtpFromVolumeLocal( cur, ts.charge()));
00173       return GlobalParametersWithPath( res, stot);
00174     }
00175 
00176     start = rkresult;
00177     if (remainingZ * startZ > 0) {
00178       LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in same direction again " 
00179                 << remainingZ ;
00180     }
00181     else {
00182       LogDebug("RKPropagatorInS")  << "Accuracy not reached yet, trying in opposite direction " 
00183                 << remainingZ ;
00184       currentDirection = invertDirection( currentDirection);
00185     }
00186     startZ = remainingZ;
00187   }
00188 }

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInS::propagateWithPath	(	const FreeTrajectoryState &	fts,
		const Cylinder &	cyl
	)			const `[virtual]`

Implements Propagator.

Definition at line 40 of file RKPropagatorInS.cc.

References SurfaceSideDefinition::afterSurface, alongMomentum, SurfaceSideDefinition::beforeSurface, GlobalParametersWithPath::parameters(), propagateParametersOnCylinder(), Propagator::propagationDirection(), and GlobalParametersWithPath::s().

00041 {
00042   GlobalParametersWithPath gp =  propagateParametersOnCylinder( fts, cyl);
00043   if (!gp) return TsosWP(TrajectoryStateOnSurface(),0.);
00044   else {
00045     SurfaceSideDefinition::SurfaceSide side = PropagationDirectionFromPath()(gp.s(),propagationDirection())==alongMomentum 
00046       ? SurfaceSideDefinition::beforeSurface : SurfaceSideDefinition::afterSurface;
00047     AnalyticalErrorPropagation errorprop;
00048     return errorprop( fts, cyl, side, gp.parameters(),gp.s());
00049   }
00050   
00051 }

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInS::propagateWithPath	(	const FreeTrajectoryState &	fts,
		const Plane &	plane
	)			const `[virtual]`

Implements Propagator.

Definition at line 26 of file RKPropagatorInS.cc.

References SurfaceSideDefinition::afterSurface, alongMomentum, SurfaceSideDefinition::beforeSurface, GlobalParametersWithPath::parameters(), propagateParametersOnPlane(), Propagator::propagationDirection(), and GlobalParametersWithPath::s().

Referenced by propagate(), and NavPropagator::propagateInVolume().

00028 {
00029   GlobalParametersWithPath gp =  propagateParametersOnPlane( fts, plane);
00030   if (!gp) return TsosWP(TrajectoryStateOnSurface(),0.);
00031   else {
00032     SurfaceSideDefinition::SurfaceSide side = PropagationDirectionFromPath()(gp.s(),propagationDirection())==alongMomentum 
00033       ? SurfaceSideDefinition::beforeSurface : SurfaceSideDefinition::afterSurface;
00034     AnalyticalErrorPropagation errorprop;
00035     return errorprop( fts, plane, side, gp.parameters(),gp.s());
00036   }
00037 }

Basic3DVector< double > RKPropagatorInS::rkMomentum ( const GlobalVector & mom ) const [private]

Definition at line 365 of file RKPropagatorInS.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), theVolume, and GloballyPositioned< T >::toLocal().

Referenced by propagateParametersOnPlane().

00366 {
00367   if (theVolume != 0) return theVolume->toLocal( mom).basicVector();
00368   else return mom.basicVector();
00369 }

Basic3DVector< double > RKPropagatorInS::rkPosition ( const GlobalPoint & pos ) const [private]

Definition at line 359 of file RKPropagatorInS.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), theVolume, and GloballyPositioned< T >::toLocal().

Referenced by propagateParametersOnPlane().

00360 {
00361   if (theVolume != 0) return theVolume->toLocal( pos).basicVector();
00362   else return pos.basicVector();
00363 }

double RKPropagatorInS::theTolerance [private]

Definition at line 63 of file RKPropagatorInS.h.

Referenced by propagateParametersOnCylinder(), and propagateParametersOnPlane().

const MagVolume* RKPropagatorInS::theVolume [private]

Definition at line 62 of file RKPropagatorInS.h.

Referenced by fieldProvider(), globalMomentum(), globalPosition(), gtpFromLocal(), gtpFromVolumeLocal(), magneticField(), propagateParametersOnCylinder(), propagateParametersOnPlane(), rkMomentum(), and rkPosition().


Public Member Functions
virtual Propagator *	clone () const
virtual const MagneticField *	magneticField () const
TrajectoryStateOnSurface	propagate (const TrajectoryStateOnSurface &ts, const Plane &plane) const
virtual TrajectoryStateOnSurface	propagate (const FreeTrajectoryState &, const Cylinder &) const
virtual TrajectoryStateOnSurface	propagate (const FreeTrajectoryState &, const Plane &) const
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Cylinder &) const
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Plane &) const
	RKPropagatorInS (const MagVolume &vol, PropagationDirection dir=alongMomentum, double tolerance=5.e-5)
	~RKPropagatorInS ()
Private Types
typedef std::pair < TrajectoryStateOnSurface, double >	TsosWP
Private Member Functions
RKLocalFieldProvider	fieldProvider (const Cylinder &cyl) const
RKLocalFieldProvider	fieldProvider () const
GlobalVector	globalMomentum (const Basic3DVector< double > &mom) const
GlobalPoint	globalPosition (const Basic3DVector< double > &pos) const
GlobalTrajectoryParameters	gtpFromLocal (const Basic3DVector< double > &lpos, const Basic3DVector< double > &lmom, TrackCharge ch, const Surface &surf) const
GlobalTrajectoryParameters	gtpFromVolumeLocal (const CartesianStateAdaptor &state, TrackCharge charge) const
PropagationDirection	invertDirection (PropagationDirection dir) const
GlobalParametersWithPath	propagateParametersOnCylinder (const FreeTrajectoryState &ts, const Cylinder &cyl) const
GlobalParametersWithPath	propagateParametersOnPlane (const FreeTrajectoryState &ts, const Plane &plane) const
Basic3DVector< double >	rkMomentum (const GlobalVector &mom) const
Basic3DVector< double >	rkPosition (const GlobalPoint &pos) const
Private Attributes
double	theTolerance
const MagVolume *	theVolume

RKPropagatorInS Class Reference

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation