WireGroupDesign.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_WIREGROUPDESIGN_H
#define MUONREADOUTGEOMETRYR4_WIREGROUPDESIGN_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 WireGroupDesign;
using WireDesignPtr = GeoModel::TransientConstSharedPtr<WireGroupDesign>;
    
class WireGroupDesign: public StripDesign {
    public:
        WireGroupDesign() = default;
        bool operator<(const WireGroupDesign& other) const;
        void declareGroup(const unsigned int x);
        unsigned int numWiresInGroup(unsigned int groupNum) const;
        unsigned int numPitchesToGroup(unsigned int groupNum) const;
        unsigned int nAllWires() const;
        CheckVector2D wirePosition(unsigned int groupNum, 
                                   unsigned int wireNum) const;
    
        CheckVector2D leftWireEdge(unsigned int groupNum,
                                   unsigned int wireNum) const;
        
        CheckVector2D rightWireEdge(unsigned int groupNum,
                                    unsigned int wireNum) const;
    
        int stripNumber(const Amg::Vector2D& pos) const override;

        std::pair<int, int> wireNumber(const Amg::Vector2D& extPos) const;

        double wireLength(unsigned int groupNum,
                          unsigned int wireNum) const;

        double wireCutout() const; 
        void defineWireCutout(const double wireCutout);
        
    private:
        void print(std::ostream& ostr) const override final; 
        Amg::Vector2D stripPosition(int stripNum) const override final;
        struct wireGroup{
            wireGroup(unsigned int nWires, unsigned int accWires):
                numWires{nWires}, accumlWires{accWires} {} 

            unsigned int numWires{0};
            unsigned int accumlWires{0};
        };
        using wireGrpVector = std::vector<wireGroup>;
        using wireGrpVectorItr = wireGrpVector::const_iterator;
        wireGrpVector m_groups{};
        double m_wireCutout{0};
 
 
};
 
struct WireDesignSorter{
    bool operator()(const WireDesignPtr&a, const WireDesignPtr& b) const {
            return (*a) < (*b);
    }
    bool operator()(const WireGroupDesign&a ,const WireGroupDesign& b) const {
        return a < b;
    }
};
 
using WireGroupDesignSet = std::set<WireDesignPtr, WireDesignSorter>;
 
}
#include <MuonReadoutGeometryR4/WireGroupDesign.icc>
#endif