MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/TgcReadoutElement.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUONREADOUTGEOMETRY_TGCREADOUTELEMENT_H
#define MUONREADOUTGEOMETRY_TGCREADOUTELEMENT_H
 
 
#include "Identifier/Identifier.h"
#include "MuonReadoutGeometry/MuonClusterReadoutElement.h"
#include "MuonReadoutGeometry/MuonDetectorManager.h"
#include "MuonReadoutGeometry/TgcReadoutParams.h"
#include "CxxUtils/ArrayHelper.h"
 
 
namespace MuonGM {

 
    class TgcReadoutElement final : public MuonClusterReadoutElement {
 
    public:
        TgcReadoutElement(GeoVFullPhysVol* pv, const std::string& stName, MuonDetectorManager* mgr);
 
        virtual ~TgcReadoutElement() = default;

        virtual double distanceToReadout(const Amg::Vector2D& pos, const Identifier& id) const override;

        virtual int stripNumber(const Amg::Vector2D& pos, const Identifier& id) const override final;

        virtual bool stripPosition(const Identifier& id, Amg::Vector2D& pos) const override;

        virtual  int layerHash(const Identifier& id) const override;

        virtual  int surfaceHash(const Identifier& id) const override;

        virtual  int boundaryHash(const Identifier& id) const override;

        virtual  bool measuresPhi(const Identifier& id) const override;

        virtual int numberOfLayers(bool isStrip) const override;

        virtual int numberOfStrips(const Identifier& layerId) const override;
        virtual int numberOfStrips(int layer, bool isStrip) const override;

        virtual  bool spacePointPosition(const Identifier& phiId, const Identifier& etaId, Amg::Vector2D& pos) const override;

        virtual  bool spacePointPosition(const Identifier& phiId, const Identifier& etaId, Amg::Vector3D& pos) const override;
 
        Amg::Vector3D localSpacePoint(const Identifier& stripId, 
                                      const Amg::Vector3D& etaHitPos,
                                      const Amg::Vector3D& phiHitPos) const;
        
        virtual bool containsId(const Identifier& id) const override;
        
        
        
        int chamberType() const;
        int nGasGaps() const;
        bool validGap(int gasGap) const;
 

        int nStrips(int gasGap) const;
        bool validStrip(int gasGap, int strip) const;
        double stripWidth(int gasGap, int strip) const;
        double stripLength() const;
        double stripPitch(int gasGap, int strip) const;
        double stripPitch(int gasGap, int strip, double radialPos) const; 
        double stripLowEdgeLocX(int gasGap, int strip, double radialPos) const; 
        double stripHighEdgeLocX(int gasGap, int strip , double radialPos) const;
        double stripCenterLocX(int gasGap, int strip, double radialPos) const;
 
        double physicalDistanceFromBase() const;
        double stripPosOnLargeBase(int strip, int gasGap) const;
        double stripPosOnShortBase(int strip, int gasGap) const;
 
        
        double stripDeltaPhi() const;
        double stripDeltaPhi(int gasGap, int strip) const;
 
        
        bool validGang(int gasGap, int wireGang) const;
        int nWires(int gasGap) const;
        int nWires(int gasGap, int gang) const;
        int nWireGangs(int gasGap) const;
        double nPitchesToGang(int gasGap, int gang) const;
        double gangRadialLength(int gasGap, int gang) const;
        double gangShortWidth(int gasGap, int gang) const;
        double gangCentralWidth(int gasGap, int gang) const;
        double gangLongWidth(int gasGap, int gang) const;
        double gangThickness() const;
        double wireLength(int wire) const;
        double wirePitch() const;
 
 
 
        double stripShortWidth(int, int) const;
        double stripLongWidth(int, int) const;

        const Amg::Transform3D& localToGlobalTransf(const Identifier& id) const;
        Amg::Transform3D globalToLocalTransf(const Identifier& id) const;

        Amg::Vector3D channelPos(const Identifier& id) const;
        Amg::Vector3D channelPos(int gasGap, bool isStrip, int channel) const;

        Amg::Vector3D wireGangPos(const Identifier& id) const;
        Amg::Vector3D wireGangPos(int gasGap, int gang) const;

        Amg::Vector3D stripPos(int gasGap, int strip) const;
        Amg::Vector3D stripPos(const Identifier& id) const;
        Amg::Vector3D stripDir(int gasGap, int strip) const;
        Amg::Vector3D stripDir(const Identifier& id) const;
 
        
        bool isEndcap() const;
        bool isForward() const;
        bool isDoublet() const;
        bool isTriplet() const;
        double chamberLocPhiMin() const;
        double chamberLocPhiMax() const;
        
 
        double length() const;
 
        double frameZwidth() const;
        double frameXwidth() const;
        double chamberWidth(double z) const;
 
 
        int nPhiChambers() const;
        int nPhiSectors() const;
        
        void setPlaneZ(double value, int gasGap);
        
        int findGang(int gasGap, const Amg::Vector3D& extPos) const;
        int findStrip(int gasGap, const Amg::Vector3D& extPos) const;
        
 
 
        void setReadOutName(const std::string& rName);
        void setReadOutParams(GeoModel::TransientConstSharedPtr<TgcReadoutParams> pars);
        
        void setFrameThickness(const double frameH, const double frameAB);
 
       // Access to readout parameters
        const std::string& readOutName() const;
 
        const TgcReadoutParams* getReadoutParams() const;
 
        virtual void fillCache() override;
 
        
    private:
        Amg::Vector3D localChannelPos(const Identifier& id) const;
        Amg::Vector3D localChannelPos(int gasGap, bool isStrip, int channel) const;
        // Returns the local position of a strip
        Amg::Vector3D localStripPos(int gasGap, int strip) const;
        Amg::Vector3D localStripPos(const Identifier& id) const;
        Amg::Vector3D localStripDir(int gasGap, int strip) const;
        Amg::Vector3D localStripDir(const Identifier& id) const;

        Amg::Vector3D localWireGangPos(int gasGap, int gang) const;
        Amg::Vector3D localWireGangPos(const Identifier& id) const;
        Amg::Vector3D localGasGapPos(int gg) const;
 
        static int surfaceHash(int GasGap, bool isStrip);

        double stripLocalX(const int gasGap,
                           const int stripNum,
                           const double locY,
                           const double refPoint) const;
        
        double wireGangLocalX(const int gasGap,
                              const int gangNum) const;

        double wireGangBottomX(int gasGap, int gangNum) const;
        double wireGangTopX(int gasGap, int gangNum) const;
        
        std::pair<double, int> stripNumberToFetch(int gasGap, int inStrip) const;
 
        const TgcIdHelper& m_idHelper{idHelperSvc()->tgcIdHelper()};
 
        static constexpr int s_maxGasGap{3};
        std::array<double, s_maxGasGap> m_gasPlaneZ{make_array<double, s_maxGasGap>(-9999.)};
 
        std::string m_readout_name{};
        GeoModel::TransientConstSharedPtr <TgcReadoutParams> m_readoutParams{nullptr};
        double m_stripSlope{0.};
        double m_locMinPhi{0.};
        double m_locMaxPhi{0.};
 
        int m_stIdxT4E{m_idHelper.stationNameIndex("T4E")};
 
        double m_frameH{0.};
        double m_frameAB{0.};
 
    };
 
}  // namespace MuonGM
#include <MuonReadoutGeometry/TgcReadoutElement.icc>
#endif  // MUONREADOUTGEOMETRY_TGCREADOUTELEMENT_H