RpcTrackAnaAlg Class Reference

#include <RpcTrackAnaAlg.h>

Inheritance diagram for RpcTrackAnaAlg:

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

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

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

Private Types
enum	BarrelDL {   BI = 1 , BM1 , BM2 , BO1 ,   BO2 , OUT }
enum	MuonSource { AllMuon = 0 , ZCand }
typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > >	MonVarVec_t
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	readElIndexFromXML ()
StatusCode	initRpcPanel ()
StatusCode	setPanelIndex (std::shared_ptr< RpcPanel > panel)
StatusCode	initTrigTag ()
StatusCode	initArrayHistosMap ()
StatusCode	fillMuonExtrapolateEff (const EventContext &ctx) const
StatusCode	fillHistPRD (const EventContext &ctx) const
StatusCode	triggerMatching (const xAOD::Muon *, const std::vector< TagDef > &) const
StatusCode	extrapolate2RPC (const xAOD::TrackParticle *track, const Trk::PropDirection direction, std::vector< GasGapResult > &results, BarrelDL barrelDL) const
std::unique_ptr< Trk::TrackParameters >	computeTrackIntersectionWithGasGap (ExResult &result, const xAOD::TrackParticle *track_particle, const std::shared_ptr< GasGapData > &gap) const
StatusCode	extrapolate2RPC (std::unique_ptr< Trk::TrackParameters > trackParam, const Trk::PropDirection direction, std::vector< GasGapResult > &results, BarrelDL barrelDL) const
std::unique_ptr< Trk::TrackParameters >	computeTrackIntersectionWithGasGap (ExResult &result, const Trk::TrackParameters *trackParam, const std::shared_ptr< GasGapData > &gap) const
StatusCode	readHitsPerGasgap (const EventContext &ctx, std::vector< GasGapResult > &results, MuonSource muon_source) const
StatusCode	fillClusterSize (std::vector< const Muon::RpcPrepData * > &view_hits, const int panel_index, int LB, int phiSector, int isPhi) const
bool	IsNearbyHit (const std::vector< const Muon::RpcPrepData * > &cluster_hits, const Muon::RpcPrepData *hit) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
BooleanProperty	m_plotMuonEff
BooleanProperty	m_plotPRD
BooleanProperty	m_useAODParticle
FloatProperty	m_lbDuraThr
StringProperty	m_packageName
StringProperty	m_elementsFileName
StringProperty	m_trigTagList
FloatProperty	m_trigMatchWindow
BooleanProperty	m_TagAndProbe
BooleanProperty	m_TagAndProbeZmumu
FloatProperty	m_minPt {this, "minPt", 25.0e3, "minmum pT of muon"}
FloatProperty	m_maxEta
FloatProperty	m_barrelMinPt {this, "barrelMinPt", 2.0e3}
FloatProperty	m_barrelMinEta {this, "barrelMinEta", 0.1}
FloatProperty	m_barrelMaxEta {this, "barrelMaxEta", 1.05}
FloatProperty	m_muonMass
FloatProperty	m_zMass_lowLimit
FloatProperty	m_zMass_upLimit
FloatProperty	m_isolationWindow
FloatProperty	m_l1trigMatchWindow
FloatProperty	m_minDRTrackToGasGap
FloatProperty	m_boundsToleranceReadoutElement
FloatProperty	m_boundsToleranceReadoutElementTight
FloatProperty	m_diffHitTrackPostion
FloatProperty	m_outtime {this, "outtime", 12.5, "the out-time time"}
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc
const RpcIdHelper *	m_rpcIdHelper {nullptr}
ToolHandle< Trk::IExtrapolator >	m_extrapolator
SG::ReadHandleKey< xAOD::MuonRoIContainer >	m_MuonRoIContainerKey
SG::ReadHandleKey< xAOD::MuonContainer >	m_MuonContainerKey
SG::ReadHandleKey< Muon::RpcPrepDataContainer >	m_rpcPrdKey
SG::ReadHandleKey< xAOD::VertexContainer >	m_PrimaryVertexContainerKey
SG::ReadDecorHandleKey< xAOD::EventInfo >	m_beamSigmaX
SG::ReadDecorHandleKey< xAOD::EventInfo >	m_beamSigmaY
SG::ReadDecorHandleKey< xAOD::EventInfo >	m_beamSigmaXY
RpcPanelMap	m_rpcPanelMap {}
std::vector< TagDef >	m_trigTagDefs {}
std::map< std::pair< int, int >, std::vector< std::shared_ptr< GasGapData > > >	m_gasGapData {}
std::map< BarrelDL, std::vector< int > >	m_StationNames {}
std::map< std::string, int >	m_elementIndex {}
std::map< std::string, int >	m_SectorGroup {}
std::map< std::string, int >	m_TriggerThrGroup {}
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 39 of file RpcTrackAnaAlg.h.

Member Typedef Documentation

GasGapResult

typedef std::pair<ExResult, const std::shared_ptr<GasGapData> > RpcTrackAnaAlg::GasGapResult

Definition at line 48 of file RpcTrackAnaAlg.h.

MonVarVec_t

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

privateinherited

Definition at line 370 of file AthMonitorAlgorithm.h.

RpcPanelMap

typedef std::map<Identifier, std::shared_ptr<RpcPanel> > RpcTrackAnaAlg::RpcPanelMap

Definition at line 47 of file RpcTrackAnaAlg.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

BarrelDL

enum RpcTrackAnaAlg::BarrelDL

private

Enumerator
BI
BM1
BM2
BO1
BO2
OUT

Definition at line 51 of file RpcTrackAnaAlg.h.

                  {
        BI = 1,
        BM1,
        BM2,
        BO1,
        BO2,
        OUT
    };  // Barrel doublet: BM_dbR

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

MuonSource

enum RpcTrackAnaAlg::MuonSource

private

Enumerator
AllMuon
ZCand

Definition at line 59 of file RpcTrackAnaAlg.h.

{ AllMuon = 0, ZCand };

Constructor & Destructor Documentation

RpcTrackAnaAlg()

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

Definition at line 30 of file RpcTrackAnaAlg.cxx.

    : AthMonitorAlgorithm(name, pSvcLocator) {}

~RpcTrackAnaAlg()

RpcTrackAnaAlg::~RpcTrackAnaAlg ( )

virtual

Definition at line 35 of file RpcTrackAnaAlg.cxx.

{}

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

computeTrackIntersectionWithGasGap() [1/2]

std::unique_ptr< Trk::TrackParameters > RpcTrackAnaAlg::computeTrackIntersectionWithGasGap ( ExResult & result, const Trk::TrackParameters * trackParam, const std::shared_ptr< GasGapData > & gap ) const

private

Definition at line 1035 of file RpcTrackAnaAlg.cxx.

                                                {
    const EventContext& ctx = Gaudi::Hive::currentContext();
    /*
      This function:
      - constructs Identifier for specific gasgap
      - extrapolates muon to this gas gap
    */
 
    // Get surface of this gas gap and extrapolate track to this surface
    const Trk::SurfaceBounds& bounds = gap->readoutEl->bounds(gap->gapid);
    const Trk::PlaneSurface& gapSurface = gap->readoutEl->surface(gap->gapid);
    auto detParameters = m_extrapolator->extrapolate(
        ctx, *trackParam, gapSurface, result.direction, true, Trk::muon);
 
    if (!detParameters) {
        return detParameters;
    }
 
    //
    // Transform global extrapolated track position on surface to local
    // coordinates
    //
    const Amg::Vector3D local3dTrackPosition =
        gap->readoutEl->globalToLocalCoords(detParameters->position(),
                                            gap->gapid);
    const Amg::Vector2D local2dTrackPosition(local3dTrackPosition.y(),
                                             local3dTrackPosition.z());
 
    //
    // Check if the track position on surface is within tolerable bounds
    //
    const bool inbounds =
        bounds.inside(local2dTrackPosition, m_boundsToleranceReadoutElement,
                      m_boundsToleranceReadoutElement);
    const bool inbounds_tight = bounds.inside(
        local2dTrackPosition, m_boundsToleranceReadoutElementTight,
        m_boundsToleranceReadoutElementTight);
 
    result.localTrackPosInBounds = inbounds;
    result.localTrackPosInBoundsTight = inbounds_tight;
    result.localPos = local3dTrackPosition;
    result.globalPos = detParameters->position();
 
    return detParameters;
}

computeTrackIntersectionWithGasGap() [2/2]

std::unique_ptr< Trk::TrackParameters > RpcTrackAnaAlg::computeTrackIntersectionWithGasGap ( ExResult & result, const xAOD::TrackParticle * track_particle, const std::shared_ptr< GasGapData > & gap ) const

private

Definition at line 866 of file RpcTrackAnaAlg.cxx.

                                                {
    const EventContext& ctx = Gaudi::Hive::currentContext();
    /*
      This function:
      - constructs Identifier for specific gasgap
      - extrapolates muon to this gas gap
    */
 
    // Get surface of this gas gap and extrapolate track to this surface
    const Trk::SurfaceBounds& bounds = gap->readoutEl->bounds(gap->gapid);
    const Trk::PlaneSurface& gapSurface = gap->readoutEl->surface(gap->gapid);
    std::unique_ptr<Trk::TrackParameters> detParameters{};
 
    if (m_useAODParticle) {
        detParameters = m_extrapolator->extrapolate(
            ctx, track_particle->perigeeParameters(), gapSurface,
            result.direction, false, Trk::muon);
    } else if (track_particle->track()) {
        detParameters = m_extrapolator->extrapolateTrack(
            ctx, *(track_particle->track()), gapSurface, result.direction, true,
            Trk::muon);
    } else {
        return detParameters;
    }
 
    if (!detParameters) {
        return detParameters;
    }
 
    //
    // Transform global extrapolated track position on surface to local
    // coordinates
    //
    const Amg::Vector3D local3dTrackPosition =
        gap->readoutEl->globalToLocalCoords(detParameters->position(),
                                            gap->gapid);
    const Amg::Vector2D local2dTrackPosition(local3dTrackPosition.y(),
                                             local3dTrackPosition.z());
 
    //
    // Check if the track position on surface is within tolerable bounds
    //
    const bool inbounds =
        bounds.inside(local2dTrackPosition, m_boundsToleranceReadoutElement,
                      m_boundsToleranceReadoutElement);
    const bool inbounds_tight = bounds.inside(
        local2dTrackPosition, m_boundsToleranceReadoutElementTight,
        m_boundsToleranceReadoutElementTight);
 
    result.localTrackPosInBounds = inbounds;
    result.localTrackPosInBoundsTight = inbounds_tight;
    result.localPos = local3dTrackPosition;
    result.globalPos = detParameters->position();
 
    return detParameters;
}

dataType()

DataType_t AthMonitorAlgorithm::dataType ( ) const

inlineinherited

Accessor functions for the data type.

Returns

the current value of the class's DataType_t instance.

Definition at line 224 of file AthMonitorAlgorithm.h.

{ return m_dataType; }

dataTypeStringToEnum()

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

inherited

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

Returns

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

Definition at line 144 of file AthMonitorAlgorithm.cxx.

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

declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

environment()

Environment_t AthMonitorAlgorithm::environment ( ) const

inlineinherited

Accessor functions for the environment.

Returns

the current value of the class's Environment_t instance.

Definition at line 208 of file AthMonitorAlgorithm.h.

{ return m_environment; }

envStringToEnum()

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

inherited

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

Returns

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

Definition at line 116 of file AthMonitorAlgorithm.cxx.

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

evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

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

overridevirtualinherited

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters

ctx	event context for reentrant Athena call

Returns

StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

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

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

extraOutputDeps()

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

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

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

extrapolate2RPC() [1/2]

StatusCode RpcTrackAnaAlg::extrapolate2RPC ( const xAOD::TrackParticle * track, const Trk::PropDirection direction, std::vector< GasGapResult > & results, BarrelDL barrelDL ) const

private

Definition at line 756 of file RpcTrackAnaAlg.cxx.

                                                                    {
    /*
      get intersections of the muon with the RPC planes
 
      Iterate over all valid RPCDetectorElements and RPCReadoutElements:
       1) compute DR distance between track and center of ReadoutElement
          if this distance within tolerance - proceed
       2) Next, compute:
            -- min DR distance between track and strips within this gas gap
            -- number of valid eta and phi strips within this gas gap
          if both results within their tolerances - proceed
       3) Extrapolate track to the surface of this gas gap:
            -- Check that extrapolation result valid
            -- Check that extrapolated position is in the gas gap surface bounds
          if both within checks valid - then save RPC extrapolation result
    */
    int doubletR = 1;
    if (barrelDL >= OUT) {
        return StatusCode::SUCCESS;
    } else if (barrelDL == BM2 || barrelDL == BO2) {
        doubletR = 2;
    }
 
    using namespace Monitored;
    auto tool = getGroup(m_packageName);
 
    if (!track) {
        return StatusCode::FAILURE;
    }
 
    std::map<BarrelDL, std::vector<int>>::const_iterator dl_vec_it =
        m_StationNames.find(barrelDL);
    if (dl_vec_it == m_StationNames.end()) {
        return StatusCode::FAILURE;
    }
 
    std::unique_ptr<Trk::TrackParameters> trackParamLayer{};
    double minDR = 1.0;  // A intial value
 
    const std::vector<int> dl_vec = dl_vec_it->second;
    std::vector<int>::const_iterator it_dl = dl_vec.begin();
    for (; it_dl != dl_vec.end(); ++it_dl) {
        int stName = *it_dl;
        std::pair<int, int> st_dbR = std::make_pair(stName, doubletR);
        std::map<std::pair<int, int>,
                 std::vector<std::shared_ptr<GasGapData>>>::const_iterator
            gasgapIt = m_gasGapData.find(st_dbR);
        if (gasgapIt == m_gasGapData.end()) {
            continue;
        }
 
        //
        // Iterate over RPC readout elements and compute intersections with each
        // gas gap
        //
        for (const std::shared_ptr<GasGapData>& gap : gasgapIt->second) {
            ExResult result(gap->gapid, direction);
 
            // Compute track distance to the gas gap surface
            gap->computeTrackDistanceToGasGap(result, *track);
 
            if (result.minTrackGasGapDR > m_minDRTrackToGasGap) {
                continue;
            }
 
            //
            // Extrapolate track to the gas gap surface and check whether the
            // track position is in bounds returns new object
            auto trackParamInGap =
                computeTrackIntersectionWithGasGap(result, track, gap);
 
            if (!trackParamInGap) {
                continue;
            }
 
            if (!result.localTrackPosInBoundsTight) {
                continue;
            }
 
            if (result.minTrackGasGapDR < minDR) {
                minDR = result.minTrackGasGapDR;
                // new object moved to trackParamLayer; previous trackParamLayer
                // gets destroyed
                trackParamLayer = std::move(trackParamInGap);
            }
            ATH_MSG_DEBUG(name() << " extrapolated gasgap: " << gap->gapid_str);
 
            results.push_back(std::make_pair(result, gap));
        }
    }
 
    // Go to next layer of doublet
    BarrelDL nextDL = BarrelDL(barrelDL + 1);
 
    // propgate the track parameter of the last doublet
    // if no track paramater, use the input track
    if (trackParamLayer != nullptr) {
        // trackParamLayer used and then goes out of scope to destroy the object
        return extrapolate2RPC(std::move(trackParamLayer), direction, results,
                               nextDL);
    } else {
        return extrapolate2RPC(track, direction, results, nextDL);
    }
}

extrapolate2RPC() [2/2]

StatusCode RpcTrackAnaAlg::extrapolate2RPC ( std::unique_ptr< Trk::TrackParameters > trackParam, const Trk::PropDirection direction, std::vector< GasGapResult > & results, BarrelDL barrelDL ) const

private

Definition at line 926 of file RpcTrackAnaAlg.cxx.

                             {
    /*
      get intersections of the muon with the RPC planes
 
      Iterate over all valid RPCDetectorElements and RPCReadoutElements:
       1) compute DR distance between track and center of ReadoutElement
          if this distance within tolerance - proceed
       2) Next, compute:
            -- min DR distance between track and strips within this gas gap
            -- number of valid eta and phi strips within this gas gap
          if both results within their tolerances - proceed
       3) Extrapolate track to the surface of this gas gap:
            -- Check that extrapolation result valid
            -- Check that extrapolated position is in the gas gap surface bounds
          if both within checks valid - then save RPC extrapolation result
    */
    int doubletR = 1;
    if (barrelDL >= OUT) {
        return StatusCode::SUCCESS;
    } else if (barrelDL == BM2 || barrelDL == BO2) {
        doubletR = 2;
    }
 
    using namespace Monitored;
    auto tool = getGroup(m_packageName);
 
    if (!trackParam) {
        return StatusCode::FAILURE;
    }
 
    std::map<BarrelDL, std::vector<int>>::const_iterator dl_vec_it =
        m_StationNames.find(barrelDL);
    if (dl_vec_it == m_StationNames.end()) {
        return StatusCode::FAILURE;
    }
 
    std::unique_ptr<Trk::TrackParameters> trackParamLayer{};
    double minDR = 1.0;  // A intial value
 
    const std::vector<int> dl_vec = dl_vec_it->second;
 
    std::vector<int>::const_iterator it_dl = dl_vec.begin();
    for (; it_dl != dl_vec.end(); ++it_dl) {
        int stName = *it_dl;
        std::pair<int, int> st_dbR = std::make_pair(stName, doubletR);
        std::map<std::pair<int, int>,
                 std::vector<std::shared_ptr<GasGapData>>>::const_iterator
            gasgapIt = m_gasGapData.find(st_dbR);
 
        if (gasgapIt == m_gasGapData.end()) {
            continue;
        }
 
        //
        // Iterate over RPC readout elements and compute intersections with each
        // gas gap
        //
        for (const std::shared_ptr<GasGapData>& gap : gasgapIt->second) {
            ExResult result(gap->gapid, direction);
 
            // Compute track distance to the gas gap surface; doesnt take
            // ownership
            gap->computeTrackDistanceToGasGap(result, trackParam.get());
 
            if (result.minTrackGasGapDR > m_minDRTrackToGasGap) {
                continue;
            }
 
            //
            // Extrapolate track to the gas gap surface and check whether the
            // track position is in bounds; doesnt take ownership but returns a
            // new object
            auto trackParamInGap = computeTrackIntersectionWithGasGap(
                result, trackParam.get(), gap);
            if (trackParamInGap == nullptr) {
                continue;
            }
 
            if (!result.localTrackPosInBoundsTight) {
                continue;
            }
 
            ATH_MSG_DEBUG(name() << " extrapolated gasgap: " << gap->gapid_str);
 
            if (result.minTrackGasGapDR < minDR) {
                minDR = result.minTrackGasGapDR;
                // previously created trackParamInGap moved to trackParamLayer;
                // previous trackParamLayer deleted
                trackParamLayer = std::move(trackParamInGap);
            }
 
            results.push_back(std::make_pair(result, gap));
        }
    }
 
    if (trackParamLayer == nullptr) {
        trackParamLayer = std::move(trackParam);
    }
    BarrelDL nextDL = BarrelDL(barrelDL + 1);
    // trackParamLayer used and then goes out of scope, destroying the object
    return extrapolate2RPC(std::move(trackParamLayer), direction, results,
                           nextDL);
}

fillClusterSize()

StatusCode RpcTrackAnaAlg::fillClusterSize ( std::vector< const Muon::RpcPrepData * > & view_hits, const int panel_index, int LB, int phiSector, int isPhi ) const

private

Definition at line 1343 of file RpcTrackAnaAlg.cxx.

                                            {
    using namespace Monitored;
 
    auto tool = getGroup(m_packageName);
    auto i_LB = Scalar<int>("LB_nrpchit", LB);
 
    // Make clusters from hits that are close together in space and time
    std::vector<const Muon::RpcPrepData*> cluster_hits;
    while (!view_hits.empty()) {
        cluster_hits.clear();
 
        // Seed cluster with first (random) hit
        cluster_hits.push_back(view_hits.back());
 
        // Erase the selected first hit from the list
        view_hits.pop_back();
 
        // Collect all other hits which are close to the selected hit in time
        // and space
        std::vector<const Muon::RpcPrepData*>::const_iterator hit =
            view_hits.begin();
 
        while (hit != view_hits.end()) {
            const Muon::RpcPrepData* hit_ptr = *hit;
 
            if (IsNearbyHit(cluster_hits, hit_ptr)) {
                cluster_hits.push_back(*hit);
                view_hits.erase(hit);
 
                // Start loop from the beginning since we have increased cluster
                // size
                hit = view_hits.begin();
            } else {
                ++hit;
            }
        }
 
        int cluster_size = cluster_hits.size();
        for (int i_hit = 0; i_hit < cluster_size; i_hit++) {
            auto i_phiSector = Scalar<int>("PhiSector", phiSector);
            fill(tool, i_LB, i_phiSector);
        }
 
        auto i_panelIndex = Scalar<int>("panelInd_clust", panel_index);
        auto i_clusterSize = Scalar<int>("clusterSize", cluster_size);
        fill(tool, i_panelIndex, i_clusterSize);
 
        auto i_cs_sec = Scalar<int>("cs_sec", cluster_size);
        fill(m_tools[m_SectorGroup.at("sector" +
                                      std::to_string(std::abs(phiSector)))],
             i_cs_sec);
 
        if (isPhi == 1) {
            auto i_clusterSize_view =
                Scalar<int>("clusterSize_phi", cluster_size);
            fill(tool, i_clusterSize_view);
        } else {
            auto i_clusterSize_view =
                Scalar<int>("clusterSize_eta", cluster_size);
            fill(tool, i_clusterSize_view);
        }
    }
 
    return StatusCode::SUCCESS;
}

fillHistograms()

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

overridevirtual

adds event to the monitoring histograms

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

Parameters

ctx	forwarded from execute

Returns

StatusCode

Implements AthMonitorAlgorithm.

Definition at line 306 of file RpcTrackAnaAlg.cxx.

                                                                       {
    using namespace Monitored;
 
    if (m_plotMuonEff) {
        ATH_CHECK(fillMuonExtrapolateEff(ctx));
    }
 
    if (m_plotPRD) {
        ATH_CHECK(fillHistPRD(ctx));
    }
 
    auto tool = getGroup(m_packageName);
    auto evtLB = Scalar<int>("evtLB", ctx.eventID().lumi_block());
    auto run = Scalar<int>("run", ctx.eventID().run_number());
    fill(tool, evtLB, run);
 
    return StatusCode::SUCCESS;
}

fillHistPRD()

StatusCode RpcTrackAnaAlg::fillHistPRD ( const EventContext & ctx ) const

private

Definition at line 616 of file RpcTrackAnaAlg.cxx.

                                                                    {
    using namespace Monitored;
    //
    // Read RPC Prepare data
    //
 
    SG::ReadHandle<Muon::RpcPrepDataContainer> rpcContainer(m_rpcPrdKey, ctx);
    const RpcIdHelper& rpcIdHelper = m_idHelperSvc->rpcIdHelper();
 
    const int i_lb = ctx.eventID().lumi_block();
    std::vector<double> v_prdTime = {};
 
    auto prd_sec_all = Scalar<int>("prd_sec", 0);
    auto prd_layer_all = Scalar<int>("prd_layer", 0);
    auto prd_sec_1214 = Scalar<int>("prd_sec_1214", 0);
    auto prd_layer_1214 = Scalar<int>("prd_layer_1214", 0);
 
    auto prd_sec_all_eta = Scalar<int>("prd_sec_eta", 0);
    auto prd_layer_all_eta = Scalar<int>("prd_layer_eta", 0);
    auto prd_sec_all_phi = Scalar<int>("prd_sec_phi", 0);
    auto prd_layer_all_phi = Scalar<int>("prd_layer_phi", 0);
 
    auto i_prd_LB = Scalar<int>("LB", i_lb);
    auto i_panelIndex = Scalar<int>("panelInd", -1);
 
    auto tool = getGroup(m_packageName);
 
    int panel_index;
    std::pair<int, int> sec_layer;
 
    //
    // loop on RpcPrepData container
    //
    for (const Muon::RpcPrepDataCollection* rpcCollection : *rpcContainer) {
        //
        // loop on RpcPrepData
        //
        for (const Muon::RpcPrepData* rpcData : *rpcCollection) {
 
            Identifier id = rpcData->identify();
            const int measphi = rpcIdHelper.measuresPhi(id);
 
            auto temp_panel = std::make_unique<RpcPanel>(
                    *m_idHelperSvc, rpcData->detectorElement(), 
                    rpcIdHelper.doubletZ(id), rpcIdHelper.doubletPhi(id), 
                    rpcIdHelper.gasGap(id), rpcIdHelper.measuresPhi(id));
 
            std::map<Identifier, std::shared_ptr<RpcPanel>>::const_iterator
                i_panel = m_rpcPanelMap.find(temp_panel->panelId);
            if (i_panel == m_rpcPanelMap.end()) {
                ATH_MSG_WARNING(
                    "The panelID corresponding prd hit does NOT link to a "
                    "known Panel "
                    << m_idHelperSvc->toString(id) << " <"
                    << m_idHelperSvc->toString(
                           rpcData->detectorElement()->identify())
                    << "> !!!");
                continue;
            } else {
                panel_index = i_panel->second->panel_index;
            }
 
            sec_layer = temp_panel->getSectorLayer();
            prd_sec_all = sec_layer.first;
            prd_layer_all = sec_layer.second;
 
            if (std::abs(sec_layer.first) == 12 ||
                std::abs(sec_layer.first) == 14) {
                prd_sec_1214 = sec_layer.first;
                prd_layer_1214 = sec_layer.second;
            }
 
            fill(tool, prd_sec_all, prd_layer_all, prd_sec_1214,
                 prd_layer_1214);
 
            if (measphi == 0) {
                prd_sec_all_eta = sec_layer.first;
                prd_layer_all_eta = sec_layer.second;
                fill(tool, prd_sec_all_eta, prd_layer_all_eta);
            } else {
                prd_sec_all_phi = sec_layer.first;
                prd_layer_all_phi = sec_layer.second;
                fill(tool, prd_sec_all_phi, prd_layer_all_phi);
            }
 
            i_panelIndex = panel_index;
            fill(tool, i_prd_LB, i_panelIndex);
 
            v_prdTime.push_back(rpcData->time());
        }  // loop on RpcPrepData
    }      // loop on RpcPrepData container
 
    auto prdTimeCollection = Collection("prdTime", v_prdTime);
    fill(tool, prdTimeCollection);
 
    ATH_MSG_DEBUG(" fillHistPRD finished ");
    return StatusCode::SUCCESS;
}

fillMuonExtrapolateEff()

StatusCode RpcTrackAnaAlg::fillMuonExtrapolateEff ( const EventContext & ctx ) const

private

Definition at line 326 of file RpcTrackAnaAlg.cxx.

                                   {
    using namespace Monitored;
    auto tool = getGroup(m_packageName);
 
    //
    // read PrimaryVertex Z
    //
    const xAOD::Vertex* primVertex = nullptr;
    if (!m_PrimaryVertexContainerKey.empty()) {
        SG::ReadHandle<xAOD::VertexContainer> primVtxContainer(
            m_PrimaryVertexContainerKey, ctx);
        if (primVtxContainer.isValid()) {
            for (const auto vtx : *primVtxContainer) {
                if (vtx->vertexType() == xAOD::VxType::VertexType::PriVtx) {
                    primVertex = vtx;
                    break;
                }
            }
        }
    }
    const double primaryVertexZ = primVertex ? (primVertex->z()) : (-999);
 
    //
    // read beam position sigma from eventinfo
    //
    SG::ReadDecorHandle<xAOD::EventInfo, float> beamSigmaX(m_beamSigmaX, ctx);
    SG::ReadDecorHandle<xAOD::EventInfo, float> beamSigmaY(m_beamSigmaY, ctx);
    SG::ReadDecorHandle<xAOD::EventInfo, float> beamSigmaXY(m_beamSigmaXY, ctx);
 
    const float beamPosSigmaX = beamSigmaX(0);
    const float beamPosSigmaY = beamSigmaY(0);
    const float beamPosSigmaXY = beamSigmaXY(0);
 
    //
    // read muon
    //
    SG::ReadHandle<xAOD::MuonContainer> muons(m_MuonContainerKey, ctx);
    ATH_MSG_DEBUG(" muons size = " << muons->size());
 
    if (!muons.isValid()) {
        ATH_MSG_ERROR("evtStore() does not contain muon Collection with name "
                      << m_MuonContainerKey);
        return StatusCode::FAILURE;
    }
 
    //
    // select out tag and probe muons
    //
    std::vector<std::shared_ptr<MyMuon>> tagmuons;
    std::vector<std::shared_ptr<MyMuon>> probemuons;
    for (const xAOD::Muon* muon : *muons) {
        if (std::abs(muon->eta()) > m_maxEta)
            continue;
 
        auto mymuon = std::make_shared<MyMuon>();
        mymuon->muon = muon;
        mymuon->fourvec.SetPtEtaPhiM(muon->pt(), muon->eta(), muon->phi(),
                                     m_muonMass.value());
        mymuon->tagged =
            triggerMatching(muon, m_trigTagDefs) == StatusCode::SUCCESS;
 
        // muon quality
        if (muon->quality() > xAOD::Muon::Medium)
            continue;
 
        //
        // calculate muon z0sin(\theta) and d0significance
        auto track = muon->primaryTrackParticle();
        const double z0 = track->z0() + track->vz() - primaryVertexZ;
 
        double z0sin = z0 * std::sin(track->theta());
        double d0sig = xAOD::TrackingHelpers::d0significance(
            track, beamPosSigmaX, beamPosSigmaY, beamPosSigmaXY);
 
        // select probe muons
        if (std::abs(z0sin) < 0.5 && std::abs(d0sig) < 3.0) {
            mymuon->tagProbeOK = true;
        }
        probemuons.push_back(mymuon);
 
        // select tag muons
        if (muon->pt() > 27.0e3 && std::abs(z0sin) < 1.0 &&
            std::abs(d0sig) < 5.0) {
            tagmuons.push_back(mymuon);
        }
    }
 
    //
    // Z tag & probe
    //
    for (const auto& tag_muon : tagmuons) {
        if (!(tag_muon->tagged))
            continue;
 
        for (const auto& probe_muon : probemuons) {
            if (tag_muon->muon == probe_muon->muon)
                continue;
 
            // probe muon
            if (!probe_muon->tagProbeOK)
                continue;
 
            // Opposite charge
            if (tag_muon->muon->charge() == probe_muon->muon->charge())
                continue;
 
            // Z mass window
            float dimuon_mass = (tag_muon->fourvec + probe_muon->fourvec).M();
            if (dimuon_mass < m_zMass_lowLimit || dimuon_mass > m_zMass_upLimit)
                continue;
 
            // angular separation
            float dr = (tag_muon->fourvec).DeltaR(probe_muon->fourvec);
            if (dr < m_isolationWindow)
                continue;
 
            probe_muon->tagProbeAndZmumuOK = true;
        }
    }
 
    //
    // read rois: raw LVL1MuonRoIs
    //
    std::vector<const xAOD::MuonRoI*> roisBarrel;
    if (!m_MuonRoIContainerKey.empty()) {
        SG::ReadHandle<xAOD::MuonRoIContainer> muonRoIs(m_MuonRoIContainerKey,
                                                        ctx);
 
        if (!muonRoIs.isValid()) {
            ATH_MSG_ERROR(
                "evtStore() does not contain muon L1 ROI Collection with name "
                << m_MuonRoIContainerKey);
            return StatusCode::FAILURE;
        }
        const xAOD::MuonRoIContainer* rois = muonRoIs.cptr();
 
        std::vector<float> roiEtaVec{}, roiBarrelEtaVec{};
        std::vector<int> roiBarrelThrVec{};
 
        roiEtaVec.reserve(muonRoIs->size());
        roiBarrelEtaVec.reserve(muonRoIs->size());
        roiBarrelThrVec.reserve(muonRoIs->size());
        roisBarrel.reserve(muonRoIs->size());
        for (const xAOD::MuonRoI* roi : *rois) {
            roiEtaVec.push_back(roi->eta());
            if (roi->getSource() != xAOD::MuonRoI::RoISource::Barrel)
                continue;
 
            roiBarrelEtaVec.push_back(roi->eta());
            roiBarrelThrVec.push_back(roi->getThrNumber());
            roisBarrel.push_back(roi);
        }
 
        auto roiEtaCollection = Collection("roiEta", roiEtaVec);
        auto roiBarrelEtaCollection =
            Collection("roiBarrelEta", roiBarrelEtaVec);
        auto roiBarrelThrCollection =
            Collection("roiBarrelThr", roiBarrelThrVec);
        fill(tool, roiEtaCollection);
        fill(tool, roiBarrelEtaCollection);
        fill(tool, roiBarrelThrCollection);
    }
 
    //
    // Fill Lv1 trigger efficiency
    // Fill muon detection efficiency
    auto i_pt_allMu = Scalar<double>("muPt_allMu", 0.);
    auto i_eta_allMu = Scalar<double>("muEta_allMu", 0.);
    auto i_phi_allMu = Scalar<double>("muPhi_allMu", 0.);
    auto i_pt_zMu = Scalar<double>("muPt_MuonFromZ", 0.);
    auto i_eta_zMu = Scalar<double>("muEta_MuonFromZ", 0.);
    auto i_phi_zMu = Scalar<double>("muPhi_MuonFromZ", 0.);
    auto i_eta_zMu_p =
        Scalar<double>("muEta_p_MuonFromZ",
                       0.);  // kine variable at the plateau of pt turn-on curve
    auto i_phi_zMu_p =
        Scalar<double>("muPhi_p_MuonFromZ",
                       0.);  // kine variable at the plateau of pt turn-on curve
 
    std::vector<GasGapResult> results;
    int nmuon{0}, nmuon_barrel{0};
    for (const auto& probe_muon : probemuons) {
        nmuon++;
        double pt = probe_muon->muon->pt();
        double eta = probe_muon->muon->eta();
        double phi = probe_muon->muon->phi();
 
        // barrel muon
        if (std::abs(eta) < m_barrelMinEta || std::abs(eta) > m_barrelMaxEta)
            continue;
        nmuon_barrel++;
 
        // has MuonSpectrometerTrackParticle
        const xAOD::TrackParticle* track = probe_muon->muon->trackParticle(
            xAOD::Muon::MuonSpectrometerTrackParticle);
        if (track) {
            results.clear();
 
            ATH_CHECK(extrapolate2RPC(track, Trk::anyDirection, results, BI));
            ATH_CHECK(readHitsPerGasgap(ctx, results, AllMuon));
        }
 
        // fill kine variables for probe muons
        i_pt_allMu = pt;
        i_eta_allMu = eta;
        i_phi_allMu = phi;
        fill(tool, i_pt_allMu, i_eta_allMu, i_phi_allMu);
 
        //
        // Probe muons with Z tag&probe method
        // 1) Plot L1 trigger efficiency use probe muons
        // 2) Plot muon detection efficiency use probe muons
        if (probe_muon->tagProbeAndZmumuOK) {
            //
            // do match muon with ROI
            std::vector<bool> isMatcheds(6, 0);
            for (const xAOD::MuonRoI* roi : roisBarrel) {
                const double dphi = TVector2::Phi_mpi_pi(roi->phi() - phi);
                const double deta = roi->eta() - eta;
                const double dr = std::hypot(dphi, deta);
 
                if (dr > m_l1trigMatchWindow)
                    continue;
 
                int thr = std::min(7, roi->getThrNumber());
                for (int i_thr = 1; i_thr <= thr; i_thr++) {
                    isMatcheds[i_thr - 1] = true;
                }
            }
 
            //
            // fill L1 trigger Efficiency
            auto i_pt = Scalar<double>("muPt_l1", pt);
            auto i_eta = Scalar<double>("muEta_l1", eta);
            auto i_phi = Scalar<double>("muPhi_l1", phi);
            auto i_passTrigger = Scalar<bool>("passTrigger", false);
            auto i_passTrigger_1 = Scalar<bool>("passTrigger_plateau", false);
 
            for (int i_thr = 1; i_thr < 7; i_thr++) {
                i_passTrigger = isMatcheds[i_thr - 1];
                i_passTrigger_1 = isMatcheds[i_thr - 1];
 
                // plot L1 trigger pt turn-on curve
                fill(m_tools[m_TriggerThrGroup.at("thr" +
                                                  std::to_string(i_thr))],
                     i_pt, i_passTrigger);
 
                // plot L1 trigger efficiency on the plateau of pt turn-on curve
                if (std::abs(pt) > m_minPt) {
                    fill(m_tools[m_TriggerThrGroup.at("thr" +
                                                      std::to_string(i_thr))],
                         i_eta, i_phi, i_passTrigger_1);
                }
            }
 
            //
            // muon track
            // fill muon detection Efficiency
            //
            if (track) {
                ATH_CHECK(readHitsPerGasgap(ctx, results, ZCand));
            }
 
            // fill kine variables for probe muons in ZTP
            i_pt_zMu = pt;
            i_eta_zMu = eta;
            i_phi_zMu = phi;
            fill(tool, i_pt_zMu, i_eta_zMu, i_phi_zMu);
 
            if (std::abs(pt) > m_minPt) {
                i_eta_zMu_p = eta;
                i_phi_zMu_p = phi;
                fill(tool, i_eta_zMu_p, i_phi_zMu_p);
            }
        }  // tagProbeAndZmumuOK
    }      // probemuons
 
    auto Nmuon = Scalar<int>("nMu", nmuon);
    auto Nmuon_barrel = Scalar<int>("nMuBarrel", nmuon_barrel);
 
    //
    // Fill histograms
    //
    fill(tool, Nmuon, Nmuon_barrel);
 
    return StatusCode::SUCCESS;
}

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

getTrigDecisionTool()

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

inherited

Get the trigger decision tool member.

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

Returns

m_trigDecTool

Definition at line 198 of file AthMonitorAlgorithm.cxx.

                                                                                       {
    return m_trigDecTool;
}

initArrayHistosMap()

StatusCode RpcTrackAnaAlg::initArrayHistosMap ( )

private

initialize()

StatusCode RpcTrackAnaAlg::initialize ( )

overridevirtual

initialize

Returns

StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 38 of file RpcTrackAnaAlg.cxx.

                                      {
    ATH_MSG_INFO(" RpcTrackAnaAlg initialize begin ");
    ATH_CHECK(AthMonitorAlgorithm::initialize());
 
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(m_beamSigmaX.initialize());
    ATH_CHECK(m_beamSigmaY.initialize());
    ATH_CHECK(m_beamSigmaXY.initialize());
 
    ATH_CHECK(m_MuonRoIContainerKey.initialize(SG::AllowEmpty));
    ATH_CHECK(m_MuonContainerKey.initialize());
    ATH_CHECK(m_rpcPrdKey.initialize());
    ATH_CHECK(m_PrimaryVertexContainerKey.initialize(SG::AllowEmpty));
 
    ATH_CHECK(readElIndexFromXML());
    ATH_CHECK(initRpcPanel());
 
    ATH_CHECK(initTrigTag());
 
    std::vector<std::string> sectorStr = {
        "sector1",  "sector2",  "sector3",  "sector4",  "sector5",  "sector6",
        "sector7",  "sector8",  "sector9",  "sector10", "sector11", "sector12",
        "sector13", "sector14", "sector15", "sector16"};
    m_SectorGroup =
        Monitored::buildToolMap<int>(m_tools, "RpcTrackAnaAlg", sectorStr);
 
    std::vector<std::string> triggerThrs = {"thr1", "thr2", "thr3",
                                            "thr4", "thr5", "thr6"};
    m_TriggerThrGroup =
        Monitored::buildToolMap<int>(m_tools, "RpcTrackAnaAlg", triggerThrs);
 
    ATH_MSG_INFO(" initialize extrapolator ");
    ATH_CHECK(m_extrapolator.retrieve());
 
    ATH_MSG_INFO(" RpcTrackAnaAlg initialize END ");
    return StatusCode::SUCCESS;
}

initRpcPanel()

StatusCode RpcTrackAnaAlg::initRpcPanel ( )

private

Definition at line 77 of file RpcTrackAnaAlg.cxx.

                                        {
    /*
      Iterate over all RpcDetectorElements and RpcReadoutElements
      and cache locally all panel
    */
    ATH_MSG_INFO(name() << " - RpcTrackAnaAlg::initRpcPanel - start");
 
    unsigned int nValidPanel = 0, nTriedPanel{0};
    const RpcIdHelper& rpcIdHelper = m_idHelperSvc->rpcIdHelper();
    const MuonGM::MuonDetectorManager* muonMgr{nullptr};
    ATH_CHECK(detStore()->retrieve(muonMgr));
 
    ATH_MSG_INFO("MuonGM::MuonDetectorManager::RpcDetElMaxHash= "
                 << rpcIdHelper.module_hash_max());
 
    m_StationNames[BI] = {};
    m_StationNames[BM1] = {2, 3, 8, 53};  // doubletR = 1
    m_StationNames[BM2] = {2, 3, 8, 53};  // doubletR = 2
    m_StationNames[BO1] = {4, 5, 9, 10};  // doubletR = 1
    m_StationNames[BO2] = {4, 5, 9, 10};  // doubletR = 2
 
    std::vector<int> BMBO_StationNames = {2, 3, 4, 5, 8, 9, 10, 53};
    for (auto idetEl = rpcIdHelper.detectorElement_begin();
              idetEl != rpcIdHelper.detectorElement_end(); ++idetEl) {
        const MuonGM::RpcReadoutElement* readoutEl = muonMgr->getRpcReadoutElement(*idetEl);
        if (!readoutEl)
            continue;
 
        const unsigned int ngasgap = 2;
        const unsigned int doubletZ = readoutEl->getDoubletZ();
        if (readoutEl->getDoubletPhi() != 1) {
            continue;
        }
 
        const Identifier readEl_id = readoutEl->identify();
 
        int stName = rpcIdHelper.stationName(readEl_id);
 
        if (!std::count(BMBO_StationNames.begin(), BMBO_StationNames.end(),
                        stName)) {
            continue;  // Will be changed to include BIS
        }
 
        if (!std::count(BMBO_StationNames.begin(), BMBO_StationNames.end(),
                        stName)) {
            continue;  // Will be changed to include BIS
        }
 
        for (unsigned doubletPhi = 1; doubletPhi <= 2; ++doubletPhi) {
            for (unsigned gasgap = 1; gasgap <= ngasgap; ++gasgap) {
                std::shared_ptr<GasGapData> gap = std::make_shared<GasGapData>(
                    *m_idHelperSvc, readoutEl, doubletZ, doubletPhi, gasgap);
 
                std::pair<int, int> st_dbR = std::make_pair(stName, gap->doubletR);
                m_gasGapData[st_dbR].push_back(gap);
 
                std::shared_ptr<RpcPanel> rpcPanel_eta = std::make_shared<RpcPanel>(
                    *m_idHelperSvc, readoutEl, doubletZ, doubletPhi, gasgap, 0);
                ATH_CHECK(setPanelIndex(rpcPanel_eta));
 
                std::shared_ptr<RpcPanel> rpcPanel_phi = std::make_shared<RpcPanel>(
                    *m_idHelperSvc, readoutEl, doubletZ, doubletPhi, gasgap, 1);
                ATH_CHECK(setPanelIndex(rpcPanel_phi));
 
                if (rpcPanel_eta->panel_valid) {
                    ATH_MSG_DEBUG(" Panel  stationName:"
                                  << rpcPanel_eta->stationName
                                  << " stationEta:" << rpcPanel_eta->stationEta
                                  << " stationPhi:" << rpcPanel_eta->stationPhi
                                  << " doubletR:" << rpcPanel_eta->doubletR
                                  << "doubletZ:" << rpcPanel_eta->doubletZ
                                  << " doubletPhi:" << rpcPanel_eta->doubletPhi
                                  << " gasGap:" << rpcPanel_eta->gasGap
                                  << " measPhi:" << rpcPanel_eta->measPhi);
 
                    m_rpcPanelMap.insert(
                        std::map<Identifier, std::shared_ptr<RpcPanel>>::value_type(
                            rpcPanel_eta->panelId, rpcPanel_eta));
                    nValidPanel++;
                }
 
                if (rpcPanel_phi->panel_valid) {
                    ATH_MSG_DEBUG(" Panel  stationName:"
                                  << rpcPanel_phi->stationName
                                  << " stationEta:" << rpcPanel_phi->stationEta
                                  << " stationPhi:" << rpcPanel_phi->stationPhi
                                  << " doubletR:" << rpcPanel_phi->doubletR
                                  << "doubletZ:" << rpcPanel_phi->doubletZ
                                  << " doubletPhi:" << rpcPanel_phi->doubletPhi
                                  << " gasGap:" << rpcPanel_phi->gasGap
                                  << " measPhi:" << rpcPanel_phi->measPhi);
 
                    m_rpcPanelMap.insert(
                        std::map<Identifier, std::shared_ptr<RpcPanel>>::value_type(
                            rpcPanel_phi->panelId, rpcPanel_phi));
                    nValidPanel++;
                }
 
                gap->RpcPanel_eta_phi = std::make_pair(rpcPanel_eta, rpcPanel_phi);
                ++nTriedPanel;
            }
        }
    }
 
    ATH_MSG_INFO("Number of valid panels = " << nValidPanel << " tried panels "
                                             << nTriedPanel);
    return StatusCode::SUCCESS;
}

initTrigTag()

StatusCode RpcTrackAnaAlg::initTrigTag ( )

private

Definition at line 280 of file RpcTrackAnaAlg.cxx.

                                       {
    TString trigStr = m_trigTagList.value();
    std::unique_ptr<TObjArray> tagList(trigStr.Tokenize(";"));
 
    // TObjArray* tagList = TString(m_trigTagList.value()).Tokenize(",");
    std::set<TString> alllist;
    for (int i = 0; i < tagList->GetEntries(); i++) {
        TString tagTrig = tagList->At(i)->GetName();
        if (alllist.find(tagTrig) != alllist.end())
            continue;
        alllist.insert(tagTrig);
        std::unique_ptr<TObjArray> arr(tagTrig.Tokenize(";"));
        if (arr->GetEntries() == 0)
            continue;
        TagDef def;
        def.eventTrig = TString(arr->At(0)->GetName());
        def.tagTrig = def.eventTrig;
        if (arr->GetEntries() == 2)
            def.tagTrig = TString(arr->At(1)->GetName());
        m_trigTagDefs.push_back(def);
    }
 
    return StatusCode::SUCCESS;
}

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

isClonable()

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

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

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

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

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

IsNearbyHit()

bool RpcTrackAnaAlg::IsNearbyHit ( const std::vector< const Muon::RpcPrepData * > & cluster_hits, const Muon::RpcPrepData * hit ) const

private

Definition at line 1412 of file RpcTrackAnaAlg.cxx.

                                      {
    const RpcIdHelper& rpcIdHelper = m_idHelperSvc->rpcIdHelper();
 
    // Check whether this hit is close to any hits in the cluster
    for (const Muon::RpcPrepData* it_hit : cluster_hits) {
        if (abs(rpcIdHelper.strip(it_hit->identify()) -
                rpcIdHelper.strip(hit->identify())) < 2 &&
            std::abs(it_hit->time() - hit->time()) < 6.5) {
            return true;
        }
    }
 
    return false;
}

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

parseList()

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

virtualinherited

Parse a string into a vector.

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

Parameters

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

Returns

StatusCode

Definition at line 345 of file AthMonitorAlgorithm.cxx.

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

readElIndexFromXML()

StatusCode RpcTrackAnaAlg::readElIndexFromXML ( )

private

Definition at line 216 of file RpcTrackAnaAlg.cxx.

                                              {
    /*
      Read xml file in share/Element.xml to give each element(in the afterwards,
      for panel) a index
    */
    ATH_MSG_INFO(name() << " - read element index ");
 
    //
    // Read xml file
    //
    const std::string xml_file =
        PathResolver::find_file(m_elementsFileName, "DATAPATH");
    if (xml_file.empty()) {
        ATH_MSG_ERROR(m_elementsFileName << " not found!");
        return StatusCode::FAILURE;
    }
 
    ATH_MSG_INFO(name() << " - read xml file: " << xml_file);
 
    boost::property_tree::ptree pt;
    read_xml(xml_file, pt);
 
    for (boost::property_tree::ptree::value_type& child:
                   pt.get_child("Elements")) {
        if (child.first == "Element") {
            int index = child.second.get<int>("<xmlattr>.index");
            int stName = child.second.get<int>("<xmlattr>.stationName");
            int stEta = child.second.get<int>("<xmlattr>.stationEta");
            int stPhi = child.second.get<int>("<xmlattr>.stationPhi");
            int dbR = child.second.get<int>("<xmlattr>.doubletR");
            int dbZ = child.second.get<int>("<xmlattr>.doubletZ");
 
            ATH_MSG_DEBUG(" index = " << index << ", stationName = " << stName
                                      << ", stationEta = " << stEta
                                      << ", stationPhi = " << stPhi
                                      << ", doubletR = " << dbR
                                      << ", doubletZ = " << dbZ);
            std::ostringstream ele_key;
 
            ele_key << stName << "_" << stEta << "_" << stPhi << "_" << dbR
                    << "_" << dbZ;
 
            std::map<std::string, int>::iterator it_eleInd =
                m_elementIndex.find(ele_key.str());
            if (it_eleInd == m_elementIndex.end()) {
                m_elementIndex.insert(std::map<std::string, int>::value_type(
                    ele_key.str(), index));
            } else {
                ATH_MSG_ERROR("These is duplicated elements");
            }
        } else {
            ATH_MSG_ERROR(" node key != Element");
        }
    }
 
    ATH_MSG_INFO(
        "Total number of elements: "
        << m_elementIndex.size()
        << ", should be consistent with the number of elements in xml file!");
 
    return StatusCode::SUCCESS;
}

readHitsPerGasgap()

StatusCode RpcTrackAnaAlg::readHitsPerGasgap ( const EventContext & ctx, std::vector< GasGapResult > & results, MuonSource muon_source ) const

private

Definition at line 1084 of file RpcTrackAnaAlg.cxx.

                                                                           {
    using namespace Monitored;
    auto tool = getGroup(m_packageName);
 
    int lumiBlock = ctx.eventID().lumi_block();
 
    SG::ReadHandle<Muon::RpcPrepDataContainer> rpcContainer(m_rpcPrdKey, ctx);
    const RpcIdHelper& rpcIdHelper = m_idHelperSvc->rpcIdHelper();
 
    ATH_MSG_DEBUG(" RpcPrepDataContainer size = " << rpcContainer->size());
    ATH_MSG_DEBUG(" results size = " << results.size());
 
    auto i_hitTime_sec = Scalar<int>("hitTime_sec", 0);
 
    auto isOutTime_prd = Scalar<bool>("isOutTime_prd", false);
    auto isOutTime_onTrack = Scalar<bool>("isOutTime_prd_onTrack", false);
    auto i_panelIndex = Scalar<int>("panelInd_prd", -1);
    auto i_panelIndex_onTrack = Scalar<int>("panelInd_prd_onTrack", -1);
 
    auto res_eta = Scalar<int>("residual_eta", 0);
    auto res_phi = Scalar<int>("residual_phi", 0);
    auto closest_res_eta = Scalar<int>("closest_residual_eta", 0);
    auto closest_res_phi = Scalar<int>("closest_residual_phi", 0);
 
    auto res_panel = Scalar<int>("residual_panel", 0);
    auto i_panelInd_res = Scalar<int>("panelInd_res_inTime", -1);
 
    for (GasGapResult& exr : results) {
        const std::shared_ptr<GasGapData> gap = exr.second;
 
        int sector = (gap->RpcPanel_eta_phi.first->getSectorLayer()).first;
        
        float clo_res_eta = 1001.; // initial value; no sense
        float clo_res_phi = 1001.;
        int NHitwithCut_perMuon_eta = 0;
        int NHitwithCut_perMuon_phi = 0;
        int NHitnoCut_perMuon_eta = 0;
        int NHitnoCut_perMuon_phi = 0;
        std::vector<const Muon::RpcPrepData*> view_hits_eta;
        std::vector<const Muon::RpcPrepData*> view_hits_phi;
 
        // loop on RpcPrepDataCollection
        for (const Muon::RpcPrepDataCollection* rpcCollection : *rpcContainer) {
            if (!rpcCollection) {
                continue;
            }
 
            // loop on RpcPrepData
            for (const Muon::RpcPrepData* rpcData : *rpcCollection) {
                if (!rpcData) {
                    continue;
                }
 
                const Identifier id = rpcData->identify();
                const int stationName = rpcIdHelper.stationName(id);
                const int stationEta = rpcIdHelper.stationEta(id);
                const int stationPhi = rpcIdHelper.stationPhi(id);
 
                const int doubletR = rpcIdHelper.doubletR(id);
                const int doubletZ = rpcIdHelper.doubletZ(id);
                const int doubletPhi = rpcIdHelper.doubletPhi(id);
                const unsigned gasGap = rpcIdHelper.gasGap(id);
                const int measuresPhi = rpcIdHelper.measuresPhi(id);
 
                // match hit to the gasgap
                if (stationName == gap->stationName &&
                    stationPhi == gap->stationPhi &&
                    stationEta == gap->stationEta &&
                    doubletR == gap->doubletR && gasGap == gap->gasgap &&
                    doubletPhi == gap->doubletPhi &&
                    doubletZ == gap->doubletZ) {
 
 
                    i_hitTime_sec = rpcData->time();
                    fill(m_tools[m_SectorGroup.at(
                             "sector" + std::to_string(std::abs(sector)))],
                         i_hitTime_sec);
 
                    // ---------------------------------
                    // Calculate distance between extrapolated muon track position and hit
                    float hit_local_x = rpcData->localPosition().x();
                    float trackPos_localY = exr.first.localPos.y();
                    float trackPos_localZ = exr.first.localPos.z();
 
                    float residual_phi = trackPos_localY-hit_local_x;
                    float residual_eta = trackPos_localZ-hit_local_x;
 
                    // If hit is out-of-time
                    bool isOutTime = (std::abs(i_hitTime_sec+50.) > m_outtime); // for run 3
                    int i_panel = measuresPhi ? gap->RpcPanel_eta_phi.second->panel_index :
                        gap->RpcPanel_eta_phi.first->panel_index;
                    
                    // Fill histograms of out-of-time hit fraction
                    if (muon_source == ZCand) {
                        
                        isOutTime_prd = isOutTime;
                        isOutTime_onTrack = isOutTime;
                        
                        if (measuresPhi) {
                            i_panelIndex = i_panel;
                            fill(tool, i_panelIndex, isOutTime_prd);
                            if (std::abs(residual_phi) < m_diffHitTrackPostion) {
                                i_panelIndex_onTrack = i_panel;
                                fill(tool, i_panelIndex_onTrack, isOutTime_onTrack);
                            }
                        } else {
                            i_panelIndex = i_panel;
                            fill(tool, i_panelIndex, isOutTime_prd);
                            if (std::abs(residual_eta) < m_diffHitTrackPostion) {
                                i_panelIndex_onTrack = i_panel;
                                fill(tool, i_panelIndex_onTrack, isOutTime_onTrack);
                            }
                        }
                    }
 
                    if (measuresPhi) {
                        NHitnoCut_perMuon_phi++;
                    } else {
                        NHitnoCut_perMuon_eta++;
                    }
 
                    // process hit within |time| < 12.5 ns && on Track
                    if (isOutTime){
                        continue;
                    }
                    i_panelInd_res = i_panel;
 
                    if (measuresPhi) {
                        res_phi = residual_phi;
                        res_panel = residual_phi;
                        clo_res_phi = std::min(clo_res_phi, residual_phi);
 
                        fill(tool, res_phi, res_panel, i_panelInd_res);
 
                        if (std::abs(residual_phi) < m_diffHitTrackPostion) {
                            NHitwithCut_perMuon_phi++;
                            view_hits_phi.push_back(rpcData);
                        }
                    } else {
                        res_eta = residual_eta;
                        res_panel = residual_eta;
                        clo_res_eta = std::min(clo_res_eta, residual_eta);
 
                        fill(tool, res_eta, res_panel, i_panelInd_res);
 
                        if (std::abs(residual_eta) < m_diffHitTrackPostion){
                            NHitwithCut_perMuon_eta++;
                            view_hits_eta.push_back(rpcData);
                        }
                    }
                }
            }
        }
 
        int etaPanel_ind = -1;
        int phiPanel_ind = -1;
 
        etaPanel_ind = gap->RpcPanel_eta_phi.first->panel_index;
        phiPanel_ind = gap->RpcPanel_eta_phi.second->panel_index;
 
        // Declare the quantities which should be monitored
        if (muon_source == ZCand) {
            auto hitMulti_eta = Scalar<int>("hitMulti_eta", NHitwithCut_perMuon_eta);
            auto hitMulti_phi = Scalar<int>("hitMulti_phi", NHitwithCut_perMuon_phi);
            auto hitMulti = Scalar<int>("hitMulti", 0);
            auto i_panelIndex = Scalar<int>("panelInd_hM", -1);
            auto i_passExtrap = Scalar<bool>("muon_passExtrap", false);
            auto i_passExtrap_or  = Scalar<bool>("muon_passExtrap_or", false);
            auto i_passExtrap_and = Scalar<bool>("muon_passExtrap_and", false);
            auto i_passExtrap_sig_or  = Scalar<bool>("muon_passExtrap_signalhit_or", false);
            auto i_passExtrap_sig_and = Scalar<bool>("muon_passExtrap_signalhit_and", false);
            auto i_passExtrap_sig     = Scalar<bool>("muon_passExtrap_signalhit", false);
            auto i_LB = Scalar<int>("LB_detEff", lumiBlock);
 
            fill(tool, hitMulti_eta, hitMulti_phi);
 
            //
            // Eta OR Phi panel
            if (NHitnoCut_perMuon_eta > 0 || NHitnoCut_perMuon_phi > 0)
                i_passExtrap_or = true;
            if (NHitwithCut_perMuon_eta > 0 || NHitwithCut_perMuon_phi > 0)
                i_passExtrap_sig_or = true;
                       
            //
            // Eta AND Phi panel
            if (NHitnoCut_perMuon_eta > 0 && NHitnoCut_perMuon_phi > 0)
                i_passExtrap_and = true;
            if (NHitwithCut_perMuon_eta > 0 && NHitwithCut_perMuon_phi > 0)
                i_passExtrap_sig_and = true;
            //
            // Eta panel
            hitMulti = NHitwithCut_perMuon_eta;
            i_panelIndex = etaPanel_ind;
 
            if (NHitnoCut_perMuon_eta > 0)
                i_passExtrap = true;
 
            if (NHitwithCut_perMuon_eta > 0)
                i_passExtrap_sig = true;
 
            fill(tool, hitMulti, i_panelIndex, i_passExtrap, i_passExtrap_sig, 
                    i_passExtrap_or, i_passExtrap_and, i_passExtrap_sig_or, i_passExtrap_sig_and, i_LB);
            ATH_CHECK(fillClusterSize(view_hits_eta, etaPanel_ind, lumiBlock,
                                      sector, 0));  // isPhi = 0
 
            if (clo_res_eta < 1000.) {
                closest_res_eta = clo_res_eta;
                fill(tool, closest_res_eta);
            }
 
            //
            // Phi panel
            hitMulti = NHitwithCut_perMuon_phi;
            i_panelIndex = phiPanel_ind;
            if (NHitnoCut_perMuon_phi > 0)
                i_passExtrap = true;
 
            if (NHitwithCut_perMuon_phi > 0)
                i_passExtrap_sig = true; 
 
            fill(tool, hitMulti, i_panelIndex, i_passExtrap, i_passExtrap_sig, 
                    i_passExtrap_or, i_passExtrap_and, i_passExtrap_sig_or, i_passExtrap_sig_and, i_LB);
            ATH_CHECK(fillClusterSize(view_hits_phi, phiPanel_ind, lumiBlock,
                                      sector, 1));  // isPhi = 1
 
            if (clo_res_phi < 1000.) {
                closest_res_phi = clo_res_phi;
                fill(tool, closest_res_phi);
            }
        }
 
        // 
        // All muon
        // 
        auto i_panelIndex_allMu = Scalar<int>("panelInd_hM_allMu", -1);
        auto i_passExtrap_allMu = Scalar<bool>("muon_passExtrap_allMu", false);
 
        //
        // Eta panel
        i_panelIndex_allMu = etaPanel_ind;
        if (NHitwithCut_perMuon_eta > 0)
            i_passExtrap_allMu = true;
        fill(tool, i_panelIndex_allMu, i_passExtrap_allMu);
 
        //
        // Phi panel
        i_panelIndex_allMu = phiPanel_ind;
        if (NHitwithCut_perMuon_phi > 0)
            i_passExtrap_allMu = true;
 
        fill(tool, i_panelIndex_allMu, i_passExtrap_allMu);
    }
 
    return StatusCode::SUCCESS;
}

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

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

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

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

setPanelIndex()

StatusCode RpcTrackAnaAlg::setPanelIndex ( std::shared_ptr< RpcPanel > panel )

private

Definition at line 187 of file RpcTrackAnaAlg.cxx.

                                                                      {
    /*
      Specify every panel with a unique index according the information in xml
    */
    std::string ele_str = panel->getElementStr();
 
    std::map<std::string, int>::iterator it_eleInd =
        m_elementIndex.find(ele_str);
 
    if (it_eleInd == m_elementIndex.end()) {
        ATH_MSG_ERROR("Can not find element: stationName = "
                      << panel->stationName
                      << ", stationEta = " << panel->stationEta
                      << ", stationPhi = " << panel->stationPhi
                      << ", doubletR = " << panel->doubletR
                      << ", doubletZ = " << panel->doubletZ);
        return StatusCode::FAILURE;
    }
 
    int ele_index = it_eleInd->second;
    int panel_index = (ele_index - 1) * 8 + (panel->doubletPhi - 1) * 4 +
                      (panel->gasGap - 1) * 2 + panel->measPhi;
 
    panel->SetPanelIndex(panel_index);
 
    return StatusCode::SUCCESS;
}

sysExecute()

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

overridevirtualinherited

Execute an algorithm.

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

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

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

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

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

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

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

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

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

sysStart()

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

overridevirtualinherited

Handle START transition.

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

trigChainsArePassed()

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

inherited

Check whether triggers are passed.

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

Parameters

vTrigNames

List of trigger names.

Returns

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

Definition at line 203 of file AthMonitorAlgorithm.cxx.

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

triggerMatching()

StatusCode RpcTrackAnaAlg::triggerMatching ( const xAOD::Muon * offline_muon, const std::vector< TagDef > & list_of_triggers ) const

private

Definition at line 716 of file RpcTrackAnaAlg.cxx.

                                                     {
    if (!m_TagAndProbe.value())
        return StatusCode::SUCCESS;
    if (getTrigDecisionTool().empty())
        return StatusCode::SUCCESS;
    TVector3 muonvec;
    muonvec.SetPtEtaPhi(offline_muon->pt(), offline_muon->eta(),
                        offline_muon->phi());
 
    for (const auto& tagTrig : list_of_triggers) {
        if (!getTrigDecisionTool()->isPassed(tagTrig.eventTrig.Data()))
            continue;
 
        ATH_MSG_DEBUG("tagTrig.eventTrig = " << tagTrig.eventTrig
                                             << ";  tagTrig.tagTrig = "
                                             << tagTrig.tagTrig);
        // ATH_MSG_DEBUG("m_MuonEFContainerName.value() = "<<
        // m_MuonEFContainerName.value());
 
        std::vector<TrigCompositeUtils::LinkInfo<xAOD::MuonContainer>>
            features = getTrigDecisionTool()->features<xAOD::MuonContainer>(
                tagTrig.tagTrig.Data(), TrigDefs::Physics);
 
        for (const auto& aaa : features) {
            ATH_CHECK(aaa.isValid());
            auto trigmuon_link = aaa.link;
            auto trigmuon = *trigmuon_link;
            TVector3 trigvec;
            trigvec.SetPtEtaPhi(trigmuon->pt(), trigmuon->eta(),
                                trigmuon->phi());
            if (trigvec.DeltaR(muonvec) < m_trigMatchWindow.value())
                return StatusCode::SUCCESS;
        }
    }
    return StatusCode::FAILURE;
}

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_barrelMaxEta

FloatProperty RpcTrackAnaAlg::m_barrelMaxEta {this, "barrelMaxEta", 1.05}

private

Definition at line 136 of file RpcTrackAnaAlg.h.

{this, "barrelMaxEta", 1.05};

m_barrelMinEta

FloatProperty RpcTrackAnaAlg::m_barrelMinEta {this, "barrelMinEta", 0.1}

private

Definition at line 135 of file RpcTrackAnaAlg.h.

{this, "barrelMinEta", 0.1};

m_barrelMinPt

FloatProperty RpcTrackAnaAlg::m_barrelMinPt {this, "barrelMinPt", 2.0e3}

private

Definition at line 134 of file RpcTrackAnaAlg.h.

{this, "barrelMinPt", 2.0e3};

m_beamSigmaX

SG::ReadDecorHandleKey<xAOD::EventInfo> RpcTrackAnaAlg::m_beamSigmaX

private

Initial value:

{
        this, "beamPosSigmaX", "EventInfo.beamPosSigmaX",
        "Beam spot position sigma in X"}

Definition at line 193 of file RpcTrackAnaAlg.h.

                                                    {
        this, "beamPosSigmaX", "EventInfo.beamPosSigmaX",
        "Beam spot position sigma in X"};

m_beamSigmaXY

SG::ReadDecorHandleKey<xAOD::EventInfo> RpcTrackAnaAlg::m_beamSigmaXY

private

Initial value:

{
        this, "beamPosSigmaXY", "EventInfo.beamPosSigmaXY",
        "Beam spot covariance in XY"}

Definition at line 201 of file RpcTrackAnaAlg.h.

                                                     {
        this, "beamPosSigmaXY", "EventInfo.beamPosSigmaXY",
        "Beam spot covariance in XY"};

m_beamSigmaY

SG::ReadDecorHandleKey<xAOD::EventInfo> RpcTrackAnaAlg::m_beamSigmaY

private

Initial value:

{
        this, "beamPosSigmaY", "EventInfo.beamPosSigmaY",
        "Beam spot position sigma in Y"}

Definition at line 196 of file RpcTrackAnaAlg.h.

                                                    {
        this, "beamPosSigmaY", "EventInfo.beamPosSigmaY",
        "Beam spot position sigma in Y"};

m_boundsToleranceReadoutElement

FloatProperty RpcTrackAnaAlg::m_boundsToleranceReadoutElement

private

Initial value:

{
        this, "boundsToleranceReadoutElement", 100.0,
        "boundsToleranceReadoutElement"}

Definition at line 160 of file RpcTrackAnaAlg.h.

                                                 {
        this, "boundsToleranceReadoutElement", 100.0,
        "boundsToleranceReadoutElement"};

m_boundsToleranceReadoutElementTight

FloatProperty RpcTrackAnaAlg::m_boundsToleranceReadoutElementTight

private

Initial value:

{
        this, "boundsToleranceReadoutElementTight", 20.0,
        "boundsToleranceReadoutElementTight"}

Definition at line 163 of file RpcTrackAnaAlg.h.

                                                      {
        this, "boundsToleranceReadoutElementTight", 20.0,
        "boundsToleranceReadoutElementTight"};

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_diffHitTrackPostion

FloatProperty RpcTrackAnaAlg::m_diffHitTrackPostion

private

Initial value:

{
        this, "diffHitTrackPostion", 30.0,
        "the largest distance between hit and tracks local postion"}

Definition at line 167 of file RpcTrackAnaAlg.h.

                                       {
        this, "diffHitTrackPostion", 30.0,
        "the largest distance between hit and tracks local postion"};

m_DQFilterTools

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

protectedinherited

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

{this,"FilterTools",{}}; 

m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy

privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

m_elementIndex

std::map<std::string, int> RpcTrackAnaAlg::m_elementIndex {}

private

Definition at line 214 of file RpcTrackAnaAlg.h.

{};

m_elementsFileName

StringProperty RpcTrackAnaAlg::m_elementsFileName

private

Initial value:

{this, "ElementsFileName", "RpcRawDataMonitoring/Element.xml",
                                      "Elements xml file"}

Definition at line 114 of file RpcTrackAnaAlg.h.

                                     {this, "ElementsFileName", "RpcRawDataMonitoring/Element.xml",
                                      "Elements xml file"};

m_enforceExpressTriggers

Gaudi::Property<bool> AthMonitorAlgorithm::m_enforceExpressTriggers

privateinherited

Initial value:

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

Definition at line 377 of file AthMonitorAlgorithm.h.

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

m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment

protectedinherited

Instance of the Environment_t enum.

Definition at line 355 of file AthMonitorAlgorithm.h.

m_environmentStr

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

protectedinherited

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

Definition at line 357 of file AthMonitorAlgorithm.h.

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

m_EventInfoKey

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

protectedinherited

Key for retrieving EventInfo from StoreGate.

Definition at line 367 of file AthMonitorAlgorithm.h.

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

m_evtStore

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

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_extrapolator

ToolHandle<Trk::IExtrapolator> RpcTrackAnaAlg::m_extrapolator

private

Initial value:

{
        this, "TrackExtrapolator", "Trk::Extrapolator/AtlasExtrapolator",
        "Track extrapolator"}

Definition at line 178 of file RpcTrackAnaAlg.h.

                                               {
        this, "TrackExtrapolator", "Trk::Extrapolator/AtlasExtrapolator",
        "Track extrapolator"};

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_gasGapData

std::map<std::pair<int, int>, std::vector<std::shared_ptr<GasGapData> > > RpcTrackAnaAlg::m_gasGapData {}

private

Definition at line 210 of file RpcTrackAnaAlg.h.

{};

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> RpcTrackAnaAlg::m_idHelperSvc

private

Initial value:

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

Definition at line 173 of file RpcTrackAnaAlg.h.

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

m_isolationWindow

FloatProperty RpcTrackAnaAlg::m_isolationWindow

private

Initial value:

{
        this, "IsolationWindow", 0.1,
        "Window size in R for isolation with other muons"}

Definition at line 147 of file RpcTrackAnaAlg.h.

                                   {
        this, "IsolationWindow", 0.1,
        "Window size in R for isolation with other muons"};

m_l1trigMatchWindow

FloatProperty RpcTrackAnaAlg::m_l1trigMatchWindow

private

Initial value:

{
        this, "L1TrigMatchingWindow", 0.3,
        "Window size in R for L1 trigger matching"}

Definition at line 150 of file RpcTrackAnaAlg.h.

                                     {
        this, "L1TrigMatchingWindow", 0.3,
        "Window size in R for L1 trigger matching"};

m_lbDuraThr

FloatProperty RpcTrackAnaAlg::m_lbDuraThr

private

Initial value:

{this, "lbDuraThr", 10.,
                              "Thrshold of luminosity block deruation"}

Definition at line 109 of file RpcTrackAnaAlg.h.

                             {this, "lbDuraThr", 10.,
                              "Thrshold of luminosity block deruation"};

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_maxEta

FloatProperty RpcTrackAnaAlg::m_maxEta

private

Initial value:

{this, "maxEta", 2.5,
                           "max eta absolute value of muon"}

Definition at line 130 of file RpcTrackAnaAlg.h.

                          {this, "maxEta", 2.5,
                           "max eta absolute value of muon"};

m_minDRTrackToGasGap

FloatProperty RpcTrackAnaAlg::m_minDRTrackToGasGap

private

Initial value:

{
        this, "minDRTrackToGasGap", 0.02,
        "minimum of DR between track and gasgap"}

Definition at line 156 of file RpcTrackAnaAlg.h.

                                      {
        this, "minDRTrackToGasGap", 0.02,
        "minimum of DR between track and gasgap"};

m_minPt

FloatProperty RpcTrackAnaAlg::m_minPt {this, "minPt", 25.0e3, "minmum pT of muon"}

private

Definition at line 129 of file RpcTrackAnaAlg.h.

{this, "minPt", 25.0e3, "minmum pT of muon"};

m_MuonContainerKey

SG::ReadHandleKey<xAOD::MuonContainer> RpcTrackAnaAlg::m_MuonContainerKey

private

Initial value:

{
        this, "MuonContainerKey", "Muons",
        "Key for Offline muon track Containers"}

Definition at line 183 of file RpcTrackAnaAlg.h.

                                                         {
        this, "MuonContainerKey", "Muons",
        "Key for Offline muon track Containers"};

m_muonMass

FloatProperty RpcTrackAnaAlg::m_muonMass

private

Initial value:

{this, "MuonMass", ParticleConstants::muonMassInMeV,
                             "muon invariant mass in MeV"}

Definition at line 138 of file RpcTrackAnaAlg.h.

                            {this, "MuonMass", ParticleConstants::muonMassInMeV,
                             "muon invariant mass in MeV"};

m_MuonRoIContainerKey

SG::ReadHandleKey<xAOD::MuonRoIContainer> RpcTrackAnaAlg::m_MuonRoIContainerKey

private

Initial value:

{
        this, "MuonRoIContainerName", "LVL1MuonRoIs", "Key for L1 ROIs"}

Definition at line 181 of file RpcTrackAnaAlg.h.

                                                               {
        this, "MuonRoIContainerName", "LVL1MuonRoIs", "Key for L1 ROIs"};

m_name

std::string AthMonitorAlgorithm::m_name

privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

m_outtime

FloatProperty RpcTrackAnaAlg::m_outtime {this, "outtime", 12.5, "the out-time time"}

private

Definition at line 170 of file RpcTrackAnaAlg.h.

{this, "outtime", 12.5, "the out-time time"};

m_packageName

StringProperty RpcTrackAnaAlg::m_packageName

private

Initial value:

{this, "PackageName", "RpcTrackAnaAlg",
                                 "group name for histograming"}

Definition at line 111 of file RpcTrackAnaAlg.h.

                                {this, "PackageName", "RpcTrackAnaAlg",
                                 "group name for histograming"};

m_plotMuonEff

BooleanProperty RpcTrackAnaAlg::m_plotMuonEff

private

Initial value:

{
        this, "plotMuonEff", false,
        "switch to plot histograms for Muon Efficiency"}

Definition at line 100 of file RpcTrackAnaAlg.h.

                                 {
        this, "plotMuonEff", false,
        "switch to plot histograms for Muon Efficiency"};

m_plotPRD

BooleanProperty RpcTrackAnaAlg::m_plotPRD

private

Initial value:

{
        this, "plotPRD", false,
        "switch to plot histograms for Prepare Data objects"}

Definition at line 103 of file RpcTrackAnaAlg.h.

                             {
        this, "plotPRD", false,
        "switch to plot histograms for Prepare Data objects"};

m_PrimaryVertexContainerKey

SG::ReadHandleKey<xAOD::VertexContainer> RpcTrackAnaAlg::m_PrimaryVertexContainerKey

private

Initial value:

{
        this, "PrimaryVertexContainerName", "PrimaryVertices",
        "Primary Vertex Container"}

Definition at line 188 of file RpcTrackAnaAlg.h.

                                                                    {
        this, "PrimaryVertexContainerName", "PrimaryVertices",
        "Primary Vertex Container"};

m_rpcIdHelper

const RpcIdHelper* RpcTrackAnaAlg::m_rpcIdHelper {nullptr}

private

Definition at line 176 of file RpcTrackAnaAlg.h.

{nullptr};

m_rpcPanelMap

RpcPanelMap RpcTrackAnaAlg::m_rpcPanelMap {}

private

Definition at line 205 of file RpcTrackAnaAlg.h.

{};

m_rpcPrdKey

SG::ReadHandleKey<Muon::RpcPrepDataContainer> RpcTrackAnaAlg::m_rpcPrdKey

private

Initial value:

{
        this, "RpcPrepDataContainer", "RPC_Measurements", "RPC PRDs"}

Definition at line 186 of file RpcTrackAnaAlg.h.

                                                         {
        this, "RpcPrepDataContainer", "RPC_Measurements", "RPC PRDs"};

m_SectorGroup

std::map<std::string, int> RpcTrackAnaAlg::m_SectorGroup {}

private

Definition at line 216 of file RpcTrackAnaAlg.h.

{};

m_StationNames

std::map<BarrelDL, std::vector<int> > RpcTrackAnaAlg::m_StationNames {}

private

Definition at line 213 of file RpcTrackAnaAlg.h.

{};

m_TagAndProbe

BooleanProperty RpcTrackAnaAlg::m_TagAndProbe

private

Initial value:

{this, "TagAndProbe", false,
                                  "switch to perform tag-and-probe method"}

Definition at line 122 of file RpcTrackAnaAlg.h.

                                 {this, "TagAndProbe", false,
                                  "switch to perform tag-and-probe method"};

m_TagAndProbeZmumu

BooleanProperty RpcTrackAnaAlg::m_TagAndProbeZmumu

private

Initial value:

{
        this, "TagAndProbeZmumu", false,
        "switch to perform tag-and-probe method Z->mumu"}

Definition at line 124 of file RpcTrackAnaAlg.h.

                                      {
        this, "TagAndProbeZmumu", false,
        "switch to perform tag-and-probe method Z->mumu"};

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_TriggerThrGroup

std::map<std::string, int> RpcTrackAnaAlg::m_TriggerThrGroup {}

private

Definition at line 217 of file RpcTrackAnaAlg.h.

{};

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_trigMatchWindow

FloatProperty RpcTrackAnaAlg::m_trigMatchWindow

private

Initial value:

{this, "TrigMatchingWindow", 0.005,
                                    "Window size in R for trigger matching"}

Definition at line 120 of file RpcTrackAnaAlg.h.

                                   {this, "TrigMatchingWindow", 0.005,
                                    "Window size in R for trigger matching"};

m_trigTagDefs

std::vector<TagDef> RpcTrackAnaAlg::m_trigTagDefs {}

private

Definition at line 207 of file RpcTrackAnaAlg.h.

{};

m_trigTagList

StringProperty RpcTrackAnaAlg::m_trigTagList

private

Initial value:

{
        this, "TagTrigList", "HLT_mu_ivarmedium;HLT_mu50",
        "list of triggers to be used for trigger matching"}

Definition at line 117 of file RpcTrackAnaAlg.h.

                                {
        this, "TagTrigList", "HLT_mu_ivarmedium;HLT_mu50",
        "list of triggers to be used for trigger matching"};

m_useAODParticle

BooleanProperty RpcTrackAnaAlg::m_useAODParticle

private

Initial value:

{this, "useAODParticle", false,
                                     "use AOD Particle"}

Definition at line 106 of file RpcTrackAnaAlg.h.

                                    {this, "useAODParticle", false,
                                     "use AOD Particle"};

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.

m_zMass_lowLimit

FloatProperty RpcTrackAnaAlg::m_zMass_lowLimit

private

Initial value:

{
        this, "zMass_lowLimit", 50000.,
        "2 muon invariant mass low limit in Zmumu event"}

Definition at line 140 of file RpcTrackAnaAlg.h.

                                  {
        this, "zMass_lowLimit", 50000.,
        "2 muon invariant mass low limit in Zmumu event"};

m_zMass_upLimit

FloatProperty RpcTrackAnaAlg::m_zMass_upLimit

private

Initial value:

{
        this, "zMass_upLimit", 150000.,
        "2 muon invariant mass up limit in Zmumu event"}

Definition at line 143 of file RpcTrackAnaAlg.h.

                                 {
        this, "zMass_upLimit", 150000.,
        "2 muon invariant mass up limit in Zmumu event"};

---

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

• RpcTrackAnaAlg.h
• RpcTrackAnaAlg.cxx