MdtCalibInputSvc Class Reference

#include <MdtCalibInputSvc.h>

Inheritance diagram for MdtCalibInputSvc:

Collaboration diagram for MdtCalibInputSvc:

Public Member Functions
	MdtCalibInputSvc (const std::string &name, ISvcLocator *svc_locator)
	Service constructor. More...
virtual	~MdtCalibInputSvc ()
	destructor More...
virtual StatusCode	initialize () override
	service initalizer - reads files More...
const MuonCalib::MdtStationT0Container *	GetT0 (const MuonCalib::NtupleStationId &id) const
	Get t0 container for Station. More...
const MuonCalib::IRtRelation *	GetRtRelation (const MuonCalib::NtupleStationId &id) const
	Get rt relation container. More...
const MuonCalib::BFieldCorFunc *	GetBCorr (const MuonCalib::NtupleStationId &id)
	Get B field correction. More...
const MuonCalib::IRtResolution *	GetResolution (const MuonCalib::NtupleStationId &id) const
	Get Resolution. More...
const MuonCalib::IRtRelation *	GetRtRelation () const
	Get rt-Relation for calibration region. More...
const MuonCalib::BFieldCorFunc *	GetBCorr ()
	Get b-field correction for calibratino region. More...
const MuonCalib::IRtResolution *	GetResolution () const
	Get resolution for calibration region. More...
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Private Member Functions
bool	create_b_field_correction (const MuonCalib::NtupleStationId &id)
	create the b-field correction More...
const MuonCalib::BFieldCorFunc *	findbfieldfun (const MuonCalib::NtupleStationId &id)
void	create_mean_rts ()
	create mean rt relations, and resolutions for the selected calibration region More...
StatusCode	read_calib_input ()

Private Attributes
ToolHandle< MuonCalib::ICalibrationIOTool >	m_calib_input_tool {this, "CalibrationInputTool", "MuonCalib::CalibrationDummyIOTool"}
	calibration io tool to be used More...
std::map< MuonCalib::NtupleStationId, MuonCalib::MdtStationT0Container * >	m_t0
	calibration data sorted by station id More...
std::map< MuonCalib::NtupleStationId, MuonCalib::IRtRelation * >	m_rt_relation
std::map< MuonCalib::NtupleStationId, MuonCalib::BFieldCorFunc * >	m_B_corr
std::map< MuonCalib::NtupleStationId, MuonCalib::IRtResolution * >	m_spat_res
ServiceHandle< RegionSelectionSvc >	m_reg_sel_svc
	pointer to region selection service More...
const MuonCalib::IRtRelation *	m_sel_region_rt {nullptr}
	rt relation - resolution - and correction function for the selected region - is average of all matching rt relations More...
const MuonCalib::BFieldCorFunc *	m_sel_region_b {nullptr}
const MuonCalib::IRtResolution *	m_sel_region_res {nullptr}
MuonCalib::NtupleStationId	m_mean_station_id
	station id for mean rt More...
std::set< MuonCalib::NtupleStationId > m_t0_warned	ATLAS_THREAD_SAFE
	give warnings about missing calibration only once per chamber More...
std::set< MuonCalib::NtupleStationId > m_rt_warned	ATLAS_THREAD_SAFE

Detailed Description

Athena service which read calibration from text files and sorts them by station id

Author

rausc.nosp@m.her@.nosp@m.cern..nosp@m.ch

Definition at line 30 of file MdtCalibInputSvc.h.

Constructor & Destructor Documentation

MdtCalibInputSvc()

MdtCalibInputSvc::MdtCalibInputSvc	(	const std::string &	name,
		ISvcLocator *	svc_locator
	)

Service constructor.

Definition at line 22 of file MdtCalibInputSvc.cxx.

                                                                                  :
    AthService(name, svc_locator), m_reg_sel_svc("RegionSelectionSvc", name) {
    declareProperty("RegionSelectionSvc", m_reg_sel_svc);
}

~MdtCalibInputSvc()

MdtCalibInputSvc::~MdtCalibInputSvc ( )

virtual

destructor

Definition at line 27 of file MdtCalibInputSvc.cxx.

                                    {
    std::map<MuonCalib::NtupleStationId, MuonCalib::MdtStationT0Container *>::iterator it1;
    for (it1 = m_t0.begin(); it1 != m_t0.end(); ++it1) delete it1->second;
    std::map<MuonCalib::NtupleStationId, MuonCalib::IRtRelation *>::iterator it2;
    for (it2 = m_rt_relation.begin(); it2 != m_rt_relation.end(); ++it2) delete it2->second;
    std::map<MuonCalib::NtupleStationId, MuonCalib::BFieldCorFunc *>::iterator it3;
    for (it3 = m_B_corr.begin(); it3 != m_B_corr.end(); ++it3) delete it3->second;
    std::map<MuonCalib::NtupleStationId, MuonCalib::IRtResolution *>::iterator it4;
    for (it4 = m_spat_res.begin(); it4 != m_spat_res.end(); ++it4) delete it4->second;
}

Member Function Documentation

create_b_field_correction()

bool MdtCalibInputSvc::create_b_field_correction ( const MuonCalib::NtupleStationId &  id )

private

create the b-field correction

Definition at line 114 of file MdtCalibInputSvc.cxx.

                                                                                   {
    std::map<MuonCalib::NtupleStationId, MuonCalib::IRtRelation *>::const_iterator it(m_rt_relation.find(id));
    if (it == m_rt_relation.end()) return false;
    ATH_MSG_INFO("Initiailizing B-Field correction for " << id.regionId());
    const MuonCalib::IRtRelation *rt_rel(it->second);
    // create magnetic field correction
    std::vector<double> corr_params(2);
    corr_params[0] = 3080.0;  // high voltage
    corr_params[1] = 0.11;    // epsilon parameter
    m_B_corr[id] = new MuonCalib::BFieldCorFunc(corr_params, rt_rel);
    return true;
}

create_mean_rts()

void MdtCalibInputSvc::create_mean_rts ( )

private

create mean rt relations, and resolutions for the selected calibration region

Definition at line 127 of file MdtCalibInputSvc.cxx.

                                       {
    ATH_MSG_INFO("MdtCalibInputSvc::create_mean_rts()");
    std::list<const MuonCalib::IRtRelation *> matching_relations;
    std::list<MuonCalib::NtupleStationId> matching_ids;
    // loop over all rts
    for (std::map<MuonCalib::NtupleStationId, MuonCalib::IRtRelation *>::const_iterator it = m_rt_relation.begin();
         it != m_rt_relation.end(); ++it) {
        MuonCalib::MuonFixedId id;
        // if the rt relation is stored per chamber, check both multilayers
        if (it->first.GetMl() == 0) {
            id.set_mdt();
            for (int i = 1; i <= 2; i++) {
                id.setStationName(it->first.GetStation());
                id.setStationEta(it->first.GetEta());
                id.setStationPhi(it->first.GetPhi());
                id.setMdtMultilayer(i);
                if (m_reg_sel_svc->isInRegion(id)) {
                    matching_relations.push_back(it->second);
                    matching_ids.push_back(it->first);
                    break;
                }
            }
        } else {
            id.setStationName(it->first.GetStation());
            id.setStationEta(it->first.GetEta());
            id.setStationPhi(it->first.GetPhi());
            id.setMdtMultilayer(it->first.GetMl());
            if (m_reg_sel_svc->isInRegion(id)) {
                matching_relations.push_back(it->second);
                matching_ids.push_back(it->first);
            }
        }
    }
 
    ATH_MSG_INFO("Found " << matching_relations.size() << " rt-relations for calibration region");
    if (matching_relations.empty()) return;
    // averageing over rt relations is not yet implemented - take the first found
    if (matching_relations.size() > 1) { ATH_MSG_WARNING("More than one rt relation for this region loaded! Taking first!"); }
    m_sel_region_rt = *(matching_relations.begin());
    m_mean_station_id = *(matching_ids.begin());
    m_sel_region_res = GetResolution(*(matching_ids.begin()));
}  // end MdtCalibInputSvc::create_mean_rts

findbfieldfun()

const MuonCalib::BFieldCorFunc * MdtCalibInputSvc::findbfieldfun ( const MuonCalib::NtupleStationId &  id )

private

Definition at line 170 of file MdtCalibInputSvc.cxx.

                                                                                                {
    std::map<MuonCalib::NtupleStationId, MuonCalib::BFieldCorFunc *>::const_iterator it(m_B_corr.find(id));
    if (it == m_B_corr.end()) {
        if (!create_b_field_correction(id)) return nullptr;
        return findbfieldfun(id);
    }
    return it->second;
}

GetBCorr() [1/2]

const MuonCalib::BFieldCorFunc* MdtCalibInputSvc::GetBCorr ( )

inline

Get b-field correction for calibratino region.

Definition at line 56 of file MdtCalibInputSvc.h.

                                             {
        if (!m_sel_region_b) m_sel_region_b = GetBCorr(m_mean_station_id);
        return m_sel_region_b;
    }

GetBCorr() [2/2]

const MuonCalib::BFieldCorFunc * MdtCalibInputSvc::GetBCorr

(

const MuonCalib::NtupleStationId &

id

)

Get B field correction.

Definition at line 83 of file MdtCalibInputSvc.cxx.

                                                                                           {
    const MuonCalib::BFieldCorFunc *fun;
    if ((fun = findbfieldfun(id)) != nullptr) return fun;
    MuonCalib::NtupleStationId chamber_id(id);
    chamber_id.SetMultilayer(0);
    return findbfieldfun(chamber_id);
}

GetResolution() [1/2]

const MuonCalib::IRtResolution* MdtCalibInputSvc::GetResolution ( ) const

inline

Get resolution for calibration region.

Definition at line 61 of file MdtCalibInputSvc.h.

{ return m_sel_region_res; }

GetResolution() [2/2]

const MuonCalib::IRtResolution * MdtCalibInputSvc::GetResolution

(

const MuonCalib::NtupleStationId &

id

)

const

Get Resolution.

Definition at line 91 of file MdtCalibInputSvc.cxx.

                                                                                                      {
    std::map<MuonCalib::NtupleStationId, MuonCalib::IRtResolution *>::const_iterator it;
    if ((it = m_spat_res.find(id)) == m_spat_res.end()) {
        MuonCalib::NtupleStationId chamber_id = id;
        chamber_id.SetMultilayer(0);
        if ((it = m_spat_res.find(chamber_id)) == m_spat_res.end()) {
            if (m_rt_warned.find(chamber_id) == m_rt_warned.end()) {
                ATH_MSG_FATAL("Rt relation not loaded for station");
                m_rt_warned.insert(chamber_id);
            }
            return nullptr;
        }
    }
    return it->second;
}

GetRtRelation() [1/2]

const MuonCalib::IRtRelation* MdtCalibInputSvc::GetRtRelation ( ) const

inline

Get rt-Relation for calibration region.

Definition at line 54 of file MdtCalibInputSvc.h.

{ return m_sel_region_rt; }

GetRtRelation() [2/2]

const MuonCalib::IRtRelation * MdtCalibInputSvc::GetRtRelation

(

const MuonCalib::NtupleStationId &

id

)

const

Get rt relation container.

Definition at line 65 of file MdtCalibInputSvc.cxx.

                                                                                                    {
    static std::mutex warn_mutex;
    std::map<MuonCalib::NtupleStationId, MuonCalib::IRtRelation *>::const_iterator it;
    if ((it = m_rt_relation.find(id)) == m_rt_relation.end()) {
        MuonCalib::NtupleStationId chamber_id = id;
        chamber_id.SetMultilayer(0);
        if ((it = m_rt_relation.find(chamber_id)) == m_rt_relation.end()) {
            std::lock_guard<std::mutex> guard(warn_mutex);
            if (m_rt_warned.find(chamber_id) == m_rt_warned.end()) {
                ATH_MSG_WARNING("Rt relation not loaded for station" << chamber_id.regionId());
                m_rt_warned.insert(chamber_id);
            }
            return nullptr;
        }
    }
    return it->second;
}

GetT0()

const MuonCalib::MdtStationT0Container * MdtCalibInputSvc::GetT0

(

const MuonCalib::NtupleStationId &

id

)

const

Get t0 container for Station.

Definition at line 47 of file MdtCalibInputSvc.cxx.

                                                                                                      {
    static std::mutex warn_mutex;
    std::map<MuonCalib::NtupleStationId, MuonCalib::MdtStationT0Container *>::const_iterator it;
    if ((it = m_t0.find(id)) == m_t0.end()) {
        MuonCalib::NtupleStationId chamber_id = id;
        chamber_id.SetMultilayer(0);
        if ((it = m_t0.find(chamber_id)) == m_t0.end()) {
            std::lock_guard<std::mutex> guard(warn_mutex);
            if (m_t0_warned.find(chamber_id) == m_t0_warned.end()) {
                ATH_MSG_WARNING("T0 not loaded for station " << id.regionId());
                m_t0_warned.insert(chamber_id);
            }
            return nullptr;
        }
    }
    return it->second;
}

initialize()

StatusCode MdtCalibInputSvc::initialize ( )

overridevirtual

service initalizer - reads files

Definition at line 38 of file MdtCalibInputSvc.cxx.

                                        {
    ATH_MSG_INFO("initialize()");
    // get region selection service
    ATH_CHECK(m_reg_sel_svc.retrieve());
    ATH_CHECK(m_calib_input_tool.retrieve());
    ATH_CHECK(read_calib_input());
    return StatusCode::SUCCESS;
}

interfaceID()

static const InterfaceID& MdtCalibInputSvc::interfaceID ( )

inlinestatic

Definition at line 37 of file MdtCalibInputSvc.h.

                                            {
        // interface to enable retrieving of a pointer to the singleton //
        static const InterfaceID IID_IMdtCalibInputSvc("MdtCalibInputSvc", 1, 0);
        return IID_IMdtCalibInputSvc;
    }

msg() [1/2]

MsgStream& AthCommonMsg< Service >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< Service >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

bool AthCommonMsg< Service >::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

read_calib_input()

StatusCode MdtCalibInputSvc::read_calib_input ( )

private

Definition at line 107 of file MdtCalibInputSvc.cxx.

                                              {
    ATH_CHECK(m_calib_input_tool->LoadT0(m_t0, -1));
    ATH_CHECK(m_calib_input_tool->LoadRt(m_rt_relation, m_spat_res, -1));
    create_mean_rts();
    return StatusCode::SUCCESS;
}

Member Data Documentation

ATLAS_THREAD_SAFE [1/2]

std::set<MuonCalib::NtupleStationId> m_t0_warned MdtCalibInputSvc::ATLAS_THREAD_SAFE

mutableprivate

give warnings about missing calibration only once per chamber

Definition at line 86 of file MdtCalibInputSvc.h.

ATLAS_THREAD_SAFE [2/2]

std::set<MuonCalib::NtupleStationId> m_rt_warned MdtCalibInputSvc::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 87 of file MdtCalibInputSvc.h.

m_B_corr

std::map<MuonCalib::NtupleStationId, MuonCalib::BFieldCorFunc *> MdtCalibInputSvc::m_B_corr

private

Definition at line 69 of file MdtCalibInputSvc.h.

m_calib_input_tool

ToolHandle<MuonCalib::ICalibrationIOTool> MdtCalibInputSvc::m_calib_input_tool {this, "CalibrationInputTool", "MuonCalib::CalibrationDummyIOTool"}

private

calibration io tool to be used

Definition at line 65 of file MdtCalibInputSvc.h.

m_mean_station_id

MuonCalib::NtupleStationId MdtCalibInputSvc::m_mean_station_id

private

station id for mean rt

Definition at line 84 of file MdtCalibInputSvc.h.

m_reg_sel_svc

ServiceHandle<RegionSelectionSvc> MdtCalibInputSvc::m_reg_sel_svc

private

pointer to region selection service

Definition at line 72 of file MdtCalibInputSvc.h.

m_rt_relation

std::map<MuonCalib::NtupleStationId, MuonCalib::IRtRelation *> MdtCalibInputSvc::m_rt_relation

private

Definition at line 68 of file MdtCalibInputSvc.h.

m_sel_region_b

const MuonCalib::BFieldCorFunc* MdtCalibInputSvc::m_sel_region_b {nullptr}

private

Definition at line 81 of file MdtCalibInputSvc.h.

m_sel_region_res

const MuonCalib::IRtResolution* MdtCalibInputSvc::m_sel_region_res {nullptr}

private

Definition at line 82 of file MdtCalibInputSvc.h.

m_sel_region_rt

const MuonCalib::IRtRelation* MdtCalibInputSvc::m_sel_region_rt {nullptr}

private

rt relation - resolution - and correction function for the selected region - is average of all matching rt relations

Definition at line 80 of file MdtCalibInputSvc.h.

m_spat_res

std::map<MuonCalib::NtupleStationId, MuonCalib::IRtResolution *> MdtCalibInputSvc::m_spat_res

private

Definition at line 70 of file MdtCalibInputSvc.h.

m_t0

std::map<MuonCalib::NtupleStationId, MuonCalib::MdtStationT0Container *> MdtCalibInputSvc::m_t0

private

calibration data sorted by station id

Definition at line 67 of file MdtCalibInputSvc.h.

---

The documentation for this class was generated from the following files:

• MdtCalibInputSvc.h
• MdtCalibInputSvc.cxx