MdtTubeLayer.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONREADOUTGEOMETRYR4_MDTTUBELAYER_H
#define MUONREADOUTGEOMETRYR4_MDTTUBELAYER_H
 
#include <MuonReadoutGeometryR4/MuonDetectorDefs.h>
#include <GeoModelKernel/GeoVPhysVol.h>
#include <GeoModelKernel/GeoTransform.h>
#include <GeoModelKernel/GeoVolumeCursor.h>
#include <GeoModelUtilities/TransientConstSharedPtr.h>
 
#include <set>
namespace MuonGMR4{
    class MdtTubeLayer;
    using MdtTubeLayerPtr = GeoModel::TransientConstSharedPtr<MdtTubeLayer>;
    struct  MdtTubeLayerSorter{
        bool operator()(const MdtTubeLayer& a, const MdtTubeLayer& b) const;
        bool operator()(const MdtTubeLayerPtr&a, const MdtTubeLayerPtr& b) const;
    };
    using MdtTubeLayerSet = std::set<MdtTubeLayerPtr, MdtTubeLayerSorter>;

    class MdtTubeLayer{
    public:
        
        struct CutTubes {
            unsigned int firstTube{0};
            unsigned int lastTube{0};
            double unCutHalfLength{0.};
            bool operator<(const CutTubes& other) const {
                return lastTube <= other.firstTube;
            }
        };
        using CutTubeSet = std::set<CutTubes, std::less<>>;

        struct TubePositioner{
            unsigned int firstTube{0};
            unsigned int lastTube{0};
            PVConstLink tubeVol{};
            Amg::Vector3D firstTubePos{Amg::Vector3D::Zero()};
            bool operator<(const TubePositioner& other) const {
                return lastTube < other.firstTube;
            }
        };
        using TubePositionerSet = std::set<TubePositioner, std::less<>>;
        
        friend MdtTubeLayerSorter;
        MdtTubeLayer(const PVConstLink layer, 
                     const GeoIntrusivePtr<const GeoTransform> toLayTrf,
                     const CutTubeSet& cutTubes);
        unsigned int nTubes() const;
        const Amg::Transform3D& layerTransform() const;
        const Amg::Transform3D tubeTransform(const unsigned int tube) const;
        const Amg::Vector3D tubePosInLayer(const unsigned int tube) const;
        double tubeHalfLength(const unsigned int tube) const;
        double uncutHalfLength(const unsigned int tube) const;
        PVConstLink getTubeNode(unsigned int tube) const;
 
 
    private:
        PVConstLink m_layerNode{nullptr};
        GeoIntrusivePtr<const GeoTransform> m_layTrf{nullptr};
        CutTubeSet m_cutTubes{};
        TubePositionerSet m_tubePos{};
        Amg::Vector3D m_tubePitch{Amg::Vector3D::Zero()};
    };
 
    bool operator<(const MdtTubeLayer::CutTubes& a, const unsigned int tube);
    bool operator<(const unsigned int tube, const MdtTubeLayer::CutTubes& a);
 
    bool operator<(const MdtTubeLayer::TubePositioner&a,const unsigned int tube);
    bool operator<(const unsigned int tube, const MdtTubeLayer::TubePositioner& a);
 
}
#endif