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 constant:" << superPart;
  }
  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);