ARICHGeoHAPD.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Luka Santelj                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <arich/dbobjects/ARICHGeoBase.h>
#include <geometry/dbobjects/GeoOpticalSurface.h>
#include <string>
 
namespace Belle2 {
  class ARICHGeoHAPD: public ARICHGeoBase {
 
  public:
 
    ARICHGeoHAPD()
    {}
 
    void setWallThickness(double thickness) {m_wallThickness = thickness;}
 
    void setWinThickness(double thickness) {m_winThickness = thickness;}
 
    void setPadSize(double padSize) {m_padSize = padSize;}
 
    void setChipGap(double chipGap) {m_chipGap = chipGap;}
 
    void setNPads(int nx, int ny) {m_nPadX = nx; m_nPadY = ny; }
 
    void setWinRefIndex(double refInd) {m_winRefIndex = refInd;}
 
    void setModuleSizeZ(double modZ) {m_moduleSizeZ = modZ;}
 
    void setWallMaterial(const std::string& material) {m_wallMaterial = material;}
 
    void setWinMaterial(const std::string& material) {m_winMaterial = material;}
 
    void setAPDMaterial(const std::string& material) {m_apdMaterial = material;}
 
    void setFEBMaterial(const std::string& material) {m_febMaterial = material;}
 
    void setFillMaterial(const std::string& material) {m_fillMaterial = material;}
 
    void setFEBVolume(double sizeX, double sizeY, double sizeZ,
                      const std::string& material)
    {
      m_FEBSizeX = sizeX;
      m_FEBSizeY = sizeY;
      m_FEBSizeZ = sizeZ;
      m_febMaterial = material;
    }
 
    void setAPDVolume(double sizeX, double sizeY, double sizeZ,
                      const std::string& material, const GeoOpticalSurface& surf)
    {
      m_APDSizeX = sizeX;
      m_APDSizeY = sizeY;
      m_APDSizeZ = sizeZ;
      m_apdMaterial = material;
      m_apdSurface = surf;
    }
 
    void setHAPDVolume(double sizeX, double sizeY, double sizeZ,
                       const std::string& material)
    {
      m_HAPDSizeX = sizeX;
      m_HAPDSizeY = sizeY;
      m_HAPDSizeZ = sizeZ;
      m_fillMaterial = material;
    }
 
 
    double getSizeX() const {return m_HAPDSizeX / s_unit;}
 
    double getSizeY() const {return m_HAPDSizeY / s_unit;}
 
    double getSizeZ() const {return m_HAPDSizeZ / s_unit;}
 
    double getWallThickness() const {return m_wallThickness / s_unit;}
 
    double getWinThickness() const {return m_winThickness / s_unit;}
 
    double getAPDSizeX() const {return m_APDSizeX / s_unit;}
 
    double getAPDSizeY() const {return m_APDSizeY / s_unit;}
 
    double getAPDSizeZ() const {return m_APDSizeZ / s_unit;}
 
    double getFEBSizeX() const {return m_FEBSizeX / s_unit;}
 
    double getFEBSizeY() const {return m_FEBSizeY / s_unit;}
 
    double getFEBSizeZ() const {return m_FEBSizeZ / s_unit;}
 
    double getModuleSizeZ() const {return m_moduleSizeZ / s_unit;}
 
    double getWinRefIndex() const {return m_winRefIndex;}
 
    double getPadSize() const {return m_padSize / s_unit;}
 
    double getChipGap() const {return m_chipGap / s_unit;}
 
    unsigned getNumPadsX() const {return m_nPadX;}
 
    unsigned getNumPadsY() const {return m_nPadY;}
 
 
    const std::string& getWallMaterial() const {return m_wallMaterial;}
 
    const std::string& getFillMaterial() const {return m_fillMaterial;}
 
    const std::string& getWinMaterial() const {return m_winMaterial;}
 
    const std::string& getAPDMaterial() const {return m_apdMaterial;}
 
    const std::string& getFEBMaterial() const {return m_febMaterial;}
 
    const GeoOpticalSurface& getAPDSurface() const {return m_apdSurface;}
 
    void getXYChannel(double x, double y, int& chX, int& chY) const;
 
    unsigned getChipID(double x, double y) const;
 
    void getXYChannelPos(int chX, int chY, double& x, double& y) const;
 
 
    double getPhotocathodeApdDistance() const
    {
      return (m_HAPDSizeZ - m_winThickness - m_wallThickness - m_APDSizeZ) / s_unit;
    }
 
 
    bool isConsistent() const override;
 
    void print(const std::string& title = "HAPD module geometry parameters") const override;
 
  private:
 
    std::string m_wallMaterial; 
    std::string m_winMaterial;  
    std::string m_apdMaterial;  
    std::string m_febMaterial;  
    std::string m_fillMaterial; 
    double m_HAPDSizeX = 0;     
    double m_HAPDSizeY = 0;     
    double m_HAPDSizeZ = 0;     
    double m_wallThickness = 0; 
    double m_winThickness = 0;  
    double m_padSize = 0;       
    double m_chipGap = 0;       
    unsigned m_nPadX = 0;       
    unsigned m_nPadY = 0;       
    double m_winRefIndex = 0;   
    double m_APDSizeX = 0;      
    double m_APDSizeY = 0;      
    double m_APDSizeZ = 0;      
    double m_FEBSizeX = 0;      
    double m_FEBSizeY = 0;      
    double m_FEBSizeZ = 0;      
    double m_moduleSizeZ = 0;   
    GeoOpticalSurface m_apdSurface; 
    ClassDefOverride(ARICHGeoHAPD, 1); 
  };
 
} // end namespace Belle2