KalmanVertexTrackUpdator.h

Go to the documentation of this file.

#ifndef KalmanVertexTrackUpdator_H
#define KalmanVertexTrackUpdator_H

#include "RecoVertex/VertexPrimitives/interface/VertexTrackUpdator.h"
#include "RecoVertex/VertexPrimitives/interface/CachingVertex.h"
#include "RecoVertex/VertexPrimitives/interface/VertexTrack.h"
#include "RecoVertex/VertexTools/interface/VertexTrackFactory.h"
#include "RecoVertex/KalmanVertexFit/interface/KVFHelper.h"
#include "RecoVertex/KalmanVertexFit/interface/KalmanVertexUpdator.h"

template <unsigned int N>
class KalmanVertexTrackUpdator : public VertexTrackUpdator<N> {

public:

  typedef typename CachingVertex<N>::RefCountedVertexTrack RefCountedVertexTrack;
  typedef typename VertexTrack<N>::RefCountedLinearizedTrackState RefCountedLinearizedTrackState;
  typedef typename VertexTrack<N>::RefCountedRefittedTrackState RefCountedRefittedTrackState;


  KalmanVertexTrackUpdator(){}

  virtual ~KalmanVertexTrackUpdator(){}

  RefCountedVertexTrack update(const CachingVertex<N> & vertex,
                               RefCountedVertexTrack track) const;


  KalmanVertexTrackUpdator<N> * clone() const
  {
    return new KalmanVertexTrackUpdator(*this);
  }

  typedef ROOT::Math::SMatrix<double,3,N-2,ROOT::Math::MatRepStd<double,3,N-2> > AlgebraicMatrix3M;
  typedef ROOT::Math::SMatrix<double,N+1,N+1,ROOT::Math::MatRepSym<double,N+1> > AlgebraicSymMatrixOO;
  typedef std::pair< RefCountedRefittedTrackState, AlgebraicSymMatrixOO > trackMatrixPair; 

  trackMatrixPair trackRefit(const VertexState & vertex,
        RefCountedLinearizedTrackState linTrackState,
        float weight) const;

private:

  VertexTrackFactory<N> theVTFactory;
  KVFHelper<N> helper;
  KalmanVertexUpdator<N> updator;
};

#endif