#include <PerigeeLinearizedTrackState.h>

Inheritance diagram for PerigeeLinearizedTrackState:

Public Types
typedef ReferenceCountingPointer< LinearizedTrackState< 5 > >	RefCountedLinearizedTrackState

Public Types inherited from LinearizedTrackState< 5 >
typedef ROOT::Math::SMatrix< double, N - 2, 3, ROOT::Math::MatRepStd< double, N - 2, 3 > >	AlgebraicMatrixM3

typedef ROOT::Math::SMatrix< double, N, 3, ROOT::Math::MatRepStd< double, N, 3 > >	AlgebraicMatrixN3

typedef ROOT::Math::SMatrix< double, N, N - 2, ROOT::Math::MatRepStd< double, N, N - 2 > >	AlgebraicMatrixNM

typedef ROOT::Math::SMatrix< double, N - 2, N - 2, ROOT::Math::MatRepSym< double, N - 2 > >	AlgebraicSymMatrixMM

typedef ROOT::Math::SMatrix< double, N, N, ROOT::Math::MatRepSym< double, N > >	AlgebraicSymMatrixNN

typedef ROOT::Math::SMatrix< double, N+1, N+1, ROOT::Math::MatRepSym< double, N+1 > >	AlgebraicSymMatrixOO

typedef ROOT::Math::SVector< double, N - 2 >	AlgebraicVectorM

typedef ROOT::Math::SVector< double, N >	AlgebraicVectorN

typedef ReferenceCountingPointer< RefittedTrackState< N > >	RefCountedRefittedTrackState

Public Member Functions
void	checkParameters (AlgebraicVector5 &parameters) const override

std::vector< ReferenceCountingPointer< LinearizedTrackState< 5 > > >	components () const override

const AlgebraicVector5 &	constantTerm () const override

RefCountedRefittedTrackState	createRefittedTrackState (const GlobalPoint &vertexPosition, const AlgebraicVector3 &vectorParameters, const AlgebraicSymMatrix66 &covarianceMatrix) const override

bool	isValid () const override

const GlobalPoint &	linearizationPoint () const override

const AlgebraicMatrix53 &	momentumJacobian () const override

const AlgebraicVector5 &	parametersFromExpansion () const override

const AlgebraicMatrix53 &	positionJacobian () const override

const TrajectoryStateClosestToPoint &	predictedState () const

AlgebraicSymMatrix55	predictedStateError () const override

AlgebraicSymMatrix33	predictedStateMomentumError () const override

AlgebraicVector3	predictedStateMomentumParameters () const override

AlgebraicVector5	predictedStateParameters () const override

AlgebraicSymMatrix55	predictedStateWeight (int &error) const override

AlgebraicVector5	refittedParamFromEquation (const RefCountedRefittedTrackState &theRefittedState) const override

const TrajectoryStateOnSurface	state () const

RefCountedLinearizedTrackState	stateWithNewLinearizationPoint (const GlobalPoint &newLP) const override

reco::TransientTrack	track () const override

double	weightInMixture () const override

Public Member Functions inherited from LinearizedTrackState< 5 >
virtual TrackCharge	charge () const=0

virtual void	checkParameters (AlgebraicVectorN &parameters) const

virtual RefCountedRefittedTrackState	createRefittedTrackState (const GlobalPoint &vertexPosition, const AlgebraicVectorM &vectorParameters, const AlgebraicSymMatrixOO &covarianceMatrix) const=0

virtual bool	hasError () const=0

virtual bool	operator== (LinearizedTrackState< N > &other) const=0

virtual AlgebraicVectorN	refittedParamFromEquation (const RefCountedRefittedTrackState &theRefittedState) const=0

	~LinearizedTrackState () override

Private Member Functions
void	compute3DImpactPoint () const

void	computeChargedJacobians () const

void	computeJacobians () const

void	computeNeutralJacobians () const

	PerigeeLinearizedTrackState (const GlobalPoint &linP, const reco::TransientTrack &track, const TrajectoryStateOnSurface &tsos)

Friends
class	LinearizedTrackStateFactory

Detailed Description

Calculates and stores the ImpactPointMeasurement of the impact point (point of closest approach in 3D) to the given linearization point. (see V.Karimaki, HIP-1997-77 / EXP)

Computes the parameters of the trajectory state at the transverse point of closest approach (in the global transverse plane) to the linearization point, and the jacobiam matrices. (see R.Fruehwirth et al. Data Analysis Techniques in HEP Experiments Second Edition, Cambridge University Press 2000, or R.Fruehwirth et al. Vertex reconstruction and track bundling at the LEP collider using robust algorithms. Computer Physics Communications 96 (1996) 189-208).

Both are done ‘on-demand’ to improve CPU performance.

This particular implementation works with "perigee" helix parametrization: see Billoir et al. NIM in PR A311(1992) 139-150

Definition at line 33 of file PerigeeLinearizedTrackState.h.

Member Typedef Documentation

◆ RefCountedLinearizedTrackState

typedef ReferenceCountingPointer<LinearizedTrackState<5> > PerigeeLinearizedTrackState::RefCountedLinearizedTrackState

Definition at line 39 of file PerigeeLinearizedTrackState.h.

Constructor & Destructor Documentation

◆ PerigeeLinearizedTrackState()

PerigeeLinearizedTrackState::PerigeeLinearizedTrackState	(	const GlobalPoint &	linP,
		const reco::TransientTrack &	track,
		const TrajectoryStateOnSurface &	tsos
	)

inlineprivate

Constructor with the linearization point and the track. Private, can only be used by LinearizedTrackFactory.

Definition at line 140 of file PerigeeLinearizedTrackState.h.

143 : theLinPoint(linP), theTrack(track), theTSOS(tsos), theCharge(theTrack.charge()), jacobiansAvailable(false) {}

Referenced by stateWithNewLinearizationPoint().

Member Function Documentation

◆ checkParameters()

void PerigeeLinearizedTrackState::checkParameters ( AlgebraicVector5 & parameters ) const

inlineoverride

Definition at line 87 of file PerigeeLinearizedTrackState.cc.

                                                                                     {
   if (parameters(2) > M_PI)
     parameters(2) -= 2 * M_PI;
   if (parameters(2) < -M_PI)
     parameters(2) += 2 * M_PI;
 }

References M_PI.

◆ components()

std::vector< PerigeeLinearizedTrackState::RefCountedLinearizedTrackState > PerigeeLinearizedTrackState::components ( ) const

overridevirtual

Implements LinearizedTrackState< 5 >.

Definition at line 59 of file PerigeeLinearizedTrackState.cc.

           {
   std::vector<RefCountedLinearizedTrackState> result;
   result.reserve(1);
   result.push_back(RefCountedLinearizedTrackState(const_cast<PerigeeLinearizedTrackState*>(this)));
   return result;
 }

References mps_fire::result.

◆ compute3DImpactPoint()

void PerigeeLinearizedTrackState::compute3DImpactPoint ( ) const

private

Method calculating the 3D impact point measurement

◆ computeChargedJacobians()

void PerigeeLinearizedTrackState::computeChargedJacobians ( ) const

private

Method calculating the track parameters and the Jacobians for charged particles.

Definition at line 94 of file PerigeeLinearizedTrackState.cc.

                                                                 {
   GlobalPoint paramPt(theLinPoint);
   //tarjectory parameters
   double field = theTrack.field()->inInverseGeV(thePredState.theState().position()).z();
   double signTC = -theCharge;
  
   double thetaAtEP = thePredState.perigeeParameters().theta();
   double phiAtEP = thePredState.perigeeParameters().phi();
   double ptAtEP = thePredState.pt();
   double transverseCurvatureAtEP = field / ptAtEP * signTC;
  
   double x_v = thePredState.theState().position().x();
   double y_v = thePredState.theState().position().y();
   double z_v = thePredState.theState().position().z();
   double X = x_v - paramPt.x() - sin(phiAtEP) / transverseCurvatureAtEP;
   double Y = y_v - paramPt.y() + cos(phiAtEP) / transverseCurvatureAtEP;
   double SS = X * X + Y * Y;
   double S = sqrt(SS);
  
   // The track parameters at the expansion point
  
   theExpandedParams[0] = transverseCurvatureAtEP;
   theExpandedParams[1] = thetaAtEP;
   theExpandedParams[3] = 1 / transverseCurvatureAtEP - signTC * S;
   double phiFEP;
   if (std::abs(X) > std::abs(Y)) {
     double signX = (X > 0.0 ? +1.0 : -1.0);
     phiFEP = -signTC * signX * acos(signTC * Y / S);
   } else {
     phiFEP = asin(-signTC * X / S);
     if ((signTC * Y) < 0.0)
       phiFEP = M_PI - phiFEP;
   }
   if (phiFEP > M_PI)
     phiFEP -= 2 * M_PI;
   theExpandedParams[2] = phiFEP;
   theExpandedParams[4] =
       z_v - paramPt.z() - (phiAtEP - theExpandedParams[2]) / tan(thetaAtEP) / transverseCurvatureAtEP;
  
   // The Jacobian: (all at the expansion point)
   // [i,j]
   // i = 0: rho , 1: theta, 2: phi_p, 3: epsilon, 4: z_p
   // j = 0: x_v, 1: y_v, 2: z_v
  
   thePositionJacobian(2, 0) = -Y / (SS);
   thePositionJacobian(2, 1) = X / (SS);
   thePositionJacobian(3, 0) = -signTC * X / S;
   thePositionJacobian(3, 1) = -signTC * Y / S;
   thePositionJacobian(4, 0) = thePositionJacobian(2, 0) / tan(thetaAtEP) / transverseCurvatureAtEP;
   thePositionJacobian(4, 1) = thePositionJacobian(2, 1) / tan(thetaAtEP) / transverseCurvatureAtEP;
   thePositionJacobian(4, 2) = 1;
  
   // [i,j]
   // i = 0: rho , 1: theta, 2: phi_p, 3: epsilon, 4: z_p
   // j = 0: rho, 1: theta, 2: phi_v
   theMomentumJacobian(0, 0) = 1;
   theMomentumJacobian(1, 1) = 1;
  
   theMomentumJacobian(2, 0) =
       -(X * cos(phiAtEP) + Y * sin(phiAtEP)) / (SS * transverseCurvatureAtEP * transverseCurvatureAtEP);
  
   theMomentumJacobian(2, 2) = (Y * cos(phiAtEP) - X * sin(phiAtEP)) / (SS * transverseCurvatureAtEP);
  
   theMomentumJacobian(3, 0) =
       (signTC * (Y * cos(phiAtEP) - X * sin(phiAtEP)) / S - 1) / (transverseCurvatureAtEP * transverseCurvatureAtEP);
  
   theMomentumJacobian(3, 2) = signTC * (X * cos(phiAtEP) + Y * sin(phiAtEP)) / (S * transverseCurvatureAtEP);
  
   theMomentumJacobian(4, 0) =
       (phiAtEP - theExpandedParams[2]) / tan(thetaAtEP) / (transverseCurvatureAtEP * transverseCurvatureAtEP) +
       theMomentumJacobian(2, 0) / tan(thetaAtEP) / transverseCurvatureAtEP;
  
   theMomentumJacobian(4, 1) =
       (phiAtEP - theExpandedParams[2]) * (1 + 1 / (tan(thetaAtEP) * tan(thetaAtEP))) / transverseCurvatureAtEP;
  
   theMomentumJacobian(4, 2) = (theMomentumJacobian(2, 2) - 1) / tan(thetaAtEP) / transverseCurvatureAtEP;
  
   // And finally the residuals:
  
   auto p = thePredState.theState().position();
   AlgebraicVector3 expansionPoint(p.x(), p.y(), p.z());
   AlgebraicVector3 momentumAtExpansionPoint(transverseCurvatureAtEP, thetaAtEP, phiAtEP);
  
   theConstantTerm = AlgebraicVector5(theExpandedParams - thePositionJacobian * expansionPoint -
                                      theMomentumJacobian * momentumAtExpansionPoint);
 }

References funct::abs(), funct::cos(), reco::TransientTrack::field(), MagneticField::inInverseGeV(), M_PI, AlCaHLTBitMon_ParallelJobs::p, TrajectoryStateClosestToPoint::perigeeParameters(), PerigeeTrajectoryParameters::phi(), FreeTrajectoryState::position(), TrajectoryStateClosestToPoint::pt(), S(), funct::sin(), mathSSE::sqrt(), funct::tan(), theCharge, theConstantTerm, theExpandedParams, theLinPoint, theMomentumJacobian, thePositionJacobian, thePredState, TrajectoryStateClosestToPoint::theState(), PerigeeTrajectoryParameters::theta(), theTrack, X, PV3DBase< T, PVType, FrameType >::x(), DOFs::Y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by computeJacobians().

◆ computeJacobians()

void PerigeeLinearizedTrackState::computeJacobians ( ) const

private

Method calculating the track parameters and the Jacobians.

Definition at line 9 of file PerigeeLinearizedTrackState.cc.

                                                          {
   GlobalPoint paramPt(theLinPoint);
  
   thePredState = builder(theTSOS, paramPt);
   if
     UNLIKELY(!thePredState.isValid()) return;
  
   double field = theTrack.field()->inInverseGeV(thePredState.theState().position()).z();
  
   if ((std::abs(theCharge) < 1e-5) || (fabs(field) < 1.e-10)) {
     //neutral track
     computeNeutralJacobians();
   } else {
     //charged track
     computeChargedJacobians();
   }
  
   jacobiansAvailable = true;
 }

References funct::abs(), builder, computeChargedJacobians(), computeNeutralJacobians(), MillePedeFileConverter_cfg::e, reco::TransientTrack::field(), MagneticField::inInverseGeV(), TrajectoryStateClosestToPoint::isValid(), jacobiansAvailable, FreeTrajectoryState::position(), theCharge, theLinPoint, thePredState, TrajectoryStateClosestToPoint::theState(), theTrack, theTSOS, UNLIKELY, and z.

Referenced by constantTerm(), isValid(), momentumJacobian(), parametersFromExpansion(), positionJacobian(), predictedState(), predictedStateError(), predictedStateMomentumError(), predictedStateMomentumParameters(), predictedStateParameters(), and predictedStateWeight().

◆ computeNeutralJacobians()

void PerigeeLinearizedTrackState::computeNeutralJacobians ( ) const

private

Method calculating the track parameters and the Jacobians for neutral particles.

Definition at line 181 of file PerigeeLinearizedTrackState.cc.

                                                                 {
   GlobalPoint paramPt(theLinPoint);
  
   //tarjectory parameters
   double thetaAtEP = thePredState.theState().momentum().theta();
   double phiAtEP = thePredState.theState().momentum().phi();
   double ptAtEP = thePredState.theState().momentum().perp();
  
   double x_v = thePredState.theState().position().x();
   double y_v = thePredState.theState().position().y();
   double z_v = thePredState.theState().position().z();
   double X = x_v - paramPt.x();
   double Y = y_v - paramPt.y();
  
   // The track parameters at the expansion point
  
   theExpandedParams(0) = 1 / ptAtEP;
   theExpandedParams(1) = thetaAtEP;
   theExpandedParams(2) = phiAtEP;
   theExpandedParams(3) = X * sin(phiAtEP) - Y * cos(phiAtEP);
   theExpandedParams(4) = z_v - paramPt.z() - (X * cos(phiAtEP) + Y * sin(phiAtEP)) / tan(thetaAtEP);
  
   // The Jacobian: (all at the expansion point)
   // [i,j]
   // i = 0: rho = 1/pt , 1: theta, 2: phi_p, 3: epsilon, 4: z_p
   // j = 0: x_v, 1: y_v, 2: z_v
  
   thePositionJacobian(3, 0) = sin(phiAtEP);
   thePositionJacobian(3, 1) = -cos(phiAtEP);
   thePositionJacobian(4, 0) = -cos(phiAtEP) / tan(thetaAtEP);
   thePositionJacobian(4, 1) = -sin(phiAtEP) / tan(thetaAtEP);
   thePositionJacobian(4, 2) = 1;
  
   // [i,j]
   // i = 0: rho = 1/pt , 1: theta, 2: phi_p, 3: epsilon, 4: z_p
   // j = 0: rho = 1/pt , 1: theta, 2: phi_v
  
   theMomentumJacobian(0, 0) = 1;
   theMomentumJacobian(1, 1) = 1;
   theMomentumJacobian(2, 2) = 1;
  
   theMomentumJacobian(3, 2) = X * cos(phiAtEP) + Y * sin(phiAtEP);
  
   theMomentumJacobian(4, 1) = theMomentumJacobian(3, 2) * (1 + 1 / (tan(thetaAtEP) * tan(thetaAtEP)));
  
   theMomentumJacobian(4, 2) = (X * sin(phiAtEP) - Y * cos(phiAtEP)) / tan(thetaAtEP);
  
   // And finally the residuals:
  
   auto p = thePredState.theState().position();
   AlgebraicVector3 expansionPoint(p.x(), p.y(), p.z());
   AlgebraicVector3 momentumAtExpansionPoint(1. / ptAtEP, thetaAtEP, phiAtEP);
  
   theConstantTerm = AlgebraicVector5(theExpandedParams - thePositionJacobian * expansionPoint -
                                      theMomentumJacobian * momentumAtExpansionPoint);
 }

References funct::cos(), FreeTrajectoryState::momentum(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), FreeTrajectoryState::position(), funct::sin(), funct::tan(), theConstantTerm, theExpandedParams, theLinPoint, theMomentumJacobian, thePositionJacobian, thePredState, TrajectoryStateClosestToPoint::theState(), PV3DBase< T, PVType, FrameType >::theta(), X, PV3DBase< T, PVType, FrameType >::x(), DOFs::Y, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by computeJacobians().

◆ constantTerm()

const AlgebraicVector5 & PerigeeLinearizedTrackState::constantTerm ( ) const

inlineoverridevirtual

Method returning the constant term of the Taylor expansion of the measurement equation

Implements LinearizedTrackState< 5 >.

Definition at line 176 of file PerigeeLinearizedTrackState.h.

                                                                                {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return theConstantTerm;
 }

References computeJacobians(), jacobiansAvailable, theConstantTerm, and UNLIKELY.

Referenced by TwoBodyDecayEstimator::constructMatrices(), and refittedParamFromEquation().

◆ createRefittedTrackState()

PerigeeLinearizedTrackState::RefCountedRefittedTrackState PerigeeLinearizedTrackState::createRefittedTrackState	(	const GlobalPoint &	vertexPosition,
		const AlgebraicVector3 &	vectorParameters,
		const AlgebraicSymMatrix66 &	covarianceMatrix
	)		const

override

Method returning the impact point measurement ‍/ ImpactPointMeasurement impactPointMeasurement() const;

TrackCharge charge() const override { return theCharge; }

bool hasError() const override;

bool operator==(LinearizedTrackState<5>& other) const override;

bool operator==(ReferenceCountingPointer<LinearizedTrackState<5> >& other) const;

/** Creates the correct refitted state according to the results of the track refit.

Definition at line 50 of file PerigeeLinearizedTrackState.cc.

                                                         {
   TrajectoryStateClosestToPoint refittedTSCP = PerigeeConversions::trajectoryStateClosestToPoint(
       vectorParameters, vertexPosition, charge(), covarianceMatrix, theTrack.field());
   return RefCountedRefittedTrackState(new PerigeeRefittedTrackState(refittedTSCP, vectorParameters));
 }

References LinearizedTrackState< 5 >::charge(), reco::TransientTrack::field(), theTrack, and PerigeeConversions::trajectoryStateClosestToPoint().

◆ isValid()

bool PerigeeLinearizedTrackState::isValid ( void ) const

inlineoverridevirtual

Reimplemented from LinearizedTrackState< 5 >.

Definition at line 259 of file PerigeeLinearizedTrackState.h.

                                                        {
   if
     UNLIKELY(!theTSOS.isValid()) return false;
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return jacobiansAvailable;
 }

References computeJacobians(), TrajectoryStateOnSurface::isValid(), jacobiansAvailable, theTSOS, and UNLIKELY.

Referenced by ntupleDataFormat._Object::_checkIsValid(), and core.AutoHandle.AutoHandle::ReallyLoad().

◆ linearizationPoint()

const GlobalPoint& PerigeeLinearizedTrackState::linearizationPoint ( ) const

inlineoverridevirtual

The point at which the track state has been linearized

Implements LinearizedTrackState< 5 >.

Definition at line 50 of file PerigeeLinearizedTrackState.h.

50 { return theLinPoint; }

References theLinPoint.

Referenced by TwoBodyDecayEstimator::constructMatrices().

◆ momentumJacobian()

const AlgebraicMatrix53 & PerigeeLinearizedTrackState::momentumJacobian ( ) const

inlineoverridevirtual

Method returning the Momentum Jacobian from the Taylor expansion (Matrix B)

Method returning the Momentum Jacobian (Matrix B)

Implements LinearizedTrackState< 5 >.

Definition at line 194 of file PerigeeLinearizedTrackState.h.

                                                                                     {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return theMomentumJacobian;
 }

References computeJacobians(), jacobiansAvailable, theMomentumJacobian, and UNLIKELY.

Referenced by TwoBodyDecayEstimator::constructMatrices(), and refittedParamFromEquation().

◆ parametersFromExpansion()

const AlgebraicVector5 & PerigeeLinearizedTrackState::parametersFromExpansion ( ) const

inlineoverridevirtual

Method returning the parameters of the Taylor expansion

Implements LinearizedTrackState< 5 >.

Definition at line 202 of file PerigeeLinearizedTrackState.h.

                                                                                           {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return theExpandedParams;
 }

References computeJacobians(), jacobiansAvailable, theExpandedParams, and UNLIKELY.

◆ positionJacobian()

const AlgebraicMatrix53 & PerigeeLinearizedTrackState::positionJacobian ( ) const

inlineoverridevirtual

Method returning the Position Jacobian from the Taylor expansion (Matrix A)

Method returning the Position Jacobian (Matrix A)

Implements LinearizedTrackState< 5 >.

Definition at line 185 of file PerigeeLinearizedTrackState.h.

                                                                                     {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return thePositionJacobian;
 }

References computeJacobians(), jacobiansAvailable, thePositionJacobian, and UNLIKELY.

Referenced by TwoBodyDecayEstimator::constructMatrices(), and refittedParamFromEquation().

◆ predictedState()

const TrajectoryStateClosestToPoint & PerigeeLinearizedTrackState::predictedState ( ) const

inline

Method returning the track state at the point of closest approach to the linearization point, in the transverse plane (a.k.a. transverse impact point).

Method returning the TrajectoryStateClosestToPoint at the point of closest approch to the z-axis (a.k.a. transverse impact point)

Definition at line 212 of file PerigeeLinearizedTrackState.h.

                                                                                               {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return thePredState;
 }

References computeJacobians(), jacobiansAvailable, thePredState, and UNLIKELY.

Referenced by TwoBodyDecayLinearizationPointFinder::getLinearizationPoint(), and PerigeeMultiLTS::predictedState().

◆ predictedStateError()

AlgebraicSymMatrix55 PerigeeLinearizedTrackState::predictedStateError ( ) const

inlineoverridevirtual

Method returning the covariance matrix of the track state at the transverse impact point.

Implements LinearizedTrackState< 5 >.

Definition at line 247 of file PerigeeLinearizedTrackState.h.

                                                                                    {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return thePredState.perigeeError().covarianceMatrix();
 }

References computeJacobians(), PerigeeTrajectoryError::covarianceMatrix(), jacobiansAvailable, TrajectoryStateClosestToPoint::perigeeError(), thePredState, and UNLIKELY.

Referenced by TwoBodyDecayEstimator::constructMatrices().

◆ predictedStateMomentumError()

AlgebraicSymMatrix33 PerigeeLinearizedTrackState::predictedStateMomentumError ( ) const

inlineoverridevirtual

Method returning the momentum covariance matrix of the track state at the transverse impact point.

Implements LinearizedTrackState< 5 >.

Definition at line 253 of file PerigeeLinearizedTrackState.h.

                                                                                            {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   return thePredState.perigeeError().covarianceMatrix().Sub<AlgebraicSymMatrix33>(0, 2);
 }

References computeJacobians(), PerigeeTrajectoryError::covarianceMatrix(), jacobiansAvailable, TrajectoryStateClosestToPoint::perigeeError(), thePredState, and UNLIKELY.

◆ predictedStateMomentumParameters()

AlgebraicVector3 PerigeeLinearizedTrackState::predictedStateMomentumParameters ( ) const

inlineoverridevirtual

Method returning the momentum part of the parameters of the track state at the linearization point.

Implements LinearizedTrackState< 5 >.

Definition at line 230 of file PerigeeLinearizedTrackState.h.

                                                                                             {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   auto v = thePredState.perigeeParameters().vector();
   return AlgebraicVector3(v[0], v[1], v[2]);
 }

References computeJacobians(), jacobiansAvailable, TrajectoryStateClosestToPoint::perigeeParameters(), thePredState, UNLIKELY, findQualityFiles::v, and PerigeeTrajectoryParameters::vector().

Referenced by TwoBodyDecayEstimator::constructMatrices(), and refittedParamFromEquation().

◆ predictedStateParameters()

AlgebraicVector5 PerigeeLinearizedTrackState::predictedStateParameters ( ) const

inlineoverridevirtual

Method returning the parameters of the track state at the transverse impact point.

Implements LinearizedTrackState< 5 >.

Definition at line 224 of file PerigeeLinearizedTrackState.h.

                                                                                     {
   if (!jacobiansAvailable)
     computeJacobians();
   return thePredState.perigeeParameters().vector();
 }

References computeJacobians(), jacobiansAvailable, TrajectoryStateClosestToPoint::perigeeParameters(), thePredState, and PerigeeTrajectoryParameters::vector().

Referenced by TwoBodyDecayEstimator::constructMatrices().

◆ predictedStateWeight()

AlgebraicSymMatrix55 PerigeeLinearizedTrackState::predictedStateWeight ( int & error ) const

inlineoverridevirtual

Method returning the weight matrix of the track state at the transverse impact point. The error variable is 0 in case of success.

Implements LinearizedTrackState< 5 >.

Definition at line 237 of file PerigeeLinearizedTrackState.h.

                                                                                               {
   if
     UNLIKELY(!jacobiansAvailable) computeJacobians();
   if (!thePredState.isValid()) {
     error = 1;
     return AlgebraicSymMatrix55();
   }
   return thePredState.perigeeError().weightMatrix(error);
 }

References computeJacobians(), relativeConstraints::error, TrajectoryStateClosestToPoint::isValid(), jacobiansAvailable, TrajectoryStateClosestToPoint::perigeeError(), thePredState, UNLIKELY, and PerigeeTrajectoryError::weightMatrix().

◆ refittedParamFromEquation()

AlgebraicVector5 PerigeeLinearizedTrackState::refittedParamFromEquation ( const RefCountedRefittedTrackState & theRefittedState ) const

override

Definition at line 67 of file PerigeeLinearizedTrackState.cc.

                                                                 {
   auto p = theRefittedState->position();
   AlgebraicVector3 vertexPosition(p.x(), p.y(), p.z());
   AlgebraicVector3 momentum = theRefittedState->momentumVector();
   if ((momentum(2) * predictedStateMomentumParameters()(2) < 0) && (fabs(momentum(2)) > M_PI / 2)) {
     if (predictedStateMomentumParameters()(2) < 0.)
       momentum(2) -= 2 * M_PI;
     if (predictedStateMomentumParameters()(2) > 0.)
       momentum(2) += 2 * M_PI;
   }
   AlgebraicVectorN param = constantTerm() + positionJacobian() * vertexPosition + momentumJacobian() * momentum;
   if (param(2) > M_PI)
     param(2) -= 2 * M_PI;
   if (param(2) < -M_PI)
     param(2) += 2 * M_PI;
  
   return param;
 }

References constantTerm(), M_PI, momentumJacobian(), AlCaHLTBitMon_ParallelJobs::p, positionJacobian(), and predictedStateMomentumParameters().

◆ state()

const TrajectoryStateOnSurface PerigeeLinearizedTrackState::state ( ) const

inline

Definition at line 54 of file PerigeeLinearizedTrackState.h.

54 { return theTSOS; }

References theTSOS.

◆ stateWithNewLinearizationPoint()

PerigeeLinearizedTrackState::RefCountedLinearizedTrackState PerigeeLinearizedTrackState::stateWithNewLinearizationPoint ( const GlobalPoint & newLP ) const

overridevirtual

Returns a new linearized state with respect to a new linearization point. A new object of the same type is returned, without change to the existing one.

Implements LinearizedTrackState< 5 >.

Definition at line 45 of file PerigeeLinearizedTrackState.cc.

                                     {
   return RefCountedLinearizedTrackState(new PerigeeLinearizedTrackState(newLP, track(), theTSOS));
 }

References PerigeeLinearizedTrackState(), theTSOS, and track().

◆ track()

reco::TransientTrack PerigeeLinearizedTrackState::track ( ) const

inlineoverridevirtual

Implements LinearizedTrackState< 5 >.

Definition at line 52 of file PerigeeLinearizedTrackState.h.

52 { return theTrack; }

References theTrack.

Referenced by TwoBodyDecayEstimator::constructMatrices(), stateWithNewLinearizationPoint(), and TrimmedVertexFinder::vertices().

◆ weightInMixture()

double PerigeeLinearizedTrackState::weightInMixture ( ) const

inlineoverridevirtual

Implements LinearizedTrackState< 5 >.

Definition at line 128 of file PerigeeLinearizedTrackState.h.

128 { return theTSOS.weight(); }

References theTSOS, and TrajectoryStateOnSurface::weight().

Friends And Related Function Documentation

◆ LinearizedTrackStateFactory

friend class LinearizedTrackStateFactory

friend

Friend class properly dealing with creation of reference-counted pointers to LinearizedTrack objects

Definition at line 38 of file PerigeeLinearizedTrackState.h.

Member Data Documentation

◆ builder

TSCPBuilderNoMaterial PerigeeLinearizedTrackState::builder

private

Definition at line 167 of file PerigeeLinearizedTrackState.h.

Referenced by computeJacobians().

◆ jacobiansAvailable

bool PerigeeLinearizedTrackState::jacobiansAvailable

mutableprivate

Definition at line 170 of file PerigeeLinearizedTrackState.h.

Referenced by computeJacobians(), constantTerm(), isValid(), momentumJacobian(), parametersFromExpansion(), positionJacobian(), predictedState(), predictedStateError(), predictedStateMomentumError(), predictedStateMomentumParameters(), predictedStateParameters(), and predictedStateWeight().

◆ theCharge

TrackCharge PerigeeLinearizedTrackState::theCharge

private

Definition at line 169 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), and computeJacobians().

◆ theConstantTerm

AlgebraicVector5 PerigeeLinearizedTrackState::theConstantTerm

mutableprivate

Definition at line 164 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeNeutralJacobians(), and constantTerm().

◆ theExpandedParams

AlgebraicVector5 PerigeeLinearizedTrackState::theExpandedParams

mutableprivate

Definition at line 165 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeNeutralJacobians(), and parametersFromExpansion().

◆ theLinPoint

GlobalPoint PerigeeLinearizedTrackState::theLinPoint

private

Definition at line 158 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeJacobians(), computeNeutralJacobians(), and linearizationPoint().

◆ theMomentumJacobian

AlgebraicMatrix53 PerigeeLinearizedTrackState::theMomentumJacobian

mutableprivate

Definition at line 162 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeNeutralJacobians(), and momentumJacobian().

◆ thePositionJacobian

AlgebraicMatrix53 PerigeeLinearizedTrackState::thePositionJacobian

mutableprivate

Definition at line 162 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeNeutralJacobians(), and positionJacobian().

◆ thePredState

TrajectoryStateClosestToPoint PerigeeLinearizedTrackState::thePredState

mutableprivate

Definition at line 163 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeJacobians(), computeNeutralJacobians(), predictedState(), predictedStateError(), predictedStateMomentumError(), predictedStateMomentumParameters(), predictedStateParameters(), and predictedStateWeight().

◆ theTrack

reco::TransientTrack PerigeeLinearizedTrackState::theTrack

private

Definition at line 159 of file PerigeeLinearizedTrackState.h.

Referenced by computeChargedJacobians(), computeJacobians(), createRefittedTrackState(), and track().

◆ theTSOS

const TrajectoryStateOnSurface PerigeeLinearizedTrackState::theTSOS

private

Definition at line 160 of file PerigeeLinearizedTrackState.h.

Referenced by computeJacobians(), isValid(), state(), stateWithNewLinearizationPoint(), and weightInMixture().

Private Attributes
TSCPBuilderNoMaterial	builder

bool	jacobiansAvailable

TrackCharge	theCharge

AlgebraicVector5	theConstantTerm

AlgebraicVector5	theExpandedParams

GlobalPoint	theLinPoint

AlgebraicMatrix53	theMomentumJacobian

AlgebraicMatrix53	thePositionJacobian

TrajectoryStateClosestToPoint	thePredState

reco::TransientTrack	theTrack

const TrajectoryStateOnSurface	theTSOS

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ RefCountedLinearizedTrackState

Constructor & Destructor Documentation

◆ PerigeeLinearizedTrackState()

Member Function Documentation

◆ checkParameters()

◆ components()

◆ compute3DImpactPoint()

◆ computeChargedJacobians()

◆ computeJacobians()

◆ computeNeutralJacobians()

◆ constantTerm()

◆ createRefittedTrackState()

◆ isValid()

◆ linearizationPoint()

◆ momentumJacobian()

◆ parametersFromExpansion()

◆ positionJacobian()

◆ predictedState()

◆ predictedStateError()

◆ predictedStateMomentumError()

◆ predictedStateMomentumParameters()

◆ predictedStateParameters()

◆ predictedStateWeight()

◆ refittedParamFromEquation()

◆ state()

◆ stateWithNewLinearizationPoint()

◆ track()

◆ weightInMixture()

Friends And Related Function Documentation

◆ LinearizedTrackStateFactory

Member Data Documentation

◆ builder

◆ jacobiansAvailable

◆ theCharge

◆ theConstantTerm

◆ theExpandedParams

◆ theLinPoint

◆ theMomentumJacobian

◆ thePositionJacobian

◆ thePredState

◆ theTrack

◆ theTSOS