HBHERecHit.cc

Go to the documentation of this file.

#include "DataFormats/HcalRecHit/interface/HBHERecHit.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHitAuxSetter.h"
#include "DataFormats/HcalRecHit/interface/CaloRecHitAuxSetter.h"

HcalDetId HBHERecHit::idFront() const {
  if (auxPhase1_ & (1U << HBHERecHitAuxSetter::OFF_COMBINED)) {
    const HcalDetId myId(id());
    return HcalDetId(myId.subdet(), myId.ieta(), myId.iphi(), auxHBHE_ & 0xf);
  } else {
    return id();
  }
}

bool HBHERecHit::isMerged() const { return auxPhase1_ & (1U << HBHERecHitAuxSetter::OFF_COMBINED); }

std::ostream& operator<<(std::ostream& s, const HBHERecHit& hit) {
  s << hit.id() << ": " << hit.energy() << " GeV";
  if (hit.eraw() > -0.9e21) {
    s << ", eraw=" << hit.eraw() << " GeV";
  }
  if (hit.eaux() > -0.9e21) {
    s << ", eaux=" << hit.eaux() << " GeV";
  }
  if (hit.time() > -998) {
    s << ", t= " << hit.time() << " to " << hit.timeFalling() << " ns";
  }
  return s;
}

void HBHERecHit::getMergedIds(std::vector<HcalDetId>* ids) const {
  if (ids) {
    ids->clear();
    if (auxPhase1_ & (1U << HBHERecHitAuxSetter::OFF_COMBINED)) {
      const unsigned nMerged = CaloRecHitAuxSetter::getField(
          auxPhase1_, HBHERecHitAuxSetter::MASK_NSAMPLES, HBHERecHitAuxSetter::OFF_NSAMPLES);
      ids->reserve(nMerged);
      const HcalDetId myId(id());
      for (unsigned i = 0; i < nMerged; ++i) {
        const unsigned depth = CaloRecHitAuxSetter::getField(auxHBHE_, 0xf, i * 4);
        ids->emplace_back(myId.subdet(), myId.ieta(), myId.iphi(), depth);
      }
    }
  }
}