DD4hep_DTNumberingScheme.cc

Go to the documentation of this file.

#include "Geometry/MuonNumbering/interface/DD4hep_DTNumberingScheme.h"
#include "Geometry/MuonNumbering/interface/DD4hep_MuonNumbering.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include <cassert>

using namespace cms;

DTNumberingScheme::DTNumberingScheme(const MuonConstants& muonConstants) { initMe(muonConstants); }

void DTNumberingScheme::initMe(const MuonConstants& muonConstants) {
  int levelPart = get("level", muonConstants);
  assert(levelPart != 0);
  theRegionLevel = get("mb_region", muonConstants) / levelPart;
  theWheelLevel = get("mb_wheel", muonConstants) / levelPart;
  theStationLevel = get("mb_station", muonConstants) / levelPart;
  theSuperLayerLevel = get("mb_superlayer", muonConstants) / levelPart;
  theLayerLevel = get("mb_layer", muonConstants) / levelPart;
  theWireLevel = get("mb_wire", muonConstants) / levelPart;
}

int DTNumberingScheme::getDetId(const MuonBaseNumber& num) const {
  int wire_id(0);
  int layer_id(0);
  int superlayer_id(0);
  int sector_id(0);
  int station_id(0);
  int wheel_id(0);

  //decode significant barrel levels
  decode(num, wire_id, layer_id, superlayer_id, sector_id, station_id, wheel_id);

  DTWireId id(wheel_id, station_id, sector_id, superlayer_id, layer_id, wire_id);

  return id.rawId();
}

void DTNumberingScheme::decode(const MuonBaseNumber& num,
                               int& wire_id,
                               int& layer_id,
                               int& superlayer_id,
                               int& sector_id,
                               int& station_id,
                               int& wheel_id) const {
  for (int level = 1; level <= num.getLevels(); ++level) {
    //decode
    if (level == theWheelLevel) {
      const int copyno = num.getBaseNo(level);
      wheel_id = copyno - 2;

    } else if (level == theStationLevel) {
      const int station_tag = num.getSuperNo(level);
      const int copyno = num.getBaseNo(level);
      station_id = station_tag;
      sector_id = copyno + 1;
    } else if (level == theSuperLayerLevel) {
      const int copyno = num.getBaseNo(level);
      superlayer_id = copyno + 1;
    } else if (level == theLayerLevel) {
      const int copyno = num.getBaseNo(level);
      layer_id = copyno + 1;
    } else if (level == theWireLevel) {
      const int copyno = num.getBaseNo(level);
      wire_id = copyno + 1;
    }
  }
}

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