xMuonHitAnalysis.h

Go to the documentation of this file.

 
/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef HitAnalysis_xMuonSimHitAnalysis_H
#define HitAnalysis_xMuonSimHitAnalysis_H
 
#include "AthenaBaseComps/AthHistogramAlgorithm.h"
#include "StoreGate/ReadHandleKey.h"
 
#include "xAODMuonSimHit/MuonSimHitContainer.h"
#include "MuonReadoutGeometryR4/MuonDetectorManager.h"
#include "ActsGeometryInterfaces/GeometryContext.h"
 
namespace MuonValR4{
    class xMuonHitAnalysis : public AthHistogramAlgorithm{
        public:
            using AthHistogramAlgorithm::AthHistogramAlgorithm;
 
            virtual StatusCode initialize() override final;
            virtual StatusCode execute(const EventContext& ctx) override final;
        private:
            SG::ReadHandleKey<xAOD::MuonSimHitContainer> m_inputKey{this, "InputKey", ""};
            // ACTS geometry context
            SG::ReadHandleKey<ActsTrk::GeometryContext> m_geoCtxKey{this, "AlignmentKey", "ActsAlignment", "cond handle key"};
            ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc{this, "MuonIdHelperSvc", 
                                                                "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

            const MuonGMR4::MuonDetectorManager* m_detMgr{nullptr};
 
            Gaudi::Property<std::string> m_path{this, "HistPath", "MuonSDOAnalysis"};
 
            using StIdx_t = Muon::MuonStationIndex::StIndex;
            using TechIdx_t = Muon::MuonStationIndex::TechnologyIndex;
 
            Gaudi::Property<int> m_techIdx{this, "techIndex", -1};
 
            struct HistoSet{
                HistoSet(const Muon::IMuonIdHelperSvc* idHelperSvc,
                         const ServiceHandle<ITHistSvc>& histSvc,
                         const std::string& basePath,
                         const TechIdx_t techIdx,
                         const int stIdent);

                int identifier{-1};
                TH1* h_dirTheta{nullptr};
                TH1* h_dirPhi{nullptr};
                TH1* h_energy{nullptr};
                TH1* h_deposit{nullptr};
                TH1* h_pdgId{nullptr};
                TH1* h_stationPhi{nullptr};
                TH1* h_stationEta{nullptr};
                TH1* h_localHitPos{nullptr};
                TH1* h_localHitZ{nullptr};
                TH1* h_globHitXY{nullptr};
                TH1* h_globHitRZ{nullptr};
            };
            std::vector<HistoSet> m_histos{};
    };
}
 
#endif