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==(const 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