MuonSensitiveDetector.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONSENSITIVEDETECTORSR4_MUONSENSITIVEDETECTOR_H
#define MUONSENSITIVEDETECTORSR4_MUONSENSITIVEDETECTOR_H
 
#include <GeoPrimitives/GeoPrimitives.h>
 
#include <GaudiKernel/EventContext.h>
#include <StoreGate/ReadHandleKey.h>
#include <MuonReadoutGeometryR4/MuonDetectorManager.h>
#include <xAODMuonSimHit/MuonSimHitContainer.h>
#include <xAODMuonSimHit/MuonSimHitAuxContainer.h>
#include <HitManagement/AthenaHitsVector.h>
#include <HitManagement/HitCollectionMap.h>
#include <AthenaBaseComps/AthMessaging.h>
 
#include <G4VSensitiveDetector.hh>
 
class AtlasG4EventUserInfo;
class G4Track;
class HepMcParticleLink;

namespace MuonG4R4 {
    using MuonSimHitsVec = AthHitVec::AuxStoreHitCollection<xAOD::MuonSimHitContainer,
                                                            xAOD::MuonSimHitAuxContainer>;
 
    class MuonSensitiveDetector : public G4VSensitiveDetector, public AthMessaging {
        public:
            MuonSensitiveDetector(const std::string& name, 
                                  const std::string& output_key,
                                  const std::string& trf_storeKey,
                                  const MuonGMR4::MuonDetectorManager* detMgr);
            ~MuonSensitiveDetector() = default;

            virtual void Initialize(G4HCofThisEvent* HCE) override final;
 
        private:
            std::string m_writeKey{};
            xAOD::MuonSimHitContainer* m_outContainer{};
            AtlasG4EventUserInfo* m_g4UserEventInfo{};
            SG::ReadHandleKey<ActsTrk::DetectorAlignStore> m_trfCacheKey;

            AtlasG4EventUserInfo* eventInfo() const;

            xAOD::MuonSimHit* lastSnapShot(const Identifier& gasGapId,
                                           const G4Step* hitStep,
                                           const HepMcParticleLink& particleLink);
 
        protected:
            bool processStep(const G4Step* step) const;

            const EventContext& eventContext() const;

            ActsTrk::GeometryContext getGeoContext() const;

            HepMcParticleLink genParticleLink(const G4Track* track) const;

            xAOD::MuonSimHit* lastSnapShot(const Identifier& gasGapId,
                                           const G4Step* hitStep);

            xAOD::MuonSimHit* propagateAndSaveStrip(const Identifier& hitId,
                                                    const Amg::Transform3D& toGasGap,
                                                    const G4Step* hitStep);

            xAOD::MuonSimHit* saveHit(const Identifier& hitId,
                                      const Amg::Vector3D& hitPos,
                                      const Amg::Vector3D& hitDir,
                                      const double globTime,
                                      const G4Step* hitStep);

            const MuonGMR4::MuonDetectorManager* m_detMgr{nullptr};
 
            
    };
}
 
#endif