EMECCell.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef LARREADOUTGEOMETRY_EMECCELL_H
#define LARREADOUTGEOMETRY_EMECCELL_H
#include <iostream>
#include "LArReadoutGeometry/EMECDetDescr.h"
#include "LArReadoutGeometry/EMECDetectorManager.h"
#include "GeoModelKernel/CellBinning.h"
#include "GeoModelKernel/RCBase.h"
#include "LArHV/EMECHVElectrode.h"
#include "LArHV/EMECPresamplerHVModule.h"
#include "CxxUtils/CachedValue.h"
 
class EMECCell : public RCBase  
{
  
 public:
 
  typedef enum {FRONT,CENTER,BACK} CELLPOS;
       
  EMECCell (unsigned int endcap, const EMECDetDescr *emecDescriptor, 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;
  
  unsigned int getRadialIndex () const;
  
  const EMECDetDescr * getDescriptor () const;
  
  double getZLocal (CELLPOS position) const;
  
  double getRMinLocal (CELLPOS position) const;
  
  double getRMaxLocal (CELLPOS position) const;
  
  double getPhiLocalLower () const;
  
  double getPhiLocalUpper () const;
 
  double getPhiMinNominal() const;
  
  double getPhiMaxNominal() const;
  
  double getEtaMin () const;
  
  double getEtaMax () const;
 
 
  unsigned int getNumElectrodes() const;
  
  const EMECHVElectrode& getElectrode (unsigned int i) const;
 
 
  const EMECPresamplerHVModule& getPresamplerHVModule () const;
 
 private:
  
  EMECCell(const EMECCell &right);
 
  virtual ~EMECCell();
 
  EMECCell & operator=(const EMECCell &right);
  //    Pointer to the Descriptor.  This class responds to its
  //    interface mostly by forwarding requests to the
  //    descriptor.
 
  const EMECDetDescr *m_emecDetDescr;
 
  struct HVInfo
  {
    // The cell does NOT own the pointers to its electrodes
    std::vector<const EMECHVElectrode*> m_electrode;
    const EMECPresamplerHVModule* m_presamplerModule = nullptr;
  };
  CxxUtils::CachedValue<HVInfo> m_hvinfo;
 
  unsigned int m_clockwork;
 
 
  friend class ImaginaryFriend;
 
  const HVInfo& getHVInfo() const;
  void initHV (HVInfo& hvinfo) const;
};
 
 
 
 
inline EMECCell::EMECCell (unsigned int endcap, const EMECDetDescr *emecDescriptor, unsigned int eta, unsigned int phi)
  // 8 bits are needed for phi index, 6 for eta, 1 for endcap.
 :m_emecDetDescr(emecDescriptor),m_clockwork(phi | (eta<<8) | (endcap <<15) )
{
  
}
 
 
 
 
inline unsigned int EMECCell::getPhiIndex () const
{
  
   return m_clockwork & 0x000000FF;
  
}
 
inline unsigned int EMECCell::getEndcapIndex () const
{
  
  return  (m_clockwork & 0x00008000) >> 15;
  
}
 
inline unsigned int EMECCell::getEtaIndex () const
{
  
  return  (m_clockwork & 0x00007F00) >> 8;
  
}
 
inline unsigned int EMECCell::getSamplingIndex () const
{
  
  return getDescriptor()->getSamplingIndex();
  
}
 
inline unsigned int EMECCell::getRegionIndex () const
{
  
  return getDescriptor()->getRegionIndex();
  
}
 
inline unsigned int EMECCell::getRadialIndex () const
{
  
  return getDescriptor()->getRadialIndex();
  
}
 
inline const EMECDetDescr * EMECCell::getDescriptor () const
{
  
  return m_emecDetDescr;
  
}
 
inline double EMECCell::getZLocal (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 EMECCell::getRMinLocal (CELLPOS position) const
{
  
  double toActive = getDescriptor()->getSamplingIndex()==0 ? getDescriptor()->getManager()->getRefToPresampler() : getDescriptor()->getManager()->getRefToActive();
  return (getDescriptor()->getManager()->getFocalToRef()+toActive+getZLocal(position))/sinh(getEtaMax());
  
}
 
inline double EMECCell::getRMaxLocal (CELLPOS position) const
{
  
  double toActive = getDescriptor()->getSamplingIndex()==0 ? getDescriptor()->getManager()->getRefToPresampler() : getDescriptor()->getManager()->getRefToActive();
  return (getDescriptor()->getManager()->getFocalToRef()+toActive+getZLocal(position))/sinh(getEtaMin());
  
}
 
inline double EMECCell::getPhiLocalLower () const
{
  
  return getDescriptor()->getPhiBinning().binLower(getPhiIndex());
  
}
 
inline double EMECCell::getPhiLocalUpper () const
{
  
  return getDescriptor()->getPhiBinning().binUpper(getPhiIndex());
  
}
 
 
inline double EMECCell::getPhiMinNominal () const
{
  return getEndcapIndex()==0 ? M_PI-getPhiLocalLower() : getPhiLocalLower();
}
 
inline double EMECCell::getPhiMaxNominal () const
{
  return getEndcapIndex()==0 ? M_PI-getPhiLocalUpper() : getPhiLocalUpper();
}
 
inline double EMECCell::getEtaMin () const
{
  
  return getDescriptor()->getEtaBinning().binLower(getEtaIndex()); 
  
}
 
inline double EMECCell::getEtaMax () const
{
  
  return getDescriptor()->getEtaBinning().binUpper(getEtaIndex());
  
}
 
 
#endif