MmReadoutElement.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_MMREADOUTELEMENT_H
#define MUONREADOUTGEOMETRYR4_MMREADOUTELEMENT_H
 
#include <MuonReadoutGeometryR4/MuonReadoutElement.h>
#include <MuonReadoutGeometryR4/StripDesign.h>
#include <MuonReadoutGeometryR4/StripLayer.h>
 
#ifndef SIMULATIONBASE
#   include "Acts/Utilities/BoundFactory.hpp"
#endif
 
 
namespace MuonGMR4 {
 
class MmReadoutElement : public MuonReadoutElement {
 
   public:
    
    struct parameterBook {
      double halfThickness{0.};
      double halfShortWidth{0.};
      double halfLongWidth{0.};
      double halfHeight{0.};
      unsigned int nGasGaps{0};
      std::vector<int> readoutSide{};
      std::vector<StripLayerPtr> layers{};
 
#ifndef SIMULATIONBASE
        std::shared_ptr<Acts::SurfaceBoundFactory> layerBounds{};
#endif
 
 
    };
 
    struct defineArgs : public MuonReadoutElement::defineArgs,
                        public parameterBook {};
 
    MmReadoutElement(defineArgs&& args);
 
    const parameterBook& getParameters() const;
    ActsTrk::DetectorType detectorType() const override final {
        return ActsTrk::DetectorType::Mm;
    }
 
    
    double thickness() const override final;
    
    int multilayer() const;
    double moduleHeight() const;
    double moduleWidthS() const;
    double moduleWidthL() const;
    double moduleThickness() const;
    double gapLengthS(const IdentifierHash& layerHash) const;
    double gapLengthL(const IdentifierHash& layerHash) const;
    double gapHeight(const IdentifierHash& layerHash) const;
    unsigned int nGasGaps() const;
    unsigned int numStrips(const IdentifierHash& layerHash) const;
    unsigned int firstStrip(const IdentifierHash& layerHash) const;
    double stripLength(const IdentifierHash& measHash) const;
    int readoutSide(const IdentifierHash& measHash) const;
 
    StatusCode initElement() override final;

 
    // measurementHash : Creates a channelHash based on channel's and gasGap's number, which are
    // retrieved through the MmIdHelper by inputting the calculated channel Identifier.
    IdentifierHash measurementHash(const Identifier& measId) const override final;
    // layerHash : Creates a layerHash based solely on gasGap's number, which is
    // retrieved through the MmIdHelper by inputting the calculated channel Identifier.
    IdentifierHash layerHash(const Identifier& measId) const override final;
    // layerHash : Creates a layerHash based on the first 4 bits of the channelHash,
    // which retrieve the layer.
    IdentifierHash layerHash(const IdentifierHash& measHash) const;
    // measurementId : Retrieves the channel Identifier based on the above created channelHash (measurementHash)
    Identifier measurementId(const IdentifierHash& measHash) const override final;
    
    static IdentifierHash createHash(const int gasGap, const int strip);
    Amg::Vector3D stripPosition(const ActsTrk::GeometryContext& ctx, const Identifier& measId) const;
    Amg::Vector3D stripPosition(const ActsTrk::GeometryContext& ctx, const IdentifierHash& measHash) const;
    Amg::Vector3D leftStripEdge(const ActsTrk::GeometryContext& ctx, const Identifier& measId) const;
    Amg::Vector3D leftStripEdge(const ActsTrk::GeometryContext& ctx, const IdentifierHash& measHash) const; 
    Amg::Vector3D rightStripEdge(const ActsTrk::GeometryContext& ctx, const Identifier& measId) const;
    Amg::Vector3D rightStripEdge(const ActsTrk::GeometryContext& ctx, const IdentifierHash& measHash) const;
 
    const StripLayer& stripLayer(const Identifier& measId) const;    
    const StripLayer& stripLayer(const IdentifierHash& measHash) const;    
 
    friend ActsTrk::TransformCacheDetEle<MmReadoutElement>;
 
   private:
       
    
 
    static unsigned int gasGapNumber(const IdentifierHash& measHash);
    static unsigned int stripNumber(const IdentifierHash& measHash);
 
    parameterBook m_pars{};        
 
    const MmIdHelper& m_idHelper{idHelperSvc()->mmIdHelper()};
 
    const int m_multilayer{m_idHelper.multilayer(identify())};
 
    Amg::Transform3D fromGapToChamOrigin(const IdentifierHash& layerHash) const;
 
};
std::ostream& operator<<(std::ostream& ostr, const MmReadoutElement::parameterBook& pars);
}  // namespace MuonGMR4
 
namespace ActsTrk{
    template <> Amg::Transform3D 
        TransformCacheDetEle<MuonGMR4::MmReadoutElement>::fetchTransform(const DetectorAlignStore* store) const;
    template <> Identifier
        TransformCacheDetEle<MuonGMR4::MmReadoutElement>::identify() const;
}
 
 
#include <MuonReadoutGeometryR4/MmReadoutElement.icc>
#endif