DTRecHit1DPair.cc

Go to the documentation of this file.

#include "DataFormats/DTRecHit/interface/DTRecHit1DPair.h"
#include "DataFormats/DTDigi/interface/DTDigi.h"

#include "FWCore/Utilities/interface/Exception.h"

using namespace DTEnums;
using namespace std;

// Constructor without components: must use setPos and Err!
DTRecHit1DPair::DTRecHit1DPair(const DTWireId& wireId, const DTDigi& digi)
    : theLeftHit(wireId, Left, digi.time()), theRightHit(wireId, Right, digi.time()) {}

// Default constructor
DTRecHit1DPair::DTRecHit1DPair() : theLeftHit(), theRightHit() {}

// Destructor
DTRecHit1DPair::~DTRecHit1DPair() {}

DTRecHit1DPair* DTRecHit1DPair::clone() const { return new DTRecHit1DPair(*this); }

// Return the 3-dimensional local position.
// The average theLeftHit/theRightHit hits position, namely the wire position
// is returned.
LocalPoint DTRecHit1DPair::localPosition() const {
  return theLeftHit.localPosition() + (theRightHit.localPosition() - theLeftHit.localPosition()) / 2.;
}

// Return the 3-dimensional error on the local position.
// The error is defiened as half
// the distance between theLeftHit and theRightHit pos
LocalError DTRecHit1DPair::localPositionError() const {
  return LocalError((theRightHit.localPosition().x() - theLeftHit.localPosition().x()) / 2., 0., 0.);
}

// Access to component RecHits.
vector<const TrackingRecHit*> DTRecHit1DPair::recHits() const {
  vector<const TrackingRecHit*> result;
  result.push_back(componentRecHit(Left));
  result.push_back(componentRecHit(Right));
  return result;
}

// Non-const access to component RecHits.
vector<TrackingRecHit*> DTRecHit1DPair::recHits() {
  vector<TrackingRecHit*> result;
  result.push_back(componentRecHit(Left));
  result.push_back(componentRecHit(Right));
  return result;
}

// Return the detId of the Det (a DTLayer).
DetId DTRecHit1DPair::geographicalId() const { return wireId().layerId(); }

// Comparison operator, based on the wireId and the digi time
bool DTRecHit1DPair::operator==(const DTRecHit1DPair& hit) const {
  return wireId() == hit.wireId() && fabs(digiTime() - hit.digiTime()) < 0.1;
}

// Return position in the local (layer) coordinate system for a
// certain hypothesis about the L/R cell side
LocalPoint DTRecHit1DPair::localPosition(DTCellSide lrside) const { return componentRecHit(lrside)->localPosition(); }

// Return position error in the local (layer) coordinate system for a
// certain hypothesis about the L/R cell side
LocalError DTRecHit1DPair::localPositionError(DTCellSide lrside) const {
  return componentRecHit(lrside)->localPositionError();
}

// Set the 3-dimensional local position for the component hit
// corresponding to the given cell side. Default value is assumed for the error.
void DTRecHit1DPair::setPosition(DTCellSide lrside, const LocalPoint& point) {
  if (lrside != undefLR)
    componentRecHit(lrside)->setPosition(point);
  else
    throw cms::Exception("DTRecHit1DPair::setPosition with undefined LR");
}

// Set the 3-dimensional local position and error for the component hit
// corresponding to the given cell side. Default value is assumed for the error.
void DTRecHit1DPair::setPositionAndError(DTCellSide lrside, const LocalPoint& point, const LocalError& err) {
  if (lrside != undefLR) {
    componentRecHit(lrside)->setPosition(point);
    componentRecHit(lrside)->setError(err);
  } else
    throw cms::Exception("DTRecHit1DPair::setPosition with undefined LR");
}

// Return the left/right DTRecHit1D
const DTRecHit1D* DTRecHit1DPair::componentRecHit(DTCellSide lrSide) const {
  if (lrSide == Left) {
    return &theLeftHit;
  } else if (lrSide == Right) {
    return &theRightHit;
  } else {
    throw cms::Exception("DTRecHit1DPair::recHit with undefined LR");
  }
}

// Non const access to left/right DTRecHit1D
DTRecHit1D* DTRecHit1DPair::componentRecHit(DTCellSide lrSide) {
  if (lrSide == Left) {
    return &theLeftHit;
  } else if (lrSide == Right) {
    return &theRightHit;
  } else {
    throw cms::Exception("DTRecHit1DPair::recHit with undefined LR");
  }
}

pair<const DTRecHit1D*, const DTRecHit1D*> DTRecHit1DPair::componentRecHits() const {
  return make_pair(componentRecHit(Left), componentRecHit(Right));
}

// Ostream operator
ostream& operator<<(ostream& os, const DTRecHit1DPair& hit) {
  os << "Pos: " << hit.localPosition();
  return os;
}