MuonChamberToolTest.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONGEOMODELTESTR4_MUONCHAMBERTOOLTEST_H
#define MUONGEOMODELTESTR4_MUONCHAMBERTOOLTEST_H
 
#include <AthenaBaseComps/AthReentrantAlgorithm.h>
 
#include <MuonReadoutGeometryR4/MuonDetectorManager.h>
#include <MuonReadoutGeometryR4/Chamber.h>
#include <MuonReadoutGeometryR4/MdtReadoutElement.h>
#include <MuonReadoutGeometryR4/RpcReadoutElement.h>
#include <MuonReadoutGeometryR4/TgcReadoutElement.h>
#include <MuonReadoutGeometryR4/sTgcReadoutElement.h>
#include <MuonReadoutGeometryR4/MmReadoutElement.h>
 
#include <ActsGeometryInterfaces/GeometryContext.h>
#include <ActsGeometryInterfaces/ITrackingGeometrySvc.h>
#include <StoreGate/ReadCondHandleKey.h>
 
 
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Surfaces/StrawSurface.hpp"
#include "Acts/Surfaces/PlaneSurface.hpp"
 
namespace MuonGMR4 { 
 
 
class MuonChamberToolTest: public AthReentrantAlgorithm {
    public:
        using AthReentrantAlgorithm::AthReentrantAlgorithm;
        ~MuonChamberToolTest() = default;
 
        StatusCode execute(const EventContext& ctx) const override;        
        StatusCode initialize() override;        
 
        bool isReEntrant() const override final {return false;}   
    
    private:
        std::array<Amg::Vector3D, 8> cornerPoints(const ActsTrk::GeometryContext& gctx, const Acts::Volume& volume) const;
        std::array<Amg::Vector3D, 8> cornerPoints(const ActsTrk::GeometryContext& gctx, const Acts::StrawSurface& surface) const;
        std::array<Amg::Vector3D, 4> cornerPoints(const ActsTrk::GeometryContext& gctx, const Acts::PlaneSurface&) const;
        
        
        void saveEnvelope(const ActsTrk::GeometryContext& gctx,
                          const std::string& envName,
                          const Acts::Volume& envelopeVol,
                          const std::vector<const MuonGMR4::MuonReadoutElement*>& assocRE,
                          const std::vector<std::shared_ptr<Acts::Volume>>& subVolumes={}) const;
        StatusCode checkChambers(const ActsTrk::GeometryContext& gctx) const;
        StatusCode checkEnvelopes(const ActsTrk::GeometryContext& gctx) const;
        template <class EnvelopeType>
          StatusCode allReadoutInEnvelope(const ActsTrk::GeometryContext& ctx,
                                          const EnvelopeType& envelope) const; 
        
        template <class EnvelopeType>        
            StatusCode pointInside(const ActsTrk::GeometryContext& gctx,
                                   const EnvelopeType& envelope,
                                   const Acts::Volume& boundVol,
                                   const Amg::Vector3D& point,
                                   const std::string& descr,
                                   const Identifier& channelId) const;
        StatusCode pointInside(const ActsTrk::GeometryContext& gctx,
                               const Acts::TrackingVolume& volume,
                               const Amg::Vector3D& point,
                               const std::string& descr,
                               const Identifier& chamberId) const;           
        bool hasOverlap(const ActsTrk::GeometryContext& gctx,
                        const std::array<Amg::Vector3D, 8>& chamberEdges,
                        const Acts::Volume& volume) const;
        template <class EnvelopeType>
          StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                    const MdtReadoutElement& readOutEle,
                                    const EnvelopeType& envelope,
                                    const Acts::Volume& boundVol) const;
        template <class EnvelopeType>
          StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                    const RpcReadoutElement& readOutEle,
                                    const EnvelopeType& envelope,
                                    const Acts::Volume& boundVol) const;
        template <class EnvelopeType>
          StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                    const TgcReadoutElement& readOutEle,
                                    const EnvelopeType& envelope,
                                    const Acts::Volume& boundVol) const;
        template <class EnvelopeType>
          StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                    const sTgcReadoutElement& readOutEle,
                                    const EnvelopeType& envelope,
                                    const Acts::Volume& boundVol) const;
        template <class EnvelopeType>
          StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                    const MmReadoutElement& readOutEle,
                                    const EnvelopeType& envelope,
                                    const Acts::Volume& boundVol) const;
 
 
 
        ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "IdHelperSvc", 
                                                "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
 
        SG::ReadHandleKey<ActsTrk::GeometryContext> m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};
 
        ServiceHandle<ActsTrk::ITrackingGeometrySvc> m_trackingGeometrySvc{this, "TrackingGeometrySvc","ActsTrackingGeometrySvc"};
        Gaudi::Property<unsigned> m_overlapSamples{this, "overlapSamples", 100};
        Gaudi::Property<std::string> m_overlapChambObj{this, "chamberOverlapFile", "OverlapingChambers.obj"};
        Gaudi::Property<bool> m_ignoreOverlapCh{this, "ignoreChamberOverlap", true};
        Gaudi::Property<bool> m_dumpObjs{this, "dumpVolumes" , false};
         const MuonDetectorManager* m_detMgr{nullptr};
 
};
}
#endif