PerigeeRefittedTrackState.h

Go to the documentation of this file.

#ifndef PerigeeRefittedTrackState_H
#define PerigeeRefittedTrackState_H

#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "RecoVertex/VertexPrimitives/interface/RefittedTrackState.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateClosestToPoint.h"

class FreeTrajectoryState;
class TrajectoryStateOnSurface;
class Surface;
class Propagator;

class PerigeeRefittedTrackState : public RefittedTrackState<5> {

public:

  typedef ReferenceCountingPointer<RefittedTrackState<5> > RefCountedRefittedTrackState;

  PerigeeRefittedTrackState(const TrajectoryStateClosestToPoint & tscp,
                const AlgebraicVector3 & aMomentumAtVertex,
                const double aWeight = 1.) :
    theState(tscp), momentumAtVertex(aMomentumAtVertex), theWeight(aWeight) {}

 virtual ~PerigeeRefittedTrackState(){}

  virtual FreeTrajectoryState freeTrajectoryState() const
    {return theState.theState();}

  virtual TrajectoryStateOnSurface trajectoryStateOnSurface(
        const Surface & surface) const;

  virtual TrajectoryStateOnSurface trajectoryStateOnSurface(
        const Surface & surface, const Propagator & propagator) const;

  virtual AlgebraicVector5 parameters() const
    {return theState.perigeeParameters().vector();}

  virtual AlgebraicSymMatrix55  covariance() const
    {return theState.perigeeError().covarianceMatrix();}

  virtual GlobalPoint position() const
    {return theState.referencePoint();}

  virtual AlgebraicVector3 momentumVector() const;

  virtual double weight() const {return theWeight;}

  virtual ReferenceCountingPointer<RefittedTrackState<5> > stateWithNewWeight
    (const double newWeight) const;

  virtual std::vector<ReferenceCountingPointer<RefittedTrackState<5> > > components() const;

  virtual reco::TransientTrack transientTrack() const;

private:

  TrajectoryStateClosestToPoint theState;
  AlgebraicVector3 momentumAtVertex;
  double theWeight;
};
#endif