DigitEffiData.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#include "MuonCondData/DigitEffiData.h"
#include "GeoModelKernel/throwExcept.h"
 
namespace Muon{
    DigitEffiData::DigitEffiData(const Muon::IMuonIdHelperSvc* idHelperSvc, double defaultEffi):
        AthMessaging{"DigitEffiData"},
        m_idHelperSvc{idHelperSvc},
        m_defaultEffi{defaultEffi}{}
 
    Identifier DigitEffiData::getLookUpId(const Identifier& channelId, bool isInnerQ1 /*=false*/) const {
        Identifier lookUpId{};
        using TechIndex = Muon::MuonStationIndex::TechnologyIndex;
        switch (m_idHelperSvc->technologyIndex(channelId)){
            case TechIndex::MDT:
                lookUpId = channelId;
                break;
            case TechIndex::RPC:
            case TechIndex::TGC:
            case TechIndex::CSCI:
                lookUpId = m_idHelperSvc->gasGapId(channelId);
                break;
            case TechIndex::MM:
                lookUpId = m_idHelperSvc->mmIdHelper().febID(channelId);
                break;
            case TechIndex::STGC:
                lookUpId = m_idHelperSvc->stgcIdHelper().hvID(channelId, isInnerQ1);
                break;
            default:
                THROW_EXCEPTION("Invalid muon technology");
        };
        return lookUpId;
    }
    double DigitEffiData::getEfficiency(const Identifier& channelId, bool isInnerQ1 /*=false*/) const {
        EffiMap::const_iterator effi_itr = m_effiData.find(getLookUpId(channelId, isInnerQ1));
        if (effi_itr != m_effiData.end()) {
            ATH_MSG_VERBOSE("Channel "<<m_idHelperSvc->toString(channelId) <<" has an efficiency of "<<effi_itr->second);
            return effi_itr->second;
        }
        ATH_MSG_VERBOSE("Efficiency of channel "<<m_idHelperSvc->toString(channelId)<<" is unknown. Return 1.");
        return m_defaultEffi;
    }
 
    StatusCode DigitEffiData::setEfficiency(const Identifier& channelId, const double effi, bool isInnerQ1 /*=false*/){
        const Identifier gasGapId = getLookUpId(channelId, isInnerQ1);
        auto insert_itr = m_effiData.insert(std::make_pair(gasGapId, effi));
        if (!insert_itr.second) {
            ATH_MSG_ERROR("An efficiency for gasGap "<<m_idHelperSvc->toString(gasGapId)
            <<" has already been stored "<<m_effiData[gasGapId]<<" vs. "<<effi);
            return StatusCode::FAILURE;
        }
        return StatusCode::SUCCESS;
    }
}