MMRawDataMonAlg Class Reference

#include <MMRawDataMonAlg.h>

Inheritance diagram for MMRawDataMonAlg:

Collaboration diagram for MMRawDataMonAlg:

Public Types
enum class	Environment_t {   user = 0 , online , tier0 , tier0Raw ,   tier0ESD , AOD , altprod }
	Specifies the processing environment. More...
enum class	DataType_t {   userDefined = 0 , monteCarlo , collisions , cosmics ,   heavyIonCollisions }
	Specifies what type of input data is being monitored. More...

Public Member Functions
	MMRawDataMonAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~MMRawDataMonAlg ()=default
virtual StatusCode	initialize () override
	initialize
virtual StatusCode	fillHistograms (const EventContext &ctx) const override
	adds event to the monitoring histograms
virtual StatusCode	execute (const EventContext &ctx) const override
	Applies filters and trigger requirements.
void	fill (const ToolHandle< GenericMonitoringTool > &groupHandle, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
	Fills a vector of variables to a group by reference.
void	fill (const ToolHandle< GenericMonitoringTool > &groupHandle, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
	Fills a vector of variables to a group by reference.
template<typename... T>
void	fill (const ToolHandle< GenericMonitoringTool > &groupHandle, T &&... variables) const
	Fills a variadic list of variables to a group by reference.
void	fill (const std::string &groupName, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
	Fills a vector of variables to a group by name.
void	fill (const std::string &groupName, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
	Fills a vector of variables to a group by name.
template<typename... T>
void	fill (const std::string &groupName, T &&... variables) const
	Fills a variadic list of variables to a group by name.
Environment_t	environment () const
	Accessor functions for the environment.
Environment_t	envStringToEnum (const std::string &str) const
	Convert the environment string from the python configuration to an enum object.
DataType_t	dataType () const
	Accessor functions for the data type.
DataType_t	dataTypeStringToEnum (const std::string &str) const
	Convert the data type string from the python configuration to an enum object.
const ToolHandle< GenericMonitoringTool > &	getGroup (const std::string &name) const
	Get a specific monitoring tool from the tool handle array.
const ToolHandle< Trig::TrigDecisionTool > &	getTrigDecisionTool () const
	Get the trigger decision tool member.
bool	trigChainsArePassed (const std::vector< std::string > &vTrigNames) const
	Check whether triggers are passed.
SG::ReadHandle< xAOD::EventInfo >	GetEventInfo (const EventContext &) const
	Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
virtual float	lbAverageInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate the average mu, i.e.
virtual float	lbInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate instantaneous number of interactions, i.e.
virtual float	lbAverageLuminosity (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate average luminosity (in ub-1 s-1 => 10^30 cm-2 s-1).
virtual float	lbLuminosityPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate the instantaneous luminosity per bunch crossing.
virtual double	lbDuration (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate the duration of the luminosity block (in seconds)
virtual float	lbAverageLivefraction (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate the average luminosity livefraction.
virtual float	livefractionPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate the live fraction per bunch crossing ID.
virtual double	lbLumiWeight (const EventContext &ctx=Gaudi::Hive::currentContext()) const
	Calculate the average integrated luminosity multiplied by the live fraction.
virtual StatusCode	parseList (const std::string &line, std::vector< std::string > &result) const
	Parse a string into a vector.
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Protected Attributes
ToolHandleArray< GenericMonitoringTool >	m_tools {this,"GMTools",{}}
	Array of Generic Monitoring Tools.
PublicToolHandle< Trig::TrigDecisionTool >	m_trigDecTool
	Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterTool >	m_DQFilterTools {this,"FilterTools",{}}
	Array of Data Quality filter tools.
SG::ReadCondHandleKey< LuminosityCondData >	m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
SG::ReadCondHandleKey< LBDurationCondData >	m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
SG::ReadCondHandleKey< TrigLiveFractionCondData >	m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
AthMonitorAlgorithm::Environment_t	m_environment
	Instance of the Environment_t enum.
AthMonitorAlgorithm::DataType_t	m_dataType
	Instance of the DataType_t enum.
Gaudi::Property< std::string >	m_environmentStr {this,"Environment","user"}
	Environment string pulled from the job option and converted to enum.
Gaudi::Property< std::string >	m_dataTypeStr {this,"DataType","userDefined"}
	DataType string pulled from the job option and converted to enum.
Gaudi::Property< std::string >	m_triggerChainString {this,"TriggerChain",""}
	Trigger chain string pulled from the job option and parsed into a vector.
std::vector< std::string >	m_vTrigChainNames
	Vector of trigger chain names parsed from trigger chain string.
Gaudi::Property< std::string >	m_fileKey {this,"FileKey",""}
	Internal Athena name for file.
Gaudi::Property< bool >	m_useLumi {this,"EnableLumi",false}
	Allows use of various luminosity functions.
Gaudi::Property< float >	m_defaultLBDuration {this,"DefaultLBDuration",60.}
	Default duration of one lumi block.
Gaudi::Property< int >	m_detailLevel {this,"DetailLevel",0}
	Sets the level of detail used in the monitoring.
SG::ReadHandleKey< xAOD::EventInfo >	m_EventInfoKey {this,"EventInfoKey","EventInfo"}
	Key for retrieving EventInfo from StoreGate.

Private Types
typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > >	MonVarVec_t
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	fillMMOverviewVects (const Muon::MMPrepData *, MMOverviewHistogramStruct &vects, MMByPhiStruct(&occupancyPlots)[16][2]) const
void	fillMMOverviewHistograms (const MMOverviewHistogramStruct &vects, MMByPhiStruct(&occupancyPlots)[16][2], const int lb) const
StatusCode	fillMMSummaryVects (const Muon::MMPrepData *, MMSummaryHistogramStruct(&vects)[2][16][2][2][4]) const
StatusCode	fillMMHistograms (const Muon::MMPrepData *) const
StatusCode	fillMMSummaryHistograms (const MMSummaryHistogramStruct(&vects)[2][16][2][2][4]) const
void	clusterFromTrack (const xAOD::TrackParticleContainer *, const int lb) const
void	clusterFromSegments (const Trk::SegmentCollection *, const int lb) const
int	get_PCB_from_channel (const int channel) const
int	get_FEB_from_channel (const int channel, const int stEta) const
int	get_sectorPhi_from_stationPhi_stName (const int stationPhi, const std::string &stName) const
int	get_sectorEta_from_stationEta (const int stationEta) const
int	get_bin_for_occ_CSide_hist (const int stationEta, const int multiplet, const int gas_gap) const
int	get_bin_for_occ_ASide_hist (const int stationEta, const int multiplet, const int gas_gap) const
int	get_bin_for_occ (const int gas_gap, const int PCB) const
int	get_bin_for_feb_occ (const int gas_gap, const int FEB) const
int	get_bin_for_occ_CSide_pcb_eta2_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB) const
int	get_bin_for_occ_CSide_pcb_eta1_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB) const
int	get_bin_for_occ_ASide_pcb_eta2_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB) const
int	get_bin_for_occ_ASide_pcb_eta1_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB) const
int	get_bin_for_occ_lb_CSide_pcb_eta2_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB, const int isector) const
int	get_bin_for_occ_lb_CSide_pcb_eta1_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB, int isector) const
int	get_bin_for_occ_lb_ASide_pcb_eta1_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB, int isector) const
int	get_bin_for_occ_lb_ASide_pcb_eta2_hist (const int stationEta, const int multiplet, const int gas_gap, const int PCB, const int isector) const
int	get_bin_for_occ_lb_pcb_hist (const int multiplet, const int gas_gap, const int PCB) const
void	MMEfficiency (const xAOD::TrackParticleContainer *) const
void	fillMMTrigger (const xAOD::NSWMMTPRDOContainer *, const int) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
ToolHandle< CP::IMuonSelectionTool >	m_muonSelectionTool {this,"MuonSelectionTool","CP::MuonSelectionTool/MuonSelectionTool"}
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager >	m_DetectorManagerKey {this, "DetectorManagerKey", "MuonDetectorManager","Key of input MuonDetectorManager condition data"}
SG::ReadHandleKey< Trk::SegmentCollection >	m_segm_type {this,"Eff_segm_type","TrackMuonSegments","muon segments"}
SG::ReadHandleKey< Muon::MMPrepDataContainer >	m_MMContainerKey {this,"MMPrepDataContainerName","MM_Measurements"}
SG::ReadHandleKey< xAOD::MuonContainer >	m_muonKey {this,"MuonKey","Muons","muons"}
SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_meTrkKey {this, "METrkContainer", "ExtrapolatedMuonTrackParticles"}
SG::ReadHandleKey< xAOD::NSWMMTPRDOContainer >	m_mmtpRdoKey {this, "NSW_MMTPDataKey", "NSW_MMTrigProcessor_RDO"}
Gaudi::Property< bool >	m_doMMESD {this,"DoMMESD",true}
Gaudi::Property< bool >	m_do_mm_overview {this,"do_mm_overview",true}
Gaudi::Property< bool >	m_do_stereoCorrection {this,"do_stereoCorrection",false}
Gaudi::Property< float >	m_cut_pt {this,"cut_pt",15000}
Gaudi::Property< bool >	m_doDetailedHists {this,"doDetailedHists",true}
std::string	m_name
std::unordered_map< std::string, size_t >	m_toolLookupMap
const ToolHandle< GenericMonitoringTool >	m_dummy
Gaudi::Property< bool >	m_enforceExpressTriggers
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 50 of file MMRawDataMonAlg.h.

Member Typedef Documentation

MonVarVec_t

typedef std::vector<std::reference_wrapper<Monitored::IMonitoredVariable> > AthMonitorAlgorithm::MonVarVec_t

privateinherited

Definition at line 370 of file AthMonitorAlgorithm.h.

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

DataType_t

enum class AthMonitorAlgorithm::DataType_t

stronginherited

Specifies what type of input data is being monitored.

An enumeration of the different types of data the monitoring application may be running over. This can be used to select which histograms to produce, e.g. to prevent the production of colliding-beam histograms when running on cosmic-ray data. Strings of the same names may be given as jobOptions.

Enumerator
userDefined
monteCarlo
collisions
cosmics
heavyIonCollisions

Definition at line 194 of file AthMonitorAlgorithm.h.

                          {
        userDefined = 0, 
        monteCarlo, 
        collisions, 
        cosmics, 
        heavyIonCollisions, 
    };

Environment_t

enum class AthMonitorAlgorithm::Environment_t

stronginherited

Specifies the processing environment.

The running environment may be used to select which histograms are produced, and where they are located in the output. For example, the output paths of the histograms are different for the "user", "online" and the various offline flags. Strings of the same names may be given as jobOptions.

Enumerator
user
online
tier0
tier0Raw
tier0ESD
AOD
altprod

Definition at line 175 of file AthMonitorAlgorithm.h.

                             {
        user = 0, 
        online, 
        tier0, 
        tier0Raw, 
        tier0ESD, 
        AOD, 
        altprod, 
    };

Constructor & Destructor Documentation

MMRawDataMonAlg()

MMRawDataMonAlg::MMRawDataMonAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 98 of file MMRawDataMonAlg.cxx.

                                                                                  : 
    AthMonitorAlgorithm(name,pSvcLocator)
{ }

~MMRawDataMonAlg()

virtual MMRawDataMonAlg::~MMRawDataMonAlg ( )

virtualdefault

Member Function Documentation

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

clusterFromSegments()

void MMRawDataMonAlg::clusterFromSegments ( const Trk::SegmentCollection * segms, const int lb ) const

private

Definition at line 845 of file MMRawDataMonAlg.cxx.

{
    MMOverviewHistogramStruct overviewPlots;
    MMByPhiStruct occupancyPlots[16][2];
    MMSummaryHistogramStruct summaryPlots[2][16][2][2][4];
    int nseg=0;
 
    for (Trk::SegmentCollection::const_iterator s = segms->begin(); s != segms->end(); ++s) {
        const Muon::MuonSegment* segment = dynamic_cast<const Muon::MuonSegment*>(*s);
        if (segment == nullptr) {
            ATH_MSG_DEBUG("no pointer to segment!!!");
            break;
        }
        bool isMM=false;
        for(unsigned int irot=0;irot<segment->numberOfContainedROTs();irot++){
            const Trk::RIO_OnTrack* rot = segment->rioOnTrack(irot);
            if(!rot) continue;
            Identifier rot_id = rot->identify();
            if(!m_idHelperSvc->isMM(rot_id)) continue;
        isMM=true;
            const Muon::MMClusterOnTrack* cluster = dynamic_cast<const Muon::MMClusterOnTrack*>(rot);
            if(!cluster) continue;
 
            std::string stName = m_idHelperSvc->mmIdHelper().stationNameString(m_idHelperSvc->mmIdHelper().stationName(rot_id));
        int stEta          = m_idHelperSvc->mmIdHelper().stationEta(rot_id);
        int stPhi          = m_idHelperSvc->mmIdHelper().stationPhi(rot_id);
        int multi          = m_idHelperSvc->mmIdHelper().multilayer(rot_id);
        int gap            = m_idHelperSvc->mmIdHelper().gasGap(rot_id);
        int ch             = m_idHelperSvc->mmIdHelper().channel(rot_id);
 
            // MMS and MML phi sectors                                                                                                                                                                                                  
            int sectorPhi = get_sectorPhi_from_stationPhi_stName(stPhi,stName);
            int PCB = get_PCB_from_channel(ch);
            int iside = (stEta > 0) ? 1 : 0;
 
            const Muon::MMPrepData* prd = cluster->prepRawData();
            const std::vector<Identifier>& stripIds = prd->rdoList();
            unsigned int csize = stripIds.size();
            const std::vector<uint16_t>& stripNumbers = prd->stripNumbers();
 
            auto& pcb_vects = summaryPlots[iside][sectorPhi-1][std::abs(stEta)-1][multi-1][gap-1];
            pcb_vects.cl_size.push_back(csize);
        pcb_vects.pcb.push_back(PCB);
            std::vector<short int> s_times = prd->stripTimes();
            float c_time = 0;
            for(unsigned int sIdx=0; sIdx<csize; ++sIdx) {
          pcb_vects.strp_times.push_back(s_times.at(sIdx));
          pcb_vects.pcb_strip.push_back( get_PCB_from_channel(stripNumbers[sIdx]));
          c_time += s_times.at(sIdx);
        }
            c_time /= s_times.size();
            pcb_vects.cl_times.push_back(c_time);
 
            float charge = prd->charge()*conversion_charge;
        pcb_vects.charge.push_back(charge);
 
            auto& vects = overviewPlots;
 
            auto& thisSect = occupancyPlots[sectorPhi-1][iside];
 
        const int gap_offset=4;
        int gas_gap8 = (multi==1) ?  gap :  gap + gap_offset;
        int FEB = get_FEB_from_channel(ch, stEta);
 
        int bin=get_bin_for_feb_occ(gas_gap8,FEB);
        thisSect.sector_lb_onseg.push_back(bin);
 
        if(stEta<0) {
          vects.sector_CSide_onseg.push_back(bin);
          vects.stationPhi_CSide_onseg.push_back(sectorPhi);
        } else {
          vects.sector_ASide_onseg.push_back(bin);
          vects.stationPhi_ASide_onseg.push_back(sectorPhi);
        }
 
        } // loop on ROT container                                                                                                                                                                                                           
    if (isMM==true) ++nseg;
 
    } // loop on segment collection                                                                                                                                                                  
    auto nsegs = Monitored::Scalar<int>("nseg",nseg);
    fill("mmMonitor", nsegs);
 
    auto& vects = overviewPlots;
    auto stationPhi_CSide_onseg = Monitored::Collection("stationPhi_CSide_onseg",vects.stationPhi_CSide_onseg);
    auto stationPhi_ASide_onseg = Monitored::Collection("stationPhi_ASide_onseg",vects.stationPhi_ASide_onseg);
    auto sector_ASide_onseg = Monitored::Collection("sector_ASide_onseg",vects.sector_ASide_onseg);
    auto sector_CSide_onseg = Monitored::Collection("sector_CSide_onseg",vects.sector_CSide_onseg);
 
    auto lb_onseg = Monitored::Scalar<int>("lb_onseg", lb);
 
    fill("mmMonitor", stationPhi_CSide_onseg, stationPhi_ASide_onseg, sector_CSide_onseg, sector_ASide_onseg, lb_onseg);
 
    for(int iside = 0; iside < 2; ++iside) {
        std::string MM_sideGroup = "MM_sideGroup" + MM_Side[iside];
        for(int statPhi=0; statPhi<16; ++statPhi) {
            auto& occ_lb = occupancyPlots[statPhi][iside];
            auto sector_lb_onseg = Monitored::Collection("sector_lb_"+MM_Side[iside]+"_phi"+std::to_string(statPhi+1)+"_onseg",occ_lb.sector_lb_onseg);
            fill(MM_sideGroup, lb_onseg, sector_lb_onseg);
 
            for(int statEta = 0; statEta < 2; ++statEta) {
                for(int multiplet = 0; multiplet < 2; ++multiplet) {
                    for(int gas_gap = 0; gas_gap < 4; ++gas_gap) {
                        auto& pcb_vects = summaryPlots[iside][statPhi][statEta][multiplet][gas_gap];
 
                        if(pcb_vects.pcb.empty()) continue;
                    if(m_doDetailedHists){
                      auto pcb_mon = Monitored::Collection("pcb_mon_onseg_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), pcb_vects.pcb);
                      auto pcb_strip_mon = Monitored::Collection("pcb_strip_mon_onseg_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), pcb_vects.pcb_strip);
                      auto strip_times = Monitored::Collection("strp_time_onseg_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), pcb_vects.strp_times);
                      auto cluster_time = Monitored::Collection("cluster_time_onseg_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), pcb_vects.cl_times);
                      auto clus_size = Monitored::Collection("cluster_size_onseg_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), pcb_vects.cl_size);
                      auto charge_perPCB = Monitored::Collection("charge_perPCB_onseg_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), pcb_vects.charge);
                      
                      fill(MM_sideGroup, clus_size, strip_times, charge_perPCB, cluster_time, pcb_mon, pcb_strip_mon);
                    }
                    auto clus_size_all = Monitored::Collection("cluster_size_onseg", pcb_vects.cl_size); 
                    auto charge_all = Monitored::Collection("charge_onseg", pcb_vects.charge);
                    auto strip_times_all = Monitored::Collection("strp_time_onseg", pcb_vects.strp_times);
                    fill("mmMonitor", clus_size_all, charge_all, strip_times_all);
                    }
                }
            }
        }
    }
}

clusterFromTrack()

void MMRawDataMonAlg::clusterFromTrack ( const xAOD::TrackParticleContainer * muonContainer, const int lb ) const

private

Definition at line 497 of file MMRawDataMonAlg.cxx.

{
    MMSummaryHistogramStruct summaryPlots[2][2][4]; // side, multilayer, gas gap
    MMSummaryHistogramStruct summaryPlots_full[2][16][2][2][4]; // side, phi, eta, multilayer, gas gap
    MMSummaryHistogramStruct sumPlots[2][16][2][2][4]; // side, phi, eta, multilayer, gas gap
    MMOverviewHistogramStruct overviewPlots;
    MMByPhiStruct occupancyPlots[16][2]; // sector, side
    int ntrk=0;
    for(const xAOD::TrackParticle* meTP : *muonContainer) {
 
        if(!meTP) continue;
        auto eta_trk = Monitored::Scalar<float>("eta_trk", meTP->eta());
        auto phi_trk = Monitored::Scalar<float>("phi_trk", meTP->phi());
        auto pt_trk = Monitored::Scalar<float>("pt_trk", meTP->pt()/1000.);
 
        //retrieve the original track
        const Trk::Track* meTrack = meTP->track();
 
        if(!meTrack) continue;
        
        // get the vector of measurements on track
        const DataVector<const Trk::MeasurementBase>* meas = meTrack->measurementsOnTrack();
        bool isMM=false;
        for(const Trk::MeasurementBase* it : *meas) {
            const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(it);
            if(!rot) continue;
            Identifier rot_id = rot->identify();
            if(!m_idHelperSvc->isMM(rot_id)) continue;
            isMM=true;
            const Muon::MMClusterOnTrack* cluster = dynamic_cast<const Muon::MMClusterOnTrack*>(rot);
            if(!cluster) continue;
 
            std::string stName = m_idHelperSvc->mmIdHelper().stationNameString(m_idHelperSvc->mmIdHelper().stationName(rot_id));
            int stEta          = m_idHelperSvc->mmIdHelper().stationEta(rot_id);
            int stPhi          = m_idHelperSvc->mmIdHelper().stationPhi(rot_id);
            int multi          = m_idHelperSvc->mmIdHelper().multilayer(rot_id);
            int gap            = m_idHelperSvc->mmIdHelper().gasGap(rot_id);
            int ch             = m_idHelperSvc->mmIdHelper().channel(rot_id);
            
            // MMS and MML phi sectors
            //              int phisec = (stNumber%2==0) ? 1 : 0;
            int sectorPhi = get_sectorPhi_from_stationPhi_stName(stPhi,stName); // 1->16
            int PCB = get_PCB_from_channel(ch);
            int iside = (stEta > 0) ? 1 : 0;
            auto& vects = overviewPlots;
            auto& thisSect = occupancyPlots[sectorPhi-1][iside];
            
 
            const Muon::MMPrepData* prd = cluster->prepRawData();
            const std::vector<Identifier>& stripIds = prd->rdoList();
            unsigned int csize = stripIds.size();
            const std::vector<uint16_t>& stripNumbers = prd->stripNumbers();
            float charge = prd->charge()*conversion_charge;
            std::vector<short int> s_times = prd->stripTimes();
 
            vects.charge_all.push_back(charge);
            
            float c_time = 0;
            for(unsigned int sIdx=0; sIdx<stripIds.size(); ++sIdx){
              vects.strp_times.push_back(s_times.at(sIdx));
              c_time += s_times.at(sIdx);
            }
            c_time /= s_times.size();
            vects.cl_times.push_back(c_time);
 
            if(m_doDetailedHists){
              auto& vect = sumPlots[iside][sectorPhi-1][std::abs(stEta)-1][multi-1][gap-1];
              vect.cl_size.push_back(csize);
              vect.pcb.push_back(PCB);
              for(unsigned int sIdx=0; sIdx<stripIds.size(); ++sIdx)
                { 
                  vect.strip_number.push_back(stripNumbers[sIdx]);
                  vect.strp_times.push_back(s_times.at(sIdx));
                  vect.pcb_strip.push_back(get_PCB_from_channel(stripNumbers[sIdx]));
                }
              vect.cl_times.push_back(c_time);
              vect.charge.push_back(charge);
            }
 
 
            const int gap_offset=4;
            int gas_gap8 = (multi==1) ?  gap :  gap + gap_offset; 
 
            int FEB = get_FEB_from_channel(ch, stEta);
            int bin = get_bin_for_feb_occ(gas_gap8,FEB);
            thisSect.sector_lb_ontrack.push_back(bin);
            // Occupancy plots with FEB granularity further divided for each eta sector: -2, -1, 1, 2
            // Filling Vectors for stationEta=-1 - cluster on track
            if(stEta<0) {
              vects.sector_CSide_ontrack.push_back(bin);
              vects.stationPhi_CSide_ontrack.push_back(sectorPhi);
            } else {
              vects.sector_ASide_ontrack.push_back(bin);
              vects.stationPhi_ASide_ontrack.push_back(sectorPhi);
            }
 
            float x = cluster->localParameters()[Trk::loc1];
            for(const Trk::TrackStateOnSurface* trkState : *meTrack->trackStateOnSurfaces()) {
                    
                if(!(trkState)) continue;
                if (!trkState->type(Trk::TrackStateOnSurface::Measurement)) continue;
 
                Identifier surfaceId = (trkState)->surface().associatedDetectorElementIdentifier();
                if(!m_idHelperSvc->isMM(surfaceId)) continue;
            
                int trk_stEta = m_idHelperSvc->mmIdHelper().stationEta(surfaceId);
                int trk_stPhi = m_idHelperSvc->mmIdHelper().stationPhi(surfaceId);
                int trk_multi = m_idHelperSvc->mmIdHelper().multilayer(surfaceId);
                int trk_gap   = m_idHelperSvc->mmIdHelper().gasGap(surfaceId);
 
                if( (trk_stPhi == stPhi) && (trk_stEta == stEta) && (trk_multi == multi) && (trk_gap == gap)) {
                  double x_trk = trkState->trackParameters()->parameters()[Trk::loc1];
                  int sectorPhi = get_sectorPhi_from_stationPhi_stName(trk_stPhi,stName); // 1->16
                  int side  = (stEta > 0) ? 1 : 0;
                  float res_stereo = (x - x_trk);
                    if(m_do_stereoCorrection) {
                        float stereo_angle = ((multi == 1 && gap < 3) || (multi == 2 && gap > 2)) ? 0 : 0.02618;
                        double y_trk = trkState->trackParameters()->parameters()[Trk::locY];
                        float stereo_correction = ( (multi == 1 && gap < 3) || (multi == 2 && gap > 2) ) ? 0 : ( ((multi == 1 && gap == 3) || (multi == 2 && gap ==1 )) ? (-std::sin(stereo_angle)*y_trk) : std::sin(stereo_angle)*y_trk );
                        res_stereo = (x - x_trk)*std::cos(stereo_angle) - stereo_correction;
                    }
                    auto residual_mon = Monitored::Scalar<float>("residual", res_stereo);
                    auto stPhi_mon = Monitored::Scalar<float>("stPhi_mon",sectorPhi);
 
                    fill("mmMonitor", residual_mon, eta_trk, phi_trk, stPhi_mon);
                    int abs_stEta = get_sectorEta_from_stationEta(stEta); // 0 or 1
 
                    if(m_doDetailedHists){
                      auto& vectors = summaryPlots_full[side][sectorPhi-1][abs_stEta][multi-1][gap-1];
                      vectors.residuals.push_back(res_stereo);
                    }
                }
            }//TrackStates
 
        } // loop on meas
        if(isMM) {
          ++ntrk;
          fill("mmMonitor", pt_trk);
        }
 
        if(m_doDetailedHists){
          for(int iside = 0; iside < 2; ++iside) {
            std::string MM_sideGroup = "MM_sideGroup" + MM_Side[iside];
            for(int statPhi = 0; statPhi < 16; ++statPhi) {
              //                for(int statEta = 0; statEta < 2; ++statEta) {
              for(int multiplet = 0; multiplet < 2; ++multiplet) {
            for(int gas_gap = 0; gas_gap < 4; ++gas_gap) {
              auto layer=gas_gap+multiplet*4;
              MMSummaryHistogramStruct vects;
              for(int statEta = 0; statEta < 2; ++statEta) {
                vects = summaryPlots_full[iside][statPhi][statEta][multiplet][gas_gap];
                auto residuals_gap = Monitored::Collection("residuals_"+MM_Side[iside]+"_phi"+std::to_string(statPhi+1)+"_stationEta"+EtaSector[statEta]+"_multiplet"+std::to_string(multiplet+1)+"_gas_gap"+std::to_string(gas_gap+1),vects.residuals);
                auto residuals_layer = Monitored::Collection("residuals_"+MM_Side[iside]+"_phi"+std::to_string(statPhi+1)+"_layer"+std::to_string(layer+1),vects.residuals);
                
                fill(MM_sideGroup, residuals_gap,residuals_layer);
              }
            }
              }
            }
          }
        }
 
        for(const Trk::TrackStateOnSurface* trkState : *meTrack->trackStateOnSurfaces()) {
            if(!(trkState)) continue;
            if (!trkState->type(Trk::TrackStateOnSurface::Measurement)) continue;
            Identifier surfaceId = (trkState)->surface().associatedDetectorElementIdentifier();
            if(!m_idHelperSvc->isMM(surfaceId)) continue;
 
            const Trk::MeasurementBase* meas = trkState->measurementOnTrack() ;
            if(!meas) continue;
            
            const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(meas);
                        if(!rot) continue;
                        Identifier rot_id = rot->identify();
                        if(!m_idHelperSvc->isMM(rot_id)) continue;
            
            const Amg::Vector3D& pos = (trkState)->trackParameters()->position();
            int stEta = m_idHelperSvc->mmIdHelper().stationEta(surfaceId);
            int multi = m_idHelperSvc->mmIdHelper().multilayer(surfaceId);
            int gap   = m_idHelperSvc->mmIdHelper().gasGap(surfaceId);
 
            // CSide and ASide
            int iside = (stEta > 0) ? 1 : 0;
            auto& Vectors = summaryPlots[iside][multi-1][gap-1];
 
            // Filling x-y position vectors using the trackStateonSurface
            Vectors.x_ontrack.push_back(pos.x());
            Vectors.y_ontrack.push_back(pos.y());
        }
    } // loop on muonContainer
 
    auto ntrack = Monitored::Scalar<int>("ntrk",ntrk);
    fill("mmMonitor", ntrack);
 
    auto& vects = overviewPlots;
    auto stationPhi_CSide_ontrack = Monitored::Collection("stationPhi_CSide_ontrack",vects.stationPhi_CSide_ontrack);
    auto stationPhi_ASide_ontrack = Monitored::Collection("stationPhi_ASide_ontrack",vects.stationPhi_ASide_ontrack);
    auto sector_ASide_ontrack = Monitored::Collection("sector_ASide_ontrack",vects.sector_ASide_ontrack);
    auto sector_CSide_ontrack = Monitored::Collection("sector_CSide_ontrack",vects.sector_CSide_ontrack);
 
    auto lb_ontrack = Monitored::Scalar<int>("lb_ontrack", lb);
    auto csize = Monitored::Collection("nstrips_ontrack", vects.numberofstrips_percluster);
    auto charge = Monitored::Collection("charge_ontrack", vects.charge_all);
    auto stime = Monitored::Collection("strip_time_on_track", vects.strp_times);
    auto ctime = Monitored::Collection("cluster_time_on_track", vects.cl_times);
 
    fill("mmMonitor", csize, charge, stime, ctime, stationPhi_CSide_ontrack, stationPhi_ASide_ontrack, sector_CSide_ontrack,sector_ASide_ontrack, lb_ontrack);
 
    for(int iside = 0; iside < 2; ++iside) {
        std::string MM_sideGroup = "MM_sideGroup" + MM_Side[iside];
        for(int statPhi = 0; statPhi < 16; ++statPhi) {
            for(int statEta = 0; statEta < 2; ++statEta) {
                for(int multiplet = 0; multiplet < 2; ++multiplet) {
                    for(int gas_gap = 0; gas_gap < 4; ++gas_gap) {
                        auto& vects = sumPlots[iside][statPhi][statEta][multiplet][gas_gap];
                        if(m_doDetailedHists){
                          if(!vects.strip_number.empty())
                            {
                              auto clus_size = Monitored::Collection("cluster_size_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.cl_size);
                              auto strip_times = Monitored::Collection("strp_time_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.strp_times);
                              auto cluster_time = Monitored::Collection("cluster_time_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.cl_times);
                              auto charge_perPCB = Monitored::Collection("charge_perPCB_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.charge);
                              auto charge_perlayer = Monitored::Collection("charge_perlayer_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.charge);
                              auto clus_size_perlayer = Monitored::Collection("cluster_size_perlayer_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.cl_size);
                              auto pcb_mon = Monitored::Collection("pcb_mon_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.pcb);
                              auto pcb_strip_mon = Monitored::Collection("pcb_strip_mon_ontrack_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), vects.pcb_strip);
                              
                              fill(MM_sideGroup, clus_size, strip_times, cluster_time, charge_perPCB, pcb_mon, pcb_strip_mon,charge_perlayer,clus_size_perlayer);
                            }
                        }
 
                    }
                }
            }
            auto& occ_lb = occupancyPlots[statPhi][iside];
            auto sector_lb_ontrack = Monitored::Collection("sector_lb_"+MM_Side[iside]+"_phi"+std::to_string(statPhi+1)+"_ontrack",occ_lb.sector_lb_ontrack);
            fill(MM_sideGroup, lb_ontrack, sector_lb_ontrack);
        }
        for(int multiplet=0; multiplet<2; ++multiplet) {
            for(int gas_gap=0; gas_gap<4; ++gas_gap) {
                auto& Vectors = summaryPlots[iside][multiplet][gas_gap];
                auto x_ontrack = Monitored::Collection("x_"+MM_Side[iside]+"_multiplet"+std::to_string(multiplet+1)+"_gas_gap_"+std::to_string(gas_gap+1)+"_ontrack", Vectors.x_ontrack);
                auto y_ontrack = Monitored::Collection("y_"+MM_Side[iside]+"_multiplet"+std::to_string(multiplet+1)+"_gas_gap_"+std::to_string(gas_gap+1)+"_ontrack", Vectors.y_ontrack);
                fill(MM_sideGroup, x_ontrack, y_ontrack);
            }
        }
    }
 
}

dataType()

DataType_t AthMonitorAlgorithm::dataType ( ) const

inlineinherited

Accessor functions for the data type.

Returns

the current value of the class's DataType_t instance.

Definition at line 224 of file AthMonitorAlgorithm.h.

{ return m_dataType; }

dataTypeStringToEnum()

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::dataTypeStringToEnum ( const std::string & str ) const

inherited

Convert the data type string from the python configuration to an enum object.

Returns

a value in the DataType_t enumeration which matches the input string.

Definition at line 144 of file AthMonitorAlgorithm.cxx.

                                                                                                    {
    // convert the string to all lowercase
    std::string lowerCaseStr = str;
    std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
 
    // check if it matches one of the enum choices
    if( lowerCaseStr == "userdefined" ) {
        return DataType_t::userDefined;
    } else if( lowerCaseStr == "montecarlo" ) {
        return DataType_t::monteCarlo;
    } else if( lowerCaseStr == "collisions" ) {
        return DataType_t::collisions;
    } else if( lowerCaseStr == "cosmics" ) {
        return DataType_t::cosmics;
    } else if( lowerCaseStr == "heavyioncollisions" ) {
        return DataType_t::heavyIonCollisions;
    } else { // otherwise, warn the user and return "userDefined"
        ATH_MSG_WARNING("AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
            <<str<<", returning userDefined.");
        return DataType_t::userDefined;
    }
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

environment()

Environment_t AthMonitorAlgorithm::environment ( ) const

inlineinherited

Accessor functions for the environment.

Returns

the current value of the class's Environment_t instance.

Definition at line 208 of file AthMonitorAlgorithm.h.

{ return m_environment; }

envStringToEnum()

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::envStringToEnum ( const std::string & str ) const

inherited

Convert the environment string from the python configuration to an enum object.

Returns

a value in the Environment_t enumeration which matches the input string.

Definition at line 116 of file AthMonitorAlgorithm.cxx.

                                                                                                  {
    // convert the string to all lowercase
    std::string lowerCaseStr = str;
    std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
 
    // check if it matches one of the enum choices
    if( lowerCaseStr == "user" ) {
        return Environment_t::user;
    } else if( lowerCaseStr == "online" ) {
        return Environment_t::online;
    } else if( lowerCaseStr == "tier0" ) {
        return Environment_t::tier0;
    } else if( lowerCaseStr == "tier0raw" ) {
        return Environment_t::tier0Raw;
    } else if( lowerCaseStr == "tier0esd" ) {
        return Environment_t::tier0ESD;
    } else if( lowerCaseStr == "aod" ) {
        return Environment_t::AOD;
    } else if( lowerCaseStr == "altprod" ) {
        return Environment_t::altprod;
    } else { // otherwise, warn the user and return "user"
        ATH_MSG_WARNING("AthMonitorAlgorithm::envStringToEnum(): Unknown environment "
            <<str<<", returning user.");
        return Environment_t::user;
    }
}

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode AthMonitorAlgorithm::execute ( const EventContext & ctx ) const

overridevirtualinherited

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters

ctx	event context for reentrant Athena call

Returns

StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

                                                                       {
 
    // Checks that all of the  DQ filters are passed. If any one of the filters
    // fails, return SUCCESS code and do not fill the histograms with the event.
    for ( const auto& filterItr : m_DQFilterTools ) {
        if (!filterItr->accept()) {
            ATH_MSG_DEBUG("Event rejected due to filter tool.");
            return StatusCode::SUCCESS;
        }
    }
 
    // Trigger: If there is a decision tool and the chains fail, skip the event.
    if ( !m_trigDecTool.empty() && !trigChainsArePassed(m_vTrigChainNames) ) {
        ATH_MSG_DEBUG("Event rejected due to trigger filter.");
        return StatusCode::SUCCESS;
    }
 
    ATH_MSG_DEBUG("Event accepted!");
    return fillHistograms(ctx);
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

fillHistograms()

StatusCode MMRawDataMonAlg::fillHistograms ( const EventContext & ctx ) const

overridevirtual

adds event to the monitoring histograms

User will overwrite this function. Histogram booking is no longer done in C++. This function is called in execute once the filters are all passed.

Parameters

ctx	forwarded from execute

Returns

StatusCode

Implements AthMonitorAlgorithm.

Definition at line 129 of file MMRawDataMonAlg.cxx.

{
    int lumiblock = -1;
    lumiblock = GetEventInfo(ctx)->lumiBlock();
    ATH_MSG_DEBUG("MMRawDataMonAlg::MM RawData Monitoring Histograms being filled" );
 
    SG::ReadHandle<Muon::MMPrepDataContainer> mm_container(m_MMContainerKey,ctx);
    ATH_MSG_DEBUG("****** mmContainer->size() : " << mm_container->size());
 
    if(m_doMMESD) {
        MMOverviewHistogramStruct overviewPlots;
        MMSummaryHistogramStruct summaryPlots[2][16][2][2][4];
        MMByPhiStruct occupancyPlots[16][2];
 
        //loop in MMPrepDataContainer
        for(const Muon::MMPrepDataCollection* coll : *mm_container) {
          for(const Muon::MMPrepData* prd : *coll) {
            ATH_CHECK(fillMMOverviewVects(prd, overviewPlots, occupancyPlots));
                ATH_CHECK(fillMMSummaryVects(prd, summaryPlots));
                ATH_CHECK(fillMMHistograms(prd));
          }
        }
        
        if(m_do_mm_overview) fillMMOverviewHistograms(overviewPlots, occupancyPlots, lumiblock);
        
        ATH_CHECK(fillMMSummaryHistograms(summaryPlots));
        SG::ReadHandle<xAOD::TrackParticleContainer> meTPContainer{m_meTrkKey,ctx};
        if (!meTPContainer.isValid()) {
          ATH_MSG_FATAL("Nope. Could not retrieve "<<m_meTrkKey.fullKey());
          return StatusCode::FAILURE;
        }
        clusterFromTrack(meTPContainer.cptr(),lumiblock);
        MMEfficiency(meTPContainer.cptr());
 
 
        
        //trigger
        SG::ReadHandle<xAOD::NSWMMTPRDOContainer> rdos =  SG::ReadHandle<xAOD::NSWMMTPRDOContainer>{m_mmtpRdoKey, ctx};
        if (not rdos.isValid()) {
          ATH_MSG_INFO("NSW MMTP failed. Skipping");
          // return StatusCode::SUCCESS;
        } else fillMMTrigger(rdos.cptr(),lumiblock);
        
        SG::ReadHandle<Trk::SegmentCollection> segms(m_segm_type, ctx);
 
        if (!segms.isValid()) {
          ATH_MSG_INFO("evtStore() does not contain MM segms Collection with name " << m_segm_type);
          //          return StatusCode::FAILURE;
        }else
          clusterFromSegments(segms.cptr(),lumiblock);
    }
 
 
    return StatusCode::SUCCESS;
}

fillMMHistograms()

StatusCode MMRawDataMonAlg::fillMMHistograms ( const Muon::MMPrepData * ) const

private

Definition at line 393 of file MMRawDataMonAlg.cxx.

                                                                         {
  return StatusCode::SUCCESS;
}

fillMMOverviewHistograms()

void MMRawDataMonAlg::fillMMOverviewHistograms ( const MMOverviewHistogramStruct & vects, MMByPhiStruct(&) occupancyPlots[16][2], const int lb ) const

private

Definition at line 259 of file MMRawDataMonAlg.cxx.

{
    auto charge_all = Monitored::Collection("charge_all", vects.charge_all);
    auto numberofstrips_percluster = Monitored::Collection("numberofstrips_percluster", vects.numberofstrips_percluster);
    fill("mmMonitor", charge_all, numberofstrips_percluster);
 
    auto strip_times = Monitored::Collection("strip_times", vects.strp_times);
    auto cluster_times = Monitored::Collection("cluster_times", vects.cl_times);
    auto strip_number = Monitored::Collection("strip_number", vects.strip_number);
    auto statEta_strip = Monitored::Collection("statEta_strip", vects.statEta_strip);
    fill("mmMonitor", strip_times, cluster_times, strip_number, statEta_strip);
 
    auto x_mon = Monitored::Collection("x_mon", vects.x_mon);
    auto y_mon = Monitored::Collection("y_mon", vects.y_mon);
    auto z_mon = Monitored::Collection("z_mon", vects.z_mon);
    auto R_mon = Monitored::Collection("R_mon", vects.R_mon);
    fill("mmMonitor", x_mon, y_mon, z_mon, R_mon);
 
    auto lb_mon = Monitored::Scalar<int>("lb_mon", lb);
 
    for(int statPhi=0; statPhi<16; ++statPhi) {
        for(int iside=0; iside<2; ++iside) {
            auto& occ_lb = occupancyPlots[statPhi][iside];
            auto sector_lb = Monitored::Collection("sector_lb_"+MM_Side[iside]+"_phi"+std::to_string(statPhi+1),occ_lb.sector_lb);
            std::string MM_sideGroup = "MM_sideGroup" + MM_Side[iside];
            fill(MM_sideGroup, lb_mon, sector_lb);
        }
    }
    auto sector_CSide      = Monitored::Collection("sector_CSide",vects.sector_CSide);
    auto sector_ASide      = Monitored::Collection("sector_ASide",vects.sector_ASide);
    auto stationPhi_CSide = Monitored::Collection("stationPhi_CSide",vects.stationPhi_CSide);
    auto stationPhi_ASide = Monitored::Collection("stationPhi_ASide",vects.stationPhi_ASide);
 
    fill("mmMonitor", sector_CSide, sector_ASide, stationPhi_CSide, stationPhi_ASide );
}

fillMMOverviewVects()

StatusCode MMRawDataMonAlg::fillMMOverviewVects ( const Muon::MMPrepData * prd, MMOverviewHistogramStruct & vects, MMByPhiStruct(&) occupancyPlots[16][2] ) const

private

Definition at line 185 of file MMRawDataMonAlg.cxx.

{
    Identifier Id = prd->identify();
    const std::vector<Identifier>& stripIds = prd->rdoList();
    unsigned int nStrips = stripIds.size(); // number of strips in this cluster (cluster size)
    const std::vector<uint16_t>& stripNumbers = prd->stripNumbers();
 
    std::string stName = m_idHelperSvc->mmIdHelper().stationNameString(m_idHelperSvc->mmIdHelper().stationName(Id));
    int gas_gap        = m_idHelperSvc->mmIdHelper().gasGap(Id);
    int stationEta       = m_idHelperSvc->mmIdHelper().stationEta(Id);
    int stationPhi     = m_idHelperSvc->mmIdHelper().stationPhi(Id);
    int multiplet      = m_idHelperSvc->mmIdHelper().multilayer(Id);
    int channel        = m_idHelperSvc->mmIdHelper().channel(Id);
 
    // Returns the charge (number of electrons) converted in fC
    float charge = prd->charge()*conversion_charge;
    // Returns the times of each strip (in ns)
    std::vector<short int> strip_times = prd->stripTimes();
 
    Amg::Vector3D pos = prd->globalPosition();
    float R = std::hypot(pos.x(),pos.y());
 
    // MM gaps are back to back, so the direction of the drift (time) is different for the even and odd gaps -> flip for the even gaps
 
    vects.charge_all.push_back(charge);
    vects.numberofstrips_percluster.push_back(nStrips);
    vects.x_mon.push_back(pos.x());
    vects.y_mon.push_back(pos.y());
    vects.z_mon.push_back(pos.z());
    vects.R_mon.push_back(R);
 
    // 16 phi sectors, 8 stationPhi times 2 stName, MMS and MML
    int sectorPhi = get_sectorPhi_from_stationPhi_stName(stationPhi,stName);
 
    // Occupancy plots with PCB granularity further divided for each eta sector: -2, -1, 1, 2
    // CSide and ASide
    int iside = (stationEta>0) ? 1 : 0;
 
    auto& thisSect = occupancyPlots[sectorPhi-1][iside];
    const int gap_offset=4;
    int gas_gap8 = (multiplet==1) ?  gas_gap :  gas_gap + gap_offset; 
    int FEB = get_FEB_from_channel(channel, stationEta);
    
    int bin = get_bin_for_feb_occ(gas_gap8,FEB);
 
    thisSect.sector_lb.push_back(bin);
 
    if(stationEta<0) {
      vects.sector_CSide.push_back(bin);
      vects.stationPhi_CSide.push_back(sectorPhi);
    } else {
      vects.sector_ASide.push_back(bin);
      vects.stationPhi_ASide.push_back(sectorPhi);
    }
 
    // loop on each strip
    int sIdx = 0; // index-counter for the vector of Id's
    float cluster_time = 0;
    for(const Identifier& id: stripIds) {
        
        std::string stName_strip = m_idHelperSvc->mmIdHelper().stationNameString(m_idHelperSvc->mmIdHelper().stationName(id));
        int stationEta_strip     = m_idHelperSvc->mmIdHelper().stationEta(id);
        vects.statEta_strip.push_back(stationEta_strip);
        vects.strip_number.push_back(stripNumbers[sIdx]);
        vects.strp_times.push_back(strip_times.at(sIdx));
        cluster_time += strip_times.at(sIdx);
        ++sIdx;
    }
    cluster_time /= strip_times.size();
    vects.cl_times.push_back(cluster_time);
 
    return StatusCode::SUCCESS;
}

fillMMSummaryHistograms()

StatusCode MMRawDataMonAlg::fillMMSummaryHistograms ( const MMSummaryHistogramStruct(&) vects[2][16][2][2][4] ) const

private

Definition at line 351 of file MMRawDataMonAlg.cxx.

                                                                                                                  {
 
    for(int iside=0; iside<2; ++iside) {
        std::string MM_sideGroup = "MM_sideGroup" + MM_Side[iside];
 
 
 
        for(int statPhi=0; statPhi<16; ++statPhi) {
          for(int multiplet=0; multiplet<2; ++multiplet) {
 
 
            for(int statEta=0; statEta<2; ++statEta) {
 
              for(int gas_gap=0; gas_gap<4; ++gas_gap) {
                        auto& Vectors = vects[iside][statPhi][statEta][multiplet][gas_gap];
                        auto sector_strip = Monitored::Collection("sector_strip_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1), Vectors.sector_strip);
                        auto strip_number = Monitored::Collection("strip_number_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1), Vectors.strip_number);
                        if(m_doDetailedHists){
                          if(!Vectors.strip_number.empty())
                            {
                              auto cluster_size = Monitored::Collection("cluster_size_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.cl_size);
                              auto strip_times = Monitored::Collection("strp_time_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.strp_times);
                              auto cluster_time = Monitored::Collection("cluster_time_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.cl_times);
                              auto charge_perPCB = Monitored::Collection("charge_perPCB_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.charge);
                              auto charge_perlayer = Monitored::Collection("charge_perlayer_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.charge);
                              auto cluster_size_perlayer = Monitored::Collection("cluster_size_perlayer_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.cl_size);
                              auto pcb_mon = Monitored::Collection("pcb_mon_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.pcb);
                              auto pcb_strip_mon = Monitored::Collection("pcb_strip_mon_" + MM_Side[iside] + "_phi" + std::to_string(statPhi+1) + "_eta" + std::to_string(statEta+1) + "_ml" + std::to_string(multiplet+1) + "_gap" + std::to_string(gas_gap+1), Vectors.pcb_strip);
                              fill(MM_sideGroup, cluster_size, strip_times, cluster_time, charge_perPCB, pcb_mon, pcb_strip_mon, charge_perlayer, cluster_size_perlayer);
                            }
                        }
                    
                        fill(MM_sideGroup, strip_number, sector_strip);
                    }
                }
            }
        }
    }
 
    return StatusCode::SUCCESS;
}

fillMMSummaryVects()

StatusCode MMRawDataMonAlg::fillMMSummaryVects ( const Muon::MMPrepData * prd, MMSummaryHistogramStruct(&) vects[2][16][2][2][4] ) const

private

Definition at line 295 of file MMRawDataMonAlg.cxx.

{
  Identifier Id = prd->identify();
  const std::vector<Identifier>& stripIds = prd->rdoList();
  
  std::string stName      = m_idHelperSvc->mmIdHelper().stationNameString(m_idHelperSvc->mmIdHelper().stationName(Id));
  int thisStationEta      = m_idHelperSvc->mmIdHelper().stationEta(Id);
  int thisStationPhi       = m_idHelperSvc->mmIdHelper().stationPhi(Id);
  int thisMultiplet        = m_idHelperSvc->mmIdHelper().multilayer(Id);
  int thisGasgap          = m_idHelperSvc->mmIdHelper().gasGap(Id);
  int ch             = m_idHelperSvc->mmIdHelper().channel(Id);
  float thisCharge=prd->charge()*conversion_charge;
  std::vector<short int> strip_times = prd->stripTimes();
 
    
    int phi=get_sectorPhi_from_stationPhi_stName(thisStationPhi ,stName);
 
    // CSide and ASide
    int iside = (thisStationEta>0) ? 1 : 0;
 
    // 2 eta sectors depending on Eta=+-1 (0) and +-2 (1)
    int sectorEta=get_sectorEta_from_stationEta(thisStationEta);
    unsigned int csize = stripIds.size();
    int PCB = get_PCB_from_channel(ch);
    auto& Vectors = vects[iside][phi-1][sectorEta][thisMultiplet-1][thisGasgap-1];    
 
    // loop on strips
    int sIdx = 0;
    const std::vector<uint16_t>& stripNumbers=prd->stripNumbers();
    float cluster_time = 0;
    for ( const Identifier& id : stripIds) {
      
      int stationEta       = m_idHelperSvc->mmIdHelper().stationEta(id);
      int gas_gap          = m_idHelperSvc->mmIdHelper().gasGap(Id);
      int multiplet        = m_idHelperSvc->mmIdHelper().multilayer(Id);
      // Filling Vectors for both sides, considering each strip
      if(m_doDetailedHists){
    Vectors.strp_times.push_back(strip_times.at(sIdx));
    cluster_time += strip_times.at(sIdx);
      }
      Vectors.strip_number.push_back(stripNumbers[sIdx]);
      Vectors.pcb_strip.push_back( get_PCB_from_channel(stripNumbers[sIdx]));
      ++sIdx;
      if(iside==1)    Vectors.sector_strip.push_back(get_bin_for_occ_ASide_hist(stationEta,multiplet,gas_gap));
      if(iside==0)    Vectors.sector_strip.push_back(get_bin_for_occ_CSide_hist(stationEta,multiplet,gas_gap));
    }
    if(m_doDetailedHists){
      Vectors.cl_size.push_back(csize);
      Vectors.pcb.push_back(PCB);
      cluster_time /= strip_times.size();
      Vectors.cl_times.push_back(cluster_time);
      Vectors.charge.push_back(thisCharge);
    }      
    return StatusCode::SUCCESS;
}

fillMMTrigger()

void MMRawDataMonAlg::fillMMTrigger ( const xAOD::NSWMMTPRDOContainer * mmtp, const int lb ) const

private

Definition at line 397 of file MMRawDataMonAlg.cxx.

                                                                                           {
 
    auto lb_tri=Monitored::Scalar<int>("lb_tri",lb);
 
    for (const auto* rdo : *mmtp) {
       auto sourceID = rdo->sourceID();
       auto moduleID = rdo->moduleID();
       int s_side = (sourceID >> 16) == 107 ? 1:-1;                  // 1 for A; -1 for C
       auto side = Monitored::Scalar<int>("tri_side", s_side);
       int iside= s_side>0 ? s_side : 0;                
       uint s_sector = (moduleID & 0xF) + 1;                  //0-15 --> 1-16
       int oct = (int)((s_sector-1)/2.);
       float sector_pos=(45/180.)*M_PI*oct; //large
       if(s_sector%2==0 ) sector_pos=(45*oct+22.5)*M_PI/180.; //small
       if(sector_pos>M_PI) sector_pos=sector_pos-2*M_PI;
 
       auto trig_sector = Monitored::Scalar<int>("trig_sector", s_sector*s_side);
 
       auto event_bcid=rdo->ROD_BCID();
 
       auto bcid=Monitored::Scalar<int>("bcid",event_bcid);
 
       std::vector<short unsigned int> art_bcids = rdo->art_BCID();
       std::vector<unsigned char> layers=rdo->art_layer();
       auto channels=rdo->art_channel();
 
       fill("mmTrigger", trig_sector, lb_tri);
 
       for (long unsigned int i=0; i< rdo->art_BCID().size(); i++ ){
     auto art_layer=static_cast<unsigned int>(layers[i]);
     auto art_channel = Monitored::Scalar<int>("art_channel", channels[i]);
     auto art_sector_layer = Monitored::Scalar<int>("art_sector_layer", s_side*8*(s_sector-1)+art_layer);
 
     const int rollover=3564;
     int art_bc=art_bcids[i];
     int relative = art_bc - event_bcid;;
     if (relative > rollover / 2) {
       relative -= rollover;
     } else if (relative <= -rollover / 2) {
       relative += rollover;
     }
     if (relative > (rollover-2048) / 2) {
       relative -= (rollover-2048);
     }
     auto art_deltaBC=Monitored::Scalar<int>("art_deltaBC",relative);
     auto art_deltaBC_perSector=Monitored::Scalar<int>("art_deltaBC_"+MM_Side[iside]+"_s"+std::to_string(s_sector),relative);
     auto art_bc_mon=Monitored::Scalar<int>("art_bc",art_bc);
     
     fill("mmTrigger", art_channel, trig_sector, art_sector_layer, art_deltaBC, art_bc_mon, bcid, lb_tri, art_deltaBC_perSector);
       }
 
       auto bcids =rdo->trig_BCID();
       auto dthetas=rdo->trig_dTheta();
       auto rids=rdo->trig_ROI_rID();
       auto phiids=rdo->trig_ROI_phiID() ;
 
 
       std::unordered_map<int, int> NROIPerBC;
       for (int bc : bcids)      NROIPerBC[bc]++;
       for (const auto& [value, count] : NROIPerBC) {
     auto nROIPerBC=Monitored::Scalar<int>("nROIPerBC",count);
     fill("mmTrigger_roi", nROIPerBC, trig_sector);
       }
 
       int nROI=rdo->trig_BCID().size();
       for (int i=0; i< nROI; i++ ){
     auto phiID = (phiids[i] & 0b11111) * ((phiids[i] >> 5) ? 1 : -1);
     int sign = phiID > 0 ? 1 : -1;
     auto phi_conv = (phiID-0.5*sign)*(16./31.)*M_PI/180. + sector_pos; 
     if(phi_conv> M_PI)phi_conv = phi_conv - 2*M_PI;
 
     const float z_ref=7824.46;
     const float r_step=(5000-900)/256.;
     auto rID = static_cast<unsigned int>(rids[i]);
     float r_conv=r_step*rID+900;
     float eta_conv=-log(0.5*atan(r_conv/z_ref))*s_side;
     //  auto dTheta=Monitored::Scalar<float>("dTheta_roi",static_cast<float>(dthetas[i])); 
     auto deltaBC=   Monitored::Scalar<int>("deltaBC", bcids[i]-event_bcid);
     auto deltaBC_perSector=   Monitored::Scalar<int>("deltaBC_"+MM_Side[iside]+"_s"+std::to_string(s_sector), bcids[i]-event_bcid);
     auto rid=Monitored::Scalar<int>("rid",rID);
     auto phiid=Monitored::Scalar<int>("phiid",phiID);
     auto rid_sector=Monitored::Scalar<int>("rid_"+MM_Side[iside]+"_s"+std::to_string(s_sector),rID);
     auto phiid_sector=Monitored::Scalar<int>("phiid_"+MM_Side[iside]+"_s"+std::to_string(s_sector),phiID);
     auto r_roi=Monitored::Scalar<float>("r_roi",r_conv);
     auto phi_roi=Monitored::Scalar<float>("phi_roi",phi_conv);
     auto eta_roi=Monitored::Scalar<float>("eta_roi",eta_conv);
     auto x_roi_sideA=Monitored::Scalar<float>("x_roi_sideA", r_conv*cos(phi_roi));
     auto y_roi_sideA=Monitored::Scalar<float>("y_roi_sideA", r_conv*sin(phi_roi));
     auto x_roi_sideC=Monitored::Scalar<float>("x_roi_sideC", r_conv*cos(phi_roi) );
     auto y_roi_sideC=Monitored::Scalar<float>("y_roi_sideC", r_conv*sin(phi_roi));
 
     fill("mmTrigger_roi", deltaBC, phiid, rid, trig_sector, phi_roi, eta_roi,r_roi, lb_tri, rid_sector, phiid_sector, deltaBC_perSector);
     if(s_side>0) fill("mmTrigger_roi", x_roi_sideA, y_roi_sideA);
     if(s_side<0) fill("mmTrigger_roi", x_roi_sideC, y_roi_sideC);
 
     }
  }
  
}

filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

get_bin_for_feb_occ()

int MMRawDataMonAlg::get_bin_for_feb_occ ( const int gas_gap, const int FEB ) const

private

Definition at line 136 of file MMRawDataUtils.cxx.

                                                                   {
 
  static const int max_feb = 16;
 
  return  (gas_gap-1)*max_feb + (FEB-1);
 
}

get_bin_for_occ()

int MMRawDataMonAlg::get_bin_for_occ ( const int gas_gap, const int PCB ) const

private

Definition at line 128 of file MMRawDataUtils.cxx.

                                                               {
 
  static const int max_pcb = 8;
 
  return  (gas_gap-1)*max_pcb + (PCB-1);
 
}

get_bin_for_occ_ASide_hist()

int MMRawDataMonAlg::get_bin_for_occ_ASide_hist ( const int stationEta, const int multiplet, const int gas_gap ) const

private

Definition at line 106 of file MMRawDataUtils.cxx.

                                                                                                {
 
  static const int max_gas_gap = 4;
  static const int max_multiplet = 2;
 
  return (stationEta-1)*(max_gas_gap*max_multiplet)+(multiplet-1)*max_gas_gap +(gas_gap-1);
 
}

get_bin_for_occ_ASide_pcb_eta1_hist()

int MMRawDataMonAlg::get_bin_for_occ_ASide_pcb_eta1_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB ) const

private

Definition at line 167 of file MMRawDataUtils.cxx.

                                                                                                                  {
 
  static const int max_pcb = 5;
  static const int max_gas_gap = 4;
  if (stationEta != 1) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb+ (gas_gap-1)*max_pcb + (PCB-1);
}

get_bin_for_occ_ASide_pcb_eta2_hist()

int MMRawDataMonAlg::get_bin_for_occ_ASide_pcb_eta2_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB ) const

private

Definition at line 156 of file MMRawDataUtils.cxx.

                                                                                                                  {
 
  static const int max_pcb = 3;
  static const int max_gas_gap = 4;
  if (stationEta != 2) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb+ (gas_gap-1)*max_pcb + (PCB-1);
 
}

get_bin_for_occ_CSide_hist()

int MMRawDataMonAlg::get_bin_for_occ_CSide_hist ( const int stationEta, const int multiplet, const int gas_gap ) const

private

Definition at line 97 of file MMRawDataUtils.cxx.

                                                                                                {
 
  static const int max_gas_gap = 4;
  static const int max_multiplet = 2;
 
  return (stationEta+2)*(max_gas_gap*max_multiplet)+(multiplet-1)*max_gas_gap +(gas_gap-1);
 
}

get_bin_for_occ_CSide_pcb_eta1_hist()

int MMRawDataMonAlg::get_bin_for_occ_CSide_pcb_eta1_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB ) const

private

Definition at line 145 of file MMRawDataUtils.cxx.

                                                                                                                  {
 
  static const int max_pcb = 5;
  static const int max_gas_gap = 4;
  if (stationEta != -1) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb+ (gas_gap-1)*max_pcb + (PCB-1);
 
}

get_bin_for_occ_CSide_pcb_eta2_hist()

int MMRawDataMonAlg::get_bin_for_occ_CSide_pcb_eta2_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB ) const

private

Definition at line 116 of file MMRawDataUtils.cxx.

                                                                                                                  {
 
  static const int max_pcb = 3;
  static const int max_gas_gap = 4;
  if (stationEta != -2) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb+ (gas_gap-1)*max_pcb + (PCB-1);
 
}

get_bin_for_occ_lb_ASide_pcb_eta1_hist()

int MMRawDataMonAlg::get_bin_for_occ_lb_ASide_pcb_eta1_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB, int isector ) const

private

Definition at line 200 of file MMRawDataUtils.cxx.

                                                                                                                                 {
 
  static const int max_pcb = 5;
  static const int max_gas_gap = 4;
  static const int max_isector = 2;
  if (stationEta != 1) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb*max_isector+ (gas_gap-1)*max_pcb*max_isector + isector*max_pcb + (PCB-1);
}

get_bin_for_occ_lb_ASide_pcb_eta2_hist()

int MMRawDataMonAlg::get_bin_for_occ_lb_ASide_pcb_eta2_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB, const int isector ) const

private

Definition at line 211 of file MMRawDataUtils.cxx.

                                                                                                                                 {
 
  static const int max_pcb = 3;
  static const int max_gas_gap = 4;
  static const int max_isector = 2;
  if (stationEta != 2) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb*max_isector+ (gas_gap-1)*max_pcb*max_isector + isector*max_pcb + (PCB-1);
 
}

get_bin_for_occ_lb_CSide_pcb_eta1_hist()

int MMRawDataMonAlg::get_bin_for_occ_lb_CSide_pcb_eta1_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB, int isector ) const

private

Definition at line 189 of file MMRawDataUtils.cxx.

                                                                                                                                 {
 
  static const int max_pcb = 5;
  static const int max_gas_gap = 4;
  static const int max_isector = 2;
  if (stationEta != -1) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb*max_isector+ (gas_gap-1)*max_pcb*max_isector + isector*max_pcb + (PCB-1);
 
}

get_bin_for_occ_lb_CSide_pcb_eta2_hist()

int MMRawDataMonAlg::get_bin_for_occ_lb_CSide_pcb_eta2_hist ( const int stationEta, const int multiplet, const int gas_gap, const int PCB, const int isector ) const

private

Definition at line 177 of file MMRawDataUtils.cxx.

                                                                                                                                 {
 
  static const int max_pcb = 3;
  static const int max_gas_gap = 4;
  static const int max_isector = 2;
  if (stationEta != -2) return -1;
 
  return  (multiplet-1)*max_gas_gap*max_pcb*max_isector+ (gas_gap-1)*max_pcb*max_isector + isector*max_pcb+ (PCB-1);
 
}

get_bin_for_occ_lb_pcb_hist()

int MMRawDataMonAlg::get_bin_for_occ_lb_pcb_hist ( const int multiplet, const int gas_gap, const int PCB ) const

private

Definition at line 222 of file MMRawDataUtils.cxx.

                                                                                          {
 
  static const int max_pcb = 8;
  static const int max_gas_gap = 4;
  
  return (multiplet-1)*max_gas_gap*max_pcb + (gas_gap-1)*max_pcb + (PCB-1);
 
}

get_FEB_from_channel()

int MMRawDataMonAlg::get_FEB_from_channel ( const int channel, const int stEta ) const

private

Definition at line 34 of file MMRawDataUtils.cxx.

                                                                      {
 
  if(fabs(stEta)==1){
    if (channel>0 && channel<=512) return 1;
    if (channel>512 && channel<=1024) return 2;
    
    if (channel>1024 && channel<=1536) return 3;
    if (channel>1536 && channel<=2048) return 4;
 
    if (channel>2048 && channel<=2560) return 5;
    if (channel>2560 && channel<=3072) return 6;
    
    if (channel>3072 && channel<=3584) return 7;
    if (channel>3584 && channel<=4096) return 8;
    
    if (channel>4096 && channel<=4608) return 9;
    if (channel>4608 && channel<=5120) return 10;
  }else if(fabs(stEta)==2){
    if (channel>0 && channel<=512) return 11;
    if (channel>512 && channel<=1024) return 12;
    
    if (channel>1024 && channel<=1536) return 13;
    if (channel>1536 && channel<=2048) return 14;
 
    if (channel>2048 && channel<=2560) return 15;
    if (channel>2560 && channel<=3072) return 16;
  } 
  throw std::invalid_argument( "channel is not valid!" );
}

get_PCB_from_channel()

int MMRawDataMonAlg::get_PCB_from_channel ( const int channel ) const

private

Definition at line 23 of file MMRawDataUtils.cxx.

                                                           {
 
  if (channel>0 && channel<=1024) return 1; 
  if (channel>1024 && channel<=2048) return 2;
  if (channel>2048 && channel<=3072) return 3;
  if (channel>3072 && channel<=4096) return 4;
  if (channel>4096 && channel<=5120) return 5;
 
  throw std::invalid_argument( "channel is not valid!" );
}

get_sectorEta_from_stationEta()

int MMRawDataMonAlg::get_sectorEta_from_stationEta ( const int stationEta ) const

private

Definition at line 87 of file MMRawDataUtils.cxx.

                                                                       {
 
  //  1<-0  0-> 1
  if (std::abs(stationEta)==1) return 0;                                                            
  if (std::abs(stationEta)==2) return 1;
 
  return -1;
 
}

get_sectorPhi_from_stationPhi_stName()

int MMRawDataMonAlg::get_sectorPhi_from_stationPhi_stName ( const int stationPhi, const std::string & stName ) const

private

Definition at line 64 of file MMRawDataUtils.cxx.

                                                                                                       {
  
  if (stationPhi==1 && stName=="MML") return 1;
  if (stationPhi==1 && stName=="MMS") return 2;
  if (stationPhi==2 && stName=="MML") return 3;
  if (stationPhi==2 && stName=="MMS") return 4;
  if (stationPhi==3 && stName=="MML") return 5;
  if (stationPhi==3 && stName=="MMS") return 6;
  if (stationPhi==4 && stName=="MML") return 7;
  if (stationPhi==4 && stName=="MMS") return 8;
  if (stationPhi==5 && stName=="MML") return 9;
  if (stationPhi==5 && stName=="MMS") return 10;
  if (stationPhi==6 && stName=="MML") return 11;
  if (stationPhi==6 && stName=="MMS") return 12;
  if (stationPhi==7 && stName=="MML") return 13;
  if (stationPhi==7 && stName=="MMS") return 14;
  if (stationPhi==8 && stName=="MML") return 15;
  if (stationPhi==8 && stName=="MMS") return 16;
 
  throw std::invalid_argument( "stationPhi and stName are not valid!" );
 
}

GetEventInfo()

SG::ReadHandle< xAOD::EventInfo > AthMonitorAlgorithm::GetEventInfo ( const EventContext & ctx ) const

inherited

Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)

Parameters

ctx	EventContext for the event

Returns

a SG::ReadHandle<xAOD::EventInfo>

Definition at line 111 of file AthMonitorAlgorithm.cxx.

                                                                                             {
    return SG::ReadHandle<xAOD::EventInfo>(m_EventInfoKey, ctx);
}

getGroup()

const ToolHandle< GenericMonitoringTool > & AthMonitorAlgorithm::getGroup ( const std::string & name ) const

inherited

Get a specific monitoring tool from the tool handle array.

Finds a specific GenericMonitoringTool instance from the list of monitoring tools (a ToolHandleArray). Throws a FATAL warning if the object found is empty.

Parameters

name	string name of the desired tool

Returns

reference to the desired monitoring tool

Definition at line 168 of file AthMonitorAlgorithm.cxx.

                                                                                                    {
    // get the pointer to the tool, and check that it exists
    auto idx = m_toolLookupMap.find(name);
    if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
        return m_tools[idx->second];
    }
    else {
      // treat empty tool handle case as in Monitored::Group
      if (m_toolLookupMap.empty()) {
    return m_dummy;
      }
 
      if (!isInitialized()) {
        ATH_MSG_FATAL(
            "It seems that the AthMonitorAlgorithm::initialize was not called "
            "in derived class initialize method");
      } else {
        std::string available = std::accumulate(
            m_toolLookupMap.begin(), m_toolLookupMap.end(), std::string(""),
            [](const std::string& s, auto h) { return s + "," + h.first; });
        ATH_MSG_FATAL("The tool " << name << " could not be found in the tool array of the "
                      << "monitoring algorithm " << m_name << ". This probably reflects a discrepancy between "
                      << "your python configuration and c++ filling code. Note: your available groups are {"
                      << available << "}.");
        }
    }
    return m_dummy;
}

getTrigDecisionTool()

const ToolHandle< Trig::TrigDecisionTool > & AthMonitorAlgorithm::getTrigDecisionTool ( ) const

inherited

Get the trigger decision tool member.

The trigger decision tool is used to check whether a specific trigger is passed by an event.

Returns

m_trigDecTool

Definition at line 198 of file AthMonitorAlgorithm.cxx.

                                                                                       {
    return m_trigDecTool;
}

initialize()

StatusCode MMRawDataMonAlg::initialize ( )

overridevirtual

initialize

Returns

StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 103 of file MMRawDataMonAlg.cxx.

{
    //init message stream
    ATH_MSG_DEBUG("initialize MMRawDataMonAlg");
    ATH_MSG_DEBUG("******************");
    ATH_MSG_DEBUG("doMMESD: " << m_doMMESD );
    ATH_MSG_DEBUG("******************");
 
    ATH_CHECK(AthMonitorAlgorithm::initialize());
    ATH_CHECK(m_DetectorManagerKey.initialize());
    ATH_CHECK(m_idHelperSvc.retrieve());
 
    ATH_MSG_INFO(" Found the MuonIdHelperSvc ");
    ATH_CHECK(m_muonKey.initialize());
    ATH_CHECK(m_MMContainerKey.initialize());
    ATH_CHECK(m_meTrkKey.initialize());
        ATH_CHECK(m_segm_type.initialize());
    ATH_CHECK(m_mmtpRdoKey.initialize());
 
    ATH_MSG_DEBUG(" end of initialize " );
    ATH_MSG_INFO("MMRawDataMonAlg initialization DONE " );
 
    return StatusCode::SUCCESS;
} 

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

MMEfficiency()

void MMRawDataMonAlg::MMEfficiency ( const xAOD::TrackParticleContainer * muonContainer ) const

private

Definition at line 747 of file MMRawDataMonAlg.cxx.

{
    MMEfficiencyHistogramStruct effPlots[2][2][16][2][4];
    MMEfficiencyHistogramStruct Gaps[2][2][16][2];
 
    static const std::array<std::string,2> MM_Side = {"CSide", "ASide"};
    static const std::array<std::string,2> EtaSector = {"1","2"};
 
    for (const xAOD::TrackParticle* meTP  : *muonContainer) {
        if (!meTP) continue;
        auto eta_trk = Monitored::Scalar<float>("eta_trk", meTP->eta());
        auto phi_trk = Monitored::Scalar<float>("phi_trk", meTP->phi());
 
        float pt_trk = meTP->pt();
        if(pt_trk < m_cut_pt) continue;
        // retrieve the original track
        const Trk::Track* meTrack = meTP->track();
        if(!meTrack) continue;
        // get the vector of measurements on track
        const DataVector<const Trk::MeasurementBase>* meas = meTrack->measurementsOnTrack();
 
        for(const Trk::MeasurementBase* it: *meas) {
            const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(it);
            if (!rot) continue;
            Identifier rot_id = rot->identify();
            if (!m_idHelperSvc->isMM(rot_id)) continue;
 
            const Muon::MMClusterOnTrack* cluster = dynamic_cast<const Muon::MMClusterOnTrack*>(rot);
            if (!cluster) continue;
            std::string stName   = m_idHelperSvc->mmIdHelper().stationNameString(m_idHelperSvc->mmIdHelper().stationName(rot_id));
            int stEta= m_idHelperSvc->mmIdHelper().stationEta(rot_id);
            int stPhi= m_idHelperSvc->mmIdHelper().stationPhi(rot_id);
            int phi = get_sectorPhi_from_stationPhi_stName(stPhi,stName);
            int multi = m_idHelperSvc->mmIdHelper().multilayer(rot_id);
            int gap=  m_idHelperSvc->mmIdHelper().gasGap(rot_id);
            int ch=  m_idHelperSvc->mmIdHelper().channel(rot_id);
            int pcb=get_PCB_from_channel(ch);
            int abs_stEta= get_sectorEta_from_stationEta(stEta);
            int iside = (stEta > 0) ? 1 : 0;
                        if( ! (std::find( Gaps[iside][abs_stEta][phi-1][multi-1].nGaps.begin(), Gaps[iside][abs_stEta][phi-1][multi-1].nGaps.end(), gap ) != Gaps[iside][abs_stEta][phi-1][multi-1].nGaps.end()) )
                          Gaps[iside][abs_stEta][phi-1][multi-1].nGaps.push_back(gap);
            //numerator
            if(effPlots[iside][abs_stEta][phi-1][multi-1][gap-1].num.size()==0) effPlots[iside][abs_stEta][phi-1][multi-1][gap-1].num.push_back(pcb-1);
        }
    } // loop on tracks
 
    unsigned  int nGaptag=3;
 
    for(int s=0; s<2; ++s) {
        std::string MM_sideGroup = "MM_sideGroup"+MM_Side[s];
        for(int e=0; e<2; ++e) {
            for(int p=0; p<16; ++p) {
                for(int m=0; m<2; ++m) {
                    if(Gaps[s][e][p][m].nGaps.size()<nGaptag) continue;
                                        if(Gaps[s][e][p][m].nGaps.size()>4) continue;
                                        //find the missing gap                                                                                                                                                                         
                                        int gapsum=0;
                    for (unsigned int g=0; g<Gaps[s][e][p][m].nGaps.size(); ++g)
                      gapsum+= Gaps[s][e][p][m].nGaps.at(g);
                    int missing_gap=10-gapsum-1;
                    //if missing gap = -1 --> nGaps=4 --> all efficient                                                                                                                                           
                    if(Gaps[s][e][p][m].nGaps.size()==4){
                      for (unsigned int ga=0; ga<Gaps[s][e][p][m].nGaps.size(); ++ga){
                        for (unsigned int i=0; i<effPlots[s][e][p][m][ga].num.size(); ++i){
                          int pcb = effPlots[s][e][p][m][ga].num.at(i);
                          auto traversed_pcb = Monitored::Scalar<int>("pcb_eta"+std::to_string(e+1)+"_"+MM_Side[s]+"_phi"+std::to_string(p)+"_multiplet"+std::to_string(m+1)+"_gas_gap"+std::to_string(ga+1),pcb);
                          int layer=ga+4*m+8*p;
                          
                          auto traversed_gap = Monitored::Scalar<int>(MM_Side[s]+"_eta"+std::to_string(e+1),layer);
                          auto isHit = 1;
                          auto hitcut = Monitored::Scalar<int>("hitcut", (int)isHit);
                          fill(MM_sideGroup, traversed_pcb, traversed_gap, hitcut);
                        }
                      }
                    } else {// 3 gaps, the fourth is inefficient                                                                                                                                                  
                      int ref_gap = missing_gap+1;
                      if(missing_gap==3) ref_gap=0;
                      if (ref_gap>3){
                        ATH_MSG_FATAL("ref_gap is out of range in MMRawDataMonAlg::MMEfficiency");
                        return;
                      }
                      int ref_pcb=effPlots[s][e][p][m][ref_gap].num.at(0);
                      auto traversed_pcb = Monitored::Scalar<int>("pcb_eta"+std::to_string(e+1)+"_"+MM_Side[s]+"_phi"+std::to_string(p)+"_multiplet"+std::to_string(m+1)+"_gas_gap"+std::to_string(missing_gap+1), ref_pcb);
                      int layer=missing_gap+4*m+8*p;
                      auto traversed_gap = Monitored::Scalar<int>(MM_Side[s]+"_eta"+std::to_string(e+1),layer);
                      auto isHit = 0;
                      auto hitcut = Monitored::Scalar<int>("hitcut", (int)isHit);
                      fill(MM_sideGroup, traversed_pcb, traversed_gap, hitcut);
                    }
 
 
                }
            }
        }
    }
}

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

parseList()

StatusCode AthMonitorAlgorithm::parseList ( const std::string & line, std::vector< std::string > & result ) const

virtualinherited

Parse a string into a vector.

The input string is a single long string of all of the trigger names. It parses this string and turns it into a vector, where each element is one trigger or trigger category.

Parameters

line	The input string.
result	The parsed output vector of strings.

Returns

StatusCode

Definition at line 345 of file AthMonitorAlgorithm.cxx.

                                                                                                   {
    std::string item;
    std::stringstream ss(line);
 
    ATH_MSG_DEBUG( "AthMonitorAlgorithm::parseList()" );
 
    while ( std::getline(ss, item, ',') ) {
        std::stringstream iss(item); // remove whitespace
        iss >> item;
        result.push_back(item);
    }
 
    return StatusCode::SUCCESS;
}

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

                                                               {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

trigChainsArePassed()

bool AthMonitorAlgorithm::trigChainsArePassed ( const std::vector< std::string > & vTrigNames ) const

inherited

Check whether triggers are passed.

For the event, use the trigger decision tool to check that at least one of the triggers listed in the supplied vector is passed.

Parameters

vTrigNames

List of trigger names.

Returns

If empty input, default to true. If at least one trigger is specified, returns whether at least one trigger was passed.

Definition at line 203 of file AthMonitorAlgorithm.cxx.

                                                                                            {
 
  
  // If no triggers were given, return true.
  if (vTrigNames.empty()) return true;
  
  
  // Trigger: Check if this Algorithm is being run as an Express Stream job.
  // Events are entering the express stream are chosen randomly, and by chain,
  // Hence an additional check should be aplied to see if the chain(s)
  // monitored here are responsible for the event being selected for
  // the express stream.
 
  const auto group =  m_trigDecTool->getChainGroup(vTrigNames);
  if (m_enforceExpressTriggers){  
    const auto passedBits = m_trigDecTool->isPassedBits(group);
    bool expressPass = passedBits & TrigDefs::Express_passed; //bitwise AND
    if(!expressPass) {
      return false;
    }
  }
  
  // monitor the event if any of the chains in the chain group passes the event.
  return group->isPassed();
  
}

updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_cut_pt

Gaudi::Property<float> MMRawDataMonAlg::m_cut_pt {this,"cut_pt",15000}

private

Definition at line 109 of file MMRawDataMonAlg.h.

{this,"cut_pt",15000};

m_dataType

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::m_dataType

protectedinherited

Instance of the DataType_t enum.

Definition at line 356 of file AthMonitorAlgorithm.h.

m_dataTypeStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_dataTypeStr {this,"DataType","userDefined"}

protectedinherited

DataType string pulled from the job option and converted to enum.

Definition at line 358 of file AthMonitorAlgorithm.h.

{this,"DataType","userDefined"}; 

m_defaultLBDuration

Gaudi::Property<float> AthMonitorAlgorithm::m_defaultLBDuration {this,"DefaultLBDuration",60.}

protectedinherited

Default duration of one lumi block.

Definition at line 365 of file AthMonitorAlgorithm.h.

{this,"DefaultLBDuration",60.}; 

m_detailLevel

Gaudi::Property<int> AthMonitorAlgorithm::m_detailLevel {this,"DetailLevel",0}

protectedinherited

Sets the level of detail used in the monitoring.

Definition at line 366 of file AthMonitorAlgorithm.h.

{this,"DetailLevel",0}; 

m_DetectorManagerKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> MMRawDataMonAlg::m_DetectorManagerKey {this, "DetectorManagerKey", "MuonDetectorManager","Key of input MuonDetectorManager condition data"}

private

Definition at line 65 of file MMRawDataMonAlg.h.

{this, "DetectorManagerKey", "MuonDetectorManager","Key of input MuonDetectorManager condition data"};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_do_mm_overview

Gaudi::Property<bool> MMRawDataMonAlg::m_do_mm_overview {this,"do_mm_overview",true}

private

Definition at line 107 of file MMRawDataMonAlg.h.

{this,"do_mm_overview",true};

m_do_stereoCorrection

Gaudi::Property<bool> MMRawDataMonAlg::m_do_stereoCorrection {this,"do_stereoCorrection",false}

private

Definition at line 108 of file MMRawDataMonAlg.h.

{this,"do_stereoCorrection",false};

m_doDetailedHists

Gaudi::Property<bool> MMRawDataMonAlg::m_doDetailedHists {this,"doDetailedHists",true}

private

Definition at line 110 of file MMRawDataMonAlg.h.

{this,"doDetailedHists",true};

m_doMMESD

Gaudi::Property<bool> MMRawDataMonAlg::m_doMMESD {this,"DoMMESD",true}

private

Definition at line 106 of file MMRawDataMonAlg.h.

{this,"DoMMESD",true};

m_DQFilterTools

ToolHandleArray<IDQFilterTool> AthMonitorAlgorithm::m_DQFilterTools {this,"FilterTools",{}}

protectedinherited

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

{this,"FilterTools",{}}; 

m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy

privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

m_enforceExpressTriggers

Gaudi::Property<bool> AthMonitorAlgorithm::m_enforceExpressTriggers

privateinherited

Initial value:

{this,
                          "EnforceExpressTriggers", false,
                          "Requires that matched triggers made the event enter the express stream"}

Definition at line 377 of file AthMonitorAlgorithm.h.

                                               {this,
                          "EnforceExpressTriggers", false,
                          "Requires that matched triggers made the event enter the express stream"};

m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment

protectedinherited

Instance of the Environment_t enum.

Definition at line 355 of file AthMonitorAlgorithm.h.

m_environmentStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_environmentStr {this,"Environment","user"}

protectedinherited

Environment string pulled from the job option and converted to enum.

Definition at line 357 of file AthMonitorAlgorithm.h.

{this,"Environment","user"}; 

m_EventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> AthMonitorAlgorithm::m_EventInfoKey {this,"EventInfoKey","EventInfo"}

protectedinherited

Key for retrieving EventInfo from StoreGate.

Definition at line 367 of file AthMonitorAlgorithm.h.

{this,"EventInfoKey","EventInfo"}; 

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_fileKey

Gaudi::Property<std::string> AthMonitorAlgorithm::m_fileKey {this,"FileKey",""}

protectedinherited

Internal Athena name for file.

Definition at line 363 of file AthMonitorAlgorithm.h.

{this,"FileKey",""}; 

m_idHelperSvc

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

private

Definition at line 62 of file MMRawDataMonAlg.h.

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

m_lbDurationDataKey

SG::ReadCondHandleKey<LBDurationCondData> AthMonitorAlgorithm::m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}

protectedinherited

Definition at line 350 of file AthMonitorAlgorithm.h.

{this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"};

m_lumiDataKey

SG::ReadCondHandleKey<LuminosityCondData> AthMonitorAlgorithm::m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}

protectedinherited

Definition at line 348 of file AthMonitorAlgorithm.h.

{this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"};

m_meTrkKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> MMRawDataMonAlg::m_meTrkKey {this, "METrkContainer", "ExtrapolatedMuonTrackParticles"}

private

Definition at line 69 of file MMRawDataMonAlg.h.

{this, "METrkContainer", "ExtrapolatedMuonTrackParticles"};

m_MMContainerKey

SG::ReadHandleKey<Muon::MMPrepDataContainer> MMRawDataMonAlg::m_MMContainerKey {this,"MMPrepDataContainerName","MM_Measurements"}

private

Definition at line 67 of file MMRawDataMonAlg.h.

{this,"MMPrepDataContainerName","MM_Measurements"};

m_mmtpRdoKey

SG::ReadHandleKey<xAOD::NSWMMTPRDOContainer> MMRawDataMonAlg::m_mmtpRdoKey {this, "NSW_MMTPDataKey", "NSW_MMTrigProcessor_RDO"}

private

Definition at line 70 of file MMRawDataMonAlg.h.

{this, "NSW_MMTPDataKey", "NSW_MMTrigProcessor_RDO"};

m_muonKey

SG::ReadHandleKey<xAOD::MuonContainer> MMRawDataMonAlg::m_muonKey {this,"MuonKey","Muons","muons"}

private

Definition at line 68 of file MMRawDataMonAlg.h.

{this,"MuonKey","Muons","muons"};

m_muonSelectionTool

ToolHandle<CP::IMuonSelectionTool> MMRawDataMonAlg::m_muonSelectionTool {this,"MuonSelectionTool","CP::MuonSelectionTool/MuonSelectionTool"}

private

Definition at line 64 of file MMRawDataMonAlg.h.

{this,"MuonSelectionTool","CP::MuonSelectionTool/MuonSelectionTool"};

m_name

std::string AthMonitorAlgorithm::m_name

privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

m_segm_type

SG::ReadHandleKey<Trk::SegmentCollection> MMRawDataMonAlg::m_segm_type {this,"Eff_segm_type","TrackMuonSegments","muon segments"}

private

Definition at line 66 of file MMRawDataMonAlg.h.

{this,"Eff_segm_type","TrackMuonSegments","muon segments"};

m_toolLookupMap

std::unordered_map<std::string, size_t> AthMonitorAlgorithm::m_toolLookupMap

privateinherited

Definition at line 372 of file AthMonitorAlgorithm.h.

m_tools

ToolHandleArray<GenericMonitoringTool> AthMonitorAlgorithm::m_tools {this,"GMTools",{}}

protectedinherited

Array of Generic Monitoring Tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

{this,"GMTools",{}}; 

m_trigDecTool

PublicToolHandle<Trig::TrigDecisionTool> AthMonitorAlgorithm::m_trigDecTool

protectedinherited

Tool to tell whether a specific trigger is passed.

Definition at line 345 of file AthMonitorAlgorithm.h.

m_triggerChainString

Gaudi::Property<std::string> AthMonitorAlgorithm::m_triggerChainString {this,"TriggerChain",""}

protectedinherited

Trigger chain string pulled from the job option and parsed into a vector.

Definition at line 360 of file AthMonitorAlgorithm.h.

{this,"TriggerChain",""}; 

m_trigLiveFractionDataKey

SG::ReadCondHandleKey<TrigLiveFractionCondData> AthMonitorAlgorithm::m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}

protectedinherited

Definition at line 352 of file AthMonitorAlgorithm.h.

{this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"};

m_useLumi

Gaudi::Property<bool> AthMonitorAlgorithm::m_useLumi {this,"EnableLumi",false}

protectedinherited

Allows use of various luminosity functions.

Definition at line 364 of file AthMonitorAlgorithm.h.

{this,"EnableLumi",false}; 

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_vTrigChainNames

std::vector<std::string> AthMonitorAlgorithm::m_vTrigChainNames

protectedinherited

Vector of trigger chain names parsed from trigger chain string.

Definition at line 361 of file AthMonitorAlgorithm.h.

---

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

• MMRawDataMonAlg.h
• MMRawDataMonAlg.cxx
• MMRawDataUtils.cxx