#include <RegionSelectionSvc.h>

Inheritance diagram for RegionSelectionSvc:

Collaboration diagram for RegionSelectionSvc:

Public Member Functions
	RegionSelectionSvc (const std::string &name, ISvcLocator *svc_locator)
	contructor More...

virtual	~RegionSelectionSvc ()
	destructor More...

virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
	just some crazy atheta function More...

StatusCode	initialize ()
	initialize More...

void	Print (std::ostream &os) const
	debug print More...

bool	isInRegion (const MuonCalib::MuonFixedId &id) const
	return true if id is in selected region More...

std::string	GetRegionSelection () const
	get region selection More...

const std::vector< MuonCalib ::NtupleStationId > &	GetStationsInRegions () const
	get stations/multilayers in region More...

int	AddRegionNtuples (const std::string &infile, TChain *chain)
	add region ntuples which are relevant for this region to TChain More...

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

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
void	search_chambers_in_region ()

bool	ProcessString (const std::string &input)
	process string More...

void	print_list_of_selected_chambers () const
	print a list of selected chambers - debug function More...

void	convert_to_tower (int station, int phi, int eta, std::vector< int > &tower) const
	convert station, eta, phi in tower vector More...

Private Attributes
Gaudi::Property< std::string >	m_region_string {this, "Region", ""}
	region string - job option More...

Gaudi::Property< bool >	m_print_list_of_selected_chambers {this, "PrintList", false}
	print list of selected chambers if set to true - job option More...

std::unique_ptr< MuonCalib ::RegionSelectorBase >	m_master_region
	master region More...

std::vector< MuonCalib::NtupleStationId >	m_stations_in_region
	chambers/mutlilayers in calibraition region More...

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 More...

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.

29 : 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 67 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 71 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()

std::string RegionSelectionSvc::GetRegionSelection ( ) const

get region selection

Definition at line 65 of file RegionSelectionSvc.cxx.

65 { return m_region_string; }

◆ GetStationsInRegions()

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

get stations/multilayers in region

Definition at line 66 of file RegionSelectionSvc.cxx.

66 { return m_stations_in_region; }

◆ initialize()

StatusCode RegionSelectionSvc::initialize ( )

initialize

Definition at line 43 of file RegionSelectionSvc.cxx.

                                            {
     ATH_MSG_INFO("initialize RegionSelectionSvc");
  
     if (!ProcessString(m_region_string)) return StatusCode::FAILURE;
  
     // detector stre - id to fixed id ...
     StoreGateSvc *detStore{nullptr};
     ATH_CHECK(serviceLocator()->service("DetectorStore", detStore));
     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;
 }

◆ interfaceID()

static const InterfaceID& RegionSelectionSvc::interfaceID ( )

inlinestatic

Definition at line 41 of file RegionSelectionSvc.h.

                                             {
         static const InterfaceID IID_IRegionSelectionSvc("RegionSelectionSvc", 1, 0);
         return IID_IRegionSelectionSvc;
     }

◆ isInRegion()

bool RegionSelectionSvc::isInRegion ( const MuonCalib::MuonFixedId & id ) const

return true if id is in selected region

Definition at line 63 of file RegionSelectionSvc.cxx.

63 { return m_master_region->Result(id); }

◆ 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);
   }

◆ Print()

void RegionSelectionSvc::Print ( std::ostream & os ) const

debug print

Definition at line 98 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 149 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 144 of file RegionSelectionSvc.cxx.

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

◆ queryInterface()

StatusCode RegionSelectionSvc::queryInterface	(	const InterfaceID &	riid,
		void **	ppvUnknown
	)

virtual

just some crazy atheta function

Definition at line 31 of file RegionSelectionSvc.cxx.

                                                                                         {
     ATH_MSG_VERBOSE("StatusCode RegionSelectionSvc::queryInterface");
  
     if (interfaceID().versionMatch(riid)) {
         *ppvUnknown = (RegionSelectionSvc *)this;
     } else {
         return AthService::queryInterface(riid, ppvUnknown);
     }
  
     return StatusCode::SUCCESS;
 }

◆ search_chambers_in_region()

void RegionSelectionSvc::search_chambers_in_region ( )

private

Definition at line 106 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 71 of file RegionSelectionSvc.h.

◆ 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 73 of file RegionSelectionSvc.h.

◆ m_idToFixedIdTool

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

private

Definition at line 74 of file RegionSelectionSvc.h.

◆ m_master_region

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

private

master region

Definition at line 67 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 65 of file RegionSelectionSvc.h.

◆ m_region_string

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

private

region string - job option

Definition at line 63 of file RegionSelectionSvc.h.

◆ m_stations_in_region

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

private

chambers/mutlilayers in calibraition region

Definition at line 69 of file RegionSelectionSvc.h.

◆ m_unique_chambers

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

private

Definition at line 70 of file RegionSelectionSvc.h.

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

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ RegionSelectionSvc()

◆ ~RegionSelectionSvc()

Member Function Documentation

◆ AddRegionNtuples() [1/2]

◆ AddRegionNtuples() [2/2]

◆ convert_to_tower()

◆ GetRegionSelection()

◆ GetStationsInRegions()

◆ initialize()

◆ interfaceID()

◆ isInRegion()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ Print()

◆ print_list_of_selected_chambers()

◆ ProcessString()

◆ queryInterface()

◆ search_chambers_in_region()

Member Data Documentation

◆ m_detMgr

◆ m_idHelperSvc

◆ m_idToFixedIdTool

◆ m_master_region

◆ m_print_list_of_selected_chambers

◆ m_region_string

◆ m_stations_in_region

◆ m_unique_chambers