EMBHVManager.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArHV/EMBHVManager.h"
#include "LArHV/EMBHVModule.h"
#include "LArHV/EMBHVElectrode.h"
#include "GeoModelKernel/CellBinning.h"
#include <cmath>
#include <vector>
#include "EMBHVPayload.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"
#include "LArCabling/LArHVCablingSimTool.h"
 
#ifndef SIMULATIONBASE
#ifndef GENERATIONBASE
#include "LArRecConditions/LArHVIdMapping.h"
#endif
#endif
 
#include "Identifier/HWIdentifier.h"
#include <mutex>
#include <atomic>
 
 
namespace {
 
 
struct SimIdFunc
{
  SimIdFunc();
  std::vector<HWIdentifier> operator()(HWIdentifier id) const
  {
    return m_cablingTool->getLArElectrodeIDvec (id);
  }
  LArHVCablingSimTool* m_cablingTool;
};
 
 
SimIdFunc::SimIdFunc()
{
  ToolHandle<LArHVCablingSimTool> tool ("LArHVCablingSimTool");
  if (!tool.retrieve().isSuccess()) {
    std::abort();
  }
  m_cablingTool = tool.get();
}
 
 
} // Anonymous namespace
 
 
class EMBHVManager::Clockwork {
public:
  explicit Clockwork(const EMBHVManager* manager) {
    for(int iSide=0; iSide<2; ++iSide) {
      for(int iEta=0; iEta<8; ++iEta) {
    for(int iPhi=0; iPhi<16; ++iPhi) {
      for(int iSector=0; iSector<2; ++iSector) {
        moduleArray[iSide][iEta][iPhi][iSector] = std::make_unique<EMBHVModule>(manager,iSide,iEta,iPhi,iSector);
      }
    }
      }
    }
 
    ServiceHandle<StoreGateSvc> detStore ("DetectorStore", "HECHVManager");
    if (StatusCode::SUCCESS!=detStore->retrieve(elecId, "LArElectrodeID")) {
      throw std::runtime_error("EMBHVManager failed to retrieve LArElectrodeID");
    }
 
    if (StatusCode::SUCCESS!=detStore->retrieve(hvId,"LArHVLineID")) {
      throw std::runtime_error("EMBHVManager failed to retrieve LArHVLineID");
    }
  }
  Clockwork(const Clockwork&) = delete;
  ~Clockwork() = default;
  EMBHVDescriptor       descriptor{CellBinning(0.0, 1.4, 7, 1),CellBinning(0.0, 2*M_PI, 16)};
  std::unique_ptr<const EMBHVModule>    moduleArray[2][8][16][2];
  const LArElectrodeID* elecId = nullptr;
  const LArHVLineID* hvId = nullptr;
};
 
 
class EMBHVManager::EMBHVData::Payload
{
public:
  std::vector<EMBHVPayload> m_payloadArray;
};
 
 
EMBHVManager::EMBHVData::EMBHVData() = default;
 
 
EMBHVManager::EMBHVData::EMBHVData (std::unique_ptr<Payload> payload)
  : m_payload (std::move (payload))
{
}
 
 
EMBHVManager::EMBHVData&
EMBHVManager::EMBHVData::operator= (EMBHVData&& other) noexcept {
  if (this != &other) {
    m_payload = std::move (other.m_payload);
  }
  return *this;
}
 
 
EMBHVManager::EMBHVData::~EMBHVData()
= default;
 
 
bool EMBHVManager::EMBHVData::hvOn (const EMBHVElectrode& electrode, const int& iGap) const
{
  return voltage (electrode, iGap) > INVALID;
}
 
 
double EMBHVManager::EMBHVData::voltage (const EMBHVElectrode& electrode, const int& iGap) const
{
  return m_payload->m_payloadArray[index(electrode)].voltage[iGap];
}
 
 
double EMBHVManager::EMBHVData::current (const EMBHVElectrode& electrode, const int& iGap) const
{
  return m_payload->m_payloadArray[index(electrode)].current[iGap];
}
 
 
int  EMBHVManager::EMBHVData::hvLineNo  (const EMBHVElectrode& electrode, const int& iGap) const
{
  return m_payload->m_payloadArray[index(electrode)].hvLineNo[iGap];
}
 
 
int  EMBHVManager::EMBHVData::index  (const EMBHVElectrode& electrode)
{
  unsigned int electrodeIndex    = electrode.getElectrodeIndex();
  const EMBHVModule& module      = electrode.getModule();
  unsigned int etaIndex          = module.getEtaIndex();
  unsigned int phiIndex          = module.getPhiIndex();
  unsigned int sectorIndex       = module.getSectorIndex();
  unsigned int sideIndex         = module.getSideIndex();
  unsigned int index             = 8192*sideIndex+1024*etaIndex+64*phiIndex+32*sectorIndex+electrodeIndex;
  return index;
}
 
 
EMBHVManager::EMBHVManager()
  : m_c (std::make_unique<Clockwork> (this))
{
}
 
EMBHVManager::~EMBHVManager()
= default;
 
const EMBHVDescriptor& EMBHVManager::getDescriptor() const
{
  return m_c->descriptor;
}
 
unsigned int EMBHVManager::beginPhiIndex() const
{
  return m_c->descriptor.getPhiBinning().getFirstDivisionNumber();
}
 
unsigned int EMBHVManager::endPhiIndex() const
{
  return m_c->descriptor.getPhiBinning().getFirstDivisionNumber() + m_c->descriptor.getPhiBinning().getNumDivisions();
}
 
unsigned int EMBHVManager::beginEtaIndex() const
{
  return m_c->descriptor.getEtaBinning().getFirstDivisionNumber();
}
 
unsigned int EMBHVManager::endEtaIndex() const
{
  return m_c->descriptor.getEtaBinning().getFirstDivisionNumber() + m_c->descriptor.getEtaBinning().getNumDivisions();
}
 
const EMBHVModule& EMBHVManager::getHVModule(unsigned int iSide, unsigned int iEta,unsigned int iPhi, unsigned int iSector) const
{
  return *(m_c->moduleArray[iSide][iEta][iPhi][iSector]);
}
 
unsigned int EMBHVManager::beginSectorIndex()
{
  return 0;
}
 
unsigned int EMBHVManager::endSectorIndex()
{
  return 2;
}
 
unsigned int EMBHVManager::beginSideIndex()
{
  return 0;
}
 
unsigned int EMBHVManager::endSideIndex()
{
  return 2;
}
 
EMBHVManager::EMBHVData
EMBHVManager::getData (const idfunc_t& idfunc,
                       const std::vector<const CondAttrListCollection*>& attrLists) const
{
  auto payload = std::make_unique<EMBHVData::Payload>();
  payload->m_payloadArray.reserve(2*8*16*2*32);
 
  for (int i=0;i<16384;i++) {
    payload->m_payloadArray[i].voltage[0] = EMBHVData::INVALID;
    payload->m_payloadArray[i].voltage[1] = EMBHVData::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==0) {
 
// 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);
// phi index
//  offline 0 to 2pi in 2pi/16 bins
// this is module in the electrode numbering: on the A side 0 to 15, 0 is halfway around phi=0 (FT-1 (hv_phi=1 is a lower phi)
//   offline phi      0               pi                    2pi
//   Module           M0 M1 M1                          M15 M0
//   FT               0 -1  0                           0  -1
//   hv_phi           0  1  0                           0   1
//   phiIndex         0  0  1                          15   15
// sector Index       0  1  0                           0   1
          unsigned int phiIndex;
          unsigned int sectorIndex;
          if (sideIndex==1) {
            phiIndex=m_c->elecId->module(elecHWID);
            sectorIndex=m_c->elecId->hv_phi(elecHWID);
          }
// module numbering on the C side 0 around phi=pi, running backwards
//   offline phi     0               pi                      2pi
//   Module          P8 P7 P7      P0 P0                 P9 P8
//   FT              -1 0 -1       0  -1                -1  0
//   hv_phi           1 0  1       0   1                 1   0
//  phiIndex          0 0  1       7   8                15  15
// sectorIndex        0 1  0       1   0                 0  1
          else {
            int imodule=m_c->elecId->module(elecHWID);
            if (imodule<9) phiIndex = 8 - imodule;
            else           phiIndex = 24 - imodule;
            sectorIndex = 1-m_c->elecId->hv_phi(elecHWID);
          }
 
          if (sectorIndex==1) {
            if (phiIndex>0) phiIndex = phiIndex - 1;
            else phiIndex=15;
          }
 
          unsigned int electrodeIndex=m_c->elecId->electrode(elecHWID);
          if (sideIndex==0) {
            if (m_c->elecId->hv_phi(elecHWID)==1) electrodeIndex=31-electrodeIndex;  // FT-1 change 0->31 to 31->0
            else                             electrodeIndex=63-electrodeIndex;  // FT 0 change 32->63 to 31-0
          }
          else {
            if (m_c->elecId->hv_phi(elecHWID)==0) electrodeIndex=electrodeIndex-32;  // FT 0 change 31-63 to 0-31
          }
 
          unsigned int index             = 8192*sideIndex+1024*etaIndex+64*phiIndex+32*sectorIndex+electrodeIndex;
 
          unsigned int gapIndex=m_c->elecId->gap(elecHWID);
          if (sideIndex==0) gapIndex=1-gapIndex;
 
          float voltage = EMBHVData::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)};
}
 
 
EMBHVManager::EMBHVData
EMBHVManager::getDataSim() const
{
  std::vector<const CondAttrListCollection*> attrLists;
  ServiceHandle<StoreGateSvc> detStore ("DetectorStore", "EMBHVManager");
  const CondAttrListCollection* atrlistcol = nullptr;
  // Not a typo --- this folder has a lower-case l in the database...
  if (detStore->retrieve(atrlistcol, "/LAR/DCS/HV/BARREl/I16").isSuccess()) {
    attrLists.push_back (atrlistcol);
  }
  if (detStore->retrieve(atrlistcol, "/LAR/DCS/HV/BARREL/I8").isSuccess()) {
    attrLists.push_back (atrlistcol);
  }
  return getData (SimIdFunc(), attrLists);
}
 
 
#ifndef SIMULATIONBASE
#ifndef GENERATIONBASE
EMBHVManager::EMBHVData
EMBHVManager::getData (const LArHVIdMapping& hvIdMapping,
                       const std::vector<const CondAttrListCollection*>& attrLists) const
{
  auto idfunc = [&] (HWIdentifier id) { return hvIdMapping.getLArElectrodeIDvec(id); };
  return getData (idfunc, attrLists);
}
 
 
int EMBHVManager::hvLineNo(const EMBHVElectrode& electrode
               , int gap
               , const LArHVIdMapping* hvIdMapping
                           , HWIdentifier *hvlId) const
{
  const EMBHVModule& module      = electrode.getModule();
  int etaIndex          = module.getEtaIndex();
  int phiIndex          = module.getPhiIndex();
  int sectorIndex       = module.getSectorIndex();
  int sideIndex         = module.getSideIndex();
  int electrodeIndex    = electrode.getElectrodeIndex();
 
  // ________________________ Construct ElectrodeID ________________________________
  int id_detector = 0;
  int id_zside = 1 - sideIndex;
  int id_hv_phi = (sideIndex==1 ? sectorIndex : 1-sectorIndex);
  int id_hv_eta = etaIndex;
  int id_gap = (sideIndex==0 ? 1-gap : gap);
  int tmpPhi = phiIndex;
  if(sectorIndex==1) {
    if(tmpPhi==15) {
      tmpPhi=0;
    }
    else {
      tmpPhi += 1;
    }
  }
 
  int id_module;
  if(sideIndex==1) {
    id_module = tmpPhi;
  }
  else {
    if(tmpPhi<9) {
      id_module = 8-tmpPhi;
    }
    else {
      id_module = 24-tmpPhi;
    }
  }
 
  int id_electrode;
  if(sideIndex==0) {
    if(id_hv_phi==1) {
      id_electrode = 31-electrodeIndex;
    }
    else {
      id_electrode = 63-electrodeIndex;
    }
  }
  else {
    if(id_hv_phi==0) {
      id_electrode = electrodeIndex+32;
    }
    else {
      id_electrode = electrodeIndex;
    }
  }
 
  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