ETLNumberingScheme.cc

Go to the documentation of this file.

#include "Geometry/MTDCommonData/interface/ETLNumberingScheme.h"
#include "DataFormats/ForwardDetId/interface/ETLDetId.h"

#include <iostream>

//#define EDM_ML_DEBUG

ETLNumberingScheme::ETLNumberingScheme() : MTDNumberingScheme() {
#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "Creating ETLNumberingScheme";
#endif
}

ETLNumberingScheme::~ETLNumberingScheme() {
#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "Deleting ETLNumberingScheme";
#endif
}

uint32_t ETLNumberingScheme::getUnitID(const MTDBaseNumber& baseNumber) const {
  const uint32_t nLevels(baseNumber.getLevels());

#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "ETLNumberingScheme geometry levels = " << nLevels;
#endif

  if (11 > nLevels) {
    edm::LogWarning("MTDGeom") << "ETLNumberingScheme::getUnitID(): "
                               << "Not enough levels found in MTDBaseNumber ( " << nLevels << ") Returning 0";
    return 0;
  }

  const uint32_t modCopy(baseNumber.getCopyNumber(2));

  const std::string& ringName(baseNumber.getLevelName(3));  // name of ring volume
  const int modtyp(0);
  std::string baseName = ringName.substr(ringName.find(":") + 1);
  const int ringCopy(::atoi(baseName.c_str() + 4));

  // Side choice: up to scenario D38 is given by level 7 (HGCal v9)
  int nSide(7);
  const std::string& sideName(baseNumber.getLevelName(nSide));
  // Side choice: from scenario D41 is given by level 8 (HGCal v10)
  if (sideName.find("CALOECTSFront") != std::string::npos) {
    nSide = 8;
  }
  const uint32_t sideCopy(baseNumber.getCopyNumber(nSide));
  const uint32_t zside(sideCopy == 1 ? 1 : 0);

#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << baseNumber.getLevelName(0) << ", " << baseNumber.getLevelName(1) << ", "
                          << baseNumber.getLevelName(2) << ", " << baseNumber.getLevelName(3) << ", "
                          << baseNumber.getLevelName(4) << ", " << baseNumber.getLevelName(5) << ", "
                          << baseNumber.getLevelName(6) << ", " << baseNumber.getLevelName(7) << ", "
                          << baseNumber.getLevelName(8) << ", " << baseNumber.getLevelName(9) << ", "
                          << baseNumber.getLevelName(10) << ", " << baseNumber.getLevelName(11);
#endif

  // error checking

  if (modtyp != 0) {
    edm::LogWarning("MTDGeom") << "ETLNumberingScheme::getUnitID(): "
                               << "****************** Bad module name = " << modtyp
                               << ", Volume Name = " << baseNumber.getLevelName(4);
    return 0;
  }

  if (1 > modCopy || 176 < modCopy) {
    edm::LogWarning("MTDGeom") << "ETLNumberingScheme::getUnitID(): "
                               << "****************** Bad module copy = " << modCopy
                               << ", Volume Number = " << baseNumber.getCopyNumber(4);
    return 0;
  }

  if (1 > ringCopy || 11 < ringCopy) {
    edm::LogWarning("MTDGeom") << "ETLNumberingScheme::getUnitID(): "
                               << "****************** Bad ring copy = " << ringCopy
                               << ", Volume Number = " << baseNumber.getCopyNumber(3);
    return 0;
  }

  // all inputs are fine. Go ahead and decode

  ETLDetId thisETLdetid(zside, ringCopy, modCopy, modtyp);
  const int32_t intindex = thisETLdetid.rawId();

#ifdef EDM_ML_DEBUG
  edm::LogInfo("MTDGeom") << "ETL Numbering scheme: "
                          << " ring = " << ringCopy << " zside = " << zside << " module = " << modCopy
                          << " modtyp = " << modtyp << " Raw Id = " << intindex << thisETLdetid;
#endif

  return intindex;
}