ExtendedPerigeeTrajectoryError.h

Go to the documentation of this file.

#ifndef ExtendedPerigeeTrajectoryError_H
#define ExtendedPerigeeTrajectoryError_H
 
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
class ExtendedPerigeeTrajectoryError {
public:
  ExtendedPerigeeTrajectoryError() : weightAvailable(false), vl(false) {}
 
  ExtendedPerigeeTrajectoryError(const AlgebraicSymMatrix66& covariance)
      : cov(covariance), weightAvailable(false), vl(true) {}
 
  bool isValid() const { return vl; }
 
  bool weightIsAvailable() const { return weightAvailable; }
 
  const AlgebraicSymMatrix66& covarianceMatrix() const { return cov; }
 
  const AlgebraicSymMatrix66& weightMatrix(int& error) const {
    error = 0;
    if (!weightIsAvailable()) {
      weight = cov.Inverse(error);
      if (error != 0)
        LogDebug("RecoVertex/ExtendedPerigeeTrajectoryError") << "unable to invert covariance matrix\n";
      weightAvailable = true;
    }
    return weight;
  }
 
private:
  AlgebraicSymMatrix66 cov;
  mutable AlgebraicSymMatrix66 weight;
  mutable bool weightAvailable;
  mutable bool vl;
};
#endif