EMBDetDescr.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArReadoutGeometry/EMBDetDescr.h"
#include "LArReadoutGeometry/EMBBasicReadoutNumbers.h"
#include "LArReadoutGeometry/EMBDetectorManager.h"
#include "GeoModelKernel/CellBinning.h"
 
 
namespace {
  
 
struct Boundary {
  double min;
  double max;
  unsigned int nDiv;
  unsigned int firstDiv;
};
 
const Boundary
etaBoundariesBarrel [4][2] = {{{0.000, 1.525,  61,0},{0.000, 0.000,  0,0}},    // 
                {{0.003125, 1.400, 447,1},{1.400, 1.475,  3,0}},    // 
                {{0.000, 1.400,  56,0},{1.400, 1.475,  1,0}},    // 
                {{0.000, 1.350,  27,0},{0.000, 0.000,  0,0}}};   //
 
} // anonymous namespace
 
 
 
 
EMBDetDescr::EMBDetDescr (const EMBDetectorManager *detManager, unsigned int sampling, unsigned int region, const CellBinning &phiBinning)
  :m_manager(detManager),m_samplingIndex(sampling),m_regionIndex(region),m_phiBinning(phiBinning),m_etaBinning(etaBoundariesBarrel[sampling][region].min,
                                                     etaBoundariesBarrel[sampling][region].max,
                                                     etaBoundariesBarrel[sampling][region].nDiv,
                                                     etaBoundariesBarrel[sampling][region].firstDiv)
{
  
  for (unsigned int i=0;i<m_etaBinning.getFirstDivisionNumber()+m_etaBinning.getNumDivisions();i++) {
    if (i<m_etaBinning.getFirstDivisionNumber()) {
      m_offset.push_back(0);
      m_halfLength.push_back(0);
    }
    else {
      double front, back;
      switch (sampling) {
      case 0:
    front = m_manager->getBasicReadoutNumbers()->getPresamplerRadius();
    back  = m_manager->getBasicReadoutNumbers()->getRinAc();
    break;
      case 1:
    front  = m_manager->getBasicReadoutNumbers()->getRinAc();
    back   = region==0 ? m_manager->getBasicReadoutNumbers()->getEMBSamplingSepInnerRMax(i) : m_manager->getBasicReadoutNumbers()->getEMBSamplingSepInnerRMax(447);
    break;
      case 2:
    front   = region==0 ? m_manager->getBasicReadoutNumbers()->getEMBSamplingSepInnerRMax(std::max((unsigned int) 0,8*i+7)) : m_manager->getBasicReadoutNumbers()->getEMBSamplingSepInnerRMax(447);
    back  = m_manager->getBasicReadoutNumbers()->getRMX23(std::min((unsigned int) 52,i));
    break;
      case 3:
    front  = (m_manager->getBasicReadoutNumbers()->getRMX23(std::min((unsigned int)52,2*i))+m_manager->getBasicReadoutNumbers()->getRMX23(std::min((unsigned int)52,2*i+1)))/2.0;
    back  = m_manager->getBasicReadoutNumbers()->getRoutAc();
    break;
      default:
    throw std::range_error("Attempt to construct EMEB DetDescr with illegal sampling index");
    
      }
      m_offset.push_back((front+back)/2.);
      m_halfLength.push_back((back-front)/2.);
    }
  }
}
 
 
EMBDetDescr::~EMBDetDescr()
= default;