HGCScintillatorDetId.cc

Go to the documentation of this file.

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

const HGCScintillatorDetId HGCScintillatorDetId::Undefined(0, 0, 0, 0, false);

HGCScintillatorDetId::HGCScintillatorDetId() : DetId() {}

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

HGCScintillatorDetId::HGCScintillatorDetId(int type, int layer, int radius, int phi, bool trigger)
    : DetId(HGCalHSc, ForwardEmpty) {
  int zside = (radius < 0) ? 1 : 0;
  int itrig = trigger ? 1 : 0;
  int radiusAbs = std::abs(radius);
  id_ |= (((type & kHGCalTypeMask) << kHGCalTypeOffset) | ((zside & kHGCalZsideMask) << kHGCalZsideOffset) |
          ((itrig & kHGCalTriggerMask) << kHGCalTriggerOffset) | ((layer & kHGCalLayerMask) << kHGCalLayerOffset) |
          ((radiusAbs & kHGCalRadiusMask) << kHGCalRadiusOffset) | ((phi & kHGCalPhiMask) << kHGCalPhiOffset));
}

HGCScintillatorDetId::HGCScintillatorDetId(const DetId& gen) {
  if (!gen.null()) {
    if (gen.det() != HGCalHSc) {
      throw cms::Exception("Invalid DetId")
          << "Cannot initialize HGCScintillatorDetId from " << std::hex << gen.rawId() << std::dec;
    }
  }
  id_ = gen.rawId();
}

HGCScintillatorDetId& HGCScintillatorDetId::operator=(const DetId& gen) {
  if (!gen.null()) {
    if (gen.det() != HGCalHSc) {
      throw cms::Exception("Invalid DetId")
          << "Cannot assign HGCScintillatorDetId from " << std::hex << gen.rawId() << std::dec;
    }
  }
  id_ = gen.rawId();
  return (*this);
}

int HGCScintillatorDetId::iradiusAbs() const {
  if (trigger())
    return (2 * ((id_ >> kHGCalRadiusOffset) & kHGCalRadiusMask));
  else
    return ((id_ >> kHGCalRadiusOffset) & kHGCalRadiusMask);
}

int HGCScintillatorDetId::iradiusTriggerAbs() const {
  if (trigger())
    return ((iradiusAbs() + 1) / 2);
  else
    return iradiusAbs();
}

int HGCScintillatorDetId::iphi() const {
  if (trigger())
    return (2 * ((id_ >> kHGCalPhiOffset) & kHGCalPhiMask));
  else
    return ((id_ >> kHGCalPhiOffset) & kHGCalPhiMask);
}

int HGCScintillatorDetId::iphiTrigger() const {
  if (trigger())
    return ((iphi() + 1) / 2);
  else
    return iphi();
}

std::vector<HGCScintillatorDetId> HGCScintillatorDetId::detectorCells() const {
  std::vector<HGCScintillatorDetId> cells;
  int irad = iradiusAbs();
  int ifi = iphi();
  int iz = zside();
  if (trigger()) {
    cells.emplace_back(HGCScintillatorDetId(type(), layer(), (2 * irad - 1) * iz, 2 * ifi - 1, false));
    cells.emplace_back(HGCScintillatorDetId(type(), layer(), 2 * irad * iz, 2 * ifi - 1, false));
    cells.emplace_back(HGCScintillatorDetId(type(), layer(), (2 * irad - 1) * iz, 2 * ifi, false));
    cells.emplace_back(HGCScintillatorDetId(type(), layer(), 2 * irad * iz, 2 * ifi, false));
  } else {
    cells.emplace_back(HGCScintillatorDetId(type(), layer(), irad * iz, ifi, false));
  }
  return cells;
}

HGCScintillatorDetId HGCScintillatorDetId::geometryCell() const {
  if (trigger()) {
    return HGCScintillatorDetId(type(), layer(), iradiusTrigger(), iphiTrigger(), false);
  } else {
    return HGCScintillatorDetId(type(), layer(), iradius(), iphi(), false);
  }
}

HGCScintillatorDetId HGCScintillatorDetId::triggerCell() const {
  if (trigger())
    return HGCScintillatorDetId(type(), layer(), iradius(), iphi(), true);
  else
    return HGCScintillatorDetId(type(), layer(), iradiusTrigger(), iphiTrigger(), true);
}

std::ostream& operator<<(std::ostream& s, const HGCScintillatorDetId& id) {
  return s << " HGCScintillatorDetId::EE:HE= " << id.isEE() << ":" << id.isHE() << " trigger= " << id.trigger()
           << " type= " << id.type() << " layer= " << id.layer() << " radius= " << id.iradius() << ":"
           << id.iradiusTrigger() << " phi= " << id.iphi() << ":" << id.iphiTrigger();
}