GlobalTrajectoryParameters Class Reference

#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 &aP, TrackCharge aCharge, const MagneticField *fieldProvider, GlobalVector fieldValue)
	GlobalTrajectoryParameters (const GlobalPoint &aX, const GlobalVector &direction, float transverseCurvature, int, const MagneticField *fieldProvider)
	GlobalTrajectoryParameters (const GlobalPoint &aX, const GlobalVector &direction, float transverseCurvature, int, const MagneticField *fieldProvider, GlobalVector fieldValue)
const MagneticField &	magneticField () const
GlobalVector	magneticFieldInInverseGeV (const GlobalPoint &x) const
GlobalVector	magneticFieldInInverseGeV () const
GlobalVector	magneticFieldInTesla () const
GlobalVector	momentum () const
GlobalPoint	position () const
float	signedInverseMomentum () const
float	signedInverseTransverseMomentum () const
float	transverseCurvature () const
AlgebraicVector6	vector () const

Private Member Functions
void	setCache ()

Private Attributes
GlobalVector	cachedMagneticField
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 15 of file GlobalTrajectoryParameters.h.

Constructor & Destructor Documentation

GlobalTrajectoryParameters::GlobalTrajectoryParameters ( )

inline

Definition at line 18 of file GlobalTrajectoryParameters.h.

Referenced by GlobalTrajectoryParameters().

                               :
    theField(nullptr), 
    theX(), theP(), 
    theCharge(0)
  {}  // 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 27 of file GlobalTrajectoryParameters.h.

References setCache().

                                                                :
    theField(fieldProvider),
    theX(aX), theP(aP),
    theCharge(aCharge)
    {setCache();}

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

inline

Definition at line 37 of file GlobalTrajectoryParameters.h.

References direction(), GlobalTrajectoryParameters(), and transverseCurvature().

                                                     :
    theField(fieldProvider),
    theX(aX), theP(aP),
    cachedMagneticField(fieldValue),
    theCharge(aCharge)
    {}

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 funct::abs(), cachedMagneticField, f, MagneticField::inTesla(), PV3DBase< T, PVType, FrameType >::perp(), theCharge, theField, theP, theX, transverseCurvature(), and PV3DBase< T, PVType, FrameType >::z().

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

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

Definition at line 20 of file GlobalTrajectoryParameters.cc.

References funct::abs(), cachedMagneticField, f, PV3DBase< T, PVType, FrameType >::perp(), theCharge, theP, transverseCurvature(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                :
  theField(fieldProvider),
  theX(aX),    
  cachedMagneticField(fieldValue)
{
  float bza = -2.99792458e-3f * cachedMagneticField.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 89 of file GlobalTrajectoryParameters.h.

References theCharge.

Referenced by TwoTrackMinimumDistanceLineLine::calculate(), TwoTrackMinimumDistance::calculate(), FreeTrajectoryState::charge(), MuonPatternRecoDumper::dumpFTS(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), SeedForPhotonConversion1Leg::initialKinematic(), SeedFromConsecutiveHitsCreator::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), ClosestApproachInRPhi::newTrajectory(), TwoTrackMinimumDistance::pointsHelixHelix(), ConvBremSeedProducer::produce(), Geant4ePropagator::propagateGeneric(), TangentApproachInRPhi::trajectoryParameters(), and TwoTrackMinimumDistanceHelixHelix::updateCoeffs().

                             {
    return theCharge;
  }

GlobalVector GlobalTrajectoryParameters::direction ( ) const

inline

Definition at line 82 of file GlobalTrajectoryParameters.h.

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

Referenced by GlobalTrajectoryParameters().

                                 {
    return theP.unit();
  }

const MagneticField& GlobalTrajectoryParameters::magneticField ( ) const

inline

Definition at line 134 of file GlobalTrajectoryParameters.h.

References setCache(), and theField.

Referenced by TangentApproachInRPhi::calculate(), TwoTrackMinimumDistanceLineLine::calculate(), ClosestApproachInRPhi::calculate(), TwoTrackMinimumDistance::calculate(), TrajectoryStateCombiner::combine(), TSCPBuilderNoMaterial::createFTSatTransverseImpactPointCharged(), TSCPBuilderNoMaterial::createFTSatTransverseImpactPointNeutral(), fastProp(), SeedForPhotonConversion1Leg::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), KinematicState::magneticField(), SeedFromConsecutiveHitsCreator::makeSeed(), ClosestApproachInRPhi::newTrajectory(), TSCBLBuilderNoMaterial::operator()(), TrajectoryStateWithArbitraryError::operator()(), RecHitPropagator::propagate(), Geant4ePropagator::propagateGeneric(), NuclearInteractionFinder::rescaleError(), FreeTrajectoryState::rescaleError(), TangentApproachInRPhi::trajectoryParameters(), TrajectoryStateClosestToPoint::TrajectoryStateClosestToPoint(), PerigeeRefittedTrackState::trajectoryStateOnSurface(), GsfMultiStateUpdator::update(), KFStrip1DUpdator::update(), TwoTrackMinimumDistanceHelixHelix::updateCoeffs(), GsfMaterialEffectsUpdator::updateState(), and VolumeMaterialEffectsUpdator::updateState().

{return *theField;}

GlobalVector GlobalTrajectoryParameters::magneticFieldInInverseGeV

(

const GlobalPoint &

x

)

const

Definition at line 41 of file GlobalTrajectoryParameters.cc.

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

Referenced by ReferenceTrajectory::addMaterialEffectsBrl(), AnalyticalCurvilinearJacobian::AnalyticalCurvilinearJacobian(), PerigeeConversions::ftsToPerigeeParameters(), PerigeeConversions::jacobianCurvilinear2Perigee(), JacobianCurvilinearToLocal::JacobianCurvilinearToLocal(), JacobianLocalToCurvilinear::JacobianLocalToCurvilinear(), PerigeeConversions::jacobianPerigee2Curvilinear(), KinematicState::magneticFieldInInverseGeV(), operator<<(), and TrackKinematicStatePropagator::propagateToTheTransversePCA().

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

GlobalVector GlobalTrajectoryParameters::magneticFieldInInverseGeV ( ) const

inline

Definition at line 126 of file GlobalTrajectoryParameters.h.

References cachedMagneticField, and f.

Referenced by vector().

                                                 {
    return 2.99792458e-3f * cachedMagneticField;
  }

GlobalVector GlobalTrajectoryParameters::magneticFieldInTesla ( ) const

inline

Definition at line 130 of file GlobalTrajectoryParameters.h.

References cachedMagneticField.

Referenced by BasicTrajectoryState::notValid(), and BasicTrajectoryState::update().

                                            {
    return cachedMagneticField;
  }

GlobalVector GlobalTrajectoryParameters::momentum ( ) const

inline

Global momentum.

Definition at line 78 of file GlobalTrajectoryParameters.h.

References theP.

Referenced by ReferenceTrajectory::addMaterialEffectsBrl(), analyticalErrorPropagation(), MuonTrackLoader::buildTrackExtra(), TwoTrackMinimumDistanceLineLine::calculate(), TkGluedMeasurementDet::checkHitProjection(), VZeroFinder::checkTrackPair(), TkGluedMeasurementDet::collectRecHits(), PixelClusterShapeSeedComparitor::compatible(), StripSubClusterShapeSeedFilter::compatible(), AnalyticalCurvilinearJacobian::computeFullJacobian(), AnalyticalCurvilinearJacobian::computeInfinitesimalJacobian(), AnalyticalCurvilinearJacobian::computeStraightLineJacobian(), NuclearTrackCorrector::getNewTrackExtra(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), SeedForPhotonConversion1Leg::initialKinematic(), SeedFromConsecutiveHitsCreator::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), JacobianCurvilinearToLocal::JacobianCurvilinearToLocal(), PerigeeConversions::jacobianPerigee2Curvilinear(), SeedFromConsecutiveHitsCreator::makeSeed(), TkClonerImpl::makeShared(), FreeTrajectoryState::momentum(), ClosestApproachInRPhi::newTrajectory(), TwoTrackMinimumDistance::pointsHelixHelix(), ConvBremSeedProducer::produce(), TkGluedMeasurementDet::projectOnGluedDet(), AnalyticalPropagator::propagatedStateWithPath(), KfTrackProducerBase::putInEvt(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), TkGluedMeasurementDet::testStrips(), TangentApproachInRPhi::trajectoryParameters(), SeedForPhotonConversion1Leg::trajectorySeed(), and TwoTrackMinimumDistanceHelixHelix::updateCoeffs().

                                {
    return theP;
  }

GlobalPoint GlobalTrajectoryParameters::position ( ) const

inline

Global position.

Definition at line 71 of file GlobalTrajectoryParameters.h.

References theX.

Referenced by ReferenceTrajectory::addMaterialEffectsBrl(), analyticalErrorPropagation(), MuonTrackLoader::buildTrackExtra(), TwoTrackMinimumDistanceHelixHelix::calculate(), TwoTrackMinimumDistanceLineLine::calculate(), TwoTrackMinimumDistance::calculate(), PixelClusterShapeSeedComparitor::compatible(), StripSubClusterShapeSeedFilter::compatible(), PixelHitMatcher::compatibleSeeds(), AnalyticalCurvilinearJacobian::computeFullJacobian(), ConversionForwardEstimator::estimate(), ConversionBarrelEstimator::estimate(), ForwardMeasurementEstimator::estimate(), BarrelMeasurementEstimator::estimate(), AnalyticalImpactPointExtrapolator::extrapolateSingleState(), AnalyticalTrajectoryExtrapolatorToLine::extrapolateSingleState(), TwoTrackMinimumDistanceHelixHelix::finalPoints(), NuclearTrackCorrector::getNewTrackExtra(), SeedFromConsecutiveHitsTripletOnlyCreator::initialKinematic(), SeedForPhotonConversion1Leg::initialKinematic(), SeedFromConsecutiveHitsCreator::initialKinematic(), SeedForPhotonConversionFromQuadruplets::initialKinematic(), ClosestApproachInRPhi::newTrajectory(), TwoTrackMinimumDistance::pointsHelixHelix(), FreeTrajectoryState::position(), ConvBremSeedProducer::produce(), AnalyticalPropagator::propagatedStateWithPath(), TauTagTools::propagTrackECALSurfContactPoint(), KfTrackProducerBase::putInEvt(), GsfTrackProducerBase::putInEvt(), TrackProducerWithSCAssociation::putInEvt(), TrajSeedMatcher::SCHitMatch::SCHitMatch(), TangentApproachInRPhi::trajectoryParameters(), and TwoTrackMinimumDistanceHelixHelix::updateCoeffs().

                               {
    return theX;
  }

void GlobalTrajectoryParameters::setCache ( )

private

Definition at line 35 of file GlobalTrajectoryParameters.cc.

References cachedMagneticField, MagneticField::inTesla(), theField, and theX.

Referenced by GlobalTrajectoryParameters(), and magneticField().

                                          {
    cachedMagneticField = theField ? theField->inTesla(theX) : GlobalVector(0, 0, 0);
} // we must initialize cache to non-NAN to avoid FPE

float GlobalTrajectoryParameters::signedInverseMomentum ( ) const

inline

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

Definition at line 96 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 103 of file GlobalTrajectoryParameters.h.

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

Referenced by transverseCurvature().

                                                {
    return theCharge/theP.perp();
  }

float GlobalTrajectoryParameters::transverseCurvature ( ) const

inline

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 112 of file GlobalTrajectoryParameters.h.

References cachedMagneticField, f, signedInverseTransverseMomentum(), and PV3DBase< T, PVType, FrameType >::z().

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

                                    {
    return -2.99792458e-3f * signedInverseTransverseMomentum() * cachedMagneticField.z();
  }

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 120 of file GlobalTrajectoryParameters.h.

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

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

                                  {
    return AlgebraicVector6(theX.x(),theX.y(),theX.z(),theP.x(),theP.y(),theP.z());
  }

Member Data Documentation

GlobalVector GlobalTrajectoryParameters::cachedMagneticField

private

Definition at line 143 of file GlobalTrajectoryParameters.h.

Referenced by GlobalTrajectoryParameters(), magneticFieldInInverseGeV(), magneticFieldInTesla(), setCache(), and transverseCurvature().

signed char GlobalTrajectoryParameters::theCharge

private

Definition at line 144 of file GlobalTrajectoryParameters.h.

Referenced by charge(), GlobalTrajectoryParameters(), signedInverseMomentum(), and signedInverseTransverseMomentum().

const MagneticField* GlobalTrajectoryParameters::theField

private

Definition at line 140 of file GlobalTrajectoryParameters.h.

Referenced by GlobalTrajectoryParameters(), magneticField(), magneticFieldInInverseGeV(), and setCache().

GlobalVector GlobalTrajectoryParameters::theP

private

Definition at line 142 of file GlobalTrajectoryParameters.h.

Referenced by direction(), GlobalTrajectoryParameters(), momentum(), signedInverseMomentum(), signedInverseTransverseMomentum(), and vector().

GlobalPoint GlobalTrajectoryParameters::theX

private

Definition at line 141 of file GlobalTrajectoryParameters.h.

Referenced by GlobalTrajectoryParameters(), position(), setCache(), and vector().