#include <GlobalTrajectoryParameters.h>

Public Member Functions
TrackCharge	charge () const

GlobalVector	direction () const

	GlobalTrajectoryParameters ()

	GlobalTrajectoryParameters (const GlobalPoint &aX, const GlobalVector &aP, TrackCharge aCharge, const MagneticField *fieldProvider)

	GlobalTrajectoryParameters (const GlobalPoint &aX, const GlobalVector &direction, float transverseCurvature, int, const MagneticField *fieldProvider)

const MagneticField &	magneticField () const

GlobalVector	magneticFieldInInverseGeV (const GlobalPoint &x) const

GlobalVector	momentum () const

GlobalPoint	position () const

float	signedInverseMomentum () const

float	signedInverseTransverseMomentum () const

float	transverseCurvature () const

AlgebraicVector6	vector () const

Private Attributes
float	cachedCurvature_

bool	hasCurvature_

signed char	theCharge

const MagneticField *	theField

GlobalVector	theP

GlobalPoint	theX

Detailed Description

Class providing access to a set of relevant parameters of a trajectory in the global, Cartesian frame. The basic data members used to calculate these parameters are the charge and global position and momentum.

Definition at line 16 of file GlobalTrajectoryParameters.h.

Constructor & Destructor Documentation

GlobalTrajectoryParameters::GlobalTrajectoryParameters ( )

inline

Definition at line 19 of file GlobalTrajectoryParameters.h.

                                :
     theField(0), 
     theX(), theP(), 
     cachedCurvature_(0.0),
     theCharge(0),
     hasCurvature_(false)  {}  // we must initialize cache to non-NAN to avoid FPE

GlobalTrajectoryParameters::GlobalTrajectoryParameters	(	const GlobalPoint &	aX,
		const GlobalVector &	aP,
		TrackCharge	aCharge,
		const MagneticField *	fieldProvider
	)

inline

Constructing class from global position, global momentum and charge.

Definition at line 29 of file GlobalTrajectoryParameters.h.

                                                      :
     theField(fieldProvider),
     theX(aX), theP(aP),     
     cachedCurvature_(aCharge),
     theCharge(aCharge),
     hasCurvature_(false) {
   } // we must initialize cache to non-NAN to avoid FPE

GlobalTrajectoryParameters::GlobalTrajectoryParameters	(	const GlobalPoint &	aX,
		const GlobalVector &	direction,
		float	transverseCurvature,
		int	,
		const MagneticField *	fieldProvider
	)

Constructing class from global position, direction (unit length) and transverse curvature. The fourth int argument is dummy, it serves only to distinguish this constructor from the one above.

Definition at line 5 of file GlobalTrajectoryParameters.cc.

References abs, f, MagneticField::inTesla(), PV3DBase< T, PVType, FrameType >::perp(), theCharge, theField, theP, theX, transverseCurvature(), and PV3DBase< T, PVType, FrameType >::z().

                                                                    :
   theField(fieldProvider),
   theX(aX), cachedCurvature_(transverseCurvature),  hasCurvature_(true)
 {
   float bza = -2.99792458e-3f * theField->inTesla(theX).z();
   float qbpi = bza/(direction.perp()*transverseCurvature);
   theP = direction*std::abs(qbpi);
   theCharge = qbpi > 0.f ? 1 : -1;
 }

Member Function Documentation

TrackCharge GlobalTrajectoryParameters::charge ( void ) const

inline

Charge q of particle, either +1 or -1.

Definition at line 71 of file GlobalTrajectoryParameters.h.

References theCharge.

Referenced by TwoTrackMinimumDistanceLineLine::calculate(), FreeTrajectoryState::charge(), MuonPatternRecoDumper::dumpFTS(), SeedFromConsecutiveHitsCreator::initialKinematic(), SeedForPhotonConversion1Leg::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), JacobianCartesianToCurvilinear::JacobianCartesianToCurvilinear(), JacobianCurvilinearToCartesian::JacobianCurvilinearToCartesian(), ConvBremSeedProducer::produce(), TangentApproachInRPhi::trajectoryParameters(), and TwoTrackMinimumDistanceHelixHelix::updateCoeffs().

                              {
     return theCharge;
   }

GlobalVector GlobalTrajectoryParameters::direction ( ) const

inline

Definition at line 64 of file GlobalTrajectoryParameters.h.

References theP, and Vector3DBase< T, FrameTag >::unit().

                                  {
     return theP.unit();
   }

const MagneticField& GlobalTrajectoryParameters::magneticField ( ) const

inline

Definition at line 113 of file GlobalTrajectoryParameters.h.

References theField.

Referenced by TransientTrackKinematicStateBuilder::buildState(), TangentApproachInRPhi::calculate(), TwoTrackMinimumDistanceLineLine::calculate(), TrajectoryStateCombiner::combine(), fastProp(), PerigeeConversions::ftsToPerigeeParameters(), SeedForPhotonConversion1Leg::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), PerigeeConversions::jacobianCurvilinear2Perigee(), PerigeeConversions::jacobianPerigee2Curvilinear(), BasicTrajectoryState::magneticField(), TrajectoryStateWithArbitraryError::operator()(), TrackKinematicStatePropagator::planeCrossing(), RecHitPropagator::propagate(), NuclearInteractionFinder::rescaleError(), FreeTrajectoryState::rescaleError(), TangentApproachInRPhi::trajectoryParameters(), TrajectoryStateClosestToPoint::TrajectoryStateClosestToPoint(), PerigeeRefittedTrackState::trajectoryStateOnSurface(), GsfMultiStateUpdator::update(), TwoTrackMinimumDistanceHelixHelix::updateCoeffs(), GsfMaterialEffectsUpdator::updateState(), VolumeMaterialEffectsUpdator::updateState(), and MaterialEffectsUpdator::updateStateInPlace().

113 {return *theField;}

GlobalTrajectoryParameters::theField

const MagneticField * theField

Definition: GlobalTrajectoryParameters.h:117

GlobalVector GlobalTrajectoryParameters::magneticFieldInInverseGeV ( const GlobalPoint & x ) const

Definition at line 28 of file GlobalTrajectoryParameters.cc.

References f, MagneticField::inTesla(), and theField.

Referenced by ReferenceTrajectory::addMaterialEffectsBrl(), AnalyticalCurvilinearJacobian::AnalyticalCurvilinearJacobian(), and operator<<().

 {
   return 2.99792458e-3f * theField->inTesla(x);
 }

GlobalVector GlobalTrajectoryParameters::momentum ( ) const

inline

Global momentum.

Definition at line 60 of file GlobalTrajectoryParameters.h.

References theP.

Referenced by ReferenceTrajectory::addMaterialEffectsBrl(), MuonTrackLoader::buildTrackExtra(), TwoTrackMinimumDistanceLineLine::calculate(), VZeroFinder::checkTrackPair(), PixelClusterShapeSeedComparitor::compatible(), AnalyticalCurvilinearJacobian::computeFullJacobian(), AnalyticalCurvilinearJacobian::computeInfinitesimalJacobian(), AnalyticalCurvilinearJacobian::computeStraightLineJacobian(), NuclearTrackCorrector::getNewTrackExtra(), SeedFromConsecutiveHitsCreator::initialKinematic(), SeedForPhotonConversion1Leg::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), JacobianCartesianToCurvilinear::JacobianCartesianToCurvilinear(), JacobianCurvilinearToCartesian::JacobianCurvilinearToCartesian(), JacobianCurvilinearToLocal::JacobianCurvilinearToLocal(), PerigeeConversions::jacobianPerigee2Curvilinear(), SeedFromConsecutiveHitsCreator::makeSeed(), FreeTrajectoryState::momentum(), AnalyticalErrorPropagation::operator()(), ConvBremSeedProducer::produce(), TrackingRecHitProjector< ResultingHit >::project(), KfTrackProducerBase::putInEvt(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), KalmanAlignmentTrackRefitter::refitSingleTracklet(), TangentApproachInRPhi::trajectoryParameters(), SeedForPhotonConversion1Leg::trajectorySeed(), and TwoTrackMinimumDistanceHelixHelix::updateCoeffs().

                                 {
     return theP;
   }

GlobalPoint GlobalTrajectoryParameters::position ( ) const

inline

Global position.

Definition at line 53 of file GlobalTrajectoryParameters.h.

References theX.

Referenced by ReferenceTrajectory::addMaterialEffectsBrl(), AnalyticalCurvilinearJacobian::AnalyticalCurvilinearJacobian(), MuonTrackLoader::buildTrackExtra(), TwoTrackMinimumDistanceHelixHelix::calculate(), TwoTrackMinimumDistanceLineLine::calculate(), PixelClusterShapeSeedComparitor::compatible(), PixelHitMatcher::compatibleSeeds(), AnalyticalCurvilinearJacobian::computeFullJacobian(), ConversionForwardEstimator::estimate(), ConversionBarrelEstimator::estimate(), ForwardMeasurementEstimator::estimate(), BarrelMeasurementEstimator::estimate(), AnalyticalImpactPointExtrapolator::extrapolateSingleState(), AnalyticalTrajectoryExtrapolatorToLine::extrapolateSingleState(), TwoTrackMinimumDistanceHelixHelix::finalPoints(), NuclearTrackCorrector::getNewTrackExtra(), SeedFromConsecutiveHitsCreator::initialKinematic(), SeedForPhotonConversion1Leg::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), JacobianCurvilinearToLocal::JacobianCurvilinearToLocal(), JacobianLocalToCurvilinear::JacobianLocalToCurvilinear(), PerigeeConversions::jacobianPerigee2Curvilinear(), AnalyticalErrorPropagation::operator()(), FreeTrajectoryState::position(), ConvBremSeedProducer::produce(), TauTagTools::propagTrackECALSurfContactPoint(), KfTrackProducerBase::putInEvt(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), TangentApproachInRPhi::trajectoryParameters(), and TwoTrackMinimumDistanceHelixHelix::updateCoeffs().

                                {
     return theX;
   }

float GlobalTrajectoryParameters::signedInverseMomentum ( ) const

inline

Charge divided by (magnitude of) momentum, i.e. q/p.

Definition at line 78 of file GlobalTrajectoryParameters.h.

References PV3DBase< T, PVType, FrameType >::mag(), theCharge, and theP.

Referenced by AnalyticalCurvilinearJacobian::computeFullJacobian(), AnalyticalCurvilinearJacobian::computeInfinitesimalJacobian(), PerigeeConversions::jacobianPerigee2Curvilinear(), and FreeTrajectoryState::signedInverseMomentum().

                                       {
     return theCharge/theP.mag();
   }

float GlobalTrajectoryParameters::signedInverseTransverseMomentum ( ) const

inline

Charge divided by transverse momentum, i.e. q/p_T.

Definition at line 85 of file GlobalTrajectoryParameters.h.

References PV3DBase< T, PVType, FrameType >::perp(), theCharge, and theP.

Referenced by transverseCurvature().

                                                 {
     return theCharge/theP.perp();
   }

float GlobalTrajectoryParameters::transverseCurvature ( ) const

Transverse curvature kappa (which is the inverse radius of curvature in the transverse plane) in cm^{-1}. Sign convention is such that positive kappa means counterclockwise rotation of the track with respect to the global z-axis.

Definition at line 18 of file GlobalTrajectoryParameters.cc.

References cachedCurvature_, f, hasCurvature_, MagneticField::inTesla(), signedInverseTransverseMomentum(), theField, theX, unlikely, and PV3DBase< T, PVType, FrameType >::z().

Referenced by AnalyticalCurvilinearJacobian::AnalyticalCurvilinearJacobian(), GlobalTrajectoryParameters(), SeedForPhotonConversion1Leg::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), PerigeeConversions::jacobianPerigee2Curvilinear(), and FreeTrajectoryState::transverseCurvature().

 {
   if unlikely(!hasCurvature_) {
       float bza = -2.99792458e-3f * theField->inTesla(theX).z();
       cachedCurvature_ = bza*signedInverseTransverseMomentum();
       hasCurvature_ = true;
   }
   return cachedCurvature_;
 }

AlgebraicVector6 GlobalTrajectoryParameters::vector ( ) const

inline

Vector whose first three elements are the global position coordinates and whose last three elements are the global momentum coordinates.

Definition at line 100 of file GlobalTrajectoryParameters.h.

References theP, theX, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by GlobalTrackerMuonAlignment::debugTrajectorySOS(), GlobalTrackerMuonAlignment::debugTrajectorySOSv(), operator<<(), MinPtTrajectoryFilter::test(), and ThresholdPtTrajectoryFilter::test().

                                   {
     //AlgebraicVector6 v;
     //v[0] = theX.x();
     //v[1] = theX.y();
     //v[2] = theX.z();
     //v[3] = theP.x();
     //v[4] = theP.y();
     //v[5] = theP.z();
     return AlgebraicVector6(theX.x(),theX.y(),theX.z(),theP.x(),theP.y(),theP.z());
   }