RegionSelectionSvc Class Reference

Seolect calibration region. More...

#include <RegionSelectionSvc.h>

Inheritance diagram for RegionSelectionSvc:

Collaboration diagram for RegionSelectionSvc:

Public Member Functions
	RegionSelectionSvc (const std::string &name, ISvcLocator *svc_locator)
	contructor
virtual	~RegionSelectionSvc ()
	destructor
virtual StatusCode	initialize () override
	initialize
void	Print (std::ostream &os) const
	debug print
bool	isInRegion (const MuonCalib::MuonFixedId &id) const
	return true if id is in selected region
const std::string &	GetRegionSelection () const
	get region selection
const std::vector< MuonCalib ::NtupleStationId > &	GetStationsInRegions () const
	get stations/multilayers in region
int	AddRegionNtuples (const std::string &infile, TChain *chain)
	add region ntuples which are relevant for this region to TChain
int	AddRegionNtuples (const std::string &infile, TChain *chain, std::list< std::string > &dirnames)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Private Member Functions
void	search_chambers_in_region ()
bool	ProcessString (const std::string &input)
	process string
void	print_list_of_selected_chambers () const
	print a list of selected chambers - debug function
void	convert_to_tower (int station, int phi, int eta, std::vector< int > &tower) const
	convert station, eta, phi in tower vector

Private Attributes
Gaudi::Property< std::string >	m_region_string {this, "Region", ""}
	region string - job option
Gaudi::Property< bool >	m_print_list_of_selected_chambers {this, "PrintList", false}
	print list of selected chambers if set to true - job option
std::unique_ptr< MuonCalib ::RegionSelectorBase >	m_master_region
	master region
std::vector< MuonCalib::NtupleStationId >	m_stations_in_region
	chambers/mutlilayers in calibraition region
std::set< MuonCalib::NtupleStationId >	m_unique_chambers
const MuonGM::MuonDetectorManager *	m_detMgr {nullptr}
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
	search for chambers and multilayers in selected region towers in selected region
ToolHandle< MuonCalib::IIdToFixedIdTool >	m_idToFixedIdTool {this, "IdToFixedIdTool", "MuonCalib::IdToFixedIdTool"}

Detailed Description

Seolect calibration region.

Author

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

Definition at line 35 of file RegionSelectionSvc.h.

Constructor & Destructor Documentation

RegionSelectionSvc()

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

contructor

Definition at line 29 of file RegionSelectionSvc.cxx.

: AthService(name, svc_locator) {}

~RegionSelectionSvc()

RegionSelectionSvc::~RegionSelectionSvc ( )

virtualdefault

destructor

Member Function Documentation

AddRegionNtuples() [1/2]

int RegionSelectionSvc::AddRegionNtuples	(	const std::string &	infile,
		TChain *	chain )

add region ntuples which are relevant for this region to TChain

Definition at line 55 of file RegionSelectionSvc.cxx.

                                                                                {
    std::list<std::string> dirnames;
    return AddRegionNtuples(infile, chain, dirnames);
}

AddRegionNtuples() [2/2]

int RegionSelectionSvc::AddRegionNtuples	(	const std::string &	infile,
		TChain *	chain,
		std::list< std::string > &	dirnames )

Definition at line 59 of file RegionSelectionSvc.cxx.

                                                                                                              {
    int n_trees(0);
    if (dirnames.empty()) {
        // open file
        std::unique_ptr<TFile> inf{TFile::Open(infile.c_str(), "READ")};
        TIter nextkey(inf->GetListOfKeys());
        TKey *key;
        // take a list here because erase does not invalidate iterators
        while ((key = (TKey *)nextkey())) { dirnames.push_back(key->GetName()); }
    }
    // loop over found ntuples and selected region
    for (std::set<MuonCalib::NtupleStationId>::const_iterator it = m_unique_chambers.begin(); it != m_unique_chambers.end(); ++it)
        for (std::list<std::string>::iterator nt_it = dirnames.begin(); nt_it != dirnames.end(); ++nt_it) {
            // truncate region name
            std::string regname(it->regionId(), 0, nt_it->size());
            if (regname == (*nt_it)) {
                chain->AddFile(infile.c_str(), TChain::kBigNumber, ((*nt_it) + "/" + "Segments").c_str());
                std::list<std::string>::iterator new_it = nt_it;
                --new_it;
                dirnames.erase(nt_it);
                nt_it = new_it;
                ++n_trees;
            }
        }
    return n_trees;
}

convert_to_tower()

void RegionSelectionSvc::convert_to_tower ( int station, int phi, int eta, std::vector< int > & tower ) const

private

convert station, eta, phi in tower vector

GetRegionSelection()

const std::string & RegionSelectionSvc::GetRegionSelection

(

)

const

get region selection

Definition at line 53 of file RegionSelectionSvc.cxx.

{ return m_region_string; }

GetStationsInRegions()

const std::vector< MuonCalib::NtupleStationId > & RegionSelectionSvc::GetStationsInRegions

(

)

const

get stations/multilayers in region

Definition at line 54 of file RegionSelectionSvc.cxx.

{ return m_stations_in_region; }

initialize()

StatusCode RegionSelectionSvc::initialize ( )

overridevirtual

initialize

Definition at line 31 of file RegionSelectionSvc.cxx.

                                           {
    ATH_MSG_INFO("initialize RegionSelectionSvc");
 
    if (!ProcessString(m_region_string)) return StatusCode::FAILURE;
 
    // detector stre - id to fixed id ...
    SmartIF<StoreGateSvc> detStore{service("DetectorStore")};
    ATH_CHECK(detStore.isValid());
    ATH_MSG_DEBUG("Retrieved DetectorStore");
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(detStore->retrieve(m_detMgr));
    ATH_CHECK(m_idToFixedIdTool.retrieve());
    ATH_MSG_INFO("Retrieved " << m_idToFixedIdTool);
 
    search_chambers_in_region();
    if (m_print_list_of_selected_chambers) print_list_of_selected_chambers();
 
    return StatusCode::SUCCESS;
}

isInRegion()

bool RegionSelectionSvc::isInRegion

(

const MuonCalib::MuonFixedId &

id

)

const

return true if id is in selected region

Definition at line 51 of file RegionSelectionSvc.cxx.

{ return m_master_region->Result(id); }

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

Print()

void RegionSelectionSvc::Print

(

std::ostream &

os

)

const

debug print

Definition at line 86 of file RegionSelectionSvc.cxx.

                                                    {
    if (!m_master_region) {
        os << "ERROR";
        return;
    }
    m_master_region->Print(os);
}

print_list_of_selected_chambers()

void RegionSelectionSvc::print_list_of_selected_chambers ( ) const

private

print a list of selected chambers - debug function

Definition at line 137 of file RegionSelectionSvc.cxx.

                                                                {
    ATH_MSG_INFO("Selected regions: " << m_stations_in_region.size());
    msg(MSG::INFO) << "Selected regions: ";
    for (std::vector<MuonCalib::NtupleStationId>::const_iterator it = m_stations_in_region.begin(); it != m_stations_in_region.end();
         ++it) {
        msg(MSG::INFO) << " " << it->regionId();
    }
    msg(MSG::INFO) << endmsg;
}

ProcessString()

bool RegionSelectionSvc::ProcessString ( const std::string & input )

private

process string

Definition at line 132 of file RegionSelectionSvc.cxx.

                                                              {
    m_master_region = RegionSelectorBase::GetRegion(input);
    return (m_master_region != nullptr);
}

search_chambers_in_region()

void RegionSelectionSvc::search_chambers_in_region ( )

private

Definition at line 94 of file RegionSelectionSvc.cxx.

                                                    {
    MdtIdHelper::const_id_iterator it = m_idHelperSvc->mdtIdHelper().module_begin();
    MdtIdHelper::const_id_iterator it_end = m_idHelperSvc->mdtIdHelper().module_end();
    for (; it != it_end; ++it) {
        const MuonGM::MdtReadoutElement *detEl = m_detMgr->getMdtReadoutElement(m_idHelperSvc->mdtIdHelper().channelID(*it, 1, 1, 1));
        if (!detEl) continue;
        // get number of mls;
        int n_mls = m_idHelperSvc->mdtIdHelper().numberOfMultilayers(*it);
        // fixed id
        MuonCalib::MuonFixedId fixed_id(m_idToFixedIdTool->idToFixedId(*it));
        std::vector<MuonCalib::NtupleStationId> the_ids;
        for (int i = 1; i <= n_mls; i++) {
            fixed_id.setMdtMultilayer(i);
            if (isInRegion(fixed_id)) {
                NtupleStationId the_id(fixed_id);
                the_ids.push_back(the_id);
            }
        }
        // all multilayers are in the region - inser the whole station
        if (static_cast<int>(the_ids.size()) == n_mls) {
            NtupleStationId full_station(fixed_id);
            full_station.SetMultilayer(0);
            m_stations_in_region.push_back(full_station);
        } else {
            for (std::vector<MuonCalib::NtupleStationId>::iterator it2 = the_ids.begin(); it2 != the_ids.end(); ++it2) {
                m_stations_in_region.push_back(*it2);
            }
        }
    }
    // create a set of chambers (not multilayers)
    std::vector<int> tower(3);
    for (std::vector<MuonCalib::NtupleStationId>::iterator it2 = m_stations_in_region.begin(); it2 != m_stations_in_region.end(); ++it2) {
        MuonCalib::NtupleStationId station(*it2);
        station.SetMultilayer(0);
        m_unique_chambers.insert(station);
    }
}

Member Data Documentation

m_detMgr

const MuonGM::MuonDetectorManager* RegionSelectionSvc::m_detMgr {nullptr}

private

Definition at line 65 of file RegionSelectionSvc.h.

{nullptr};  

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> RegionSelectionSvc::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

search for chambers and multilayers in selected region towers in selected region

Definition at line 67 of file RegionSelectionSvc.h.

{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

m_idToFixedIdTool

ToolHandle<MuonCalib::IIdToFixedIdTool> RegionSelectionSvc::m_idToFixedIdTool {this, "IdToFixedIdTool", "MuonCalib::IdToFixedIdTool"}

private

Definition at line 68 of file RegionSelectionSvc.h.

{this, "IdToFixedIdTool", "MuonCalib::IdToFixedIdTool"};

m_master_region

std::unique_ptr<MuonCalib ::RegionSelectorBase> RegionSelectionSvc::m_master_region

private

master region

Definition at line 61 of file RegionSelectionSvc.h.

m_print_list_of_selected_chambers

Gaudi::Property<bool> RegionSelectionSvc::m_print_list_of_selected_chambers {this, "PrintList", false}

private

print list of selected chambers if set to true - job option

Definition at line 59 of file RegionSelectionSvc.h.

{this, "PrintList", false};

m_region_string

Gaudi::Property<std::string> RegionSelectionSvc::m_region_string {this, "Region", ""}

private

region string - job option

Definition at line 57 of file RegionSelectionSvc.h.

{this, "Region", ""};

m_stations_in_region

std::vector<MuonCalib::NtupleStationId> RegionSelectionSvc::m_stations_in_region

private

chambers/mutlilayers in calibraition region

Definition at line 63 of file RegionSelectionSvc.h.

m_unique_chambers

std::set<MuonCalib::NtupleStationId> RegionSelectionSvc::m_unique_chambers

private

Definition at line 64 of file RegionSelectionSvc.h.

---

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

• RegionSelectionSvc.h
• RegionSelectionSvc.cxx