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