RadialStripDesign.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_RADIALSTRIPDESIGN_H
#define MUONREADOUTGEOMETRYR4_RADIALSTRIPDESIGN_H
 
#include <MuonReadoutGeometryR4/StripDesign.h>
 
namespace MuonGMR4{
 
/* The phi readout channels of the Tgc chambers are an essamble of wires that are readout 
 * together. The wires are equi-distant, but not every wire group in a gas gap contains the
 * same number of wires. Hence, the channel width & pitch vary across the board. The wire group 
 * design accounts for this detector design feature. It inherits from the StripDesign class and 
 * overloads the stripPosition feature to place the measurement onto the middle wire of each group.
 * Nevertheless, it's important to know about the position and the lenghts of each wire seperately
 * as these quantities are fed into the digitization. Therefore, extra methods are added to the design
 * to provide this information as well.
*/
class RadialStripDesign;
using RadialStripDesignPtr = GeoModel::TransientConstSharedPtr<RadialStripDesign>;
    
class RadialStripDesign: public StripDesign {
    public:
        RadialStripDesign() = default;
        bool operator<(const RadialStripDesign& other) const;
        
        void addStrip(const double posOnBottom,
                      const double posOnTop);
        
        Amg::Vector2D stripDir(int stripNumber) const;
        Amg::Vector2D stripLeftEdge(int stripNumber) const; 
        Amg::Vector2D stripRightEdge(int stripNumber) const;
        Amg::Vector2D stripNormal(int stripNumber) const;
        Amg::Vector2D stripLeftBottom(int stripNumber) const;
        Amg::Vector2D stripRightBottom(int stripNumber) const;
        Amg::Vector2D stripLeftTop(int stripNumber) const;
        Amg::Vector2D stripRightTop(int stripNumber) const;

        int numStrips() const override;
        int stripNumber(const Amg::Vector2D& extPos) const override final;
        double stripLength(int stripNumb) const override;
        double stripPitch(int stripNumb) const;
    private:
        CheckVector2D leftInterSect(int stripNum, bool uncapped = false) const override final;
        CheckVector2D rightInterSect(int stripNum, bool uncapped = false) const override final;
        struct stripEdges{
            stripEdges(double dBot, double dTop, 
                       const RadialStripDesign& parent_):
                parent{parent_},
                distOnBottom{dBot},
                distOnTop{dTop} {}
 
            /* Returns the bottom mounting point */
            Amg::Vector2D bottomMounting() const;
            /* Returns the top mounting point */
            Amg::Vector2D topMounting() const;
            /* Returns the center point between bottom & top*/
            Amg::Vector2D center() const;
            /* Returns the connecting vector from bottom top */
            Amg::Vector2D fromBottomToTop() const;
 
            const RadialStripDesign& parent;
            double distOnBottom{0.};
            double distOnTop{0.};
        };
        using stripEdgeVec = std::vector<stripEdges>;
        using stripEdgeVecItr = stripEdgeVec::const_iterator;
        stripEdgeVec m_strips{};
        
        bool m_reversedStripOrder{false};
 
};
 
struct RadialDesignSorter{
    bool operator()(const RadialStripDesignPtr& a, const RadialStripDesignPtr& b) const {
            return (*a) < (*b);
    }
    bool operator()(const RadialStripDesign&a ,const RadialStripDesign& b) const {
        return a < b;
    }
};
 
using RadialStripDesignSet = std::set<RadialStripDesignPtr, RadialDesignSorter>;
 
}
#include <MuonReadoutGeometryR4/RadialStripDesign.icc>
#endif