ParticleKinematicLinearizedTrackState.h

Go to the documentation of this file.

#ifndef ParticleKinematicLinearizedTrackState_H
#define ParticleKinematicLinearizedTrackState_H

#include "RecoVertex/VertexPrimitives/interface/LinearizedTrackState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/RefCountedKinematicParticle.h"
#include "RecoVertex/VertexPrimitives/interface/RefittedTrackState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/PerigeeKinematicState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/TransientTrackKinematicStateBuilder.h"
#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"
#include "DataFormats/CLHEP/interface/Migration.h"
#include "RecoVertex/KinematicFitPrimitives/interface/Matrices.h"

class ParticleKinematicLinearizedTrackState : public LinearizedTrackState<6> {

public:
 
 friend class ParticleKinematicLinearizedTrackStateFactory;
 typedef ReferenceCountingPointer<LinearizedTrackState<6> > RefCountedLinearizedTrackState;

 ParticleKinematicLinearizedTrackState()
 {jacobiansAvailable = false;}
 
 virtual ReferenceCountingPointer<LinearizedTrackState<6> > stateWithNewLinearizationPoint
    (const GlobalPoint & newLP) const; 
    
 const GlobalPoint & linearizationPoint() const { return theLinPoint; }
 

 const AlgebraicVector6 & constantTerm() const;

 virtual const AlgebraicMatrix63 & positionJacobian() const;

 virtual const AlgebraicMatrix64 & momentumJacobian() const;

 const AlgebraicVector6 & parametersFromExpansion() const;

 AlgebraicVector6 predictedStateParameters() const;
    
 AlgebraicVectorM predictedStateMomentumParameters() const; 

 AlgebraicSymMatrix44 predictedStateMomentumError() const; 

 AlgebraicSymMatrix66 predictedStateWeight(int & error) const;
  
 AlgebraicSymMatrix66 predictedStateError() const; 

// /** 
//  * Method returning the impact point measurement     
//  */      
//  ImpactPointMeasurement impactPointMeasurement() const;
  
 TrackCharge charge() const;
 
 RefCountedKinematicParticle particle() const;
   
 bool operator ==(LinearizedTrackState<6>& other)const;
 
 bool hasError() const;
 
 RefCountedRefittedTrackState createRefittedTrackState(
                     const GlobalPoint & vertexPosition, 
             const AlgebraicVectorM & vectorParameters,
             const AlgebraicSymMatrix77 & covarianceMatrix) const;
             
  virtual AlgebraicVectorN refittedParamFromEquation(
    const RefCountedRefittedTrackState & theRefittedState) const;

  virtual inline void checkParameters(AlgebraicVectorN & parameters) const;
             
 double weightInMixture() const;

 std::vector<ReferenceCountingPointer<LinearizedTrackState<6> > > components()
                                    const;
 
 virtual reco::TransientTrack track() const;

 
private:

  ParticleKinematicLinearizedTrackState(const GlobalPoint & linP, RefCountedKinematicParticle & prt) 
              : theLinPoint(linP), part(prt), jacobiansAvailable(false),
          theCharge(prt->currentState().particleCharge())  ,impactPointAvailable(false)
       
    {}
    
  void computeJacobians() const;
  void computeChargedJacobians() const;
 
  void computeNeutralJacobians() const;   
     
     
  GlobalPoint theLinPoint;
  RefCountedKinematicParticle part;
  TransientTrackKinematicStateBuilder builder;
  
  mutable bool errorAvailable;
  mutable bool jacobiansAvailable;
  mutable AlgebraicMatrix63 thePositionJacobian;
  mutable AlgebraicMatrix64 theMomentumJacobian;
  mutable PerigeeKinematicState thePredState;
  mutable AlgebraicVector6 theConstantTerm;
  mutable AlgebraicVector6 theExpandedParams;
  
  TrackCharge theCharge;
//   ImpactPointMeasurementExtractor theIPMExtractor;
  mutable bool impactPointAvailable; 
};
#endif