EMBPresamplerHVManager.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArHV/EMBPresamplerHVManager.h"
#include "LArHV/EMBPresamplerHVDescriptor.h"
#include "GeoModelKernel/CellBinning.h"
#include "GeoModelKernel/CellPartitioning.h"
#include <cmath>
#include <vector>
#include <fstream>
#include "EMBPresamplerHVPayload.h"
 
#include "StoreGate/StoreGateSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/IToolSvc.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ServiceHandle.h"
 
#include "AthenaPoolUtilities/AthenaAttributeList.h"
#include "AthenaPoolUtilities/CondAttrListCollection.h"
 
#include "LArIdentifier/LArElectrodeID.h"
#include "LArIdentifier/LArHVLineID.h"
 
#ifndef SIMULATIONBASE
#ifndef GENERATIONBASE
#include "LArRecConditions/LArHVIdMapping.h"
#endif
#endif
 
#include "Identifier/HWIdentifier.h"
 
#include <atomic>
 
 
class EMBPresamplerHVManager::Clockwork {
public:
  explicit Clockwork(const EMBPresamplerHVManager* manager) {
    CellPartitioning etaPartitioning;
    for (unsigned int i= 0; i<4; i++)  etaPartitioning.addValue(i*0.4);
    etaPartitioning.addValue(1.525);
    descriptor = std::make_unique<EMBPresamplerHVDescriptor>(etaPartitioning,CellBinning(0.0, 2*M_PI, 32));
 
    for(int iSide=0; iSide<2; ++iSide) {
      for(int iEta=0; iEta<4; ++iEta) {
    for(int iPhi=0; iPhi<32; ++iPhi) {
      moduleArray[iSide][iEta][iPhi] = std::make_unique<EMBPresamplerHVModule>(manager, iSide, iEta,iPhi);
    }
      }
    }
 
    ServiceHandle<StoreGateSvc> detStore ("DetectorStore", "HECHVManager");
    if (StatusCode::SUCCESS!=detStore->retrieve(elecId, "LArElectrodeID")) {
      throw std::runtime_error("EMBPresamplerHVManager failed to retrieve LArElectrodeID");
    }
 
    if (StatusCode::SUCCESS!=detStore->retrieve(hvId,"LArHVLineID")) {
      throw std::runtime_error("EMBPresamplerHVManager failed to retrieve LArHVLineID");
    }
  }
  ~Clockwork() = default;
  std::unique_ptr<EMBPresamplerHVDescriptor>     descriptor;
  std::unique_ptr<const EMBPresamplerHVModule>   moduleArray[2][4][32];
  const LArElectrodeID* elecId = nullptr;
  const LArHVLineID* hvId = nullptr;
};
 
 
class EMBPresamplerHVManager::EMBPresamplerHVData::Payload
{
public:
  std::vector<EMBPresamplerHVPayload> m_payloadArray;
};
 
 
EMBPresamplerHVManager::EMBPresamplerHVData::EMBPresamplerHVData() = default;
 
 
EMBPresamplerHVManager::EMBPresamplerHVData::EMBPresamplerHVData
  (std::unique_ptr<Payload> payload)
  : m_payload (std::move (payload))
{
}
 
 
EMBPresamplerHVManager::EMBPresamplerHVData&
EMBPresamplerHVManager::EMBPresamplerHVData::operator= (EMBPresamplerHVData&& other) noexcept {
  if (this != &other) {
    m_payload = std::move (other.m_payload);
  }
  return *this;
}
 
 
EMBPresamplerHVManager::EMBPresamplerHVData::~EMBPresamplerHVData()
= default;
 
 
bool EMBPresamplerHVManager::EMBPresamplerHVData::hvOn
  (const EMBPresamplerHVModule& module, const int& iGap) const
{
  return voltage (module, iGap) > INVALID;
}
 
 
double EMBPresamplerHVManager::EMBPresamplerHVData::voltage
  (const EMBPresamplerHVModule& module, const int& iGap) const
{
  return m_payload->m_payloadArray[index(module)].voltage[iGap];
}
 
 
double EMBPresamplerHVManager::EMBPresamplerHVData::current
  (const EMBPresamplerHVModule& module, const int& iGap) const
{
  return m_payload->m_payloadArray[index(module)].current[iGap];
}
 
 
int EMBPresamplerHVManager::EMBPresamplerHVData::hvLineNo
  (const EMBPresamplerHVModule& module, const int& iGap) const
{
  return m_payload->m_payloadArray[index(module)].hvLineNo[iGap];
}
 
 
int EMBPresamplerHVManager::EMBPresamplerHVData::index
  (const EMBPresamplerHVModule& module)
{
  unsigned int sideIndex         = module.getSideIndex();
  unsigned int phiIndex          = module.getPhiIndex();
  unsigned int etaIndex          = module.getEtaIndex();
  unsigned int index             =  128*sideIndex+32*etaIndex+phiIndex;
  return index;
}
 
 
EMBPresamplerHVManager::EMBPresamplerHVManager()
  : m_c (std::make_unique<Clockwork> (this))
{
}
 
EMBPresamplerHVManager::~EMBPresamplerHVManager()
= default;
 
const EMBPresamplerHVDescriptor* EMBPresamplerHVManager::getDescriptor() const
{
  return m_c->descriptor.get();
}
 
unsigned int EMBPresamplerHVManager::beginPhiIndex() const
{
  return m_c->descriptor->getPhiBinning().getFirstDivisionNumber();
}
 
unsigned int EMBPresamplerHVManager::endPhiIndex() const
{
  return m_c->descriptor->getPhiBinning().getFirstDivisionNumber() + m_c->descriptor->getPhiBinning().getNumDivisions();
}
 
unsigned int EMBPresamplerHVManager::beginEtaIndex() const
{
  return m_c->descriptor->getEtaPartitioning().getFirstDivisionNumber();
}
 
unsigned int EMBPresamplerHVManager::endEtaIndex() const
{
  return m_c->descriptor->getEtaPartitioning().getFirstDivisionNumber() + m_c->descriptor->getEtaPartitioning().getNumDivisions();
}
 
const EMBPresamplerHVModule& EMBPresamplerHVManager::getHVModule(unsigned int iSide, unsigned int iEta,unsigned int iPhi) const
{
  return *(m_c->moduleArray[iSide][iEta][iPhi]);
}
 
unsigned int EMBPresamplerHVManager::beginSideIndex()
{
  return 0;
}
 
unsigned int EMBPresamplerHVManager::endSideIndex()
{
  return 2;
}
 
 
EMBPresamplerHVManager::EMBPresamplerHVData
EMBPresamplerHVManager::getData (const idfunc_t& idfunc,
                                 const std::vector<const CondAttrListCollection*>& attrLists) const
{
  auto payload = std::make_unique<EMBPresamplerHVData::Payload>();
  payload->m_payloadArray.reserve(2*4*32);
 
  for (int i=0;i<256;i++) {
    payload->m_payloadArray[i].voltage[0] = EMBPresamplerHVData::INVALID;
    payload->m_payloadArray[i].voltage[1] = EMBPresamplerHVData::INVALID;
  }
 
  for (const CondAttrListCollection* atrlistcol : attrLists) {
 
    for (CondAttrListCollection::const_iterator citr=atrlistcol->begin(); citr!=atrlistcol->end();++citr) {
 
      // Construct HWIdentifier
      // 1. decode COOL Channel ID
      unsigned int chanID = (*citr).first;
      int cannode = chanID/1000;
      int line = chanID%1000;
 
      // 2. Construct the identifier
      HWIdentifier id = m_c->hvId->HVLineId(1,1,cannode,line);
 
      std::vector<HWIdentifier> electrodeIdVec = idfunc(id);
 
 
      for(size_t i=0;i<electrodeIdVec.size();i++)
      {
        HWIdentifier& elecHWID = electrodeIdVec[i];
        int detector = m_c->elecId->detector(elecHWID);
        if (detector==1) {
 
// side  in standard offline 0 for z<0 (C) 1 for z>0 (A)
//  in electrode numbering, this is the opposite (0 for A and 1 for C)
          unsigned int sideIndex=1-m_c->elecId->zside(elecHWID);
// eta index, no trouble
          unsigned int etaIndex=m_c->elecId->hv_eta(elecHWID)-1;
// phi index
//  offline 0 to 2pi in 2pi/32 bins  (2 presampler cells per phi_HV)
//  module from elecID : 0 to 31
//    phi        0                                     pi                                            2pi
//               M0 M1  M2  M3  M4  M5   M6    M7    M8    M9    M10   M11   M12   M13   M14   M15   M0
//     FT        0 -1 0                                                                        -1 0  -1
//   Module      0  1 2 3 4 5 6 7 9 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 39 30 31
//  offline phi  0  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
          unsigned int phiIndex;
          if (sideIndex==1) {
            phiIndex=m_c->elecId->module(elecHWID);
          }
// module numbering on the C side 0 around phi=pi, running backwards
//   phi           0                                     pi                                            2pi
//                 P8 P7  P6  P5  P4  P3   P2    P1    P0    P15   P14   P13   P12   P11   P10   P9    P8
//     FT          -1 0                                0  -1                                      0 -1  0
// Module          15 14        9 8 7 6 5  4  3  2  1   0 31 30 29 38 27 26 25 24 23 22 21 20 19 18 17 16
//  offline phi    0  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
          else {
            int imodule=m_c->elecId->module(elecHWID);
            if (imodule<16) phiIndex = 15 - imodule;
            else           phiIndex =  47 - imodule;
          }
 
          unsigned int index             =  128*sideIndex+32*etaIndex+phiIndex;
 
          unsigned int gapIndex=m_c->elecId->gap(elecHWID);
          if (sideIndex==0) gapIndex=1-gapIndex;
 
 
          float voltage = EMBPresamplerHVData::INVALID;
          if (!((*citr).second)["R_VMEAS"].isNull()) voltage = ((*citr).second)["R_VMEAS"].data<float>();
          float current = 0.;
          if (!((*citr).second)["R_IMEAS"].isNull()) current = ((*citr).second)["R_IMEAS"].data<float>();
 
          payload->m_payloadArray[index].voltage[gapIndex]=voltage;
          payload->m_payloadArray[index].current[gapIndex]=current;
          payload->m_payloadArray[index].hvLineNo[gapIndex]=chanID;
        }
      }
    }
  }
 
  return {std::move (payload)};
}
 
#ifndef SIMULATIONBASE
#ifndef GENERATIONBASE
EMBPresamplerHVManager::EMBPresamplerHVData
EMBPresamplerHVManager::getData (const LArHVIdMapping& hvIdMapping,
                                 const std::vector<const CondAttrListCollection*>& attrLists) const
{
  auto idfunc = [&] (HWIdentifier id) { return hvIdMapping.getLArElectrodeIDvec(id); };
  return getData (idfunc, attrLists);
}
 
 
int EMBPresamplerHVManager::hvLineNo(const EMBPresamplerHVModule& module
                     , int gap
                     , const LArHVIdMapping* hvIdMapping
                                     , HWIdentifier *hvlId ) const {
 
  int sideIndex         = module.getSideIndex();
  int phiIndex          = module.getPhiIndex();
  int etaIndex          = module.getEtaIndex();
 
  // ________________________ Construct ElectrodeID ________________________________
  int id_detector = 1;
  int id_zside = 1 - sideIndex;
  int id_module{0};
  if(sideIndex==1){
    id_module = phiIndex;
  }
  else {
    if(phiIndex<16) {
      id_module = 15 - phiIndex;
    }
    else {
      id_module = 47 - phiIndex;
    }
  }
  int id_hv_phi{0};
  int id_hv_eta = etaIndex + 1;
  int id_gap = sideIndex==0 ? 1-gap : gap;
  int id_electrode{0};
 
  HWIdentifier elecHWID = m_c->elecId->ElectrodeId(id_detector
                           ,id_zside
                           ,id_module
                           ,id_hv_phi
                           ,id_hv_eta
                           ,id_gap
                           ,id_electrode);
  // ________________________  ________________________________
 
  // Get LArHVLineID corresponding to a given LArElectrodeId
  HWIdentifier id = hvIdMapping->getLArHVLineID(elecHWID);
 
  if(hvlId) *hvlId = id;
 
  // Extract HV Line No
  return m_c->hvId->can_node(id)*1000 + m_c->hvId->hv_line(id);
}
#endif
#endif