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 <boost/shared_ptr.hpp>

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() != 0) && (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 boost::shared_ptr<AlgebraicSymMatrix55> theWeightMatrixPtr;
};

#endif