DD4hep_MuonNumbering.cc

Go to the documentation of this file.

#include "Geometry/MuonNumbering/interface/DD4hep_MuonNumbering.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace cms;
using namespace edm;

const MuonBaseNumber MuonNumbering::geoHistoryToBaseNumber(const cms::ExpandedNodes& nodes) const {
  MuonBaseNumber num;

  int levelPart = get("level");
  int superPart = get("super");
  int basePart = get("base");
  int startCopyNo = get("xml_starts_with_copyno");

  // some consistency checks
  if (basePart != 1) {
    edm::LogError("Geometry") << "MuonNumbering finds unusual base constant: " << basePart;
  }
  if (superPart < 100) {
    edm::LogError("Geometry") << "MuonNumbering finds unusual super part: " << superPart
                              << " -- resetting super part to 100";
    superPart = 100;  // Reset to sensible value so calculations below don't go out of range
  }
  if (levelPart < 10 * superPart) {
    edm::LogError("Geometry") << "MuonNumbering finds unusual level constant: " << levelPart;
  }
  if ((startCopyNo != 0) && (startCopyNo != 1)) {
    edm::LogError("Geometry") << "MuonNumbering finds unusual start value for copy numbers: " << startCopyNo;
  }
  int ctr(0);
  for (auto const& it : nodes.tags) {
    int tag = it / levelPart;
    if (tag > 0) {
      int offset = nodes.offsets[ctr];
      int copyno = nodes.copyNos[ctr] + offset % superPart;
      int super = offset / superPart;
      num.addBase(tag, super, copyno - startCopyNo);
    }
    ++ctr;
  }
  return num;
}

const int MuonNumbering::get(const char* key) const {
  int result(0);
  auto const& it = values_.find(key);
  if (it != end(values_))
    result = it->second;
  return result;
}

void MuonNumbering::put(std::string_view str, int num) { values_.emplace(str, num); }

#include "FWCore/Utilities/interface/typelookup.h"

TYPELOOKUP_DATA_REG(MuonNumbering);