#include <LocalTrajectoryParameters.h>

Public Member Functions
TrackCharge	charge () const
	Charge (-1, 0 or 1)
LocalVector	direction () const
	Momentum vector unit in the local frame.
double	dxdz () const
double	dydz () const
	LocalTrajectoryParameters (const LocalPoint &pos, const LocalVector &p, TrackCharge charge)
	Constructor from local position, momentum and charge.
	LocalTrajectoryParameters (const AlgebraicVector5 &v, double aPzSign, bool charged=true)
	LocalTrajectoryParameters (double aQbp, double aDxdz, double aDydz, double aX, double aY, double aPzSign, bool charged=true)
	LocalTrajectoryParameters ()
AlgebraicVector5	mixedFormatVector () const
LocalVector	momentum () const
	Momentum vector in the local frame.
LocalPoint	position () const
	Local x and y position coordinates.
double	pzSign () const
	Sign of the z-component of the momentum in the local frame.
double	qbp () const
double	signedInverseMomentum () const
	Signed inverse momentum q/p (zero for neutrals).
bool	updateP (double dP)
	Update of momentum by a scalar dP.
AlgebraicVector5	vector () const
Private Attributes
TrackCharge	theCharge
	charge
double	theDxdz
	tangent of direction in local x vs. z
double	theDydz
	tangent of direction in local y vs. z
float	thePzSign
	sign of local pz
double	theQbp
	q/p (charged) or 1/p (neutral)
double	theX
	local x position
double	theY
	local y position

Detailed Description

Class providing access to a set of relevant parameters of a trajectory in a local, Cartesian frame. The set consists of the following parameters:

q/p : charged particles: charge (plus or minus one) divided by magnitude of momentum
neutral particles: inverse magnitude of momentum
dxdz : direction tangent in local xz-plane
dydz : direction tangent in local yz-plane
x : local x-coordinate
y : local y-coordinate

In addition, the sign of local p_z is needed to fully define the direction of the track in this local frame.

Definition at line 25 of file LocalTrajectoryParameters.h.

Constructor & Destructor Documentation

LocalTrajectoryParameters::LocalTrajectoryParameters ( ) [inline]

Definition at line 29 of file LocalTrajectoryParameters.h.

{}

LocalTrajectoryParameters::LocalTrajectoryParameters	(	const AlgebraicVector5 &	v,
		double	aPzSign,
		bool	charged = `true`
	)		`[inline]`

Constructor from vector of parameters.

Expects a vector of parameters as defined above, plus the sign of p_z. For charged particles the charge will be determined by the sign of the first element. For neutral particles the last argument should be false, in which case the charge of the first element will be neglected.

Definition at line 38 of file LocalTrajectoryParameters.h.

References theCharge, theDxdz, theDydz, thePzSign, theQbp, theX, and theY.

                                                                                            {
    theQbp  = v[0];
    theDxdz = v[1];
    theDydz = v[2];
    theX    = v[3];
    theY    = v[4];
    thePzSign = aPzSign;
    if ( charged )
      theCharge = theQbp>0 ? 1 : -1;
    else
      theCharge = 0;
  }

LocalTrajectoryParameters::LocalTrajectoryParameters	(	double	aQbp,
		double	aDxdz,
		double	aDydz,
		double	aX,
		double	aY,
		double	aPzSign,
		bool	charged = `true`
	)		`[inline]`

Constructor from individual parameters.

Expects parameters as defined above, plus the sign of p_z. For charged particles the charge will be determined by the sign of the first argument. For neutral particles the last argument should be false, in which case the charge of the first argument will be neglected.

Definition at line 59 of file LocalTrajectoryParameters.h.

References theCharge, and theQbp.

                                                                                       :
    theDxdz(aDxdz), theDydz(aDydz),
    theX(aX), theY(aY), thePzSign(aPzSign) {
    if ( charged ) {
      theQbp = aQbp;
      theCharge = theQbp>0 ? 1 : -1;
    }
    else {
      theQbp = aQbp;
      theCharge = 0;
    }
  }

LocalTrajectoryParameters::LocalTrajectoryParameters	(	const LocalPoint &	pos,
		const LocalVector &	p,
		TrackCharge	charge
	)		`[inline]`

Constructor from local position, momentum and charge.

Definition at line 74 of file LocalTrajectoryParameters.h.

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

                                                 :
    theQbp( charge/p.mag()), theDxdz( p.x()/p.z()), theDydz( p.y()/p.z()), 
    theX( pos.x()), theY(pos.y()), thePzSign( p.z()>0. ? 1.:-1.), theCharge(charge) {
    if ( charge==0 )  theQbp = 1./p.mag();
  }

Member Function Documentation

TrackCharge LocalTrajectoryParameters::charge ( ) const [inline]

Charge (-1, 0 or 1)

Definition at line 108 of file LocalTrajectoryParameters.h.

References theCharge.

Referenced by BzeroReferenceTrajectory::BzeroReferenceTrajectory(), BasicTrajectoryState::checkGlobalParameters(), display_seed(), reco::ElectronSeed::getCharge(), JacobianCartesianToLocal::JacobianCartesianToLocal(), JacobianLocalToCartesian::JacobianLocalToCartesian(), operator<<(), PrintRecoObjects::print(), FakeTrackProducer< T >::produce(), signedInverseMomentum(), and InOutConversionTrackFinder::tracks().

{return theCharge;}

LocalVector LocalTrajectoryParameters::direction ( ) const [inline]

Momentum vector unit in the local frame.

Definition at line 99 of file LocalTrajectoryParameters.h.

References mathSSE::sqrt(), theDxdz, theDydz, and thePzSign.

Referenced by JacobianCurvilinearToLocal::JacobianCurvilinearToLocal(), and JacobianLocalToCurvilinear::JacobianLocalToCurvilinear().

                                {
    double dz = thePzSign/sqrt(1. + theDxdz*theDxdz + theDydz*theDydz);
    double dx = dz*theDxdz;
    double dy = dz*theDydz;
    return LocalVector(dx, dy, dz);
  }

double LocalTrajectoryParameters::dxdz ( ) const [inline]

Definition at line 161 of file LocalTrajectoryParameters.h.

References theDxdz.

Referenced by TrackCandidateProducer::isDuplicateCandidate(), JacobianLocalToCartesian::JacobianLocalToCartesian(), PixelCPEBase::measurementError(), and PixelCPEBase::measurementPosition().

{ return theDxdz;}

double LocalTrajectoryParameters::dydz ( ) const [inline]

Definition at line 162 of file LocalTrajectoryParameters.h.

References theDydz.

Referenced by JacobianLocalToCartesian::JacobianLocalToCartesian(), PixelCPEBase::measurementError(), and PixelCPEBase::measurementPosition().

{ return theDydz;}

AlgebraicVector5 LocalTrajectoryParameters::mixedFormatVector ( ) const [inline]

Vector of parameters in internal representation.

Vector of parameters as defined above, with the first element = q/p for charged and = 1/p for neutral.

Definition at line 135 of file LocalTrajectoryParameters.h.

References theDxdz, theDydz, theQbp, theX, theY, and v.

Referenced by BzeroReferenceTrajectory::BzeroReferenceTrajectory(), PixelCPEBase::computeAnglesFromTrajectory(), TwoBowedSurfacesAlignmentParameters::derivatives(), DualReferenceTrajectory::extractParameters(), DualBzeroReferenceTrajectory::extractParameters(), MuonKinkFinder::getChi2(), MatcherUsingTracksAlgorithm::getChi2(), SegmentAlignmentDerivatives4D::operator()(), BowedSurfaceAlignmentDerivatives::operator()(), KarimakiAlignmentDerivatives::operator()(), operator<<(), and ReferenceTrajectory::ReferenceTrajectory().

                                             {
    AlgebraicVector5 v;
    v[0] = theQbp;    // signed in case of charged particles, 1/p for neutrals
    v[1] = theDxdz;
    v[2] = theDydz;
    v[3] = theX;
    v[4] = theY;
    return v;
  }

LocalVector LocalTrajectoryParameters::momentum ( ) const [inline]

Momentum vector in the local frame.

Definition at line 89 of file LocalTrajectoryParameters.h.

References alignCSCRings::e, mathSSE::sqrt(), theDxdz, theDydz, thePzSign, and theQbp.

Referenced by SiPixelTrackResidualSource::analyze(), SiPixelErrorEstimation::analyze(), BasicTrajectoryState::checkGlobalParameters(), MuonSeedSimpleCleaner::checkPt(), PixelCPEBase::computeAnglesFromTrajectory(), TwoBodyDecayTrajectory::constructSingleTsosWithErrors(), display_seed(), SiPixelHitEfficiencyModule::fill(), JacobianCartesianToLocal::JacobianCartesianToLocal(), BasicTrajectoryState::localMomentum(), StripCPEgeometric::localParameters(), StripCPEfromTemplate::localParameters(), StripCPEfromTrackAngle::localParameters(), reco::TrackInfo::localTrackMomentum(), PrintRecoObjects::print(), TrackCandidateProducer::produce(), KalmanAlignmentTrackRefitter::refitSingleTracklet(), TrackInfoProducerAlgorithm::run(), TrackClusterSplitter::splitCluster(), and InOutConversionTrackFinder::tracks().

                               {
    double op = fabs(theQbp);
    if ( op<1.e-9 )  op = 1.e-9;
    double pz = thePzSign/(op*sqrt(1. + theDxdz*theDxdz + theDydz*theDydz));
    double px = pz*theDxdz;
    double py = pz*theDydz;
    return LocalVector(px, py, pz);
  }

LocalPoint LocalTrajectoryParameters::position ( ) const [inline]

double LocalTrajectoryParameters::pzSign ( ) const [inline]

Sign of the z-component of the momentum in the local frame.

Definition at line 146 of file LocalTrajectoryParameters.h.

References thePzSign.

Referenced by BzeroReferenceTrajectory::BzeroReferenceTrajectory(), TrajectoryStateCombiner::combine(), JacobianLocalToCartesian::JacobianLocalToCartesian(), GlobalTrackerMuonAlignment::misalignMuonL(), MultiTrajectoryStateMode::momentumFromModeLocal(), GsfTrackProducerBase::putInEvt(), GaussianStateConversions::tsosFromMultiGaussianState(), MultiGaussianStateTransform::tsosFromSingleState(), KFStripUpdator::update(), and KFStrip1DUpdator::update().

                        {
    return thePzSign;
  }

double LocalTrajectoryParameters::qbp ( ) const [inline]

Definition at line 160 of file LocalTrajectoryParameters.h.

References theQbp.

Referenced by TrackCandidateProducer::isDuplicateCandidate(), and JacobianLocalToCartesian::JacobianLocalToCartesian().

{ return theQbp;}

double LocalTrajectoryParameters::signedInverseMomentum ( ) const [inline]

Signed inverse momentum q/p (zero for neutrals).

Definition at line 111 of file LocalTrajectoryParameters.h.

References charge(), and theQbp.

Referenced by ReferenceTrajectory::construct(), TwoBodyDecayTrajectory::constructSingleTsosWithErrors(), JacobianCurvilinearToLocal::JacobianCurvilinearToLocal(), JacobianLocalToCurvilinear::JacobianLocalToCurvilinear(), NuclearInteractionFinder::rescaleError(), and vector().

                                       {
    return charge()==0 ? 0. : theQbp;
  }

bool LocalTrajectoryParameters::updateP ( double dP ) [inline]

Update of momentum by a scalar dP.

Definition at line 151 of file LocalTrajectoryParameters.h.

References AlCaHLTBitMon_ParallelJobs::p, and theQbp.

Referenced by VolumeMaterialEffectsUpdator::updateState(), GsfMaterialEffectsUpdator::updateState(), and MaterialEffectsUpdator::updateStateInPlace().

                          {
    double p = 1./fabs(theQbp);
    if ((p += dP) <= 0.) return false;
    double newQbp = theQbp > 0. ? 1./p : -1./p;
    theQbp = newQbp;
    return true;
  }

AlgebraicVector5 LocalTrajectoryParameters::vector ( ) const [inline]

Vector of parameters with signed inverse momentum.

Vector of parameters as defined above, with the first element = q/p .

Definition at line 120 of file LocalTrajectoryParameters.h.

References signedInverseMomentum(), theDxdz, theDydz, theX, theY, and v.

Referenced by GlobalTrackerMuonAlignment::analyzeTrackTrack(), GlobalTrackerMuonAlignment::analyzeTrackTrajectory(), CollinearFitAtTM2::CollinearFitAtTM2(), TrajectoryStateCombiner::combine(), GsfTrackProducerBase::computeModeAtTM(), PFGsfHelper::computeQpMode(), GlobalTrackerMuonAlignment::debugTrajectorySOS(), GlobalTrackerMuonAlignment::debugTrajectorySOSv(), Chi2MeasurementEstimatorForTrackerHits::estimate(), CollinearFitAtTM::fit(), StripCPEgeometric::localParameters(), StripCPEfromTemplate::localParameters(), StripCPEfromTrackAngle::localParameters(), GsfTrackProducerBase::localParametersFromQpMode(), GlobalMuonTrackMatcher::match_Chi2(), MeasurementExtractor::measuredParameters(), GlobalTrackerMuonAlignment::misalignMuonL(), TrackCandidateProducer::produce(), TrackAssociatorByPosition::quality(), ChargeSignificanceTrajectoryFilter::TBC(), CkfDebugger::testSeed(), Strip1DMeasurementTransformator::trajectoryParameters(), StripMeasurementTransformator::trajectoryParameters(), Tsos2DPhi::Tsos2DPhi(), Tsos2DZed::Tsos2DZed(), and Tsos4D::Tsos4D().

                                  {
    AlgebraicVector5 v;
    v[0] = signedInverseMomentum();
    v[1] = theDxdz;
    v[2] = theDydz;
    v[3] = theX;
    v[4] = theY;
    return v;
  }