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