RpcFastDigiTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#ifndef RPC_DIGITIZATIONR4_RPCFASTDIGITOOL_H
#define RPC_DIGITIZATIONR4_RPCFASTDIGITOOL_H
 
 
#include "MuonDigitizationR4/MuonDigitizationTool.h"
#include "MuonDigitContainer/RpcDigitContainer.h"
#include "MuonCondData/DigitEffiData.h"
#include "MuonReadoutGeometryR4/RpcReadoutElement.h"
 
#include "GaudiKernel/PhysicalConstants.h"
 
namespace MuonR4{
    class RpcFastDigiTool final: public MuonDigitizationTool {
        public:
            RpcFastDigiTool(const std::string& type, const std::string& name, const IInterface* pIID);
 
            StatusCode initialize() override final;
            StatusCode finalize() override final;
        protected:
            StatusCode digitize(const EventContext& ctx,
                                const TimedHits& hitsToDigit,
                                xAOD::MuonSimHitContainer* sdoContainer) const override final; 
        
 
        private:
            using EdgeSide = MuonGMR4::RpcReadoutElement::EdgeSide;
            int m_stIdxBIL{-1}; // Station name index of the BIL stations
        int m_stIdxBIS{-1}; // Station name index of the BIS stations
            
            bool digitizeHit(const Identifier& gasGapId,
                             const bool measuresPhi,
                             const MuonGMR4::RpcReadoutElement& reEle,
                             const double hitTime,
                             const Amg::Vector2D& locPos,
                             const Muon::DigitEffiData* effiMap,
                             RpcDigitCollection& outContainer,
                             CLHEP::HepRandomEngine* rndEngine,
                             DeadTimeMap& deadTimes) const;
            
            bool digitizeHitBI(const Identifier& gasGapId,
                               const MuonGMR4::RpcReadoutElement& reEle,
                               const double hitTime,
                               const Amg::Vector2D& locPos,
                               const Muon::DigitEffiData* effiMap,
                               RpcDigitCollection& outContainer,
                               CLHEP::HepRandomEngine* rndEngine,
                               DeadTimeMap& deadTimes) const;
 
            static double timeOverThreshold(CLHEP::HepRandomEngine* rndmEngine) ;
 
            using DigiCache = OutDigitCache_t<RpcDigitCollection>;
            SG::WriteHandleKey<RpcDigitContainer> m_writeKey{this, "OutputObjectName", "RPC_DIGITS"};
 
            SG::ReadCondHandleKey<Muon::DigitEffiData> m_effiDataKey{this, "EffiDataKey", "RpcDigitEff",
                                                                    "Efficiency constants of the individual Rpc gasGaps"};
 
            mutable std::array<std::atomic<unsigned>, 2> m_allHits ATLAS_THREAD_SAFE{};
            mutable std::array<std::atomic<unsigned>, 2> m_acceptedHits ATLAS_THREAD_SAFE{};
 
            Gaudi::Property<double> m_propagationVelocity{this, "propSpeed", 0.5 * Gaudi::Units::c_light,
                                                         "Propagation speed of the signal inside the strip"}; // in mm/ns
    
            Gaudi::Property<double> m_stripTimeResolution{this, "timeResolution", 0.6 * Gaudi::Units::nanosecond,
                                                          "Estimated time resolution of the strip readout"};
 
            Gaudi::Property<double> m_deadTime{this, "deadTime", 100.*Gaudi::Units::nanosecond};
 
            Gaudi::Property<bool> m_digitizeMuonOnly{this, "ProcessTrueMuonsOnly", false, 
                                                     "If set to true hit with pdgId != 13 are skipped"};
 
    };
}
#endif