PadDesign.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#include <MuonReadoutGeometryR4/PadDesign.h>
#include <GaudiKernel/SystemOfUnits.h>
namespace {
    constexpr double tolerance = 0.001 * Gaudi::Units::mm;
}
#define ORDER_PROP(PROP)                                            \
      {                                                             \
        if (std::abs(PROP - other.PROP) > tolerance) {              \
            return PROP < other.PROP;                               \
        }                                                           \
      }
 
namespace MuonGMR4{
    void PadDesign::print(std::ostream& ostr) const {
        ostr<<"Dimension  -- width x height [mm]: "<<halfWidth() * Gaudi::Units::mm<<" x ";
        ostr<<shortHalfHeight()<<"/"<<longHalfHeight()<<" [mm], ";
        if (hasStereoAngle()) ostr<<"stereo angle: "<<stereoAngle() / Gaudi::Units::deg<<", ";
        ostr<<"position first pad "<<Amg::toString(stripPosition(firstStripNumber()),1);
        ostr<<" *** Trapezoid edges "<<Amg::toString(cornerBotLeft(),1)<<" - "<<Amg::toString(cornerBotRight(), 1)<<" --- ";
        ostr<<Amg::toString(cornerTopLeft(), 1)<<" - "<<Amg::toString(cornerTopRight(), 1);
        ostr<<" -- numPadEta: "<<numPadEta()<<", numPadPhi: "<<numPadPhi()<<", pad height: "<<padHeight();
        ostr<<" -- firstPadPhiDiv: "<<firstPadPhiDiv()<<", anglePadPhi: "<<anglePadPhi()<<", padPhiShift: "<<padPhiShift();
        ostr<<" -- firstPadHeight: "<<firstPadHeight()<<", beamlineRadius: "<<beamlineRadius();
    }
 
    bool PadDesign::operator<(const PadDesign& other) const {
        ORDER_PROP(numPadPhi());
        ORDER_PROP(numPadEta());
        ORDER_PROP(firstPadPhiDiv());
        ORDER_PROP(padPhiShift());
        return static_cast<const StripDesign&>(*this) < other;
    }
 
    void PadDesign::definePadRow(const double firstPadPhiDiv,
                                const int numPadPhi,
                                const double anglePadPhi,
                                const int padPhiShift) {
 
        m_firstPadPhiDiv = firstPadPhiDiv;
        m_numPadPhi = numPadPhi;
        m_anglePadPhi = anglePadPhi;
        m_padPhiShift = static_cast<double>(padPhiShift);
    };
 
    void PadDesign::definePadColumn(const double firstPadHeight,
                                const int numPadEta,
                                const double padHeight,
                                const int maxPadEta) {
 
        m_firstPadHeight = firstPadHeight;
        m_numPadEta = numPadEta;
        m_padHeight = padHeight;
        m_maxPadEta = maxPadEta;
    };
 
    void PadDesign::defineBeamlineRadius(const double radius) {
        m_radius = radius;
    };
 
    localCornerArray PadDesign::padCorners(const std::pair<int, int>& padEtaPhi) const {
        int padEta = padEtaPhi.first;
        int padPhi = padEtaPhi.second;
        localCornerArray padCorners = {make_array<Amg::Vector2D, 4>(Amg::Vector2D::Zero())};
        double botEdge{0.};
        double topEdge{0.};
        double botLeftPoint{0.};
        double botRightPoint{0.};
        double topLeftPoint{0.};
        double topRightPoint{0.};
        double maxBottom = -halfWidth();
        double maxTop = halfWidth();
        if(padEta == 1) {
            botEdge = maxBottom;
            topEdge = botEdge + firstPadHeight();
        }
        else if(padEta > 1 && padEta <= numPadEta()) {
            botEdge = maxBottom + firstPadHeight() + (padEta - 2) * padHeight();
            topEdge = botEdge + padHeight();
            if(padEta == numPadEta()) topEdge = maxTop;
        }
        else {
            ATH_MSG_WARNING(__FILE__<<":"<<__LINE__<<" The given pad Eta " << padEta << " is out of range. Maximum range is " << numPadEta());
            return padCorners;
        }

        double phiRight = firstPadPhiDiv() + (padPhi - 2) * anglePadPhi();
        double phiLeft =  firstPadPhiDiv() + (padPhi - 1) * anglePadPhi();
        double tanRight = std::tan(phiRight * Gaudi::Units::deg);
        double tanLeft = std::tan(phiLeft * Gaudi::Units::deg);
        double botBase = botEdge + beamlineRadius();
        double topBase = topEdge + beamlineRadius();
        botLeftPoint = -botBase * tanLeft;
        botRightPoint = -botBase * tanRight;
        topLeftPoint = -topBase * tanLeft;
        topRightPoint = -topBase * tanRight;
        double cosRight = botBase / std::hypot(botRightPoint, botBase);
        double cosLeft = botBase / std::hypot(botLeftPoint, botBase);
        botLeftPoint += padPhiShift() * cosLeft;
        botRightPoint += padPhiShift() * cosRight;
        topLeftPoint += padPhiShift() * cosLeft;
        topRightPoint += padPhiShift() * cosRight;
        double adjHeight = yCutout()? 2 * halfWidth() - yCutout() : 2 * halfWidth();
        if(padPhi == 1) {
            botRightPoint = shortHalfHeight() + ((longHalfHeight() - shortHalfHeight()) * (botEdge - maxBottom) / adjHeight);
            topRightPoint = shortHalfHeight() + ((longHalfHeight() - shortHalfHeight()) * (topEdge - maxBottom) / adjHeight);
        }
        if(padPhi == numPadPhi()) {
            botLeftPoint = -shortHalfHeight() - ((longHalfHeight() - shortHalfHeight()) * (botEdge - maxBottom) / adjHeight);
            topLeftPoint = -shortHalfHeight() - ((longHalfHeight() - shortHalfHeight()) * (topEdge - maxBottom) / adjHeight);
        }
 
        if (yCutout() && topEdge > (halfWidth() - yCutout())) {
            if (padPhi == 1) {
                topRightPoint = longHalfHeight();
                if (botEdge > (halfWidth() - yCutout())) botRightPoint = longHalfHeight();
            }
            if (padPhi == numPadPhi()) {
                topLeftPoint = -longHalfHeight();
                if (botEdge > (halfWidth() - yCutout())) botLeftPoint = -longHalfHeight();
            }
        }

        if (botEdge > topEdge) {
        ATH_MSG_VERBOSE("Swap top and bottom side "<<padEtaPhi.first<<"/"<<padEtaPhi.second);
        std::swap(botEdge, topEdge);
        }
        if (botLeftPoint > botRightPoint) {
            ATH_MSG_VERBOSE("Swap bottom left and right points "<<padEtaPhi.first<<"/"<<padEtaPhi.second);
            std::swap(botLeftPoint, botRightPoint);
        }
        if (topLeftPoint > topRightPoint) {
            ATH_MSG_VERBOSE("Swap top left and right points "<<padEtaPhi.first<<"/"<<padEtaPhi.second);
            std::swap(topLeftPoint, topRightPoint);
        }
        padCorners[botLeft] = Amg::Vector2D(botLeftPoint, botEdge);
        padCorners[botRight] = Amg::Vector2D(botRightPoint, botEdge);
        padCorners[topLeft] = Amg::Vector2D(topLeftPoint, topEdge);
        padCorners[topRight] = Amg::Vector2D(topRightPoint, topEdge);
        return padCorners; 
    };
}