MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/sTgcReadoutElement.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONREADOUTGEOMETRY_STGCREADOUTELEMENT_H
#define MUONREADOUTGEOMETRY_STGCREADOUTELEMENT_H
 
#include <string>
#include <utility>
 
#include "MuonIdHelpers/sTgcIdHelper.h"
#include "MuonReadoutGeometry/MuonChannelDesign.h"
#include "MuonReadoutGeometry/MuonClusterReadoutElement.h"
#include "MuonReadoutGeometry/MuonDetectorManager.h"
#include "MuonReadoutGeometry/MuonPadDesign.h"
#include "GeoModelInterfaces/IGeoDbTagSvc.h"
#include "GeoPrimitives/GeoPrimitivesToStringConverter.h"
 
 
namespace MuonGMR4 {
    class ReadoutGeomCnvAlg;
}
 
namespace MuonGM {
    class sTgcReadoutElement final : public MuonClusterReadoutElement {
    public:
        friend class MuonGMR4::ReadoutGeomCnvAlg;
   
        sTgcReadoutElement(GeoVFullPhysVol* pv, 
                           const std::string& stName, 
                           int zi, int fi, int mL, MuonDetectorManager* mgr);
 
        ~sTgcReadoutElement();
 
        virtual bool containsId(const Identifier& id) const override final;
 
        virtual double distanceToReadout(const Amg::Vector2D& pos, const Identifier& id) const override final;
 
        virtual int stripNumber(const Amg::Vector2D& pos, const Identifier& id) const override final;
 
        double channelPitch(const Identifier& id) const;
 
        virtual bool stripPosition(const Identifier& id, Amg::Vector2D& pos) const override final;
 
        bool stripGlobalPosition(const Identifier& id, Amg::Vector3D& gpos) const;
 
        int padNumber(const Amg::Vector2D& pos, const Identifier& id) const;
 
        int wireNumber(const Amg::Vector2D& pos, const Identifier& id) const;
 
        double wirePitch(int gas_gap = 1) const;
 
        double positionFirstWire(const Identifier& id) const;
 
        int numberOfWires(const Identifier& id) const;
 
        bool padPosition(const Identifier& id, Amg::Vector2D& pos) const;
 
        bool padGlobalPosition(const Identifier& id, Amg::Vector3D& gpos) const;
 
        bool padCorners(const Identifier& id, std::array<Amg::Vector2D,4>& corners) const;
 
        bool padGlobalCorners(const Identifier& id, std::array<Amg::Vector3D, 4>& gcorners) const;
 
        bool isEtaZero(const Identifier& id, const Amg::Vector2D& localPosition) const;
 
        virtual int numberOfLayers(bool) const override final;
 
        virtual int numberOfStrips(const Identifier& layerId) const override final;
        virtual int numberOfStrips(int, bool measuresPhi) const override final;
 
        int numberOfPads(const Identifier& layerId) const;
 
        int maxPadNumber(const Identifier& layerId) const;
 
        virtual bool spacePointPosition(const Identifier& phiId, const Identifier& etaId, Amg::Vector2D& pos) const override final;
 
        virtual bool spacePointPosition(const Identifier& phiId, const Identifier& etaId, Amg::Vector3D& pos) const override final;
 
        void spacePointPosition(const Amg::Vector2D& phiPos, const Amg::Vector2D& etaPos, Amg::Vector2D& pos) const;
 
        void spacePointPosition(const Identifier& layerId, double locXpos, double locYpos, Amg::Vector3D& pos) const;
 
        Amg::Vector3D localToGlobalCoords(const Amg::Vector3D& locPos, Identifier id) const;
 
        virtual void fillCache() override final;
 
        virtual int surfaceHash(const Identifier& id) const override final;
 
        int surfaceHash(int gasGap, int channelType) const;
 
        virtual int layerHash(const Identifier& id) const override;
 
        // int layerHash( int gasGap) const;        // does not fit in the scheme ( layer hash needs to follow surface hash )
 
        virtual int boundaryHash(const Identifier& id) const override final;
 
        virtual bool measuresPhi(const Identifier& id) const override final;
 
        void initDesign(double thickness);
 
        const MuonChannelDesign* getDesign(const Identifier& id) const;
 
        const MuonChannelDesign* getDesign(int gasGap, int channelType) const;
 
        const MuonPadDesign* getPadDesign(const Identifier& id) const;
        MuonPadDesign* getPadDesign(const Identifier& id);
 
        const MuonPadDesign* getPadDesign(int gasGap) const;
 
        void setChamberLayer(int ml) { m_ml = ml; }
 
        // double getSectorOpeningAngle(bool isLargeSector);
 
        const Amg::Transform3D& getDelta() const { return m_delta; }
        void setDelta(const ALinePar& aline);
       
        void setBLinePar(const BLinePar& bLine);
    
        void posOnDefChamber(Amg::Vector3D& locPosML) const;
 
        bool  has_ALines() const { return (m_ALinePar != nullptr); }
        bool  has_BLines() const { return (m_BLinePar != nullptr); }
        const ALinePar* getALinePar() const { return m_ALinePar; }
        const BLinePar* getBLinePar() const { return m_BLinePar; }
        void  clearALinePar();
        void  clearBLinePar() { m_BLinePar = nullptr; }
 
        // Amdb local (szt) to global coord
        virtual Amg::Vector3D AmdbLRSToGlobalCoords(const Amg::Vector3D& x) const override final { return AmdbLRSToGlobalTransform()*x; }
        virtual Amg::Transform3D AmdbLRSToGlobalTransform() const override final { return absTransform()*Amg::getTranslateZ3D(m_offset)*getDelta(); }
        // Global to Amdb local (szt) coord
        virtual Amg::Vector3D GlobalToAmdbLRSCoords(const Amg::Vector3D& x) const override final { return GlobalToAmdbLRSTransform()*x; }
        virtual Amg::Transform3D GlobalToAmdbLRSTransform() const override final { return AmdbLRSToGlobalTransform().inverse(); }
 
        static double triggerBandIdToRadius(bool isLarge, int triggerBand); 
        
    private:
        
        void initDesignFromSQLite(double thickness);
        
        void initDesignFromAGDD(double thickness);
 
        const sTgcIdHelper& m_idHelper{idHelperSvc()->stgcIdHelper()};
        std::array<MuonChannelDesign,4> m_phiDesign{};
        std::array<MuonChannelDesign,4> m_etaDesign{};
        std::array<MuonPadDesign,4> m_padDesign{};
 
        static constexpr int m_nlayers{4};
        int    m_ml{0};
        double m_offset{0.};
        
        bool   m_diamondShape{false};
        const ALinePar*  m_ALinePar{nullptr};
        const BLinePar*  m_BLinePar{nullptr};
        Amg::Transform3D m_delta{Amg::Transform3D::Identity()};
 
        // transforms (RE->layer)
        std::array<Amg::Transform3D, 4> m_Xlg{make_array<Amg::Transform3D, 4>(Amg::Transform3D::Identity())};
 
 
        // The radial positions of the trigger bands cannot be derived from the readout geometry, therefore hard-coding them for now
        //Position of Band IDs for Large sectors
        static constexpr std::array<double, 94>  LBANDIDSP = {
          10, 10, 10, 10, 10, 10, 1156.72, 1187.98, 1231.87, 1271.34, 1312.87, 1354.56, 1396.38, 1438.07,
          1479.73, 1521.44, 1563.11, 1604.8, 1646.48, 1688.02, 1729.84, 1771.51, 1813.2, 1854.89, 1896.57, 1938.26, 1979.93, 2021.61,
          2063.14, 2104.98, 2146.55, 2181.64, 2209.01, 2251.65, 2282.54, 2313.27, 2356.24, 2396.73, 2438.29, 2480.09, 2521.75, 2563.46,
          2605.11, 2646.85, 2688.48, 2730.19, 2771.86, 2813.41, 2855.21, 2896.93, 2938.61, 2980.26, 3021.95, 3063.63, 3105.31, 3146.98,
          3188.85, 3230.37, 3272.05, 3313.77, 3353.77, 3376.19, 3426.09, 3464.49, 3506.78, 3563.91, 3589.03, 3626.17, 3667.84, 3709.56,
          3751.33, 3792.92, 3834.58, 3876.27, 3917.9, 3959.62, 4001.29, 4043.03, 4084.66, 4126.39, 4168.05, 4209.74, 4251.38, 4293.16,
          4334.72, 4376.47, 4417.54, 4459.75, 4496.31, 4543.27, 4584.77, 4626.47, 4668.25, 4701.14
        };
 
        //Position of Band IDs for Small Sector
        static constexpr std::array<double, 92> SBANDIDSP  = {
          10.0, 10.0, 958.077, 998.248, 1037.405, 1076.535, 1115.69, 1154.82, 1193.97, 1233.135, 1272.265, 1311.395,
          1350.59, 1389.705, 1428.865, 1468.01, 1507.175, 1546.305, 1585.435, 1624.58, 1663.71, 1702.895, 1742.055,
          1781.165, 1820.315, 1859.44, 1898.575, 1937.75, 1976.885, 2016.04, 2055.15, 2094.345, 2136.125, 2172.61,
          2217.68, 2255.125, 2316.115, 2348.91, 2388.06, 2427.245, 2466.385, 2505.515, 2544.69, 2583.8, 2622.99,
          2662.115, 2701.31, 2740.395, 2779.55, 2818.715, 2857.905, 2897.0, 2936.185, 2975.315, 3014.47, 3053.615,
          3092.775, 3131.895, 3171.075, 3210.225, 3249.375, 3288.485, 3317.74, 3347.075, 3396.65, 3440.175, 3475.575,
          3540.81, 3581.97, 3621.13, 3660.285, 3699.41, 3738.535, 3777.73, 3816.89, 3856.055, 3895.105, 3934.3, 3974.34,
          4012.565, 4051.71, 4090.865, 4130.04, 4169.145, 4208.285, 4247.55, 4286.65, 4320.075, 4364.84, 4404.12, 4443.14, 4482.29
        };
 
    };
 
    inline int sTgcReadoutElement::surfaceHash(const Identifier& id) const {
        return surfaceHash(m_idHelper.gasGap(id), m_idHelper.channelType(id));
    }
 
    inline int sTgcReadoutElement::surfaceHash(int gasGap, int channelType) const {
        return (gasGap - 1) * 3 + (2 - channelType);  // assumes channelType=2 (wires), 1(strips), 0(pads)
    }
 
    inline int sTgcReadoutElement::layerHash(const Identifier& id) const {
        return surfaceHash(id);  // don't have a choice here : rewrite MuonClusterReadoutElement first
    }
 
    inline int sTgcReadoutElement::boundaryHash(const Identifier& id) const {
        int iphi = m_idHelper.channelType(id) != sTgcIdHelper::sTgcChannelTypes::Strip;  // wires and pads have locX oriented along phi
        if (std::abs(getStationEta()) < 3) iphi += 2 * (m_idHelper.gasGap(id) - 1);
        return iphi;
    }
 
    inline bool sTgcReadoutElement::measuresPhi(const Identifier& id) const { 
        return (m_idHelper.channelType(id) != sTgcIdHelper::sTgcChannelTypes::Strip); 
    }
 
    inline const MuonChannelDesign* sTgcReadoutElement::getDesign(const Identifier& id) const {
        if (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Strip) return &(m_etaDesign[m_idHelper.gasGap(id) - 1]);
        if (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Wire) return &(m_phiDesign[m_idHelper.gasGap(id) - 1]);
        return nullptr;
    }
 
    inline const MuonPadDesign* sTgcReadoutElement::getPadDesign(const Identifier& id) const {
        if (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Pad) return &(m_padDesign[m_idHelper.gasGap(id) - 1]);
        return nullptr;
    }
 
    inline MuonPadDesign* sTgcReadoutElement::getPadDesign(const Identifier& id) {
        if (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Pad) return &(m_padDesign[m_idHelper.gasGap(id) - 1]);
        return nullptr;
    }
 
    inline const MuonChannelDesign* sTgcReadoutElement::getDesign(int gasGap, int channelType) const {
        if (channelType == 1) return &(m_etaDesign[gasGap - 1]);
        if (channelType == 2) return &(m_phiDesign[gasGap - 1]);
        return 0;
    }
 
    inline const MuonPadDesign* sTgcReadoutElement::getPadDesign(int gasGap) const { return &(m_padDesign[gasGap - 1]); }
 
    inline double sTgcReadoutElement::distanceToReadout(const Amg::Vector2D& pos, const Identifier& id) const {
        const MuonChannelDesign* design = getDesign(id);
        if (!design) return -10000.;
        return design->distanceToReadout(pos);
    }
 
    inline int sTgcReadoutElement::stripNumber(const Amg::Vector2D& pos, const Identifier& id) const {
        if (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Pad) return padNumber(pos, id);
 
        const MuonChannelDesign* design = getDesign(id);
        if (!design) {
            ATH_MSG_WARNING("Cannot associate the strip number for "<<Amg::toString(pos)<<" in layer "
                            <<idHelperSvc()->toStringGasGap(id));
            return -1;
        }
        return design->channelNumber(pos);
    }
 
    inline bool sTgcReadoutElement::stripPosition(const Identifier& id, Amg::Vector2D& pos) const {
        if (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Pad) return padPosition(id, pos);
 
        const MuonChannelDesign* design = getDesign(id);
        if (!design) {
            ATH_MSG_WARNING("Cannot determine the strip postion for "<<idHelperSvc()->toString(id));
            return false;
        }
        return design->center(m_idHelper.channel(id), pos);
    }
 
    inline bool sTgcReadoutElement::stripGlobalPosition(const Identifier& id, Amg::Vector3D& gpos) const {
        Amg::Vector2D lpos{Amg::Vector2D::Zero()};
        if (!stripPosition(id, lpos)) return false;
        surface(id).localToGlobal(lpos, Amg::Vector3D::Zero(), gpos);
        return true;
    }
 
    inline bool sTgcReadoutElement::padPosition(const Identifier& id, Amg::Vector2D& pos) const {
        const MuonPadDesign* design = getPadDesign(id);
        if (!design){
            ATH_MSG_WARNING("Cannot determine the pad position for "<<idHelperSvc()->toString(id));
            return false;
        }
        int padEta = m_idHelper.padEta(id);
        int padPhi = m_idHelper.padPhi(id);
 
        return design->channelPosition(std::make_pair(padEta, padPhi), pos);
    }
 
    inline bool sTgcReadoutElement::padGlobalPosition(const Identifier& id, Amg::Vector3D& gpos) const {
        Amg::Vector2D lpos{Amg::Vector2D::Zero()};
        if (!padPosition(id, lpos)) return false;
        surface(id).localToGlobal(lpos, Amg::Vector3D::Zero(), gpos);
        return true;
    }
 
    inline bool sTgcReadoutElement::padCorners(const Identifier& id, std::array<Amg::Vector2D, 4>& corners) const {
        const MuonPadDesign* design = getPadDesign(id);
        if (!design) {
            ATH_MSG_WARNING("Cannot find the pad corners for "<<idHelperSvc()->toString(id));
            return false;
        }
        int padEta = m_idHelper.padEta(id);
        int padPhi = m_idHelper.padPhi(id);
 
        return design->channelCorners(std::make_pair(padEta, padPhi), corners);
    }
 
    inline bool sTgcReadoutElement::padGlobalCorners(const Identifier& id, std::array<Amg::Vector3D, 4>& gcorners) const {
        std::array<Amg::Vector2D,4> lcorners{make_array<Amg::Vector2D, 4>(Amg::Vector2D::Zero())};
        if (!padCorners(id, lcorners)) {
            return false;
        }
        for (size_t c = 0; c < lcorners.size() ; ++c) {
            surface(id).localToGlobal(lcorners[c], Amg::Vector3D::Zero(), gcorners[c]);
        }
        return true;
    }
 
    // This function returns true if we are in the eta 0 region of QL1/QS1
    inline bool sTgcReadoutElement::isEtaZero(const Identifier& id, const Amg::Vector2D& localPosition) const {
        // False if not a QL1 or QS1 quadruplet
        if (std::abs(m_idHelper.stationEta(id)) != 1) return false;
        const MuonChannelDesign* wireDesign = (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Wire) ?
                                                getDesign(id) :
                                                getDesign(m_idHelper.channelID(id,
                                                                            m_idHelper.multilayer(id),
                                                                            m_idHelper.gasGap(id),
                                                                            sTgcIdHelper::sTgcChannelTypes::Wire,
                                                                            1));
        if (!wireDesign) {
            ATH_MSG_WARNING("Cannot determine whether the pos "<<Amg::toString(localPosition)<<" is etaZero() for "
                          <<idHelperSvc()->toString(id));
            return false;
        } 
 
        // Require the x coordinate for strips, and the y coordinate for wires and pads
        double lpos = (m_idHelper.channelType(id) == sTgcIdHelper::sTgcChannelTypes::Strip) ?
                      localPosition.x() : localPosition.y();
        if (lpos < 0.5 * wireDesign->xSize() - wireDesign->wireCutout) return true;
 
        return false;
    }
 
    inline int sTgcReadoutElement::numberOfLayers(bool) const { return m_nlayers; }
 
    inline int sTgcReadoutElement::numberOfStrips(const Identifier& layerId) const {
        return numberOfStrips(m_idHelper.gasGap(layerId) - 1, 
                              m_idHelper.channelType(layerId) == sTgcIdHelper::sTgcChannelTypes::Wire);
    }
 
    inline int sTgcReadoutElement::numberOfStrips(int lay, bool measPhi) const {
        if (lay > -1 && lay < m_nlayers) { return !measPhi ? m_etaDesign[lay].nch : m_phiDesign[lay].nGroups; }
        return -1;
    }
 
    inline int sTgcReadoutElement::numberOfPads(const Identifier& layerId) const {
        const MuonPadDesign* design = getPadDesign(layerId);
        if (!design) {
            ATH_MSG_WARNING("no pad design found when trying to get the number of pads "<<idHelperSvc()->toString(layerId));
            return 0;
        }
        return design->nPadColumns * design->nPadH;
    }
 
    inline int sTgcReadoutElement::maxPadNumber(const Identifier& layerId) const {
        const MuonPadDesign* design = getPadDesign(layerId);
        if (!design) {
            ATH_MSG_WARNING("no pad design found when trying to get the largest pad number "<<idHelperSvc()->toString(layerId));
            return 0;
        }
        return (design->nPadColumns - 1) * m_idHelper.padEtaMax() + design->nPadH;
    }
 
    inline bool sTgcReadoutElement::spacePointPosition(const Identifier& phiId, const Identifier& etaId, Amg::Vector2D& pos) const {
        Amg::Vector2D phiPos{Amg::Vector2D::Zero()}, etaPos{Amg::Vector2D::Zero()};
        if (!stripPosition(phiId, phiPos) || !stripPosition(etaId, etaPos)) return false;
        spacePointPosition(phiPos, etaPos, pos);
        return true;
    }
 
    inline bool sTgcReadoutElement::spacePointPosition(const Identifier& phiId, const Identifier& etaId, Amg::Vector3D& pos) const {
        Amg::Vector2D lpos{Amg::Vector2D::Zero()};
        spacePointPosition(phiId, etaId, lpos);
        surface(phiId).localToGlobal(lpos, pos, pos);
        return true;
    }
 
    inline void sTgcReadoutElement::spacePointPosition(const Amg::Vector2D& phiPos, const Amg::Vector2D& etaPos, Amg::Vector2D& pos) const {
        pos[0] = phiPos.x();
        pos[1] = etaPos.x();
    }
 
 
    inline double sTgcReadoutElement::triggerBandIdToRadius(bool isLarge, int triggerBand) {
        if(isLarge){
            return LBANDIDSP[triggerBand];
        } else {
            return SBANDIDSP[triggerBand]; 
        }
    }; 
 
}  // namespace MuonGM
 
#endif  // MUONREADOUTGEOMETRY_STGCREADOUTELEMENT_H