LinearizedTrackState.h

Go to the documentation of this file.

#ifndef LinearizedTrackState_H
#define LinearizedTrackState_H

//#include "CommonReco/CommonVertex/interface/ImpactPointMeasurement.h"

#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/TrajectoryState/interface/TrackCharge.h"
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"

#define SMATRIX_USE_CONSTEXPR
#include "Math/SMatrix.h"

#include "RecoVertex/VertexPrimitives/interface/RefittedTrackState.h"
#include <vector>

template <unsigned int N>
class LinearizedTrackState : public ReferenceCounted {
public:
  typedef ROOT::Math::SVector<double, N> AlgebraicVectorN;
  typedef ROOT::Math::SVector<double, N - 2> AlgebraicVectorM;
  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, 3, ROOT::Math::MatRepStd<double, N - 2, 3> > AlgebraicMatrixM3;
  typedef ROOT::Math::SMatrix<double, N, N, ROOT::Math::MatRepSym<double, N> > AlgebraicSymMatrixNN;
  typedef ROOT::Math::SMatrix<double, N - 2, N - 2, ROOT::Math::MatRepSym<double, N - 2> > AlgebraicSymMatrixMM;
  typedef ROOT::Math::SMatrix<double, N + 1, N + 1, ROOT::Math::MatRepSym<double, N + 1> > AlgebraicSymMatrixOO;

  typedef ReferenceCountingPointer<RefittedTrackState<N> > RefCountedRefittedTrackState;

  ~LinearizedTrackState() override {}

  virtual ReferenceCountingPointer<LinearizedTrackState<N> > stateWithNewLinearizationPoint(
      const GlobalPoint& newLP) const = 0;

  virtual const GlobalPoint& linearizationPoint() const = 0;

  virtual const AlgebraicVectorN& constantTerm() const = 0;

  virtual const AlgebraicMatrixN3& positionJacobian() const = 0;

  virtual const AlgebraicMatrixNM& momentumJacobian() const = 0;

  virtual const AlgebraicVectorN& parametersFromExpansion() const = 0;

  virtual AlgebraicVectorN predictedStateParameters() const = 0;

  virtual AlgebraicVectorM predictedStateMomentumParameters() const = 0;

  virtual AlgebraicSymMatrixNN predictedStateWeight(int& error) const = 0;

  virtual AlgebraicSymMatrixMM predictedStateMomentumError() const = 0;

  virtual AlgebraicSymMatrixNN predictedStateError() const = 0;

  virtual bool hasError() const = 0;

  virtual TrackCharge charge() const = 0;

  //   virtual ImpactPointMeasurement impactPointMeasurement() const = 0;

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

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

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

  virtual inline void checkParameters(AlgebraicVectorN& parameters) const {}

  virtual double weightInMixture() const = 0;

  virtual std::vector<ReferenceCountingPointer<LinearizedTrackState<N> > > components() const = 0;

  virtual reco::TransientTrack track() const = 0;

  virtual bool isValid() const { return true; }
};

#endif