CaloDetectorElements.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/

 
#ifndef CALODETDESCR_CALODETECTORELEMENTS_H
#define CALODETDESCR_CALODETECTORELEMENTS_H
 
#include "CaloDetDescr/CaloDetDescrElement.h"
 
#include "LArReadoutGeometry/EMBCell.h"
#include "LArReadoutGeometry/EMBDetectorRegion.h"
#include "LArReadoutGeometry/EMECCell.h"
#include "LArReadoutGeometry/EMECDetectorRegion.h"
#include "LArReadoutGeometry/HECCell.h"
#include "LArReadoutGeometry/HECDetectorRegion.h"
#include "LArReadoutGeometry/FCALTile.h"
#include "LArReadoutGeometry/FCALModule.h"
 
// MBTS
#include "Identifier/Identifier.h"
 
class CaloDetDescriptor;
class GeoAlignmentStore;

 
typedef struct 
{
  float dx;
  float dy;
  float dz;
} CaloElementPositionShift;

class EMBDetectorElement : public CaloDetDescrElement
{
 public:
  EMBDetectorElement(const IdentifierHash subcaloHash
             , const IdentifierHash onl1
             , const IdentifierHash onl2
             , const CaloDetDescriptor* descriptor
             , EMBCellConstLink& embCell
             , const EMBDetectorRegion* embRegion
             , bool isTestBeam
             , const GeoAlignmentStore* geoAlignStore = nullptr
                     , const CaloElementPositionShift* posShift = nullptr);
  
  inline EMBCellConstLink getEMBCell() const {return m_cell;}
  
  virtual int getLayer() const override;
 
 private:
  EMBCellConstLink m_cell;
  const EMBDetectorRegion* m_region;
  
  EMBDetectorElement() = delete;
  EMBDetectorElement & operator=(const EMBDetectorElement &right) = delete;
  
  void init_description(const GeoAlignmentStore* geoAlignStore
            , const CaloElementPositionShift* posShift);
};

class EMECDetectorElement : public CaloDetDescrElement
{
 public:
  EMECDetectorElement(const IdentifierHash subcaloHash
              , const IdentifierHash onl1
              , const IdentifierHash onl2
              , const CaloDetDescriptor* descriptor
              , EMECCellConstLink& emecCell
              , const EMECDetectorRegion* emecRegion
              , bool isTestBeam
              , const GeoAlignmentStore* geoAlignStore = nullptr
                      , const CaloElementPositionShift* posShift = nullptr);

  inline EMECCellConstLink getEMECCell() const {return m_cell;}

  virtual int getLayer() const override;
 
 private:
  EMECCellConstLink m_cell;
  const EMECDetectorRegion* m_region;
 
  EMECDetectorElement() = delete;
  EMECDetectorElement & operator=(const EMECDetectorElement &right) = delete;

  void init_description(bool isTestBeam
            , const GeoAlignmentStore* geoAlignStore
            , const CaloElementPositionShift* posShift);

  void init_interpretation();
 
};
 

class HECDetectorElement : public CaloDetDescrElement
{
 public:
  HECDetectorElement(const IdentifierHash subcaloHash
             , const IdentifierHash onl1
             , const IdentifierHash onl2
             , const CaloDetDescriptor* descriptor
             , HECCellConstLink& hecCell
             , const HECDetectorRegion* hecRegion
             , bool isTestBeam
             , const GeoAlignmentStore* geoAlignStore = nullptr
             , const CaloElementPositionShift* posShift = nullptr);

  inline HECCellConstLink getHECCell() const {return m_cell;}

  virtual int getLayer() const override;
 
 private:
  HECCellConstLink m_cell;
  const HECDetectorRegion* m_region;
  
  HECDetectorElement() = delete;
  HECDetectorElement & operator=(const HECDetectorElement &right) = delete;

  void init_description(bool isTestBeam
            , const GeoAlignmentStore* geoAlignStore
            , const CaloElementPositionShift* posShift);

  void init_interpretation();
 
};

class FCALDetectorElement : public CaloDetDescrElement
{
 public:
  FCALDetectorElement(const IdentifierHash subcaloHash
              , const IdentifierHash onl1
              , const IdentifierHash onl2
              , const CaloDetDescriptor* descriptor
              , const FCALTile* fcalTile
              , const FCALModule* fcalModule
              , bool isTestBeam
              , const GeoAlignmentStore* geoAlignStore = nullptr
              , const CaloElementPositionShift* posShift = nullptr);

  inline const FCALTile* getFCALTile() const {return m_tile;}

  virtual int getLayer() const override;
 
 private:
  const FCALTile*   m_tile;
  const FCALModule* m_module;
  
  FCALDetectorElement() = delete;
  FCALDetectorElement & operator=(const FCALDetectorElement &right) = delete;

  void init_description(bool isTestBeam
            , const GeoAlignmentStore* geoAlignStore
            , const CaloElementPositionShift* posShift);
 
};

class TileDetectorElement : public CaloDetDescrElement
{
 public:
  TileDetectorElement(const IdentifierHash subcaloHash
              , const IdentifierHash onl1
              , const IdentifierHash onl2
              , const CaloDetDescriptor* descriptor);

  void set_cylindric_raw(double eta_raw
             , double phi_raw
             , double r_raw);

  void set_cylindric(double eta
                     , double phi
                     , double r);

  inline void set_deta(double deta) {m_deta = static_cast<float> (deta);}
  inline void set_dphi(double dphi) {m_dphi = static_cast<float> (dphi);}
  inline void set_dz(double dz) {m_dz = static_cast<float> (dz);}
  inline void set_dr(double dr) {m_dr = static_cast<float> (dr);}
  inline void set_z(double z) {m_z = static_cast<float> (z);}
  inline void set_r(double r) {
    m_r = static_cast<float> (r);
    m_x = static_cast<float> (m_cosPhi * r);
    m_y = static_cast<float> (m_sinPhi * r);
  }
};

class MbtsDetectorElement : public CaloDetDescrElement
{
 public:
  MbtsDetectorElement();
 
  // Initialize data members
  inline void set_id(const Identifier& id) {m_id = id;}
  inline void set_z(double z) {m_z = static_cast<float> (z);}
  inline void set_dz(double dz) {m_dz = static_cast<float> (dz);}
  inline void set_r(double r) {m_r = static_cast<float> (r);}
  inline void set_dr(double dr) {m_dr = static_cast<float> (dr);}
  inline void set_eta(double eta) {m_eta = static_cast<float> (eta);}
  inline void set_deta(double deta) {m_deta = static_cast<float> (deta);}
  inline void set_phi(double phi) {m_phi = static_cast<float> (phi);}
  inline void set_dphi(double dphi) {m_dphi = static_cast<float> (dphi);}
 
  // Compute derived quantities
  void compute_derived();
 
 private:
  Identifier m_id;
 
  virtual Identifier customID() const;
};
 
 

class CaloSuperCellDetectorElement
  : public CaloDetDescrElement
{
public:
  CaloSuperCellDetectorElement(const IdentifierHash subcaloHash,
                               const CaloDetDescriptor* descriptor)
    : CaloDetDescrElement (subcaloHash, 0, 0, descriptor) {}
 

  StatusCode update (const std::vector<const CaloDetDescrElement*>& fromelts);
 
 
private:
  StatusCode updateNull();
 

  StatusCode updateSingle (const CaloDetDescrElement* fromelt);
 

  StatusCode
  updateBE (const std::vector<const CaloDetDescrElement*>& fromelts);
 

  StatusCode
  updateFCAL (const std::vector<const CaloDetDescrElement*>& fromelts);
};
 

class DummyDetDescrElement : public CaloDetDescrElement
{
 public:
  DummyDetDescrElement(const IdentifierHash subcaloHash,
                       const IdentifierHash onl1,
                       const IdentifierHash onl2,
                       const CaloDetDescriptor* descriptor);

  void set_cylindric(double eta, 
             double phi, 
             double r);

  void set_cylindric_raw(double eta_raw, 
             double phi_raw, 
             double r_raw);

  void set_cylindric_size(double deta, 
              double dphi,  
              double dr);
};
 
#endif