/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/RecoVertex/VertexTools/interface/PerigeeLinearizedTrackState.h

Go to the documentation of this file.

#ifndef PerigeeLinearizedTrackState_H
#define PerigeeLinearizedTrackState_H

#include "RecoVertex/VertexPrimitives/interface/LinearizedTrackState.h"
#include "TrackingTools/PatternTools/interface/TSCPBuilderNoMaterial.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "RecoVertex/VertexPrimitives/interface/LinearizedTrackState.h"
#include "Math/SMatrix.h"
#include "DataFormats/CLHEP/interface/Migration.h"

class PerigeeLinearizedTrackState : public LinearizedTrackState<5> {


public:

  friend class LinearizedTrackStateFactory;
  typedef ReferenceCountingPointer<LinearizedTrackState<5> > RefCountedLinearizedTrackState;

  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; }

  const AlgebraicVector5 & constantTerm() const;

  const AlgebraicMatrix53 & positionJacobian() const;

  const AlgebraicMatrix53 & momentumJacobian() const;

  const AlgebraicVector5 & parametersFromExpansion() const;

  const TrajectoryStateClosestToPoint & predictedState() const;

  AlgebraicVector5 predictedStateParameters() const;

  virtual AlgebraicVector3 predictedStateMomentumParameters() const;

  AlgebraicSymMatrix55 predictedStateWeight(int & error) const;

  AlgebraicSymMatrix55 predictedStateError() const;

  AlgebraicSymMatrix33 predictedStateMomentumError() const;

//   /** Method returning the impact point measurement
//    */
//   ImpactPointMeasurement impactPointMeasurement() const;

  TrackCharge charge() const {return theCharge;}

  bool hasError() const;

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

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

  virtual RefCountedRefittedTrackState createRefittedTrackState(
        const GlobalPoint & vertexPosition,
        const AlgebraicVector3 & vectorParameters,
        const AlgebraicSymMatrix66 & covarianceMatrix) const;


  virtual AlgebraicVector5 refittedParamFromEquation(
        const RefCountedRefittedTrackState & theRefittedState) const;

  virtual double weightInMixture() const {return theTSOS.weight();}

  virtual void inline checkParameters(AlgebraicVector5 & parameters) const;

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

  virtual bool isValid() const;

private:

   PerigeeLinearizedTrackState(const GlobalPoint & linP, const reco::TransientTrack & track, 
        const TrajectoryStateOnSurface& tsos)
     : theLinPoint(linP), theTrack(track), jacobiansAvailable(false),
       // impactPointAvailable(false), 
       theCharge(theTrack.charge()), theTSOS(tsos) {}

  void computeJacobians() const;
  void computeChargedJacobians() const;
  void computeNeutralJacobians() const;
  void compute3DImpactPoint() const;

  GlobalPoint theLinPoint;
  reco::TransientTrack theTrack;

  mutable bool jacobiansAvailable;
  mutable AlgebraicMatrix53 thePositionJacobian, theMomentumJacobian;
  mutable TrajectoryStateClosestToPoint thePredState;
  mutable AlgebraicVector5 theConstantTerm;
  mutable AlgebraicVector5 theExpandedParams;

//   ImpactPointMeasurementExtractor theIPMExtractor;
  TSCPBuilderNoMaterial builder;
//   mutable bool impactPointAvailable;
//   mutable ImpactPointMeasurement the3DImpactPoint;
  TrackCharge theCharge;
  const TrajectoryStateOnSurface theTSOS;
};

#endif