PerigeeKinematicState.h

Go to the documentation of this file.

#ifndef PerigeeKinematicState_H
#define PerigeeKinematicState_H

#include "RecoVertex/KinematicFitPrimitives/interface/KinematicState.h"
#include "RecoVertex/KinematicFitPrimitives/interface/ExtendedPerigeeTrajectoryParameters.h"
#include "RecoVertex/KinematicFitPrimitives/interface/ExtendedPerigeeTrajectoryError.h"

class PerigeeKinematicState
{
public:

 PerigeeKinematicState()
 {
  vl = false;
  errorIsAvailable = false;
 }

 virtual ~PerigeeKinematicState(){}
 bool hasError() const
 {
  if(!(isValid()))throw VertexException("PerigeeKinematicState::error is requested for the invalid state");
  return errorIsAvailable;
 }
 
 bool isValid() const
 {return vl;}

 const KinematicState theState() const
 {
  if(!isValid()) throw VertexException("PerigeeKinematicState::initial state is requested for the invalid state");
  return inState;
 }

 const GlobalPoint referencePoint() const
 {
  if(!isValid()) throw VertexException("PerigeeKinematicState::point is requested for the invalid state");
  return point;
 }

 const ExtendedPerigeeTrajectoryError& perigeeError() const
 {
  if(!(isValid()))  throw VertexException("PerigeeKinematicState::requesting perigee error for invalid state");
  if(!(hasError())) throw VertexException("PerigeeKinematicState::requesting perigee error when none available");
  return cov;
 }

 const ExtendedPerigeeTrajectoryParameters & perigeeParameters() const
 {
  if(!(isValid()))  throw VertexException("PerigeeKinematicState::requesting perigee parameters for invalid state");
  return par;
 }

private:

 friend class TransientTrackKinematicStateBuilder;


 PerigeeKinematicState(const KinematicState& state, const GlobalPoint& pt);
/*
 AlgebraicMatrix jacobianKinematicToExPerigee(const KinematicState& state, 
                                              const GlobalPoint& pt)const;
 AlgebraicMatrix jacobianExPerigeeToKinematic(const ExtendedPerigeeTrajectoryParameters& state,
                                              const GlobalPoint& point)const;
*/                        
 AlgebraicMatrix jacobianCurvilinear2Perigee(const FreeTrajectoryState& fts) const;
 
private:
 GlobalPoint point;
 ExtendedPerigeeTrajectoryParameters par;
 ExtendedPerigeeTrajectoryError cov;
 KinematicState inState;
 bool errorIsAvailable;
 bool vl;

};
#endif