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;
}