ATLAS Offline Software
Loading...
Searching...
No Matches
L1CaloL1TopoMonitorAlgorithm Class Reference

Monitor legacy CMX to L1Topo transmission of TOBs and timing of the L1Topo trigger bits with respect to L1 accept. More...

#include <L1CaloL1TopoMonitorAlgorithm.h>

Inheritance diagram for L1CaloL1TopoMonitorAlgorithm:
Collaboration diagram for L1CaloL1TopoMonitorAlgorithm:

Public Types

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

Public Member Functions

 L1CaloL1TopoMonitorAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~L1CaloL1TopoMonitorAlgorithm ()=default
virtual StatusCode initialize () override
 initialize
virtual StatusCode fillHistograms (const EventContext &ctx) const override
 adds event to the monitoring histograms
virtual StatusCode execute (const EventContext &ctx) const override
 Applies filters and trigger requirements.
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
 Fills a vector of variables to a group by reference.
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
 Fills a vector of variables to a group by reference.
template<typename... T>
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, T &&... variables) const
 Fills a variadic list of variables to a group by reference.
void fill (const std::string &groupName, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
 Fills a vector of variables to a group by name.
void fill (const std::string &groupName, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
 Fills a vector of variables to a group by name.
template<typename... T>
void fill (const std::string &groupName, T &&... variables) const
 Fills a variadic list of variables to a group by name.
Environment_t environment () const
 Accessor functions for the environment.
Environment_t envStringToEnum (const std::string &str) const
 Convert the environment string from the python configuration to an enum object.
DataType_t dataType () const
 Accessor functions for the data type.
DataType_t dataTypeStringToEnum (const std::string &str) const
 Convert the data type string from the python configuration to an enum object.
const ToolHandle< GenericMonitoringTool > & getGroup (const std::string &name) const
 Get a specific monitoring tool from the tool handle array.
const ToolHandle< Trig::TrigDecisionTool > & getTrigDecisionTool () const
 Get the trigger decision tool member.
bool trigChainsArePassed (const std::vector< std::string > &vTrigNames) const
 Check whether triggers are passed.
SG::ReadHandle< xAOD::EventInfoGetEventInfo (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

Public Attributes

 inputs = glob.glob('/eos/atlas/atlascerngroupdisk/data-art/build-output/master/Athena/x86_64-centos7-gcc8-opt/2020-04-06T2139/TrigP1Test/test_trigP1_v1PhysP1_T0Mon_build/ESD.pool.root')
 flags = initConfigFlags()
 Files
 HISTFileName
 cfg = MainServicesSerialCfg()
 L1CaloL1TopoMonitorCfg = L1CaloL1TopoMonitoringConfig(flags)
 OutputLevel
 withDetails
 False
 summariseProps
int nevents = -1

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< GenericMonitoringToolm_tools {this,"GMTools",{}}
 Array of Generic Monitoring Tools.
PublicToolHandle< Trig::TrigDecisionToolm_trigDecTool
 Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterToolm_DQFilterTools {this,"FilterTools",{}}
 Array of Data Quality filter tools.
SG::ReadCondHandleKey< LuminosityCondDatam_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
SG::ReadCondHandleKey< LBDurationCondDatam_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
SG::ReadCondHandleKey< TrigLiveFractionCondDatam_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::EventInfom_EventInfoKey {this,"EventInfoKey","EventInfo"}
 Key for retrieving EventInfo from StoreGate.

Private Types

enum class  SummaryErrors : int {
  CALO_CONV =0 , NO_CMX , DAQ_CONV , NO_DAQ ,
  ROI_CONV , F_CRC , PAYL_CRC , CMX_MATCH ,
  NUMBEROFBITS
}
enum  TOB_TYPE {
  JETL_TOB =0 , JETS_TOB , TAU_TOB , EM_TOB ,
  MU_TOB , TOB_TYPES
}
typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > MonVarVec_t
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

void jem2Coord (const int crate, const int jem, const int frame, const int location, int &x, int &y, double &eta, double &phi) const
void cpm2Coord (const int crate, const int cpm, const int chip, const int location, int &x, int &y, double &eta, double &phi) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

StringProperty m_packageName
SG::ReadHandleKey< CTP_RDOm_ctpRdoKey
SG::ReadHandleKey< xAOD::CMXCPTobContainerm_cmxCpTobLocation
SG::ReadHandleKey< xAOD::CMXJetTobContainerm_cmxJetTobLocation
SG::ReadHandleKey< xAOD::L1TopoRawDataContainerm_l1TopoKey
std::string m_name
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringToolm_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

Static Private Attributes

static const int MAXTOBS =30

Detailed Description

Monitor legacy CMX to L1Topo transmission of TOBs and timing of the L1Topo trigger bits with respect to L1 accept.

Author
Joergen Sjoelin

Definition at line 29 of file L1CaloL1TopoMonitorAlgorithm.h.

Member Typedef Documentation

◆ MonVarVec_t

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

Definition at line 370 of file AthMonitorAlgorithm.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ DataType_t

enum class AthMonitorAlgorithm::DataType_t
stronginherited

Specifies what type of input data is being monitored.

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

Enumerator
userDefined 
monteCarlo 
collisions 
cosmics 
heavyIonCollisions 

Definition at line 194 of file AthMonitorAlgorithm.h.

194 {
195 userDefined = 0,
196 monteCarlo,
197 collisions,
198 cosmics,
199 heavyIonCollisions,
200 };
AthConfigFlags cosmics(AthConfigFlags flags, str instanceName, str recoMode)

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

175 {
176 user = 0,
177 online,
178 tier0,
179 tier0Raw,
180 tier0ESD,
181 AOD,
182 altprod,
183 };

◆ SummaryErrors

enum class L1CaloL1TopoMonitorAlgorithm::SummaryErrors : int
strongprivate
Enumerator
CALO_CONV 
NO_CMX 
DAQ_CONV 
NO_DAQ 
ROI_CONV 
F_CRC 
PAYL_CRC 
CMX_MATCH 
NUMBEROFBITS 

Definition at line 46 of file L1CaloL1TopoMonitorAlgorithm.h.

46 : int {
47 CALO_CONV=0, NO_CMX, DAQ_CONV,
48 NO_DAQ, ROI_CONV, F_CRC, PAYL_CRC, CMX_MATCH, NUMBEROFBITS
49 };

◆ TOB_TYPE

enum L1CaloL1TopoMonitorAlgorithm::TOB_TYPE
private
Enumerator
JETL_TOB 
JETS_TOB 
TAU_TOB 
EM_TOB 
MU_TOB 
TOB_TYPES 

Definition at line 50 of file L1CaloL1TopoMonitorAlgorithm.h.

50{JETL_TOB=0, JETS_TOB, TAU_TOB, EM_TOB, MU_TOB, TOB_TYPES};

Constructor & Destructor Documentation

◆ L1CaloL1TopoMonitorAlgorithm()

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

Definition at line 26 of file L1CaloL1TopoMonitorAlgorithm.cxx.

28 : AthMonitorAlgorithm(name,pSvcLocator)
29{
30}
Base class for Athena Monitoring Algorithms.

◆ ~L1CaloL1TopoMonitorAlgorithm()

virtual L1CaloL1TopoMonitorAlgorithm::~L1CaloL1TopoMonitorAlgorithm ( )
virtualdefault

Member Function Documentation

◆ cardinality()

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

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

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ cpm2Coord()

void L1CaloL1TopoMonitorAlgorithm::cpm2Coord ( const int crate,
const int cpm,
const int chip,
const int location,
int & x,
int & y,
double & eta,
double & phi ) const
private

Definition at line 512 of file L1CaloL1TopoMonitorAlgorithm.cxx.

515{
516 x = crate*14 + cpm - 1;
517 y = chip*4 + location;
518 const uint32_t roiWord = // formula from CPSimMon
519 ((((((crate << 4) + cpm) << 4) + chip) << 2)
520 + location) << 18;
521 LVL1::CPRoIDecoder decoder;
522 const LVL1::CoordinateRange coord(decoder.coordinate(roiWord));
523 eta = coord.eta();
524 phi = coord.phi();
525}
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
double coord
Type of coordination system.
#define y
#define x
A level 1 calorimeter trigger conversion service: returns the Coordinate represented by a RoI word.
CoordinateRange class declaration.

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

224{ return m_dataType; }
AthMonitorAlgorithm::DataType_t m_dataType
Instance of the DataType_t enum.

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

144 {
145 // convert the string to all lowercase
146 std::string lowerCaseStr = str;
147 std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
148
149 // check if it matches one of the enum choices
150 if( lowerCaseStr == "userdefined" ) {
152 } else if( lowerCaseStr == "montecarlo" ) {
154 } else if( lowerCaseStr == "collisions" ) {
156 } else if( lowerCaseStr == "cosmics" ) {
157 return DataType_t::cosmics;
158 } else if( lowerCaseStr == "heavyioncollisions" ) {
160 } else { // otherwise, warn the user and return "userDefined"
161 ATH_MSG_WARNING("AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
162 <<str<<", returning userDefined.");
164 }
165}
#define ATH_MSG_WARNING(x)
void tolower(std::string &s)

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

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

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

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

208{ return m_environment; }
AthMonitorAlgorithm::Environment_t m_environment
Instance of the Environment_t enum.

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

116 {
117 // convert the string to all lowercase
118 std::string lowerCaseStr = str;
119 std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
120
121 // check if it matches one of the enum choices
122 if( lowerCaseStr == "user" ) {
123 return Environment_t::user;
124 } else if( lowerCaseStr == "online" ) {
126 } else if( lowerCaseStr == "tier0" ) {
128 } else if( lowerCaseStr == "tier0raw" ) {
130 } else if( lowerCaseStr == "tier0esd" ) {
132 } else if( lowerCaseStr == "aod" ) {
133 return Environment_t::AOD;
134 } else if( lowerCaseStr == "altprod" ) {
136 } else { // otherwise, warn the user and return "user"
137 ATH_MSG_WARNING("AthMonitorAlgorithm::envStringToEnum(): Unknown environment "
138 <<str<<", returning user.");
139 return Environment_t::user;
140 }
141}

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

◆ execute()

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

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters
ctxevent context for reentrant Athena call
Returns
StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

77 {
78
79 // Checks that all of the DQ filters are passed. If any one of the filters
80 // fails, return SUCCESS code and do not fill the histograms with the event.
81 for ( const auto& filterItr : m_DQFilterTools ) {
82 if (!filterItr->accept()) {
83 ATH_MSG_DEBUG("Event rejected due to filter tool.");
84 return StatusCode::SUCCESS;
85 }
86 }
87
88 // Trigger: If there is a decision tool and the chains fail, skip the event.
90 ATH_MSG_DEBUG("Event rejected due to trigger filter.");
91 return StatusCode::SUCCESS;
92 }
93
94 ATH_MSG_DEBUG("Event accepted!");
95 return fillHistograms(ctx);
96}
#define ATH_MSG_DEBUG(x)
virtual StatusCode fillHistograms(const EventContext &ctx) const =0
adds event to the monitoring histograms
bool trigChainsArePassed(const std::vector< std::string > &vTrigNames) const
Check whether triggers are passed.
std::vector< std::string > m_vTrigChainNames
Vector of trigger chain names parsed from trigger chain string.
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool
Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterTool > m_DQFilterTools
Array of Data Quality filter tools.

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

90{
91 // If we didn't find any symlinks to add, just return the collection
92 // from the base class. Otherwise, return the extended collection.
93 if (!m_extendedExtraObjects.empty()) {
95 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ fillHistograms()

StatusCode L1CaloL1TopoMonitorAlgorithm::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
ctxforwarded from execute
Returns
StatusCode

Implements AthMonitorAlgorithm.

Definition at line 46 of file L1CaloL1TopoMonitorAlgorithm.cxx.

46 {
47
48 ATH_MSG_DEBUG("L1CaloL1TopoMonitorAlgorithm::fillHistograms");
49 ATH_MSG_DEBUG("m_ctpRdoKey=" << m_ctpRdoKey);
50 ATH_MSG_DEBUG("m_cmxCpTobLocation=" << m_cmxCpTobLocation);
51 ATH_MSG_DEBUG("m_cmxJetTobLocation=" << m_cmxJetTobLocation);
52 ATH_MSG_DEBUG("m_l1TopoKey=" << m_l1TopoKey);
53
54 typedef std::tuple<int,int,int,int,int,int> TobKey;
55 std::set<TobKey> cmxKeys[TOB_TYPES],topoKeys[TOB_TYPES],
56 keyDiff[TOB_TYPES],keyIntersect[TOB_TYPES];
57
58 const int LUMI_BLOCK=GetEventInfo(ctx)->lumiBlock();
59 const int NFPGA=4;
60 const int MAXTOBS=20;
61
62 // monitored variables
63 std::vector<int> errors;
64 std::vector<std::reference_wrapper<Monitored::IMonitoredVariable>> variables;
65
66 // 1D
67 auto run = Monitored::Scalar<int>("run",GetEventInfo(ctx)->runNumber());
68 auto nJetTobs = Monitored::Scalar<int>("nJetTobs",0);
69 auto nTauTobs = Monitored::Scalar<int>("nTauTobs",0);
70 auto nEMTobs = Monitored::Scalar<int>("nEMTobs",0);
71 auto nMuonTobs = Monitored::Scalar<int>("nMuonTobs",0);
72
73 variables.push_back(run);
74 variables.push_back(nJetTobs);
75 variables.push_back(nTauTobs);
76 variables.push_back(nEMTobs);
77 variables.push_back(nMuonTobs);
78
79 auto jetTobsEnergyL = Monitored::Scalar<int>("jetTobsEnergyL",0);
80 auto wFPGA = Monitored::Scalar<float>("wFPGA",1.0/NFPGA);
81 auto lbErrors = Monitored::Scalar<int>("lbErrors",0);
82 auto sumErrors = Monitored::Collection("sumErrors",errors,[](const auto &_e){return _e;});
83
84 // 2D Hitmaps
85 auto etaJetSTobs_match = Monitored::Scalar<double>("etaJetSTobs_match", 0.);
86 auto phiJetSTobs_match = Monitored::Scalar<double>("phiJetSTobs_match", 0.);
87 auto etaJetLTobs_match = Monitored::Scalar<double>("etaJetLTobs_match", 0.);
88 auto phiJetLTobs_match = Monitored::Scalar<double>("phiJetLTobs_match", 0.);
89 auto etaTauTobs_match = Monitored::Scalar<double>("etaTauTobs_match", 0.);
90 auto phiTauTobs_match = Monitored::Scalar<double>("phiTauTobs_match", 0.);
91 auto etaEMTobs_match = Monitored::Scalar<double>("etaEMTobs_match", 0.);
92 auto phiEMTobs_match = Monitored::Scalar<double>("phiEMTobs_match", 0.);
93 auto xJetSTobs_match = Monitored::Scalar<double>("xJetSTobs_match", 0.);
94 auto yJetSTobs_match = Monitored::Scalar<double>("yJetSTobs_match", 0.);
95 auto xJetLTobs_match = Monitored::Scalar<double>("xJetLTobs_match", 0.);
96 auto yJetLTobs_match = Monitored::Scalar<double>("yJetLTobs_match", 0.);
97 auto xTauTobs_match = Monitored::Scalar<double>("xTauTobs_match", 0.);
98 auto yTauTobs_match = Monitored::Scalar<double>("yTauTobs_match", 0.);
99 auto xEMTobs_match = Monitored::Scalar<double>("xEMTobs_match", 0.);
100 auto yEMTobs_match = Monitored::Scalar<double>("yEMTobs_match", 0.);
101
102 auto etaJetSTobs_mismatch = Monitored::Scalar<double>("etaJetSTobs_mismatch", 0.);
103 auto phiJetSTobs_mismatch = Monitored::Scalar<double>("phiJetSTobs_mismatch", 0.);
104 auto etaJetLTobs_mismatch = Monitored::Scalar<double>("etaJetLTobs_mismatch", 0.);
105 auto phiJetLTobs_mismatch = Monitored::Scalar<double>("phiJetLTobs_mismatch", 0.);
106 auto etaTauTobs_mismatch = Monitored::Scalar<double>("etaTauTobs_mismatch", 0.);
107 auto phiTauTobs_mismatch = Monitored::Scalar<double>("phiTauTobs_mismatch", 0.);
108 auto etaEMTobs_mismatch = Monitored::Scalar<double>("etaEMTobs_mismatch", 0.);
109 auto phiEMTobs_mismatch = Monitored::Scalar<double>("phiEMTobs_mismatch", 0.);
110 auto xJetSTobs_mismatch = Monitored::Scalar<double>("xJetSTobs_mismatch", 0.);
111 auto yJetSTobs_mismatch = Monitored::Scalar<double>("yJetSTobs_mismatch", 0.);
112 auto xJetLTobs_mismatch = Monitored::Scalar<double>("xJetLTobs_mismatch", 0.);
113 auto yJetLTobs_mismatch = Monitored::Scalar<double>("yJetLTobs_mismatch", 0.);
114 auto xTauTobs_mismatch = Monitored::Scalar<double>("xTauTobs_mismatch", 0.);
115 auto yTauTobs_mismatch = Monitored::Scalar<double>("yTauTobs_mismatch", 0.);
116 auto xEMTobs_mismatch = Monitored::Scalar<double>("xEMTobs_mismatch", 0.);
117 auto yEMTobs_mismatch = Monitored::Scalar<double>("yEMTobs_mismatch", 0.);
118
119 // 2D Topo items vs BC
120 auto item0 = Monitored::Scalar<int>("item0", 0.);
121 auto item1 = Monitored::Scalar<int>("item1", 0.);
122 auto item2 = Monitored::Scalar<int>("item2", 0.);
123 auto item3 = Monitored::Scalar<int>("item3", 0.);
124 auto bc0 = Monitored::Scalar<int>("bc0", 0.);
125 auto bc1 = Monitored::Scalar<int>("bc1", 0.);
126 auto bc2 = Monitored::Scalar<int>("bc2", 0.);
127 auto bc3 = Monitored::Scalar<int>("bc3", 0.);
128
129 // read all objects needed
130 ATH_MSG_DEBUG("Begin to fetch store gate objects ..");
131 SG::ReadHandle<CTP_RDO> const_ctpRdo(m_ctpRdoKey,ctx);
132 ATH_CHECK(const_ctpRdo.isValid());
133 SG::ReadHandle<xAOD::CMXCPTobContainer> cmxCpTob(m_cmxCpTobLocation,ctx);
134 ATH_CHECK(cmxCpTob.isValid());
135 SG::ReadHandle<xAOD::CMXJetTobContainer> cmxJetTob(m_cmxJetTobLocation,ctx);
136 ATH_CHECK(cmxJetTob.isValid());
137 SG::ReadHandle<xAOD::L1TopoRawDataContainer> l1TopoRawData{m_l1TopoKey,ctx};
138 if (!l1TopoRawData.isValid()) {
139 ATH_MSG_DEBUG ("L1Topo DAQ raw data collection not valid");
140 errors.push_back(static_cast<int>(SummaryErrors::NO_DAQ));
141 lbErrors=LUMI_BLOCK;
142 fill(m_packageName,lbErrors,sumErrors);
143 return StatusCode::SUCCESS;
144 }
145 if (l1TopoRawData->empty()) {
146 ATH_MSG_DEBUG ("L1Topo DAQ raw data collection empty");
147 errors.push_back(static_cast<int>(SummaryErrors::NO_DAQ));
148 lbErrors=LUMI_BLOCK;
149 fill(m_packageName,lbErrors,sumErrors);
150 return StatusCode::SUCCESS;
151 }
152 ATH_MSG_DEBUG("Done fetching from store gate.");
153
154 // Make a writable copy and Set CTP version number to 4 when reading persistified data
155 std::vector<uint32_t> ctp_data=const_ctpRdo->getDataWords();
156 CTP_RDO ctpRdo(4,std::move(ctp_data));
157 ctpRdo.setL1AcceptBunchPosition(const_ctpRdo->getL1AcceptBunchPosition());
158 ctpRdo.setTurnCounter(const_ctpRdo->getTurnCounter());
159
160 // Approximate timing test of L1Topo by comparing TIP bits to L1A
161 if (ctpRdo.getCTPVersionNumber()==0) {
162 ATH_MSG_DEBUG("CTP version number not set, skipping CTP test");
163 }
164 else {
165 CTP_Decoder ctp;
166 ctp.setRDO(&ctpRdo);
167 const uint16_t l1aPos = ctpRdo.getL1AcceptBunchPosition();
168 const uint32_t bx=ctp.getBunchCrossings().size();
169 const int dbx=bx/2;
170 if (l1aPos >= bx) {
171 ATH_MSG_DEBUG( "CTP_RDO gave invalid l1aPos." );
172 }
173 else {
174 ATH_MSG_DEBUG( "CTP l1aPos, size : " << l1aPos << ", " << bx );
175 if ( bx>=3 && l1aPos>0) {
176 for (int bc=-dbx; bc<=dbx; ++bc) {
177 const CTP_BC& ctp_bc = ctp.getBunchCrossings().at(l1aPos+bc);
178 std::bitset<512> tip = ctp_bc.getTIP();
179 ATH_MSG_DEBUG( "got CTP TIP bits: " << tip.to_string() );
180 const unsigned int topoTipStart(384);
181 const unsigned int nTopoCTPOutputs(128);
182 for (unsigned int item=0; item<nTopoCTPOutputs;++item) {
183 // Check if algorithms fired
184 if (tip.test(item+topoTipStart)) {
185 if (item/32==0) {item0=item; bc0=bc; fill(m_packageName,item0,bc0);}
186 else if (item/32==1) {item1=item; bc1=bc; fill(m_packageName,item1,bc1);}
187 else if (item/32==2) {item2=item; bc2=bc; fill(m_packageName,item2,bc2);}
188 else if (item/32==3) {item3=item; bc3=bc; fill(m_packageName,item3,bc3);}
189 ATH_MSG_DEBUG(" Fired (item, bc) =" << item << ", " << bc);
190 }
191 }
192 }
193 }
194 }
195 }
196
197 // Retrieve CMX CP tobs
198 if (cmxCpTob->empty()) {
199 ATH_MSG_DEBUG ("No CMX CP TOBs found");
200 }
201 else {
202 ATH_MSG_DEBUG( "Found CMXCPTobCollection, looping on TOBs ..." );
203 for (const xAOD::CMXCPTob* t : *cmxCpTob) {
204 if (t->energy()) {
205 for (int clone=0;clone<4;++clone) {
206 if (t->cmx()==0)
207 cmxKeys[TAU_TOB].insert(std::make_tuple(t->crate(),t->cpm(),
208 t->chip(),t->location(),
209 t->energy(),clone));
210 else
211 cmxKeys[EM_TOB].insert(std::make_tuple(t->crate(),t->cpm(),
212 t->chip(),t->location(),
213 t->energy(),clone));
214 }
215 }
216 }
217 }
218
219 // Retrieve CMX jet tobs
220 if (cmxJetTob->empty()) {
221 ATH_MSG_DEBUG ("No CMX jet TOBs found");
222 }
223 else {
224 ATH_MSG_DEBUG( "Found CMXJetTobCollection, looping on TOBs ..." );
225 for (const xAOD::CMXJetTob* t : *cmxJetTob) {
226 for (int clone=0;clone<4;++clone) {
227 if (t->energyLarge())
228 cmxKeys[JETL_TOB].insert(std::make_tuple(t->crate(),t->jem(),
229 t->frame(),t->location(),
230 t->energyLarge(),clone));
231 if (t->energySmall())
232 cmxKeys[JETS_TOB].insert(std::make_tuple(t->crate(),t->jem(),
233 t->frame(),t->location(),
234 t->energySmall(),clone));
235 }
236 }
237 }
238
239 // analyse and register the L1Topo DAQ TOBs
240 int daqJetTobs=0;
241 int daqTauTobs=0;
242 int daqEMTobs=0;
243 int daqMuonTobs=0;
244 for (const xAOD::L1TopoRawData* rawdata : *l1TopoRawData) {
245 ATH_MSG_DEBUG("Found L1Topo raw DAQ data with source ID "
246 << L1Topo::formatHex8(rawdata->sourceID()));
247 //rdo.setStatusWords(rawdata->statusWords());
248 if (rawdata->errorWord()) {
249 ATH_MSG_INFO("DAQ Converter errors reported: " << rawdata->errorWord());
250 errors.push_back(static_cast<int>(SummaryErrors::DAQ_CONV));
251 }
252 if (rawdata->dataWords().empty()) {
253 ATH_MSG_INFO("L1Topo raw data DAQ payload is empty");
254 errors.push_back(static_cast<int>(SummaryErrors::NO_DAQ));
255 }
256 else {
257 std::vector<L1Topo::L1TopoTOB> daqTobs;
258
259 // initialise header: beware, this can make a valid-looking header
260 // set version 15, BCN -7, which is unlikely:
261 L1Topo::Header header(0xf, 0, 0, 0, 0, 1, 0x7);
262 int i_fpga=-1;
263 for (auto word : rawdata->dataWords()) {
264 switch (L1Topo::blockType(word)) {
265
267 {
268 header = L1Topo::Header(word);
269 if (header.payload_crc()!=0) {
270 errors.push_back(static_cast<int>(SummaryErrors::PAYL_CRC));
271 }
272 i_fpga=(((rawdata->sourceID())>>3)&2)+header.fpga();
273 break;
274 }
276 {
277 auto fibreBlock = L1Topo::Fibre(word);
278 break;
279 }
281 {
282 auto status = L1Topo::Status(word);
283 ATH_MSG_WARNING( "fibre overflow: " << status.overflow()
284 << " fibre crc: " << status.crc() );
285 if (status.crc()) {
286 errors.push_back(static_cast<int>(SummaryErrors::F_CRC));
287 }
288 break;
289 }
291 {
292 if (header.bcn_offset()==0){
293 auto tob = L1Topo::L1TopoTOB(word);
294 ATH_MSG_DEBUG(tob);
295 daqTobs.push_back(tob);
296 // not reimplemented yet
297 //auto index = L1Topo::triggerBitIndex(rawdata->sourceID(),tob);
298 for (unsigned int i=0; i<8; ++i){
299 if ((tob.trigger_bits() >>i)&1) {
300 //m_h_l1topo_1d_DAQTriggerBits->Fill(index+i);
301 }
302 if ((tob.overflow_bits()>>i)&1) {
303 //m_h_l1topo_1d_DAQOverflowBits->Fill(index+i);
304 }
305 }
306 }
307 break;
308 }
311 {
312 if (header.bcn_offset()==0) {
313 const int crate = (word >> 28) & 0x1;
314 const int jem = (word >> 24) & 0xF;
315 const int frame = (word >> 21) & 0x7;
316 const int location = (word >> 19) & 0x3;
317 const int energyS = (word >> 10) & 0x1FF;
318 const int energyL = (word & 0x3FF);
319 if (energyL) {
320 topoKeys[JETL_TOB].
321 insert(std::make_tuple(crate,jem,frame,
322 location,energyL,i_fpga));
323 ++daqJetTobs;
324 jetTobsEnergyL=energyL;
325 fill(m_packageName,jetTobsEnergyL,wFPGA);
326 }
327 if (energyS) {
328 topoKeys[JETS_TOB].
329 insert(std::make_tuple(crate,jem,frame,
330 location,energyS,i_fpga));
331 }
332 }
333 break;
334 }
336 {
337 if (header.bcn_offset()==0) {
338 ++daqTauTobs;
339 const int crate = (word >> 26) & 0x3;
340 const int cpm = (word >> 20) & 0xF;
341 const int chip = (word >> 15) & 0xF;
342 const int location = (word >> 13) & 0x3;
343 const int energy = (word & 0xFF);
344 if (energy)
345 topoKeys[TAU_TOB].
346 insert(std::make_tuple(crate,cpm,chip,
347 location,energy,i_fpga));
348 }
349 break;
350 }
352 {
353 if (header.bcn_offset()==0) {
354 ++daqEMTobs;
355 const int crate = (word >> 26) & 0x3;
356 const int cpm = (word >> 20) & 0xF;
357 const int chip = (word >> 15) & 0xF;
358 const int location = (word >> 13) & 0x3;
359 const int energy = (word & 0xFF);
360 if (energy)
361 topoKeys[EM_TOB].
362 insert(std::make_tuple(crate,cpm,chip,
363 location,energy,i_fpga));
364 }
365 break;
366 }
368 {
369 if (header.bcn_offset()==0){
370 ++daqMuonTobs;
371 }
372 break;
373 }
375 {
376 if (header.bcn_offset()==0){
377 }
378 break;
379 }
380 default:
381 {
382 break;
383 }
384
385 }
386 } // for word
387 } // L1Topo raw data
388 } // L1 Topo collection
389
390
391 // look in both directions for CMX and Topo TOBs (mis)matches
392 bool mismatch=false;
393 for (int t=0; t<TOB_TYPES; ++t) {
394 // create set of mismatches
395 set_symmetric_difference(cmxKeys[t].begin(),cmxKeys[t].end(),
396 topoKeys[t].begin(),topoKeys[t].end(),
397 inserter(keyDiff[t],keyDiff[t].begin()));
398 if (keyDiff[t].size()>0) {
399 ATH_MSG_DEBUG("Error: CMX <-> L1Topo TOB mismatches = "
400 << keyDiff[t].size());
401 mismatch=true;
402 }
403 for (auto& tob : keyDiff[t]) {
404 int x,y;
405 double eta,phi;
406 if (t==JETS_TOB || t==JETL_TOB) {
407 jem2Coord(std::get<0>(tob),std::get<1>(tob),
408 std::get<2>(tob),std::get<3>(tob),x,y,eta,phi);
409 if (t==JETS_TOB) {
410 xJetSTobs_mismatch=x; yJetSTobs_mismatch=y;
411 etaJetSTobs_mismatch=eta; phiJetSTobs_mismatch=phi;
412 fill(m_packageName, xJetSTobs_mismatch, yJetSTobs_mismatch,
413 etaJetSTobs_mismatch, phiJetSTobs_mismatch);
414 }
415 else { // JETL_TOB
416 xJetLTobs_mismatch=x; yJetLTobs_mismatch=y;
417 etaJetLTobs_mismatch=eta; phiJetLTobs_mismatch=phi;
418 fill(m_packageName, xJetLTobs_mismatch, yJetLTobs_mismatch,
419 etaJetLTobs_mismatch, phiJetLTobs_mismatch);
420 }
421 }
422 else {
423 cpm2Coord(std::get<0>(tob),std::get<1>(tob),
424 std::get<2>(tob),std::get<3>(tob),x,y,eta,phi);
425 if (t==TAU_TOB) {
426 xTauTobs_mismatch=x; yTauTobs_mismatch=y;
427 etaTauTobs_mismatch=eta; phiTauTobs_mismatch=phi;
428 fill(m_packageName, xTauTobs_mismatch, yTauTobs_mismatch,
429 etaTauTobs_mismatch, phiTauTobs_mismatch);
430 }
431 else { // EM_TOB
432 xEMTobs_mismatch=x; yEMTobs_mismatch=y;
433 etaEMTobs_mismatch=eta; phiEMTobs_mismatch=phi;
434 fill(m_packageName, xEMTobs_mismatch, yEMTobs_mismatch,
435 etaEMTobs_mismatch, phiEMTobs_mismatch);
436 }
437 }
438 }
439 // create set of matches
440 set_intersection(cmxKeys[t].begin(),cmxKeys[t].end(),
441 topoKeys[t].begin(),topoKeys[t].end(),
442 inserter(keyIntersect[t],keyIntersect[t].begin()));
443 for (auto& tob : keyIntersect[t]) {
444 int x,y;
445 double eta,phi;
446 if (t==JETS_TOB || t==JETL_TOB) {
447 jem2Coord(std::get<0>(tob),std::get<1>(tob),
448 std::get<2>(tob),std::get<3>(tob),x,y,eta,phi);
449 if (t==JETS_TOB) {
450 xJetSTobs_match=x; yJetSTobs_match=y;
451 etaJetSTobs_match=eta; phiJetSTobs_match=phi;
452 fill(m_packageName, xJetSTobs_match, yJetSTobs_match,
453 etaJetSTobs_match, phiJetSTobs_match);
454 }
455 else { // JETL_TOB
456 xJetLTobs_match=x; yJetLTobs_match=y;
457 etaJetLTobs_match=eta; phiJetLTobs_match=phi;
458 fill(m_packageName, xJetLTobs_match, yJetLTobs_match,
459 etaJetLTobs_match, phiJetLTobs_match);
460 }
461 }
462 else {
463 cpm2Coord(std::get<0>(tob),std::get<1>(tob),
464 std::get<2>(tob),std::get<3>(tob),x,y,eta,phi);
465 if (t==TAU_TOB) {
466 xTauTobs_match=x; yTauTobs_match=y;
467 etaTauTobs_match=eta; phiTauTobs_match=phi;
468 fill(m_packageName, xTauTobs_match, yTauTobs_match,
469 etaTauTobs_match, phiTauTobs_match);
470 }
471 else { // EM_TOB
472 xEMTobs_match=x; yEMTobs_match=y;
473 etaEMTobs_match=eta; phiEMTobs_match=phi;
474 fill(m_packageName, xEMTobs_match, yEMTobs_match,
475 etaEMTobs_match, phiEMTobs_match);
476 }
477 }
478 }
479 }
480
481 if (mismatch) {
482 errors.push_back(static_cast<int>(SummaryErrors::CMX_MATCH));
483 }
484 if (!errors.empty()) {
485 lbErrors=LUMI_BLOCK;
486 fill(m_packageName,lbErrors,sumErrors);
487 }
488
489 nJetTobs=std::min(daqJetTobs/NFPGA,MAXTOBS-1);
490 nTauTobs=std::min(daqTauTobs/NFPGA,MAXTOBS-1);
491 nEMTobs=std::min(daqEMTobs/NFPGA,MAXTOBS-1);
492 nMuonTobs=std::min(daqMuonTobs/NFPGA,MAXTOBS-1);
493
494 fill(m_packageName,variables);
495 return StatusCode::SUCCESS;
496}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_INFO(x)
const std::bitset< 512 > & getTIP() const
get bitset of TIP words
Definition CTP_Decoder.h:80
void cpm2Coord(const int crate, const int cpm, const int chip, const int location, int &x, int &y, double &eta, double &phi) const
void jem2Coord(const int crate, const int jem, const int frame, const int location, int &x, int &y, double &eta, double &phi) const
Declare a monitored scalar variable.
virtual bool isValid() override final
Can the handle be successfully dereferenced?
L1Topo::BlockTypes blockType(const uint32_t word, uint32_t offset=28, uint32_t size=0x0f)
Function to return the block type of a data word from L1Topo.
Definition BlockTypes.cxx:9
std::string formatHex8(uint32_t word)
Helper function to format a 32-bit integer as an 8-digit hex number for printing.
ValuesCollection< T > Collection(std::string name, const T &collection)
Declare a monitored (double-convertible) collection.
ctp
Retrieve trigger EDM objects.
status
Definition merge.py:16
Definition run.py:1
L1TopoRawData_v1 L1TopoRawData
Define the latest version of the L1TopoRawData class.
CMXJetTob_v1 CMXJetTob
Define the latest version of the CMXJetTob class.
CMXCPTob_v1 CMXCPTob
Define the latest version of the CMXCPTob class.
setWord1 uint16_t
setEventNumber uint32_t
Set * set_intersection(Set *set1, Set *set2)
Perform an intersection of two sets.
void fill(H5::Group &out_file, size_t iterations)

◆ filterPassed()

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

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ 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
ctxEventContext for the event
Returns
a SG::ReadHandle<xAOD::EventInfo>

Definition at line 111 of file AthMonitorAlgorithm.cxx.

111 {
112 return SG::ReadHandle<xAOD::EventInfo>(m_EventInfoKey, ctx);
113}
SG::ReadHandleKey< xAOD::EventInfo > m_EventInfoKey
Key for retrieving EventInfo from StoreGate.

◆ 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
namestring name of the desired tool
Returns
reference to the desired monitoring tool

Definition at line 168 of file AthMonitorAlgorithm.cxx.

168 {
169 // get the pointer to the tool, and check that it exists
170 auto idx = m_toolLookupMap.find(name);
171 if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
172 return m_tools[idx->second];
173 }
174 else {
175 // treat empty tool handle case as in Monitored::Group
176 if (m_toolLookupMap.empty()) {
177 return m_dummy;
178 }
179
180 if (!isInitialized()) {
182 "It seems that the AthMonitorAlgorithm::initialize was not called "
183 "in derived class initialize method");
184 } else {
185 std::string available = std::accumulate(
186 m_toolLookupMap.begin(), m_toolLookupMap.end(), std::string(""),
187 [](const std::string& s, auto h) { return s + "," + h.first; });
188 ATH_MSG_FATAL("The tool " << name << " could not be found in the tool array of the "
189 << "monitoring algorithm " << m_name << ". This probably reflects a discrepancy between "
190 << "your python configuration and c++ filling code. Note: your available groups are {"
191 << available << "}.");
192 }
193 }
194 return m_dummy;
195}
#define ATH_MSG_FATAL(x)
#define ATH_LIKELY(x)
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringTool > m_dummy
ToolHandleArray< GenericMonitoringTool > m_tools
Array of Generic Monitoring Tools.

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

198 {
199 return m_trigDecTool;
200}

◆ initialize()

StatusCode L1CaloL1TopoMonitorAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 32 of file L1CaloL1TopoMonitorAlgorithm.cxx.

32 {
33
34 ATH_MSG_DEBUG("L1CaloL1TopoMonitorAlgorith::initialize");
35 ATH_MSG_DEBUG("Package Name "<< m_packageName);
36
37 // we initialise all the containers that we need
38 ATH_CHECK(m_ctpRdoKey.initialize());
39 ATH_CHECK(m_cmxJetTobLocation.initialize());
40 ATH_CHECK(m_cmxCpTobLocation.initialize());
41 ATH_CHECK(m_l1TopoKey.initialize());
42
44}
virtual StatusCode initialize() override
initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.

◆ 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()

◆ jem2Coord()

void L1CaloL1TopoMonitorAlgorithm::jem2Coord ( const int crate,
const int jem,
const int frame,
const int location,
int & x,
int & y,
double & eta,
double & phi ) const
private

Definition at line 498 of file L1CaloL1TopoMonitorAlgorithm.cxx.

501{
502 x = crate*16 + jem;
503 y = frame*4 + location;
504 const uint32_t roiWord = // formula from JEPSimMon
505 ((((((crate << 4) + jem) << 3) + frame) << 2) + location) << 19;
506 LVL1::JEPRoIDecoder decoder;
507 const LVL1::CoordinateRange coord(decoder.coordinate(roiWord));
508 eta = coord.eta();
509 phi = coord.phi();
510}
A level 1 calorimeter trigger conversion service: returns the Coordinate represented by a RoI word.

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ 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
lineThe input string.
resultThe parsed output vector of strings.
Returns
StatusCode

Definition at line 345 of file AthMonitorAlgorithm.cxx.

345 {
346 std::string item;
347 std::stringstream ss(line);
348
349 ATH_MSG_DEBUG( "AthMonitorAlgorithm::parseList()" );
350
351 while ( std::getline(ss, item, ',') ) {
352 std::stringstream iss(item); // remove whitespace
353 iss >> item;
354 result.push_back(item);
355 }
356
357 return StatusCode::SUCCESS;
358}
static Double_t ss

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

381 {
382 h.renounce();
384 }
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)

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

364 {
366 }

◆ setFilterPassed()

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

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

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

77{
78 return BaseAlg::sysExecute (ctx);
79}

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

107 {
109
110 if (sc.isFailure()) {
111 return sc;
112 }
113
114 ServiceHandle<ICondSvc> cs("CondSvc",name());
115 for (auto h : outputHandles()) {
116 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117 // do this inside the loop so we don't create the CondSvc until needed
118 if ( cs.retrieve().isFailure() ) {
119 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120 return StatusCode::SUCCESS;
121 }
122 if (cs->regHandle(this,*h).isFailure()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
#define ATH_MSG_ERROR(x)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

◆ 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
vTrigNamesList 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.

203 {
204
205
206 // If no triggers were given, return true.
207 if (vTrigNames.empty()) return true;
208
209
210 // Trigger: Check if this Algorithm is being run as an Express Stream job.
211 // Events are entering the express stream are chosen randomly, and by chain,
212 // Hence an additional check should be aplied to see if the chain(s)
213 // monitored here are responsible for the event being selected for
214 // the express stream.
215
216 const auto group = m_trigDecTool->getChainGroup(vTrigNames);
218 const auto passedBits = m_trigDecTool->isPassedBits(group);
219 bool expressPass = passedBits & TrigDefs::Express_passed; //bitwise AND
220 if(!expressPass) {
221 return false;
222 }
223 }
224
225 // monitor the event if any of the chains in the chain group passes the event.
226 return group->isPassed();
227
228}
Gaudi::Property< bool > m_enforceExpressTriggers

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }

Member Data Documentation

◆ cfg

L1CaloL1TopoMonitorAlgorithm.cfg = MainServicesSerialCfg()

Definition at line 215 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ False

L1CaloL1TopoMonitorAlgorithm.False

Definition at line 224 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ Files

L1CaloL1TopoMonitorAlgorithm.Files

Definition at line 207 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ flags

L1CaloL1TopoMonitorAlgorithm.flags = initConfigFlags()

Definition at line 206 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ HISTFileName

L1CaloL1TopoMonitorAlgorithm.HISTFileName

Definition at line 208 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ inputs

L1CaloL1TopoMonitorAlgorithm.inputs = glob.glob('/eos/atlas/atlascerngroupdisk/data-art/build-output/master/Athena/x86_64-centos7-gcc8-opt/2020-04-06T2139/TrigP1Test/test_trigP1_v1PhysP1_T0Mon_build/ESD.pool.root')

Definition at line 204 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ L1CaloL1TopoMonitorCfg

L1CaloL1TopoMonitorAlgorithm.L1CaloL1TopoMonitorCfg = L1CaloL1TopoMonitoringConfig(flags)

Definition at line 218 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ m_cmxCpTobLocation

SG::ReadHandleKey<xAOD::CMXCPTobContainer> L1CaloL1TopoMonitorAlgorithm::m_cmxCpTobLocation
private
Initial value:
{this,
"CMXCPTob container"}
static const std::string CMXCPTobLocation

Definition at line 56 of file L1CaloL1TopoMonitorAlgorithm.h.

56 {this,
58 "CMXCPTob container"};

◆ m_cmxJetTobLocation

SG::ReadHandleKey<xAOD::CMXJetTobContainer> L1CaloL1TopoMonitorAlgorithm::m_cmxJetTobLocation
private
Initial value:
{this,
"CMXJetTob container"}
static const std::string CMXJetTobLocation

Definition at line 59 of file L1CaloL1TopoMonitorAlgorithm.h.

59 {this,
60 "CMXJetTobLocation", LVL1::TrigT1CaloDefs::CMXJetTobLocation,
61 "CMXJetTob container"};

◆ m_ctpRdoKey

SG::ReadHandleKey<CTP_RDO> L1CaloL1TopoMonitorAlgorithm::m_ctpRdoKey
private
Initial value:
{this, "CTPRDOLocation",
LVL1CTP::DEFAULT_RDOOutputLocation,"Key of the CTP RDO object"}
static const std::string DEFAULT_RDOOutputLocation
location of CTP RDO output in StoreGate

Definition at line 54 of file L1CaloL1TopoMonitorAlgorithm.h.

54 {this, "CTPRDOLocation",
55 LVL1CTP::DEFAULT_RDOOutputLocation,"Key of the CTP RDO object"};

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

358{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.

365{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.

366{this,"DetailLevel",0};

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_DQFilterTools

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

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

346{this,"FilterTools",{}};

◆ m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy
privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

◆ m_enforceExpressTriggers

Gaudi::Property<bool> AthMonitorAlgorithm::m_enforceExpressTriggers
privateinherited
Initial value:
{this,
"EnforceExpressTriggers", false,
"Requires that matched triggers made the event enter the express stream"}

Definition at line 377 of file AthMonitorAlgorithm.h.

377 {this,
378 "EnforceExpressTriggers", false,
379 "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.

357{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.

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

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_fileKey

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

Internal Athena name for file.

Definition at line 363 of file AthMonitorAlgorithm.h.

363{this,"FileKey",""};

◆ m_l1TopoKey

SG::ReadHandleKey<xAOD::L1TopoRawDataContainer> L1CaloL1TopoMonitorAlgorithm::m_l1TopoKey
private
Initial value:
{this,
"L1TopoRawData", "L1TopoRawData",
"Key of the L1Topo DAQ object"}

Definition at line 62 of file L1CaloL1TopoMonitorAlgorithm.h.

62 {this,
63 "L1TopoRawData", "L1TopoRawData",
64 "Key of the L1Topo DAQ object"};

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

351{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.

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

◆ m_name

std::string AthMonitorAlgorithm::m_name
privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

◆ m_packageName

StringProperty L1CaloL1TopoMonitorAlgorithm::m_packageName
private
Initial value:
{this,"PackageName","L1CaloL1TopoMonitor",
"group name for histograming"}

Definition at line 43 of file L1CaloL1TopoMonitorAlgorithm.h.

43 {this,"PackageName","L1CaloL1TopoMonitor",
44 "group name for histograming"};

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

341{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.

360{this,"TriggerChain",""};

◆ m_trigLiveFractionDataKey

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

Definition at line 352 of file AthMonitorAlgorithm.h.

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

◆ m_useLumi

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

Allows use of various luminosity functions.

Definition at line 364 of file AthMonitorAlgorithm.h.

364{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.

◆ MAXTOBS

const int L1CaloL1TopoMonitorAlgorithm::MAXTOBS =30
staticprivate

Definition at line 51 of file L1CaloL1TopoMonitorAlgorithm.h.

◆ nevents

int L1CaloL1TopoMonitorAlgorithm.nevents = -1

Definition at line 226 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ OutputLevel

L1CaloL1TopoMonitorAlgorithm.OutputLevel

Definition at line 222 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ summariseProps

L1CaloL1TopoMonitorAlgorithm.summariseProps

Definition at line 224 of file L1CaloL1TopoMonitorAlgorithm.py.

◆ withDetails

L1CaloL1TopoMonitorAlgorithm.withDetails

Definition at line 224 of file L1CaloL1TopoMonitorAlgorithm.py.


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