MuonChamberToolTest.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 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"
#include "Acts/Geometry/TrackingVolume.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::vector<Amg::Vector3D> cornerPoints(const ActsTrk::GeometryContext& gctx, const Acts::Volume& volume) const;
        std::vector<Amg::Vector3D> cornerPoints(const ActsTrk::GeometryContext& gctx, const Acts::Surface& surface) const;
        
        void saveEnvelope(const ActsTrk::GeometryContext& gctx,
                          const std::string& envName,
                          const Acts::Volume& envelopeVol,
                          const std::vector<const Acts::Surface*>& assocSurfaces,
                          const std::vector<const Acts::Volume*>& subVolumes={}) const;
 
        StatusCode checkPortals(const ActsTrk::GeometryContext& gctx,
                                const Acts::TrackingVolume& volume) const;
                                
        StatusCode checkChambers(const ActsTrk::GeometryContext& gctx) const;
        StatusCode checkEnvelopes(const ActsTrk::GeometryContext& gctx) const;
        StatusCode checkTrackingGeometry(const ActsTrk::GeometryContext& gctx, 
                                         const Acts::TrackingGeometry& trackingGeometry) const;
        StatusCode allReadoutInEnvelope(const ActsTrk::GeometryContext& ctx,
                                        const Chamber& envelope) const; 
        
        StatusCode pointInside(const ActsTrk::GeometryContext& gctx,
                                const Acts::Volume& boundVol,
                                const Amg::Vector3D& point,
                                const std::string& descr,
                                const Identifier& channelId) const;           
        bool hasOverlap(const ActsTrk::GeometryContext& gctx,
                        const std::vector<Amg::Vector3D>& chamberEdges,
                        const Acts::Volume& volume) const;
        
        StatusCode testReadoutSurfaces(const ActsTrk::GeometryContext& gctx,
                                       const MuonReadoutElement& element,
                                       const Acts::Volume& boundVol) const;

        StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                  const MdtReadoutElement& readOutEle,
                                  const Acts::Volume& boundVol) const;
        StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                  const RpcReadoutElement& readOutEle,
                                  const Acts::Volume& boundVol) const;
        StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                  const TgcReadoutElement& readOutEle,
                                  const Acts::Volume& boundVol) const;
        StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                  const sTgcReadoutElement& readOutEle,
                                  const Acts::Volume& boundVol) const;
        StatusCode testReadoutEle(const ActsTrk::GeometryContext& gctx,
                                  const MmReadoutElement& readOutEle,
                                  const Acts::Volume& boundVol) const;
 
 
 
        ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "IdHelperSvc", 
                                                "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
 
        ActsTrk::GeoContextReadKey_t m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};
 
        ServiceHandle<ActsTrk::ITrackingGeometrySvc> m_trackingGeometrySvc{this, "TrackingGeometrySvc","ActsTrackingGeometrySvc"};
        Gaudi::Property<unsigned> m_overlapSamples{this, "overlapSamples", 50};
        Gaudi::Property<std::string> m_overlapChambObj{this, "chamberOverlapFile", "OverlapingChambers.obj"};
        Gaudi::Property<bool> m_ignoreOverlapCh{this, "ignoreChamberOverlap", true};
        Gaudi::Property<bool> m_ignoreOutsideSurf{this, "ignoreOutsideSurface", true};
        Gaudi::Property<bool> m_dumpObjs{this, "dumpVolumes" , false};
         const MuonDetectorManager* m_detMgr{nullptr};
 
};
}
#endif