PerigeeMultiLTS.h

Go to the documentation of this file.

#ifndef PerigeeMultiLTS_H
#define PerigeeMultiLTS_H

#include "RecoVertex/VertexPrimitives/interface/LinearizedTrackState.h"
#include "RecoVertex/VertexTools/interface/LinearizedTrackStateFactory.h"
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include <vector>

class PerigeeMultiLTS : public LinearizedTrackState<5> {


public:

  typedef ReferenceCountingPointer<LinearizedTrackState<5> > RefCountedLinearizedTrackState;

  friend class MultiPerigeeLTSFactory;

  virtual  RefCountedLinearizedTrackState stateWithNewLinearizationPoint
    (const GlobalPoint & newLP) const;


  const GlobalPoint & linearizationPoint() const { return theLinPoint; }

  virtual reco::TransientTrack track() const { return theTrack; }

  const TrajectoryStateOnSurface state() const { return theTSOS; }

  virtual const AlgebraicVectorN & constantTerm() const;

  virtual const AlgebraicMatrixN3 & positionJacobian() const;

  virtual const AlgebraicMatrixNM & momentumJacobian() const;

  virtual const AlgebraicVectorN & parametersFromExpansion() const;

  const TrajectoryStateClosestToPoint & predictedState() const;

  AlgebraicVectorN predictedStateParameters() const;

  virtual AlgebraicVectorM predictedStateMomentumParameters() const;

  AlgebraicSymMatrixNN predictedStateWeight(int & error) const;

  AlgebraicSymMatrixNN predictedStateError() const;

  AlgebraicSymMatrixMM predictedStateMomentumError() const;

  TrackCharge charge() const {return theCharge;}

  bool hasError() const;

  bool operator ==(LinearizedTrackState<5>& other)const;

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

  virtual AlgebraicVector5 refittedParamFromEquation(
    const RefCountedRefittedTrackState & theRefittedState) const;

  virtual void inline checkParameters(AlgebraicVector5 & parameters) const;
  virtual double weightInMixture() const {return theTSOS.weight();}

  virtual std::vector<ReferenceCountingPointer<LinearizedTrackState<5> > > components()
                const {return ltComp;}

private:

   PerigeeMultiLTS(const GlobalPoint & linP, const reco::TransientTrack & track,
    const TrajectoryStateOnSurface& tsos);

  void prepareCollapsedState() const;

  GlobalPoint theLinPoint;
  reco::TransientTrack theTrack;


  const TrajectoryStateOnSurface theTSOS;
  std::vector<RefCountedLinearizedTrackState> ltComp;
  mutable RefCountedLinearizedTrackState collapsedStateLT;
  LinearizedTrackStateFactory theLTSfactory;
  TrackCharge theCharge;
  mutable bool collapsedStateAvailable;
};

#endif