SpectrometerSector.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_SPECTROMETERSECTOR_H
#define MUONREADOUTGEOMETRYR4_SPECTROMETERSECTOR_H
 
#include "AthenaBaseComps/AthMessaging.h"
#ifndef SIMULATIONBASE
 
#include <MuonReadoutGeometryR4/MuonReadoutElement.h>
#include <MuonReadoutGeometryR4/Chamber.h>
#include <ActsGeometryInterfaces/GeometryDefs.h>
#include <Acts/Geometry/TrapezoidVolumeBounds.hpp>
 
namespace Acts {
    class Volume;
    class PlaneSurface;
}
namespace MuonGMR4{
    class SpectrometerSector;
    class Chamber;
}
 
namespace MuonGMR4 {
    class SpectrometerSector : public AthMessaging {   
        public:
            using ChamberPtr = GeoModel::TransientConstSharedPtr<Chamber>;
            using ChamberSet = std::vector<ChamberPtr>;
 
            struct chamberLocation{
                public:
                    using enum Acts::TrapezoidVolumeBounds::BoundValues;
                    using BoundPtr_t = std::shared_ptr<const Acts::TrapezoidVolumeBounds>;
                    chamberLocation(const Amg::Vector3D& origin,
                                    const MuonReadoutElement* reEle,
                                    BoundPtr_t bounds):
                        m_origin{origin}, m_reEle{reEle}, m_bounds{std::move(bounds)}{}
                    chamberLocation(const chamberLocation& other) = default;
                    chamberLocation(chamberLocation&& other) = default;
                    chamberLocation& operator=(const chamberLocation& other) = default;
                    chamberLocation& operator=(chamberLocation&& other) = default;
                    bool insideYZ(const Amg::Vector3D& pos) const {
                      return minY() <= pos.y() &&  maxY() >= pos.y() &&
                             minZ() <= pos.z() &&  maxZ() >= pos.z();
                    }
                     double minY() const { return m_origin.y() - m_bounds->get(eHalfLengthY); }
                     double maxY() const { return m_origin.y() + m_bounds->get(eHalfLengthY); }
                     double minZ() const { return m_origin.z() - m_bounds->get(eHalfLengthZ); }
                     double maxZ() const { return m_origin.z() + m_bounds->get(eHalfLengthZ); }
                     double width(const double y0) const {
                        const double tanPhiHalf = 0.5*(m_bounds->get(eHalfLengthXposY) - m_bounds->get(eHalfLengthXnegY)) 
                                                / m_bounds->get(eHalfLengthY); 
                        return m_bounds->get(eHalfLengthXnegY) + tanPhiHalf * (y0 - minY());
                     }
                     const MuonReadoutElement* readoutEle() const { return m_reEle; }
                    const BoundPtr_t& bounds() const { return m_bounds; }
                    const Amg::Vector3D& location() const { return m_origin; }
                private:
                    Amg::Vector3D m_origin{Amg::Vector3D::Zero()};
                    const MuonReadoutElement* m_reEle{nullptr};
                    BoundPtr_t m_bounds{};
            };
 
            struct defineArgs{
                ChamberSet chambers{};
                std::shared_ptr<Acts::TrapezoidVolumeBounds> bounds{};
                std::shared_ptr<const Acts::PlaneSurface> surface{};
                
                std::vector<chamberLocation> detectorLocs{}; 
            };
 
            SpectrometerSector(defineArgs&& args);
            SpectrometerSector(const SpectrometerSector& other) = delete;
            SpectrometerSector& operator=(const SpectrometerSector& other) = delete;
 
            ~SpectrometerSector() = default;
 
            bool operator<(const SpectrometerSector& other) const;
 
            std::string identString() const;
            const Muon::IMuonIdHelperSvc* idHelperSvc() const;
            Muon::MuonStationIndex::ChIndex chamberIndex() const;
            int stationPhi() const;
            int sector() const;
            int8_t side() const;
            bool barrel() const;
            const Amg::Transform3D& localToGlobalTrans(const ActsGeometryContext& gctx) const;
            Amg::Transform3D globalToLocalTrans(const ActsGeometryContext& gctx) const;
            const Acts::PlaneSurface& surface() const;
            const ChamberSet& chambers() const;
            double halfXLong() const;
            double halfXShort() const;
            double halfY() const;
            double halfZ() const;  
            const defineArgs& parameters() const;
            std::shared_ptr<Acts::Volume> boundingVolume(const ActsGeometryContext& gctx) const;
            std::shared_ptr<Acts::TrapezoidVolumeBounds> bounds() const;
            Chamber::ReadoutSet readoutEles() const;
            const std::vector<chamberLocation> & chamberLocations() const; 
            const std::vector<unsigned int>& logicalLayerIdx(const MuonReadoutElement* reEle) const;
 
        private:
            defineArgs m_args{};
 
            std::unordered_map<const MuonReadoutElement*, std::vector<unsigned int>> fillDetLayIdCache() const;
            const std::unordered_map<const MuonReadoutElement*, std::vector<unsigned int>> m_detLayIdCache{fillDetLayIdCache()};
            Identifier computeDetLayerId(const MuonReadoutElement* rele) const;
            unsigned int nLayerPerReadout (const MuonReadoutElement* rele) const;
    };
    
    std::ostream& operator<<(std::ostream& ostr,
                             const SpectrometerSector::defineArgs& args);
 
    std::ostream& operator<<(std::ostream& ostr,
                             const SpectrometerSector& chamber);
 
}
 
#endif
#endif