ParticleKinematicLinearizedTrackState Class Reference

#include <ParticleKinematicLinearizedTrackState.h>

Inheritance diagram for ParticleKinematicLinearizedTrackState:

Public Types
typedef ReferenceCountingPointer< LinearizedTrackState< 6 > >	RefCountedLinearizedTrackState
Public Types inherited from LinearizedTrackState< 6 >
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 (AlgebraicVectorN &parameters) const override
std::vector< ReferenceCountingPointer< LinearizedTrackState< 6 > > >	components () const override
const AlgebraicVector6 &	constantTerm () const override
const GlobalPoint &	linearizationPoint () const override
const AlgebraicMatrix64 &	momentumJacobian () const override
const AlgebraicVector6 &	parametersFromExpansion () const override
	ParticleKinematicLinearizedTrackState ()
const AlgebraicMatrix63 &	positionJacobian () const override
AlgebraicSymMatrix66	predictedStateError () const override
AlgebraicSymMatrix44	predictedStateMomentumError () const override
AlgebraicVectorM	predictedStateMomentumParameters () const override
AlgebraicVector6	predictedStateParameters () const override
AlgebraicSymMatrix66	predictedStateWeight (int &error) const override
AlgebraicVectorN	refittedParamFromEquation (const RefCountedRefittedTrackState &theRefittedState) const override
ReferenceCountingPointer< LinearizedTrackState< 6 > >	stateWithNewLinearizationPoint (const GlobalPoint &newLP) const override
reco::TransientTrack	track () const override
double	weightInMixture () const override
Public Member Functions inherited from LinearizedTrackState< 6 >
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	isValid () const
virtual bool	operator== (LinearizedTrackState< N > &other) const=0
virtual AlgebraicVectorN	refittedParamFromEquation (const RefCountedRefittedTrackState &theRefittedState) const=0
	~LinearizedTrackState () override

Private Member Functions
void	computeChargedJacobians () const
void	computeJacobians () const
void	computeNeutralJacobians () const
	ParticleKinematicLinearizedTrackState (const GlobalPoint &linP, RefCountedKinematicParticle &prt)

Private Attributes
TransientTrackKinematicStateBuilder	builder
bool	errorAvailable
bool	impactPointAvailable
bool	jacobiansAvailable
RefCountedKinematicParticle	part
TrackCharge	theCharge
AlgebraicVector6	theConstantTerm
AlgebraicVector6	theExpandedParams
GlobalPoint	theLinPoint
AlgebraicMatrix64	theMomentumJacobian
AlgebraicMatrix63	thePositionJacobian
PerigeeKinematicState	thePredState

Friends
class	ParticleKinematicLinearizedTrackStateFactory

Detailed Description

Definition at line 13 of file ParticleKinematicLinearizedTrackState.h.

Member Typedef Documentation

RefCountedLinearizedTrackState

typedef ReferenceCountingPointer<LinearizedTrackState<6> > ParticleKinematicLinearizedTrackState::RefCountedLinearizedTrackState

Definition at line 16 of file ParticleKinematicLinearizedTrackState.h.

Constructor & Destructor Documentation

ParticleKinematicLinearizedTrackState() [1/2]

ParticleKinematicLinearizedTrackState::ParticleKinematicLinearizedTrackState ( )

inline

Definition at line 18 of file ParticleKinematicLinearizedTrackState.h.

{ jacobiansAvailable = false; }

References jacobiansAvailable.

Referenced by stateWithNewLinearizationPoint().

ParticleKinematicLinearizedTrackState() [2/2]

ParticleKinematicLinearizedTrackState::ParticleKinematicLinearizedTrackState ( const GlobalPoint &  linP, RefCountedKinematicParticle &  prt  )

inlineprivate

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

Definition at line 115 of file ParticleKinematicLinearizedTrackState.h.

      : theLinPoint(linP),
        part(prt),
        jacobiansAvailable(false),
        theCharge(prt->currentState().particleCharge()),
        impactPointAvailable(false)
 
  {}

Member Function Documentation

checkParameters()

void ParticleKinematicLinearizedTrackState::checkParameters ( AlgebraicVectorN &  parameters ) const

inlineoverride

Definition at line 172 of file ParticleKinematicLinearizedTrackState.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< 6 > > > ParticleKinematicLinearizedTrackState::components ( ) const

overridevirtual

Implements LinearizedTrackState< 6 >.

Definition at line 140 of file ParticleKinematicLinearizedTrackState.cc.

          {
  std::vector<ReferenceCountingPointer<LinearizedTrackState<6> > > res;
  res.reserve(1);
  res.push_back(RefCountedLinearizedTrackState(const_cast<ParticleKinematicLinearizedTrackState*>(this)));
  return res;
}

computeChargedJacobians()

void ParticleKinematicLinearizedTrackState::computeChargedJacobians ( ) const

private

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

Definition at line 179 of file ParticleKinematicLinearizedTrackState.cc.

                                                                          {
  GlobalPoint paramPt(theLinPoint);
 
  double field = part->currentState().magneticField()->inInverseGeV(thePredState.theState().globalPosition()).z();
  double signTC = -part->currentState().particleCharge();
 
  double thetaAtEP = thePredState.theState().globalMomentum().theta();
  double phiAtEP = thePredState.theState().globalMomentum().phi();
  double ptAtEP = thePredState.theState().globalMomentum().perp();
  double transverseCurvatureAtEP = field / ptAtEP * signTC;
 
  double x_v = thePredState.theState().globalPosition().x();
  double y_v = thePredState.theState().globalPosition().y();
  double z_v = thePredState.theState().globalPosition().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;
 
  theExpandedParams[5] = part->currentState().mass();
 
  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;
 
  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;
 
  //debug code - to be removed later
  //   cout<<"parameters for momentum jacobian: X "<<X<<endl;
  //   cout<<"parameters for momentum jacobian: Y "<<Y<<endl;
  //   cout<<"parameters for momentum jacobian: SS "<<SS<<endl;
  //   cout<<"parameters for momentum jacobian: PhiAtEP "<<phiAtEP<<endl;
  //   cout<<"parameters for momentum jacobian: curv "<<transverseCurvatureAtEP<<endl;
  //   cout<<"sin phi Atep "<<sin(phiAtEP)<<endl;
  //   cout<<"cos phi at EP "<<cos(phiAtEP)<<endl;
  //   cout<<"upper  part is "<<X*cos(phiAtEP) + Y*sin(phiAtEP) <<endl;
  //   cout<<"lower part is"<<SS*transverseCurvatureAtEP*transverseCurvatureAtEP<<endl;
 
  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;
 
  theMomentumJacobian(5, 3) = 1;
 
  // And finally the residuals:
  AlgebraicVector3 expansionPoint;
  expansionPoint[0] = thePredState.theState().globalPosition().x();
  expansionPoint[1] = thePredState.theState().globalPosition().y();
  expansionPoint[2] = thePredState.theState().globalPosition().z();
  AlgebraicVector4 momentumAtExpansionPoint;
  momentumAtExpansionPoint[0] = transverseCurvatureAtEP;  // Transverse Curv
  momentumAtExpansionPoint[1] = thetaAtEP;
  momentumAtExpansionPoint[2] = phiAtEP;
  momentumAtExpansionPoint[3] = theExpandedParams[5];
 
  theConstantTerm = AlgebraicVector6(theExpandedParams - thePositionJacobian * expansionPoint -
                                     theMomentumJacobian * momentumAtExpansionPoint);
}

References funct::abs(), funct::cos(), KinematicState::globalMomentum(), KinematicState::globalPosition(), M_PI, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), S(), funct::sin(), mathSSE::sqrt(), funct::tan(), theConstantTerm, theExpandedParams, theLinPoint, theMomentumJacobian, thePositionJacobian, thePredState, PerigeeKinematicState::theState(), PV3DBase< T, PVType, FrameType >::theta(), 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 ParticleKinematicLinearizedTrackState::computeJacobians ( ) const

private

Method calculating the track parameters and the Jacobians.

Definition at line 79 of file ParticleKinematicLinearizedTrackState.cc.

                                                                   {
  GlobalPoint paramPt(theLinPoint);
  thePredState = builder(part->currentState(), paramPt);
  //  bool valid = thePredState.isValid();
  //  if (!valid) std::cout <<"Help!!!!!!!!! State is invalid\n";
  //  if (!valid) return;
  if (std::abs(theCharge) < 1e-5) {
    //neutral track
    computeNeutralJacobians();
  } else {
    //charged track
    computeChargedJacobians();
  }
  jacobiansAvailable = true;
}

References funct::abs(), builder, computeChargedJacobians(), computeNeutralJacobians(), MillePedeFileConverter_cfg::e, jacobiansAvailable, theCharge, theLinPoint, and thePredState.

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

computeNeutralJacobians()

void ParticleKinematicLinearizedTrackState::computeNeutralJacobians ( ) const

private

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

Definition at line 278 of file ParticleKinematicLinearizedTrackState.cc.

                                                                          {
  GlobalPoint paramPt(theLinPoint);
  double thetaAtEP = thePredState.theState().globalMomentum().theta();
  double phiAtEP = thePredState.theState().globalMomentum().phi();
  double ptAtEP = thePredState.theState().globalMomentum().perp();
 
  double x_v = thePredState.theState().globalPosition().x();
  double y_v = thePredState.theState().globalPosition().y();
  double z_v = thePredState.theState().globalPosition().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);
  theExpandedParams[5] = part->currentState().mass();
 
  // 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);
  theMomentumJacobian(5, 3) = 1;
 
  // And finally the residuals:
  AlgebraicVector3 expansionPoint;
  expansionPoint[0] = thePredState.theState().globalPosition().x();
  expansionPoint[1] = thePredState.theState().globalPosition().y();
  expansionPoint[2] = thePredState.theState().globalPosition().z();
  AlgebraicVector4 momentumAtExpansionPoint;
  momentumAtExpansionPoint[0] = 1 / ptAtEP;
  momentumAtExpansionPoint[1] = thetaAtEP;
  momentumAtExpansionPoint[2] = phiAtEP;
  momentumAtExpansionPoint[3] = theExpandedParams[5];
 
  theConstantTerm = AlgebraicVector6(theExpandedParams - thePositionJacobian * expansionPoint -
                                     theMomentumJacobian * momentumAtExpansionPoint);
}

References funct::cos(), KinematicState::globalMomentum(), KinematicState::globalPosition(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), funct::sin(), funct::tan(), theConstantTerm, theExpandedParams, theLinPoint, theMomentumJacobian, thePositionJacobian, thePredState, PerigeeKinematicState::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 AlgebraicVector6 & ParticleKinematicLinearizedTrackState::constantTerm ( ) const

overridevirtual

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

Implements LinearizedTrackState< 6 >.

Definition at line 5 of file ParticleKinematicLinearizedTrackState.cc.

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

References computeJacobians(), jacobiansAvailable, and theConstantTerm.

Referenced by refittedParamFromEquation().

linearizationPoint()

const GlobalPoint& ParticleKinematicLinearizedTrackState::linearizationPoint ( ) const

inlineoverridevirtual

The point at which the track state has been linearized

Implements LinearizedTrackState< 6 >.

Definition at line 30 of file ParticleKinematicLinearizedTrackState.h.

{ return theLinPoint; }

References theLinPoint.

momentumJacobian()

const AlgebraicMatrix64 & ParticleKinematicLinearizedTrackState::momentumJacobian ( ) const

overridevirtual

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

Implements LinearizedTrackState< 6 >.

Definition at line 17 of file ParticleKinematicLinearizedTrackState.cc.

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

References computeJacobians(), jacobiansAvailable, and theMomentumJacobian.

Referenced by refittedParamFromEquation().

parametersFromExpansion()

const AlgebraicVector6 & ParticleKinematicLinearizedTrackState::parametersFromExpansion ( ) const

overridevirtual

Method returning the parameters of the Taylor expansion

Implements LinearizedTrackState< 6 >.

Definition at line 23 of file ParticleKinematicLinearizedTrackState.cc.

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

References computeJacobians(), jacobiansAvailable, and theExpandedParams.

positionJacobian()

const AlgebraicMatrix63 & ParticleKinematicLinearizedTrackState::positionJacobian ( ) const

overridevirtual

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

Implements LinearizedTrackState< 6 >.

Definition at line 11 of file ParticleKinematicLinearizedTrackState.cc.

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

References computeJacobians(), jacobiansAvailable, and thePositionJacobian.

Referenced by refittedParamFromEquation().

predictedStateError()

AlgebraicSymMatrix66 ParticleKinematicLinearizedTrackState::predictedStateError ( ) const

overridevirtual

Full predicted error matrix

Implements LinearizedTrackState< 6 >.

Definition at line 45 of file ParticleKinematicLinearizedTrackState.cc.

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

References computeJacobians(), ExtendedPerigeeTrajectoryError::covarianceMatrix(), jacobiansAvailable, PerigeeKinematicState::perigeeError(), and thePredState.

predictedStateMomentumError()

AlgebraicSymMatrix44 ParticleKinematicLinearizedTrackState::predictedStateMomentumError ( ) const

overridevirtual

4x4 error matrix ofe xtended perigee mometum components

Implements LinearizedTrackState< 6 >.

Definition at line 124 of file ParticleKinematicLinearizedTrackState.cc.

                                                                                              {
  int error;
  if (!jacobiansAvailable)
    computeJacobians();
  AlgebraicSymMatrix44 res;
  AlgebraicSymMatrix33 m3 = thePredState.perigeeError().weightMatrix(error).Sub<AlgebraicSymMatrix33>(0, 0);
  res.Place_at(m3, 0, 0);
  res(3, 0) = thePredState.perigeeError().weightMatrix(error)(5, 5);
  res(3, 1) = thePredState.perigeeError().weightMatrix(error)(5, 0);
  res(3, 2) = thePredState.perigeeError().weightMatrix(error)(5, 1);
  res(3, 3) = thePredState.perigeeError().weightMatrix(error)(5, 2);
  return res;
}

References computeJacobians(), relativeConstraints::error, jacobiansAvailable, PerigeeKinematicState::perigeeError(), thePredState, and ExtendedPerigeeTrajectoryError::weightMatrix().

predictedStateMomentumParameters()

AlgebraicVector4 ParticleKinematicLinearizedTrackState::predictedStateMomentumParameters ( ) const

overridevirtual

returns predicted 4-momentum in extended perigee parametrization

Implements LinearizedTrackState< 6 >.

Definition at line 113 of file ParticleKinematicLinearizedTrackState.cc.

                                                                                               {
  if (!jacobiansAvailable)
    computeJacobians();
  AlgebraicVector4 res;
  res[0] = thePredState.perigeeParameters().vector()(0);
  res[1] = thePredState.perigeeParameters().vector()(1);
  res[2] = thePredState.perigeeParameters().vector()(2);
  res[3] = thePredState.perigeeParameters().vector()(5);
  return res;
}

References computeJacobians(), jacobiansAvailable, PerigeeKinematicState::perigeeParameters(), thePredState, and ExtendedPerigeeTrajectoryParameters::vector().

Referenced by refittedParamFromEquation().

predictedStateParameters()

AlgebraicVector6 ParticleKinematicLinearizedTrackState::predictedStateParameters ( ) const

overridevirtual

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). extended perigee predicted parameters

Implements LinearizedTrackState< 6 >.

Definition at line 29 of file ParticleKinematicLinearizedTrackState.cc.

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

References computeJacobians(), jacobiansAvailable, PerigeeKinematicState::perigeeParameters(), thePredState, and ExtendedPerigeeTrajectoryParameters::vector().

predictedStateWeight()

AlgebraicSymMatrix66 ParticleKinematicLinearizedTrackState::predictedStateWeight ( int &  error ) const

overridevirtual

Full predicted weight matrix

Implements LinearizedTrackState< 6 >.

Definition at line 35 of file ParticleKinematicLinearizedTrackState.cc.

                                                                                                 {
  if (!jacobiansAvailable)
    computeJacobians();
  int i = 0;
  AlgebraicSymMatrix66 z = thePredState.perigeeError().weightMatrix(i);
  error = i;
  return z;
  //   return thePredState.perigeeError().weightMatrix(error);
}

References computeJacobians(), relativeConstraints::error, mps_fire::i, jacobiansAvailable, PerigeeKinematicState::perigeeError(), thePredState, ExtendedPerigeeTrajectoryError::weightMatrix(), and z.

refittedParamFromEquation()

AlgebraicVector6 ParticleKinematicLinearizedTrackState::refittedParamFromEquation ( const RefCountedRefittedTrackState &  theRefittedState ) const

override

    Method returning the impact point measurement
   ‍/

ImpactPointMeasurement impactPointMeasurement() const;

TrackCharge charge() const override;

RefCountedKinematicParticle particle() const;

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

bool hasError() const override;

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

/** Method returning the parameters of the Taylor expansion evaluated with the refitted state.

Definition at line 148 of file ParticleKinematicLinearizedTrackState.cc.

                                                                {
  AlgebraicVectorM momentum = theRefittedState->momentumVector();
  if ((momentum(2) * predictedStateMomentumParameters()(2) < 0) && (std::fabs(momentum(2)) > M_PI / 2)) {
    if (predictedStateMomentumParameters()(2) < 0.)
      momentum(2) -= 2 * M_PI;
    if (predictedStateMomentumParameters()(2) > 0.)
      momentum(2) += 2 * M_PI;
  }
 
  AlgebraicVector3 vertexPosition;
  vertexPosition(0) = theRefittedState->position().x();
  vertexPosition(1) = theRefittedState->position().y();
  vertexPosition(2) = theRefittedState->position().z();
  AlgebraicVector6 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(), positionJacobian(), and predictedStateMomentumParameters().

stateWithNewLinearizationPoint()

ReferenceCountingPointer< LinearizedTrackState< 6 > > ParticleKinematicLinearizedTrackState::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< 6 >.

Definition at line 96 of file ParticleKinematicLinearizedTrackState.cc.

                                                                                                    {
  RefCountedKinematicParticle pr = part;
  return new ParticleKinematicLinearizedTrackState(newLP, pr);
}

References part, and ParticleKinematicLinearizedTrackState().

track()

reco::TransientTrack ParticleKinematicLinearizedTrackState::track ( ) const

overridevirtual

Implements LinearizedTrackState< 6 >.

Definition at line 336 of file ParticleKinematicLinearizedTrackState.cc.

                                                                      {
  throw VertexException(" ParticleKinematicLinearizedTrackState:: no TransientTrack to return");
}

weightInMixture()

double ParticleKinematicLinearizedTrackState::weightInMixture ( ) const

overridevirtual

Implements LinearizedTrackState< 6 >.

Definition at line 138 of file ParticleKinematicLinearizedTrackState.cc.

{ return 1.; }

Friends And Related Function Documentation

ParticleKinematicLinearizedTrackStateFactory

friend class ParticleKinematicLinearizedTrackStateFactory

friend

Definition at line 15 of file ParticleKinematicLinearizedTrackState.h.

Member Data Documentation

builder

TransientTrackKinematicStateBuilder ParticleKinematicLinearizedTrackState::builder

private

Definition at line 140 of file ParticleKinematicLinearizedTrackState.h.

Referenced by computeJacobians().

errorAvailable

bool ParticleKinematicLinearizedTrackState::errorAvailable

mutableprivate

Definition at line 142 of file ParticleKinematicLinearizedTrackState.h.

impactPointAvailable

bool ParticleKinematicLinearizedTrackState::impactPointAvailable

mutableprivate

Definition at line 152 of file ParticleKinematicLinearizedTrackState.h.

jacobiansAvailable

bool ParticleKinematicLinearizedTrackState::jacobiansAvailable

mutableprivate

Definition at line 143 of file ParticleKinematicLinearizedTrackState.h.

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

part

RefCountedKinematicParticle ParticleKinematicLinearizedTrackState::part

private

Definition at line 139 of file ParticleKinematicLinearizedTrackState.h.

Referenced by stateWithNewLinearizationPoint().

theCharge

TrackCharge ParticleKinematicLinearizedTrackState::theCharge

private

Definition at line 150 of file ParticleKinematicLinearizedTrackState.h.

Referenced by computeJacobians().

theConstantTerm

AlgebraicVector6 ParticleKinematicLinearizedTrackState::theConstantTerm

mutableprivate

Definition at line 147 of file ParticleKinematicLinearizedTrackState.h.

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

theExpandedParams

AlgebraicVector6 ParticleKinematicLinearizedTrackState::theExpandedParams

mutableprivate

Definition at line 148 of file ParticleKinematicLinearizedTrackState.h.

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

theLinPoint

GlobalPoint ParticleKinematicLinearizedTrackState::theLinPoint

private

Definition at line 138 of file ParticleKinematicLinearizedTrackState.h.

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

theMomentumJacobian

AlgebraicMatrix64 ParticleKinematicLinearizedTrackState::theMomentumJacobian

mutableprivate

Definition at line 145 of file ParticleKinematicLinearizedTrackState.h.

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

thePositionJacobian

AlgebraicMatrix63 ParticleKinematicLinearizedTrackState::thePositionJacobian

mutableprivate

Definition at line 144 of file ParticleKinematicLinearizedTrackState.h.

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

thePredState

PerigeeKinematicState ParticleKinematicLinearizedTrackState::thePredState

mutableprivate

Definition at line 146 of file ParticleKinematicLinearizedTrackState.h.

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