EmulatedME0Segment.cc

Go to the documentation of this file.

 
#include <DataFormats/MuonReco/interface/EmulatedME0Segment.h>
#include <iostream>
 
EmulatedME0Segment::EmulatedME0Segment(const LocalPoint& origin,
                                       const LocalVector& direction,
                                       const AlgebraicSymMatrix& errors,
                                       const double chi2)
    : theOrigin(origin), theLocalDirection(direction), theCovMatrix(errors), theChi2(chi2) {}
 
EmulatedME0Segment::~EmulatedME0Segment() {}
 
LocalError EmulatedME0Segment::localPositionError() const {
  return LocalError(theCovMatrix[2][2], theCovMatrix[2][3], theCovMatrix[3][3]);
}
 
LocalError EmulatedME0Segment::localDirectionError() const {
  return LocalError(theCovMatrix[0][0], theCovMatrix[0][1], theCovMatrix[1][1]);
}
 
AlgebraicVector EmulatedME0Segment::parameters() const {
  // For consistency with DT and what we require for the TrackingRecHit interface,
  // the order of the parameters in the returned vector should be (dx/dz, dy/dz, x, z)
 
  AlgebraicVector result(4);
 
  result[0] = theLocalDirection.x() / theLocalDirection.z();
  result[1] = theLocalDirection.y() / theLocalDirection.z();
  result[2] = theOrigin.x();
  result[3] = theOrigin.y();
 
  return result;
}
 
namespace {
  AlgebraicMatrix createStaticMatrix() {
    AlgebraicMatrix m(4, 5, 0);
    m[0][1] = 1;
    m[1][2] = 1;
    m[2][3] = 1;
    m[3][4] = 1;
    return m;
  }
};  // namespace
 
namespace {
  const AlgebraicMatrix theProjectionMatrix = createStaticMatrix();
};
 
AlgebraicMatrix EmulatedME0Segment::projectionMatrix() const { return theProjectionMatrix; }
 
//
void EmulatedME0Segment::print() const { std::cout << *this << std::endl; }
 
std::ostream& operator<<(std::ostream& os, const EmulatedME0Segment& seg) {
  os << "EmulatedME0Segment: local pos = " << seg.localPosition() << " posErr = ("
     << sqrt(seg.localPositionError().xx()) << "," << sqrt(seg.localPositionError().yy()) << "0,)\n"
     << "            dir = " << seg.localDirection() << " dirErr = (" << sqrt(seg.localDirectionError().xx()) << ","
     << sqrt(seg.localDirectionError().yy()) << "0,)\n";
  return os;
}