EcalHodoscopeNumberingScheme.cc

Go to the documentation of this file.

// File: EcalHodoscopeNumberingScheme.cc
// Description: Numbering scheme for TB H4 hodoscope detector

#include "Geometry/EcalTestBeam/interface/EcalHodoscopeNumberingScheme.h"
#include "SimDataFormats/EcalTestBeam/interface/HodoscopeDetId.h"

#include <iostream>

EcalHodoscopeNumberingScheme::EcalHodoscopeNumberingScheme() : EcalNumberingScheme() {
  edm::LogInfo("EcalTBGeom") << "Creating EcalHodoscopeNumberingScheme";
}

EcalHodoscopeNumberingScheme::~EcalHodoscopeNumberingScheme() {
  edm::LogInfo("EcalTBGeom") << "Deleting EcalHodoscopeNumberingScheme";
}

uint32_t EcalHodoscopeNumberingScheme::getUnitID(const EcalBaseNumber& baseNumber) const {
  int level = baseNumber.getLevels();
  uint32_t intIndex = 0;
  if (level > 0) {
    // depth index - plans and fibers
    if (baseNumber.getLevelName(0) == "FIBR") {
      uint32_t iFibr;

      if (baseNumber.getCopyNumber(0) > 32)
        iFibr = 2 * (baseNumber.getCopyNumber(0) - 33);
      else
        iFibr = 2 * baseNumber.getCopyNumber(0) - 1;

      const uint32_t iPlane = baseNumber.getCopyNumber(1) - 1;

      LogDebug("EcalTBGeom") << "Fibr/plane " << iFibr << " " << iPlane;

      intIndex = HodoscopeDetId(iPlane, iFibr).rawId();

      LogDebug("EcalTBGeom") << "Index for fiber volume " << baseNumber.getLevelName(0) << " in plane "
                             << baseNumber.getLevelName(1) << " = " << intIndex;
    }
  }
  return intIndex;
}