PerigeeTrajectoryError.h

Go to the documentation of this file.

#ifndef PerigeeTrajectoryError_H
#define PerigeeTrajectoryError_H

#include "TrackingTools/TrajectoryParametrization/interface/TrajectoryStateExceptions.h"
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"
#include "FWCore/Utilities/interface/Likely.h"
class PerigeeTrajectoryError {
public:
  PerigeeTrajectoryError() {}
  // ~PerigeeTrajectoryError() {}

  /*
  PerigeeTrajectoryError(AlgebraicSymMatrix aPerigeeError):
    thePerigeeError(asSMatrix<5>(aPerigeeError)), weightIsAvailable(false) {}
  */
  PerigeeTrajectoryError(const AlgebraicSymMatrix55 &aPerigeeError)
      : thePerigeeError(aPerigeeError), weightIsAvailable(false) {}

  const AlgebraicSymMatrix55 &covarianceMatrix() const { return thePerigeeError; }

  const AlgebraicSymMatrix55 &weightMatrix(int &error) const {
    if
      UNLIKELY(!weightIsAvailable) calculateWeightMatrix();
    error = inverseError;
    return thePerigeeWeight;
  }

  void calculateWeightMatrix() const;

  double transverseCurvatureError() const { return sqrt(thePerigeeError(0, 0)); }

  double thetaError() const { return sqrt(thePerigeeError(1, 1)); }

  double phiError() const { return sqrt(thePerigeeError(2, 2)); }

  double transverseImpactParameterError() const { return sqrt(thePerigeeError(3, 3)); }

  double longitudinalImpactParameterError() const { return sqrt(thePerigeeError(4, 4)); }

private:
  AlgebraicSymMatrix55 thePerigeeError;
  mutable AlgebraicSymMatrix55 thePerigeeWeight;
  mutable int inverseError;
  mutable bool weightIsAvailable;
};
#endif