HGCHEDetId.cc

Go to the documentation of this file.

#include "DataFormats/ForwardDetId/interface/HGCHEDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <ostream>

const HGCHEDetId HGCHEDetId::Undefined(ForwardEmpty, 0, 0, 0, 0, 0);

HGCHEDetId::HGCHEDetId() : DetId() {}

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

HGCHEDetId::HGCHEDetId(ForwardSubdetector subdet, int zp, int lay, int sec, int subsec, int cell)
    : DetId(Forward, subdet) {
  id_ |= ((cell & kHGCHECellMask) << kHGCHECellOffset);
  id_ |= ((sec & kHGCHESectorMask) << kHGCHESectorOffset);
  if (subsec < 0)
    subsec = 0;
  id_ |= ((subsec & kHGCHESubSectorMask) << kHGCHESubSectorOffset);
  id_ |= ((lay & kHGCHELayerMask) << kHGCHELayerOffset);
  if (zp > 0)
    id_ |= ((zp & kHGCHEZsideMask) << kHGCHEZsideOffset);
}

HGCHEDetId::HGCHEDetId(const DetId& gen) {
  if (!gen.null()) {
    ForwardSubdetector subdet = (ForwardSubdetector(gen.subdetId()));
    if ((gen.det() != Forward) || (subdet != HGCHEF && subdet != HGCHEB)) {
      throw cms::Exception("Invalid DetId") << "Cannot initialize HGCHEDetId from " << std::hex << gen.rawId()
                                            << std::dec << " Det|SubDet " << gen.det() << "|" << subdet;
    }
  }
  id_ = gen.rawId();
}

HGCHEDetId& HGCHEDetId::operator=(const DetId& gen) {
  if (!gen.null()) {
    ForwardSubdetector subdet = (ForwardSubdetector(gen.subdetId()));
    if ((gen.det() != Forward) || (subdet != HGCHEF && subdet != HGCHEB)) {
      throw cms::Exception("Invalid DetId") << "Cannot assign HGCHEDetId from " << std::hex << gen.rawId() << std::dec
                                            << " Det|SubDet " << gen.det() << "|" << subdet;
    }
  }
  id_ = gen.rawId();
  return (*this);
}

HGCHEDetId HGCHEDetId::geometryCell() const {
  int sub = ((subdet() == HGCHEF) ? 0 : ((id_ >> kHGCHESubSectorOffset) & kHGCHESubSectorMask));
  return HGCHEDetId(subdet(), zside(), layer(), sector(), sub, 0);
}

std::ostream& operator<<(std::ostream& s, const HGCHEDetId& id) {
  if (id.subdet() == HGCHEF || id.subdet() == HGCHEB) {
    return s << "isHE=" << id.isHE() << " zpos=" << id.zside() << " layer=" << id.layer()
             << " phi subSector=" << id.subsector() << " sector=" << id.sector() << " cell=" << id.cell();
  } else {
    return s << std::hex << id.rawId() << std::dec;
  }
}