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