ParticleKinematicLinearizedTrackState.h

Go to the documentation of this file.

#ifndef ParticleKinematicLinearizedTrackState_H
#define ParticleKinematicLinearizedTrackState_H
 
#include "RecoVertex/VertexPrimitives/interface/LinearizedTrackState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/RefCountedKinematicParticle.h"
#include "RecoVertex/VertexPrimitives/interface/RefittedTrackState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/PerigeeKinematicState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/TransientTrackKinematicStateBuilder.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "DataFormats/CLHEP/interface/Migration.h"
#include "RecoVertex/KinematicFitPrimitives/interface/Matrices.h"
 
class ParticleKinematicLinearizedTrackState : public LinearizedTrackState<6> {
public:
  friend class ParticleKinematicLinearizedTrackStateFactory;
  typedef ReferenceCountingPointer<LinearizedTrackState<6> > RefCountedLinearizedTrackState;
 
  ParticleKinematicLinearizedTrackState() { jacobiansAvailable = false; }
 
  ReferenceCountingPointer<LinearizedTrackState<6> > stateWithNewLinearizationPoint(
      const GlobalPoint& newLP) const override;
 
  const GlobalPoint& linearizationPoint() const override { return theLinPoint; }
 
  const AlgebraicVector6& constantTerm() const override;
 
  const AlgebraicMatrix63& positionJacobian() const override;
 
  const AlgebraicMatrix64& momentumJacobian() const override;
 
  const AlgebraicVector6& parametersFromExpansion() const override;
 
  AlgebraicVector6 predictedStateParameters() const override;
 
  AlgebraicVectorM predictedStateMomentumParameters() const override;
 
  AlgebraicSymMatrix44 predictedStateMomentumError() const override;
 
  AlgebraicSymMatrix66 predictedStateWeight(int& error) const override;
 
  AlgebraicSymMatrix66 predictedStateError() 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;
 
  AlgebraicVectorN refittedParamFromEquation(const RefCountedRefittedTrackState& theRefittedState) const override;
 
  inline void checkParameters(AlgebraicVectorN& parameters) const override;
 
  double weightInMixture() const override;
 
  std::vector<ReferenceCountingPointer<LinearizedTrackState<6> > > components() const override;
 
  reco::TransientTrack track() const override;
 
private:
  ParticleKinematicLinearizedTrackState(const GlobalPoint& linP, RefCountedKinematicParticle& prt)
      : theLinPoint(linP),
        part(prt),
        jacobiansAvailable(false),
        theCharge(prt->currentState().particleCharge()),
        impactPointAvailable(false)
 
  {}
 
  void computeJacobians() const;
  void computeChargedJacobians() const;
 
  void computeNeutralJacobians() const;
 
  GlobalPoint theLinPoint;
  RefCountedKinematicParticle part;
  TransientTrackKinematicStateBuilder builder;
 
  mutable bool errorAvailable;
  mutable bool jacobiansAvailable;
  mutable AlgebraicMatrix63 thePositionJacobian;
  mutable AlgebraicMatrix64 theMomentumJacobian;
  mutable PerigeeKinematicState thePredState;
  mutable AlgebraicVector6 theConstantTerm;
  mutable AlgebraicVector6 theExpandedParams;
 
  TrackCharge theCharge;
  //   ImpactPointMeasurementExtractor theIPMExtractor;
  mutable bool impactPointAvailable;
};
#endif