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

Class for Tile Jet based monitoring. More...

#include <TileJetMonitorAlgorithm.h>

Inheritance diagram for TileJetMonitorAlgorithm:
Collaboration diagram for TileJetMonitorAlgorithm:

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

 TileJetMonitorAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~TileJetMonitorAlgorithm ()
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

 flags = initConfigFlags()
 Files
 HISTFileName
 useTrigger
 enableLumiAccess
 MaxEvents
 cfg = MainServicesCfg(flags)
 tileJetMonitorAccumulator
 sc = cfg.run()

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

typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > MonVarVec_t
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

StatusCode fillTimeHistograms (const xAOD::Jet &jet, uint32_t lumiBlock, std::set< Identifier > &usedCells) const
unsigned int findIndex (const int gain, const float energy) const
bool isGoodChannel (int part, int module, int channel, uint32_t bad, unsigned int qbit, Identifier id) const
bool passesJvt (const xAOD::Jet &jet) const
bool isGoodJet (const xAOD::Jet &jet) const
bool isGoodEvent (const EventContext &ctx) const
std::string sampleName (const int ros, const int sample, const int tower) const
std::string cellName (const int ros, const int sample, const int tower, const int module) const
bool matchesEnergyRange (const int sample, const int tower, const float energy, const int gain) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

Gaudi::Property< float > m_jetPtMin {this, "JetPtMin", 20000., "Threshold in MeV"}
Gaudi::Property< float > m_jetPtMax {this, "JetPtMax", 120000, ""}
Gaudi::Property< float > m_jetEtaMax {this, "JetEtaMax", 1.6, ""}
Gaudi::Property< float > m_energyChanMin {this, "ChannelEnergyMin", 2000, ""}
Gaudi::Property< float > m_energyChanMax {this, "ChannelEnergyMax", 4000, ""}
Gaudi::Property< float > m_gain {this, "Gain", 1, ""}
Gaudi::Property< float > m_energyE1Min {this, "E1EnergyMin", -1, ""}
Gaudi::Property< float > m_energyE1Max {this, "E1EnergyMax", -1, ""}
Gaudi::Property< float > m_gainE1 {this, "GainE1", -1, ""}
Gaudi::Property< float > m_energyE2Min {this, "E2EnergyMin", -1, ""}
Gaudi::Property< float > m_energyE2Max {this, "E2EnergyMax", -1, ""}
Gaudi::Property< float > m_gainE2 {this, "GainE2", -1, ""}
Gaudi::Property< float > m_energyE3Min {this, "E3EnergyMin", -1, ""}
Gaudi::Property< float > m_energyE3Max {this, "E3EnergyMax", -1, ""}
Gaudi::Property< float > m_gainE3 {this, "GainE3", -1, ""}
Gaudi::Property< float > m_energyE4Min {this, "E4EnergyMin", -1, ""}
Gaudi::Property< float > m_energyE4Max {this, "E4EnergyMax", -1, ""}
Gaudi::Property< float > m_gainE4 {this, "GainE4", -1, ""}
Gaudi::Property< bool > m_do1DHistograms {this, "Do1DHistograms", false, ""}
Gaudi::Property< bool > m_do2DHistograms {this, "Do2DHistograms", false, ""}
Gaudi::Property< bool > m_doCellHistograms {this, "DoCellHistograms", false, ""}
Gaudi::Property< bool > m_doEnergyDiffHistograms {this, "DoEnergyDiffHistograms", false, ""}
Gaudi::Property< float > m_energyDiffThreshold {this, "EnergyDiffThreshold", 2000, ""}
Gaudi::Property< bool > m_doEnergyProfiles {this, "DoEnergyProfiles", true, ""}
Gaudi::Property< bool > m_doEventCleaning {this, "DoEventCleaning", true, ""}
Gaudi::Property< bool > m_doJetCleaning {this, "DoJetCleaning", false, ""}
Gaudi::Property< float > m_jetTrackingEtaLimit {this, "JetTrackingEtaLimit", 2.4, ""}
Gaudi::Property< float > m_jvtThreshold {this, "JvtThreshold", 0.59, ""}
Gaudi::Property< std::vector< float > > m_cellEnergyUpperLimitsHG
Gaudi::Property< std::vector< float > > m_cellEnergyUpperLimitsLG
ToolHandle< ITileBadChanToolm_tileBadChanTool
ToolHandle< IJetUpdateJvtm_jvt {this, "JVT", "", ""}
ToolHandle< IJetSelectorm_jetCleaningTool {this, "JetCleaningTool", "", ""}
ToolHandle< ECUtils::IEventCleaningToolm_eventCleaningTool {this, "EventCleaningTool", "", ""}
SG::ReadHandleKey< xAOD::JetContainerm_jetContainerKey
SG::ReadHandleKey< CaloCellContainerm_caloCellContainerKey
const TileIDm_tileID
const TileHWIDm_tileHWID
const TileCablingServicem_cabling
 TileCabling instance.
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

Detailed Description

Class for Tile Jet based monitoring.

Definition at line 33 of file TileJetMonitorAlgorithm.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 };

Constructor & Destructor Documentation

◆ TileJetMonitorAlgorithm()

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

Definition at line 22 of file TileJetMonitorAlgorithm.cxx.

23 : AthMonitorAlgorithm(name, pSvcLocator)
24 , m_tileID{nullptr}
25 , m_tileHWID{nullptr}
26 , m_cabling{nullptr}
27{}
Base class for Athena Monitoring Algorithms.
const TileCablingService * m_cabling
TileCabling instance.

◆ ~TileJetMonitorAlgorithm()

TileJetMonitorAlgorithm::~TileJetMonitorAlgorithm ( )
virtual

Definition at line 30 of file TileJetMonitorAlgorithm.cxx.

30{}

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}

◆ cellName()

std::string TileJetMonitorAlgorithm::cellName ( const int ros,
const int sample,
const int tower,
const int module ) const
private

Definition at line 479 of file TileJetMonitorAlgorithm.cxx.

479 {
480/*---------------------------------------------------------*/
481 // names compatible with names in TileMonitoringCfgHelper.py
482 std::array<std::string, 3> sample_Name_LB{"A", "B", "D"};
483 std::array<std::string, 4> sample_Name_EB{"A", "B", "D", "E"};
484 std::string c_name;
485 if (ros < 3) { // LBA, LBC
486 c_name = sample_Name_LB[sample] + std::to_string(tower + 1); // default, standard cells
487 if ((sample == TileID::SAMP_BC) && (tower == 8)) { // cell B9
488 c_name = "B9";
489 }
490 if ((sample == TileID::SAMP_D) && (tower == 0)) { // cell D0
491 c_name = "D0";
492 }
493 if ((sample == TileID::SAMP_D) && (tower == 2)) { // cell D1
494 c_name = "D1";
495 }
496 if ((sample == TileID::SAMP_D) && (tower == 4)) { // cell D2
497 c_name = "D2";
498 }
499 if ((sample == TileID::SAMP_D) && (tower == 6)) { // cell D3
500 c_name = "D3";
501 }
502 } else { // EBA, EBC
503 c_name = sample_Name_EB[sample] + std::to_string(tower + 1); // default, standard cells
504 if ((sample == TileID::SAMP_C) && (tower == 9)) { // cell C10
505 c_name = "C10";
506 }
507 if ((sample == TileID::SAMP_D) && (tower == 8)) { // cell D4
508 c_name = "D4";
509 }
510 if ((sample == TileID::SAMP_D) && (tower == 10)) { // cell D5
511 c_name = "D5";
512 }
513 if ((sample == TileID::SAMP_D) && (tower == 12)) { // cell D6
514 c_name = "D6";
515 }
516 if ((sample == TileID::SAMP_E) && (tower == 10)) { // cell E1
517 c_name = "E1";
518 }
519 if ((sample == TileID::SAMP_E) && (tower == 11)) { // cell E2
520 c_name = "E2";
521 }
522 if ((sample == TileID::SAMP_E) && (tower == 13)) { // cell E3
523 if (((ros == 3) && ((module + 1) == 15)) || ((ros == 4) && ((module + 1) == 18))) { // special modules EBA15 & EBC18
524 c_name = "E3*";
525 } else {
526 c_name = "E3";
527 }
528 }
529 if ((sample == TileID::SAMP_E) && (tower == 15)) { // cell E4
530 if (((ros == 3) && ((module + 1) == 15)) || ((ros == 4) && ((module + 1) == 18))) { // special modules EBA15 & EBC18
531 c_name = "E4*";
532 } else {
533 c_name = "E4";
534 }
535 }
536 }
537 return c_name;
538}

◆ 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 TileJetMonitorAlgorithm::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 99 of file TileJetMonitorAlgorithm.cxx.

99 {
100
101 // In case you want to measure the execution time
102 auto timer = Monitored::Timer("TIME_execute");
103
104 if (!isGoodEvent(ctx)) {
105 ATH_MSG_DEBUG("::fillHistograms(), event skipped ");
106 return StatusCode::SUCCESS;
107 } else {
108 ATH_MSG_DEBUG("::fillHistograms(), event accepted ");
109 }
110
111 SG::ReadHandle<xAOD::JetContainer> jetContainer(m_jetContainerKey, ctx);
112 if (!jetContainer.isValid()) {
113 ATH_MSG_WARNING("Can't retrieve Jet Container: " << m_jetContainerKey.key());
114 return StatusCode::SUCCESS;
115 }
116
117 uint32_t lumiBlock = GetEventInfo(ctx)->lumiBlock();
118
119 ATH_MSG_VERBOSE("::fillHistograms(), lumiblock " << lumiBlock);
120
121 std::set<Identifier> usedCells; // cells already used in the given event
122
123 int iJet = 0;
124 for (const xAOD::Jet* jet : *jetContainer) {
125 if ((jet->pt() > m_jetPtMin) && (fabs(jet->eta()) < m_jetEtaMax)) {
126 if (isGoodJet(*jet)) {
127 ATH_MSG_DEBUG("::fillHistograms, jet " << iJet
128 << ", eta " << jet->eta()
129 << ", phi " << jet->phi()
130 << ", constituents " << jet->numConstituents());
131 CHECK(fillTimeHistograms(*jet, lumiBlock, usedCells));
132 } else {
133 ATH_MSG_DEBUG("::fillHistogram, BAD jet " << iJet
134 << ", eta " << jet->eta()
135 << ", phi " << jet->phi()
136 << ", constituents " << jet->numConstituents());
137 }
138 }
139 ++iJet;
140 }
141
142
143 fill("TileJetMonExecuteTime", timer);
144
145 ATH_MSG_VERBOSE("::fillHistograms(), end-of-loops ");
146
147 return StatusCode::SUCCESS;
148}
#define ATH_MSG_VERBOSE(x)
#define CHECK(...)
Evaluate an expression and check for errors.
SG::ReadHandle< xAOD::EventInfo > GetEventInfo(const EventContext &) const
Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
A monitored timer.
StatusCode fillTimeHistograms(const xAOD::Jet &jet, uint32_t lumiBlock, std::set< Identifier > &usedCells) const
SG::ReadHandleKey< xAOD::JetContainer > m_jetContainerKey
bool isGoodJet(const xAOD::Jet &jet) const
Jet_v1 Jet
Definition of the current "jet version".
setTeId lumiBlock
setEventNumber uint32_t
void fill(H5::Group &out_file, size_t iterations)

◆ fillTimeHistograms()

StatusCode TileJetMonitorAlgorithm::fillTimeHistograms ( const xAOD::Jet & jet,
uint32_t lumiBlock,
std::set< Identifier > & usedCells ) const
private

Definition at line 152 of file TileJetMonitorAlgorithm.cxx.

152 {
153/*---------------------------------------------------------*/
154
155 ATH_MSG_VERBOSE( "in fillTimeHistograms()" );
156
157 if( jet.numConstituents() == 0 || !jet.getConstituents().isValid()) return StatusCode::SUCCESS;
158
159 int cellIndex(-1);
160
161 ToolHandle<GenericMonitoringTool> tileJetChannTimeDQTool = getGroup("TileJetChanTimeDQ");
162
163 std::array<std::string, 2> gainName{"LG", "HG"};
164 std::array<std::string, 5> partitionName{"AUX", "LBA", "LBC", "EBA", "EBC"};
165
166 for (const xAOD::JetConstituent* jet_constituent : jet.getConstituents()) {
167 if( jet_constituent->type() == xAOD::Type::CaloCluster ){
168 const xAOD::CaloCluster* calo_cluster = static_cast<const xAOD::CaloCluster*>(jet_constituent->rawConstituent());
169 if (calo_cluster && calo_cluster->getCellLinks()) {
170 for (const CaloCell* cell : *calo_cluster) {
171 cellIndex++;
172 if (cell->caloDDE()->getSubCalo() == CaloCell_ID::TILE) { // a Tile Cell
173 ATH_MSG_DEBUG("Cell " << cellIndex << " IS TILECAL !!");
174 const TileCell *tilecell = static_cast<const TileCell*> (cell);
175 Identifier id = tilecell->ID();
176 if (usedCells.find(id) == usedCells.end()) {
177 usedCells.insert(id);
178 } else {
179 continue;
180 }
181 // int section= m_tileID->section(id);
182 // int module = m_tileID->module(id); // ranges 0..63
183 auto module = Monitored::Scalar<int>("module", m_tileID->module(id));
184 int sample = m_tileID->sample(id); // ranges 0..3 (A, BC, D, E)
185 int tower = m_tileID->tower(id);
186 int ros1 = 0;
187 int ros2 = 0;
188 int chan1 = -1;
189 int chan2 = -1;
190 uint32_t bad1 = 0;
191 uint32_t bad2 = 0;
192 int gain1 = tilecell->gain1();
193 int gain2 = tilecell->gain2();
194 unsigned int qbit1 = tilecell->qbit1();
195 unsigned int qbit2 = tilecell->qbit2();
196
197 const CaloDetDescrElement * caloDDE = tilecell->caloDDE();
198 IdentifierHash hash1 = caloDDE->onl1();
199 if (hash1 != TileHWID::NOT_VALID_HASH) {
200 HWIdentifier adc_id = m_tileHWID->adc_id(hash1, gain1);
201 ros1 = m_tileHWID->ros(adc_id);
202 chan1 = m_tileHWID->channel(adc_id);
203 bad1 = m_tileBadChanTool->encodeStatus(m_tileBadChanTool->getAdcStatus(adc_id));
204 }
205
206 // How is it here with partition? D0 spans two partitions....
207 // It should be ok to treat it in this way:
208 IdentifierHash hash2 = caloDDE->onl2();
209 if (hash2 != TileHWID::NOT_VALID_HASH) {
210 HWIdentifier adc_id = m_tileHWID->adc_id(hash2, gain2);
211 ros2 = m_tileHWID->ros(adc_id);
212 chan2 = m_tileHWID->channel(adc_id);
213 bad2 = m_tileBadChanTool->encodeStatus(m_tileBadChanTool->getAdcStatus(adc_id));
214 }
215
216 bool is_good1 = isGoodChannel(ros1, module, chan1, bad1, qbit1, id);
217 bool is_good2 = isGoodChannel(ros2, module, chan2, bad2, qbit2, id);
218 float ene1 = is_good1 ? tilecell->ene1() : -1;
219 float ene2 = is_good2 ? tilecell->ene2() : -1;
220
221 ATH_MSG_DEBUG(".... " << TileCalibUtils::getDrawerString(ros1, module)
222 << ", ch1 " << chan1
223 << ", ch2 " << chan2
224 << ", qbit " << qbit1 << "/" << qbit2
225 << ", is_bad " << bad1 << "/" << bad2
226 << ", isGood " << is_good1
227 << "/" << is_good2
228 << ", ene " << tilecell->energy());
229 /*
230 Now really fill the histograms time vs lumiblock and 1dim time
231 */
232
233 // first channel
234 if (is_good1 && matchesEnergyRange(sample, tower, ene1, gain1)) {
235 if (m_do1DHistograms) {
236 std::string name = TileCalibUtils::getDrawerString(ros1, module) + "_ch_" + std::to_string(chan1) + "_1d";
237 auto channelTime = Monitored::Scalar<float>(std::move(name), tilecell->time1());
238 fill("TileJetChanTime1D", channelTime);
239 }
240
241 if (m_do2DHistograms) {
242 ATH_MSG_WARNING("These histograms are not implemented yet!");
243 }
244
245 // info for DQ histograms
246 auto moduleDQ = Monitored::Scalar<int>("module" + partitionName[ros1], module + 1);
247 auto channelDQ = Monitored::Scalar<int>("channel" + partitionName[ros1], chan1);
248 auto timeDQ = Monitored::Scalar<float>("time" + partitionName[ros1], tilecell->time1());
249 Monitored::fill(tileJetChannTimeDQTool, moduleDQ, channelDQ, timeDQ);
250
251 // general histograms, only require non-affected channels
252 if (bad1 < 2) {
253 ATH_MSG_DEBUG( "Filling in time1 for " << TileCalibUtils::getDrawerString(ros1, module)
254 << ", ch " << chan1
255 << ", ene " << ene1
256 << ", LB " << lumiBlock
257 << ", time: " << tilecell->time1());
258
259 auto channelTime = Monitored::Scalar<float>("channelTime" + partitionName[ros1], tilecell->time1());
260 fill("TileJetChanTime", channelTime);
261
262 if ((ros1 > 2) && (sample < TileID::SAMP_E)) {
263 std::string nameNoScint("channelTime" + partitionName[ros1] + "_NoScint");
264 auto channelTimeNoScint = Monitored::Scalar<float>(std::move(nameNoScint), tilecell->time1());
265 fill("TileJetChanTime", channelTimeNoScint);
266 }
267 }
268 }
269
270 // second channel
271 if (is_good2 && matchesEnergyRange(sample, tower, ene2, gain2)) {
272 if (m_do1DHistograms) {
273 std::string name = TileCalibUtils::getDrawerString(ros2, module) + "_ch_" + std::to_string(chan2) + "_1d";
274 auto channelTime = Monitored::Scalar<float>(std::move(name), tilecell->time2());
275 fill("TileJetChanTime1D", channelTime);
276 }
277
278 if (m_do2DHistograms) {
279 ATH_MSG_WARNING("This histograms are not implemented yet!");
280 }
281
282 // info for DQ histograms
283 auto moduleDQ = Monitored::Scalar<int>("module" + partitionName[ros2], module + 1);
284 auto channelDQ = Monitored::Scalar<int>("channel" + partitionName[ros2], chan2);
285 auto timeDQ = Monitored::Scalar<float>("time" + partitionName[ros2], tilecell->time2());
286 Monitored::fill(tileJetChannTimeDQTool, moduleDQ, channelDQ, timeDQ);
287
288 // general histograms, only require non-affected channels
289 if (bad2 < 2) {
290 ATH_MSG_DEBUG( "Filling in time2 for " << TileCalibUtils::getDrawerString(ros2, module)
291 << ", ch " << chan2
292 << ", ene " << ene2
293 << ", LB " << lumiBlock
294 << ", time: " << tilecell->time2()
295 <<" (qbit2 " << qbit2 << ", ch1 " << chan1 << ", ene1 " << ene1 << ", bad1 " << bad1 << ", qbit1 " << qbit1 << ")" );
296
297 auto channelTime = Monitored::Scalar<float>("channelTime" + partitionName[ros2], tilecell->time2());
298 fill("TileJetChanTime", channelTime);
299
300 if ((ros2 > 2) && (sample < TileID::SAMP_E)) {
301 std::string nameNoScint("channelTime" + partitionName[ros2] + "_NoScint");
302 auto channelTimeNoScint = Monitored::Scalar<float>(std::move(nameNoScint), tilecell->time2());
303 fill("TileJetChanTime", channelTimeNoScint);
304 }
305 }
306 }
307
308 /*
309 Now filling the cell-based histograms,
310 HG-HG and LG-LG combinations only for normal cells,
311 also include E-cells
312 */
313 if ((is_good1) && (((is_good2) && (gain1 == gain2)) || (sample == TileID::SAMP_E))) {
314 // E-cells are read-out by one channel only, so is_good2 = false for E-cells
315
316 if (m_doEnergyDiffHistograms && (tilecell->energy() > m_energyDiffThreshold)) {
317 // EneDiff histograms
318
319 int evenChannnel = (chan1 % 2 == 0) ? chan1 : chan2;
320 std::string name = TileCalibUtils::getDrawerString(ros1, module) + "_enediff_"
321 + gainName[gain1] + "_ch1_" + std::to_string(evenChannnel);
322 auto energyDifference = Monitored::Scalar<float>(std::move(name), tilecell->eneDiff() / tilecell->energy());
323 fill("TileJetEnergyDiff", energyDifference);
324 }
325
326 if ((bad1 < 2) && (bad2 < 2)) {
327
328 // cell-time histograms, only overall, require not affected channels
329 int index = findIndex(gain1, tilecell->energy());
330 ATH_MSG_DEBUG( "Filling in cell-time for " << TileCalibUtils::getDrawerString(ros1, module)
331 << ", ch1 " << chan1
332 << ", ch2 " << chan2
333 << ", ene " << tilecell->energy()
334 << ", index " << index
335 << ", time: " << tilecell->time());
336
337 // TD adding histograms per partition and per radial sampling
338 std::string name1("Cell_time_" + partitionName[ros1] + "_" + sampleName(ros1, sample, tower) + "_" + gainName[gain1] + "_slice_" + std::to_string(index));
339 auto cellTime1 = Monitored::Scalar<float>(std::move(name1), tilecell->time());
340 fill("TileJetCellTime", cellTime1);
341
342 // Adding histograms per selected individual cell
343 int index_ch1 = -1;
344 int index_ch2 = -1;
345 if (m_doCellHistograms)
346 {
347 std::string name_sel("Cell_time_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_" + gainName[gain1] + "_slice_" + std::to_string(index));
348 auto cellTime_sel = Monitored::Scalar<float>(std::move(name_sel), tilecell->time());
349 fill("TileJetSelCellTime", cellTime_sel);
350
351 // Adding histograms per channels of selected individual cell
352 // First channel
353 index_ch1 = findIndex(gain1, tilecell->ene1() * 2); // Using twice the channel energy to find the correct index
354 std::string name_selCh1("Cell_time_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan1) + "_" + gainName[gain1] + "_slice_" + std::to_string(index_ch1));
355 auto cellTime_selCh1 = Monitored::Scalar<float>(std::move(name_selCh1), tilecell->time1());
356 fill("TileJetSelChanTime", cellTime_selCh1);
357
358 // Second channel
359 index_ch2 = findIndex(gain2, tilecell->ene2() * 2); // Using twice the channel energy to find the correct index
360 std::string name_selCh2("Cell_time_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan2) + "_" + gainName[gain2] + "_slice_" + std::to_string(index_ch2));
361 auto cellTime_selCh2 = Monitored::Scalar<float>(std::move(name_selCh2), tilecell->time2());
362 fill("TileJetSelChanTime", cellTime_selCh2);
363 }
364
365 if (m_doEnergyProfiles) {
366 // TD adding energy profiles per partition and per radial sampling
367 std::string indexName1("index_" + partitionName[ros1] + "_" + sampleName(ros1, sample, tower) + "_" + gainName[gain1]);
368 auto energyIndex1 = Monitored::Scalar<float>(std::move(indexName1), index);
369
370 std::string energyName1("energy_" + partitionName[ros1] + "_" + sampleName(ros1, sample, tower) + "_" + gainName[gain1]);
371 auto cellEnergy1 = Monitored::Scalar<float>(std::move(energyName1), tilecell->energy());
372
373 fill("TileJetCellEnergyProfile", energyIndex1, cellEnergy1);
374
375 // Adding energy profiles per selected individual cell
376 if (m_doCellHistograms)
377 {
378 std::string indexname_sel("index_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_" + gainName[gain1]);
379 auto energyIndex_sel = Monitored::Scalar<float>(std::move(indexname_sel), index);
380
381 std::string energyname_sel("energy_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_" + gainName[gain1]);
382 auto cellEnergy_sel = Monitored::Scalar<float>(std::move(energyname_sel), tilecell->energy());
383
384 fill("TileJetSelCellEnergyProfile", energyIndex_sel, cellEnergy_sel);
385
386 // Adding histograms per channels of selected individual cell
387 // First channel
388 std::string indexname_selCh1("index_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan1) + "_" + gainName[gain1]);
389 auto energyIndex_selCh1 = Monitored::Scalar<float>(std::move(indexname_selCh1), index_ch1);
390
391 std::string energyname_selCh1("energy_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan1) + "_" + gainName[gain1]);
392 auto cellEnergy_selCh1 = Monitored::Scalar<float>(std::move(energyname_selCh1), tilecell->ene1());
393
394 fill("TileJetSelChanEnergyProfile", energyIndex_selCh1, cellEnergy_selCh1);
395
396 // Second channel
397 std::string indexname_selCh2("index_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan2) + "_" + gainName[gain2]);
398 auto energyIndex_selCh2 = Monitored::Scalar<float>(std::move(indexname_selCh2), index_ch2);
399
400 std::string energyname_selCh2("energy_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan2) + "_" + gainName[gain2]);
401 auto cellEnergy_selCh2 = Monitored::Scalar<float>(std::move(energyname_selCh2), tilecell->ene2());
402
403 fill("TileJetSelChanEnergyProfile", energyIndex_selCh2, cellEnergy_selCh2);
404 }
405 } else {
406 // TD adding energy histograms per partition and per radial sampling
407 std::string name1("Cell_ene_" + partitionName[ros1] + "_" + sampleName(ros1, sample, tower) + "_" + gainName[gain1] + "_slice_" + std::to_string(index));
408 auto cellEnergy1 = Monitored::Scalar<float>(std::move(name1), tilecell->energy());
409 fill("TileJetCellEnergy", cellEnergy1);
410
411 // Adding energy histograms per selected individual cell
412 if (m_doCellHistograms)
413 {
414 std::string name_sel("Cell_ene_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_" + gainName[gain1] + "_slice_" + std::to_string(index));
415 auto cellEnergy_sel = Monitored::Scalar<float>(std::move(name_sel), tilecell->energy());
416 fill("TileJetSelCellEnergy", cellEnergy_sel);
417
418 // Adding histograms per channels of selected individual cell
419 // First channel
420 std::string name_selCh1("Cell_ene_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan1) + "_" + gainName[gain1] + "_slice_" + std::to_string(index_ch1));
421 auto cellEnergy_selCh1 = Monitored::Scalar<float>(std::move(name_selCh1), tilecell->ene1());
422 fill("TileJetSelChanEnergy", cellEnergy_selCh1);
423
424 // Second channel
425 std::string name_selCh2("Cell_ene_" + TileCalibUtils::getDrawerString(ros1, module) + "_" + cellName(ros1, sample, tower, module) + "_ch" + std::to_string(chan2) + "_" + gainName[gain2] + "_slice_" + std::to_string(index_ch2));
426 auto cellEnergy_selCh2 = Monitored::Scalar<float>(std::move(name_selCh2), tilecell->ene2());
427 fill("TileJetSelChanEnergy", cellEnergy_selCh2);
428 }
429 }
430 }
431 }
432 } else {
433 ATH_MSG_DEBUG("Cell " << cellIndex << " is NOT Tilecal");
434 }
435 }
436 }
437 }
438 }
439
440 return StatusCode::SUCCESS;
441}
float time() const
get time (data member)
Definition CaloCell.h:368
double energy() const
get energy (data member)
Definition CaloCell.h:327
const CaloDetDescrElement * caloDDE() const
get pointer to CaloDetDescrElement (data member)
Definition CaloCell.h:321
Identifier ID() const
get ID (from cached data member) non-virtual and inline for fast access
Definition CaloCell.h:295
IdentifierHash onl2() const
cell online identifier 2
IdentifierHash onl1() const
cell online identifier 1
static std::string getDrawerString(unsigned int ros, unsigned int drawer)
Return the drawer name, e.g.
float time1(void) const
get time of first PMT
Definition TileCell.h:192
float eneDiff(void) const
all get methods
Definition TileCell.h:182
int gain2(void) const
get gain of second PMT
Definition TileCell.cxx:175
uint8_t qbit2(void) const
get quality bits of second PMT (data member)
Definition TileCell.h:206
int gain1(void) const
get gain of first PMT
Definition TileCell.cxx:168
float ene1(void) const
get energy of first PMT
Definition TileCell.h:187
float time2(void) const
get time of second PMT
Definition TileCell.h:194
uint8_t qbit1(void) const
get quality bits of first PMT (data member)
Definition TileCell.h:203
float ene2(void) const
get energy of second PMT
Definition TileCell.h:189
int channel(const HWIdentifier &id) const
extract channel field from HW identifier
Definition TileHWID.h:189
@ NOT_VALID_HASH
Definition TileHWID.h:314
HWIdentifier adc_id(int ros, int drawer, int channel, int adc) const
adc HWIdentifer
Definition TileHWID.cxx:229
int ros(const HWIdentifier &id) const
extract ros field from HW identifier
Definition TileHWID.h:167
std::string cellName(const int ros, const int sample, const int tower, const int module) const
ToolHandle< ITileBadChanTool > m_tileBadChanTool
bool matchesEnergyRange(const int sample, const int tower, const float energy, const int gain) const
bool isGoodChannel(int part, int module, int channel, uint32_t bad, unsigned int qbit, Identifier id) const
std::string sampleName(const int ros, const int sample, const int tower) const
unsigned int findIndex(const int gain, const float energy) const
int sample(const Identifier &id) const
int tower(const Identifier &id) const
const CaloClusterCellLink * getCellLinks() const
Get a pointer to the CaloClusterCellLink object (const version)
JetConstituentVector getConstituents() const
Return a vector of consituents. The object behaves like vector<const IParticle*>. See JetConstituentV...
Definition Jet_v1.cxx:147
str index
Definition DeMoScan.py:362
void fill(const ToolHandle< GenericMonitoringTool > &tool, T &&... variables)
@ CaloCluster
The object is a calorimeter cluster.
Definition ObjectType.h:39
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.

◆ 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

◆ findIndex()

unsigned int TileJetMonitorAlgorithm::findIndex ( const int gain,
const float energy ) const
private

Definition at line 680 of file TileJetMonitorAlgorithm.cxx.

680 {
681
682 if (gain == 1) {
683 return (std::upper_bound(m_cellEnergyUpperLimitsHG.begin(), m_cellEnergyUpperLimitsHG.end(), energy)
684 - m_cellEnergyUpperLimitsHG.begin());
685 } else {
686 return (std::upper_bound(m_cellEnergyUpperLimitsLG.begin(), m_cellEnergyUpperLimitsLG.end(), energy)
687 - m_cellEnergyUpperLimitsLG.begin());
688 }
689
690}

◆ 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 TileJetMonitorAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 33 of file TileJetMonitorAlgorithm.cxx.

33 {
34
36
37 ATH_MSG_INFO("in initialize()");
38
39 ATH_CHECK( detStore()->retrieve(m_tileID) );
40 ATH_CHECK( detStore()->retrieve(m_tileHWID) );
41
42 //=== get TileCablingSvc
43 ServiceHandle<TileCablingSvc> cablingSvc("TileCablingSvc", name());
44 ATH_CHECK( cablingSvc.retrieve() );
45
46 //=== cache pointers to cabling helpers
47 m_cabling = cablingSvc->cablingService();
48
49 if (!m_cabling) {
50 ATH_MSG_ERROR( "Pointer to TileCablingService is zero: " << m_cabling);
51 return StatusCode::FAILURE;
52 }
53
54 ATH_MSG_INFO("value of m_doJetCleaning: " << m_doJetCleaning);
55
56 //=== get TileBadChanTool
57 ATH_MSG_DEBUG("::Retrieving tile bad channel tool");
58 ATH_CHECK(m_tileBadChanTool.retrieve());
59 ATH_MSG_DEBUG("::Retrieved tile bad channel tool");
60
61 if (m_doJetCleaning) {
62 ATH_MSG_DEBUG("::initializing JVT updater");
63 ATH_CHECK(m_jvt.retrieve());
64 ATH_MSG_DEBUG("::initialized JVT updater");
65
66 ATH_MSG_DEBUG("::initializing JetCleaningTool");
67 ATH_CHECK(m_jetCleaningTool.retrieve());
68 ATH_CHECK(m_eventCleaningTool.retrieve());
69 ATH_MSG_DEBUG("::initialized JetCleaningTool");
70 } else {
71 m_jvt.disable();
72 m_jetCleaningTool.disable();
73 m_eventCleaningTool.disable();
74 }
75
76 ATH_CHECK( m_jetContainerKey.initialize() );
78
79 /* Initialize Gain and Min/Max cell energies for E-cells:
80 if they are not explicitly given, set them as for the ordinary cells
81 */
82 if (m_energyE1Min < 0) m_energyE1Min = m_energyChanMin;
83 if (m_energyE1Max < 0) m_energyE1Max = m_energyChanMax;
84 if (m_gainE1 < 0) m_gainE1 = m_gain;
85 if (m_energyE2Min < 0) m_energyE2Min = m_energyE1Min;
86 if (m_energyE2Max < 0) m_energyE2Max = m_energyE1Max;
87 if (m_gainE2 < 0) m_gainE2 = m_gainE1;
88 if (m_energyE3Min < 0) m_energyE3Min = m_energyE2Min;
89 if (m_energyE3Max < 0) m_energyE3Max = m_energyE2Max;
90 if (m_gainE3 < 0) m_gainE3 = m_gainE2;
91 if (m_energyE4Min < 0) m_energyE4Min = m_energyE3Min;
92 if (m_energyE4Max < 0) m_energyE4Max = m_energyE3Max;
93 if (m_gainE4 < 0) m_gainE4 = m_gainE3;
94
95 return StatusCode::SUCCESS;
96}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_ERROR(x)
#define ATH_MSG_INFO(x)
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.
Gaudi::Property< float > m_gainE2
Gaudi::Property< float > m_gainE4
Gaudi::Property< float > m_gainE1
Gaudi::Property< float > m_energyE3Min
Gaudi::Property< float > m_energyE2Min
SG::ReadHandleKey< CaloCellContainer > m_caloCellContainerKey
Gaudi::Property< float > m_energyE4Max
Gaudi::Property< float > m_energyE1Min
ToolHandle< IJetSelector > m_jetCleaningTool
ToolHandle< IJetUpdateJvt > m_jvt
Gaudi::Property< float > m_energyE2Max
Gaudi::Property< float > m_gainE3
Gaudi::Property< float > m_energyChanMax
Gaudi::Property< float > m_energyE3Max
Gaudi::Property< float > m_gain
Gaudi::Property< float > m_energyE1Max
ToolHandle< ECUtils::IEventCleaningTool > m_eventCleaningTool
Gaudi::Property< float > m_energyChanMin
Gaudi::Property< float > m_energyE4Min

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

◆ isGoodChannel()

bool TileJetMonitorAlgorithm::isGoodChannel ( int part,
int module,
int channel,
uint32_t bad,
unsigned int qbit,
Identifier id ) const
private

Definition at line 571 of file TileJetMonitorAlgorithm.cxx.

571 {
572
573 if ((ros < 1) || (ros >= (int) TileCalibUtils::MAX_ROS)) return false; // invalid partition
574
575 if ((module < 0) || (module >= (int) TileCalibUtils::MAX_DRAWER)) return false; // invalid module number
576
577 if (m_cabling->isDisconnected(ros, module, channel)) {
578 return false; // non-existing PMT (empty hole)
579 }
580
581 if (((qbit & TileCell::MASK_BADCH) != 0) || // masked as bad
582 ((qbit & TileCell::MASK_TIME) != TileCell::MASK_TIME) || // flagged
583 ((qbit & TileCell::MASK_ALGO) == TileFragHash::OptFilterDsp)) // in DSP
584
585 return false;
586 /*
587 bad is the status in the DB (see http://alxr.usatlas.bnl.gov/lxr-stb6/source/atlas/TileCalorimeter/TileConditions/src/TileBadChanTool.cxx#390).
588 Meaning:
589 0 = good, 1 = noisy, 2 = affected, 3 = bad, 4 = otherwise
590 */
591 if (bad > 2) return false;
592
593 /*
594 Now check for special C10, merged E1, E4'
595 C10 spec is ok only if channel = 5 (i.e. pmt=6). The other is pmt=5
596 E1 merged and E4' should be dropped if channel = 12 (i.e. pmt=13)
597 */
598 return ((( channel != 4) && (channel != 12)) || m_cabling->TileGap_connected(id));
599}
bool isDisconnected(int ros, int drawer, int channel) const
bool TileGap_connected(const Identifier &id) const
static const unsigned int MAX_ROS
Number of ROSs.
static const unsigned int MAX_DRAWER
Number of drawers in ROS 1-4.
@ MASK_TIME
Definition TileCell.h:67
@ MASK_ALGO
Definition TileCell.h:62
@ MASK_BADCH
Definition TileCell.h:63

◆ isGoodEvent()

bool TileJetMonitorAlgorithm::isGoodEvent ( const EventContext & ctx) const
private

Definition at line 602 of file TileJetMonitorAlgorithm.cxx.

602 {
603 /* check for errors in LAr and Tile, see https://twiki.cern.ch/twiki/bin/viewauth/Atlas/DataPreparationCheckListForPhysicsAnalysis
604 */
605 if (! m_doEventCleaning) return true;
606
607 ATH_MSG_DEBUG("::isGoodEvent()....");
608
609 SG::ReadHandle<xAOD::EventInfo> eventInfo = GetEventInfo(ctx);
610
611 if (eventInfo->errorState(xAOD::EventInfo::LAr) == xAOD::EventInfo::Error) return false;
612 if (eventInfo->errorState(xAOD::EventInfo::Tile) == xAOD::EventInfo::Error) return false;
613
614 /* see https://twiki.cern.ch/twiki/bin/view/AtlasProtected/HowToCleanJets2017
615 */
616 if (! m_doJetCleaning) return true;
617
618 const xAOD::JetContainer* jetContainer = SG::get(m_jetContainerKey);
619 if (! jetContainer){
620 ATH_MSG_INFO("Cannot retrieve " << m_jetContainerKey << ". However, returning true.");
621 return true;
622 }
623
624 auto jetsSC = xAOD::shallowCopyContainer(*jetContainer);
625 std::unique_ptr< xAOD::JetContainer > jetsCopy(jetsSC.first);
626 std::unique_ptr< xAOD::ShallowAuxContainer > jetsCopyAux(jetsSC.second);
627
628 // We're attaching decorations here to a temporary object that
629 // is not recorded, so we shouldn't use a WriteDecorHandle.
630 static const SG::AuxElement::Decorator<char> passOR ("passOR");
631 static const SG::AuxElement::Decorator<char> passJvt ("passJvt");
632
633 int iJet = 0;
634 for (auto jet : *jetsCopy) {
635 ATH_MSG_DEBUG("Jet " << iJet << ", pT " << jet->pt()/1000.0 << " GeV, eta "
636 << jet->eta());
637 passJvt(*jet) = passesJvt(*jet);
638 passOR(*jet) = true;
639 ATH_MSG_DEBUG("... done with jet " << iJet);
640 ++iJet;
641 }
642
643 bool accept = m_eventCleaningTool->acceptEvent(&*jetsCopy);
644
645 return accept;
646}
SG::Decorator< T, ALLOC > Decorator
Definition AuxElement.h:575
@ Tile
The Tile calorimeter.
@ LAr
The LAr calorimeter.
@ Error
The sub-detector issued an error.
StatusCode accept(const xAOD::Muon *mu)
const T * get(const ReadCondHandleKey< T > &key, const EventContext &ctx)
Convenience function to retrieve an object given a ReadCondHandleKey.
std::pair< std::unique_ptr< T >, std::unique_ptr< ShallowAuxContainer > > shallowCopyContainer(const T &cont, const EventContext &ctx)
Function making a shallow copy of a constant container.
JetContainer_v1 JetContainer
Definition of the current "jet container version".

◆ isGoodJet()

bool TileJetMonitorAlgorithm::isGoodJet ( const xAOD::Jet & jet) const
private

Definition at line 664 of file TileJetMonitorAlgorithm.cxx.

664 {
665
666 if (m_doJetCleaning) {
667
668 if (jet.pt() >= m_jetPtMin && passesJvt(jet) && m_jetCleaningTool->keep(jet)) {
669 return true;
670 } else {
671 return false;
672 }
673
674 } else {
675 return true;
676 }
677
678}

◆ matchesEnergyRange()

bool TileJetMonitorAlgorithm::matchesEnergyRange ( const int sample,
const int tower,
const float energy,
const int gain ) const
private

Definition at line 541 of file TileJetMonitorAlgorithm.cxx.

541 {
542/*---------------------------------------------------------*/
543 /* Want to separate E-cells, D4 and C10 from the rest:
544 ros: 1-2 LBA/LBC, 3-4 EBA/EBC
545 sample: 0 = A, 1 = B/BC/C, D = 2, E = 3
546 tower: 10 = E1, 11 = E2, 13 = E3, 14 = E4, 8 = D4, 9 = C10
547 Nevertheless, C10 and D4 have the same limits, since these are ordinary cells
548 */
549 if (sample != TileID::SAMP_E) {
550 return((energy > m_energyChanMin) && (energy < m_energyChanMax) && (gain == m_gain));
551 } else {
552 switch (tower) {
553 case 10:
554 return((energy > m_energyE1Min) && (energy < m_energyE1Max) && (gain == m_gainE1));
555 break;
556 case 11:
557 return((energy > m_energyE2Min) && (energy < m_energyE2Max) && (gain == m_gainE2));
558 break;
559 case 13:
560 return((energy > m_energyE3Min) && (energy < m_energyE3Max) && (gain == m_gainE3));
561 break;
562 case 15:
563 return((energy > m_energyE4Min) && (energy < m_energyE4Max) && (gain == m_gainE4));
564 break;
565 default:
566 return false;
567 }
568 }
569}

◆ 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

◆ passesJvt()

bool TileJetMonitorAlgorithm::passesJvt ( const xAOD::Jet & jet) const
private

Definition at line 649 of file TileJetMonitorAlgorithm.cxx.

649 {
650
651 if (jet.pt() > m_jetPtMin
652 && jet.pt() < m_jetPtMax
653 && fabs(jet.getAttribute<float>("DetectorEta")) < m_jetTrackingEtaLimit
654 && m_jvt->updateJvt(jet) < m_jvtThreshold) {
655
656 return false;
657
658 } else {
659 return true;
660 }
661
662}

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

◆ sampleName()

std::string TileJetMonitorAlgorithm::sampleName ( const int ros,
const int sample,
const int tower ) const
private

Definition at line 444 of file TileJetMonitorAlgorithm.cxx.

444 {
445/*---------------------------------------------------------*/
446 std::array<std::string, 3> sample_Name_LB{"A", "B", "D"};
447 std::array<std::string, 4> sample_Name_EB{"A", "B", "D", "E"};
448 std::string s_name;
449 if (ros < 3) { // LBA, LBC
450 s_name = sample_Name_LB[sample]; // default, standard cells
451 if ((sample == TileID::SAMP_BC) && (tower == 8)) { // cell B9
452 s_name = "B9";
453 }
454 } else { // EBA, EBC
455 s_name = sample_Name_EB[sample]; //default, standard cells
456 if ((sample == TileID::SAMP_C) && (tower == 9)) { // cell C10
457 s_name = "C10";
458 }
459 if ((sample == TileID::SAMP_D) && (tower == 8)) { // cell D4
460 s_name = "D4";
461 }
462 if ((sample == TileID::SAMP_E) && (tower == 10)) { // cell E1
463 s_name = "E1";
464 }
465 if ((sample == TileID::SAMP_E) && (tower == 11)) { // cell E2
466 s_name = "E2";
467 }
468 if ((sample == TileID::SAMP_E) && (tower == 13)) { // cell E3
469 s_name = "E3";
470 }
471 if ((sample == TileID::SAMP_E) && (tower == 15)) { // cell E4
472 s_name = "E4";
473 }
474 }
475 return s_name;
476}
std::atomic< const char * > ClassID_traits< ViewVector< DV > >::s_name
Definition ViewVector.h:339

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

TileJetMonitorAlgorithm.cfg = MainServicesCfg(flags)

Definition at line 400 of file TileJetMonitorAlgorithm.py.

◆ enableLumiAccess

TileJetMonitorAlgorithm.enableLumiAccess

Definition at line 392 of file TileJetMonitorAlgorithm.py.

◆ Files

TileJetMonitorAlgorithm.Files

Definition at line 389 of file TileJetMonitorAlgorithm.py.

◆ flags

TileJetMonitorAlgorithm.flags = initConfigFlags()

Definition at line 388 of file TileJetMonitorAlgorithm.py.

◆ HISTFileName

TileJetMonitorAlgorithm.HISTFileName

Definition at line 390 of file TileJetMonitorAlgorithm.py.

◆ m_cabling

const TileCablingService* TileJetMonitorAlgorithm::m_cabling
private

TileCabling instance.

Definition at line 108 of file TileJetMonitorAlgorithm.h.

◆ m_caloCellContainerKey

SG::ReadHandleKey<CaloCellContainer> TileJetMonitorAlgorithm::m_caloCellContainerKey
private
Initial value:
{this,
"CaloCellContainer", "AllCalo", "Calo cell container name"}

Definition at line 102 of file TileJetMonitorAlgorithm.h.

102 {this,
103 "CaloCellContainer", "AllCalo", "Calo cell container name"};

◆ m_cellEnergyUpperLimitsHG

Gaudi::Property<std::vector<float> > TileJetMonitorAlgorithm::m_cellEnergyUpperLimitsHG
private
Initial value:
{this,
"CellEnergyUpperLimitsHG", {}, "Energy upper limits of the HG cell-time histograms"}

Definition at line 83 of file TileJetMonitorAlgorithm.h.

83 {this,
84 "CellEnergyUpperLimitsHG", {}, "Energy upper limits of the HG cell-time histograms"};

◆ m_cellEnergyUpperLimitsLG

Gaudi::Property<std::vector<float> > TileJetMonitorAlgorithm::m_cellEnergyUpperLimitsLG
private
Initial value:
{this,
"CellEnergyUpperLimitsLG", {}, "Energy upper limits of the LG cell-time histograms"}

Definition at line 85 of file TileJetMonitorAlgorithm.h.

85 {this,
86 "CellEnergyUpperLimitsLG", {}, "Energy upper limits of the LG cell-time histograms"};

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

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_do1DHistograms {this, "Do1DHistograms", false, ""}
private

Definition at line 73 of file TileJetMonitorAlgorithm.h.

73{this, "Do1DHistograms", false, ""};

◆ m_do2DHistograms

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_do2DHistograms {this, "Do2DHistograms", false, ""}
private

Definition at line 74 of file TileJetMonitorAlgorithm.h.

74{this, "Do2DHistograms", false, ""};

◆ m_doCellHistograms

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_doCellHistograms {this, "DoCellHistograms", false, ""}
private

Definition at line 75 of file TileJetMonitorAlgorithm.h.

75{this, "DoCellHistograms", false, ""};

◆ m_doEnergyDiffHistograms

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_doEnergyDiffHistograms {this, "DoEnergyDiffHistograms", false, ""}
private

Definition at line 76 of file TileJetMonitorAlgorithm.h.

76{this, "DoEnergyDiffHistograms", false, ""};

◆ m_doEnergyProfiles

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_doEnergyProfiles {this, "DoEnergyProfiles", true, ""}
private

Definition at line 78 of file TileJetMonitorAlgorithm.h.

78{this, "DoEnergyProfiles", true, ""};

◆ m_doEventCleaning

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_doEventCleaning {this, "DoEventCleaning", true, ""}
private

Definition at line 79 of file TileJetMonitorAlgorithm.h.

79{this, "DoEventCleaning", true, ""};

◆ m_doJetCleaning

Gaudi::Property<bool> TileJetMonitorAlgorithm::m_doJetCleaning {this, "DoJetCleaning", false, ""}
private

Definition at line 80 of file TileJetMonitorAlgorithm.h.

80{this, "DoJetCleaning", false, ""};

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

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyChanMax {this, "ChannelEnergyMax", 4000, ""}
private

Definition at line 58 of file TileJetMonitorAlgorithm.h.

58{this, "ChannelEnergyMax", 4000, ""};

◆ m_energyChanMin

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyChanMin {this, "ChannelEnergyMin", 2000, ""}
private

Definition at line 57 of file TileJetMonitorAlgorithm.h.

57{this, "ChannelEnergyMin", 2000, ""};

◆ m_energyDiffThreshold

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyDiffThreshold {this, "EnergyDiffThreshold", 2000, ""}
private

Definition at line 77 of file TileJetMonitorAlgorithm.h.

77{this, "EnergyDiffThreshold", 2000, ""};

◆ m_energyE1Max

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE1Max {this, "E1EnergyMax", -1, ""}
private

Definition at line 61 of file TileJetMonitorAlgorithm.h.

61{this, "E1EnergyMax", -1, ""};

◆ m_energyE1Min

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE1Min {this, "E1EnergyMin", -1, ""}
private

Definition at line 60 of file TileJetMonitorAlgorithm.h.

60{this, "E1EnergyMin", -1, ""};

◆ m_energyE2Max

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE2Max {this, "E2EnergyMax", -1, ""}
private

Definition at line 64 of file TileJetMonitorAlgorithm.h.

64{this, "E2EnergyMax", -1, ""};

◆ m_energyE2Min

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE2Min {this, "E2EnergyMin", -1, ""}
private

Definition at line 63 of file TileJetMonitorAlgorithm.h.

63{this, "E2EnergyMin", -1, ""};

◆ m_energyE3Max

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE3Max {this, "E3EnergyMax", -1, ""}
private

Definition at line 67 of file TileJetMonitorAlgorithm.h.

67{this, "E3EnergyMax", -1, ""};

◆ m_energyE3Min

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE3Min {this, "E3EnergyMin", -1, ""}
private

Definition at line 66 of file TileJetMonitorAlgorithm.h.

66{this, "E3EnergyMin", -1, ""};

◆ m_energyE4Max

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE4Max {this, "E4EnergyMax", -1, ""}
private

Definition at line 70 of file TileJetMonitorAlgorithm.h.

70{this, "E4EnergyMax", -1, ""};

◆ m_energyE4Min

Gaudi::Property<float> TileJetMonitorAlgorithm::m_energyE4Min {this, "E4EnergyMin", -1, ""}
private

Definition at line 69 of file TileJetMonitorAlgorithm.h.

69{this, "E4EnergyMin", -1, ""};

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

ToolHandle<ECUtils::IEventCleaningTool> TileJetMonitorAlgorithm::m_eventCleaningTool {this, "EventCleaningTool", "", ""}
private

Definition at line 96 of file TileJetMonitorAlgorithm.h.

96{this, "EventCleaningTool", "", ""};

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

Gaudi::Property<float> TileJetMonitorAlgorithm::m_gain {this, "Gain", 1, ""}
private

Definition at line 59 of file TileJetMonitorAlgorithm.h.

59{this, "Gain", 1, ""};

◆ m_gainE1

Gaudi::Property<float> TileJetMonitorAlgorithm::m_gainE1 {this, "GainE1", -1, ""}
private

Definition at line 62 of file TileJetMonitorAlgorithm.h.

62{this, "GainE1", -1, ""};

◆ m_gainE2

Gaudi::Property<float> TileJetMonitorAlgorithm::m_gainE2 {this, "GainE2", -1, ""}
private

Definition at line 65 of file TileJetMonitorAlgorithm.h.

65{this, "GainE2", -1, ""};

◆ m_gainE3

Gaudi::Property<float> TileJetMonitorAlgorithm::m_gainE3 {this, "GainE3", -1, ""}
private

Definition at line 68 of file TileJetMonitorAlgorithm.h.

68{this, "GainE3", -1, ""};

◆ m_gainE4

Gaudi::Property<float> TileJetMonitorAlgorithm::m_gainE4 {this, "GainE4", -1, ""}
private

Definition at line 71 of file TileJetMonitorAlgorithm.h.

71{this, "GainE4", -1, ""};

◆ m_jetCleaningTool

ToolHandle<IJetSelector> TileJetMonitorAlgorithm::m_jetCleaningTool {this, "JetCleaningTool", "", ""}
private

Definition at line 95 of file TileJetMonitorAlgorithm.h.

95{this, "JetCleaningTool", "", ""};

◆ m_jetContainerKey

SG::ReadHandleKey<xAOD::JetContainer> TileJetMonitorAlgorithm::m_jetContainerKey
private
Initial value:
{this,
"JetContainer", "AntiKt4EMTopoJets", "Jet container for monitoring"}

Definition at line 99 of file TileJetMonitorAlgorithm.h.

99 {this,
100 "JetContainer", "AntiKt4EMTopoJets", "Jet container for monitoring"};

◆ m_jetEtaMax

Gaudi::Property<float> TileJetMonitorAlgorithm::m_jetEtaMax {this, "JetEtaMax", 1.6, ""}
private

Definition at line 56 of file TileJetMonitorAlgorithm.h.

56{this, "JetEtaMax", 1.6, ""};

◆ m_jetPtMax

Gaudi::Property<float> TileJetMonitorAlgorithm::m_jetPtMax {this, "JetPtMax", 120000, ""}
private

Definition at line 55 of file TileJetMonitorAlgorithm.h.

55{this, "JetPtMax", 120000, ""};

◆ m_jetPtMin

Gaudi::Property<float> TileJetMonitorAlgorithm::m_jetPtMin {this, "JetPtMin", 20000., "Threshold in MeV"}
private

Definition at line 54 of file TileJetMonitorAlgorithm.h.

54{this, "JetPtMin", 20000., "Threshold in MeV"};

◆ m_jetTrackingEtaLimit

Gaudi::Property<float> TileJetMonitorAlgorithm::m_jetTrackingEtaLimit {this, "JetTrackingEtaLimit", 2.4, ""}
private

Definition at line 81 of file TileJetMonitorAlgorithm.h.

81{this, "JetTrackingEtaLimit", 2.4, ""};

◆ m_jvt

ToolHandle<IJetUpdateJvt> TileJetMonitorAlgorithm::m_jvt {this, "JVT", "", ""}
private

Definition at line 92 of file TileJetMonitorAlgorithm.h.

92{this, "JVT", "", ""};

◆ m_jvtThreshold

Gaudi::Property<float> TileJetMonitorAlgorithm::m_jvtThreshold {this, "JvtThreshold", 0.59, ""}
private

Definition at line 82 of file TileJetMonitorAlgorithm.h.

82{this, "JvtThreshold", 0.59, ""};

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

ToolHandle<ITileBadChanTool> TileJetMonitorAlgorithm::m_tileBadChanTool
private
Initial value:
{this,
"TileBadChanTool", "TileBadChanTool", "Tile bad channel tool"}

Definition at line 88 of file TileJetMonitorAlgorithm.h.

88 {this,
89 "TileBadChanTool", "TileBadChanTool", "Tile bad channel tool"};

◆ m_tileHWID

const TileHWID* TileJetMonitorAlgorithm::m_tileHWID
private

Definition at line 106 of file TileJetMonitorAlgorithm.h.

◆ m_tileID

const TileID* TileJetMonitorAlgorithm::m_tileID
private

Definition at line 105 of file TileJetMonitorAlgorithm.h.

◆ m_toolLookupMap

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

Definition at line 372 of file AthMonitorAlgorithm.h.

◆ m_tools

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

Array of Generic Monitoring Tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

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.

◆ MaxEvents

TileJetMonitorAlgorithm.MaxEvents

Definition at line 393 of file TileJetMonitorAlgorithm.py.

◆ sc

TileJetMonitorAlgorithm.sc = cfg.run()

Definition at line 413 of file TileJetMonitorAlgorithm.py.

◆ tileJetMonitorAccumulator

TileJetMonitorAlgorithm.tileJetMonitorAccumulator
Initial value:
= TileJetMonitoringConfig(flags,
Do1DHistograms = True,
DoEnergyDiffHistograms = True,
DoCellHistograms = False)

Definition at line 403 of file TileJetMonitorAlgorithm.py.

◆ useTrigger

TileJetMonitorAlgorithm.useTrigger

Definition at line 391 of file TileJetMonitorAlgorithm.py.


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