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;

   RefCountedLinearizedTrackState stateWithNewLinearizationPoint
    (const GlobalPoint & newLP) const override;


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

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

  const TrajectoryStateOnSurface state() const { return theTSOS; }

  const AlgebraicVectorN & constantTerm() const override;

  const AlgebraicMatrixN3 & positionJacobian() const override;

  const AlgebraicMatrixNM & momentumJacobian() const override;

  const AlgebraicVectorN & parametersFromExpansion() const override;

  const TrajectoryStateClosestToPoint & predictedState() const;

  AlgebraicVectorN predictedStateParameters() const override;

  AlgebraicVectorM predictedStateMomentumParameters() const override;

  AlgebraicSymMatrixNN predictedStateWeight(int & error) const override;

  AlgebraicSymMatrixNN predictedStateError() const override;

  AlgebraicSymMatrixMM predictedStateMomentumError() const override;

  TrackCharge charge() const override {return theCharge;}

  bool hasError() const override;

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

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

  AlgebraicVector5 refittedParamFromEquation(
    const RefCountedRefittedTrackState & theRefittedState) const override;

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

  std::vector<ReferenceCountingPointer<LinearizedTrackState<5> > > components()
                const override {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