PadDesign.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_PADDESIGN_H
#define MUONREADOUTGEOMETRYR4_PADDESIGN_H
 
#include <MuonReadoutGeometryR4/StripDesign.h>
#include "CxxUtils/ArrayHelper.h"
#include <GaudiKernel/SystemOfUnits.h>
 
namespace MuonGMR4{
 
/* The pad readout channels in sTgc chambers consists of 2D structure. Chambers may contain
 * different number of pads depending on the gasGap. Measurements in the eta direction for pads
 * are in millimeters, whereas for the phi direction, we have angular measurements in degrees.
 * same number of wires. Hence, the channel width & pitch vary across the board. The wire group 
 * PadDesign inherits from the StripDesign class and contains functions that can define the geometry
 * of pad layers. There are several functions to get the corners and center of each individual pads
 * w.r.t. the center of the chamber. Extra methods will be added according to the needs.
*/
class PadDesign;
using PadDesignPtr = GeoModel::TransientConstSharedPtr<PadDesign>;
    
class PadDesign: public StripDesign {
    public:
        PadDesign() = default;
        bool operator<(const PadDesign& other) const;
 
        void definePadRow(const double firstPadPhiDiv,
                                const int numPadPhi,
                                const double anglePadPhi,
                                const int padPhiShift);
        void definePadColumn(const double firstPadHeight,
                                const int numPadEta,
                                const double padHeight,
                                const int maxPadEta = 18);
        int numPads() const;                        
        double firstPadPhiDiv() const;
        int numPadPhi() const;
        double anglePadPhi() const;
        double padPhiShift() const;
        double firstPadHeight() const;
        int numPadEta() const;
        double padHeight() const;
        int maxPadEta() const;
        int padNumber(const int channel) const;
        std::pair<int, int> padEtaPhi(const int channel) const;
        int padEta(const int channel) const;
        int padPhi(const int channel) const;
        void defineBeamlineRadius(const double radius);
        double beamlineRadius() const;
        using localCornerArray = std::array<Amg::Vector2D, 4>;   
        enum  padCorners{ botLeft=0, botRight, topLeft, topRight };
        localCornerArray padCorners(const int channel) const;
        localCornerArray padCorners(const std::pair<int, int>& padEtaPhi) const;
        Amg::Vector2D stripPosition(int stripNum) const override final;
        std::pair<int, int> channelNumber(const Amg::Vector2D& hitPos) const;
        double maxActiveX(const double locY) const;
        double sectorAngle() const;
 
    private:
        void print(std::ostream& ostr) const override final; 
        double m_firstPadPhiDiv{0.};
        int m_numPadPhi{0};
        double m_anglePadPhi{0.};
        double m_padPhiShift{0.};
        double m_firstPadHeight{0.};
        int m_numPadEta{0};
        double m_padHeight{0.};
        int m_maxPadEta{18};
        double m_radius{0.};
        static constexpr double m_largeSectorAngle{28.};
        static constexpr double m_smallSectorAngle{17.};       
};
 
struct PadDesignSorter{
    bool operator()(const PadDesignPtr&a, const PadDesignPtr& b) const {
            return (*a) < (*b);
    }
    bool operator()(const PadDesign&a ,const PadDesign& b) const {
        return a < b;
    }
};
 
}
#include <MuonReadoutGeometryR4/PadDesign.icc>
#endif