HBHELinearMap.cc

Go to the documentation of this file.

#include <cassert>
#include "FWCore/Utilities/interface/Exception.h"
#include <algorithm>

#include "CondFormats/HcalObjects/interface/HBHELinearMap.h"

void HBHELinearMap::getChannelTriple(const unsigned index, unsigned* depth, int* ieta, unsigned* iphi) const {
  if (index >= ChannelCount)
    throw cms::Exception(
        "In HBHELinearMap::getChannelTriple: "
        "input index out of range");
  const HBHEChannelId& id = lookup_[index];
  if (depth)
    *depth = id.depth();
  if (ieta)
    *ieta = id.ieta();
  if (iphi)
    *iphi = id.iphi();
}

unsigned HBHELinearMap::find(const unsigned depth, const int ieta, const unsigned iphi) const {
  const HBHEChannelId id(depth, ieta, iphi);
  const unsigned loc = std::lower_bound(inverse_.begin(), inverse_.end(), MapPair(id, 0U)) - inverse_.begin();
  if (loc < ChannelCount)
    if (inverse_[loc].first == id)
      return inverse_[loc].second;
  return ChannelCount;
}

bool HBHELinearMap::isValidTriple(const unsigned depth, const int ieta, const unsigned iphi) const {
  const unsigned ind = find(depth, ieta, iphi);
  return ind < ChannelCount;
}

unsigned HBHELinearMap::linearIndex(const unsigned depth, const int ieta, const unsigned iphi) const {
  const unsigned ind = find(depth, ieta, iphi);
  if (ind >= ChannelCount)
    throw cms::Exception(
        "In HBHELinearMap::linearIndex: "
        "invalid channel triple");
  return ind;
}

HBHELinearMap::HBHELinearMap() {
  unsigned l = 0;
  unsigned depth = 1;

  for (int ieta = -29; ieta <= -21; ++ieta)
    for (unsigned iphi = 1; iphi < 72; iphi += 2)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = -20; ieta <= 20; ++ieta)
    if (ieta)
      for (unsigned iphi = 1; iphi <= 72; ++iphi)
        lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = 21; ieta <= 29; ++ieta)
    for (unsigned iphi = 1; iphi < 72; iphi += 2)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  depth = 2;

  for (int ieta = -29; ieta <= -21; ++ieta)
    for (unsigned iphi = 1; iphi < 72; iphi += 2)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = -20; ieta <= -18; ++ieta)
    for (unsigned iphi = 1; iphi <= 72; ++iphi)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = -16; ieta <= -15; ++ieta)
    for (unsigned iphi = 1; iphi <= 72; ++iphi)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = 15; ieta <= 16; ++ieta)
    for (unsigned iphi = 1; iphi <= 72; ++iphi)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = 18; ieta <= 20; ++ieta)
    for (unsigned iphi = 1; iphi <= 72; ++iphi)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = 21; ieta <= 29; ++ieta)
    for (unsigned iphi = 1; iphi < 72; iphi += 2)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  depth = 3;

  for (int ieta = -28; ieta <= -27; ++ieta)
    for (unsigned iphi = 1; iphi < 72; iphi += 2)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = -16; ieta <= -16; ++ieta)
    for (unsigned iphi = 1; iphi <= 72; ++iphi)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = 16; ieta <= 16; ++ieta)
    for (unsigned iphi = 1; iphi <= 72; ++iphi)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  for (int ieta = 27; ieta <= 28; ++ieta)
    for (unsigned iphi = 1; iphi < 72; iphi += 2)
      lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

  assert(l == ChannelCount);

  inverse_.reserve(ChannelCount);
  for (unsigned i = 0; i < ChannelCount; ++i)
    inverse_.push_back(MapPair(lookup_[i], i));
  std::sort(inverse_.begin(), inverse_.end());
}

HcalSubdetector HBHELinearMap::getSubdetector(const unsigned depth, const int ieta) {
  const unsigned abseta = std::abs(ieta);

  // Make sure the arguments are in range
  if (!(abseta <= 29U))
    throw cms::Exception(
        "In HBHELinearMap::getSubdetector: "
        "eta argument out of range");
  if (!(depth > 0U && depth < 4U))
    throw cms::Exception(
        "In HBHELinearMap::getSubdetector: "
        "depth argument out of range");
  if (abseta == 29U)
    if (!(depth <= 2U))
      throw cms::Exception(
          "In HBHELinearMap::getSubdetector: "
          "depth argument out of range "
          "for |ieta| = 29");
  if (abseta <= 15U)
    return HcalBarrel;
  else if (abseta == 16U)
    return depth <= 2U ? HcalBarrel : HcalEndcap;
  else
    return HcalEndcap;
}

const HBHELinearMap& hbheChannelMap() {
  static const HBHELinearMap chMap;
  return chMap;
}