sTgcReadoutElement.icc

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_STGCREADOUTELEMENT_ICC
#define MUONREADOUTGEOMETRYR4_STGCREADOUTELEMENT_ICC
 
 
namespace ActsTrk{
    template <> inline Amg::Transform3D 
        TransformCacheDetEle<MuonGMR4::sTgcReadoutElement>::fetchTransform(const DetectorAlignStore* store) const{
            return  m_parent->toStation(store) * m_parent->fromGapToChamOrigin(hash());
        }
    template <> inline Identifier
        TransformCacheDetEle<MuonGMR4::sTgcReadoutElement>::identify() const {
            using RE = MuonGMR4::sTgcReadoutElement;
            return m_parent->measurementId(RE::createHash(RE::gasGapNumber(hash()) + 1,
                                                          RE::chType(hash()), 1));
    }
}
 
 
namespace MuonGMR4 {
    namespace sTgcIdMeasHashFields {
        constexpr unsigned minusOne = -1;
        /// Hash field layout
        //// (wireInGrp)  | (channel)  [1-512]  | IsWireInGrp (0,1) |  (sTgcChannelType) [0-2] | (gasGap -1) [1-4]
        constexpr unsigned gasGapShift = 0;
        constexpr unsigned chTypeShift = 2;        
        constexpr unsigned wireInGrpBit = chTypeShift+ 2;
        constexpr unsigned chanShift = chTypeShift + 3;        
        constexpr unsigned wireInGrpShift = chanShift + 9;
    }
    
inline double sTgcReadoutElement::chamberHeight() const { return 2.* m_pars.halfChamberHeight; }
inline double sTgcReadoutElement::sChamberLength() const { return 2.* m_pars.sHalfChamberLength; }
inline double sTgcReadoutElement::lChamberLength() const { return 2.* m_pars.lHalfChamberLength; }
inline double sTgcReadoutElement::thickness() const { return 2.* m_pars.halfChamberTck; }
inline double sTgcReadoutElement::sFrameWidth() const { return m_pars.sFrameWidth; }
inline double sTgcReadoutElement::lFrameWidth() const { return m_pars.lFrameWidth; }
inline int sTgcReadoutElement::multilayer() const { return m_multiLayer; }
 
inline unsigned sTgcReadoutElement::numLayers() const { return m_pars.numLayers; }
inline unsigned sTgcReadoutElement::nChTypes() const { return m_pars.nChTypes; }
inline double sTgcReadoutElement::gasGapThickness() const { return m_pars.gasTck; }
inline double sTgcReadoutElement::gasGapPitch() const { return m_gasGapPitch; }
 
inline double sTgcReadoutElement::firstStripPitch(const IdentifierHash& measHash) const { 
    return m_pars.firstStripPitch[gasGapNumber(measHash)]; }
 
inline unsigned sTgcReadoutElement::numChannels(const IdentifierHash& measHash) const {
    switch (chType(measHash)) {
        case ReadoutChannelType::Strip:
            return stripDesign(measHash).numStrips();
        case ReadoutChannelType::Pad:
            return padDesign(measHash).numPads();
        case ReadoutChannelType::Wire:
            return wireDesign(measHash).numStrips();
        case ReadoutChannelType::WireInGrp:
            return wireDesign(measHash).nAllWires(); 
        default:
            ATH_MSG_WARNING("Channel type "<<chType(measHash)<<" is unknown.");
    }
    return 0;
}
 
inline double sTgcReadoutElement::stripLength(const IdentifierHash& measHash) const { 
    return stripDesign(measHash).stripLength(channelNumber(measHash)); 
}
 
inline unsigned sTgcReadoutElement::numPadEta(const IdentifierHash& measHash) const { return padDesign(measHash).numPadEta(); }
inline unsigned sTgcReadoutElement::numPadPhi(const IdentifierHash& measHash) const { return padDesign(measHash).numPadPhi(); }
 
inline double sTgcReadoutElement::padHeight(const IdentifierHash& measHash) const { return (padEta(measHash) == 1) ? 
                                                                                    padDesign(measHash).firstPadHeight() : padDesign(measHash).padHeight(); }
 
inline double sTgcReadoutElement::padPhiShift(const IdentifierHash& measHash) const { return padDesign(measHash).padPhiShift(); }
 
inline double sTgcReadoutElement::anglePadPhi(const IdentifierHash& measHash) const { return padDesign(measHash).anglePadPhi(); }
inline unsigned sTgcReadoutElement::maxPadEta(const IdentifierHash& measHash) const { return padDesign(measHash).anglePadPhi(); }
 
inline unsigned sTgcReadoutElement::padNumberSeq(const IdentifierHash& measHash) const { 
    int padEta = padEtaPhi(measHash).first;
    int padPhi = padEtaPhi(measHash).second;
    int numPadEta = padDesign(measHash).numPadEta();
    int numPadPhi = padDesign(measHash).numPadPhi();
    if((padEta > numPadEta) || (padPhi > numPadPhi)) {
        ATH_MSG_WARNING(__FILE__<<":"<<__LINE__<<" Pad Channel: "<< idHelperSvc()->toString(measurementId(measHash))
                        << " resulted in invalid (eta, phi): ("<<padEta<< ", " <<padPhi <<"), (maxEta, maxPhi): "<<numPadEta<<", "<<numPadPhi
                        <<"). total pads available: "<<padDesign(measHash).numPads());
        return 0;
    }
    return (padPhi - 1) * numPadEta + padEta; 
}
 
inline std::pair<unsigned, unsigned> sTgcReadoutElement::padEtaPhi(const IdentifierHash& measHash) const { return padDesign(measHash).padEtaPhi(channelNumber(measHash)); }
inline unsigned sTgcReadoutElement::padEta(const IdentifierHash& measHash) const { return padDesign(measHash).padEta(channelNumber(measHash)); }
inline unsigned sTgcReadoutElement::padPhi(const IdentifierHash& measHash) const { return padDesign(measHash).padPhi(channelNumber(measHash)); }
inline double sTgcReadoutElement::beamlineRadius(const IdentifierHash& measHash) const { return padDesign(measHash).beamlineRadius(); }
inline sTgcReadoutElement::localCornerArray sTgcReadoutElement::localPadCorners(const IdentifierHash& measHash) const { return padDesign(measHash).padCorners(channelNumber(measHash)); }
inline const StripDesign& sTgcReadoutElement::stripDesign(const IdentifierHash& measHash) const {
    return m_pars.stripLayers[gasGapNumber(measHash)].design();
 }
 
inline const WireGroupDesign& sTgcReadoutElement::wireDesign(const IdentifierHash& measHash) const {
    return static_cast<const WireGroupDesign&>(m_pars.stripLayers[gasGapNumber(measHash)].design(true));
 }
inline const PadDesign& sTgcReadoutElement::padDesign(const IdentifierHash& measHash) const {
    return static_cast<const PadDesign&>(m_pars.padLayers[gasGapNumber(measHash)].design(true));
 }
inline IdentifierHash sTgcReadoutElement::measurementHash(const Identifier& measId) const {
    if (idHelperSvc()->detElId(measId) != identify()) {
        ATH_MSG_WARNING("The measurement " << idHelperSvc()->toString(measId)
                        << " picks the wrong readout element " << idHelperSvc()->toStringDetEl(identify()));
    }
    return createHash(m_idHelper.gasGap(measId),
                      m_idHelper.channelType(measId),
                      m_idHelper.channel(measId));
}
inline IdentifierHash sTgcReadoutElement::createHash(const unsigned gasGap, 
                                                     const unsigned channelType, 
                                                     const unsigned channel,
                                                     const unsigned wireInGrp) {
    using namespace sTgcIdMeasHashFields;
    /// Construct the Hash such that  (channel) | WireInGrpBit (0,1) | (channelType) | (gasGap -1) 
    if (channelType == ReadoutChannelType::WireInGrp) {
        const unsigned readOutHash = static_cast<unsigned>(createHash(gasGap, ReadoutChannelType::Wire, channel));
        return IdentifierHash {wireInGrp << wireInGrpShift |  (1<<wireInGrpBit) | readOutHash};
    }
    return IdentifierHash{ channel << chanShift | channelType << chTypeShift | (gasGap -1) };
}
 
inline IdentifierHash sTgcReadoutElement::layerHash(const IdentifierHash& measHash) {
    return createHash(gasGapNumber(measHash)+1, ReadoutChannelType::Strip,0);
}
 
inline  unsigned sTgcReadoutElement::channelNumber(const IdentifierHash& measHash) {
    using namespace sTgcIdMeasHashFields;
    constexpr unsigned mask = (minusOne << wireInGrpShift);
    const unsigned stripedHash = (~mask) & static_cast<unsigned>(measHash);
    return stripedHash >> chanShift;
}
inline unsigned sTgcReadoutElement::chType(const IdentifierHash& measHash) {
    using namespace sTgcIdMeasHashFields;
    constexpr unsigned mask = (minusOne << chanShift);
    const unsigned stripedHash = (~mask) & static_cast<unsigned>(measHash);
    return stripedHash >> chTypeShift;
    
}
inline unsigned sTgcReadoutElement::gasGapNumber(const IdentifierHash& measHash) {
    using namespace sTgcIdMeasHashFields;
    constexpr unsigned mask = (minusOne << chTypeShift);
    return (static_cast<unsigned>(measHash) &(~mask) );
}
inline Identifier sTgcReadoutElement::measurementId(const IdentifierHash& measHash) const {
    return m_idHelper.channelID(identify(), multilayer(), gasGapNumber(measHash) + 1, chType(measHash), channelNumber(measHash));  
}
inline IdentifierHash sTgcReadoutElement::layerHash(const Identifier& measId) const {
    if (m_idHelper.elementID(measId) != m_idHelper.elementID(identify()) ) {
        ATH_MSG_WARNING("The measurement " << idHelperSvc()->toString(measId)
                        << " picks the wrong readout element " << idHelperSvc()->toStringDetEl(identify()));
    }
    return createHash(m_idHelper.gasGap(measId), ReadoutChannelType::Strip, 0);
}
 
inline const StripLayer& sTgcReadoutElement::stripLayer(const IdentifierHash& measHash) const {
    unsigned gasGap = gasGapNumber(measHash);
    switch (chType(measHash)) {
        case ReadoutChannelType::Strip:
        case ReadoutChannelType::Wire:
        case ReadoutChannelType::WireInGrp:
            return m_pars.stripLayers.at(gasGap);
        case ReadoutChannelType::Pad:
            return m_pars.padLayers.at(gasGap);
    }
    THROW_EXCEPTION("Invalid channel parsed");
    return stripLayer(0);
}
 
}  // namespace MuonGMR4
#endif