LocalTrajectoryError.h

Go to the documentation of this file.

#ifndef _TRACKER_LOCALTRAJECTORYERROR_H_
#define _TRACKER_LOCALTRAJECTORYERROR_H_

#include "DataFormats/GeometrySurface/interface/LocalError.h"
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"
#include <memory>

class LocalTrajectoryError {
public:
  // construct
  LocalTrajectoryError() {}

  LocalTrajectoryError(InvalidError) : theCovarianceMatrix(ROOT::Math::SMatrixNoInit()) {
    theCovarianceMatrix(0, 0) = -99999.e10;
  }
  // destruct
  ~LocalTrajectoryError() {}

  bool invalid() const { return theCovarianceMatrix(0, 0) < -1.e10; }
  bool valid() const { return !invalid(); }

  // not really full check of posdef
  bool posDef() const {
    return (theCovarianceMatrix(0, 0) >= 0) && (theCovarianceMatrix(1, 1) >= 0) && (theCovarianceMatrix(2, 2) >= 0) &&
           (theCovarianceMatrix(3, 3) >= 0) && (theCovarianceMatrix(4, 4) >= 0);
  }

  LocalTrajectoryError(const AlgebraicSymMatrix55 &aCovarianceMatrix)
      : theCovarianceMatrix(aCovarianceMatrix), theWeightMatrixPtr() {}

  LocalTrajectoryError(float dx, float dy, float dxdir, float dydir, float dpinv);

  // access

  const AlgebraicSymMatrix55 &matrix() const { return theCovarianceMatrix; }

  const AlgebraicSymMatrix55 &weightMatrix() const;

  LocalTrajectoryError &operator*=(double factor) {
    theCovarianceMatrix *= factor;
    if ((theWeightMatrixPtr.get() != nullptr) && (factor != 0.0)) {
      (*theWeightMatrixPtr) /= factor;
    }
    return *this;
  }

  LocalError positionError() const {
    return LocalError(theCovarianceMatrix(3, 3), theCovarianceMatrix(3, 4), theCovarianceMatrix(4, 4));
  }

private:
  AlgebraicSymMatrix55 theCovarianceMatrix;
  mutable std::shared_ptr<AlgebraicSymMatrix55> theWeightMatrixPtr;
};

#endif