MuonPhaseII/MuonDetDescr/MuonReadoutGeometryR4/MuonReadoutGeometryR4/MdtReadoutElement.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_MDTREADOUTELEMENT_H
#define MUONREADOUTGEOMETRYR4_MDTREADOUTELEMENT_H
 
#include <MuonReadoutGeometryR4/MuonReadoutElement.h>
#include <MuonReadoutGeometryR4/MdtTubeLayer.h>
 
#ifndef SIMULATIONBASE
#   include "Acts/Surfaces/TrapezoidBounds.hpp"
#endif
namespace MuonGMR4 {
 
class MdtReadoutElement : public MuonReadoutElement {
 
   public:
    
    struct parameterBook {
        // the number of tubes of the first layer
        std::vector<MdtTubeLayerPtr> tubeLayers{};
        
        std::set<IdentifierHash> removedTubes{};
 
        double tubeWall{0.};
        double tubeInnerRad{0.};
        double tubePitch{0.};
        double wireTension{0.};
        double endPlugLength{0.};
        double deadLength{0.};
        double radLengthX0{0.};
        double shortHalfX{0.};
        double longHalfX{0.};
        double halfY{0.};
        double halfHeight{0.};
        double readoutSide{1.};
#ifndef SIMULATIONBASE
        ActsTrk::SurfaceBoundSetPtr<Acts::LineBounds> tubeBounds{};
        ActsTrk::SurfaceBoundSetPtr<Acts::TrapezoidBounds> layerBounds{};
#endif
 
    };
 
    struct defineArgs : public MuonReadoutElement::defineArgs,
                        public parameterBook {};
 
    MdtReadoutElement(defineArgs&& args);
 
    const parameterBook& getParameters() const;
    ActsTrk::DetectorType detectorType() const override final {
        return ActsTrk::DetectorType::Mdt;
    }
    double thickness() const override final;
 
    StatusCode initElement() override final;
    unsigned int multilayer() const;
    unsigned int numLayers() const;
    unsigned int numTubesInLay() const;
 
    static unsigned int tubeNumber(const IdentifierHash& hash);
    static  unsigned int layerNumber(const IdentifierHash& hash);
    static IdentifierHash measurementHash(unsigned int layerNumber, unsigned int tubeNumber);
 
 
    IdentifierHash measurementHash(const Identifier& measId) const override final;
    
    IdentifierHash layerHash(const Identifier& measId) const override final;
    static IdentifierHash layerHash(const IdentifierHash& measHash);
    Identifier measurementId(const IdentifierHash& measHash) const override final;
 
    bool isValid(const IdentifierHash& measHash) const;
 
    bool isBarrel() const;
    
    double tubePitch() const;
    double innerTubeRadius() const;
    double tubeRadius() const;
    double moduleWidthS() const;
    double moduleWidthL() const;
    double moduleHeight() const;
    double moduleThickness() const;
 
    
    Amg::Vector3D globalTubePos(const ActsGeometryContext& ctx, 
                                const Identifier& measId) const;
    
    Amg::Vector3D globalTubePos(const ActsGeometryContext& ctx,
                                const IdentifierHash& hash) const;
    
    Amg::Vector3D readOutPos(const ActsGeometryContext& ctx,
                             const Identifier& measId) const;
 
    Amg::Vector3D readOutPos(const ActsGeometryContext& ctx,
                             const IdentifierHash& measId) const;
    Amg::Vector3D highVoltPos(const ActsGeometryContext& ctx,
                              const Identifier& measId) const;
    Amg::Vector3D highVoltPos(const ActsGeometryContext& ctx,
                              const IdentifierHash& measId) const;
    double distanceToReadout(const ActsGeometryContext& ctx,
                             const Identifier& measId,
                             const Amg::Vector3D& globPoint) const;
    double distanceToReadout(const ActsGeometryContext& ctx,
                             const IdentifierHash& measHash,
                             const Amg::Vector3D& globPoint) const;
    
    double activeTubeLength(const IdentifierHash& hash) const;
        
    double tubeLength(const IdentifierHash& hash) const;
    
    double wireLength(const IdentifierHash& hash) const;
 
#ifndef SIMULATIONBASE
    std::map<Identifier, std::shared_ptr<Acts::Surface>> getSurfaces() const override final;
#endif
 
        friend ActsTrk::TransformCacheDetEle<MdtReadoutElement>;
   private:
        Amg::Vector3D localTubePos(const IdentifierHash& hash) const;
        Amg::Transform3D toChamberLayer(const IdentifierHash& hash) const;
        Amg::Transform3D toTubeFrame(const IdentifierHash& hash) const;
 
 
        parameterBook m_pars{};
        const MdtIdHelper& m_idHelper{idHelperSvc()->mdtIdHelper()};
        int m_stML{m_idHelper.multilayer(identify())};
        bool m_isBarrel{m_idHelper.isBarrel(identify())};
};
 
std::ostream& operator<<(std::ostream& ostr, const MdtReadoutElement::parameterBook& pars);
}  // namespace MuonGMR4
 
namespace ActsTrk{
    template <> Amg::Transform3D 
        TransformCacheDetEle<MuonGMR4::MdtReadoutElement>::fetchTransform(const DetectorAlignStore* store) const;
}
 
#include <MuonReadoutGeometryR4/MdtReadoutElement.icc>
#endif