AHCalDetId.cc

Go to the documentation of this file.

#include "SimG4CMS/HGCalTestBeam/interface/AHCalDetId.h"
#include <ostream>
#include "FWCore/Utilities/interface/Exception.h"

const AHCalDetId AHCalDetId::Undefined(0, 0, 0);

AHCalDetId::AHCalDetId() : DetId() {}

AHCalDetId::AHCalDetId(uint32_t rawid) : DetId(rawid) {}

AHCalDetId::AHCalDetId(int row, int col, int depth) : DetId(Hcal, HcalOther) {
  int icol = (col > 0) ? col : 10 - col;
  int irow = (row > 0) ? row : 10 - row;
  id_ |=
      (HcalDetId::kHcalIdFormat2) |
      ((depth & HcalDetId::kHcalDepthMask2) << HcalDetId::kHcalDepthOffset2) |
      (HcalDetId::kHcalZsideMask2) |
      ((irow & HcalDetId::kHcalEtaMask2) << HcalDetId::kHcalEtaOffset2) |
      (icol & HcalDetId::kHcalPhiMask2);
}

AHCalDetId::AHCalDetId(const DetId& gen) {
  if (!gen.null()) {
    HcalSubdetector subdet = (HcalSubdetector(gen.subdetId()));
    if (gen.det() != Hcal || subdet != HcalOther) {
      throw cms::Exception("Invalid DetId")
          << "Cannot initialize AHCalDetId from " << std::hex << gen.rawId()
          << std::dec;
    }
  }
  id_ = gen.rawId();
}

int AHCalDetId::irow() const {
  int value = ((id_ >> HcalDetId::kHcalEtaOffset2) & HcalDetId::kHcalEtaMask2);
  if (value >= 10) value = -(value % 10);
  return value;
}

int AHCalDetId::icol() const {
  int value = (id_ & HcalDetId::kHcalPhiMask2);
  if (value >= 10) value = -(value % 10);
  return value;
}

int AHCalDetId::depth() const {
  return ((id_ >> HcalDetId::kHcalDepthOffset2) & HcalDetId::kHcalDepthMask2);
}

std::pair<double, double> AHCalDetId::getXY() const {
  int row = irow();
  int col = icol();
  double shiftx = (col > 0) ? -0.5 * deltaX_ : 0.5 * deltaX_;
  double shifty = (row > 0) ? -0.5 * deltaY_ : 0.5 * deltaY_;
  return std::pair<double, double>(col * deltaX_ + shiftx,
                                   row * deltaY_ + shifty);
}

double AHCalDetId::getZ() const {
  int lay = depth();
  return (zFirst_ + (lay - 1) * deltaZ_);
}

std::ostream& operator<<(std::ostream& s, const AHCalDetId& id) {
  return s << "(AHCal " << id.irow() << ',' << id.icol() << ',' << id.depth()
           << ')';
}