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

#include <TrigTauMonitorSingleAlgorithm.h>

Inheritance diagram for TrigTauMonitorSingleAlgorithm:
Collaboration diagram for TrigTauMonitorSingleAlgorithm:

Public Types

enum  TauID { None = 0 , RNN = 1 , GNTau = 2 }
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

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

Protected Member Functions

std::map< std::string, TrigTauInfo > & getTrigInfoMap ()
const TrigTauInfogetTrigInfo (const std::string &trigger) const
std::vector< const xAOD::TauJet * > getOnlineTausAll (const std::string &trigger, bool include_0P=true) const
std::tuple< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > getOnlineTaus (const std::string &trigger) const
std::vector< const xAOD::TauJet * > getOfflineTausAll (const EventContext &ctx, const float threshold=20.0) const
std::pair< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > getOfflineTaus (const EventContext &ctx, const float threshold=20.0, const TauID tau_id=TauID::None) const
std::vector< const xAOD::eFexTauRoI * > getL1eTAUs (const EventContext &ctx, const std::string &l1_item) const
std::vector< const xAOD::jFexTauRoI * > getL1jTAUs (const EventContext &ctx, const std::string &l1_item) const
std::vector< std::pair< const xAOD::eFexTauRoI *, const xAOD::jFexTauRoI * > > getL1cTAUs (const EventContext &ctx, const std::string &l1_item) const
double dR (const double eta1, const double phi1, const double eta2, const double phi2) const
template<typename T1 = xAOD::IParticle, typename T2 = xAOD::IParticle>
bool matchObjects (const T1 *tau, const std::vector< const T2 * > &tau_vec, float threshold) const
bool matchObjects (const TLorentzVector &tau, const std::vector< TLorentzVector > &tau_vec, float threshold) const
template<typename T1 = xAOD::IParticle, typename T2 = xAOD::eFexTauRoI>
bool matchObjects (const T1 *tau_1, const T2 *tau_2, float threshold) const
template<typename T1 = xAOD::IParticle, typename T2 = xAOD::IParticle>
bool matchTruthObjects (const T1 *true_tau, const std::vector< const T2 * > &tau_vec, float threshold) const
std::vector< const xAOD::TauJet * > classifyTausAll (const std::vector< const xAOD::TauJet * > &taus, const float threshold=0.0, const TauID tau_id=TauID::None) const
std::pair< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > classifyOfflineTaus (const std::vector< const xAOD::TauJet * > &taus, const float threshold=0.0, const TauID tau_id=TauID::None) const
std::tuple< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > classifyOnlineTaus (const std::vector< const xAOD::TauJet * > &taus, const float threshold=0.0) const
const SG::ReadHandleKey< xAOD::TauJetContainer > & getOnlineContainerKey (const std::string &sequence) const
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

Gaudi::Property< std::vector< std::string > > m_triggers {this, "TriggerList", {}}
Gaudi::Property< std::map< std::string, float > > m_L1_Phase1_thresholds {this, "L1Phase1Thresholds", {}}
Gaudi::Property< std::map< std::string, uint64_t > > m_L1_Phase1_threshold_patterns {this, "L1Phase1ThresholdPatterns", {}}
Gaudi::Property< bool > m_L1_select_by_et_only {this, "SelectL1ByETOnly", false}
Gaudi::Property< bool > m_do_efficiency_plots {this, "DoEfficiencyPlots", true}
Gaudi::Property< bool > m_do_variable_plots {this, "DoVariablePlots", true}
SG::ReadHandleKey< xAOD::TauJetContainerm_offlineTauJetKey {this, "OfflineTauJetKey", "TauJets", "Offline taujet container key"}
SG::ReadHandleKey< xAOD::TauJetContainerm_hltTauJetKey {this, "HLTTauJetKey", "HLT_TrigTauRecMerged_MVA", "HLT tracktwoMVA taujet container key"}
SG::ReadHandleKey< xAOD::TauJetContainerm_hltTauJetLLPKey {this, "HLTTauJetLLPKey", "HLT_TrigTauRecMerged_LLP", "HLT tracktwoLLP taujet container key"}
SG::ReadHandleKey< xAOD::TauJetContainerm_hltTauJetLRTKey {this, "HLTTauJetLRTKey", "HLT_TrigTauRecMerged_LRT", "HLT trackLRT taujet container key"}
SG::ReadHandleKey< xAOD::TauJetContainerm_hltTauJetCaloMVAOnlyKey {this, "HLTTauJetCaloMVAOnlyKey", "HLT_TrigTauRecMerged_CaloMVAOnly", "HLT ptonly taujet container key"}
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

virtual StatusCode processEvent (const EventContext &ctx) const override
void fillHLTEfficiencies (const EventContext &ctx, const std::string &trigger, const bool l1_accept_flag, const std::vector< const xAOD::TauJet * > &offline_tau_vec, const std::vector< const xAOD::TauJet * > &online_tau_vec, const std::string &nProng) const
void fillIDInputVars (const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, const std::string &nProng, bool online) const
void fillIDTrack (const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, bool online) const
void fillIDCluster (const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, bool online) const
void fillBasicVars (const EventContext &ctx, const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, const std::string &nProng, bool online) const
void fillIDScores (const EventContext &ctx, const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, const std::string &nProng, bool online) const
std::vector< TLorentzVector > getRoIsVector (const EventContext &ctx, const std::string &trigger) 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< bool > m_doTotalEfficiency {this, "DoTotalEfficiency", false, "Do total efficiency histograms"}
Gaudi::Property< bool > m_requireOfflineTaus {this, "RequireOfflineTaus", true, "Require at leat 1 offline tau per event"}
Gaudi::Property< unsigned int > m_offline_tau_id {this, "OfflineTauID", TauID::RNN, "Offline TauID (1: RNN, 2: GNTau)"}
Gaudi::Property< bool > m_doOfflineTausDistributions {this, "DoOfflineTausDistributions", true}
Gaudi::Property< std::map< std::string, std::map< std::string, std::pair< std::string, std::string > > > > m_monitoredHLTIdScores {this, "HLTTauIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each HLT TauID algorithm to be monitored, for each reconstruction sequence (type, e.g. tracktwoMVA, tracktwoLLP, etc...)"}
std::map< std::string, std::map< std::string, std::pair< SG::ReadDecorHandleKey< xAOD::TauJetContainer >, SG::ReadDecorHandleKey< xAOD::TauJetContainer > > > > m_monitoredHLTIdDecorHandleKeys
Gaudi::Property< std::map< std::string, std::pair< std::string, std::string > > > m_monitoredOfflineIdScores {this, "OfflineTauIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each Offline TauID algorithm to be monitored"}
std::map< std::string, std::pair< SG::ReadDecorHandleKey< xAOD::TauJetContainer >, SG::ReadDecorHandleKey< xAOD::TauJetContainer > > > m_monitoredOfflineIdDecorHandleKeys
std::map< std::string, TrigTauInfom_trigInfo
SG::ReadDecorHandleKey< xAOD::EventInfom_eventInfoDecorKey {this, "LArStatusFlag", "EventInfo.larFlags", "Key for EventInfo object"}
SG::ReadDecorHandleKey< xAOD::TauJetContainerm_offlineGNTauDecorKey {this, "OfflineGNTauDecorKey", "", "Offline GNTau decoration key"}
SG::ReadHandleKey< xAOD::eFexTauRoIContainerm_phase1l1eTauRoIKey {this, "Phase1L1eTauRoIKey", "L1_eTauRoI", "eTau Phase1 L1 RoI key"}
SG::ReadDecorHandleKey< xAOD::eFexTauRoIContainerm_phase1l1eTauRoIThresholdPatternsKey {this, "Phase1L1eTauRoIThresholdPatternsKey", "L1_eTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the eTau RoIs"}
SG::ReadHandleKey< xAOD::jFexTauRoIContainerm_phase1l1jTauRoIKey {this, "Phase1L1jTauRoIKey", "L1_jFexTauRoI", "jTau Phase1 L1 RoI key"}
SG::ReadDecorHandleKey< xAOD::jFexTauRoIContainerm_phase1l1jTauRoIThresholdPatternsKey {this, "Phase1L1jTauRoIThresholdPatternsKey", "L1_jFexTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the jTau RoIs"}
SG::ReadHandleKey< xAOD::eFexTauRoIContainerm_phase1l1cTauRoIKey {this, "Phase1L1cTauRoIKey", "L1_cTauRoI", "cTau Phase1 L1 RoI key"}
SG::ReadDecorHandleKey< xAOD::eFexTauRoIContainerm_phase1l1cTauRoIDecorKey {this, "Phase1L1cTauRoIjTauRoILinkKey", "L1_cTauRoI.jTauLink", "Decoration for the link from eTau to the matching jTau"}
SG::ReadDecorHandleKey< xAOD::eFexTauRoIContainerm_phase1l1cTauRoIThresholdPatternsKey {this, "Phase1L1cTauRoIThresholdPatternsKey", "L1_cTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the cTau RoIs"}
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

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

◆ TauID

Constructor & Destructor Documentation

◆ TrigTauMonitorSingleAlgorithm()

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

Definition at line 11 of file TrigTauMonitorSingleAlgorithm.cxx.

12 : TrigTauMonitorBaseAlgorithm(name, pSvcLocator)
13{}
TrigTauMonitorBaseAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)

Member Function Documentation

◆ cardinality()

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

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

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ classifyOfflineTaus()

std::pair< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > TrigTauMonitorBaseAlgorithm::classifyOfflineTaus ( const std::vector< const xAOD::TauJet * > & taus,
const float threshold = 0.0,
const TauID tau_id = TauID::None ) const
protectedinherited

Definition at line 322 of file TrigTauMonitorBaseAlgorithm.cxx.

323{
324 std::vector<const xAOD::TauJet*> tau_vec_1p, tau_vec_3p;
325
326 for(const xAOD::TauJet* const tau : classifyTausAll(taus, threshold, tau_id)) {
327 int nTracks = -1;
328 tau->detail(xAOD::TauJetParameters::nChargedTracks, nTracks);
329
330 if(nTracks == 1) tau_vec_1p.push_back(tau);
331 else if(nTracks == 3) tau_vec_3p.push_back(tau);
332 }
333
334 return {tau_vec_1p, tau_vec_3p};
335}
static Double_t taus
std::vector< const xAOD::TauJet * > classifyTausAll(const std::vector< const xAOD::TauJet * > &taus, const float threshold=0.0, const TauID tau_id=TauID::None) const
TauJet_v3 TauJet
Definition of the current "tau version".

◆ classifyOnlineTaus()

std::tuple< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > TrigTauMonitorBaseAlgorithm::classifyOnlineTaus ( const std::vector< const xAOD::TauJet * > & taus,
const float threshold = 0.0 ) const
protectedinherited

Definition at line 305 of file TrigTauMonitorBaseAlgorithm.cxx.

306{
307 std::vector<const xAOD::TauJet*> tau_vec_0p, tau_vec_1p, tau_vec_mp;
308
309 for(const xAOD::TauJet* tau : classifyTausAll(taus, threshold, TauID::None)) {
310 int nTracks = -1;
311 tau->detail(xAOD::TauJetParameters::nChargedTracks, nTracks);
312
313 if(nTracks == 0) tau_vec_0p.push_back(tau);
314 else if(nTracks == 1) tau_vec_1p.push_back(tau);
315 else tau_vec_mp.push_back(tau);
316 }
317
318 return {tau_vec_0p, tau_vec_1p, tau_vec_mp};
319}

◆ classifyTausAll()

std::vector< const xAOD::TauJet * > TrigTauMonitorBaseAlgorithm::classifyTausAll ( const std::vector< const xAOD::TauJet * > & taus,
const float threshold = 0.0,
const TauID tau_id = TauID::None ) const
protectedinherited

Definition at line 277 of file TrigTauMonitorBaseAlgorithm.cxx.

278{
279 std::vector<const xAOD::TauJet*> tau_vec;
280
281 SG::ReadDecorHandle<xAOD::TauJetContainer, char> tauid_medium{m_offlineGNTauDecorKey, Gaudi::Hive::currentContext()};
282 if(!tauid_medium.isValid()) {
283 ATH_MSG_WARNING("Cannot retrieve " << tauid_medium.key());
284 return tau_vec;
285 }
286
287 for(const xAOD::TauJet* tau : taus) {
288 if(tau->pt() < threshold*Gaudi::Units::GeV) continue;
289
290 // Consider only offline taus which pass medium ID WP
291 if(tau_id == TauID::RNN) {
292 if(!tau->isTau(xAOD::TauJetParameters::JetRNNSigMedium)) continue;
293 }
294 else if(tau_id == TauID::GNTau) {
295 if(!tauid_medium(*tau)) continue;
296 }
297
298 tau_vec.push_back(tau);
299 }
300
301 return tau_vec;
302}
#define ATH_MSG_WARNING(x)
SG::ReadDecorHandleKey< xAOD::TauJetContainer > m_offlineGNTauDecorKey

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

◆ dR()

double TrigTauMonitorBaseAlgorithm::dR ( const double eta1,
const double phi1,
const double eta2,
const double phi2 ) const
inlineprotectedinherited

Definition at line 70 of file TrigTauMonitorBaseAlgorithm.h.

71 {
72 double deta = std::fabs(eta1 - eta2);
73 double dphi = std::fabs(CxxUtils::wrapToPi(phi1-phi2));
74 return std::sqrt(deta*deta + dphi*dphi);
75 }
T wrapToPi(T phi)
Wrap angle in radians to [-pi, pi].
Definition phihelper.h:24

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

◆ fillBasicVars()

void TrigTauMonitorSingleAlgorithm::fillBasicVars ( const EventContext & ctx,
const std::string & trigger,
const std::vector< const xAOD::TauJet * > & tau_vec,
const std::string & nProng,
bool online ) const
private

Definition at line 414 of file TrigTauMonitorSingleAlgorithm.cxx.

415{
416 ATH_MSG_DEBUG("Fill Basic Variables: " << trigger);
417
418 auto monGroup = getGroup(trigger+"_"+(online ? "HLT" : "Offline")+"_basicVars_"+nProng);
419
420 auto Pt = Monitored::Collection("Pt", tau_vec, [](const xAOD::TauJet* tau){ return tau->pt()/Gaudi::Units::GeV; });
421 auto Eta = Monitored::Collection("Eta", tau_vec, [](const xAOD::TauJet* tau){ return tau->eta(); });
422 auto Phi = Monitored::Collection("Phi", tau_vec, [](const xAOD::TauJet* tau){ return tau->phi(); });
423
424 auto nTrack = Monitored::Collection("nTrack", tau_vec, [](const xAOD::TauJet* tau){
425 int nTrack = -1;
427 return nTrack;
428 });
429 auto nIsoTrack = Monitored::Collection("nIsoTrack", tau_vec, [](const xAOD::TauJet* tau){ return tau->nTracksIsolation(); });
430
431 auto averageMu = Monitored::Scalar<float>("averageMu", 0.0);
432 averageMu = lbAverageInteractionsPerCrossing(ctx);
433
434 auto TauVertexX = Monitored::Collection("TauVertexX", tau_vec, [](const xAOD::TauJet* tau){
435 double vtx = -999;
436 if(tau->vertex() != nullptr) vtx = tau->vertex()->x();
437 return vtx;
438 });
439 auto TauVertexY = Monitored::Collection("TauVertexY", tau_vec, [](const xAOD::TauJet* tau){
440 double vty = -999;
441 if(tau->vertex() != nullptr) vty = tau->vertex()->y();
442 return vty;
443 });
444 auto TauVertexZ = Monitored::Collection("TauVertexZ", tau_vec, [](const xAOD::TauJet* tau){
445 double vtz = -999;
446 if(tau->vertex() != nullptr) vtz = tau->vertex()->z();
447 return vtz;
448 });
449
450 fill(monGroup, Pt, Eta, Phi, nTrack, nIsoTrack, averageMu, TauVertexX, TauVertexY, TauVertexZ);
451
452 ATH_MSG_DEBUG("After fill Basic variables: " << trigger);
453}
@ Phi
Definition RPCdef.h:8
@ Eta
Definition RPCdef.h:8
const ToolHandle< GenericMonitoringTool > & getGroup(const std::string &name) const
Get a specific monitoring tool from the tool handle array.
virtual double phi() const
The azimuthal angle ( ) of the particle.
virtual double pt() const
The transverse momentum ( ) of the particle.
bool detail(TauJetParameters::Detail detail, int &value) const
Get and set values of common details variables via enum.
const Vertex * vertex() const
size_t nTracksIsolation() const
virtual double eta() const
The pseudorapidity ( ) of the particle.
float z() const
Returns the z position.
float y() const
Returns the y position.
float x() const
Returns the x position.
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.
virtual float lbAverageInteractionsPerCrossing(const EventContext &ctx=Gaudi::Hive::currentContext()) const
Calculate the average mu, i.e.
ValuesCollection< T > Collection(std::string name, const T &collection)
Declare a monitored (double-convertible) collection.
monGroup(analysis_chain)

◆ fillHistograms()

StatusCode TrigTauMonitorBaseAlgorithm::fillHistograms ( const EventContext & ctx) const
overridevirtualinherited

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 250 of file TrigTauMonitorBaseAlgorithm.cxx.

251{
252 ATH_MSG_DEBUG("Executing Monitoring algorithm");
253
254 // Protect against truncated events
255 // Since this happens very rarely, it won't bias the L1 distributions and efficiencies
256 if(m_trigDecTool->ExperimentalAndExpertMethods().isHLTTruncated()){
257 ATH_MSG_WARNING("HLTResult truncated, skip trigger analysis");
258 return StatusCode::SUCCESS;
259 }
260
261 // Protect against LAr noise bursts and other detector errors
262 SG::ReadHandle<xAOD::EventInfo> eventInfo(GetEventInfo(ctx));
263 ATH_CHECK(eventInfo.isValid());
264 if (eventInfo->errorState(xAOD::EventInfo::LAr) == xAOD::EventInfo::Error
265 || eventInfo->errorState(xAOD::EventInfo::Tile) == xAOD::EventInfo::Error
266 || eventInfo->errorState(xAOD::EventInfo::SCT) == xAOD::EventInfo::Error
267 || eventInfo->isEventFlagBitSet(xAOD::EventInfo::Core, 18)) {
268 return StatusCode::SUCCESS;
269 }
270
272
273 return StatusCode::SUCCESS;
274}
#define ATH_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.)
virtual StatusCode processEvent(const EventContext &ctx) const =0
@ Tile
The Tile calorimeter.
@ Core
Core flags describing the event.
@ LAr
The LAr calorimeter.
@ Error
The sub-detector issued an error.

◆ fillHLTEfficiencies()

void TrigTauMonitorSingleAlgorithm::fillHLTEfficiencies ( const EventContext & ctx,
const std::string & trigger,
const bool l1_accept_flag,
const std::vector< const xAOD::TauJet * > & offline_tau_vec,
const std::vector< const xAOD::TauJet * > & online_tau_vec,
const std::string & nProng ) const
private

Definition at line 153 of file TrigTauMonitorSingleAlgorithm.cxx.

154{
155 ATH_MSG_DEBUG("Fill HLT " << nProng << " efficiencies: " << trigger);
156
157 const TrigTauInfo& info = getTrigInfo(trigger);
158
159 // Efficiency for single leg tau triggers:
160 // denominator = offline tau + matching with L1 object with dR(offline tau,L1 item) < 0.3
161 // numerator = denominator + hlt fires + matching with HLT tau with dR(offline tau, HLT tau) < 0.2
162
163 auto monGroup = getGroup(trigger+"_HLT_Efficiency_"+nProng);
164
165 auto tauPt = Monitored::Scalar<float>("tauPt", 0.0);
166 auto tauEta = Monitored::Scalar<float>("tauEta", 0.0);
167 auto tauPhi = Monitored::Scalar<float>("tauPhi", 0.0);
168 auto averageMu = Monitored::Scalar<float>("averageMu", 0.0);
169 auto HLT_match = Monitored::Scalar<bool>("HLT_pass", false);
170 auto HLT_match_highPt = Monitored::Scalar<bool>("HLT_pass_highPt", false);
171 auto Total_match = Monitored::Scalar<bool>("Total_pass", false);
172 auto Total_match_highPt = Monitored::Scalar<bool>("Total_pass_highPt", false);
173
174 bool hlt_fires = m_trigDecTool->isPassed(trigger, TrigDefs::Physics | TrigDefs::allowResurrectedDecision);
175
176 std::vector<TLorentzVector> rois = getRoIsVector(ctx, trigger);
177 for(const auto *offline_tau : offline_tau_vec) {
178 bool L1_match = false;
179
180 // Check the matching offline tau with L1 item -> depending on the L1 type (phase-1 eTAU, jTAU, cTAU)
181 // All L1 RoIs have a core size of 3x3 TTs -> 0.3 x 0.3
182 for(const TLorentzVector& roi : rois) {
183 L1_match = offline_tau->p4().DeltaR(roi) <= 0.3;
184 if(L1_match) break;
185 }
186
187 tauPt = offline_tau->pt()/Gaudi::Units::GeV;
188 tauEta = offline_tau->eta();
189 tauPhi = offline_tau->phi();
190 averageMu = lbAverageInteractionsPerCrossing(ctx);
191
192 bool is_highPt = tauPt > info.getHLTTauThreshold() + 20.0;
193
194 // HLT matching: dR matching + HLT fires
195 HLT_match = matchObjects(offline_tau, online_tau_vec, 0.2) && hlt_fires;
196
197 // Total efficiency (without L1 matching)
199 Total_match = static_cast<bool>(HLT_match);
200 fill(monGroup, tauPt, tauEta, tauPhi, Total_match);
201
202 if(is_highPt) {
203 Total_match_highPt = static_cast<bool>(HLT_match);
204 fill(monGroup, tauEta, tauPhi, Total_match_highPt);
205 }
206 }
207
208 if(!L1_match || !l1_accept_flag) continue; // Skip this offline tau since not matched with L1 item
209
210 fill(monGroup, tauPt, tauEta, tauPhi, averageMu, HLT_match);
211
212 if(is_highPt) {
213 HLT_match_highPt = static_cast<bool>(HLT_match);
214 fill(monGroup, tauEta, tauPhi, HLT_match_highPt);
215 }
216 }
217
218 ATH_MSG_DEBUG("After fill HLT efficiencies: " << trigger);
219}
bool matchObjects(const T1 *tau, const std::vector< const T2 * > &tau_vec, float threshold) const
const TrigTauInfo & getTrigInfo(const std::string &trigger) const
std::vector< TLorentzVector > getRoIsVector(const EventContext &ctx, const std::string &trigger) const
double tauPt(const xAOD::TauJet &xTau)
return tau pt in GeV
double tauEta(const xAOD::TauJet &xTau)
return tau eta

◆ fillIDCluster()

void TrigTauMonitorSingleAlgorithm::fillIDCluster ( const std::string & trigger,
const std::vector< const xAOD::TauJet * > & tau_vec,
bool online ) const
private

Definition at line 362 of file TrigTauMonitorSingleAlgorithm.cxx.

363{
364 ATH_MSG_DEBUG("Fill ID input Cluster: " << trigger << " for online/offline " << online);
365
366 auto monGroup = getGroup(trigger+"_ID_"+(online ? "HLT" : "Offline")+"_InputCluster");
367
368 for(const auto *tau : tau_vec){
369 auto cluster_pt_jetseed_log = Monitored::Collection("cluster_pt_jetseed_log", tau_vec, [](const xAOD::TauJet* tau){ return std::log10(tau->ptJetSeed()); });
370
371 std::vector<const xAOD::CaloCluster*> clusters;
372 for(const xAOD::IParticle* particle : tau->clusters()) {
373 const xAOD::CaloCluster* cluster = static_cast<const xAOD::CaloCluster*>(particle);
374 clusters.push_back(cluster);
375 }
376 std::sort(clusters.begin(), clusters.end(), [](const xAOD::CaloCluster *lhs, const xAOD::CaloCluster *rhs){ return lhs->et() > rhs->et(); });
377
378 auto n_cluster = Monitored::Scalar<int>("n_cluster", 0);
379 n_cluster = clusters.size();
380
381 auto cluster_et_log = Monitored::Collection("cluster_et_log",clusters, [](const xAOD::CaloCluster *cluster){ return std::log10( cluster->et()); });
382 auto cluster_eta = Monitored::Collection("cluster_eta", clusters, [](const xAOD::CaloCluster *cluster){ return cluster->eta(); });
383 auto cluster_phi = Monitored::Collection("cluster_phi", clusters, [](const xAOD::CaloCluster *cluster){ return cluster->phi(); });
384 auto cluster_dEta = Monitored::Collection("cluster_dEta", clusters, [&tau](const xAOD::CaloCluster *cluster){ return cluster->eta() - tau->eta(); });
385 auto cluster_dPhi = Monitored::Collection("cluster_dPhi", clusters, [&tau](const xAOD::CaloCluster *cluster){ return cluster->p4().DeltaPhi(tau->p4()); });
386 auto cluster_SECOND_R_log10 = Monitored::Collection("cluster_SECOND_R_log10", clusters, [](const xAOD::CaloCluster *cluster){
387 double detail = -999;
388 const auto success_SECOND_R = cluster->retrieveMoment(xAOD::CaloCluster::MomentType::SECOND_R,detail);
389 if(success_SECOND_R) detail = std::log10(detail + 0.1);
390 return detail;
391 });
392
393 auto cluster_SECOND_LAMBDA_log10 = Monitored::Collection("cluster_SECOND_LAMBDA_log10", clusters, [](const xAOD::CaloCluster *cluster){
394 double detail = -999;
395 const auto success_SECOND_LAMBDA = cluster->retrieveMoment(xAOD::CaloCluster::MomentType::SECOND_LAMBDA, detail);
396 if(success_SECOND_LAMBDA) detail = std::log10(detail + 0.1);
397 return detail;
398 });
399
400 auto cluster_CENTER_LAMBDA_log10 = Monitored::Collection("cluster_CENTER_LAMBDA_log10", clusters, [](const xAOD::CaloCluster *cluster){
401 double detail = -999;
402 const auto success_CENTER_LAMBDA = cluster->retrieveMoment(xAOD::CaloCluster::MomentType::CENTER_LAMBDA, detail);
403 if(success_CENTER_LAMBDA) detail = std::log10(detail + 1e-6);
404 return detail;
405 });
406
407 fill(monGroup, n_cluster, cluster_pt_jetseed_log, cluster_et_log, cluster_eta, cluster_phi, cluster_dEta, cluster_dPhi, cluster_SECOND_R_log10, cluster_SECOND_LAMBDA_log10, cluster_CENTER_LAMBDA_log10);
408 }
409
410 ATH_MSG_DEBUG("After fill ID input Cluster: " << trigger);
411}
bool retrieveMoment(MomentType type, double &value) const
Retrieve individual moment.
virtual double eta() const
The pseudorapidity ( ) of the particle.
virtual double phi() const
The azimuthal angle ( ) of the particle.
@ SECOND_LAMBDA
Second Moment in .
@ SECOND_R
Second Moment in .
@ CENTER_LAMBDA
Shower depth at Cluster Centroid.
virtual FourMom_t p4() const
The full 4-momentum of the particle.
virtual FourMom_t p4() const
The full 4-momentum of the particle.
Definition TauJet_v3.cxx:96
std::vector< const IParticle * > clusters() const
double ptJetSeed() const
constexpr ParticleHypothesis particle[PARTICLEHYPOTHESES]
the array of masses
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.

◆ fillIDInputVars()

void TrigTauMonitorSingleAlgorithm::fillIDInputVars ( const std::string & trigger,
const std::vector< const xAOD::TauJet * > & tau_vec,
const std::string & nProng,
bool online ) const
private

Definition at line 222 of file TrigTauMonitorSingleAlgorithm.cxx.

223{
224 ATH_MSG_DEBUG("Fill ID input variables: " << trigger);
225
226 auto monGroup = getGroup(trigger+"_ID_"+(online ? "HLT" : "Offline")+"_InputScalar_"+nProng);
227
228 auto centFrac = Monitored::Collection("centFrac", tau_vec, [](const xAOD::TauJet* tau){
229 float detail = -999;
230 if(tau->detail(xAOD::TauJetParameters::centFrac, detail)) detail = std::min(detail, 1.0f);
231 return detail;
232 });
233 auto etOverPtLeadTrk = Monitored::Collection("etOverPtLeadTrk", tau_vec, [](const xAOD::TauJet* tau){
234 float detail = -999;
235 if(tau->detail(xAOD::TauJetParameters::etOverPtLeadTrk, detail)) detail = std::log10(std::max(detail, 0.1f));
236 return detail;
237 });
238 auto dRmax = Monitored::Collection("dRmax", tau_vec, [](const xAOD::TauJet* tau){
239 float detail = -999;
241 return detail;
242 });
243 auto absipSigLeadTrk = Monitored::Collection("absipSigLeadTrk", tau_vec, [](const xAOD::TauJet* tau){
244 float detail = (tau->nTracks()>0) ? std::abs(tau->track(0)->d0SigTJVA()) : 0;
245 detail = std::min(std::abs(detail), 30.0f);
246 return detail;
247 });
248 auto sumPtTrkFrac = Monitored::Collection("sumPtTrkFrac", tau_vec, [](const xAOD::TauJet* tau){
249 float detail = -999;
251 return detail;
252 });
253 auto emPOverTrkSysP = Monitored::Collection("emPOverTrkSysP", tau_vec, [](const xAOD::TauJet* tau){
254 float detail = -999;
255 if(tau->detail(xAOD::TauJetParameters::EMPOverTrkSysP, detail)) detail = std::log10(std::max(detail, 1e-3f));
256 return detail;
257 });
258 auto ptRatioEflowApprox = Monitored::Collection("ptRatioEflowApprox", tau_vec, [](const xAOD::TauJet* tau){
259 float detail = -999;
260 if(tau->detail(xAOD::TauJetParameters::ptRatioEflowApprox, detail)) detail = std::min(detail, 4.0f);
261 return detail;
262 });
263 auto mEflowApprox = Monitored::Collection("mEflowApprox", tau_vec, [](const xAOD::TauJet* tau){
264 float detail = -999;
265 if(tau->detail(xAOD::TauJetParameters::mEflowApprox, detail)) detail = std::log10(std::max(detail, 140.0f));
266 return detail;
267 });
268 auto ptDetectorAxis = Monitored::Collection("ptDetectorAxis", tau_vec, [](const xAOD::TauJet* tau){
269 float detail = -999;
270 if( tau->ptDetectorAxis() > 0) detail = std::log10(std::min(tau->ptDetectorAxis()/Gaudi::Units::GeV, 100.0));
271 return detail;
272 });
273 auto massTrkSys = Monitored::Collection("massTrkSys", tau_vec, [&nProng](const xAOD::TauJet* tau){
274 float detail = -999;
275 if( tau->detail(xAOD::TauJetParameters::massTrkSys, detail) && (nProng.find("MP") != std::string::npos || nProng.find("3P") != std::string::npos)) {
276 detail = std::log10(std::max(detail, 140.0f));
277 }
278 return detail;});
279 auto trFlightPathSig = Monitored::Collection("trFlightPathSig", tau_vec, [&nProng](const xAOD::TauJet* tau){
280 float detail = -999;
281 if(nProng.find("MP") != std::string::npos || nProng.find("3P") != std::string::npos) tau->detail(xAOD::TauJetParameters::trFlightPathSig, detail);
282 return detail;
283 });
284
285 fill(monGroup, centFrac, etOverPtLeadTrk, dRmax, absipSigLeadTrk, sumPtTrkFrac, emPOverTrkSysP, ptRatioEflowApprox, mEflowApprox, ptDetectorAxis, massTrkSys, trFlightPathSig);
286
287 ATH_MSG_DEBUG("After fill ID input variables: " << trigger);
288}
double ptDetectorAxis() const
const TauTrack * track(size_t i, TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged, int *container_index=0) const
Get the pointer to a given tauTrack associated with this tau /*container index needed by trackNonCons...
size_t nTracks(TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged) const
float d0SigTJVA() const
bool ptDetectorAxis(const xAOD::TauJet &tau, float &out)
bool absipSigLeadTrk(const xAOD::TauJet &tau, float &out)
@ centFrac
Get centrality fraction.
Definition TauDefs.h:200
@ dRmax
Get maximal dR of tracks associated to calo-seeded tau.
Definition TauDefs.h:226

◆ fillIDScores()

void TrigTauMonitorSingleAlgorithm::fillIDScores ( const EventContext & ctx,
const std::string & trigger,
const std::vector< const xAOD::TauJet * > & tau_vec,
const std::string & nProng,
bool online ) const
private

Definition at line 456 of file TrigTauMonitorSingleAlgorithm.cxx.

457{
458 ATH_MSG_DEBUG("Fill TauID Scores: " << trigger);
459
460 const TrigTauInfo& info = getTrigInfo(trigger);
461 const std::string tau_id = info.getHLTTauID();
462 bool store_all = tau_id == "idperf" || tau_id == "perf";
463
464 if(online) {
465 if(m_monitoredHLTIdDecorHandleKeys.find(info.getHLTTauType()) == m_monitoredHLTIdDecorHandleKeys.end()) return;
466 } else {
467 if(m_monitoredOfflineIdDecorHandleKeys.size() == 0) return;
468 }
469
470
471 const auto& decor_handle_keys = online ? m_monitoredHLTIdDecorHandleKeys.at(info.getHLTTauType()) : m_monitoredOfflineIdDecorHandleKeys;
472
473 // This has to be down here, because otherwise we crash on unmonitored chains
474 auto monGroup = getGroup(trigger+"_"+(online ? "HLT" : "Offline")+"_IDScores_"+nProng);
475
476 for(const auto& [key, p] : decor_handle_keys) {
477 // We will either store the TauID scores indicated in the chain name, or all of them in case of (id)perf chains
478 if(online && !store_all && tau_id != key) continue;
479
480 SG::ReadDecorHandle<xAOD::TauJetContainer, float> score_handle(p.first, ctx);
481 SG::ReadDecorHandle<xAOD::TauJetContainer, float> score_sig_trans_handle(p.second, ctx);
482
483 // Skip if both the Score and ScoreSigTrans aren't available
484 if(!score_handle.isAvailable() || !score_sig_trans_handle.isAvailable()) continue;
485
486 std::vector<float> score, score_sig_trans;
487 for(const xAOD::TauJet* tau : tau_vec) {
488 score.push_back(score_handle(*tau));
489 score_sig_trans.push_back(score_sig_trans_handle(*tau));
490 }
491
492 auto IDScore = Monitored::Collection(key + "_TauIDScore", score);
493 auto IDScoreSigTrans = Monitored::Collection(key + "_TauIDScoreSigTrans", score_sig_trans);
494 fill(monGroup, IDScore, IDScoreSigTrans);
495 }
496
497 ATH_MSG_DEBUG("After fill TauID Scores: " << trigger);
498}
std::map< std::string, std::map< std::string, std::pair< SG::ReadDecorHandleKey< xAOD::TauJetContainer >, SG::ReadDecorHandleKey< xAOD::TauJetContainer > > > > m_monitoredHLTIdDecorHandleKeys
std::map< std::string, std::pair< SG::ReadDecorHandleKey< xAOD::TauJetContainer >, SG::ReadDecorHandleKey< xAOD::TauJetContainer > > > m_monitoredOfflineIdDecorHandleKeys

◆ fillIDTrack()

void TrigTauMonitorSingleAlgorithm::fillIDTrack ( const std::string & trigger,
const std::vector< const xAOD::TauJet * > & tau_vec,
bool online ) const
private

Definition at line 291 of file TrigTauMonitorSingleAlgorithm.cxx.

292{
293 ATH_MSG_DEBUG("Fill ID input Track: " << trigger);
294
295 auto monGroup = getGroup(trigger+"_ID_"+(online ? "HLT" : "Offline")+"_InputTrack");
296
297 auto track_pt_jetseed_log = Monitored::Collection("track_pt_jetseed_log", tau_vec, [](const xAOD::TauJet* tau){ return std::log10(tau->ptJetSeed()); });
298 fill(monGroup, track_pt_jetseed_log);
299
300 for(const auto *tau : tau_vec) {
301 // Don't call ->allTracks() unless the element links are valid
302 static const SG::ConstAccessor< std::vector<ElementLink<xAOD::TauTrackContainer>> > tauTrackAcc("tauTrackLinks");
303 bool linksValid = true;
304 for(const ElementLink<xAOD::TauTrackContainer>& trackEL : tauTrackAcc(*tau)) {
305 if(!trackEL.isValid()) {
306 linksValid = false;
307 break;
308 }
309 }
310 if(!linksValid) {
311 ATH_MSG_WARNING("Invalid track element links from TauJet in " << trigger);
312 continue;
313 }
314
315 auto tracks = tau->allTracks();
316 std::sort(tracks.begin(), tracks.end(), [](const xAOD::TauTrack* lhs, const xAOD::TauTrack* rhs){ return lhs->pt() > rhs->pt(); });
317
318 auto n_track = Monitored::Scalar<int>("n_track", tracks.size());
319
320 auto track_pt_log = Monitored::Collection("track_pt_log", tracks, [](const xAOD::TauTrack *track){ return std::log10(track->pt()); });
321 auto track_eta = Monitored::Collection("track_eta", tracks, [](const xAOD::TauTrack *track){ return track->eta(); });
322 auto track_phi = Monitored::Collection("track_phi", tracks, [](const xAOD::TauTrack *track){ return track->phi(); });
323
324 auto track_dEta = Monitored::Collection("track_dEta", tracks, [&tau](const xAOD::TauTrack *track){ return track->eta() - tau->eta(); });
325 auto track_dPhi = Monitored::Collection("track_dPhi", tracks, [&tau](const xAOD::TauTrack *track){ return track->p4().DeltaPhi(tau->p4()); });
326
327 auto track_z0sinthetaTJVA_abs_log = Monitored::Collection("track_z0sinthetaTJVA_abs_log", tracks, [](const xAOD::TauTrack *track){return track->z0sinthetaTJVA(); });
328 auto track_d0_abs_log = Monitored::Collection("track_d0_abs_log", tracks, [](const xAOD::TauTrack *track){ return std::log10(std::abs(track->track()->d0()) + 1e-6); });
329
330 auto track_nIBLHitsAndExp = Monitored::Collection("track_nIBLHitsAndExp", tracks, [](const xAOD::TauTrack *track){
331 uint8_t inner_pixel_hits, inner_pixel_exp;
332 const auto success1_innerPixel_hits = track->track()->summaryValue(inner_pixel_hits, xAOD::numberOfInnermostPixelLayerHits);
333 const auto success2_innerPixel_exp = track->track()->summaryValue(inner_pixel_exp, xAOD::expectInnermostPixelLayerHit);
334 float detail = -999;
335 if(success1_innerPixel_hits && success2_innerPixel_exp) { detail = inner_pixel_exp ? inner_pixel_hits : 1.; };
336 return detail;
337 });
338 auto track_nPixelHitsPlusDeadSensors = Monitored::Collection("track_nPixelHitsPlusDeadSensors", tracks, [](const xAOD::TauTrack *track){
339 uint8_t pixel_hits, pixel_dead;
340 const auto success1_pixel_hits = track->track()->summaryValue(pixel_hits, xAOD::numberOfPixelHits);
341 const auto success2_pixel_dead = track->track()->summaryValue(pixel_dead, xAOD::numberOfPixelDeadSensors);
342 float detail = -999;
343 if(success1_pixel_hits && success2_pixel_dead) { detail = pixel_hits + pixel_dead; };
344 return detail;
345 });
346 auto track_nSCTHitsPlusDeadSensors = Monitored::Collection("track_nSCTHitsPlusDeadSensors", tracks, [](const xAOD::TauTrack *track){
347 uint8_t sct_hits, sct_dead;
348 const auto success1_sct_hits = track->track()->summaryValue(sct_hits, xAOD::numberOfSCTHits);
349 const auto success2_sct_dead = track->track()->summaryValue(sct_dead, xAOD::numberOfSCTDeadSensors);
350 float detail = -999;
351 if(success1_sct_hits && success2_sct_dead) { detail = sct_hits + sct_dead; };
352 return detail;
353 });
354
355 fill(monGroup, n_track, track_pt_log, track_eta, track_phi, track_dEta, track_dPhi, track_z0sinthetaTJVA_abs_log, track_d0_abs_log, track_nIBLHitsAndExp, track_nPixelHitsPlusDeadSensors, track_nSCTHitsPlusDeadSensors);
356 }
357
358 ATH_MSG_DEBUG("After fill ID input Track: " << trigger);
359}
std::vector< const TauTrack * > allTracks() const
Get the v<const pointer> to all tracks associated with this tau, regardless of classification.
TauTrack_v1 TauTrack
Definition of the current version.
Definition TauTrack.h:16
@ expectInnermostPixelLayerHit
Do we expect a 0th-layer barrel hit for this track?
@ numberOfSCTDeadSensors
number of dead SCT sensors crossed [unit8_t].
@ numberOfSCTHits
number of hits in SCT [unit8_t].
@ numberOfInnermostPixelLayerHits
these are the hits in the 0th pixel barrel layer
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
@ numberOfPixelDeadSensors
number of dead pixel sensors crossed [unit8_t].

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

◆ getL1cTAUs()

std::vector< std::pair< const xAOD::eFexTauRoI *, const xAOD::jFexTauRoI * > > TrigTauMonitorBaseAlgorithm::getL1cTAUs ( const EventContext & ctx,
const std::string & l1_item ) const
protectedinherited

Definition at line 193 of file TrigTauMonitorBaseAlgorithm.cxx.

194{
195 std::vector<std::pair<const xAOD::eFexTauRoI*, const xAOD::jFexTauRoI*>> roi_vec;
196
197 SG::ReadHandle<xAOD::eFexTauRoIContainer> rois(m_phase1l1cTauRoIKey, ctx);
198 if(!rois.isValid()) {
199 ATH_MSG_WARNING("Failed to retrieve the L1_cTauRoi container");
200 return roi_vec;
201 }
202 SG::ReadDecorHandle<xAOD::eFexTauRoIContainer, ElementLink<xAOD::jFexTauRoIContainer>> jTau_roi_link{m_phase1l1cTauRoIDecorKey, ctx};
203 if(!jTau_roi_link.isValid()) {
204 ATH_MSG_WARNING("Failed to create jTauLink accessor for the L1_cTauRoi container");
205 return roi_vec;
206 }
207
209 for(size_t i = 0; i < rois->size(); i++) {
210 const xAOD::eFexTauRoI* roi = (*rois)[i];
211 const xAOD::jFexTauRoI* jTau_roi = jTau_roi_link(i).isValid() ? *jTau_roi_link(i) : nullptr;
212
213 // Select by RoI ET value only
214 if(roi->et() > m_L1_Phase1_thresholds.value().at(l1_item)) roi_vec.push_back(std::make_pair(roi, jTau_roi));
215 }
216 } else {
217 SG::ReadDecorHandle<xAOD::eFexTauRoIContainer, uint64_t> thresholdPatterns(m_phase1l1cTauRoIThresholdPatternsKey, ctx);
218 if(!thresholdPatterns.isValid()) {
219 ATH_MSG_WARNING("Failed to create thresholdPatterns property accessor for the L1_cTauRoi container");
220 return roi_vec;
221 }
222
223 for(size_t i = 0; i < rois->size(); i++) {
224 const xAOD::eFexTauRoI* roi = (*rois)[i];
225 const xAOD::jFexTauRoI* jTau_roi = jTau_roi_link(i).isValid() ? *jTau_roi_link(i) : nullptr;
226
227 // Check that the RoI passed the threshold selection
228 if(thresholdPatterns(*roi) & m_L1_Phase1_threshold_patterns.value().at(l1_item)) roi_vec.push_back(std::make_pair(roi, jTau_roi));
229 }
230 }
231
232
233 return roi_vec;
234}
Gaudi::Property< std::map< std::string, float > > m_L1_Phase1_thresholds
SG::ReadDecorHandleKey< xAOD::eFexTauRoIContainer > m_phase1l1cTauRoIDecorKey
SG::ReadDecorHandleKey< xAOD::eFexTauRoIContainer > m_phase1l1cTauRoIThresholdPatternsKey
Gaudi::Property< bool > m_L1_select_by_et_only
SG::ReadHandleKey< xAOD::eFexTauRoIContainer > m_phase1l1cTauRoIKey
Gaudi::Property< std::map< std::string, uint64_t > > m_L1_Phase1_threshold_patterns
float et() const
TOB ET (decoded from TOB, stored for convenience)
jFexTauRoI_v1 jFexTauRoI
Define the latest version of the jFexSRJetRoI class.
Definition jFexTauRoI.h:13
eFexTauRoI_v1 eFexTauRoI
Define the latest version of the eFexTauRoI class.
Definition eFexTauRoI.h:16

◆ getL1eTAUs()

std::vector< const xAOD::eFexTauRoI * > TrigTauMonitorBaseAlgorithm::getL1eTAUs ( const EventContext & ctx,
const std::string & l1_item ) const
protectedinherited

Definition at line 129 of file TrigTauMonitorBaseAlgorithm.cxx.

130{
131 std::vector<const xAOD::eFexTauRoI*> roi_vec;
132
133 SG::ReadHandle<xAOD::eFexTauRoIContainer> rois(m_phase1l1eTauRoIKey, ctx);
134 if(!rois.isValid()) {
135 ATH_MSG_WARNING("Failed to retrieve the L1_eTauRoi container");
136 return roi_vec;
137 }
138
140 for(const xAOD::eFexTauRoI* roi : *rois) {
141 // Select by RoI ET value only
142 if(roi->et() > m_L1_Phase1_thresholds.value().at(l1_item)) roi_vec.push_back(roi);
143 }
144 } else {
145 SG::ReadDecorHandle<xAOD::eFexTauRoIContainer, uint64_t> thresholdPatterns(m_phase1l1eTauRoIThresholdPatternsKey, ctx);
146 if(!thresholdPatterns.isValid()) {
147 ATH_MSG_WARNING("Failed to create thresholdPatterns property accessor for the L1_eTauRoi container");
148 return roi_vec;
149 }
150
151 for(const xAOD::eFexTauRoI* roi : *rois) {
152 // Check that the RoI passed the threshold selection
153 if(thresholdPatterns(*roi) & m_L1_Phase1_threshold_patterns.value().at(l1_item)) roi_vec.push_back(roi);
154 }
155 }
156
157 return roi_vec;
158}
SG::ReadHandleKey< xAOD::eFexTauRoIContainer > m_phase1l1eTauRoIKey
SG::ReadDecorHandleKey< xAOD::eFexTauRoIContainer > m_phase1l1eTauRoIThresholdPatternsKey

◆ getL1jTAUs()

std::vector< const xAOD::jFexTauRoI * > TrigTauMonitorBaseAlgorithm::getL1jTAUs ( const EventContext & ctx,
const std::string & l1_item ) const
protectedinherited

Definition at line 161 of file TrigTauMonitorBaseAlgorithm.cxx.

162{
163 std::vector<const xAOD::jFexTauRoI*> roi_vec;
164
165 SG::ReadHandle<xAOD::jFexTauRoIContainer> rois(m_phase1l1jTauRoIKey, ctx);
166 if(!rois.isValid()) {
167 ATH_MSG_WARNING("Failed to retrieve the L1_jTauRoi container");
168 return roi_vec;
169 }
170
172 for(const xAOD::jFexTauRoI* roi : *rois) {
173 // Select by RoI ET value only
174 if(roi->et() > m_L1_Phase1_thresholds.value().at(l1_item)) roi_vec.push_back(roi);
175 }
176 } else {
177 SG::ReadDecorHandle<xAOD::jFexTauRoIContainer, uint64_t> thresholdPatterns(m_phase1l1jTauRoIThresholdPatternsKey, ctx);
178 if(!thresholdPatterns.isValid()) {
179 ATH_MSG_WARNING("Failed to create thresholdPatterns property accessor for the L1_jTauRoi container");
180 return roi_vec;
181 }
182
183 for(const xAOD::jFexTauRoI* roi : *rois) {
184 // Check that the RoI passed the threshold selection
185 if(thresholdPatterns(*roi) & m_L1_Phase1_threshold_patterns.value().at(l1_item)) roi_vec.push_back(roi);
186 }
187 }
188
189 return roi_vec;
190}
SG::ReadHandleKey< xAOD::jFexTauRoIContainer > m_phase1l1jTauRoIKey
SG::ReadDecorHandleKey< xAOD::jFexTauRoIContainer > m_phase1l1jTauRoIThresholdPatternsKey

◆ getOfflineTaus()

std::pair< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > TrigTauMonitorBaseAlgorithm::getOfflineTaus ( const EventContext & ctx,
const float threshold = 20.0,
const TauID tau_id = TauID::None ) const
protectedinherited

Definition at line 123 of file TrigTauMonitorBaseAlgorithm.cxx.

124{
125 return classifyOfflineTaus(getOfflineTausAll(ctx, threshold), threshold, tau_id);
126}
std::vector< const xAOD::TauJet * > getOfflineTausAll(const EventContext &ctx, const float threshold=20.0) const
std::pair< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > classifyOfflineTaus(const std::vector< const xAOD::TauJet * > &taus, const float threshold=0.0, const TauID tau_id=TauID::None) const

◆ getOfflineTausAll()

std::vector< const xAOD::TauJet * > TrigTauMonitorBaseAlgorithm::getOfflineTausAll ( const EventContext & ctx,
const float threshold = 20.0 ) const
protectedinherited

Definition at line 87 of file TrigTauMonitorBaseAlgorithm.cxx.

88{
89 ATH_MSG_DEBUG("Retrieving offline Taus");
90
91 std::vector<const xAOD::TauJet*> tau_vec;
92
93 SG::ReadHandle<xAOD::TauJetContainer> taus(m_offlineTauJetKey, ctx);
94 if(!taus.isValid()) {
95 ATH_MSG_WARNING("Failed to retrieve offline Taus");
96 return tau_vec;
97 }
98
99 for(const xAOD::TauJet* const tau : *taus) {
100 // Consider only offline taus with a certain minimum pT
101 if(tau->pt() < threshold*Gaudi::Units::GeV) continue;
102
103 // Consider only offline taus outside of the crack region
104 if(std::abs(tau->eta()) > 1.37 && std::abs(tau->eta()) < 1.52) continue;
105
106 // Consider only offline taus which pass RNN medium WP
107 if(!tau->isTau(xAOD::TauJetParameters::JetRNNSigMedium)) continue;
108
109 // Consider only offline taus which pass thinning
110 static const SG::ConstAccessor<char> passThinningAcc("passThinning");
111 if(!passThinningAcc.withDefault(*tau, true)) continue;
112
113 int nTracks = -1;
114 tau->detail(xAOD::TauJetParameters::nChargedTracks, nTracks);
115 ATH_MSG_DEBUG("NTracks Offline: " << nTracks);
116 if(nTracks == 1 || nTracks == 3) tau_vec.push_back(tau);
117 }
118
119 return tau_vec;
120}
SG::ReadHandleKey< xAOD::TauJetContainer > m_offlineTauJetKey

◆ getOnlineContainerKey()

const SG::ReadHandleKey< xAOD::TauJetContainer > & TrigTauMonitorBaseAlgorithm::getOnlineContainerKey ( const std::string & sequence) const
protectedinherited

Definition at line 237 of file TrigTauMonitorBaseAlgorithm.cxx.

238{
239 if(sequence == "tracktwoMVA" || sequence == "tracktwoMVABDT") return m_hltTauJetKey;
240 else if(sequence == "tracktwoLLP") return m_hltTauJetLLPKey;
241 else if(sequence == "trackLRT") return m_hltTauJetLRTKey;
242 else if(sequence == "ptonly") return m_hltTauJetCaloMVAOnlyKey;
243 else {
244 ATH_MSG_ERROR("Unknown HLT TauJet container for sequence \"" << sequence << "\". Returning the default \"" << m_hltTauJetKey.key() << "\"");
245 return m_hltTauJetKey;
246 }
247}
#define ATH_MSG_ERROR(x)
SG::ReadHandleKey< xAOD::TauJetContainer > m_hltTauJetKey
SG::ReadHandleKey< xAOD::TauJetContainer > m_hltTauJetLLPKey
SG::ReadHandleKey< xAOD::TauJetContainer > m_hltTauJetCaloMVAOnlyKey
SG::ReadHandleKey< xAOD::TauJetContainer > m_hltTauJetLRTKey

◆ getOnlineTaus()

std::tuple< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > TrigTauMonitorBaseAlgorithm::getOnlineTaus ( const std::string & trigger) const
protectedinherited

Definition at line 81 of file TrigTauMonitorBaseAlgorithm.cxx.

82{
83 return classifyOnlineTaus(getOnlineTausAll(trigger, true), 0.0);
84}
std::tuple< std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * >, std::vector< const xAOD::TauJet * > > classifyOnlineTaus(const std::vector< const xAOD::TauJet * > &taus, const float threshold=0.0) const
std::vector< const xAOD::TauJet * > getOnlineTausAll(const std::string &trigger, bool include_0P=true) const

◆ getOnlineTausAll()

std::vector< const xAOD::TauJet * > TrigTauMonitorBaseAlgorithm::getOnlineTausAll ( const std::string & trigger,
bool include_0P = true ) const
protectedinherited

Definition at line 57 of file TrigTauMonitorBaseAlgorithm.cxx.

58{
59 std::vector<const xAOD::TauJet*> tau_vec;
60
61
62 const TrigTauInfo& info = getTrigInfo(trigger);
63 const std::string tau_container_name = getOnlineContainerKey(info.getHLTTauType()).key();
64 ATH_MSG_DEBUG("Tau container name is: " << tau_container_name);
65 auto vec = m_trigDecTool->features<xAOD::TauJetContainer>(trigger, TrigDefs::Physics, tau_container_name);
66 for(auto& featLinkInfo : vec) {
67 const xAOD::TauJet* feat = *(featLinkInfo.link);
68 if(!feat) continue;
69
70 int nTracks = -1;
72 ATH_MSG_DEBUG("NTracks Online: " << nTracks);
73
74 if(include_0P && nTracks == 0) tau_vec.push_back(feat);
75 else tau_vec.push_back(feat);
76 }
77
78 return tau_vec;
79}
std::vector< size_t > vec
const SG::ReadHandleKey< xAOD::TauJetContainer > & getOnlineContainerKey(const std::string &sequence) const
TauJetContainer_v3 TauJetContainer
Definition of the current "taujet container version".

◆ getRoIsVector()

std::vector< TLorentzVector > TrigTauMonitorSingleAlgorithm::getRoIsVector ( const EventContext & ctx,
const std::string & trigger ) const
private

Definition at line 501 of file TrigTauMonitorSingleAlgorithm.cxx.

502{
503 std::vector<TLorentzVector> ret;
504
505 const TrigTauInfo& info = getTrigInfo(trigger);
506
507 TLorentzVector v;
508 if(info.getL1TauType() == "eTAU") {
509 for(const xAOD::eFexTauRoI* roi : getL1eTAUs(ctx, info.getL1TauItem())) {
510 v.SetPtEtaPhiM(roi->et(), roi->eta(), roi->phi(), 0);
511 ret.push_back(v);
512 }
513 } else if(info.getL1TauType() == "jTAU") {
514 for(const xAOD::jFexTauRoI* roi : getL1jTAUs(ctx, info.getL1TauItem())) {
515 v.SetPtEtaPhiM(roi->et(), roi->eta(), roi->phi(), 0);
516 ret.push_back(v);
517 }
518 } else if(info.getL1TauType() == "cTAU") {
519 for(const auto& [eTau_roi, jTau_roi] : getL1cTAUs(ctx, info.getL1TauItem())) {
520 v.SetPtEtaPhiM(eTau_roi->et(), eTau_roi->eta(), eTau_roi->phi(), 0);
521 ret.push_back(v);
522 }
523 }
524
525 return ret;
526}
std::vector< const xAOD::eFexTauRoI * > getL1eTAUs(const EventContext &ctx, const std::string &l1_item) const
std::vector< const xAOD::jFexTauRoI * > getL1jTAUs(const EventContext &ctx, const std::string &l1_item) const
std::vector< std::pair< const xAOD::eFexTauRoI *, const xAOD::jFexTauRoI * > > getL1cTAUs(const EventContext &ctx, const std::string &l1_item) const

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

◆ getTrigInfo()

const TrigTauInfo & TrigTauMonitorBaseAlgorithm::getTrigInfo ( const std::string & trigger) const
inlineprotectedinherited

Definition at line 51 of file TrigTauMonitorBaseAlgorithm.h.

51{ return m_trigInfo.at(trigger); }
std::map< std::string, TrigTauInfo > m_trigInfo

◆ getTrigInfoMap()

std::map< std::string, TrigTauInfo > & TrigTauMonitorBaseAlgorithm::getTrigInfoMap ( )
inlineprotectedinherited

Definition at line 50 of file TrigTauMonitorBaseAlgorithm.h.

50{ return m_trigInfo; }

◆ initialize()

StatusCode TrigTauMonitorSingleAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from TrigTauMonitorBaseAlgorithm.

Definition at line 16 of file TrigTauMonitorSingleAlgorithm.cxx.

17{
19
20 // Create the TauID score decor handle keys for the Monitoring...
21 for(const auto& [seq_name, m] : m_monitoredHLTIdScores) {
22 const std::string online_container_name = getOnlineContainerKey(seq_name).key();
23 for(const auto& [key, p] : m) {
24 if(p.first.empty() || p.second.empty()) {
25 ATH_MSG_WARNING("Invalid HLT TauID score variable names; skipping this entry for the monitoring!");
26 continue;
27 }
28
29 m_monitoredHLTIdDecorHandleKeys[seq_name].emplace(
30 key,
31 std::make_pair(
32 SG::ReadDecorHandleKey<xAOD::TauJetContainer>(online_container_name + "." + p.first),
33 SG::ReadDecorHandleKey<xAOD::TauJetContainer>(online_container_name + "." + p.second)
34 )
35 );
36 ATH_CHECK(m_monitoredHLTIdDecorHandleKeys.at(seq_name).at(key).first.initialize());
37 ATH_CHECK(m_monitoredHLTIdDecorHandleKeys.at(seq_name).at(key).second.initialize());
38 }
39 }
40
41 for(const auto& [key, p] : m_monitoredOfflineIdScores) {
42 if(p.first.empty() || p.second.empty()) {
43 ATH_MSG_WARNING("Invalid Offline TauID score variable names; skipping this entry for the monitoring!");
44 continue;
45 }
46
48 key,
49 std::make_pair(
50 SG::ReadDecorHandleKey<xAOD::TauJetContainer>(m_offlineTauJetKey.key() + "." + p.first),
51 SG::ReadDecorHandleKey<xAOD::TauJetContainer>(m_offlineTauJetKey.key() + "." + p.second)
52 )
53 );
54 ATH_CHECK(m_monitoredOfflineIdDecorHandleKeys.at(key).first.initialize());
55 ATH_CHECK(m_monitoredOfflineIdDecorHandleKeys.at(key).second.initialize());
56 }
57
58 return StatusCode::SUCCESS;
59}
virtual StatusCode initialize() override
initialize
Gaudi::Property< std::map< std::string, std::map< std::string, std::pair< std::string, std::string > > > > m_monitoredHLTIdScores
Gaudi::Property< std::map< std::string, std::pair< std::string, std::string > > > m_monitoredOfflineIdScores

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

◆ matchObjects() [1/3]

template<typename T1 = xAOD::IParticle, typename T2 = xAOD::IParticle>
bool TrigTauMonitorBaseAlgorithm::matchObjects ( const T1 * tau,
const std::vector< const T2 * > & tau_vec,
float threshold ) const
inlineprotectedinherited

Definition at line 78 of file TrigTauMonitorBaseAlgorithm.h.

79 {
80 for(auto tau_2 : tau_vec) {
81 if(tau->p4().DeltaR(tau_2->p4()) < threshold) return true;
82 }
83 return false;
84 }

◆ matchObjects() [2/3]

template<typename T1 = xAOD::IParticle, typename T2 = xAOD::eFexTauRoI>
bool TrigTauMonitorBaseAlgorithm::matchObjects ( const T1 * tau_1,
const T2 * tau_2,
float threshold ) const
inlineprotectedinherited

Definition at line 95 of file TrigTauMonitorBaseAlgorithm.h.

96 {
97 return dR(tau_1->eta(), tau_1->phi(), tau_2->eta(), tau_2->phi()) < threshold;
98 }
double dR(const double eta1, const double phi1, const double eta2, const double phi2) const

◆ matchObjects() [3/3]

bool TrigTauMonitorBaseAlgorithm::matchObjects ( const TLorentzVector & tau,
const std::vector< TLorentzVector > & tau_vec,
float threshold ) const
inlineprotectedinherited

Definition at line 86 of file TrigTauMonitorBaseAlgorithm.h.

87 {
88 for(auto& tau_2 : tau_vec) {
89 if(tau.DeltaR(tau_2) < threshold) return true;
90 }
91 return false;
92 }

◆ matchTruthObjects()

template<typename T1 = xAOD::IParticle, typename T2 = xAOD::IParticle>
bool TrigTauMonitorBaseAlgorithm::matchTruthObjects ( const T1 * true_tau,
const std::vector< const T2 * > & tau_vec,
float threshold ) const
inlineprotectedinherited

Definition at line 101 of file TrigTauMonitorBaseAlgorithm.h.

102 {
103 static const SG::ConstAccessor<double> acc_ptvis("pt_vis");
104 static const SG::ConstAccessor<double> acc_etavis("eta_vis");
105 static const SG::ConstAccessor<double> acc_phivis("phi_vis");
106 static const SG::ConstAccessor<double> acc_mvis("mvis");
107 TLorentzVector true_tau_p4;
108 true_tau_p4.SetPtEtaPhiM(acc_ptvis(*true_tau), acc_etavis(*true_tau), acc_phivis(*true_tau), acc_mvis(*true_tau));
109
110 for(auto tau : tau_vec) {
111 if(true_tau_p4.DeltaR(tau->p4()) < threshold) return true;
112 }
113 return false;
114 }

◆ 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

◆ processEvent()

StatusCode TrigTauMonitorSingleAlgorithm::processEvent ( const EventContext & ctx) const
overrideprivatevirtual

Implements TrigTauMonitorBaseAlgorithm.

Definition at line 62 of file TrigTauMonitorSingleAlgorithm.cxx.

63{
64 constexpr float threshold_offset = 10.0;
65
66 // Offline taus
67 auto offline_taus_all = getOfflineTausAll(ctx, 0.0);
68 if(m_requireOfflineTaus && offline_taus_all.empty()) return StatusCode::SUCCESS;
69
70 for(const std::string& trigger : m_triggers) {
71 const TrigTauInfo& info = getTrigInfo(trigger);
72
73 if(!info.isHLTSingleTau()) {
74 ATH_MSG_WARNING("Chain \"" << trigger << "\" is not a single tau trigger. Skipping...");
75 continue;
76 }
77
78 const auto passBits = m_trigDecTool->isPassedBits(trigger);
79 const bool l1_accept_flag = passBits & TrigDefs::L1_isPassedAfterVeto;
80 const bool hlt_not_prescaled_flag = (passBits & TrigDefs::EF_prescaled) == 0;
81
82 // Offline tau requirement check
83 const std::vector<const xAOD::TauJet*> offline_taus_with_id = classifyTausAll(offline_taus_all, 0, static_cast<TauID>(m_offline_tau_id.value()));
84 if(m_requireOfflineTaus && offline_taus_with_id.empty()) continue;
85
86 // Filter offline taus
87 auto offline_taus = classifyOfflineTaus(offline_taus_with_id, info.getHLTTauThreshold() - threshold_offset);
88 std::vector<const xAOD::TauJet*> offline_taus_1p = offline_taus.first;
89 std::vector<const xAOD::TauJet*> offline_taus_3p = offline_taus.second;
90
91 // Online taus
92 std::vector<const xAOD::TauJet*> hlt_taus_all = getOnlineTausAll(trigger, true);
93 auto hlt_taus = classifyOnlineTaus(hlt_taus_all);
94 std::vector<const xAOD::TauJet*> hlt_taus_0p = std::get<0>(hlt_taus);
95 std::vector<const xAOD::TauJet*> hlt_taus_1p = std::get<1>(hlt_taus);
96 std::vector<const xAOD::TauJet*> hlt_taus_mp = std::get<2>(hlt_taus);
97
99 // Offline variables:
100 if(m_doOfflineTausDistributions && !offline_taus_1p.empty()) {
101 fillBasicVars(ctx, trigger, offline_taus_1p, "1P", false);
102 fillIDScores(ctx, trigger, offline_taus_1p, "1P", false);
103 fillIDInputVars(trigger, offline_taus_1p, "1P", false);
104 fillIDTrack(trigger, offline_taus_1p, false);
105 fillIDCluster(trigger, offline_taus_1p, false);
106 }
107 if(m_doOfflineTausDistributions && !offline_taus_3p.empty()) {
108 fillBasicVars(ctx, trigger, offline_taus_3p, "3P", false);
109 fillIDScores(ctx, trigger, offline_taus_3p, "3P", false);
110 fillIDInputVars(trigger, offline_taus_3p, "3P", false);
111 fillIDTrack(trigger, offline_taus_3p, false);
112 fillIDCluster(trigger, offline_taus_3p, false);
113 }
114
115 // Fill information for online 0 prong taus
116 if(!hlt_taus_0p.empty()) {
117 fillBasicVars(ctx, trigger, hlt_taus_0p, "0P", true);
118 fillIDScores(ctx, trigger, hlt_taus_0p, "0P", true);
119 fillIDInputVars(trigger, hlt_taus_0p, "0P", true);
120 fillIDTrack(trigger, hlt_taus_0p, true);
121 fillIDCluster(trigger, hlt_taus_0p, true);
122 }
123
124 // Fill information for online 1 prong taus
125 if(!hlt_taus_1p.empty()) {
126 fillBasicVars(ctx, trigger, hlt_taus_1p, "1P", true);
127 fillIDScores(ctx, trigger, hlt_taus_1p, "1P", true);
128 fillIDInputVars(trigger, hlt_taus_1p, "1P", true);
129 fillIDTrack(trigger, hlt_taus_1p, true);
130 fillIDCluster(trigger, hlt_taus_1p, true);
131 }
132
133 // Fill information for online multiprong prong taus
134 if(!hlt_taus_mp.empty()) {
135 fillBasicVars(ctx, trigger, hlt_taus_mp, "MP", true);
136 fillIDScores(ctx, trigger, hlt_taus_mp, "MP", true);
137 fillIDInputVars(trigger, hlt_taus_mp, "MP", true);
138 fillIDTrack(trigger, hlt_taus_mp, true);
139 fillIDCluster(trigger, hlt_taus_mp, true);
140 }
141 }
142
143 if(m_do_efficiency_plots && hlt_not_prescaled_flag) {
144 fillHLTEfficiencies(ctx, trigger, l1_accept_flag, offline_taus_1p, hlt_taus_all, "1P");
145 fillHLTEfficiencies(ctx, trigger, l1_accept_flag, offline_taus_3p, hlt_taus_all, "3P");
146 }
147 }
148
149 return StatusCode::SUCCESS;
150}
Gaudi::Property< bool > m_do_efficiency_plots
Gaudi::Property< std::vector< std::string > > m_triggers
Gaudi::Property< bool > m_do_variable_plots
void fillIDScores(const EventContext &ctx, const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, const std::string &nProng, bool online) const
void fillIDInputVars(const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, const std::string &nProng, bool online) const
void fillHLTEfficiencies(const EventContext &ctx, const std::string &trigger, const bool l1_accept_flag, const std::vector< const xAOD::TauJet * > &offline_tau_vec, const std::vector< const xAOD::TauJet * > &online_tau_vec, const std::string &nProng) const
Gaudi::Property< unsigned int > m_offline_tau_id
void fillIDCluster(const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, bool online) const
Gaudi::Property< bool > m_doOfflineTausDistributions
void fillBasicVars(const EventContext &ctx, const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, const std::string &nProng, bool online) const
void fillIDTrack(const std::string &trigger, const std::vector< const xAOD::TauJet * > &tau_vec, bool online) const

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

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

Gaudi::Property<bool> TrigTauMonitorBaseAlgorithm::m_do_efficiency_plots {this, "DoEfficiencyPlots", true}
protectedinherited

Definition at line 46 of file TrigTauMonitorBaseAlgorithm.h.

46{this, "DoEfficiencyPlots", true};

◆ m_do_variable_plots

Gaudi::Property<bool> TrigTauMonitorBaseAlgorithm::m_do_variable_plots {this, "DoVariablePlots", true}
protectedinherited

Definition at line 47 of file TrigTauMonitorBaseAlgorithm.h.

47{this, "DoVariablePlots", true};

◆ m_doOfflineTausDistributions

Gaudi::Property<bool> TrigTauMonitorSingleAlgorithm::m_doOfflineTausDistributions {this, "DoOfflineTausDistributions", true}
private

Definition at line 28 of file TrigTauMonitorSingleAlgorithm.h.

28{this, "DoOfflineTausDistributions", true};

◆ m_doTotalEfficiency

Gaudi::Property<bool> TrigTauMonitorSingleAlgorithm::m_doTotalEfficiency {this, "DoTotalEfficiency", false, "Do total efficiency histograms"}
private

Definition at line 21 of file TrigTauMonitorSingleAlgorithm.h.

21{this, "DoTotalEfficiency", false, "Do total efficiency histograms"};

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

SG::ReadDecorHandleKey<xAOD::EventInfo> TrigTauMonitorBaseAlgorithm::m_eventInfoDecorKey {this, "LArStatusFlag", "EventInfo.larFlags", "Key for EventInfo object"}
privateinherited

Definition at line 135 of file TrigTauMonitorBaseAlgorithm.h.

135{this, "LArStatusFlag", "EventInfo.larFlags", "Key for EventInfo object"}; // To get data-dependencies right

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

SG::ReadHandleKey<xAOD::TauJetContainer> TrigTauMonitorBaseAlgorithm::m_hltTauJetCaloMVAOnlyKey {this, "HLTTauJetCaloMVAOnlyKey", "HLT_TrigTauRecMerged_CaloMVAOnly", "HLT ptonly taujet container key"}
protectedinherited

Definition at line 126 of file TrigTauMonitorBaseAlgorithm.h.

126{this, "HLTTauJetCaloMVAOnlyKey", "HLT_TrigTauRecMerged_CaloMVAOnly", "HLT ptonly taujet container key"};

◆ m_hltTauJetKey

SG::ReadHandleKey<xAOD::TauJetContainer> TrigTauMonitorBaseAlgorithm::m_hltTauJetKey {this, "HLTTauJetKey", "HLT_TrigTauRecMerged_MVA", "HLT tracktwoMVA taujet container key"}
protectedinherited

Definition at line 123 of file TrigTauMonitorBaseAlgorithm.h.

123{this, "HLTTauJetKey", "HLT_TrigTauRecMerged_MVA", "HLT tracktwoMVA taujet container key"};

◆ m_hltTauJetLLPKey

SG::ReadHandleKey<xAOD::TauJetContainer> TrigTauMonitorBaseAlgorithm::m_hltTauJetLLPKey {this, "HLTTauJetLLPKey", "HLT_TrigTauRecMerged_LLP", "HLT tracktwoLLP taujet container key"}
protectedinherited

Definition at line 124 of file TrigTauMonitorBaseAlgorithm.h.

124{this, "HLTTauJetLLPKey", "HLT_TrigTauRecMerged_LLP", "HLT tracktwoLLP taujet container key"};

◆ m_hltTauJetLRTKey

SG::ReadHandleKey<xAOD::TauJetContainer> TrigTauMonitorBaseAlgorithm::m_hltTauJetLRTKey {this, "HLTTauJetLRTKey", "HLT_TrigTauRecMerged_LRT", "HLT trackLRT taujet container key"}
protectedinherited

Definition at line 125 of file TrigTauMonitorBaseAlgorithm.h.

125{this, "HLTTauJetLRTKey", "HLT_TrigTauRecMerged_LRT", "HLT trackLRT taujet container key"};

◆ m_L1_Phase1_threshold_patterns

Gaudi::Property<std::map<std::string, uint64_t> > TrigTauMonitorBaseAlgorithm::m_L1_Phase1_threshold_patterns {this, "L1Phase1ThresholdPatterns", {}}
protectedinherited

Definition at line 42 of file TrigTauMonitorBaseAlgorithm.h.

42{this, "L1Phase1ThresholdPatterns", {}};

◆ m_L1_Phase1_thresholds

Gaudi::Property<std::map<std::string, float> > TrigTauMonitorBaseAlgorithm::m_L1_Phase1_thresholds {this, "L1Phase1Thresholds", {}}
protectedinherited

Definition at line 41 of file TrigTauMonitorBaseAlgorithm.h.

41{this, "L1Phase1Thresholds", {}};

◆ m_L1_select_by_et_only

Gaudi::Property<bool> TrigTauMonitorBaseAlgorithm::m_L1_select_by_et_only {this, "SelectL1ByETOnly", false}
protectedinherited

Definition at line 43 of file TrigTauMonitorBaseAlgorithm.h.

43{this, "SelectL1ByETOnly", false};

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

std::map<std::string, std::map<std::string, std::pair<SG::ReadDecorHandleKey<xAOD::TauJetContainer>, SG::ReadDecorHandleKey<xAOD::TauJetContainer> > > > TrigTauMonitorSingleAlgorithm::m_monitoredHLTIdDecorHandleKeys
private

Definition at line 32 of file TrigTauMonitorSingleAlgorithm.h.

◆ m_monitoredHLTIdScores

Gaudi::Property<std::map<std::string, std::map<std::string, std::pair<std::string, std::string> > > > TrigTauMonitorSingleAlgorithm::m_monitoredHLTIdScores {this, "HLTTauIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each HLT TauID algorithm to be monitored, for each reconstruction sequence (type, e.g. tracktwoMVA, tracktwoLLP, etc...)"}
private

Definition at line 31 of file TrigTauMonitorSingleAlgorithm.h.

31{this, "HLTTauIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each HLT TauID algorithm to be monitored, for each reconstruction sequence (type, e.g. tracktwoMVA, tracktwoLLP, etc...)"};

◆ m_monitoredOfflineIdDecorHandleKeys

std::map<std::string, std::pair<SG::ReadDecorHandleKey<xAOD::TauJetContainer>, SG::ReadDecorHandleKey<xAOD::TauJetContainer> > > TrigTauMonitorSingleAlgorithm::m_monitoredOfflineIdDecorHandleKeys
private

Definition at line 36 of file TrigTauMonitorSingleAlgorithm.h.

◆ m_monitoredOfflineIdScores

Gaudi::Property<std::map<std::string, std::pair<std::string, std::string> > > TrigTauMonitorSingleAlgorithm::m_monitoredOfflineIdScores {this, "OfflineTauIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each Offline TauID algorithm to be monitored"}
private

Definition at line 35 of file TrigTauMonitorSingleAlgorithm.h.

35{this, "OfflineTauIDScores", {}, "Pairs of the TauID score and signal-transformed scores for each Offline TauID algorithm to be monitored"};

◆ m_name

std::string AthMonitorAlgorithm::m_name
privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

◆ m_offline_tau_id

Gaudi::Property<unsigned int> TrigTauMonitorSingleAlgorithm::m_offline_tau_id {this, "OfflineTauID", TauID::RNN, "Offline TauID (1: RNN, 2: GNTau)"}
private

Definition at line 25 of file TrigTauMonitorSingleAlgorithm.h.

25{this, "OfflineTauID", TauID::RNN, "Offline TauID (1: RNN, 2: GNTau)"};

◆ m_offlineGNTauDecorKey

SG::ReadDecorHandleKey<xAOD::TauJetContainer> TrigTauMonitorBaseAlgorithm::m_offlineGNTauDecorKey {this, "OfflineGNTauDecorKey", "", "Offline GNTau decoration key"}
privateinherited

Definition at line 137 of file TrigTauMonitorBaseAlgorithm.h.

137{this, "OfflineGNTauDecorKey", "", "Offline GNTau decoration key"};

◆ m_offlineTauJetKey

SG::ReadHandleKey<xAOD::TauJetContainer> TrigTauMonitorBaseAlgorithm::m_offlineTauJetKey {this, "OfflineTauJetKey", "TauJets", "Offline taujet container key"}
protectedinherited

Definition at line 121 of file TrigTauMonitorBaseAlgorithm.h.

121{this, "OfflineTauJetKey", "TauJets", "Offline taujet container key"};

◆ m_phase1l1cTauRoIDecorKey

SG::ReadDecorHandleKey<xAOD::eFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1cTauRoIDecorKey {this, "Phase1L1cTauRoIjTauRoILinkKey", "L1_cTauRoI.jTauLink", "Decoration for the link from eTau to the matching jTau"}
privateinherited

Definition at line 144 of file TrigTauMonitorBaseAlgorithm.h.

144{this, "Phase1L1cTauRoIjTauRoILinkKey", "L1_cTauRoI.jTauLink", "Decoration for the link from eTau to the matching jTau"};

◆ m_phase1l1cTauRoIKey

SG::ReadHandleKey<xAOD::eFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1cTauRoIKey {this, "Phase1L1cTauRoIKey", "L1_cTauRoI", "cTau Phase1 L1 RoI key"}
privateinherited

Definition at line 143 of file TrigTauMonitorBaseAlgorithm.h.

143{this, "Phase1L1cTauRoIKey", "L1_cTauRoI", "cTau Phase1 L1 RoI key"};

◆ m_phase1l1cTauRoIThresholdPatternsKey

SG::ReadDecorHandleKey<xAOD::eFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1cTauRoIThresholdPatternsKey {this, "Phase1L1cTauRoIThresholdPatternsKey", "L1_cTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the cTau RoIs"}
privateinherited

Definition at line 145 of file TrigTauMonitorBaseAlgorithm.h.

145{this, "Phase1L1cTauRoIThresholdPatternsKey", "L1_cTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the cTau RoIs"};

◆ m_phase1l1eTauRoIKey

SG::ReadHandleKey<xAOD::eFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1eTauRoIKey {this, "Phase1L1eTauRoIKey", "L1_eTauRoI", "eTau Phase1 L1 RoI key"}
privateinherited

Definition at line 139 of file TrigTauMonitorBaseAlgorithm.h.

139{this, "Phase1L1eTauRoIKey", "L1_eTauRoI", "eTau Phase1 L1 RoI key"};

◆ m_phase1l1eTauRoIThresholdPatternsKey

SG::ReadDecorHandleKey<xAOD::eFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1eTauRoIThresholdPatternsKey {this, "Phase1L1eTauRoIThresholdPatternsKey", "L1_eTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the eTau RoIs"}
privateinherited

Definition at line 140 of file TrigTauMonitorBaseAlgorithm.h.

140{this, "Phase1L1eTauRoIThresholdPatternsKey", "L1_eTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the eTau RoIs"};

◆ m_phase1l1jTauRoIKey

SG::ReadHandleKey<xAOD::jFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1jTauRoIKey {this, "Phase1L1jTauRoIKey", "L1_jFexTauRoI", "jTau Phase1 L1 RoI key"}
privateinherited

Definition at line 141 of file TrigTauMonitorBaseAlgorithm.h.

141{this, "Phase1L1jTauRoIKey", "L1_jFexTauRoI", "jTau Phase1 L1 RoI key"};

◆ m_phase1l1jTauRoIThresholdPatternsKey

SG::ReadDecorHandleKey<xAOD::jFexTauRoIContainer> TrigTauMonitorBaseAlgorithm::m_phase1l1jTauRoIThresholdPatternsKey {this, "Phase1L1jTauRoIThresholdPatternsKey", "L1_jFexTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the jTau RoIs"}
privateinherited

Definition at line 142 of file TrigTauMonitorBaseAlgorithm.h.

142{this, "Phase1L1jTauRoIThresholdPatternsKey", "L1_jFexTauRoI.thresholdPatterns", "Decoration for the threshold patterns for the jTau RoIs"};

◆ m_requireOfflineTaus

Gaudi::Property<bool> TrigTauMonitorSingleAlgorithm::m_requireOfflineTaus {this, "RequireOfflineTaus", true, "Require at leat 1 offline tau per event"}
private

Definition at line 24 of file TrigTauMonitorSingleAlgorithm.h.

24{this, "RequireOfflineTaus", true, "Require at leat 1 offline tau per event"};

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

Gaudi::Property<std::vector<std::string> > TrigTauMonitorBaseAlgorithm::m_triggers {this, "TriggerList", {}}
protectedinherited

Definition at line 38 of file TrigTauMonitorBaseAlgorithm.h.

38{this, "TriggerList", {}};

◆ m_trigInfo

std::map<std::string, TrigTauInfo> TrigTauMonitorBaseAlgorithm::m_trigInfo
privateinherited

Definition at line 131 of file TrigTauMonitorBaseAlgorithm.h.

◆ m_trigLiveFractionDataKey

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

Definition at line 352 of file AthMonitorAlgorithm.h.

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.


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