MuonDDDNumbering.cc

Go to the documentation of this file.

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

//#define LOCAL_DEBUG

MuonDDDNumbering::MuonDDDNumbering(const MuonDDDConstants &muonConstants) {
  //  MuonDDDConstants muonConstants;
  theLevelPart = muonConstants.getValue("level");
  theSuperPart = muonConstants.getValue("super");
  theBasePart = muonConstants.getValue("base");
  theStartCopyNo = muonConstants.getValue("xml_starts_with_copyno");

  // some consistency checks

  if (theBasePart != 1) {
    edm::LogWarning("Geometry") << "MuonDDDNumbering finds unusual base constant:" << theBasePart;
  }
  if (theSuperPart < 100) {
    edm::LogWarning("Geometry") << "MuonDDDNumbering finds unusual super constant:" << theSuperPart;
  }
  if (theLevelPart < 10 * theSuperPart) {
    edm::LogWarning("Geometry") << "MuonDDDNumbering finds unusual level constant:" << theLevelPart;
  }
  if ((theStartCopyNo != 0) && (theStartCopyNo != 1)) {
    edm::LogWarning("Geometry") << "MuonDDDNumbering finds unusual start value for copy numbers:" << theStartCopyNo;
  }

#ifdef LOCAL_DEBUG
  edm::LogVerbatim("Geometry") << "MuonDDDNumbering configured with"
                               << " Level = " << theLevelPart << " Super = " << theSuperPart
                               << " Base = " << theBasePart << " StartCopyNo = " << theStartCopyNo;
#endif
}

MuonBaseNumber MuonDDDNumbering::geoHistoryToBaseNumber(const DDGeoHistory &history) {
  MuonBaseNumber num;

#ifdef LOCAL_DEBUG
  edm::LogVerbatim("Geometry") << "MuonDDDNumbering create MuonBaseNumber for " << history;
#endif

  //loop over all parents and check
  DDGeoHistory::const_iterator cur = history.begin();
  DDGeoHistory::const_iterator end = history.end();
  while (cur != end) {
    const DDLogicalPart &ddlp = cur->logicalPart();
    const int tag = getInt("CopyNoTag", ddlp) / theLevelPart;
    if (tag > 0) {
      const int offset = getInt("CopyNoOffset", ddlp);
      const int copyno = (cur->copyno()) + offset % theSuperPart;
      const int super = offset / theSuperPart;
      num.addBase(tag, super, copyno - theStartCopyNo);
    }
    cur++;
  }

#ifdef LOCAL_DEBUG
  edm::LogVerbatim("Geometry") << "MuonDDDNumbering::" << num.getLevels();
  for (int i = 1; i <= num.getLevels(); i++) {
    edm::LogVerbatim("Geometry") << "[" << i << "] " << num.getSuperNo(i) << " " << num.getBaseNo(i);
  }
#endif

  return num;
}

int MuonDDDNumbering::getInt(const std::string &s, const DDLogicalPart &part) {
  DDValue val(s);
  std::vector<const DDsvalues_type *> result = part.specifics();
  std::vector<const DDsvalues_type *>::iterator it = result.begin();
  bool foundIt = false;
  for (; it != result.end(); ++it) {
    foundIt = DDfetch(*it, val);
    if (foundIt)
      break;
  }
  if (foundIt) {
    std::vector<double> temp = val.doubles();
    if (temp.size() != 1) {
      edm::LogError("Geometry") << "MuonDDDNumbering:: ERROR: I need only 1 " << s << " in DDLogicalPart "
                                << part.name();
      abort();
    }
    return int(temp[0]);
  } else
    return 0;
}