L2MuonSAIOMon Class Reference

#include <L2MuonSAIOMon.h>

Inheritance diagram for L2MuonSAIOMon:

Collaboration diagram for L2MuonSAIOMon:

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
	L2MuonSAIOMon (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
	initialize
virtual StatusCode	fillHistograms (const EventContext &ctx) const override
	Function that steers anlayses.
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) const
	Calculate the average mu, i.e.
virtual float	lbInteractionsPerCrossing (const EventContext &ctx) const
	Calculate instantaneous number of interactions, i.e.
virtual float	lbAverageLuminosity (const EventContext &ctx) const
	Calculate average luminosity (in ub-1 s-1 => 10^30 cm-2 s-1).
virtual float	lbLuminosityPerBCID (const EventContext &ctx) const
	Calculate the instantaneous luminosity per bunch crossing.
virtual double	lbDuration (const EventContext &ctx) const
	Calculate the duration of the luminosity block (in seconds).
virtual float	lbAverageLivefraction (const EventContext &ctx) const
	Calculate the average luminosity livefraction.
virtual float	livefractionPerBCID (const EventContext &ctx) const
	Calculate the live fraction per bunch crossing ID.
virtual double	lbLumiWeight (const EventContext &ctx) 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
virtual StatusCode	fillVariablesPerChain (const EventContext &ctx, const std::string &chain) const override
	Function that fills variables of trigger objects associated to specified trigger chains.
virtual StatusCode	fillVariablesPerOfflineMuonPerChain (const EventContext &ctx, const xAOD::Muon *mu, const std::string &chain) const override
	Function that fills variables of trigger objects associated to specified trigger chains comparing offline muons.
virtual bool	selectEvents () const
	Function that defines the event selection for anlayses User should reimlement in a subclass if needed.
virtual StatusCode	selectMuons (const EventContext &ctx, SG::ReadHandle< xAOD::MuonContainer > &muons, std::vector< const xAOD::Muon * > &probes) const
	Function that defines the event selection for anlayses Users should reimlement in a subclass if needed.
virtual StatusCode	fillVariables (const EventContext &ctx) const
	Function that fills variables by just retrieving containers of trigger objects.
virtual StatusCode	fillVariablesPerOfflineMuon (const EventContext &ctx, const xAOD::Muon *mu) const
	Function that fills variables that are compared to offline muons but the trigger chains are not specified.
template<class T, class FUNCT>
StatusCode	fillVariablesRatioPlots (const EventContext &ctx, const xAOD::Muon *mu, std::string &&trigstep, xAOD::Muon::TrackParticleType type, FUNCT matchFunc) const
	Function that fills variables of ratio plots.
template<class T>
StatusCode	fillVariableEtaPhi (const EventContext &ctx, SG::ReadHandleKey< DataVector< T > > ReadHandleKey, std::string &&trigstep, std::tuple< bool, double, double >(*PosFunc)(const T *)=&TrigMuonMonitorAlgorithm::defaultPosFunc< T >) const
	Function that fills variables of etaphi2D plots.
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.

Static Protected Member Functions
template<class T>
static std::tuple< bool, double, double >	defaultPosFunc (const T *trig)

Protected Attributes
ToolHandle< MuonMatchingTool >	m_matchTool {this, "MuonMatchingTool", "MuonMatchingTool", "Tool for matching offline and online objects"}
ToolHandle< CP::IMuonSelectionTool >	m_muonSelectionTool {this,"MuonSelectionTool","CP::MuonSelectionTool/MuonSelectionTool","MuonSelectionTool"}
SG::ReadHandleKey< xAOD::MuonContainer >	m_MuonContainerKey {this, "MuonContainerName", "Muons", "Offline muon container"}
Gaudi::Property< std::vector< std::string > >	m_monitored_chains {this, "MonitoredChains", {}, "Trigger chains that are monitored"}
	List of trigger chains that are monitored in fillVariablesPerChain and fillVariablesPerOfflineMuonPerChain.
Gaudi::Property< int >	m_muontype
	Requirement for the offline muon type considered in analyses.
Gaudi::Property< std::string >	m_group {this, "Group", "", "Histogram group"}
	Name of monitored group.
const float	m_ratio_measurement_threshold = 4
	Threshold for ratio measurement.
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	matchL2IO_wContainer (const EventContext &ctx, const std::string &chain, std::vector< const xAOD::L2CombinedMuon * > &Trig_L2IOobjects) const
StatusCode	L2OverlapRemover (const std::vector< const xAOD::L2CombinedMuon * > &matchSA_L2IOobjects, std::vector< bool > &isoverlap, std::vector< bool > &passOR) const
bool	isOverlap (const xAOD::L2CombinedMuon *matchSA_L2IOobject1, const xAOD::L2CombinedMuon *matchSA_L2IOobject2) const
StatusCode	chooseBestMuon (const std::vector< const xAOD::L2CombinedMuon * > &matchSA_L2IOobjects, std::vector< bool > &passOR, std::vector< unsigned int > &mucombResult) const
StatusCode	muCombHypo_TDTworkaround (const std::string &chain, const std::vector< const xAOD::L2CombinedMuon * > &Trig_L2IOobjects, std::vector< bool > &pass_muCombHypo) const
StatusCode	isPassedmuCombHypo (const std::string &chain, const xAOD::L2CombinedMuon *Trig_L2IOobjects, bool &pass_muCombHypo) const
StatusCode	decision_ptthreshold (const std::string &chain, std::vector< float > &my_EtaBins, std::vector< float > &my_muCombThres, bool &my_pikCuts, float &my_maxPtToApplyPik, float &my_chi2MaxID) const
const xAOD::L2CombinedMuon *	searchL2InsideOut (const EventContext &ctx, const xAOD::Muon *mu, const std::string &trigger) const
const xAOD::Muon *	searchTagOfflineMuon (const EventContext &ctx, const xAOD::Muon *probe) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Static Private Member Functions
static std::tuple< float, float, float >	L2ORPosForMatchFunc (const xAOD::L2StandAloneMuon *trig)
static std::tuple< float, float, float >	L2ORPosForMatchFunc (const xAOD::L2CombinedMuon *trig)

Private Attributes
Gaudi::Property< std::vector< float > >	m_monitored_chains_plateau {this, "Plateaus", {}, "List of plateaus of measured trigger"}
std::map< std::string, double >	m_plateaus {}
SG::ReadHandleKey< xAOD::L2CombinedMuonContainer >	m_L2MuonCBIOContainerKey {this, "L2CombinedMuonContainerName_IO", "HLT_MuonL2CBInfoIOmode", "L2MuonCBIO container"}
Gaudi::Property< std::vector< float > >	m_etaBins {this, "EtaBins", {0, 0.9, 1.1, 1.9, 2.1, 9.9}, "eta bins of DR and mass thresholds"}
Gaudi::Property< std::vector< float > >	m_dRCBThres {this, "DRThres", {0.002, 0.001, 0.002, 0.002, 0.002}, "DR threshold of L2CB"}
Gaudi::Property< std::vector< float > >	m_dRbySAThres {this, "dRbySAThres", {0, 0, 0, 0, 0}, "mufast DR threshold of L2CB"}
Gaudi::Property< std::vector< float > >	m_massCBThres {this, "MassThres", {0.004, 0.002, 0.006, 0.006, 0.006}, "mass threshold of L2CB"}
bool	m_RequireDR = true
bool	m_RequireDRbySA = false
bool	m_RequireMass = true
bool	m_RequireSameSign = 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 15 of file L2MuonSAIOMon.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

L2MuonSAIOMon()

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

Definition at line 11 of file L2MuonSAIOMon.cxx.

  : TrigMuonMonitorAlgorithm(name, pSvcLocator)
{}

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

chooseBestMuon()

StatusCode L2MuonSAIOMon::chooseBestMuon ( const std::vector< const xAOD::L2CombinedMuon * > & matchSA_L2IOobjects, std::vector< bool > & passOR, std::vector< unsigned int > & mucombResult ) const

private

Definition at line 858 of file L2MuonSAIOMon.cxx.

                                                                                                                                                                                       {
 
  const double ZERO_LIMIT = 1e-4;
  unsigned int i,j,k;
 
  ATH_MSG_DEBUG( "--- choose best among overlaps & disable EVs (muComb based) ---" );
  for(i=0; i<matchSA_L2IOobjects.size(); i++) {
    ATH_MSG_DEBUG( "++ i=" << i << ": result=" << mucombResult[i] );
    if( mucombResult[i] != i ) {
      ATH_MSG_DEBUG( "   overlap to some one. skip." );
      continue;
    }
    std::vector<unsigned int> others;
    for(j=0; j<matchSA_L2IOobjects.size(); j++) {
      if( mucombResult[j] == mucombResult[i] ) others.emplace_back(j);
    }
    if( others.size() == 1 ) {
      ATH_MSG_DEBUG( "   unique object. keep it active." );
      continue;
    }
    else { // must choose one best
      ATH_MSG_DEBUG( "   overlapped objects among: " << others );
      unsigned int bestMuon = 0;
      float maxPtCombMf  = 0.;
      float mindRRoadRoI = 999.;
      for(k=0; k<others.size(); k++) {
        j=others[k];
 
        float ptCombMf  = std::abs(matchSA_L2IOobjects.at(j)->pt()/1e3);
 
    const float roadPhiP = std::atan2(matchSA_L2IOobjects.at(j)->muSATrack()->dirPhiMS(),1.);
        const float roadPhiM = std::atan2(-1*matchSA_L2IOobjects.at(j)->muSATrack()->dirPhiMS(),-1.);
        const float roadPhi = (std::abs(xAOD::P4Helpers::deltaPhi(roadPhiP, matchSA_L2IOobjects.at(j)->muSATrack()->roiPhi()))
                   < std::abs(xAOD::P4Helpers::deltaPhi(roadPhiM, matchSA_L2IOobjects.at(j)->muSATrack()->roiPhi())))? roadPhiP : roadPhiM;
    float roadAw = 0.;
    if(std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roiEta()) < 1.05) { // barrel
      if( std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(1,0)) > ZERO_LIMIT )      roadAw = matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(1,0);
          else if( std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(2,0)) > ZERO_LIMIT ) roadAw = matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(2,0);
          else if( std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(0,0)) > ZERO_LIMIT ) roadAw = matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(0,0);
        }
    else { // endcap
      if( std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(4,0)) > ZERO_LIMIT )      roadAw = matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(4,0);
          else if( std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(5,0)) > ZERO_LIMIT ) roadAw = matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(5,0);
          else if( std::abs(matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(3,0)) > ZERO_LIMIT ) roadAw = matchSA_L2IOobjects.at(j)->muSATrack()->roadAw(3,0);
        }
    float roadEta = 999.;
        if(std::abs(roadAw) > ZERO_LIMIT) roadEta = -std::log(std::tan(0.5*std::atan(std::abs(roadAw))));
        if(roadAw < 0) roadEta *= -1.;
        float detaRoadRoI = roadEta -  matchSA_L2IOobjects.at(j)->muSATrack()->roiEta();
        float dphiRoadRoI = xAOD::P4Helpers::deltaPhi(roadPhi, matchSA_L2IOobjects.at(j)->muSATrack()->roiPhi());
        float dRRoadRoI = std::sqrt(detaRoadRoI*detaRoadRoI + dphiRoadRoI*dphiRoadRoI);
    ATH_MSG_DEBUG("     j="<< j << " , ptCombMf=" << ptCombMf << ", dRRoadRoI=" << dRRoadRoI);
 
 
    if( (ptCombMf > maxPtCombMf) ||
            (std::abs(ptCombMf - maxPtCombMf) < ZERO_LIMIT &&
             dRRoadRoI < mindRRoadRoI) ) {
          maxPtCombMf  = ptCombMf;
          mindRRoadRoI = dRRoadRoI;
          bestMuon  = j;
        }
      }
      ATH_MSG_DEBUG( "      best is: bestMuon/maxPtCombMf=" << bestMuon << " / " <<  maxPtCombMf );
 
      for(k=0; k<others.size(); k++) {
        j=others[k];
        if( j != bestMuon ) {
          ATH_MSG_DEBUG( "      EventView( j=" << j << " ) is not active" );
 
          passOR.at(j) = false;
        }
        else{
          ATH_MSG_DEBUG( "      EventView( j=" << j << " ) is best one" );
        }
      }
    }
  }
 
 
  return StatusCode::SUCCESS;
}

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

decision_ptthreshold()

StatusCode L2MuonSAIOMon::decision_ptthreshold ( const std::string & chain, std::vector< float > & my_EtaBins, std::vector< float > & my_muCombThres, bool & my_pikCuts, float & my_maxPtToApplyPik, float & my_chi2MaxID ) const

private

Definition at line 1003 of file L2MuonSAIOMon.cxx.

                                                                                                                            {
 
  my_maxPtToApplyPik = 25.;
  my_chi2MaxID = 3.5;
  ATH_MSG_DEBUG("this chain is" << chain);
  if(chain == "HLT_mu4_l2io_L1MU3V"){
    my_EtaBins = {0, 1.05, 1.5, 2.0, 9.9}; //4GeV_v15a
    my_muCombThres = {3.86, 3.77, 3.69, 3.70}; //4GeV_v15a
    my_pikCuts = false;
  }else if(chain == "HLT_mu24_ivarmedium_mu6_l2io_probe_L1MU14FCH"){
    my_EtaBins = {0, 1.05, 1.5, 2.0, 9.9}; //6GeV_v15a
    my_muCombThres = {5.87, 5.79, 5.70, 5.62}; //6GeV_v15a
    my_pikCuts = false;
  }else if(chain == "HLT_mu24_ivarmedium_mu6_l2io_probe_L1MU18VFCH"){
    my_EtaBins = {0, 1.05, 1.5, 2.0, 9.9}; //6GeV_v15a
    my_muCombThres = {5.87, 5.79, 5.70, 5.62}; //6GeV_v15a
    my_pikCuts = false;
  }else{
    ATH_MSG_ERROR("muCombHypo config is NOT defined in this package:chain =   " << chain);
  }
  return StatusCode::SUCCESS;
}

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

defaultPosFunc()

template<class T>

std::tuple< bool, double, double > TrigMuonMonitorAlgorithm::defaultPosFunc ( const T * trig )

inlinestaticprotectedinherited

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(ctx)) {
            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 TrigMuonMonitorAlgorithm::fillHistograms ( const EventContext & ctx ) const

overridevirtualinherited

Function that steers anlayses.

It currently calles four types of analyses, fillVariables, fillVariablesPerOfflineMuon, fillVariablesPerChain and fillVariablesPerOfflineMuonPerChain that can be overridden in subclasses to do specific analyses.

See also

fillVariables, fillVariablesPerOfflineMuon, fillVariablesPerChain and fillVariablesPerOfflineMuonPerChain

Parameters

ctx	EventContext provided by athenaMT

Implements AthMonitorAlgorithm.

Definition at line 22 of file TrigMuonMonitorAlgorithm.cxx.

                                                                                   {
 
 
  bool passToMonitor = selectEvents();
 
  if(passToMonitor){

    ATH_CHECK(fillVariables(ctx));

    for(const std::string& chain : m_monitored_chains){
      ATH_CHECK( fillVariablesPerChain(ctx, chain) );
    }
 

    SG::ReadHandle<xAOD::MuonContainer> muons(m_MuonContainerKey, ctx);
    if (! muons.isValid() ) {
      ATH_MSG_ERROR("evtStore() does not contain muon Collection with name "<< m_MuonContainerKey);
      return StatusCode::FAILURE;
    }
    std::vector<const xAOD::Muon*> probes;
    ATH_CHECK(selectMuons(ctx, muons, probes));
 
    for(const xAOD::Muon* mu : probes){
      ATH_CHECK( fillVariablesPerOfflineMuon(ctx, mu) );

      for(const std::string& chain : m_monitored_chains){
    ATH_CHECK( fillVariablesPerOfflineMuonPerChain(ctx, mu, chain) );
      }
    }
 
  }
 
 
 
  return StatusCode::SUCCESS;
}

fillVariableEtaPhi()

template<class T>

StatusCode TrigMuonMonitorAlgorithm::fillVariableEtaPhi ( const EventContext & ctx, SG::ReadHandleKey< DataVector< T > > ReadHandleKey, std::string && trigstep, std::tuple< bool, double, double >(* PosFunc )(const T *) = &TrigMuonMonitorAlgorithm::defaultPosFunc< T > ) const

protectedinherited

Function that fills variables of etaphi2D plots.

See also

TrigMuonMonitorAlgorithm.icc for the implementation

Parameters

ctx	EventContext provided by athenaMT
ReadHandleKey	SG::ReadHandleKey of online muon.
trigstep	trigger step
PosFunc	Function pointer that implements cuts for the online muon candidates.

fillVariables()

StatusCode TrigMuonMonitorAlgorithm::fillVariables ( const EventContext & ctx ) const

protectedvirtualinherited

Function that fills variables by just retrieving containers of trigger objects.

Users should reimlement in a subclass if needed.

See also

fillHistograms

Parameters

ctx	EventContext provided by athenaMT

Reimplemented in EFMuonMon, L1MuonMon, L2muCombMon, L2MuonSAMon, and MuonTriggerCount.

Definition at line 81 of file TrigMuonMonitorAlgorithm.cxx.

                                                                               {
  return StatusCode::SUCCESS;
}

fillVariablesPerChain()

StatusCode L2MuonSAIOMon::fillVariablesPerChain ( const EventContext & ctx, const std::string & chain ) const

overrideprotectedvirtual

Function that fills variables of trigger objects associated to specified trigger chains.

This is called in the for loop of trigger chains in fillHistograms. Users should reimlement in a subclass if needed.

See also

fillHistograms

Parameters

ctx	EventContext provided by athenaMT
chain	Trigger chain provided in @cfillHistograms

Reimplemented from TrigMuonMonitorAlgorithm.

Definition at line 34 of file L2MuonSAIOMon.cxx.

                                                                                                       {
 
  ATH_MSG_DEBUG ("Filling histograms for " << name() << "...");
 
  if( chain.find("probe") != std::string::npos ) return StatusCode::SUCCESS; // don't use TagAndProbe chains
 
 
  const float ZERO_LIMIT = 0.00001;
 
 
  //TDT workaround
  std::vector< const xAOD::L2CombinedMuon* > Trig_L2IOobjects;
  //std::vector< bool > pass_muCombHypo;
  ATH_CHECK( matchL2IO_wContainer(ctx, chain, Trig_L2IOobjects) );
  //bool pass_muCombHypo_evt = muCombHypo_TDTworkaround(chain, Trig_L2IOobjects, pass_muCombHypo);
  ATH_MSG_DEBUG(" Trig_L2IOobjects.size(): " << Trig_L2IOobjects.size() );
 
  // check basic EDM variables using single muon chain
  for(const auto& Trig_L2IOobject : Trig_L2IOobjects){
    ATH_MSG_DEBUG(" Trig_L2IOobject->muSATrack()->roiWord()/Trig_L2IOobject->pt(): " << Trig_L2IOobject->muSATrack()->roiWord() << "/" << Trig_L2IOobject->pt() );
    // basic EDM variables
    bool mf_failure = false;
    auto cbioPt = Monitored::Scalar<float>(chain+"_Pt",-999.);
    auto cbioEta = Monitored::Scalar<float>(chain+"_Eta",-999.);
    auto cbioPhi = Monitored::Scalar<float>(chain+"_Phi",-999.);
    cbioPt  = Trig_L2IOobject->pt()/1e3 * Trig_L2IOobject->charge(); // convert to GeV
    cbioEta  = Trig_L2IOobject->eta();
    cbioPhi  = Trig_L2IOobject->phi();
    ATH_MSG_DEBUG("cbioPt = " << cbioPt << ", cbioEta =" << cbioEta << ", cbioPhi = " << cbioPhi);
    if(std::abs(cbioPt) < ZERO_LIMIT || std::abs(Trig_L2IOobject->muSATrack()->pt()) < ZERO_LIMIT) mf_failure = true;
 
    if( mf_failure ) continue;
 
    // region variables
    auto isBarrel = Monitored::Scalar<bool>(chain+"_isBarrel",false);
    auto isEndcap = Monitored::Scalar<bool>(chain+"_isEndcap",false);
 
    // define regions
    int saddr = Trig_L2IOobject->muSATrack()->sAddress();
    if(saddr == -1)  isEndcap = true;
    else isBarrel = true;;
 
 
    fill(m_group+"_"+chain, cbioPt, isBarrel, isEndcap);
    fill(m_group+"_"+chain, cbioEta);
    fill(m_group+"_"+chain, cbioPhi, isBarrel, isEndcap);
 
  }
 
  return StatusCode::SUCCESS;
}

fillVariablesPerOfflineMuon()

StatusCode TrigMuonMonitorAlgorithm::fillVariablesPerOfflineMuon ( const EventContext & ctx, const xAOD::Muon * mu ) const

protectedvirtualinherited

Function that fills variables that are compared to offline muons but the trigger chains are not specified.

This is called in the for loop of offline muons in fillHistograms. Users should reimlement in a subclass if needed.

See also

fillHistograms

Parameters

ctx	EventContext provided by athenaMT
mu	Pointer to an offline muon provided in fillHistograms

Reimplemented in EFMuonMon, L2muCombMon, and L2MuonSAMon.

Definition at line 86 of file TrigMuonMonitorAlgorithm.cxx.

                                                                                                              {
  return StatusCode::SUCCESS;
}

fillVariablesPerOfflineMuonPerChain()

StatusCode L2MuonSAIOMon::fillVariablesPerOfflineMuonPerChain ( const EventContext & ctx, const xAOD::Muon * mu, const std::string & chain ) const

overrideprotectedvirtual

Function that fills variables of trigger objects associated to specified trigger chains comparing offline muons.

This is called in the for loop of trigger chains and offline muons in fillHistograms. Users should reimlement in a subclass if needed.

See also

fillHistograms

Parameters

ctx	EventContext provided by athenaMT
mu	Pointer to an offline muon provided in fillHistograms
chain	Trigger chain provided in fillHistograms

Reimplemented from TrigMuonMonitorAlgorithm.

Definition at line 87 of file L2MuonSAIOMon.cxx.

                                                                                                                                         {
 
  ATH_MSG_DEBUG ("Filling histograms for " << name() << "...");
 
  const float ZERO_LIMIT = 0.00001;
 
  // get the best L2Inside-Out object matched to offline muon(dR between L2Inside-Out object and offline muon is minimum)
  const xAOD::L2CombinedMuon* Trig_L2IOobject = searchL2InsideOut(ctx, mu, chain);
 
  // offline muon variables
  auto offPt = Monitored::Scalar<float>(chain+"_offPt",-999.);
  auto offEta = Monitored::Scalar<float>(chain+"_offEta",-999.);
  auto offPhi = Monitored::Scalar<float>(chain+"_offPhi",-999.);
  offPt = mu->pt()/1e3 * mu->charge(); // convert to GeV
  offEta = mu->eta();
  offPhi = mu->phi();
  float offCharge = mu->charge();
 
 
  if( chain.find("probe") != std::string::npos ){ // L2Inside-Out efficiency using Tag&Probe chain
    if(chain.find("L1MU14FCH") != std::string::npos){
      if ( !getTrigDecisionTool()->isPassed("HLT_mu24_ivarmedium_L1MU14FCH", TrigDefs::requireDecision) ) return StatusCode::SUCCESS; // impose trigger pass in order to eliminate bias
    }
    else if(chain.find("L1MU18VFCH") != std::string::npos){
      if ( !getTrigDecisionTool()->isPassed("HLT_mu24_ivarmedium_L1MU18VFCH", TrigDefs::requireDecision) ) return StatusCode::SUCCESS; // impose trigger pass in order to eliminate bias
    }
    else
    {
      return StatusCode::SUCCESS;
    }
 
    // search tag offline muon
    const xAOD::Muon* tag = searchTagOfflineMuon( ctx, mu );
    if( tag == nullptr ) return StatusCode::SUCCESS; // mu doesn't have no tag muons  
 
 
    // efficiency variables
    auto passL2InsideOut = Monitored::Scalar<bool>(chain+"_passL2InsideOut",false);
    auto passL2SA = Monitored::Scalar<bool>(chain+"_passL2SA",false);
    auto offdR = Monitored::Scalar<float>(chain+"_offdR",1000.);
    const auto* tag_ms_track = tag->trackParticle(xAOD::Muon::TrackParticleType::MuonSpectrometerTrackParticle);
    if( !tag_ms_track) return StatusCode::SUCCESS; // tag muon dosen't have ms track
    const auto* probe_ms_track = mu->trackParticle(xAOD::Muon::TrackParticleType::MuonSpectrometerTrackParticle);
    if( !probe_ms_track ) return StatusCode::SUCCESS; // probe muon dosen't have ms track
    offdR = xAOD::P4Helpers::deltaR(tag_ms_track, probe_ms_track);
 
    passL2InsideOut = false;
    passL2SA = false;
 
    // retrieve probe l2SA objects
    int legIndex_probe = 1; // probe
    std::vector< TrigCompositeUtils::LinkInfo<xAOD::L2StandAloneMuonContainer> > featureCont = getTrigDecisionTool()->features<xAOD::L2StandAloneMuonContainer>( chain,
                                                                               TrigDefs::includeFailedDecisions,
                                                                               "HLT_MuonL2SAInfo",
                                                                               TrigDefs::lastFeatureOfType,
                                                                               TrigCompositeUtils::featureString(),
                                                                               legIndex_probe );
 
    for(const TrigCompositeUtils::LinkInfo<xAOD::L2StandAloneMuonContainer>& probe_L2SALinkInfo : featureCont){
        ATH_CHECK( probe_L2SALinkInfo.isValid() );
        const ElementLink<xAOD::L2StandAloneMuonContainer> probe_L2SAobject = probe_L2SALinkInfo.link;
        if( m_matchTool->isMatchedL2SA( (*probe_L2SAobject), mu ) ){
          if( probe_L2SALinkInfo.state == TrigCompositeUtils::ActiveState::ACTIVE ) passL2SA = true;
        }
    }
 
    if(passL2SA == true){
      if(m_matchTool->isMatchedL2InsideOut( Trig_L2IOobject, mu ))
      {
        bool isPass = false;
        ATH_CHECK(isPassedmuCombHypo( chain, Trig_L2IOobject ,isPass));
        passL2InsideOut = isPass;
      }
    }
    else{
      return StatusCode::SUCCESS;
    }
   
 
    fill(m_group+"_"+chain, passL2InsideOut, passL2SA, offPt);
 
    if( mu->pt()/1e3 > m_plateaus.at(chain) ){
      fill(m_group+"_"+chain, passL2InsideOut, passL2SA, offEta, offPhi, offdR);
    }
  }
  else{ // make detail histograms using normal L2Inside-Out chain
    if( Trig_L2IOobject == nullptr ) return StatusCode::SUCCESS; // no L2Inside-Out objects
    if( ! m_matchTool->isMatchedL2InsideOut(Trig_L2IOobject, mu) ) return StatusCode::SUCCESS; // offline muons is not matched to any L2Inside-Out objects
 
 
    //L2Muon chamberID index
    enum chamberID {
      Inn_Barrel = 0,
      Mid_Barrel,
      Out_Barrel,
      Inn_Endcap,
      Mid_Endcap,
      Out_Endcap
    };
    std::vector< int > L2Muon_chamberID_index;
    auto mon_L2Muon_chamberID_index = Monitored::Collection(chain+"_L2Muon_chamberID_index",L2Muon_chamberID_index);
    for( int i = 0; i < 6; i++){
      L2Muon_chamberID_index.push_back(i);
    }
 
    // dR wrt offline
    auto dRmin = Monitored::Scalar<float>(chain+"_dRmin",1000.);
    dRmin = xAOD::P4Helpers::deltaR(mu, Trig_L2IOobject, false);
 
    // region variables
    const float ETA_OF_BARREL   = 1.05;
    auto regionBE = Monitored::Scalar<int>(chain+"_regionBE",-999);
    auto isBarrel = Monitored::Scalar<bool>(chain+"_isBarrel",false);
    auto isBarrelA = Monitored::Scalar<bool>(chain+"_isBarrelA",false);
    auto isBarrelC = Monitored::Scalar<bool>(chain+"_isBarrelC",false);
    auto isEndcap = Monitored::Scalar<bool>(chain+"_isEndcap",false);
    auto isEndcapA = Monitored::Scalar<bool>(chain+"_isEndcapA",false);
    auto isEndcapC = Monitored::Scalar<bool>(chain+"_isEndcapC",false);
 
    // offline pt variables
    auto pt4to6 = Monitored::Scalar<bool>(chain+"_pt4to6",false);
    auto pt6to8 = Monitored::Scalar<bool>(chain+"_pt6to8",false);
    auto ptover8 = Monitored::Scalar<bool>(chain+"_ptover8",false);
 
    // define region
    if( std::abs(offEta) < ETA_OF_BARREL ) {
      regionBE = 0;
      isBarrel = true;
      if( offEta > 0. ) isBarrelA = true;
      else isBarrelC = true;
    }
    else{
      regionBE = 1;
      isEndcap = true;
      if( offEta > 0. ) isEndcapA = true;
      else isEndcapC = true;
    }
 
    if( std::abs(offPt) > 4 ){
      if( std::abs(offPt) < 6 ) pt4to6 = true;
      else if( std::abs(offPt) < 8 ) pt6to8 = true;
      else ptover8 = true;
    }
 
    // basic variables
    auto cbioPt = Monitored::Scalar<float>(chain+"_Pt_wrt_offline",-999.);
    auto cbioEta = Monitored::Scalar<float>(chain+"_Eta_wrt_offline",-999.);
    auto cbioPhi = Monitored::Scalar<float>(chain+"_Phi_wrt_offline",-999.);
    cbioPt  = Trig_L2IOobject->pt()/1e3 * Trig_L2IOobject->charge(); // convert to GeV
    cbioEta  = Trig_L2IOobject->eta();
    cbioPhi  = Trig_L2IOobject->phi();
 
    // L2Inside-Out track multiplicity per L2SA track
    auto L2InsideOut_multiplicity =  Monitored::Scalar<float>(chain+"_L2InsideOut_track_multiplicity",-999.);
    L2InsideOut_multiplicity = 0;
    std::vector< const xAOD::L2CombinedMuon* > Trig_L2IOobjects_tmp;
    ATH_CHECK( matchL2IO_wContainer(ctx, chain, Trig_L2IOobjects_tmp) );
    for(const auto& Trig_L2IOobject_tmp : Trig_L2IOobjects_tmp){
      if( Trig_L2IOobject_tmp->muSATrack()->roiWord() == Trig_L2IOobject->muSATrack()->roiWord() ) L2InsideOut_multiplicity++;
    }
 
    // pt resolution
    auto ptresol = Monitored::Scalar<float>(chain+"_ptresol",-999.);
    if ( std::abs(offPt) > ZERO_LIMIT && std::abs(cbioPt) > ZERO_LIMIT ){
      ptresol = std::abs(cbioPt)/std::abs(offPt) - 1.;
    }
 
    std::vector< float > ptresol_pos, ptresol_neg;
    auto mon_ptresol_pos = Monitored::Collection(chain+"_ptresol_pos",ptresol_pos);
    auto mon_ptresol_neg = Monitored::Collection(chain+"_ptresol_neg",ptresol_neg);
    if( offCharge > 0. ) ptresol_pos.push_back(ptresol);
    else ptresol_neg.push_back(ptresol);
 
    // distance bw FTFroad and offlinesegment
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_Inn_Barrel;
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_Mid_Barrel;
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_Out_Barrel;
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_Inn_Endcap;
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_Mid_Endcap;
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_Out_Endcap;
    auto mon_distance_bw_FTFroad_and_offlinesegment_Inn_Barrel = Monitored::Collection(chain+"_distance_bw_FTFroad_and_offlinesegment_Inn_Barrel",distance_bw_FTFroad_and_offlinesegment_Inn_Barrel);
    auto mon_distance_bw_FTFroad_and_offlinesegment_Mid_Barrel = Monitored::Collection(chain+"_distance_bw_FTFroad_and_offlinesegment_Mid_Barrel",distance_bw_FTFroad_and_offlinesegment_Mid_Barrel);
    auto mon_distance_bw_FTFroad_and_offlinesegment_Out_Barrel = Monitored::Collection(chain+"_distance_bw_FTFroad_and_offlinesegment_Out_Barrel",distance_bw_FTFroad_and_offlinesegment_Out_Barrel);
    auto mon_distance_bw_FTFroad_and_offlinesegment_Inn_Endcap = Monitored::Collection(chain+"_distance_bw_FTFroad_and_offlinesegment_Inn_Endcap",distance_bw_FTFroad_and_offlinesegment_Inn_Endcap);
    auto mon_distance_bw_FTFroad_and_offlinesegment_Mid_Endcap = Monitored::Collection(chain+"_distance_bw_FTFroad_and_offlinesegment_Mid_Endcap",distance_bw_FTFroad_and_offlinesegment_Mid_Endcap);
    auto mon_distance_bw_FTFroad_and_offlinesegment_Out_Endcap = Monitored::Collection(chain+"_distance_bw_FTFroad_and_offlinesegment_Out_Endcap",distance_bw_FTFroad_and_offlinesegment_Out_Endcap);
 
    std::vector< float > distance_bw_FTFroad_and_offlinesegment_vec;
    std::vector< float > FTFroad_Aw;
    std::vector< float > FTFroad_Bw;
    std::vector< bool > FTFroad_fill;
    for(int i=0; i<6; i++){
      distance_bw_FTFroad_and_offlinesegment_vec.push_back(10000.);
      FTFroad_Aw.push_back(Trig_L2IOobject->muSATrack()->roadAw(i, 0));
      FTFroad_Bw.push_back(Trig_L2IOobject->muSATrack()->roadBw(i, 0));
      FTFroad_fill.push_back(false);
    }
 
    // MDT hits residual
    std::vector<float> res_Inn_Barrel, res_Mid_Barrel, res_Out_Barrel, res_Inn_Endcap, res_Mid_Endcap, res_Out_Endcap;
    auto mon_res_Inn_Barrel = Monitored::Collection(chain+"_MDT_residual_Inn_Barrel",res_Inn_Barrel);
    auto mon_res_Mid_Barrel = Monitored::Collection(chain+"_MDT_residual_Mid_Barrel",res_Mid_Barrel);
    auto mon_res_Out_Barrel = Monitored::Collection(chain+"_MDT_residual_Out_Barrel",res_Out_Barrel);
    auto mon_res_Inn_Endcap = Monitored::Collection(chain+"_MDT_residual_Inn_Endcap",res_Inn_Endcap);
    auto mon_res_Mid_Endcap = Monitored::Collection(chain+"_MDT_residual_Mid_Endcap",res_Mid_Endcap);
    auto mon_res_Out_Endcap = Monitored::Collection(chain+"_MDT_residual_Out_Endcap",res_Out_Endcap);
 
    // distance bw MDT hits and offlinesegment
    std::vector< float > distance_bw_MDT_and_offlinesegment_Inn_Barrel;
    std::vector< float > distance_bw_MDT_and_offlinesegment_Mid_Barrel;
    std::vector< float > distance_bw_MDT_and_offlinesegment_Out_Barrel;
    std::vector< float > distance_bw_MDT_and_offlinesegment_Inn_Endcap;
    std::vector< float > distance_bw_MDT_and_offlinesegment_Mid_Endcap;
    std::vector< float > distance_bw_MDT_and_offlinesegment_Out_Endcap;
    auto mon_distance_bw_MDT_and_offlinesegment_Inn_Barrel = Monitored::Collection(chain+"_distance_bw_MDT_and_offlinesegment_Inn_Barrel",distance_bw_MDT_and_offlinesegment_Inn_Barrel);
    auto mon_distance_bw_MDT_and_offlinesegment_Mid_Barrel = Monitored::Collection(chain+"_distance_bw_MDT_and_offlinesegment_Mid_Barrel",distance_bw_MDT_and_offlinesegment_Mid_Barrel);
    auto mon_distance_bw_MDT_and_offlinesegment_Out_Barrel = Monitored::Collection(chain+"_distance_bw_MDT_and_offlinesegment_Out_Barrel",distance_bw_MDT_and_offlinesegment_Out_Barrel);
    auto mon_distance_bw_MDT_and_offlinesegment_Inn_Endcap = Monitored::Collection(chain+"_distance_bw_MDT_and_offlinesegment_Inn_Endcap",distance_bw_MDT_and_offlinesegment_Inn_Endcap);
    auto mon_distance_bw_MDT_and_offlinesegment_Mid_Endcap = Monitored::Collection(chain+"_distance_bw_MDT_and_offlinesegment_Mid_Endcap",distance_bw_MDT_and_offlinesegment_Mid_Endcap);
    auto mon_distance_bw_MDT_and_offlinesegment_Out_Endcap = Monitored::Collection(chain+"_distance_bw_MDT_and_offlinesegment_Out_Endcap",distance_bw_MDT_and_offlinesegment_Out_Endcap);
    std::vector< float > distance_bw_MDT_and_offlinesegment_vec;
    std::vector< int > MDTHitChamber_fill;
    std::vector< int > MDTHitChamber;
    std::vector< float > MDTHitR;
    std::vector< float > MDTHitZ;
 
    // # of MDT hits
    std::vector<int> MDT_N_Inn_Barrel, MDT_N_Mid_Barrel, MDT_N_Out_Barrel, MDT_N_Inn_Endcap, MDT_N_Mid_Endcap, MDT_N_Out_Endcap;
    auto mon_MDT_N_Inn_Barrel = Monitored::Collection(chain+"_MDT_N_Inn_Barrel",MDT_N_Inn_Barrel);
    auto mon_MDT_N_Mid_Barrel = Monitored::Collection(chain+"_MDT_N_Mid_Barrel",MDT_N_Mid_Barrel);
    auto mon_MDT_N_Out_Barrel = Monitored::Collection(chain+"_MDT_N_Out_Barrel",MDT_N_Out_Barrel);
    auto mon_MDT_N_Inn_Endcap = Monitored::Collection(chain+"_MDT_N_Inn_Endcap",MDT_N_Inn_Endcap);
    auto mon_MDT_N_Mid_Endcap = Monitored::Collection(chain+"_MDT_N_Mid_Endcap",MDT_N_Mid_Endcap);
    auto mon_MDT_N_Out_Endcap = Monitored::Collection(chain+"_MDT_N_Out_Endcap",MDT_N_Out_Endcap);
    int n_mdthits_BI = 0;
    int n_mdthits_BM = 0;
    int n_mdthits_BO = 0;
    int n_mdthits_EI = 0;
    int n_mdthits_EM = 0;
    int n_mdthits_EO = 0;
 
    int n_mdt_hits = Trig_L2IOobject->muSATrack()->nMdtHits();
    for(int i_tube=0; i_tube<n_mdt_hits; i_tube++){
      if( Trig_L2IOobject->muSATrack()->mdtHitIsOutlier(i_tube) != 0 ) continue;
      float res = Trig_L2IOobject->muSATrack()->mdtHitResidual(i_tube);
      int imr = Trig_L2IOobject->muSATrack()->mdtHitChamber(i_tube);
      MDTHitChamber.push_back(imr);
      MDTHitR.push_back(Trig_L2IOobject->muSATrack()->mdtHitR(i_tube));
      MDTHitZ.push_back(Trig_L2IOobject->muSATrack()->mdtHitZ(i_tube));
 
      if( imr == Inn_Barrel ){
    n_mdthits_BI++;
    res_Inn_Barrel.push_back(res);
      }
      else if( imr == Mid_Barrel ){
    n_mdthits_BM++;
    res_Mid_Barrel.push_back(res);
      }
      else if( imr == Out_Barrel ){
    n_mdthits_BO++;
    res_Out_Barrel.push_back(res);
      }
      else if( imr == Inn_Endcap ){
    n_mdthits_EI++;
    res_Inn_Endcap.push_back(res);
      }
      else if( imr == Mid_Endcap ){
    n_mdthits_EM++;
    res_Mid_Endcap.push_back(res);
      }
      else if( imr == Out_Endcap ){
    n_mdthits_EO++;
    res_Out_Endcap.push_back(res);
      }
    }
 
    // reconstruction efficiency of superpoint
    std::vector<bool> superpoint_exist_pt4to6, superpoint_exist_pt6to8, superpoint_exist_ptover8;
    auto mon_superpoint_exist_pt4to6 = Monitored::Collection(chain+"_superpoint_pt4to6",superpoint_exist_pt4to6);
    auto mon_superpoint_exist_pt6to8 = Monitored::Collection(chain+"_superpoint_pt6to8",superpoint_exist_pt6to8);
    auto mon_superpoint_exist_ptover8 = Monitored::Collection(chain+"_superpoint_ptover8",superpoint_exist_ptover8);
    std::vector<bool> segment_superpoint_exist(6, false);
    std::vector<bool> offlinesegment_exist_pt4to6, offlinesegment_exist_pt6to8, offlinesegment_exist_ptover8;
    auto mon_offlinesegment_exist_pt4to6 = Monitored::Collection(chain+"_offlinesegment_pt4to6",offlinesegment_exist_pt4to6);
    auto mon_offlinesegment_exist_pt6to8 = Monitored::Collection(chain+"_offlinesegment_pt6to8",offlinesegment_exist_pt6to8);
    auto mon_offlinesegment_exist_ptover8 = Monitored::Collection(chain+"_offlinesegment_ptover8",offlinesegment_exist_ptover8);
 
    // # of superpoint
    auto superpoint_multiplicity = Monitored::Scalar<int>(chain+"_superpoint_multiplicity",0);
    std::vector< bool > superpoint_exist;
    std::vector< float > superpointR;
    int Num_L2Muon_chamberID = 12;
    for( int i_chamber = 0; i_chamber < Num_L2Muon_chamberID; i_chamber++){
      if( Trig_L2IOobject->muSATrack()->superPointR(i_chamber) < ZERO_LIMIT ){
    superpoint_exist.push_back(false);
      }
      else{
    superpoint_exist.push_back(true);
    superpoint_multiplicity++;
      }
      superpointR.push_back(Trig_L2IOobject->muSATrack()->superPointR(i_chamber));
    }
 
 
    std::vector< bool > segment_exist(6, false);
    for(unsigned int i_seg = 0; i_seg < mu->nMuonSegments(); i_seg++){
      const xAOD::MuonSegment* segment = mu->muonSegment(i_seg);
      if(!segment) continue;
      float segmentX = segment->x();
      float segmentY = segment->y();
      float segmentZ = segment->z();
      float segmentR = std::sqrt( std::pow(segmentX, 2.0) + std::pow(segmentY, 2.0) );
      float segmentPx = segment->px();
      float segmentPy = segment->py();
      float segmentPz = segment->pz();
      float segmentSector = segment->sector();
      using namespace Muon::MuonStationIndex;
      int segmentChamberIndex = toInt(segment->chamberIndex());
      float distance_bw_FTFroad_and_offlinesegment = 99999.;
      float distance_bw_MDT_and_offlinesegment = 99999.;
      int roadChamberIndex = -1;
      int MDTChamberIndex = -1;
      if( segmentChamberIndex == 0 || segmentChamberIndex == 1 ){ // Inner Barrel
    segment_exist.at(Inn_Barrel) = true;
    roadChamberIndex = Inn_Barrel;
    MDTChamberIndex = Inn_Barrel;
      }
      else if( segmentChamberIndex == 2 || segmentChamberIndex == 3 ){ // Middle Barrel
    segment_exist.at(Mid_Barrel) = true;
    roadChamberIndex = Mid_Barrel;
    MDTChamberIndex = Mid_Barrel;
      }
      else if( segmentChamberIndex == 4 || segmentChamberIndex == 5 ){ // Outer Barrel
    segment_exist.at(Out_Barrel) = true;
    roadChamberIndex = Out_Barrel;
    MDTChamberIndex = Out_Barrel;
      }
      else if( segmentChamberIndex == 7 || segmentChamberIndex == 8 ){ // Inner Endcap
    segment_exist.at(Inn_Endcap) = true;
    roadChamberIndex = Inn_Endcap;
    MDTChamberIndex = Inn_Endcap;
      }
      else if( segmentChamberIndex == 9 || segmentChamberIndex == 10 ){ // Middle Endcap
    segment_exist.at(Mid_Endcap) = true;
    roadChamberIndex = Mid_Endcap;
    MDTChamberIndex = Mid_Endcap;
      }
      else if( segmentChamberIndex == 11 || segmentChamberIndex == 12 ){ // Outer Endcap
    segment_exist.at(Out_Endcap) = true;
    roadChamberIndex = Out_Endcap;
    MDTChamberIndex = Out_Endcap;
      }
 
      // Calc distance bw FTFroad and offlinesegment
      if( roadChamberIndex != -1 ){
    if( FTFroad_Aw.at(roadChamberIndex) > ZERO_LIMIT || FTFroad_Bw.at(roadChamberIndex) > ZERO_LIMIT ){
      FTFroad_fill.at(roadChamberIndex) = true;
      if( FTFroad_Aw.at(roadChamberIndex) < ZERO_LIMIT) distance_bw_FTFroad_and_offlinesegment = segmentR - FTFroad_Bw.at(roadChamberIndex);
      else{
        float ia = 1.0/FTFroad_Aw.at(roadChamberIndex);
        float iaq = ia * ia;
        distance_bw_FTFroad_and_offlinesegment = (segmentZ - ia * (segmentR - FTFroad_Bw.at(roadChamberIndex)))/std::sqrt(1.0 + iaq);
      }
      if( std::abs(distance_bw_FTFroad_and_offlinesegment) < std::abs(distance_bw_FTFroad_and_offlinesegment_vec.at(roadChamberIndex)) )
        distance_bw_FTFroad_and_offlinesegment_vec.at(roadChamberIndex) = distance_bw_FTFroad_and_offlinesegment;
    }
      }
 
      // Calc distance bw MDT hits and offlinesegment
      if( MDTChamberIndex != -1 ){
    float sector_phi = M_PI*(segmentSector - 1.0)/8.0;
    float segmentR_projection = segmentX * std::cos(sector_phi) + segmentY * std::sin(sector_phi);
    float segmentPr_projection = segmentPx * std::cos(sector_phi) + segmentPy * std::sin(sector_phi);
    
    for(unsigned int i_tube=0; i_tube<MDTHitChamber.size(); i_tube++){
      if( MDTHitChamber.at(i_tube) != MDTChamberIndex ) continue;
      if( MDTChamberIndex < 3 ){ //Barrel
        if( std::abs(segmentPz) < ZERO_LIMIT ) distance_bw_MDT_and_offlinesegment = MDTHitZ.at(i_tube) - segmentZ;
        else{
          float denominator = segmentPr_projection/segmentPz;
          if( std::abs(denominator) < ZERO_LIMIT ) continue;
          distance_bw_MDT_and_offlinesegment = MDTHitZ.at(i_tube) - ((MDTHitR.at(i_tube) - segmentR_projection)/denominator  + segmentZ);
        }
      }
      else{ //Endcap
        if( std::abs(segmentPz) < ZERO_LIMIT ){ 
          distance_bw_MDT_and_offlinesegment = MDTHitR.at(i_tube) - segmentR;
        }
        else{
          float coeffi = (MDTHitZ.at(i_tube) - segmentZ)/segmentPz;
          float segmentR_extrapolated = std::sqrt(std::pow(segmentX + coeffi * segmentPx, 2.0) + std::pow(segmentY + coeffi * segmentPy, 2.0));
          distance_bw_MDT_and_offlinesegment = MDTHitR.at(i_tube) - segmentR_extrapolated;
        }
      }
      distance_bw_MDT_and_offlinesegment_vec.push_back(distance_bw_MDT_and_offlinesegment);
      MDTHitChamber_fill.push_back(MDTHitChamber.at(i_tube));
    }
      }
    }
 
    if( FTFroad_fill.at(Inn_Barrel) ){
      distance_bw_FTFroad_and_offlinesegment_Inn_Barrel.push_back(distance_bw_FTFroad_and_offlinesegment_vec.at(Inn_Barrel));
    }
    if( FTFroad_fill.at(Mid_Barrel) ){
      distance_bw_FTFroad_and_offlinesegment_Mid_Barrel.push_back(distance_bw_FTFroad_and_offlinesegment_vec.at(Mid_Barrel));
    }
    if( FTFroad_fill.at(Out_Barrel) ){
      distance_bw_FTFroad_and_offlinesegment_Out_Barrel.push_back(distance_bw_FTFroad_and_offlinesegment_vec.at(Out_Barrel));
    }
    if( FTFroad_fill.at(Inn_Endcap) ){
      distance_bw_FTFroad_and_offlinesegment_Inn_Endcap.push_back(distance_bw_FTFroad_and_offlinesegment_vec.at(Inn_Endcap));
    }
    if( FTFroad_fill.at(Mid_Endcap) ){
      distance_bw_FTFroad_and_offlinesegment_Mid_Endcap.push_back(distance_bw_FTFroad_and_offlinesegment_vec.at(Mid_Endcap));
    }
    if( FTFroad_fill.at(Out_Endcap) ){
      distance_bw_FTFroad_and_offlinesegment_Out_Endcap.push_back(distance_bw_FTFroad_and_offlinesegment_vec.at(Out_Endcap));
    }
 
    for( unsigned int i = 0; i < distance_bw_MDT_and_offlinesegment_vec.size(); i++ ){
      if( MDTHitChamber_fill.at(i) == Inn_Barrel ){
    distance_bw_MDT_and_offlinesegment_Inn_Barrel.push_back(distance_bw_MDT_and_offlinesegment_vec.at(i));
      }
      else if( MDTHitChamber_fill.at(i) == Mid_Barrel ){
    distance_bw_MDT_and_offlinesegment_Mid_Barrel.push_back(distance_bw_MDT_and_offlinesegment_vec.at(i));
      }
      else if( MDTHitChamber_fill.at(i) == Out_Barrel ){
    distance_bw_MDT_and_offlinesegment_Out_Barrel.push_back(distance_bw_MDT_and_offlinesegment_vec.at(i));
      }
      else if( MDTHitChamber_fill.at(i) == Inn_Endcap ){
    distance_bw_MDT_and_offlinesegment_Inn_Endcap.push_back(distance_bw_MDT_and_offlinesegment_vec.at(i));
      }
      else if( MDTHitChamber_fill.at(i) == Mid_Endcap ){
    distance_bw_MDT_and_offlinesegment_Mid_Endcap.push_back(distance_bw_MDT_and_offlinesegment_vec.at(i));
      }
      else if( MDTHitChamber_fill.at(i) == Out_Endcap ){
    distance_bw_MDT_and_offlinesegment_Out_Endcap.push_back(distance_bw_MDT_and_offlinesegment_vec.at(i));
      }
      else{
    ATH_MSG_WARNING( "undefined chamberID is pushed back into MDTHitChamber_fill" );
      }
    }
 
    if( segment_exist.at(Inn_Barrel) ){
      MDT_N_Inn_Barrel.push_back(n_mdthits_BI);
      segment_superpoint_exist.at(Inn_Barrel) = superpoint_exist.at(Inn_Barrel);
    }
    if( segment_exist.at(Mid_Barrel) ){
      MDT_N_Mid_Barrel.push_back(n_mdthits_BM);
      segment_superpoint_exist.at(Mid_Barrel) = superpoint_exist.at(Mid_Barrel);
    }
    if( segment_exist.at(Out_Barrel) ){
      MDT_N_Out_Barrel.push_back(n_mdthits_BO);
      segment_superpoint_exist.at(Out_Barrel) = superpoint_exist.at(Out_Barrel);
    }
    if( segment_exist.at(Inn_Endcap) ){
      MDT_N_Inn_Endcap.push_back(n_mdthits_EI);
      segment_superpoint_exist.at(Inn_Endcap) = superpoint_exist.at(Inn_Endcap);
    }
    if( segment_exist.at(Mid_Endcap) ){
      MDT_N_Mid_Endcap.push_back(n_mdthits_EM);
      segment_superpoint_exist.at(Mid_Endcap) = superpoint_exist.at(Mid_Endcap);
    }
    if( segment_exist.at(Out_Endcap) ){
      MDT_N_Out_Endcap.push_back(n_mdthits_EO);
      segment_superpoint_exist.at(Out_Endcap) = superpoint_exist.at(Out_Endcap);
    }
 
    if( pt4to6 ){
      offlinesegment_exist_pt4to6 = segment_exist;
      superpoint_exist_pt4to6 = segment_superpoint_exist;
    }
    else if( pt6to8 ){
      offlinesegment_exist_pt6to8 = segment_exist;
      superpoint_exist_pt6to8 = segment_superpoint_exist;
    }
    else if( ptover8 ){
      offlinesegment_exist_ptover8 = segment_exist;
      superpoint_exist_ptover8 = segment_superpoint_exist;
    }
 
    fill(m_group+"_"+chain, dRmin, isBarrel, isEndcap);
    fill(m_group+"_"+chain, cbioPt, isBarrel, isEndcap);
    fill(m_group+"_"+chain, cbioEta);
    fill(m_group+"_"+chain, cbioPhi, isBarrel, isEndcap);
    fill(m_group+"_"+chain, L2InsideOut_multiplicity, offPt, isBarrel, isEndcap);
    fill(m_group+"_"+chain, ptresol, offEta, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, ptresol, offPt, isBarrelA, isBarrelC, isEndcapA, isEndcapC);
    fill(m_group+"_"+chain, mon_ptresol_pos, mon_ptresol_neg, isBarrelA, isBarrelC, isEndcapA, isEndcapC);
    fill(m_group+"_"+chain, mon_distance_bw_FTFroad_and_offlinesegment_Inn_Barrel, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_distance_bw_FTFroad_and_offlinesegment_Mid_Barrel, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_distance_bw_FTFroad_and_offlinesegment_Out_Barrel, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_distance_bw_FTFroad_and_offlinesegment_Inn_Endcap, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_distance_bw_FTFroad_and_offlinesegment_Mid_Endcap, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_distance_bw_FTFroad_and_offlinesegment_Out_Endcap, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_res_Inn_Barrel);
    fill(m_group+"_"+chain, mon_res_Mid_Barrel);
    fill(m_group+"_"+chain, mon_res_Out_Barrel);
    fill(m_group+"_"+chain, mon_res_Inn_Endcap);
    fill(m_group+"_"+chain, mon_res_Mid_Endcap);
    fill(m_group+"_"+chain, mon_res_Out_Endcap);
    fill(m_group+"_"+chain, mon_distance_bw_MDT_and_offlinesegment_Inn_Barrel);
    fill(m_group+"_"+chain, mon_distance_bw_MDT_and_offlinesegment_Mid_Barrel);
    fill(m_group+"_"+chain, mon_distance_bw_MDT_and_offlinesegment_Out_Barrel);
    fill(m_group+"_"+chain, mon_distance_bw_MDT_and_offlinesegment_Inn_Endcap);
    fill(m_group+"_"+chain, mon_distance_bw_MDT_and_offlinesegment_Mid_Endcap);
    fill(m_group+"_"+chain, mon_distance_bw_MDT_and_offlinesegment_Out_Endcap);
    fill(m_group+"_"+chain, mon_MDT_N_Inn_Barrel);
    fill(m_group+"_"+chain, mon_MDT_N_Mid_Barrel);
    fill(m_group+"_"+chain, mon_MDT_N_Out_Barrel);
    fill(m_group+"_"+chain, mon_MDT_N_Inn_Endcap);
    fill(m_group+"_"+chain, mon_MDT_N_Mid_Endcap);
    fill(m_group+"_"+chain, mon_MDT_N_Out_Endcap);
    fill(m_group+"_"+chain, superpoint_multiplicity, regionBE, pt4to6, pt6to8, ptover8);
    fill(m_group+"_"+chain, mon_L2Muon_chamberID_index, mon_superpoint_exist_pt4to6, mon_offlinesegment_exist_pt4to6);
    fill(m_group+"_"+chain, mon_L2Muon_chamberID_index, mon_superpoint_exist_pt6to8, mon_offlinesegment_exist_pt6to8);
    fill(m_group+"_"+chain, mon_L2Muon_chamberID_index, mon_superpoint_exist_ptover8, mon_offlinesegment_exist_ptover8);
  }
 
  return StatusCode::SUCCESS;
}

fillVariablesRatioPlots()

template<class T, class FUNCT>

StatusCode TrigMuonMonitorAlgorithm::fillVariablesRatioPlots ( const EventContext & ctx, const xAOD::Muon * mu, std::string && trigstep, xAOD::Muon::TrackParticleType type, FUNCT matchFunc ) const

protectedinherited

Function that fills variables of ratio plots.

See also

TrigMuonMonitorAlgorithm.icc for the implementation

Parameters

ctx	EventContext provided by athenaMT
mu	Pointer to an offline muon provided in fillHistograms
trigstep	trigger step
type	xAOD::Muon::TrackParticleType of offline muon
matchFunc	Function pointer that implements cuts for the online muon candidates gotten by ReadHandle.

filterPassed()

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

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

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);
    //Check if the tool exists in the map
    if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
        return m_tools[idx->second];
    }
    else {
      //Check if the map is empty
     if (m_toolLookupMap.empty()) {
            ATH_MSG_FATAL("The m_toolLookupMap is empty. The tool " << name << " cannot be found in an empty map.");
            return m_dummy;
      } 
      //If the map is not empty and the tool was not found, print a fatal error
      if (!m_toolLookupMap.empty()) {
            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;
      } 
      if (!isInitialized()) {
            ATH_MSG_FATAL(
            "It seems that the AthMonitorAlgorithm::initialize was not called "
            "in derived class initialize method, group name: " << name);
      } 
    }
    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 203 of file AthMonitorAlgorithm.cxx.

                                                                                       {
    return m_trigDecTool;
}

initialize()

StatusCode L2MuonSAIOMon::initialize ( )

overridevirtual

initialize

Returns

StatusCode

Reimplemented from TrigMuonMonitorAlgorithm.

Definition at line 15 of file L2MuonSAIOMon.cxx.

                                      {
 
  ATH_CHECK( TrigMuonMonitorAlgorithm::initialize() );
  ATH_CHECK( m_L2MuonCBIOContainerKey.initialize() );
 
  unsigned int nchains = m_monitored_chains.size();
  if(nchains!=m_monitored_chains_plateau.size()){
    ATH_MSG_ERROR("Configuration seems to be wrong. The size of \"Plateaus\" should be same as \"MonitoredChainds\".");
    return StatusCode::FAILURE;
  } else {
    for(unsigned int ichain=0; ichain<nchains; ++ichain){
      m_plateaus[m_monitored_chains[ichain]] = m_monitored_chains_plateau[ichain] + 2.;
    }
  }
  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 HiveAlgBase, 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;
}

isOverlap()

bool L2MuonSAIOMon::isOverlap ( const xAOD::L2CombinedMuon * matchSA_L2IOobject1, const xAOD::L2CombinedMuon * matchSA_L2IOobject2 ) const

private

Definition at line 753 of file L2MuonSAIOMon.cxx.

                                                                                                                                    {
 
  ATH_MSG_DEBUG( "   ...matchSA_L2IOobject1: pt/eta/phi=" << matchSA_L2IOobject1->pt()/Gaudi::Units::GeV << " / " << matchSA_L2IOobject1->eta() << " / " << matchSA_L2IOobject1->phi() );
  ATH_MSG_DEBUG( "   ...matchSA_L2IOobject2: pt/eta/phi=" << matchSA_L2IOobject2->pt()/Gaudi::Units::GeV << " / " << matchSA_L2IOobject2->eta() << " / " << matchSA_L2IOobject2->phi() );
 
  const auto [mu1Pt, mu1Eta, mu1Phi] = L2ORPosForMatchFunc(matchSA_L2IOobject1);
  const auto [mu2Pt, mu2Eta, mu2Phi] = L2ORPosForMatchFunc(matchSA_L2IOobject2);
 
  // if dR or invMass is necessary but (eta,phi) info is not avaiable
  // (i.e. eta,phi=0,0; rec failed)
  const double ZERO_LIMIT_FOR_ETAPHI = 1e-4;
  if( (std::abs(matchSA_L2IOobject1->eta()) <ZERO_LIMIT_FOR_ETAPHI && std::abs(matchSA_L2IOobject1->phi()) < ZERO_LIMIT_FOR_ETAPHI) ||
      (std::abs(matchSA_L2IOobject2->eta()) <ZERO_LIMIT_FOR_ETAPHI && std::abs(matchSA_L2IOobject2->phi()) < ZERO_LIMIT_FOR_ETAPHI) ) {
    ATH_MSG_DEBUG( "   ...-> (eta,phi) info not available (rec at (eta,phi)=(0,0))" );
    if( m_RequireDR || m_RequireMass ) {
      ATH_MSG_DEBUG( "   ...-> but dR of invMass check is required. cannot judge overlap -> return with false" );
      return false;
    }
  }
 
  // if charge or invMass is necessary but charge(=pT) info is not avaiable
  const double ZERO_LIMIT_FOR_PT = 1e-4;
  if( (std::abs(matchSA_L2IOobject1->pt()) <ZERO_LIMIT_FOR_PT) || (std::abs(matchSA_L2IOobject2->pt()) < ZERO_LIMIT_FOR_PT) ) {
    ATH_MSG_DEBUG( "   ...-> pT info not available (rec at pT=0)" );
    if( m_RequireSameSign || m_RequireMass ) {
      ATH_MSG_DEBUG( "   ...-> but same sign or invMass check is required. cannot judge overlap -> return with false" );
      return false;
    }
  }
 
 
  double absEta = (std::abs(mu1Pt) > std::abs(mu2Pt)) ? std::abs(mu1Eta) : std::abs(mu2Eta);
  unsigned int iThres=0;
  for(unsigned int i=0; i<(m_etaBins.size()-1); i++) {
    if ( m_etaBins[i] <= absEta && absEta < m_etaBins[i+1] ) iThres = i;
  }
  float dRThres   = m_dRCBThres[iThres];
  float dRbySAThres = m_dRbySAThres[iThres];
  float massThres = m_massCBThres[iThres];
  ATH_MSG_DEBUG( "   ...iThres=" << iThres );
  if(m_RequireDR) ATH_MSG_DEBUG( "   ...dR       threshold=" << dRThres     );
  if(m_RequireDRbySA) ATH_MSG_DEBUG( "   ...dR(byMF) threshold=" << dRbySAThres );
  if(m_RequireMass) ATH_MSG_DEBUG( "   ...mass     threshold=" << massThres   );
 
 
 
  // same sign cut
  bool sameSign = false;
  if( m_RequireSameSign ) {
    sameSign = ((mu1Pt*mu2Pt) > 0);
    ATH_MSG_DEBUG( "   ...-> sameSign=" << sameSign );
  }
 
  // dR cut
  bool dRisClose = false;
  float deta = mu1Eta - mu2Eta;
  float dphi = xAOD::P4Helpers::deltaPhi(mu1Phi, mu2Phi);
  float dR = std::sqrt(deta*deta + dphi*dphi);
  if( m_RequireDR ) {
    if( dR < dRThres ) dRisClose = true;
    ATH_MSG_DEBUG( "   ...-> dR=" << dR << " : dRisClose=" << dRisClose );
  }
 
  // dR(by L2SA) cut
  bool dRbySAisClose = false;
  const xAOD::L2StandAloneMuon* muSA1 = matchSA_L2IOobject1->muSATrack();
  const xAOD::L2StandAloneMuon* muSA2 = matchSA_L2IOobject2->muSATrack();
  if( m_RequireDRbySA ) {
    // here, we do not check (eta,phi) of mF is not (0,0)
    // (i.e. we apply muComb based cut even if muFast rec is failed)
    float deta = muSA1->etaMS() - muSA2->etaMS();
    float dphi = xAOD::P4Helpers::deltaPhi(muSA1->phiMS(), muSA2->phiMS());
    float dRBySA = std::sqrt(deta*deta + dphi*dphi);
    if( dRBySA < dRbySAThres ) dRbySAisClose = true;
    ATH_MSG_DEBUG( "   ...-> dR(by MF)=" << dRBySA << " : dRbySAisClose=" << dRbySAisClose );
  }
 
  // mass cut
  const double TRACK_MASS = 0.;  // just assume zero mass
  bool massIsClose = false;
  TLorentzVector lvioobj1, lvioobj2;
  lvioobj1.SetPtEtaPhiM(std::abs(mu1Pt), mu1Eta, mu1Phi, TRACK_MASS);
  lvioobj2.SetPtEtaPhiM(std::abs(mu2Pt), mu2Eta, mu2Phi, TRACK_MASS);
  TLorentzVector lvsum = lvioobj1 + lvioobj2;
  float invMass = lvsum.M();
  if( m_RequireMass ) {
    if( invMass < massThres ) massIsClose = true;
    ATH_MSG_DEBUG( "   ...-> invMass=" << invMass << " : massIsClose=" << massIsClose );
  }
 
 
  // total judge
  bool overlap = false;
  if( ((m_RequireSameSign &&   sameSign)     || (! m_RequireSameSign)) &&
      ((m_RequireDR       &&  dRisClose)     || (! m_RequireDR))       &&
      ((m_RequireDRbySA &&  dRbySAisClose) || (!  m_RequireDRbySA)) &&
      ((m_RequireMass     &&  massIsClose)   || (! m_RequireMass)) ) {
    overlap = true;
  }
  ATH_MSG_DEBUG( "   ...=> isOverlap=" << overlap );
 
  return overlap;
}

isPassedmuCombHypo()

StatusCode L2MuonSAIOMon::isPassedmuCombHypo ( const std::string & chain, const xAOD::L2CombinedMuon * Trig_L2IOobjects, bool & pass_muCombHypo ) const

private

Definition at line 962 of file L2MuonSAIOMon.cxx.

                                                                                                                                              {
  pass_muCombHypo = false;
 
  // config
  std::vector< float > my_EtaBins = {0, 1.05, 1.5, 2.0, 9.9};
  std::vector< float > my_muCombThres = {0., 0., 0., 0.};
  bool my_pikCuts = true;
  float my_maxPtToApplyPik = 25.;
  float my_chi2MaxID = 3.5;
  ATH_CHECK( decision_ptthreshold( chain, my_EtaBins, my_muCombThres, my_pikCuts, my_maxPtToApplyPik, my_chi2MaxID ) );
  bool pikCut = true;
  bool stdCut = true;
 
  auto ptValue = Trig_L2IOobject->pt() * Trig_L2IOobject->charge()/1e3;
  float fexPt = ptValue;
  if(my_pikCuts && (std::abs(fexPt) < my_maxPtToApplyPik)){
    if(Trig_L2IOobject->idTrack()->chiSquared() > my_chi2MaxID){
      ATH_MSG_DEBUG("this obj failed at Kpi rejection:idTrack_chiSquared =   " << Trig_L2IOobject->idTrack()->chiSquared() );
      pikCut = false;
    }
  }
 
  float absEta = std::abs(Trig_L2IOobject->eta());
  unsigned int iThres = 0;
  for(unsigned int i=0; i<(my_EtaBins.size()-1); i++) {
    if ( my_EtaBins[i] <= absEta && absEta < my_EtaBins[i+1] ) iThres = i;
  }
  const float muCombThres = my_muCombThres[iThres];
  if(Trig_L2IOobject->pt()/1e3 < muCombThres){
    ATH_MSG_DEBUG("this obj failed at std Pt cut:muCombThres =   " << muCombThres);
    stdCut = false;
  }
 
  if(stdCut && pikCut){
    ATH_MSG_DEBUG("this obj passed muCombhypo");
    pass_muCombHypo = true;
  }
  return StatusCode::SUCCESS;
}

L2ORPosForMatchFunc() [1/2]

std::tuple< float, float, float > L2MuonSAIOMon::L2ORPosForMatchFunc ( const xAOD::L2CombinedMuon * trig )

inlinestaticprivate

Definition at line 1033 of file L2MuonSAIOMon.cxx.

                                                                                              {
  return std::forward_as_tuple( (trig->pt()/1e3 * trig->charge() ), trig->eta(), trig->phi());
}

L2ORPosForMatchFunc() [2/2]

std::tuple< float, float, float > L2MuonSAIOMon::L2ORPosForMatchFunc ( const xAOD::L2StandAloneMuon * trig )

inlinestaticprivate

Definition at line 1028 of file L2MuonSAIOMon.cxx.

                                                                                                {
  return std::forward_as_tuple(trig->pt(), trig->etaMS(), trig->phiMS());
}

L2OverlapRemover()

StatusCode L2MuonSAIOMon::L2OverlapRemover ( const std::vector< const xAOD::L2CombinedMuon * > & matchSA_L2IOobjects, std::vector< bool > & isoverlap, std::vector< bool > & passOR ) const

private

Definition at line 676 of file L2MuonSAIOMon.cxx.

                                                                                                                                                                               {
 
  ATH_MSG_DEBUG ("L2OverlapRemover ..." );
 
  const size_t numMuon = matchSA_L2IOobjects.size();
  bool errorWhenIdentifyingOverlap = false;
 
  if(numMuon > 1){
    std::vector<unsigned int> mucombResult;
    //unsigned int i,j;
    for(unsigned int i=0; i<numMuon; i++) {mucombResult.emplace_back(i); }
    for(unsigned int i=0; i<numMuon-1; i++){
      for(unsigned int j=i+1; j<numMuon; j++){
    ATH_MSG_DEBUG("++ i=" << i << " vs j=" << j);
    bool overlapped = isOverlap(matchSA_L2IOobjects.at(i), matchSA_L2IOobjects.at(j));
    ATH_MSG_DEBUG("matchSA_L2IOobjects: i/j/Overlap = " << i << "/" << j << "/" << overlapped );
    if( ! overlapped ){ // judged as different
      ATH_MSG_DEBUG("   judged as: differenr objects");
      if( mucombResult[i] == mucombResult[j] ) { // but marked as same by someone
        ATH_MSG_DEBUG( "inconsistentency in muComb overlap removal for more than two objects" );
        ATH_MSG_DEBUG( "two objects are judged as different but both were already marked as identical by someone else as: " );
        ATH_MSG_DEBUG( "i/j/result[i]/result[j]=" << i << " / " << j << " / " << mucombResult[i] << " / "  << mucombResult[j] );
        errorWhenIdentifyingOverlap = true;
      }
    }
    else{ // judged as overlap
      if( (mucombResult[j] != j && mucombResult[i] != mucombResult[j]) || (mucombResult[j] == j && mucombResult[i] != i) ){
        ATH_MSG_DEBUG( "inconsistentency in muComb based overlap removal for more than two objects" );
        ATH_MSG_DEBUG( "two objects are judged as overlap but only either was already marked as overlap to someone else: " );
        ATH_MSG_DEBUG( "i/j/result[i]/result[j]=" << i << " / " << j << " / " << mucombResult[i] << " / "  << mucombResult[j] );
        errorWhenIdentifyingOverlap = true;
      }
      ATH_MSG_DEBUG("   judged as: overlapped objects");
      if( mucombResult[i] == i ) {
        ATH_MSG_DEBUG( "   i is not yet marked as overlap. so, it is a newly found overlap" );
        ATH_MSG_DEBUG( "   -> marking mucombResult[j] as i..." );
        mucombResult[j] = i;
        isoverlap[i] = true;
            isoverlap[j] = true;
      } else {
        ATH_MSG_DEBUG( "   both i/j already marked as overlap by: mucombResult[i]=" << mucombResult[i] );
        ATH_MSG_DEBUG( "   -> do nothing..." );
      }
    }
      }
    }
 
 
    if( errorWhenIdentifyingOverlap ) {
      ATH_MSG_WARNING( "error when resolving overlap. exitting with all EVs active..." );
    } else {
 
      unsigned int n_uniqueMuon = 0;
      for(unsigned int i=0; i<numMuon; i++) {
    ATH_MSG_DEBUG( "muComb based results: i=" << i << ": ");
    if( mucombResult[i] != i ) {
      ATH_MSG_DEBUG( "      overlap to j=" << mucombResult[i] );
    } else {
          n_uniqueMuon++;
      ATH_MSG_DEBUG( "      unique" );
    }
      }
 
      ATH_MSG_DEBUG( "nr of unique Muons after muComb-based removal=" << n_uniqueMuon );
 
      if( numMuon != n_uniqueMuon ){
        ATH_CHECK( chooseBestMuon(matchSA_L2IOobjects, passOR, mucombResult) );
      } else {
    ATH_MSG_DEBUG( "no overlap identified. exitting with all EventViews active" );
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

matchL2IO_wContainer()

StatusCode L2MuonSAIOMon::matchL2IO_wContainer ( const EventContext & ctx, const std::string & chain, std::vector< const xAOD::L2CombinedMuon * > & Trig_L2IOobjects ) const

private

Definition at line 610 of file L2MuonSAIOMon.cxx.

                                                                                                                                                                {
 
  ATH_MSG_DEBUG ("matchL2IO_wContainer ..." );
 
  // retrieve l2io objects
  SG::ReadHandle<xAOD::L2CombinedMuonContainer> L2IOobjects( m_L2MuonCBIOContainerKey, ctx );
 
  // retrieve l2SA objects
  std::vector< TrigCompositeUtils::LinkInfo<xAOD::L2StandAloneMuonContainer> > featureCont;
  if( chain.find("probe") != std::string::npos ){ // if tag & probe chain, retrieve probe L2SA objects
    int legIndex_probe = 1; // probe
    featureCont = getTrigDecisionTool()->features<xAOD::L2StandAloneMuonContainer>( chain,
                                            TrigDefs::includeFailedDecisions,
                                            "HLT_MuonL2SAInfo",
                                            TrigDefs::lastFeatureOfType,
                                            TrigCompositeUtils::featureString(),
                                            legIndex_probe );
  }
  else{
    featureCont = getTrigDecisionTool()->features<xAOD::L2StandAloneMuonContainer>( chain,
                                            TrigDefs::includeFailedDecisions,
                                            "HLT_MuonL2SAInfo" );
  }
 
  // match l2io objects to l2sa objects using roiWord
  std::vector< const xAOD::L2CombinedMuon* > matchSA_L2IOobjects;
  for(const auto L2IOobject : *L2IOobjects){
    ATH_MSG_DEBUG(" L2IOobject->muSATrack()->roiWord()/L2IOobject->pt(): " << L2IOobject->muSATrack()->roiWord() << "/" << L2IOobject->pt() );
    for(const TrigCompositeUtils::LinkInfo<xAOD::L2StandAloneMuonContainer>& L2SALinkInfo : featureCont){
      ATH_CHECK( L2SALinkInfo.isValid() );
      const ElementLink<xAOD::L2StandAloneMuonContainer> L2SAobject = L2SALinkInfo.link;
      if( !L2SAobject.isValid() ) continue;
      ATH_MSG_DEBUG(" L2SAobject->roiWord()/L2SALinkInfo.state: " << (*L2SAobject)->roiWord() << "/" << L2SALinkInfo.state );
 
      if( L2IOobject->muSATrack()->roiWord() != (*L2SAobject)->roiWord() ) continue;
      if( L2SALinkInfo.state != TrigCompositeUtils::ActiveState::ACTIVE ){
        break;
      }else{
        matchSA_L2IOobjects.push_back(L2IOobject);
        ATH_MSG_DEBUG(" matchSA_L2IOobject->muSATrack()->roiWord()/matchSA_L2IOobject->pt(): " << L2IOobject->muSATrack()->roiWord() << "/" << L2IOobject->pt() );
        break;
      }
    }
  }
  ATH_MSG_DEBUG(" matchSA_L2IOobjects.size(): " << matchSA_L2IOobjects.size() );
 
  const size_t num_matchSAMuon = matchSA_L2IOobjects.size();
  if( num_matchSAMuon == 0 ){
    ATH_MSG_DEBUG(" NO matchSA_L2IOobjects ");
    return StatusCode::SUCCESS;
  }
 
  std::vector< bool > isoverlap( num_matchSAMuon, false );
  std::vector< bool > passOR( num_matchSAMuon, true );
 
  // L2CBOverlapRemover
  ATH_CHECK( L2OverlapRemover( matchSA_L2IOobjects, isoverlap, passOR ) );
  for(unsigned int i=0; i<num_matchSAMuon; i++) { // push back trig Inside-Out objects passing L2CBOverlapRemover
    if( isoverlap[i] && !passOR[i] ) continue;
    Trig_L2IOobjects.push_back(matchSA_L2IOobjects.at(i));
  }
 
  return StatusCode::SUCCESS;
}

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

muCombHypo_TDTworkaround()

StatusCode L2MuonSAIOMon::muCombHypo_TDTworkaround ( const std::string & chain, const std::vector< const xAOD::L2CombinedMuon * > & Trig_L2IOobjects, std::vector< bool > & pass_muCombHypo ) const

private

Definition at line 941 of file L2MuonSAIOMon.cxx.

                                                                                                                                                                                      {
 
 
  int requireMuonNum = 1;
 
  int passHypo_MuonNum = 0;
  for(auto &Trig_L2IOobject : Trig_L2IOobjects){
    bool isPass_muCombHypo = false;
    ATH_CHECK(isPassedmuCombHypo( chain, Trig_L2IOobject ,isPass_muCombHypo));
    bool pass_muCombHypo_obj = isPass_muCombHypo;
    pass_muCombHypo.push_back(pass_muCombHypo_obj);
    if( pass_muCombHypo_obj ) passHypo_MuonNum++; 
  }
 
  if( passHypo_MuonNum >= requireMuonNum ){
    ATH_MSG_DEBUG("this evt passed muCombhypo");
  }
  return StatusCode::SUCCESS;
}

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 350 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();
  }

searchL2InsideOut()

const xAOD::L2CombinedMuon * L2MuonSAIOMon::searchL2InsideOut ( const EventContext & ctx, const xAOD::Muon * mu, const std::string & trigger ) const

private

Definition at line 1038 of file L2MuonSAIOMon.cxx.

                                                                                                                                      {
  ATH_MSG_DEBUG("MuonMonitoring::searchL2InsideOut()");
 
  const xAOD::L2CombinedMuon* offlinematched_L2IOobject = nullptr;
 
  //TDT workaround
  std::vector< const xAOD::L2CombinedMuon* > Trig_L2IOobjects;
  if( !matchL2IO_wContainer(ctx, trig, Trig_L2IOobjects).isSuccess() ) {
    ATH_MSG_WARNING("matchL2IO_wContainer failed, returning nullptr");
    return offlinematched_L2IOobject;
  }
  if( Trig_L2IOobjects.empty() ) {
    return offlinematched_L2IOobject;
  }
 
  float reqdR = 1000.;
 
  double offlEta = mu->eta();
  double offlPhi = mu->phi();
 
  int loop_counter = 0;
  int match_index = 0;
  for(auto Trig_L2IOobject : Trig_L2IOobjects){
    double trigEta = Trig_L2IOobject->eta();
    double trigPhi = Trig_L2IOobject->phi();
    double deta = offlEta - trigEta;
    double dphi = xAOD::P4Helpers::deltaPhi(offlPhi, trigPhi);
    double dR = std::sqrt(deta*deta + dphi*dphi);
 
    ATH_MSG_VERBOSE("Trigger muon candidate eta=" << trigEta << " phi=" << trigPhi  << " pt=" << Trig_L2IOobject->pt() << " dR=" << dR);
    if( dR<reqdR ){
      reqdR = dR;
      match_index = loop_counter;
      ATH_MSG_DEBUG("* Trigger muon eta=" << trigEta << " phi=" << trigPhi  << " pt=" << Trig_L2IOobject->pt() << " dR=" << dR );
    }
    loop_counter++;
  }
 
  offlinematched_L2IOobject = Trig_L2IOobjects.at(match_index);
  return offlinematched_L2IOobject;
}

searchTagOfflineMuon()

const xAOD::Muon * L2MuonSAIOMon::searchTagOfflineMuon ( const EventContext & ctx, const xAOD::Muon * probe ) const

private

Definition at line 1081 of file L2MuonSAIOMon.cxx.

                                                                                                           {
  ATH_MSG_DEBUG("MuonMonitoring::searchTagOfflineMuon()");
 
  const double ZERO_LIMIT = 0.00001;
 
  double Jpsimass = 3.0969;
  double Zmass = 91.1876;
  double my_Jpsimass_lowlim = 81.;
  double my_Jpsimass_highlim = 101.;
  double my_Zmass_lowlim = 2.7;
  double my_Zmass_highlim = 3.5;
 
  const xAOD::Muon *tag = nullptr;
 
  SG::ReadHandle<xAOD::MuonContainer> muons(m_MuonContainerKey, ctx);
  if (! muons.isValid() ) {
    ATH_MSG_ERROR("evtStore() does not contain muon Collection with name "<< m_MuonContainerKey);
    return tag;
  }
 
  double mass_diff_min = 999.;
  double tpdR_min = 999.;
  bool tpfromZ = false;
  for( const xAOD::Muon* mu : *muons ){
    if( mu->quality() != xAOD::Muon::Quality::Medium && 
        mu->quality() != xAOD::Muon::Quality::Tight ) continue;
    if( mu->charge()*probe->charge() > 0 ) continue;
    const auto* tag_ms_track = mu->trackParticle(xAOD::Muon::TrackParticleType::MuonSpectrometerTrackParticle);
    if( !tag_ms_track) continue; // tag muon dosen't have ms track 
    TLorentzVector lvmu = mu->p4();
    TLorentzVector lvprobe = probe->p4();
    double dimu_mass = (lvmu+lvprobe).M()/1.e3;
    double tpdR = lvmu.DeltaR(lvprobe);
    if( dimu_mass > my_Jpsimass_lowlim && dimu_mass < my_Jpsimass_highlim ){
      if( tpfromZ ) continue; // Z has higher priority than Jpsi
      double mass_diff = std::abs(dimu_mass - Jpsimass);
      if( mass_diff - mass_diff_min < -1.*ZERO_LIMIT ){
    mass_diff_min = mass_diff;
    tpdR_min =tpdR;
    tag = mu;
      }
      else if( std::abs(mass_diff - mass_diff_min) < ZERO_LIMIT){
    if( tpdR - tpdR_min < 0. ){
      mass_diff_min = mass_diff;
      tpdR_min = tpdR;
      tag = mu;
    }
      }
    }
    else if( dimu_mass > my_Zmass_lowlim && dimu_mass < my_Zmass_highlim ){
      tpfromZ = true;
      double mass_diff = std::abs(dimu_mass - Zmass);
      if( mass_diff - mass_diff_min < -1.*ZERO_LIMIT ){
    mass_diff_min = mass_diff;
        tpdR_min =tpdR;
    tag = mu;
      } 
      else if( std::abs(mass_diff - mass_diff_min) < ZERO_LIMIT){
        if( tpdR - tpdR_min < 0. ){
          mass_diff_min = mass_diff;
          tpdR_min =tpdR;
          tag = mu;
        }
      }
    }
  }
 
  return tag;
}

selectEvents()

bool TrigMuonMonitorAlgorithm::selectEvents ( ) const

protectedvirtualinherited

Function that defines the event selection for anlayses User should reimlement in a subclass if needed.

Returns

True if the event is used for an analysis.

Reimplemented in TrigMuonEfficiencyMon.

Definition at line 65 of file TrigMuonMonitorAlgorithm.cxx.

                                                    {
  return true;
}

selectMuons()

StatusCode TrigMuonMonitorAlgorithm::selectMuons ( const EventContext & ctx, SG::ReadHandle< xAOD::MuonContainer > & muons, std::vector< const xAOD::Muon * > & probes ) const

protectedvirtualinherited

Function that defines the event selection for anlayses Users should reimlement in a subclass if needed.

Parameters

ctx	EventContext provided by athenaMT
muons	Offline muons in the MuonContainer
probes	List of offline muons that are used in analyses

Reimplemented in TrigMuonEfficiencyMon.

Definition at line 70 of file TrigMuonMonitorAlgorithm.cxx.

                                                                                                                                                          {
  for (const xAOD::Muon* mu : *muons) {
    if(mu->author()==xAOD::Muon::Author::MuidCo && 
       m_muonSelectionTool->getQuality(*mu)<=xAOD::Muon::Quality::Medium ){
      probes.push_back(mu);
    }
  }
  return StatusCode::SUCCESS;
}

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 208 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_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_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_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_dRbySAThres

Gaudi::Property< std::vector<float> > L2MuonSAIOMon::m_dRbySAThres {this, "dRbySAThres", {0, 0, 0, 0, 0}, "mufast DR threshold of L2CB"}

private

Definition at line 48 of file L2MuonSAIOMon.h.

{this, "dRbySAThres", {0, 0, 0, 0, 0}, "mufast DR threshold of L2CB"};

m_dRCBThres

Gaudi::Property< std::vector<float> > L2MuonSAIOMon::m_dRCBThres {this, "DRThres", {0.002, 0.001, 0.002, 0.002, 0.002}, "DR threshold of L2CB"}

private

Definition at line 47 of file L2MuonSAIOMon.h.

{this, "DRThres", {0.002, 0.001, 0.002, 0.002, 0.002}, "DR threshold of L2CB"};

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_etaBins

Gaudi::Property< std::vector<float> > L2MuonSAIOMon::m_etaBins {this, "EtaBins", {0, 0.9, 1.1, 1.9, 2.1, 9.9}, "eta bins of DR and mass thresholds"}

private

Definition at line 46 of file L2MuonSAIOMon.h.

{this, "EtaBins", {0, 0.9, 1.1, 1.9, 2.1, 9.9}, "eta bins of DR and mass thresholds"};

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_group

Gaudi::Property<std::string> TrigMuonMonitorAlgorithm::m_group {this, "Group", "", "Histogram group"}

protectedinherited

Name of monitored group.

Definition at line 144 of file TrigMuonMonitorAlgorithm.h.

{this, "Group", "", "Histogram group"};

m_L2MuonCBIOContainerKey

SG::ReadHandleKey<xAOD::L2CombinedMuonContainer> L2MuonSAIOMon::m_L2MuonCBIOContainerKey {this, "L2CombinedMuonContainerName_IO", "HLT_MuonL2CBInfoIOmode", "L2MuonCBIO container"}

private

Definition at line 29 of file L2MuonSAIOMon.h.

{this, "L2CombinedMuonContainerName_IO", "HLT_MuonL2CBInfoIOmode", "L2MuonCBIO container"};

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_massCBThres

Gaudi::Property< std::vector<float> > L2MuonSAIOMon::m_massCBThres {this, "MassThres", {0.004, 0.002, 0.006, 0.006, 0.006}, "mass threshold of L2CB"}

private

Definition at line 49 of file L2MuonSAIOMon.h.

{this, "MassThres", {0.004, 0.002, 0.006, 0.006, 0.006}, "mass threshold of L2CB"};

m_matchTool

ToolHandle<MuonMatchingTool> TrigMuonMonitorAlgorithm::m_matchTool {this, "MuonMatchingTool", "MuonMatchingTool", "Tool for matching offline and online objects"}

protectedinherited

Definition at line 130 of file TrigMuonMonitorAlgorithm.h.

{this, "MuonMatchingTool", "MuonMatchingTool", "Tool for matching offline and online objects"};

m_monitored_chains

Gaudi::Property<std::vector<std::string> > TrigMuonMonitorAlgorithm::m_monitored_chains {this, "MonitoredChains", {}, "Trigger chains that are monitored"}

protectedinherited

List of trigger chains that are monitored in fillVariablesPerChain and fillVariablesPerOfflineMuonPerChain.

Definition at line 138 of file TrigMuonMonitorAlgorithm.h.

{this, "MonitoredChains", {}, "Trigger chains that are monitored"};

m_monitored_chains_plateau

Gaudi::Property<std::vector<float> > L2MuonSAIOMon::m_monitored_chains_plateau {this, "Plateaus", {}, "List of plateaus of measured trigger"}

private

Definition at line 27 of file L2MuonSAIOMon.h.

{this, "Plateaus", {}, "List of plateaus of measured trigger"};

m_MuonContainerKey

SG::ReadHandleKey<xAOD::MuonContainer> TrigMuonMonitorAlgorithm::m_MuonContainerKey {this, "MuonContainerName", "Muons", "Offline muon container"}

protectedinherited

Definition at line 134 of file TrigMuonMonitorAlgorithm.h.

{this, "MuonContainerName", "Muons", "Offline muon container"};

m_muonSelectionTool

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

protectedinherited

Definition at line 131 of file TrigMuonMonitorAlgorithm.h.

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

m_muontype

Gaudi::Property<int> TrigMuonMonitorAlgorithm::m_muontype

protectedinherited

Initial value:

{this, "MuonType", 
                                   static_cast<int>(xAOD::Muon::MuonType::Combined), 
                                   "MuonType used for monitoring"}

Requirement for the offline muon type considered in analyses.

Definition at line 140 of file TrigMuonMonitorAlgorithm.h.

                                {this, "MuonType", 
                                   static_cast<int>(xAOD::Muon::MuonType::Combined), 
                                   "MuonType used for monitoring"};

m_name

std::string AthMonitorAlgorithm::m_name

privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

m_plateaus

std::map<std::string, double> L2MuonSAIOMon::m_plateaus {}

private

Definition at line 28 of file L2MuonSAIOMon.h.

{};

m_ratio_measurement_threshold

const float TrigMuonMonitorAlgorithm::m_ratio_measurement_threshold = 4

protectedinherited

Threshold for ratio measurement.

Definition at line 147 of file TrigMuonMonitorAlgorithm.h.

m_RequireDR

bool L2MuonSAIOMon::m_RequireDR = true

private

Definition at line 52 of file L2MuonSAIOMon.h.

m_RequireDRbySA

bool L2MuonSAIOMon::m_RequireDRbySA = false

private

Definition at line 53 of file L2MuonSAIOMon.h.

m_RequireMass

bool L2MuonSAIOMon::m_RequireMass = true

private

Definition at line 54 of file L2MuonSAIOMon.h.

m_RequireSameSign

bool L2MuonSAIOMon::m_RequireSameSign = true

private

Definition at line 55 of file L2MuonSAIOMon.h.

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:

• L2MuonSAIOMon.h
• L2MuonSAIOMon.cxx