HECCell.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef LARREADOUTGEOMETRY_HECCELL_H
#define LARREADOUTGEOMETRY_HECCELL_H
 
#include "LArReadoutGeometry/HECDetDescr.h"
#include "GeoModelKernel/CellBinning.h"
#include "GeoModelKernel/RCBase.h"
#include "LArHV/HECHVSubgap.h"
#include "CxxUtils/CachedValue.h"
#include <mutex>

 
class HECCell : public RCBase  
{  
 public:
 
  typedef enum {FRONT,CENTER,BACK} CELLPOS;
  

  HECCell (unsigned int endcap, const HECDetDescr *hecDescriptor, 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 getNumBlocks () const;
      
  unsigned int getBlockIndex (unsigned int block) const;

  double getZMin (unsigned int block) const;
      
  double getZMax (unsigned int block) const;
      
  double getPhiLocalLower () const;
      
  double getPhiLocalUpper () const;
      
  double getPhiMinNominal() const;
  
  double getPhiMaxNominal() const;

  double getEtaMinNominal () const;

  double getEtaMaxNominal () const;

  double getZLocal (CELLPOS position) const;
      
  double getRMinLocalNominal (CELLPOS position) const;
      
  double getRMaxLocalNominal (CELLPOS position) const;
      
  double getRMin (unsigned int block) const;
      
  double getRMax (unsigned int block) const;
      
  const HECDetDescr * getDescriptor () const;
      
  bool occupiesBlock (unsigned int block) const;

  unsigned int getNumSubgaps() const;
  
  const HECHVSubgap& getSubgap (unsigned int i) const;
      
 
 
 private:
    
  HECCell(const HECCell &right);
 
  virtual ~HECCell();
 
  HECCell & operator=(const HECCell &right);
      
  const HECDetDescr* m_hecDetDescr;
      
  unsigned int m_clockwork;

  CxxUtils::CachedValue<std::vector<const HECHVSubgap*> > m_subgaps;

  friend class ImaginaryFriend;
 
  const std::vector<const HECHVSubgap*>& getSubgaps() const;
  void initHV (std::vector<const HECHVSubgap*>& subgaps) const;
};
 
 
 
 
 
inline HECCell::HECCell (unsigned int endcap, const HECDetDescr *hecDescriptor, unsigned int eta, unsigned int phi)
  
  
  
  :m_hecDetDescr(hecDescriptor),m_clockwork(phi | (eta<<6) | (endcap <<10))
  
{
  
  
}
 
 
 
 
inline unsigned int HECCell::getPhiIndex () const
{
  
  return m_clockwork & 0x0000003F;
  
}
 
inline unsigned int HECCell::getEndcapIndex () const
{
  
  return (m_clockwork & 0x00000400) >> 10;
  
}
 
inline unsigned int HECCell::getEtaIndex () const
{
  
  // 4 Bits for eta:
  return (m_clockwork & 0x000003C0) >> 6;
  
}
 
inline unsigned int HECCell::getSamplingIndex () const
{
  
  return getDescriptor()->getSamplingIndex();
  
}
 
inline unsigned int HECCell::getRegionIndex () const
{
  
  return getDescriptor()->getRegionIndex();
  
}
 
inline unsigned int HECCell::getNumBlocks () const
{
  
  return getDescriptor()->getNumBlocks();
  
}
 
inline unsigned int HECCell::getBlockIndex (unsigned int block) const
{
  
  return getDescriptor()->getBlockIndex(block);
  
}
 
inline double HECCell::getZMin (unsigned int block) const
{
  
  return getDescriptor()->getZMin(block);
  
}
 
inline double HECCell::getZMax (unsigned int block) const
{
  
  return getDescriptor()->getZMax(block);
  
}
 
inline double HECCell::getPhiLocalLower () const
{
  
  return getDescriptor()->getPhiBinning().binLower(getPhiIndex());
  
}
 
inline double HECCell::getPhiLocalUpper () const
{
  
  return getDescriptor()->getPhiBinning().binUpper(getPhiIndex());
  
}
 
inline double HECCell::getEtaMinNominal () const
{
  
  return getDescriptor()->getEtaBinning().binLower(getEtaIndex());
  
}
 
inline double HECCell::getEtaMaxNominal () const
{
  
  return getDescriptor()->getEtaBinning().binUpper(getEtaIndex());
  
}
 
inline double HECCell::getZLocal (CELLPOS position) const
{
  
  switch (position) {
  case FRONT:
    return getZMin(0);
  case CENTER:
    return (getZLocal(FRONT) + getZLocal(BACK))/2.0;
  case BACK:
    return getZMax(getNumBlocks()-1);
  default :
    return 0;
  };
  
}
 
inline double HECCell::getRMinLocalNominal (CELLPOS position) const
{
  bool hec2 = getSamplingIndex() > 1;
  if (hec2) {
    return (getDescriptor()->getManager()->getFocalToRef2() + getZLocal(position))/sinh(getEtaMaxNominal());
  }
  else {
    return (getDescriptor()->getManager()->getFocalToRef1() + getZLocal(position))/sinh(getEtaMaxNominal());
  }
 
}
 
inline double HECCell::getRMaxLocalNominal (CELLPOS position) const
{
  bool hec2 = getSamplingIndex() > 1;
  if (hec2) {
    return (getDescriptor()->getManager()->getFocalToRef2() + getZLocal(position))/sinh(getEtaMinNominal());
  }
  else {
    return (getDescriptor()->getManager()->getFocalToRef1() + getZLocal(position))/sinh(getEtaMinNominal());
  }
      
}
 
inline double HECCell::getRMin (unsigned int block) const
{
  
  return getDescriptor()->getRMin(block,getEtaIndex());
  
}
 
inline double HECCell::getRMax (unsigned int block) const
{
  
  return getDescriptor()->getRMax(block,getEtaIndex());
  
}
 
inline const HECDetDescr * HECCell::getDescriptor () const
{
  
  return m_hecDetDescr;
  
}
 
inline bool HECCell::occupiesBlock (unsigned int block) const
{
  
  return getDescriptor()->getRadialSegment(block,getEtaIndex())!=NULL;
  
}
 
 
inline double HECCell::getPhiMinNominal () const
{
  return getEndcapIndex()==0 ? M_PI-getPhiLocalUpper() : getPhiLocalLower();
}
 
inline double HECCell::getPhiMaxNominal () const
{
  return getEndcapIndex()==0 ? M_PI-getPhiLocalLower() : getPhiLocalUpper();
}
 
 
 
#endif