EMBCell.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef LARREADOUTGEOMETRY_EMBCELL_H
#define LARREADOUTGEOMETRY_EMBCELL_H
 
#include "LArReadoutGeometry/EMBDetDescr.h"
#include "LArReadoutGeometry/EMBHVPathologies.h"
#include "LArHV/EMBHVElectrode.h"
#include "GeoModelKernel/RCBase.h"
#include "CxxUtils/CachedValue.h"
#include <cmath>
#include <vector>
 
class EMBPresamplerHVModule;

 
class EMBCell : public RCBase  
{
 public:
  typedef enum {FRONT,CENTER,BACK} CELLPOS;
 
  public:
      EMBCell (unsigned int side, const EMBDetDescr *embDescriptor, unsigned int eta, unsigned int phi);
 

      unsigned int getPhiIndex () const;

      unsigned int getEndcapIndex () const;

      unsigned int getEtaIndex () const;

      unsigned int getSamplingIndex () const;

      unsigned int getRegionIndex () const;

      const EMBDetDescr * getDescriptor () const;

      double getRLocal (CELLPOS position) const;

      double getZMinLocal (CELLPOS position) const;

      double getZMaxLocal (CELLPOS position) const;

      double getPhiLocalLower () const;

      double getPhiLocalUpper () const;

      double getPhiMinNominal() const;

      double getPhiMaxNominal() const;

      double getEtaMin () const;

      double getEtaMax () const;
 

      double getPhiLocalLower (double r) const;
      double getPhiLocalUpper (double r) const;

      unsigned int getNumElectrodes() const;

      const EMBHVElectrode & getElectrode (unsigned int i) const;

      unsigned int getNumHVPathologies ()const ;

      const EMBHVPathologiesConstLink & getHVPathologies (unsigned int i) const;

      const EMBPresamplerHVModule& getPresamplerHVModule () const;
 

      void addHVPathologies (const EMBHVPathologiesConstLink&);
 
 
  private:
 
      EMBCell(const EMBCell &right);
 
      virtual ~EMBCell();
 
      EMBCell & operator=(const EMBCell &right);
 
      const EMBDetDescr *m_embDetDescr;
 
      struct HVInfo
      {
        // The cell does NOT own the pointers to its electrodes
        std::vector<const EMBHVElectrode*> m_electrode;
        const EMBPresamplerHVModule* m_presamplerModule = nullptr;
      };
      CxxUtils::CachedValue<HVInfo> m_hvinfo;
 
      std::vector<EMBHVPathologiesConstLink> m_hvPathologies;
 
      unsigned int m_clockwork;
 
      friend class ImaginaryFriend;
 
      const HVInfo& getHVInfo() const;
      void initHV (HVInfo& hvinfo) const;
};
 
 
// Class EMBCell 
 
inline EMBCell::EMBCell (unsigned int side, const EMBDetDescr *embDescriptor, unsigned int eta, unsigned int phi)
  :m_embDetDescr(embDescriptor),m_clockwork(phi | (eta<<8) | (side <<17) )
{
}
 
 
 
inline unsigned int EMBCell::getPhiIndex () const
{
  return m_clockwork & 0x000000FF;
}
 
inline unsigned int EMBCell::getEndcapIndex () const
{
  return  (m_clockwork & 0x00020000) >> 17;
}
 
inline unsigned int EMBCell::getEtaIndex () const
{
  return  (m_clockwork & 0x0001FF00) >> 8;
}
 
inline unsigned int EMBCell::getSamplingIndex () const
{
  return getDescriptor()->getSamplingIndex();
}
 
inline unsigned int EMBCell::getRegionIndex () const
{
  return getDescriptor()->getRegionIndex();
}
 
inline const EMBDetDescr * EMBCell::getDescriptor () const
{
  return m_embDetDescr;
}
 
inline double EMBCell::getRLocal (CELLPOS position) const
{
  switch (position) {
  case FRONT:
    return getDescriptor()->getOffset(getEtaIndex())-getDescriptor()->getHalfLength(getEtaIndex());
  case CENTER:
    return getDescriptor()->getOffset(getEtaIndex());
  case BACK:
    return getDescriptor()->getOffset(getEtaIndex())+getDescriptor()->getHalfLength(getEtaIndex());
  default :
    return 0;
  };
}
 
inline double EMBCell::getZMinLocal (CELLPOS position) const
{
  return getRLocal(position)*sinh(getEtaMin());
}
 
inline double EMBCell::getZMaxLocal (CELLPOS position) const
{
  return getRLocal(position)*sinh(getEtaMax());
}
 
inline double EMBCell::getPhiLocalLower () const
{
  return getDescriptor()->getPhiBinning().binLower(getPhiIndex());
}
 
inline double EMBCell::getPhiLocalUpper () const
{
  return getDescriptor()->getPhiBinning().binUpper(getPhiIndex());
}
 
inline double EMBCell::getPhiMinNominal () const
{
  return getEndcapIndex()==0 ? M_PI-getPhiLocalLower() : getPhiLocalLower();
}
 
inline double EMBCell::getPhiMaxNominal () const
{
  return getEndcapIndex()==0 ? M_PI-getPhiLocalUpper() : getPhiLocalUpper();
}
 
inline double EMBCell::getEtaMin () const
{
  return getDescriptor()->getEtaBinning().binLower(getEtaIndex());
}
 
inline double EMBCell::getEtaMax () const
{
  return getDescriptor()->getEtaBinning().binUpper(getEtaIndex());
}
 
 
 
 
#endif