LinearizedTrackState.h

Go to the documentation of this file.

#ifndef LinearizedTrackState_H
#define LinearizedTrackState_H

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

//#include "TrackingTools/TrajectoryState/interface/ProxyBase.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "TrackingTools/TrajectoryState/interface/CopyUsingClone.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 ==(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