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