#include <PerigeeMultiLTS.h>

Inheritance diagram for PerigeeMultiLTS:

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
TrackCharge	charge () const override

void	checkParameters (AlgebraicVector5 &parameters) const override

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

const AlgebraicVectorN &	constantTerm () const override

RefCountedRefittedTrackState	createRefittedTrackState (const GlobalPoint &vertexPosition, const AlgebraicVectorM &vectorParameters, const AlgebraicSymMatrixOO &covarianceMatrix) const override

bool	hasError () const override

const GlobalPoint &	linearizationPoint () const override

const AlgebraicMatrixNM &	momentumJacobian () const override

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

const AlgebraicVectorN &	parametersFromExpansion () const override

const AlgebraicMatrixN3 &	positionJacobian () const override

const TrajectoryStateClosestToPoint &	predictedState () const

AlgebraicSymMatrixNN	predictedStateError () const override

AlgebraicSymMatrixMM	predictedStateMomentumError () const override

AlgebraicVectorM	predictedStateMomentumParameters () const override

AlgebraicVectorN	predictedStateParameters () const override

AlgebraicSymMatrixNN	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 void	checkParameters (AlgebraicVectorN &parameters) const

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

virtual bool	isValid () const

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

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

	~LinearizedTrackState () override

Private Member Functions
	PerigeeMultiLTS (const GlobalPoint &linP, const reco::TransientTrack &track, const TrajectoryStateOnSurface &tsos)

void	prepareCollapsedState () const

Private Attributes
bool	collapsedStateAvailable

RefCountedLinearizedTrackState	collapsedStateLT

std::vector< RefCountedLinearizedTrackState >	ltComp

TrackCharge	theCharge

GlobalPoint	theLinPoint

LinearizedTrackStateFactory	theLTSfactory

reco::TransientTrack	theTrack

const TrajectoryStateOnSurface	theTSOS

Friends
class	MultiPerigeeLTSFactory

Detailed Description

Multi-state version of PerigeeLinearizedTrackState

Definition at line 14 of file PerigeeMultiLTS.h.

Member Typedef Documentation

◆ RefCountedLinearizedTrackState

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

Definition at line 16 of file PerigeeMultiLTS.h.

Constructor & Destructor Documentation

◆ PerigeeMultiLTS()

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

private

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

Parameters

linP	The linearization point
track	The RecTrack
tsos	The original (multi-state) TrajectoryStateOnSurface

Definition at line 9 of file PerigeeMultiLTS.cc.

     : theLinPoint(linP), theTrack(track), theTSOS(tsos), theCharge(theTrack.charge()), collapsedStateAvailable(false) {
   //Extract the TSOS components:
  
   vector<TrajectoryStateOnSurface> tsosComp = theTSOS.components();
   //cout << "PerigeeMultiLTS components: "<<tsosComp.size()<<endl;
   ltComp.reserve(tsosComp.size());
   for (vector<TrajectoryStateOnSurface>::iterator it = tsosComp.begin(); it != tsosComp.end(); it++) {
     // cout <<(*it).globalPosition()<<endl;
     ltComp.push_back(theLTSfactory.linearizedTrackState(theLinPoint, theTrack, *it));
   }
 }

References TrajectoryStateOnSurface::components(), LinearizedTrackStateFactory::linearizedTrackState(), ltComp, theLinPoint, theLTSfactory, theTrack, and theTSOS.

Referenced by stateWithNewLinearizationPoint().

Member Function Documentation

◆ charge()

TrackCharge PerigeeMultiLTS::charge ( void ) const

inlineoverridevirtual

Implements LinearizedTrackState< 5 >.

Definition at line 95 of file PerigeeMultiLTS.h.

95 { return theCharge; }

References theCharge.

Referenced by createRefittedTrackState().

◆ checkParameters()

void PerigeeMultiLTS::checkParameters ( AlgebraicVector5 & parameters ) const

inlineoverride

Definition at line 150 of file PerigeeMultiLTS.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<ReferenceCountingPointer<LinearizedTrackState<5> > > PerigeeMultiLTS::components ( ) const

inlineoverridevirtual

Vector of individual components in the mixture.

Implements LinearizedTrackState< 5 >.

Definition at line 122 of file PerigeeMultiLTS.h.

122 { return ltComp; }

References ltComp.

◆ constantTerm()

const AlgebraicVector5 & PerigeeMultiLTS::constantTerm ( ) const

overridevirtual

Method returning the constant term of the Taylor expansion of the measurement equation for the collapsed track state

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

Implements LinearizedTrackState< 5 >.

Definition at line 34 of file PerigeeMultiLTS.cc.

                                                             {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->constantTerm();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

Referenced by refittedParamFromEquation().

◆ createRefittedTrackState()

PerigeeMultiLTS::RefCountedRefittedTrackState PerigeeMultiLTS::createRefittedTrackState	(	const GlobalPoint &	vertexPosition,
		const AlgebraicVectorM &	vectorParameters,
		const AlgebraicSymMatrixOO &	covarianceMatrix
	)		const

override

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

Definition at line 126 of file PerigeeMultiLTS.cc.

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

References charge(), reco::TransientTrack::field(), theTrack, and PerigeeConversions::trajectoryStateClosestToPoint().

◆ hasError()

bool PerigeeMultiLTS::hasError ( void ) const

overridevirtual

Implements LinearizedTrackState< 5 >.

Definition at line 77 of file PerigeeMultiLTS.cc.

                                      {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->hasError();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

◆ linearizationPoint()

const GlobalPoint& PerigeeMultiLTS::linearizationPoint ( ) const

inlineoverridevirtual

The point at which the track state has been linearized

Implements LinearizedTrackState< 5 >.

Definition at line 33 of file PerigeeMultiLTS.h.

33 { return theLinPoint; }

References theLinPoint.

◆ momentumJacobian()

const AlgebraicMatrix53 & PerigeeMultiLTS::momentumJacobian ( ) const

overridevirtual

Method returning the Momentum Jacobian from the Taylor expansion (Matrix B) for the collapsed track state

Method returning the Momentum Jacobian (Matrix B)

Implements LinearizedTrackState< 5 >.

Definition at line 52 of file PerigeeMultiLTS.cc.

                                                                  {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->momentumJacobian();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

Referenced by refittedParamFromEquation().

◆ operator==()

bool PerigeeMultiLTS::operator== ( LinearizedTrackState< 5 > & other ) const

override

Definition at line 113 of file PerigeeMultiLTS.cc.

                                                                      {
   const PerigeeMultiLTS* otherP = dynamic_cast<const PerigeeMultiLTS*>(&other);
   if (otherP == nullptr) {
     throw VertexException("PerigeeMultiLTS: don't know how to compare myself to non-perigee track state");
   }
   return (otherP->track() == theTrack);
 }

References trackingPlots::other, theTrack, and track().

◆ parametersFromExpansion()

const AlgebraicVector5 & PerigeeMultiLTS::parametersFromExpansion ( ) const

overridevirtual

Method returning the parameters of the Taylor expansion for the collapsed track state

Method returning the parameters of the Taylor expansion

Implements LinearizedTrackState< 5 >.

Definition at line 60 of file PerigeeMultiLTS.cc.

                                                                        {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->parametersFromExpansion();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

◆ positionJacobian()

const AlgebraicMatrix53 & PerigeeMultiLTS::positionJacobian ( ) const

overridevirtual

Method returning the Position Jacobian from the Taylor expansion (Matrix A) for the collapsed track state

Method returning the Position Jacobian (Matrix A)

Implements LinearizedTrackState< 5 >.

Definition at line 43 of file PerigeeMultiLTS.cc.

                                                                  {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->positionJacobian();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

Referenced by refittedParamFromEquation().

◆ predictedState()

const TrajectoryStateClosestToPoint & PerigeeMultiLTS::predictedState ( ) const

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), for the collapsed track state

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

Definition at line 70 of file PerigeeMultiLTS.cc.

                                                                            {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   const PerigeeLinearizedTrackState* otherP = dynamic_cast<const PerigeeLinearizedTrackState*>(collapsedStateLT.get());
   return otherP->predictedState();
 }

References collapsedStateAvailable, collapsedStateLT, PerigeeLinearizedTrackState::predictedState(), and prepareCollapsedState().

◆ predictedStateError()

AlgebraicSymMatrix55 PerigeeMultiLTS::predictedStateError ( ) const

overridevirtual

Method returning the covariance matrix of the track state at the transverse impact point, for the collapsed track state

Implements LinearizedTrackState< 5 >.

Definition at line 101 of file PerigeeMultiLTS.cc.

                                                                 {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->predictedStateError();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

◆ predictedStateMomentumError()

AlgebraicSymMatrix33 PerigeeMultiLTS::predictedStateMomentumError ( ) const

overridevirtual

Method returning the momentum covariance matrix of the track state at the transverse impact point, for the collapsed track state

Implements LinearizedTrackState< 5 >.

Definition at line 107 of file PerigeeMultiLTS.cc.

                                                                         {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->predictedStateMomentumError();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

◆ predictedStateMomentumParameters()

AlgebraicVector3 PerigeeMultiLTS::predictedStateMomentumParameters ( ) const

overridevirtual

Method returning the momentum part of the parameters of the track state at the linearization point, for the collapsed track state

Implements LinearizedTrackState< 5 >.

Definition at line 89 of file PerigeeMultiLTS.cc.

                                                                          {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->predictedStateMomentumParameters();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

◆ predictedStateParameters()

AlgebraicVector5 PerigeeMultiLTS::predictedStateParameters ( ) const

overridevirtual

Method returning the parameters of the track state at the transverse impact point, for the collapsed track state

Implements LinearizedTrackState< 5 >.

Definition at line 83 of file PerigeeMultiLTS.cc.

                                                                  {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->predictedStateParameters();
 }

References collapsedStateAvailable, collapsedStateLT, and prepareCollapsedState().

◆ predictedStateWeight()

AlgebraicSymMatrix55 PerigeeMultiLTS::predictedStateWeight ( int & error ) const

overridevirtual

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

Implements LinearizedTrackState< 5 >.

Definition at line 95 of file PerigeeMultiLTS.cc.

                                                                            {
   if (!collapsedStateAvailable)
     prepareCollapsedState();
   return collapsedStateLT->predictedStateWeight(error);
 }

References collapsedStateAvailable, collapsedStateLT, relativeConstraints::error, and prepareCollapsedState().

◆ prepareCollapsedState()

void PerigeeMultiLTS::prepareCollapsedState ( ) const

private

Linearize the collapsed track state.

Definition at line 24 of file PerigeeMultiLTS.cc.

                                                   {
   if (!collapsedStateAvailable) {
     collapsedStateLT = theLTSfactory.linearizedTrackState(theLinPoint, theTrack, theTSOS);
   }
   collapsedStateAvailable = true;
 }

References collapsedStateAvailable, collapsedStateLT, LinearizedTrackStateFactory::linearizedTrackState(), theLinPoint, theLTSfactory, theTrack, and theTSOS.

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

◆ refittedParamFromEquation()

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

override

Definition at line 135 of file PerigeeMultiLTS.cc.

                                                                                                                       {
   AlgebraicVector3 vertexPosition;
   vertexPosition(0) = theRefittedState->position().x();
   vertexPosition(1) = theRefittedState->position().y();
   vertexPosition(2) = theRefittedState->position().z();
   AlgebraicVectorN param =
       constantTerm() + positionJacobian() * vertexPosition + momentumJacobian() * theRefittedState->momentumVector();
   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(), and positionJacobian().

◆ state()

const TrajectoryStateOnSurface PerigeeMultiLTS::state ( ) const

inline

Returns theoriginal (multi-state) TrajectoryStateOnSurface whith which this instance has been created with.

Definition at line 41 of file PerigeeMultiLTS.h.

41 { return theTSOS; }

References theTSOS.

◆ stateWithNewLinearizationPoint()

PerigeeMultiLTS::RefCountedLinearizedTrackState PerigeeMultiLTS::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 121 of file PerigeeMultiLTS.cc.

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

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

◆ track()

reco::TransientTrack PerigeeMultiLTS::track ( ) const

inlineoverridevirtual

Implements LinearizedTrackState< 5 >.

Definition at line 35 of file PerigeeMultiLTS.h.

35 { return theTrack; }

References theTrack.

Referenced by operator==(), and stateWithNewLinearizationPoint().

◆ weightInMixture()

double PerigeeMultiLTS::weightInMixture ( ) const

inlineoverridevirtual

The weight of this state. It will be the sum of the weights of the individual components in the mixture.

Implements LinearizedTrackState< 5 >.

Definition at line 116 of file PerigeeMultiLTS.h.

116 { return theTSOS.weight(); }

References theTSOS, and TrajectoryStateOnSurface::weight().

Friends And Related Function Documentation

◆ MultiPerigeeLTSFactory

friend class MultiPerigeeLTSFactory

friend

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

Definition at line 21 of file PerigeeMultiLTS.h.

Member Data Documentation

◆ collapsedStateAvailable

bool PerigeeMultiLTS::collapsedStateAvailable

mutableprivate

Definition at line 145 of file PerigeeMultiLTS.h.

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

◆ collapsedStateLT

RefCountedLinearizedTrackState PerigeeMultiLTS::collapsedStateLT

mutableprivate

Definition at line 142 of file PerigeeMultiLTS.h.

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

◆ ltComp

std::vector<RefCountedLinearizedTrackState> PerigeeMultiLTS::ltComp

private

Definition at line 141 of file PerigeeMultiLTS.h.

Referenced by components(), and PerigeeMultiLTS().

◆ theCharge

TrackCharge PerigeeMultiLTS::theCharge

private

Definition at line 144 of file PerigeeMultiLTS.h.

Referenced by charge().

◆ theLinPoint

GlobalPoint PerigeeMultiLTS::theLinPoint

private

Definition at line 137 of file PerigeeMultiLTS.h.

Referenced by linearizationPoint(), PerigeeMultiLTS(), and prepareCollapsedState().

◆ theLTSfactory

LinearizedTrackStateFactory PerigeeMultiLTS::theLTSfactory

private

Definition at line 143 of file PerigeeMultiLTS.h.

Referenced by PerigeeMultiLTS(), and prepareCollapsedState().

◆ theTrack

reco::TransientTrack PerigeeMultiLTS::theTrack

private

Definition at line 138 of file PerigeeMultiLTS.h.

Referenced by createRefittedTrackState(), operator==(), PerigeeMultiLTS(), prepareCollapsedState(), and track().

◆ theTSOS

const TrajectoryStateOnSurface PerigeeMultiLTS::theTSOS

private

Definition at line 140 of file PerigeeMultiLTS.h.

Referenced by PerigeeMultiLTS(), prepareCollapsedState(), state(), stateWithNewLinearizationPoint(), and weightInMixture().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ RefCountedLinearizedTrackState

Constructor & Destructor Documentation

◆ PerigeeMultiLTS()

Member Function Documentation

◆ charge()

◆ checkParameters()

◆ components()

◆ constantTerm()

◆ createRefittedTrackState()

◆ hasError()

◆ linearizationPoint()

◆ momentumJacobian()

◆ operator==()

◆ parametersFromExpansion()

◆ positionJacobian()

◆ predictedState()

◆ predictedStateError()

◆ predictedStateMomentumError()

◆ predictedStateMomentumParameters()

◆ predictedStateParameters()

◆ predictedStateWeight()

◆ prepareCollapsedState()

◆ refittedParamFromEquation()

◆ state()

◆ stateWithNewLinearizationPoint()

◆ track()

◆ weightInMixture()

Friends And Related Function Documentation

◆ MultiPerigeeLTSFactory

Member Data Documentation

◆ collapsedStateAvailable

◆ collapsedStateLT

◆ ltComp

◆ theCharge

◆ theLinPoint

◆ theLTSfactory

◆ theTrack

◆ theTSOS