PerigeeTrajectoryError.h

Go to the documentation of this file.

#ifndef PerigeeTrajectoryError_H
#define PerigeeTrajectoryError_H

#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "TrackingTools/TrajectoryParametrization/interface/TrajectoryStateExceptions.h"

class PerigeeTrajectoryError
{

public:

  PerigeeTrajectoryError() {}
  ~PerigeeTrajectoryError() {}

  PerigeeTrajectoryError(AlgebraicSymMatrix aPerigeeError):
    thePerigeeError(asSMatrix<5>(aPerigeeError)), weightIsAvailable(false) {}

  PerigeeTrajectoryError(const AlgebraicSymMatrix55 &aPerigeeError):
    thePerigeeError(aPerigeeError), weightIsAvailable(false) {
         
  }


  const AlgebraicSymMatrix covarianceMatrix_old() const {return asHepMatrix(thePerigeeError);}
  const AlgebraicSymMatrix55 & covarianceMatrix() const {return thePerigeeError;}


  const AlgebraicSymMatrix weightMatrix_old(int & error) const {
    if (!weightIsAvailable) calculateWeightMatrix();
    error = inverseError;
    return asHepMatrix(thePerigeeWeight);
  }
  const AlgebraicSymMatrix55 &weightMatrix(int & error) const
  {
    if (!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