ATLAS Offline Software
Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | Private Types | Private Member Functions | Private Attributes | List of all members
TrigEgammaMonitorElectronAlgorithm Class Reference

#include <TrigEgammaMonitorElectronAlgorithm.h>

Inheritance diagram for TrigEgammaMonitorElectronAlgorithm:
Collaboration diagram for TrigEgammaMonitorElectronAlgorithm:

Public Types

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

Public Member Functions

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

Protected Member Functions

void fillLabel (const ToolHandle< GenericMonitoringTool > &groupHandle, const std::string &histname, const std::string &label) const
 
void fillEfficiencies (const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * > > &, const TrigInfo &) const
 
void fillDistributions (const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * > > &, const TrigInfo &) const
 
void fillResolutions (const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const TrigInfo &info) const
 ***********************************************************************************************************************8 More...
 
void fillInefficiency (const std::string &pidword, const TrigInfo &info, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const std::vector< asg::AcceptData > &acceptObjs) const
 
const std::map< std::string, TrigInfo > & getTrigInfoMap ()
 Helper methods. More...
 
bool ApplyElectronPid (const xAOD::Electron *eg, const std::string &) const
 Get offline electron decision. More...
 
bool ApplyPhotonPid (const xAOD::Photon *eg, const std::string &) const
 Get offline electron decision. More...
 
const ToolHandle< Trig::TrigDecisionTool > & tdt () const
 Get the TDT
More...
 
const ToolHandle< TrigEgammaMatchingToolMT > & match () const
 Get the e/g match tool. More...
 
asg::AcceptData setAccept (const TrigCompositeUtils::Decision *, const TrigInfo &) const
 Set the accept object for all trigger levels. More...
 
TrigInfo getTrigInfo (const std::string &) const
 Get the trigger info parsed from the chain name (only single lepton triggers) More...
 
float dR (const float, const float, const float, const float) const
 Get delta R. More...
 
std::string getL1Item (const std::string &trigger) const
 Creates static map to return L1 item from trigger name. More...
 
bool isIsolated (const xAOD::Electron *, const std::string &) const
 Check if electron fulfils isolation criteria. More...
 
bool isPrescaled (const std::string &) const
 Check if the event is prescaled. More...
 
void setTrigInfo (const std::string &)
 Set the trigger info parsed from the chain name. More...
 
bool isHLTTruncated () const
 
float getEta2 (const xAOD::Egamma *eg) const
 Features helper. More...
 
float getEt (const xAOD::Electron *eg) const
 
float getEtCluster37 (const xAOD::Egamma *eg) const
 
float getDEmaxs1 (const xAOD::Egamma *eg) const
 
float rTRT (const xAOD::Electron *eg) const
 
float getSigmaD0 (const xAOD::Electron *eg) const
 
float getD0sig (const xAOD::Electron *eg) const
 
float getEnergyBE0 (const xAOD::Egamma *eg) const
 
float getEnergyBE1 (const xAOD::Egamma *eg) const
 
float getEnergyBE2 (const xAOD::Egamma *eg) const
 
float getEnergyBE3 (const xAOD::Egamma *eg) const
 
float getEaccordion (const xAOD::Egamma *eg) const
 
float getE0Eaccordion (const xAOD::Egamma *eg) const
 
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
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. More...
 

Protected Attributes

ToolHandle< TrigEgammaMatchingToolMT > m_matchTool
 
ToolHandle< Trig::TrigEgammaEmulationToolMT > m_emulatorTool
 
ToolHandleArray< IAsgElectronIsEMSelector > m_electronIsEMTool {this,"ElectronIsEMSelector",{}}
 Offline isEM Selectors. More...
 
ToolHandleArray< IAsgElectronLikelihoodTool > m_electronLHTool {this,"ElectronLikelihoodTool",{}}
 Offline LH Selectors. More...
 
ToolHandleArray< IAsgElectronLikelihoodTool > m_electronDNNTool { this, "ElectronDNNSelectorTool", {},"DNN tools" }
 Offline DNN Selectors. More...
 
ToolHandleArray< IAsgPhotonIsEMSelector > m_photonIsEMTool {this,"PhotonIsEMSelector",{}}
 Offline isEM Photon Selectors. More...
 
Gaudi::Property< bool > m_doEmulation {this, "DoEmulation", false }
 Do emulation. More...
 
Gaudi::Property< bool > m_tp {this, "TPTrigger", false }
 TP Trigger Analysis. More...
 
Gaudi::Property< std::string > m_defaultProbePidElectron {this, "DefaultProbeSelectionElectron", "lhloose"}
 default probe pid for electron trigitems that don't have pid in their name More...
 
Gaudi::Property< std::string > m_defaultProbePidPhoton {this, "DefaultProbeSelectionPhoton", "loose"}
 default probe pid for photon trigitems that don't have pid in their name More...
 
Gaudi::Property< std::vector< std::string > > m_isemname {this, "isEMResultNames", {} }
 isem names More...
 
Gaudi::Property< std::vector< std::string > > m_lhname {this, "LHResultNames", {} }
 lh names More...
 
Gaudi::Property< std::vector< std::string > > m_dnnname {this, "DNNResultNames", {}, }
 dnn names More...
 
Gaudi::Property< bool > m_detailedHists {this, "DetailedHistograms", false}
 Include more detailed histograms. More...
 
asg::AcceptInfo m_accept
 
ToolHandleArray< GenericMonitoringTool > m_tools {this,"GMTools",{}}
 Array of Generic Monitoring Tools. More...
 
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool {this, "TrigDecisionTool",""}
 Tool to tell whether a specific trigger is passed. More...
 
ToolHandleArray< IDQFilterTool > m_DQFilterTools {this,"FilterTools",{}}
 Array of Data Quality filter tools. More...
 
SG::ReadCondHandleKey< LuminosityCondData > m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
 
SG::ReadCondHandleKey< LBDurationCondData > m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
 
SG::ReadCondHandleKey< TrigLiveFractionCondData > m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
 
AthMonitorAlgorithm::Environment_t m_environment
 Instance of the Environment_t enum. More...
 
AthMonitorAlgorithm::DataType_t m_dataType
 Instance of the DataType_t enum. More...
 
Gaudi::Property< std::string > m_environmentStr {this,"Environment","user"}
 Environment string pulled from the job option and converted to enum. More...
 
Gaudi::Property< std::string > m_dataTypeStr {this,"DataType","userDefined"}
 DataType string pulled from the job option and converted to enum. More...
 
Gaudi::Property< std::string > m_triggerChainString {this,"TriggerChain",""}
 Trigger chain string pulled from the job option and parsed into a vector. More...
 
std::vector< std::string > m_vTrigChainNames
 Vector of trigger chain names parsed from trigger chain string. More...
 
Gaudi::Property< std::string > m_fileKey {this,"FileKey",""}
 Internal Athena name for file. More...
 
Gaudi::Property< bool > m_useLumi {this,"EnableLumi",false}
 Allows use of various luminosity functions. More...
 
Gaudi::Property< float > m_defaultLBDuration {this,"DefaultLBDuration",60.}
 Default duration of one lumi block. More...
 
Gaudi::Property< int > m_detailLevel {this,"DetailLevel",0}
 Sets the level of detail used in the monitoring. More...
 
SG::ReadHandleKey< xAOD::EventInfo > m_EventInfoKey {this,"EventInfoKey","EventInfo"}
 Key for retrieving EventInfo from StoreGate. More...
 

Private Types

typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > MonVarVec_t
 
typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

StatusCode executeNavigation (const EventContext &ctx, const std::string &trigItem, float, const std::string &, std::vector< std::pair< std::shared_ptr< const xAOD::Egamma >, const TrigCompositeUtils::Decision * >> &) const
 navigation method called by executeNavigation More...
 
void fillEfficiency (const std::string &subgroup, const std::string &level, const std::string &pidword, const TrigInfo &info, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const std::vector< asg::AcceptData > &acceptObjs, const std::string &dirname) const
 
void fillL1Calo (const std::string &trigger, const std::vector< const xAOD::EmTauRoI * > &l1_vec) const
 
void fillL1eEM (const std::string &trigger, const std::vector< const xAOD::eFexEMRoI * > &l1_vec) const
 
void fillL2Calo (const std::string &trigger, const std::vector< const xAOD::TrigEMCluster * > &emCluster_vec) const
 
void fillL2Electron (const std::string &trigger, const std::vector< const xAOD::TrigElectron * > &el_vec) const
 
void fillL2Photon (const std::string &trigger, const std::vector< const xAOD::TrigPhoton * > &eg_vec) const
 
void fillEFCalo (const std::string &trigger, const std::vector< const xAOD::CaloCluster * > &clus_vec) const
 
void fillShowerShapes (const std::string &trigger, const std::vector< const xAOD::Egamma * > &eg_vec, bool online) const
 
void fillTracking (const std::string &trigger, const std::vector< const xAOD::Electron * > &eg_vec, bool online) const
 
void fillL1CaloResolution (const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const xAOD::EmTauRoI * >> &pairObjs) const
 
void fillL1CaloAbsResolution (const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const xAOD::EmTauRoI * >> &pairObjs) const
 
void fillL2CaloResolution (const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs) const
 
void fillHLTElectronResolution (const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const TrigInfo &info) const
 
void fillHLTPhotonResolution (const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const TrigInfo &info) const
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

std::vector< std::string > m_trigList
 List of triggers to study. More...
 
Gaudi::Property< std::vector< std::string > > m_trigInputList {this, "TriggerList", {}}
 List of triggers from menu. More...
 
Gaudi::Property< std::string > m_anatype {this, "Analysis", "Electron"}
 Directory name for each algorithm. More...
 
Gaudi::Property< std::string > m_offProbeIsolation { this, "OfflineProbeIsolation", "Loose"}
 Define isolation working point for Probe electron. More...
 
Gaudi::Property< bool > m_forceProbeIsolation {this, "ForceProbeIsolation", false}
 force probe isolation More...
 
Gaudi::Property< bool > m_forcePidSelection { this, "ForcePidSelection", true}
 force pid and crack selection into electron navigation More...
 
Gaudi::Property< bool > m_forceEtThr { this, "ForceEtThreshold", true}
 
Gaudi::Property< bool > m_forceVetoVeryLoose { this, "ForceVetoVeryLoose", false}
 
Gaudi::Property< bool > m_rmCrack { this, "RemoveCrack", true}
 Remove crack region for Probe default True. More...
 
SG::ReadHandleKey< xAOD::ElectronContainer > m_offElectronKey { this, "ElectronKey", "Electrons", ""}
 Event Wise offline ElectronContainer Access and end iterator. More...
 
SG::ReadDecorHandleKey< xAOD::EventInfo > m_eventInfoDecorKey {this, "LArStatusFlag", "EventInfo.larFlags", "Key for EventInfo object"}
 
SG::ReadDecorHandleKey< xAOD::ElectronContainer > m_electronIsolationKeyPtCone20 { this, "PtCone20Decoration", m_offElectronKey, "ptcone20", "Decoration key for the ptcone20 isolation decoration" }
 Ensure offline electron isolation decoration is retrieved after being created. More...
 
ToolHandle< Trig::TrigDecisionTool > m_trigdec
 Trigger decision tool. More...
 
std::map< std::string, TrigInfo > m_trigInfo
 creates map of trigger name and TrigInfo struct More...
 
std::string m_name
 
std::unordered_map< std::string, size_t > m_toolLookupMap
 
const ToolHandle< GenericMonitoringTool > m_dummy
 
Gaudi::Property< bool > m_enforceExpressTriggers
 
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 14 of file TrigEgammaMonitorElectronAlgorithm.h.

Member Typedef Documentation

◆ MonVarVec_t

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

Definition at line 365 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 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 191 of file AthMonitorAlgorithm.h.

191  {
192  userDefined = 0,
193  monteCarlo,
194  collisions,
195  cosmics,
196  heavyIonCollisions,
197  };

◆ Environment_t

enum 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 172 of file AthMonitorAlgorithm.h.

172  {
173  user = 0,
174  online,
175  tier0,
176  tier0Raw,
177  tier0ESD,
178  AOD,
179  altprod,
180  };

Constructor & Destructor Documentation

◆ TrigEgammaMonitorElectronAlgorithm()

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

Definition at line 13 of file TrigEgammaMonitorElectronAlgorithm.cxx.

13  :
14  TrigEgammaMonitorAnalysisAlgorithm( name, pSvcLocator )
15 {}

◆ ~TrigEgammaMonitorElectronAlgorithm()

TrigEgammaMonitorElectronAlgorithm::~TrigEgammaMonitorElectronAlgorithm ( )
overridevirtual

Definition at line 17 of file TrigEgammaMonitorElectronAlgorithm.cxx.

18 {}

Member Function Documentation

◆ ApplyElectronPid()

bool TrigEgammaMonitorBaseAlgorithm::ApplyElectronPid ( const xAOD::Electron *  eg,
const std::string &  pidname 
) const
protectedinherited

Get offline electron decision.

Definition at line 49 of file TrigEgammaMonitorBaseAlgorithm.cxx.

50 {
51  const auto& ctx = Gaudi::Hive::currentContext() ;
52  if (pidname == "tight"){
53  return (bool) this->m_electronIsEMTool[0]->accept(ctx,eg);
54  }
55  else if (pidname == "medium"){
56  return (bool) this->m_electronIsEMTool[1]->accept(ctx,eg);
57  }
58  else if (pidname == "loose"){
59  return (bool) this->m_electronIsEMTool[2]->accept(ctx,eg);
60  }
61  else if (pidname == "lhtight"){
62  return (bool) this->m_electronLHTool[0]->accept(ctx,eg);
63  }
64  else if (pidname == "lhmedium"){
65  return (bool) this->m_electronLHTool[1]->accept(ctx,eg);
66  }
67  else if (pidname == "lhloose"){
68  return (bool) this->m_electronLHTool[2]->accept(ctx,eg);
69  }
70  else if (pidname == "lhvloose"){
71  return (bool) this->m_electronLHTool[3]->accept(ctx,eg);
72  }
73  else if (pidname == "dnntight"){
74  return (bool) this->m_electronDNNTool[0]->accept(ctx,eg);
75  }
76  else if (pidname == "dnnmedium"){
77  return (bool) this->m_electronDNNTool[1]->accept(ctx,eg);
78  }
79  else if (pidname == "dnnloose"){
80  return (bool) this->m_electronDNNTool[2]->accept(ctx,eg);
81  }
82  else ATH_MSG_DEBUG("No Pid tool, continue without PID");
83  return false;
84 }

◆ ApplyPhotonPid()

bool TrigEgammaMonitorBaseAlgorithm::ApplyPhotonPid ( const xAOD::Photon *  eg,
const std::string &  pidname 
) const
protectedinherited

Get offline electron decision.

Definition at line 91 of file TrigEgammaMonitorBaseAlgorithm.cxx.

92 {
93  const auto& ctx = Gaudi::Hive::currentContext() ;
94  if (pidname == "tight"){
95  return (bool) this->m_photonIsEMTool[0]->accept(ctx,eg);
96  }
97  else if (pidname == "medium"){
98  return (bool) this->m_photonIsEMTool[1]->accept(ctx,eg);
99  }
100  else if (pidname == "loose"){
101  return (bool) this->m_photonIsEMTool[2]->accept(ctx,eg);
102  }
103  else ATH_MSG_DEBUG("No Pid tool, continue without PID");
104  return false;
105 }

◆ cardinality()

unsigned int AthReentrantAlgorithm::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.

Override this to return 0 for reentrant algorithms.

Definition at line 55 of file AthReentrantAlgorithm.cxx.

56 {
57  return 0;
58 }

◆ 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 221 of file AthMonitorAlgorithm.h.

221 { return m_dataType; }

◆ dataTypeStringToEnum()

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

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

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

Definition at line 140 of file AthMonitorAlgorithm.cxx.

140  {
141  // convert the string to all lowercase
142  std::string lowerCaseStr = str;
143  std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
144 
145  // check if it matches one of the enum choices
146  if( lowerCaseStr == "userdefined" ) {
147  return DataType_t::userDefined;
148  } else if( lowerCaseStr == "montecarlo" ) {
149  return DataType_t::monteCarlo;
150  } else if( lowerCaseStr == "collisions" ) {
151  return DataType_t::collisions;
152  } else if( lowerCaseStr == "cosmics" ) {
153  return DataType_t::cosmics;
154  } else if( lowerCaseStr == "heavyioncollisions" ) {
155  return DataType_t::heavyIonCollisions;
156  } else { // otherwise, warn the user and return "userDefined"
157  ATH_MSG_WARNING("AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
158  <<str<<", returning userDefined.");
159  return DataType_t::userDefined;
160  }
161 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType &   
)
inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType &   
)
inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType &   
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleBase &  hndl,
const std::string &  doc,
const SG::VarHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKey &  hndl,
const std::string &  doc,
const SG::VarHandleKeyType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray &  hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType &   
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

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

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

95 { return m_detStore; }

◆ dR()

float TrigEgammaMonitorBaseAlgorithm::dR ( const float  eta1,
const float  phi1,
const float  eta2,
const float  phi2 
) const
protectedinherited

Get delta R.

Definition at line 272 of file TrigEgammaMonitorBaseAlgorithm.cxx.

272  {
273  float deta = fabs(eta1 - eta2);
274  float dphi = fabs(phi1 - phi2) < TMath::Pi() ? fabs(phi1 - phi2) : 2*TMath:: \
275  Pi() - fabs(phi1 - phi2);
276  return sqrt(deta*deta + dphi*dphi);
277 }

◆ 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 205 of file AthMonitorAlgorithm.h.

205 { return m_environment; }

◆ envStringToEnum()

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

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

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

Definition at line 112 of file AthMonitorAlgorithm.cxx.

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

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ 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 73 of file AthMonitorAlgorithm.cxx.

73  {
74 
75  // Checks that all of the DQ filters are passed. If any one of the filters
76  // fails, return SUCCESS code and do not fill the histograms with the event.
77  for ( const auto& filterItr : m_DQFilterTools ) {
78  if (!filterItr->accept()) {
79  ATH_MSG_DEBUG("Event rejected due to filter tool.");
80  return StatusCode::SUCCESS;
81  }
82  }
83 
84  // Trigger: If there is a decision tool and the chains fail, skip the event.
85  if ( !m_trigDecTool.empty() && !trigChainsArePassed(m_vTrigChainNames) ) {
86  ATH_MSG_DEBUG("Event rejected due to trigger filter.");
87  return StatusCode::SUCCESS;
88  }
89 
90  ATH_MSG_DEBUG("Event accepted!");
91  return fillHistograms(ctx);
92 }

◆ executeNavigation()

StatusCode TrigEgammaMonitorElectronAlgorithm::executeNavigation ( const EventContext &  ctx,
const std::string &  trigItem,
float  etthr,
const std::string &  pidName,
std::vector< std::pair< std::shared_ptr< const xAOD::Egamma >, const TrigCompositeUtils::Decision * >> &  pairObjs 
) const
private

navigation method called by executeNavigation

default is true

default is true

default is true

default is false

default is Loose

default is false

Definition at line 103 of file TrigEgammaMonitorElectronAlgorithm.cxx.

108 {
109 
110  ATH_MSG_DEBUG("Apply navigation selection ");
111 
112  SG::ReadHandle<xAOD::ElectronContainer> offElectrons(m_offElectronKey, ctx);
113 
114  if(!offElectrons.isValid())
115  {
116  ATH_MSG_DEBUG("Failed to retrieve offline Electrons ");
117  return StatusCode::FAILURE;
118  }
119 
120 
121  const std::string decorName="is"+pidName;
122  for(const auto *const eg : *offElectrons ){
123 
124  const TrigCompositeUtils::Decision *dec=nullptr;
125 
126  if(!eg->trackParticle()){
127  ATH_MSG_DEBUG("No track Particle");
128  continue;
129  }
130  if(!eg->caloCluster()){
131  ATH_MSG_DEBUG("No caloCluster");
132  continue;
133  }
134 
135  if(m_forceEtThr){
136  if( !( getEt(eg) > (etthr-5.)*1.e3) ) continue;
137  }
138 
139  if(m_rmCrack){
140  if ( (fabs(eg->eta())>1.37 && fabs(eg->eta())<1.52) || fabs(eg->eta())>2.47 ) continue;
141  }
142 
143  if(m_forcePidSelection){
144  if(!ApplyElectronPid(eg,pidName)){
145  ATH_MSG_DEBUG("Fails ElectronID "<< pidName);
146  continue;
147  }
148  ATH_MSG_DEBUG("Passes ElectronID "<< pidName);
149  }
150 
151  if (m_forceProbeIsolation) {
152  if (!isIsolated(eg, m_offProbeIsolation)) continue;
153  }
154 
155  if(m_forceVetoVeryLoose){
156  bool veto = ApplyElectronPid(eg,"lhvloose");
157  if(veto) continue;
158  }
159 
160  const auto el = std::make_shared<const xAOD::Electron>(*eg);
161  SG::Decorator<bool> decor (decorName);
162  decor(*el)=static_cast<bool>(true);
163 
164  match()->match(el.get(), trigItem, dec, TrigDefs::includeFailedDecisions);
165  //match()->match(el, trigItem, dec);
166  pairObjs.emplace_back(el, dec);
167 
168  }
169 
170  ATH_MSG_DEBUG("BaseToolMT::Electron TEs " << pairObjs.size() << " found.");
171  return StatusCode::SUCCESS;
172 }

◆ 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 & AthReentrantAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 79 of file AthReentrantAlgorithm.cxx.

80 {
81  // If we didn't find any symlinks to add, just return the collection
82  // from the base class. Otherwise, return the extended collection.
83  if (!m_extendedExtraObjects.empty()) {
84  return m_extendedExtraObjects;
85  }
86  return Algorithm::extraOutputDeps();
87 }

◆ fillDistributions()

void TrigEgammaMonitorAnalysisAlgorithm::fillDistributions ( const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * > > &  ,
const TrigInfo &   
) const
protectedinherited

Definition at line 407 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

409 {
410 
411  const std::string trigger = info.trigger;
412 
413  unsigned int condition=TrigDefs::Physics;
414 
415  // Offline
416  std::vector<const xAOD::Egamma*> eg_vec;
417  std::vector<const xAOD::Electron*> el_vec;
418  for( auto pairObj: pairObjs )
419  {
420  eg_vec.push_back(pairObj.first);
421  if( xAOD::EgammaHelpers::isElectron(pairObj.first)){
422  const xAOD::Electron* elOff = static_cast<const xAOD::Electron*> (pairObj.first);
423  el_vec.push_back(elOff);
424  }
425  }
426 
427  // Offline
428  fillShowerShapes( trigger, eg_vec, false );
429  fillTracking( trigger, el_vec, false );
430 
431  // L1Calo
432  {
433  // Fill L1 features
434  auto initRois = tdt()->features<TrigRoiDescriptorCollection>(trigger,condition,"",TrigDefs::allFeaturesOfType,"initialRoI");
435 
436  if (info.L1Legacy){
437  std::vector<const xAOD::EmTauRoI*> l1_vec;
438  for( auto &initRoi: initRois ){
439  if( !initRoi.link.isValid() ) continue;
440  const auto *feat = match()->getL1Feature( initRoi.source );
441  if(feat) l1_vec.push_back(feat);
442  }
443  fillL1Calo( trigger, l1_vec );
444  }else{
445  std::vector<const xAOD::eFexEMRoI*> l1_vec;
446  for( auto &initRoi: initRois ){
447  if( !initRoi.link.isValid() ) continue;
448  const auto *feat = match()->getL1eEMFeature( initRoi.source );
449  if(feat) l1_vec.push_back(feat);
450  }
451  fillL1eEM( trigger, l1_vec );
452  }
453 
454  }
455 
456 
457  // L2Calo
458  {
459  std::vector<const xAOD::TrigEMCluster*> emCluster_vec;
460  auto vec = tdt()->features<xAOD::TrigEMClusterContainer>(trigger,condition ,match()->key("FastCalo") );
461  for(auto &featLinkInfo : vec ){
462  if(! featLinkInfo.isValid() ) continue;
463  const auto *feat = *(featLinkInfo.link);
464  if(!feat) continue;
465  emCluster_vec.push_back(feat);
466  }
467  fillL2Calo( trigger, emCluster_vec );
468  }
469  // EFCalo
470  {
471  std::string key = match()->key("PrecisionCalo_Electron");
472  if(info.signature == "Photon") key = match()->key("PrecisionCalo_Photon");
473  if(info.lrt) key = match()->key("PrecisionCalo_LRT");
474  if(info.ion) key = match()->key("PrecisionCalo_HI");
475 
476  std::vector<const xAOD::CaloCluster* > clus_vec;
477  auto vec = tdt()->features<xAOD::CaloClusterContainer>(trigger,condition,key);
478  for(auto &featLinkInfo : vec ){
479  if(! featLinkInfo.isValid() ) continue;
480  const auto *feat = *(featLinkInfo.link);
481  if(!feat) continue;
482  clus_vec.push_back(feat);
483  }
484  fillEFCalo( trigger, clus_vec );
485  }
486 
487  if ( info.signature == "Electron" ){
488 
489  // L2 Electron
490  {
491  std::string key = match()->key("FastElectrons");
492  if(info.lrt) key = match()->key("FastElectrons_LRT");
493 
494  std::vector<const xAOD::TrigElectron*> el_vec;
495  // Get only passed objects
496  auto vec = tdt()->features<xAOD::TrigElectronContainer>(trigger,condition,key );
497  for( auto &featLinkInfo : vec ){
498  if(! featLinkInfo.isValid() ) continue;
499  const auto *feat = *(featLinkInfo.link);
500  if(!feat) continue;
501  el_vec.push_back(feat);
502  }
503  fillL2Electron( trigger, el_vec );
504  }
505  // HLT Electron
506  {
507  std::string key = match()->key("Electrons_GSF");
508  if(info.nogsf) key = match()->key("Electrons");
509  if(info.lrt) key = match()->key("Electrons_LRT");
510 
511  std::vector<const xAOD::Electron*> el_vec;
512  std::vector<const xAOD::Egamma*> eg_vec;
513  auto vec = tdt()->features<xAOD::ElectronContainer>(trigger, condition ,key );
514  for( auto &featLinkInfo : vec ){
515  if(! featLinkInfo.isValid() ) continue;
516  const auto *feat = *(featLinkInfo.link);
517  if(!feat) continue;
518  el_vec.push_back(feat);
519  eg_vec.push_back(feat);
520  }
521  fillShowerShapes( trigger, eg_vec, true );
522  fillTracking( trigger, el_vec, true );
523  }
524  }else if ( info.signature == "Photon"){
525  // Fast Photon
526  {
527  std::string key = match()->key("FastPhotons");
528  std::vector<const xAOD::TrigPhoton*> ph_vec;
529  // Get only passed objects
530  auto vec = tdt()->features<xAOD::TrigPhotonContainer>(trigger,condition,key );
531  for( auto &featLinkInfo : vec ){
532  if(! featLinkInfo.isValid() ) continue;
533  const auto *feat = *(featLinkInfo.link);
534  if(!feat) continue;
535  ph_vec.push_back(feat);
536  }
537  fillL2Photon( trigger, ph_vec );
538  }
539 
540  // HLT Photon
541  {
542  std::vector<const xAOD::Egamma*> ph_vec;
543  auto vec = tdt()->features<xAOD::PhotonContainer>(trigger,condition ,match()->key("Photons") );
544  for( auto &featLinkInfo : vec ){
545  if(! featLinkInfo.isValid() ) continue;
546  const auto *feat = *(featLinkInfo.link);
547  if(!feat) continue;
548  ph_vec.push_back(feat);
549  }
550  fillShowerShapes( trigger, ph_vec, true );
551  }
552  }else{
553  ATH_MSG_INFO( "Chain type not specified" );
554  }
555 
556 }

◆ fillEFCalo()

void TrigEgammaMonitorAnalysisAlgorithm::fillEFCalo ( const std::string &  trigger,
const std::vector< const xAOD::CaloCluster * > &  clus_vec 
) const
privateinherited

Definition at line 698 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

699 {
700 
701  auto monGroup = getGroup( trigger + "_Distributions_EFCalo" );
702 
703 
704  std::vector<float> energyBE0_vec, energyBE1_vec, energyBE2_vec, energyBE3_vec,
705  energy_vec, et_vec, eta_vec, phi_vec, eta_calo_vec, phi_calo_vec, highet_vec;
706 
707 
708 
709  auto energyBE0_col = Monitored::Collection("energyBE0", energyBE0_vec);
710  auto energyBE1_col = Monitored::Collection("energyBE1", energyBE1_vec);
711  auto energyBE2_col = Monitored::Collection("energyBE2", energyBE2_vec);
712  auto energyBE3_col = Monitored::Collection("energyBE3", energyBE3_vec);
713  auto energy_col = Monitored::Collection("energy" , energy_vec );
714  auto et_col = Monitored::Collection("et" , et_vec );
715  auto highet_col = Monitored::Collection("highet" , highet_vec );
716  auto eta_col = Monitored::Collection("eta" , eta_vec );
717  auto phi_col = Monitored::Collection("phi" , phi_vec );
718  auto eta_calo_col = Monitored::Collection("eta_calo" , eta_calo_vec );
719  auto phi_calo_col = Monitored::Collection("phi_calo" , phi_calo_vec );
720 
721  for ( const auto *clus : clus_vec )
722  {
723  double tmpeta = -999.;
724  if(!clus->retrieveMoment(xAOD::CaloCluster::ETACALOFRAME,tmpeta))
725  tmpeta=-999.;
726  double tmpphi = -999.;
727  if(!clus->retrieveMoment(xAOD::CaloCluster::PHICALOFRAME,tmpphi))
728  tmpphi=-999.;
729 
730  energyBE0_vec.push_back( clus->energyBE(0)/Gaudi::Units::GeV );
731  energyBE1_vec.push_back( clus->energyBE(1)/Gaudi::Units::GeV );
732  energyBE2_vec.push_back( clus->energyBE(2)/Gaudi::Units::GeV );
733  energyBE3_vec.push_back( clus->energyBE(3)/Gaudi::Units::GeV );
734  energy_vec.push_back( clus->e()/Gaudi::Units::GeV );
735  et_vec.push_back( clus->et()/Gaudi::Units::GeV );
736  highet_vec.push_back( clus->et()/Gaudi::Units::GeV );
737  eta_vec.push_back( clus->eta() );
738  phi_vec.push_back( clus->phi() );
739  eta_calo_vec.push_back( tmpeta );
740  phi_calo_vec.push_back( tmpphi );
741 
742  }
743 
744 
745  fill( monGroup, energyBE0_col, energyBE1_col, energyBE2_col, energyBE3_col,
746  energy_col, et_col, eta_col, phi_col, eta_calo_col, phi_calo_col, highet_col);
747 }

◆ fillEfficiencies()

void TrigEgammaMonitorAnalysisAlgorithm::fillEfficiencies ( const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * > > &  ,
const TrigInfo &   
) const
protectedinherited

Definition at line 42 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

44 {
45 
46  std::vector< std::pair< const xAOD::Egamma*, const TrigCompositeUtils::Decision * >> pair_vec;
47  std::vector< std::pair< const xAOD::Egamma*, const TrigCompositeUtils::Decision * >> pair_iso_vec;
48 
49  std::vector<asg::AcceptData> accept_vec;
50  std::vector<asg::AcceptData> accept_iso_vec;
51 
52  std::vector<asg::AcceptData> emu_accept_vec;
53  std::vector<asg::AcceptData> emu_accept_iso_vec;
54 
55 
56  for( auto pairObj : pairObjs ){
57 
58  if(pairObj.first->type()==xAOD::Type::Electron){
59 
60  auto passBits=tdt()->isPassedBits(info.trigger);
61  if(!((passBits & TrigDefs::L1_isPassedAfterVeto) && ((passBits & TrigDefs::EF_prescaled)==0))){
62  ATH_MSG_DEBUG("Prescaled trigger: " << info.trigger << " Skipping to normalize efficiencies");
63  continue;
64  }
65  const xAOD::Electron* el = static_cast<const xAOD::Electron *> (pairObj.first);
66  float et = getEt(el)/Gaudi::Units::GeV;
67  if(et < info.etthr-5.0) continue;
68 
69  }else if(pairObj.first->type()==xAOD::Type::Photon){
70  float et = getCluster_et(pairObj.first)/Gaudi::Units::GeV;
71  if(et < info.etthr-5.0) continue;
72 
73  // Applying FixedCutLoose isolation on the offline photon as recommended in the twiki:
74  // https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/IsolationSelectionTool#Photons
75  bool pass_CaloIso = getIsolation_topoetcone20(pairObj.first)/getCluster_et(pairObj.first) < 0.065;
76  bool pass_trkIso = getIsolation_ptcone20(pairObj.first)/getCluster_et(pairObj.first) < 0.05;
77 
78  if (!pass_CaloIso || !pass_trkIso){
79  continue; // pass FixedCutLoose offline isolation
80  }
81  } // Offline photon
82 
83 
84  // Good pair to be measure
85  if(m_doEmulation){ // Emulation
86  bool valid=false;
87  auto acceptData = m_emulatorTool->emulate( pairObj.second, info.trigger , valid);
88  // skip this probe since the emulation is not possible. Avoid diff denominators between emulation and efficiecy
89  if(!valid) {
90  ATH_MSG_DEBUG("Emulation fail. Skip this probe...");
91  continue;
92  }
93  emu_accept_vec.push_back( acceptData );
94  static const SG::Decorator<bool> IsolatedDec("Isolated");
95  if( IsolatedDec(*pairObj.first) ){
96  emu_accept_iso_vec.push_back(acceptData);
97  }
98  }
99 
100  // Good pair to be measure
101  { // Efficiency
102  pair_vec.push_back(pairObj);
103  auto acceptData = setAccept( pairObj.second, info );
104  accept_vec.push_back(acceptData);
105  static const SG::Decorator<bool> IsolatedDec("Isolated");
106  if( IsolatedDec(*pairObj.first) ){
107  pair_iso_vec.push_back(pairObj);
108  accept_iso_vec.push_back(acceptData);
109  }
110  }
111 
112 
113  }
114 
115  std::string dirname= "Efficiency";
116  std::string l2step = "FastElectron";
117  if( info.signature == "Electron" ){
118  l2step = "FastElectron";
119  }else if( info.signature == "Photon" ){
120  l2step = "FastPhoton";
121  }
122  fillEfficiency( "L1Calo" , "L1Calo" , info.pidname, info, pair_vec , accept_vec, dirname);
123  fillEfficiency( "FastCalo" , "L2Calo" , info.pidname, info, pair_vec , accept_vec, dirname);
124  fillEfficiency( l2step , "L2" , info.pidname, info, pair_vec , accept_vec, dirname);
125  fillEfficiency( "PrecisionCalo" , "EFCalo" , info.pidname, info, pair_vec , accept_vec, dirname);
126  fillEfficiency( "HLT" , "HLT" , info.pidname, info, pair_vec , accept_vec, dirname);
127 
128 
129  if( m_detailedHists ){
130  for( const auto& pid : m_isemname ){
131  fillEfficiency( "HLT_" + pid, "HLT", pid, info, pair_vec , accept_vec, dirname);
132  fillEfficiency( "HLT_" + pid + "Iso", "HLT", pid, info, pair_iso_vec, accept_iso_vec, dirname );
133  }
134  for( const auto& pid : m_lhname ){
135  fillEfficiency( "HLT_" + pid, "HLT", pid, info, pair_vec, accept_vec, dirname );
136  fillEfficiency( "HLT_" + pid + "Iso", "HLT", pid, info, pair_iso_vec, accept_iso_vec, dirname );
137  }
138  }
139 
140  // Fill emulator efficiency plots
141  if ( m_doEmulation ){
142  dirname= "Emulation";
143  fillEfficiency( "L1Calo" , "L1Calo" , info.pidname, info, pair_vec , emu_accept_vec, dirname);
144  fillEfficiency( "FastCalo" , "L2Calo" , info.pidname, info, pair_vec , emu_accept_vec, dirname);
145  fillEfficiency( l2step , "L2" , info.pidname, info, pair_vec , emu_accept_vec, dirname);
146  fillEfficiency( "PrecisionCalo" , "EFCalo" , info.pidname, info, pair_vec , emu_accept_vec, dirname);
147  fillEfficiency( "HLT" , "HLT" , info.pidname, info, pair_vec , emu_accept_vec, dirname);
148  }
149 
150  // Fill Inefficiencies
151  fillInefficiency( info.pidname, info, pair_vec , accept_vec);
152 
153 }

◆ fillEfficiency()

void TrigEgammaMonitorAnalysisAlgorithm::fillEfficiency ( const std::string &  subgroup,
const std::string &  level,
const std::string &  pidword,
const TrigInfo &  info,
const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &  pairObjs,
const std::vector< asg::AcceptData > &  acceptObjs,
const std::string &  dirname 
) const
privateinherited

Definition at line 156 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

163 {
164  const float etthr = info.etthr;
165  const std::string trigger = info.trigger;
166  auto monGroup = getGroup( trigger + "_"+dirname+"_" + subgroup );
167 
168  std::vector<float> et_vec, highet_vec, pt_vec, eta_vec, phi_vec, avgmu_vec, npvtx_vec,et_slice0_vec,et_slice1_vec,et_slice2_vec,et_slice3_vec, ptvarcone20rel_vec, z0_vec, d0_vec;
169  std::vector<float> match_et_vec, match_highet_vec, match_pt_vec, match_eta_vec, match_phi_vec, match_avgmu_vec, match_npvtx_vec, match_ptvarcone20rel_vec, match_z0_vec, match_d0_vec;
170  std::vector<bool> et_passed_vec, et_failed_vec, highet_passed_vec, highet_failed_vec, pt_passed_vec, eta_passed_vec, eta_failed_vec, phi_passed_vec, avgmu_passed_vec, npvtx_passed_vec, ptvarcone20rel_passed_vec, z0_passed_vec, d0_passed_vec;
171  std::vector<bool> et_slice0_passed_vec,et_slice1_passed_vec,et_slice2_passed_vec,et_slice3_passed_vec;
172 
173  auto et_col = Monitored::Collection( "et" , et_vec );
174  auto highet_col = Monitored::Collection( "highet" , highet_vec );
175  auto pt_col = Monitored::Collection( "pt" , pt_vec );
176  auto eta_col = Monitored::Collection( "eta" , eta_vec );
177  auto phi_col = Monitored::Collection( "phi" , phi_vec );
178  auto avgmu_col = Monitored::Collection( "avgmu" , avgmu_vec );
179  auto npvtx_col = Monitored::Collection( "npvtx" , npvtx_vec );
180  auto ptvarcone20rel_col = Monitored::Collection( "ptvarcone20rel" , ptvarcone20rel_vec );
181  auto z0_col = Monitored::Collection( "z0" , z0_vec );
182  auto d0_col = Monitored::Collection( "d0" , d0_vec );
183 
184  auto match_et_col = Monitored::Collection( "match_et" , match_et_vec );
185  auto match_highet_col = Monitored::Collection( "match_ethigh" , match_highet_vec );
186  auto match_pt_col = Monitored::Collection( "match_pt" , match_pt_vec );
187  auto match_eta_col = Monitored::Collection( "match_eta" , match_eta_vec );
188  auto match_phi_col = Monitored::Collection( "match_phi" , match_phi_vec );
189  auto match_avgmu_col = Monitored::Collection( "match_avgmu" , match_avgmu_vec );
190  auto match_npvtx_col = Monitored::Collection( "match_npvtx" , match_npvtx_vec );
191  auto match_ptvarcone20rel_col = Monitored::Collection( "match_ptvarcone20rel" , match_ptvarcone20rel_vec );
192  auto match_z0_col = Monitored::Collection( "match_z0" , match_z0_vec );
193  auto match_d0_col = Monitored::Collection( "match_d0" , match_d0_vec );
194 
195  auto et_passed_col = Monitored::Collection( "et_passed" , et_passed_vec );
196  auto et_failed_col = Monitored::Collection( "et_failed" , et_failed_vec );
197  auto highet_passed_col = Monitored::Collection( "highet_passed" , highet_passed_vec );
198  auto highet_failed_col = Monitored::Collection( "highet_failed" , highet_failed_vec );
199  auto pt_passed_col = Monitored::Collection( "pt_passed" , pt_passed_vec );
200  auto eta_passed_col = Monitored::Collection( "eta_passed" , eta_passed_vec );
201  auto eta_failed_col = Monitored::Collection( "eta_failed" , eta_failed_vec );
202  auto phi_passed_col = Monitored::Collection( "phi_passed" , phi_passed_vec );
203  auto avgmu_passed_col = Monitored::Collection( "avgmu_passed" , avgmu_passed_vec );
204  auto npvtx_passed_col = Monitored::Collection( "npvtx_passed" , npvtx_passed_vec );
205  auto ptvarcone20rel_passed_col = Monitored::Collection( "ptvarcone20rel_passed" , ptvarcone20rel_passed_vec );
206  auto z0_passed_col = Monitored::Collection( "z0_passed" , z0_passed_vec );
207  auto d0_passed_col = Monitored::Collection( "d0_passed" , d0_passed_vec );
208 
209  // For ET efficiency analysis in eta slices
210  auto et_slice0_col = Monitored::Collection( "et_slice0" , et_slice0_vec );
211  auto et_slice1_col = Monitored::Collection( "et_slice1" , et_slice1_vec );
212  auto et_slice2_col = Monitored::Collection( "et_slice2" , et_slice2_vec );
213  auto et_slice3_col = Monitored::Collection( "et_slice3" , et_slice3_vec );
214 
215  auto et_slice0_passed_col = Monitored::Collection( "et_slice0_passed" , et_slice0_passed_vec );
216  auto et_slice1_passed_col = Monitored::Collection( "et_slice1_passed" , et_slice1_passed_vec );
217  auto et_slice2_passed_col = Monitored::Collection( "et_slice2_passed" , et_slice2_passed_vec );
218  auto et_slice3_passed_col = Monitored::Collection( "et_slice3_passed" , et_slice3_passed_vec );
219 
220  unsigned iObj=0;
221 
222  SG::Decorator<bool> pidwordDec("is"+pidword);
223  for( auto pairObj : pairObjs ){
224 
225  bool pid=true;
226  bool isPassed = acceptObjs[iObj].getCutResult( level );
227  float et=0.;
228  float ptvarcone20rel = -99.0;
229  float z0 = -99.0;
230  float d0 = -99.0;
231  const auto *eg = pairObj.first;
232  ATH_MSG_DEBUG("Default pid " << pid << " te " << isPassed);
233  if(xAOD::EgammaHelpers::isElectron(eg)){
234  ATH_MSG_DEBUG("Offline Electron with pidword " << pidword);
235  const xAOD::Electron* el =static_cast<const xAOD::Electron*> (eg);
236  pid=pidwordDec(*el);
237  ATH_MSG_DEBUG("Electron pid " << pid);
238  et = getEt(el)/Gaudi::Units::GeV;
239  if (el->pt() > 0) {
240  ptvarcone20rel = getIsolation_ptvarcone20(el)/el->pt();
241  }
242  z0 = getTrack_z0(el);
243  d0 = getTrack_d0(el);
244  }
245  else et=eg->caloCluster()->et()/Gaudi::Units::GeV;
246 
247  float eta = eg->caloCluster()->etaBE(2);
248  float phi = eg->phi();
249  float pt = eg->pt()/Gaudi::Units::GeV;
250  float avgmu=lbAverageInteractionsPerCrossing( Gaudi::Hive::currentContext() );
251  float npvtx=0.0;
252 
253  ATH_MSG_DEBUG("PID decision efficiency " << pidwordDec(*eg));
254 
255  if(pid){
256  et_vec.push_back( et );
257  pt_vec.push_back( pt );
258  highet_vec.push_back( et );
259 
260  if(et > etthr+1.0){
261  eta_vec.push_back(eta);
262  phi_vec.push_back(phi);
263  avgmu_vec.push_back(avgmu);
264  npvtx_vec.push_back(npvtx);
265  ptvarcone20rel_vec.push_back(ptvarcone20rel);
266  z0_vec.push_back(z0);
267  d0_vec.push_back(d0);
268  }
269 
270  if(abs(eta)<=0.8){
271  et_slice0_vec.push_back(et);
272  }else if( abs(eta) > 0.80 && abs(eta) <= 1.37 ){
273  et_slice1_vec.push_back(et);
274  }else if( abs(eta) > 1.37 && abs(eta) <= 1.54 ){
275  et_slice2_vec.push_back(et);
276  }else if( abs(eta) > 1.54 && abs(eta) <= 2.50 ){
277  et_slice3_vec.push_back(et);
278  }
279 
280  if(isPassed) {
281  match_et_vec.push_back( et );
282  match_pt_vec.push_back( pt );
283  match_highet_vec.push_back( et );
284 
285  if(et > etthr+1.0){
286  match_eta_vec.push_back(eta);
287  match_phi_vec.push_back(phi);
288  match_avgmu_vec.push_back(avgmu);
289  match_npvtx_vec.push_back(npvtx);
290  match_ptvarcone20rel_vec.push_back(ptvarcone20rel);
291  match_z0_vec.push_back(z0);
292  match_d0_vec.push_back(d0);
293  }
294 
295  et_passed_vec.push_back( true );
296  et_failed_vec.push_back( false );
297  pt_passed_vec.push_back( true );
298  highet_passed_vec.push_back( true );
299  highet_failed_vec.push_back( false );
300 
301  if(abs(eta)<=0.8){
302  et_slice0_passed_vec.push_back(true);
303  }else if( abs(eta) > 0.80 && abs(eta) <= 1.37 ){
304  et_slice1_passed_vec.push_back(true);
305  }else if( abs(eta) > 1.37 && abs(eta) <= 1.54 ){
306  et_slice2_passed_vec.push_back(true);
307  }else if( abs(eta) > 1.54 && abs(eta) <= 2.50 ){
308  et_slice3_passed_vec.push_back(true);
309  }
310 
311  if(et > etthr+1.0){
312  eta_passed_vec.push_back( true );
313  eta_failed_vec.push_back( false );
314  phi_passed_vec.push_back( true );
315  avgmu_passed_vec.push_back( true );
316  npvtx_passed_vec.push_back( true );
317  ptvarcone20rel_passed_vec.push_back( true );
318  z0_passed_vec.push_back( true );
319  d0_passed_vec.push_back( true );
320  }
321  } // Passes Trigger selection
322  else {
323 
324  et_passed_vec.push_back( false );
325  et_failed_vec.push_back( true );
326  pt_passed_vec.push_back( false );
327  highet_passed_vec.push_back( false );
328  highet_failed_vec.push_back( true );
329 
330  if(abs(eta)<=0.8){
331  et_slice0_passed_vec.push_back(false);
332  }else if( abs(eta) > 0.80 && abs(eta) <= 1.37 ){
333  et_slice1_passed_vec.push_back(false);
334  }else if( abs(eta) > 1.37 && abs(eta) <= 1.54 ){
335  et_slice2_passed_vec.push_back(false);
336  }else if( abs(eta) > 1.54 && abs(eta) <= 2.50 ){
337  et_slice3_passed_vec.push_back(false);
338  }
339 
340  if(et > etthr+1.0){
341  eta_passed_vec.push_back( false );
342  eta_failed_vec.push_back( true );
343  phi_passed_vec.push_back( false );
344  avgmu_passed_vec.push_back( false );
345  npvtx_passed_vec.push_back( false );
346  ptvarcone20rel_passed_vec.push_back( false );
347  z0_passed_vec.push_back( false );
348  d0_passed_vec.push_back( false );
349  }
350  } // Fails Trigger selection
351 
352  } // Passes offline pid, fill histograms
353  iObj++;
354  }
355 
356  fill( monGroup, et_col, highet_col, pt_col, eta_col, phi_col, avgmu_col, npvtx_col, ptvarcone20rel_col, z0_col, d0_col,
357  match_et_col, match_highet_col, match_pt_col, match_eta_col, match_phi_col, match_avgmu_col, match_npvtx_col, match_ptvarcone20rel_col,match_z0_col,match_d0_col,
358  et_passed_col, et_failed_col, highet_passed_col, highet_failed_col, pt_passed_col, eta_passed_col, eta_failed_col, phi_passed_col, avgmu_passed_col, npvtx_passed_col, ptvarcone20rel_passed_col, z0_passed_col, d0_passed_col,
359  et_slice0_col,et_slice1_col,et_slice2_col,et_slice3_col,et_slice0_passed_col,et_slice1_passed_col,et_slice2_passed_col,et_slice3_passed_col);
360 
361 }

◆ fillHistograms()

StatusCode TrigEgammaMonitorElectronAlgorithm::fillHistograms ( const EventContext &  ctx) const
overridevirtual

adds event to the monitoring histograms

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

Parameters
ctxforwarded from execute
Returns
StatusCode

Reimplemented from TrigEgammaMonitorBaseAlgorithm.

Definition at line 46 of file TrigEgammaMonitorElectronAlgorithm.cxx.

47 {
48 
49  ATH_MSG_DEBUG("Executing TrigEgammaMonitorElectronAlgorithm");
50 
51  if(isHLTTruncated()){
52  ATH_MSG_DEBUG("HLTResult truncated, skip trigger analysis");
53  return StatusCode::SUCCESS;
54  }
55 
56  // Noise burst protection
57  SG::ReadHandle<xAOD::EventInfo> thisEvent(GetEventInfo(ctx));
58  ATH_CHECK(thisEvent.isValid());
59  if ( thisEvent->isEventFlagBitSet(xAOD::EventInfo::LAr,LArEventBitInfo::NOISEBURSTVETO)) {
60  ATH_MSG_DEBUG("LAr Noise Burst Veto, skip trigger analysis");
61  return StatusCode::SUCCESS;
62  }
63 
64  ATH_MSG_DEBUG("Chains for Analysis " << m_trigList);
65 
66  for(const auto& trigger : m_trigList){
67 
68  const TrigInfo info = getTrigInfo(trigger);
69 
70  ATH_MSG_DEBUG("Start Chain Analysis ============================= " << trigger << " " << info.trigger);
71 
72  std::vector< std::pair<std::shared_ptr<const xAOD::Egamma>, const TrigCompositeUtils::Decision*>> pairObjs;
73  if ( executeNavigation( ctx, info.trigger,info.etthr,info.pidname, pairObjs).isFailure() )
74  {
75  ATH_MSG_DEBUG("executeNavigation Fails");
76  return StatusCode::SUCCESS;
77  }
78 
79  std::vector< std::pair<const xAOD::Egamma*, const TrigCompositeUtils::Decision*>> pairObjsRaw;
80  pairObjsRaw.reserve(pairObjs.size());
81  for (const auto& itr : pairObjs) {
82  pairObjsRaw.emplace_back(itr.first.get(), itr.second);
83  }
84 
85 
86  fillDistributions( pairObjsRaw, info );
87  fillEfficiencies( pairObjsRaw, info );
88  fillResolutions( pairObjsRaw, info );
89 
90  ATH_MSG_DEBUG("End Chain Analysis ============================= " << trigger);
91  } // End loop over trigger list
92 
93 
94  return StatusCode::SUCCESS;
95 }

◆ fillHLTElectronResolution()

void TrigEgammaMonitorAnalysisAlgorithm::fillHLTElectronResolution ( const std::string &  trigger,
const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &  pairObjs,
const TrigInfo &  info 
) const
privateinherited

Definition at line 1028 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

1031 {
1032 
1033  auto monGroup = getGroup( trigger + "_Resolutions_HLT" );
1034 
1035  std::vector<float> res_pt_vec, res_et_vec, res_phi_vec, res_eta_vec, res_deta1_vec, res_deta2_vec, res_dphi2_vec, res_dphiresc_vec,
1036  res_z0_vec, res_d0_vec, res_d0sig_vec, res_eprobht_vec, res_npixhits_vec, res_nscthits_vec, res_Rhad_vec, res_Rhad1_vec, res_Reta_vec,
1037  res_Rphi_vec, res_weta1_vec, res_weta2_vec, res_wtots1_vec, res_f1_vec, res_f3_vec, res_eratio_vec, res_ethad_vec, res_ethad1_vec,
1038  et_vec, eta_vec, mu_vec;
1039  std::vector<float> res_ptcone20_vec, res_ptcone20_rel_vec, res_ptvarcone20_vec, res_ptvarcone20_rel_vec;
1040  std::vector<float> res_etInEta0_vec, res_etInEta1_vec, res_etInEta2_vec, res_etInEta3_vec;
1041 
1042 
1043  auto et_col = Monitored::Collection( "et" , et_vec );
1044  auto eta_col = Monitored::Collection( "eta" , eta_vec );
1045  auto mu_col = Monitored::Collection( "mu" , mu_vec );
1046 
1047  // For calo
1048  auto res_et_col = Monitored::Collection( "res_et" , res_et_vec );
1049  auto res_eta_col = Monitored::Collection( "res_eta" , res_eta_vec );
1050  auto res_phi_col = Monitored::Collection( "res_phi" , res_phi_vec );
1051  auto res_ethad_col = Monitored::Collection( "res_ethad" , res_ethad_vec );
1052  auto res_ethad1_col = Monitored::Collection( "res_ethad1" , res_ethad1_vec );
1053  auto res_Rhad_col = Monitored::Collection( "res_Rhad" , res_Rhad_vec );
1054  auto res_Rhad1_col = Monitored::Collection( "res_Rhad1" , res_Rhad1_vec );
1055  auto res_Reta_col = Monitored::Collection( "res_Reta" , res_Reta_vec );
1056  auto res_Rphi_col = Monitored::Collection( "res_Rphi" , res_Rphi_vec );
1057  auto res_weta1_col = Monitored::Collection( "res_weta1" , res_weta1_vec );
1058  auto res_weta2_col = Monitored::Collection( "res_weta2" , res_weta2_vec );
1059  auto res_wtots1_col = Monitored::Collection( "res_wtots1" , res_wtots1_vec );
1060  auto res_f1_col = Monitored::Collection( "res_f1" , res_f1_vec );
1061  auto res_f3_col = Monitored::Collection( "res_f3" , res_f3_vec );
1062  auto res_eratio_col = Monitored::Collection( "res_eratio" , res_eratio_vec );
1063 
1064  auto res_etInEta0_col = Monitored::Collection( "res_etInEta0" , res_etInEta0_vec );
1065  auto res_etInEta1_col = Monitored::Collection( "res_etInEta1" , res_etInEta1_vec );
1066  auto res_etInEta2_col = Monitored::Collection( "res_etInEta2" , res_etInEta2_vec );
1067  auto res_etInEta3_col = Monitored::Collection( "res_etInEta3" , res_etInEta3_vec );
1068 
1069  // For electron
1070  auto res_pt_col = Monitored::Collection( "res_pt" , res_pt_vec );
1071  auto res_deta1_col = Monitored::Collection( "res_deta1" , res_deta1_vec );
1072  auto res_deta2_col = Monitored::Collection( "res_deta2" , res_deta2_vec );
1073  auto res_dphi2_col = Monitored::Collection( "res_dphi2" , res_dphi2_vec );
1074  auto res_dphiresc_col = Monitored::Collection( "res_dphiresc" , res_dphiresc_vec );
1075  auto res_z0_col = Monitored::Collection( "res_z0" , res_z0_vec );
1076  auto res_d0_col = Monitored::Collection( "res_d0" , res_d0_vec );
1077  auto res_d0sig_col = Monitored::Collection( "res_d0sig" , res_d0sig_vec );
1078  auto res_eprobht_col = Monitored::Collection( "res_eprobht" , res_eprobht_vec );
1079  auto res_npixhits_col = Monitored::Collection( "res_npixhits" , res_npixhits_vec );
1080  auto res_nscthits_col = Monitored::Collection( "res_nscthits" , res_nscthits_vec );
1081  auto res_ptcone20_col = Monitored::Collection( "res_ptcone20" , res_ptcone20_vec );
1082  auto res_ptcone20_rel_col = Monitored::Collection( "res_ptcone20_rel" , res_ptcone20_rel_vec );
1083  auto res_ptvarcone20_col = Monitored::Collection( "res_ptvarcone20" , res_ptvarcone20_vec );
1084  auto res_ptvarcone20_rel_col = Monitored::Collection( "res_ptvarcone20_rel" , res_ptvarcone20_rel_vec );
1085 
1086 
1087 
1088 
1089  // Check for zero before filling
1090  ATH_MSG_DEBUG("Fill Resolution");
1091 
1092  std::string key = match()->key("Electrons_GSF");
1093  if(info.nogsf) key = match()->key("Electrons");
1094  if(info.lrt) key = match()->key("Electrons_LRT");
1095 
1096  for ( const auto & pairObj : pairObjs ){
1097 
1098  const xAOD::Electron *off = static_cast<const xAOD::Electron*>(pairObj.first);
1099  const xAOD::Electron *onl=nullptr;
1100 
1101 
1102  { // Get the closest electron object from the trigger starting with deltaR = 0.15
1103  float maxDeltaR=0.05;
1104  auto vec = tdt()->features<xAOD::ElectronContainer>(trigger,TrigDefs::Physics ,key );
1105  for(auto &featLinkInfo : vec ){
1106  if(! featLinkInfo.isValid() ) continue;
1107  const auto *feat = *(featLinkInfo.link);
1108  if(!feat) continue;
1109  float deltaR = dR( off->eta(), off->phi(), feat->eta(), feat->phi() );
1110  if( deltaR < maxDeltaR){
1111  maxDeltaR=deltaR;
1112  onl=feat;
1113  }
1114  }
1115  }
1116 
1117  if(!onl) continue;
1118 
1119  float val_off=0.;
1120  const float onl_eta=onl->eta();
1121  const float feta = abs(onl_eta);
1122  const float onl_et = getEt(onl)/Gaudi::Units::GeV;
1123  const float avgmu=lbAverageInteractionsPerCrossing( Gaudi::Hive::currentContext() );
1124  const float dummy=-999;
1125 
1126  eta_vec.push_back( onl_eta );
1127  et_vec.push_back( onl_et );
1128  mu_vec.push_back( avgmu );
1129 
1130  val_off=getTrack_pt(off);
1131  if(val_off!=0.){
1132  res_pt_vec.push_back( (getTrack_pt(off)-val_off)/val_off );
1133  }else{
1134  }
1135 
1136 
1137  val_off=getEt(off);
1138  if(val_off!=0.){
1139  res_et_vec.push_back( (getEt(onl)-val_off)/val_off );
1140  if( feta < 1.37 )
1141  res_etInEta0_vec.push_back((getEt(onl)-val_off)/val_off);
1142  else if( feta >=1.37 && feta <= 1.52 )
1143  res_etInEta1_vec.push_back((getEt(onl)-val_off)/val_off);
1144  else if( feta >= 1.55 && feta < 1.8 )
1145  res_etInEta2_vec.push_back((getEt(onl)-val_off)/val_off);
1146  else if( feta >= 1.8 && feta < 2.45 )
1147  res_etInEta3_vec.push_back((getEt(onl)-val_off)/val_off);
1148  }
1149 
1150  val_off=off->eta();
1151  if(val_off!=0.){
1152  res_eta_vec.push_back( (onl_eta-val_off)/val_off );
1153  }else{
1154  res_eta_vec.push_back( dummy );
1155  }
1156 
1157  val_off=off->phi();
1158  if(val_off!=0.){
1159  res_phi_vec.push_back( (onl->phi()-val_off)/val_off );
1160  }else{
1161  res_phi_vec.push_back(dummy );
1162  }
1163 
1164  val_off=getShowerShape_ethad(off);
1165  if(val_off!=0.) {
1166  res_ethad_vec.push_back((getShowerShape_ethad(onl)-val_off)/val_off);
1167  }else{
1168  res_ethad_vec.push_back( dummy );
1169  }
1170 
1171  val_off=getShowerShape_ethad1(off);
1172  if(val_off!=0){
1173  res_ethad1_vec.push_back((getShowerShape_ethad1(onl)-val_off)/val_off);
1174  }else{
1175  res_ethad1_vec.push_back( dummy);
1176  }
1177 
1178  val_off=getShowerShape_Rhad(off);
1179  if(val_off!=0.){
1180  res_Rhad_vec.push_back( (getShowerShape_Rhad(onl)-val_off)/val_off );
1181  }else{
1182  res_Rhad_vec.push_back( dummy );
1183  }
1184 
1185  val_off=getShowerShape_Rhad1(off);
1186  if(val_off!=0.){
1187  res_Rhad1_vec.push_back( (getShowerShape_Rhad1(onl)-val_off)/val_off );
1188  }else{
1189  res_Rhad1_vec.push_back( dummy );
1190  }
1191 
1192  val_off=getShowerShape_Reta(off);
1193  if(val_off!=0.){
1194  res_Reta_vec.push_back( (getShowerShape_Reta(onl)-val_off)/val_off );
1195  }else{
1196  res_Reta_vec.push_back( dummy );
1197  }
1198 
1199  val_off=getShowerShape_Rphi(off);
1200  if(val_off!=0.){
1201  res_Rphi_vec.push_back( (getShowerShape_Rphi(onl)-val_off)/val_off );
1202  }else{
1203  res_Rphi_vec.push_back( (getShowerShape_Rphi(onl)-val_off)/val_off );
1204  }
1205 
1206  val_off=getShowerShape_weta1(off);
1207  if(val_off!=0.){
1208  res_weta1_vec.push_back( (getShowerShape_weta1(onl)-val_off)/val_off );
1209  }else{
1210  res_weta1_vec.push_back( dummy );
1211  }
1212 
1213  val_off=getShowerShape_weta2(off);
1214  if(val_off!=0.){
1215  res_weta2_vec.push_back( (getShowerShape_weta2(onl)-val_off)/val_off );
1216  }else{
1217  res_weta2_vec.push_back( dummy );
1218  }
1219 
1220  val_off=getShowerShape_wtots1(off);
1221  if(val_off!=0.){
1222  res_wtots1_vec.push_back( (getShowerShape_wtots1(onl)-val_off)/val_off );
1223  }else{
1224  res_wtots1_vec.push_back( dummy );
1225  }
1226 
1227  val_off=getShowerShape_f1(off);
1228  if(val_off!=0.){
1229  res_f1_vec.push_back( (getShowerShape_f1(onl)-val_off)/val_off );
1230  }else{
1231  res_f1_vec.push_back(dummy );
1232  }
1233 
1234  val_off=getShowerShape_f3(off);
1235  if(val_off!=0.){
1236  res_f3_vec.push_back( (getShowerShape_f3(onl)-val_off)/val_off );
1237  }else{
1238  res_f3_vec.push_back( dummy );
1239  }
1240 
1241  val_off=getShowerShape_Eratio(off);
1242  if(val_off!=0.){
1243  res_eratio_vec.push_back( (getShowerShape_Eratio(onl)-val_off)/val_off );
1244  }else{
1245  res_eratio_vec.push_back( dummy );
1246  }
1247 
1248 
1249  //
1250  // Track variables
1251  //
1252 
1253  val_off=getTrack_pt(off);
1254  if(val_off!=0.){
1255  res_pt_vec.push_back( (getTrack_pt(onl)-val_off)/val_off );
1256  }else{
1257  res_pt_vec.push_back( dummy );
1258  }
1259 
1260  val_off=getEt(off);
1261  if(val_off!=0.) {
1262  res_et_vec.push_back( (getEt(onl)-val_off)/val_off );
1263  }else{
1264  res_et_vec.push_back( dummy );
1265  }
1266 
1267  val_off=getCaloTrackMatch_deltaEta1(off);
1268  if(val_off!=0.) {
1269  res_deta1_vec.push_back( (getCaloTrackMatch_deltaEta1(onl)-val_off)/val_off );
1270  }else{
1271  res_deta1_vec.push_back( dummy );
1272  }
1273 
1274  val_off=getCaloTrackMatch_deltaEta2(off);
1275  res_deta2_vec.push_back( (getCaloTrackMatch_deltaEta2(onl)-val_off)/val_off );
1276  val_off=getCaloTrackMatch_deltaPhi2(off);
1277  if(val_off!=0.) {
1278  res_dphi2_vec.push_back( (getCaloTrackMatch_deltaPhi2(onl)-val_off)/val_off );
1279  }else{
1280  res_dphi2_vec.push_back(dummy );
1281  }
1282 
1283  val_off=getCaloTrackMatch_deltaPhiRescaled2(off);
1284  res_dphiresc_vec.push_back( (getCaloTrackMatch_deltaPhiRescaled2(onl)-val_off)/val_off );
1285  // Resolution of Z0 of the track
1286  val_off=getTrack_z0(off);
1287  if(val_off!=0.) {
1288  res_z0_vec.push_back( getTrack_z0(onl)-val_off );
1289  }else{
1290  res_z0_vec.push_back( dummy );
1291  }
1292 
1293  // Absolute resolution for impact parameter
1294  val_off=getTrack_d0(off);
1295  if(val_off!=0.) {
1296  res_d0_vec.push_back( getTrack_d0(onl)-val_off );
1297  }else{
1298  res_d0_vec.push_back( dummy );
1299  }
1300 
1301  val_off=getD0sig(off);
1302  if(val_off!=0.) {
1303  res_d0sig_vec.push_back( getD0sig(onl)-val_off );
1304  }else{
1305  res_d0sig_vec.push_back( dummy );
1306  }
1307 
1308  // Absolute resolution on track summary ints/floats
1309  val_off=getTrackSummaryFloat_eProbabilityHT(off);
1310  res_eprobht_vec.push_back( (getTrackSummaryFloat_eProbabilityHT(onl)-val_off) );
1311  res_npixhits_vec.push_back( getTrackSummary_numberOfPixelHits(onl)-getTrackSummary_numberOfPixelHits(onl) );
1312  res_nscthits_vec.push_back( getTrackSummary_numberOfSCTHits(onl)-getTrackSummary_numberOfSCTHits(onl) );
1313 
1314 
1315 
1316 
1317  if(info.isolated){
1318 
1319  float val_off=getIsolation_ptcone20(off);
1320  if (val_off > 0.) {
1321  res_ptcone20_vec.push_back((getIsolation_ptcone20(onl)-val_off)/val_off);
1322  if (getEt(onl) > 0. && getEt(off) > 0.) {
1323  const float reliso_onl=getIsolation_ptcone20(onl)/getEt(onl);
1324  const float reliso_off=getIsolation_ptcone20(off)/getEt(off);
1325  res_ptcone20_rel_vec.push_back((reliso_onl-reliso_off)/reliso_off);
1326  }else{
1327  res_ptcone20_rel_vec.push_back(dummy);
1328  }
1329  }else{
1330  res_ptcone20_vec.push_back(dummy);
1331  res_ptcone20_rel_vec.push_back(dummy);
1332  }
1333 
1334  // ptvarcone20 isolation
1335  val_off=getIsolation_ptvarcone20(off);
1336  if (val_off > 0.) {
1337  if (getEt(onl) > 0. && getEt(off) > 0.) {
1338  res_ptvarcone20_vec.push_back((getIsolation_ptvarcone20(onl)-val_off)/val_off);
1339  const float reliso_onl=getIsolation_ptvarcone20(onl)/getEt(onl);
1340  const float reliso_off=getIsolation_ptvarcone20(off)/getEt(off);
1341  res_ptvarcone20_rel_vec.push_back((reliso_onl-reliso_off)/reliso_off);
1342  }else{
1343  res_ptvarcone20_rel_vec.push_back(dummy);
1344  }
1345  }else{
1346  res_ptvarcone20_vec.push_back(dummy);
1347  res_ptvarcone20_rel_vec.push_back(dummy);
1348  }
1349 
1350  }
1351 
1352 
1353 
1354 
1355 
1356  } // Loop over all offline objects
1357 
1358  // Fill everything
1359  fill( monGroup ,
1360  et_col ,
1361  eta_col ,
1362  mu_col ,
1363  res_pt_col ,
1364  res_et_col ,
1365  res_eta_col ,
1366  res_phi_col ,
1367  res_deta1_col ,
1368  res_deta2_col ,
1369  res_dphi2_col ,
1370  res_dphiresc_col,
1371  res_z0_col ,
1372  res_d0_col ,
1373  res_d0sig_col ,
1374  res_eprobht_col ,
1375  res_npixhits_col,
1376  res_nscthits_col,
1377  res_ethad_col ,
1378  res_ethad1_col ,
1379  res_Rhad_col ,
1380  res_Rhad1_col ,
1381  res_Reta_col ,
1382  res_Rphi_col ,
1383  res_weta1_col ,
1384  res_weta2_col ,
1385  res_wtots1_col ,
1386  res_f1_col ,
1387  res_f3_col ,
1388  res_eratio_col ,
1389  res_ptcone20_col ,
1390  res_ptcone20_rel_col ,
1391  res_ptvarcone20_col ,
1392  res_ptvarcone20_rel_col ,
1393  res_etInEta0_col,
1394  res_etInEta1_col,
1395  res_etInEta2_col,
1396  res_etInEta3_col );
1397 
1398 
1399 
1400 }

◆ fillHLTPhotonResolution()

void TrigEgammaMonitorAnalysisAlgorithm::fillHLTPhotonResolution ( const std::string &  trigger,
const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &  pairObjs,
const TrigInfo &  info 
) const
privateinherited

Definition at line 1406 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

1409 {
1410 
1411  auto monGroup = getGroup( trigger + "_Resolutions_HLT" );
1412 
1413  std::vector<float> res_phi_vec, res_eta_vec, res_Rhad_vec, res_Rhad1_vec, res_Reta_vec, res_ethad_vec, res_ethad1_vec,
1414  res_Rphi_vec, res_weta1_vec, res_weta2_vec, res_wtots1_vec, res_f1_vec, res_f3_vec, res_eratio_vec, et_vec, eta_vec, mu_vec;
1415 
1416  std::vector<float> res_et_vec, res_et_cnv_vec, res_et_uncnv_vec;
1417  std::vector<float> res_etInEta0_vec, res_etInEta1_vec, res_etInEta2_vec, res_etInEta3_vec;
1418  std::vector<float> res_cnv_etInEta0_vec, res_cnv_etInEta1_vec, res_cnv_etInEta2_vec, res_cnv_etInEta3_vec;
1419  std::vector<float> res_uncnv_etInEta0_vec, res_uncnv_etInEta1_vec, res_uncnv_etInEta2_vec, res_uncnv_etInEta3_vec;
1420 
1421 
1422 
1423  std::vector<float> res_topoetcone20_vec, res_topoetcone20_rel_vec;
1424 
1425  auto et_col = Monitored::Collection( "et" , et_vec );
1426  auto eta_col = Monitored::Collection( "eta" , eta_vec );
1427  auto mu_col = Monitored::Collection( "mu" , mu_vec );
1428 
1429  // For calo
1430  auto res_eta_col = Monitored::Collection( "res_eta" , res_eta_vec );
1431  auto res_phi_col = Monitored::Collection( "res_phi" , res_phi_vec );
1432  auto res_ethad_col = Monitored::Collection( "res_ethad" , res_ethad_vec );
1433  auto res_ethad1_col = Monitored::Collection( "res_ethad1" , res_ethad1_vec );
1434  auto res_Rhad_col = Monitored::Collection( "res_Rhad" , res_Rhad_vec );
1435  auto res_Rhad1_col = Monitored::Collection( "res_Rhad1" , res_Rhad1_vec );
1436  auto res_Reta_col = Monitored::Collection( "res_Reta" , res_Reta_vec );
1437  auto res_Rphi_col = Monitored::Collection( "res_Rphi" , res_Rphi_vec );
1438  auto res_weta1_col = Monitored::Collection( "res_weta1" , res_weta1_vec );
1439  auto res_weta2_col = Monitored::Collection( "res_weta2" , res_weta2_vec );
1440  auto res_wtots1_col = Monitored::Collection( "res_wtots1" , res_wtots1_vec );
1441  auto res_f1_col = Monitored::Collection( "res_f1" , res_f1_vec );
1442  auto res_f3_col = Monitored::Collection( "res_f3" , res_f3_vec );
1443  auto res_eratio_col = Monitored::Collection( "res_eratio" , res_eratio_vec );
1444 
1445 
1446  auto res_et_col = Monitored::Collection( "res_et" , res_et_vec );
1447  auto res_et_cnv_col = Monitored::Collection( "res_et_cnv" , res_et_cnv_vec );
1448  auto res_et_uncnv_col = Monitored::Collection( "res_et_uncnv" , res_et_uncnv_vec );
1449  auto res_etInEta0_col = Monitored::Collection( "res_etInEta0" , res_etInEta0_vec );
1450  auto res_etInEta1_col = Monitored::Collection( "res_etInEta1" , res_etInEta1_vec );
1451  auto res_etInEta2_col = Monitored::Collection( "res_etInEta2" , res_etInEta2_vec );
1452  auto res_etInEta3_col = Monitored::Collection( "res_etInEta3" , res_etInEta3_vec );
1453  auto res_cnv_etInEta0_col = Monitored::Collection( "res_cnv_etInEta0" , res_cnv_etInEta0_vec );
1454  auto res_cnv_etInEta1_col = Monitored::Collection( "res_cnv_etInEta1" , res_cnv_etInEta1_vec );
1455  auto res_cnv_etInEta2_col = Monitored::Collection( "res_cnv_etInEta2" , res_cnv_etInEta2_vec );
1456  auto res_cnv_etInEta3_col = Monitored::Collection( "res_cnv_etInEta3" , res_cnv_etInEta3_vec );
1457  auto res_uncnv_etInEta0_col = Monitored::Collection( "res_uncnv_etInEta0" , res_uncnv_etInEta0_vec);
1458  auto res_uncnv_etInEta1_col = Monitored::Collection( "res_uncnv_etInEta1" , res_uncnv_etInEta1_vec);
1459  auto res_uncnv_etInEta2_col = Monitored::Collection( "res_uncnv_etInEta2" , res_uncnv_etInEta2_vec);
1460  auto res_uncnv_etInEta3_col = Monitored::Collection( "res_uncnv_etInEta3" , res_uncnv_etInEta3_vec);
1461 
1462  // For photon
1463  auto res_topoetcone20_col = Monitored::Collection( "res_topoetcone20" , res_topoetcone20_vec );
1464  auto res_topoetcone20_rel_col = Monitored::Collection( "res_topoetcone20_rel" , res_topoetcone20_rel_vec );
1465 
1466 
1467 
1468 
1469  // Check for zero before filling
1470  ATH_MSG_DEBUG("Fill Resolution");
1471 
1472 
1473 
1474  for ( const auto & pairObj : pairObjs ){
1475 
1476  const xAOD::Photon *off = static_cast<const xAOD::Photon*>(pairObj.first);
1477  const xAOD::Photon *onl=nullptr;
1478 
1479 
1480  { // Get the closest electron object from the trigger starting with deltaR = 0.15
1481  float maxDeltaR=0.05;
1482  auto vec = tdt()->features<xAOD::PhotonContainer>(trigger,TrigDefs::Physics ,match()->key("Photons") );
1483  for(auto &featLinkInfo : vec ){
1484  if(! featLinkInfo.isValid() ) continue;
1485  const auto *feat = *(featLinkInfo.link);
1486  if(!feat) continue;
1487  float deltaR = dR( off->eta(), off->phi(), feat->eta(), feat->phi() );
1488  if( deltaR < maxDeltaR){
1489  maxDeltaR=deltaR;
1490  onl=feat;
1491  }
1492  }
1493  }
1494 
1495  // If not found, skip this off object!
1496  if(!onl) continue;
1497 
1498  float val_off=0.;
1499  const float onl_eta=onl->eta();
1500  const float feta = abs(onl_eta);
1501  const float onl_et = getCluster_et(onl)/Gaudi::Units::GeV;
1502  const float dummy=-999;
1503 
1504  const float avgmu=lbAverageInteractionsPerCrossing( Gaudi::Hive::currentContext() );
1505  et_vec.push_back( onl_et );
1506  eta_vec.push_back( onl_eta );
1507  mu_vec.push_back( avgmu );
1508 
1509 
1510  val_off=getCluster_et(off);
1511  if(val_off!=0.){
1512  res_et_vec.push_back( (getCluster_et(onl)-val_off)/val_off );
1513  if( feta < 1.37 )
1514  res_etInEta0_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1515  else if( feta >=1.37 && feta <= 1.52 )
1516  res_etInEta1_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1517  else if( feta >= 1.55 && feta < 1.8 )
1518  res_etInEta2_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1519  else if( feta >= 1.8 && feta < 2.45 )
1520  res_etInEta3_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1521 
1522  if(xAOD::EgammaHelpers::isConvertedPhoton(off)){
1523  res_et_cnv_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1524  if( feta < 1.37 )
1525  res_cnv_etInEta0_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1526  else if( feta >=1.37 && feta <= 1.52 )
1527  res_cnv_etInEta1_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1528  else if( feta >= 1.55 && feta < 1.8 )
1529  res_cnv_etInEta2_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1530  else if( feta >= 1.8 && feta < 2.45 )
1531  res_cnv_etInEta3_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1532  }else{
1533  res_et_uncnv_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1534  if( feta < 1.37 )
1535  res_uncnv_etInEta0_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1536  else if( feta >=1.37 && feta <= 1.52 )
1537  res_uncnv_etInEta1_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1538  else if( feta >= 1.55 && feta < 1.8 )
1539  res_uncnv_etInEta2_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1540  else if( feta >= 1.8 && feta < 2.45 )
1541  res_uncnv_etInEta3_vec.push_back((getCluster_et(onl)-val_off)/val_off);
1542  }
1543  }
1544 
1545  val_off=off->eta();
1546  if(val_off!=0.){
1547  res_eta_vec.push_back( (onl_eta-val_off)/val_off ) ;
1548  }else{
1549  res_eta_vec.push_back( dummy ) ;
1550  }
1551 
1552  val_off=off->phi();
1553  if(val_off!=0.){
1554  res_phi_vec.push_back( (onl->phi()-val_off)/val_off );
1555  }else{
1556  res_phi_vec.push_back(dummy );
1557  }
1558 
1559  val_off=getShowerShape_ethad(off);
1560  if(val_off!=0.) {
1561  res_ethad_vec.push_back((getShowerShape_ethad(onl)-val_off)/val_off);
1562  }else{
1563  res_ethad_vec.push_back(dummy);
1564  }
1565 
1566  val_off=getShowerShape_ethad1(off);
1567  if(val_off!=0){
1568  res_ethad1_vec.push_back((getShowerShape_ethad1(onl)-val_off)/val_off);
1569  }else{
1570  res_ethad1_vec.push_back(dummy);
1571  }
1572 
1573  val_off=getShowerShape_Rhad(off);
1574  if(val_off!=0.){
1575  res_Rhad_vec.push_back( (getShowerShape_Rhad(onl)-val_off)/val_off );
1576  }else{
1577  res_Rhad_vec.push_back(dummy );
1578  }
1579 
1580  val_off=getShowerShape_Rhad1(off);
1581  if(val_off!=0.){
1582  res_Rhad1_vec.push_back( (getShowerShape_Rhad1(onl)-val_off)/val_off );
1583  }else{
1584  res_Rhad1_vec.push_back(dummy );
1585  }
1586 
1587  val_off=getShowerShape_Reta(off);
1588  if(val_off!=0.){
1589  res_Reta_vec.push_back( (getShowerShape_Reta(onl)-val_off)/val_off );
1590  }else{
1591  res_Reta_vec.push_back( dummy );
1592  }
1593 
1594  val_off=getShowerShape_Rphi(off);
1595  if(val_off!=0.){
1596  res_Rphi_vec.push_back( (getShowerShape_Rphi(onl)-val_off)/val_off );
1597  }else{
1598  res_Rphi_vec.push_back(dummy );
1599  }
1600 
1601  val_off=getShowerShape_weta1(off);
1602  if(val_off!=0.){
1603  res_weta1_vec.push_back( (getShowerShape_weta1(onl)-val_off)/val_off );
1604  }else{
1605  res_weta1_vec.push_back( dummy );
1606  }
1607 
1608  val_off=getShowerShape_weta2(off);
1609  if(val_off!=0.){
1610  res_weta2_vec.push_back( (getShowerShape_weta2(onl)-val_off)/val_off );
1611  }else{
1612  res_weta2_vec.push_back(dummy);
1613  }
1614 
1615  val_off=getShowerShape_wtots1(off);
1616  if(val_off!=0.){
1617  res_wtots1_vec.push_back( (getShowerShape_wtots1(onl)-val_off)/val_off );
1618  }else{
1619  res_wtots1_vec.push_back( dummy );
1620  }
1621 
1622  val_off=getShowerShape_f1(off);
1623  if(val_off!=0.){
1624  res_f1_vec.push_back( (getShowerShape_f1(onl)-val_off)/val_off );
1625  }else{
1626  res_f1_vec.push_back( dummy );
1627  }
1628 
1629  val_off=getShowerShape_f3(off);
1630  if(val_off!=0.){
1631  res_f3_vec.push_back( (getShowerShape_f3(onl)-val_off)/val_off );
1632  }else{
1633  res_f3_vec.push_back( dummy );
1634  }
1635 
1636  val_off=getShowerShape_Eratio(off);
1637  if(val_off!=0.){
1638  res_eratio_vec.push_back( (getShowerShape_Eratio(onl)-val_off)/val_off);
1639  }else{
1640  res_eratio_vec.push_back( dummy);
1641  }
1642 
1643 
1644  if( info.isolated ){
1645  // topoetcone20 isolation
1646  float val_off=getIsolation_topoetcone20(off);
1647  float etonl=onl->pt();
1648  float etoff=off->pt();
1649  if (val_off > 0.) {
1650  res_topoetcone20_vec.push_back((getIsolation_topoetcone20(onl)-val_off)/val_off);
1651  if (etonl > 0. && etoff > 0.) {
1652  const float reliso_onl=getIsolation_topoetcone20(onl)/etonl;
1653  const float reliso_off=getIsolation_topoetcone20(off)/etoff;
1654  res_topoetcone20_rel_vec.push_back((reliso_onl-reliso_off)/reliso_off);
1655  }else{
1656  res_topoetcone20_rel_vec.push_back(dummy);
1657  }
1658  }else{
1659  res_topoetcone20_vec.push_back(dummy);
1660  res_topoetcone20_rel_vec.push_back(dummy);
1661  }
1662  }
1663 
1664 
1665 
1666  } // Loop over all offline objects
1667 
1668  // Fill everything
1669  fill( monGroup ,
1670  et_col ,
1671  eta_col ,
1672  mu_col ,
1673  res_et_col ,
1674  res_eta_col ,
1675  res_phi_col ,
1676  res_ethad_col ,
1677  res_ethad1_col ,
1678  res_Rhad_col ,
1679  res_Rhad1_col ,
1680  res_Reta_col ,
1681  res_Rphi_col ,
1682  res_weta1_col ,
1683  res_weta2_col ,
1684  res_wtots1_col ,
1685  res_f1_col ,
1686  res_f3_col ,
1687  res_eratio_col ,
1688  res_topoetcone20_col ,
1689  res_topoetcone20_rel_col ,
1690  res_etInEta0_col,
1691  res_etInEta1_col,
1692  res_etInEta2_col,
1693  res_etInEta3_col,
1694  res_et_uncnv_col,
1695  res_cnv_etInEta0_col,
1696  res_cnv_etInEta1_col,
1697  res_cnv_etInEta2_col,
1698  res_cnv_etInEta3_col,
1699  res_et_cnv_col,
1700  res_uncnv_etInEta0_col,
1701  res_uncnv_etInEta1_col,
1702  res_uncnv_etInEta2_col,
1703  res_uncnv_etInEta3_col
1704  );
1705 
1706 
1707 
1708 }

◆ fillInefficiency()

void TrigEgammaMonitorAnalysisAlgorithm::fillInefficiency ( const std::string &  pidword,
const TrigInfo &  info,
const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &  pairObjs,
const std::vector< asg::AcceptData > &  acceptObjs 
) const
protectedinherited

Definition at line 364 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

369 {
370  auto monGroup = getGroup(info.trigger+"_Inefficiency");
371  unsigned iObj=0;
372  SG::Decorator<bool> pidwordDec("is"+pidword);
373  for( auto pairObj : pairObjs ){
374 
375  bool pid=true;
376  bool isPassedL1Calo = acceptObjs[iObj].getCutResult( "L1Calo");
377  bool isPassedL2Calo = acceptObjs[iObj].getCutResult( "L2Calo");
378  bool isPassedL2 = acceptObjs[iObj].getCutResult( "L2" );
379  bool isPassedEFCalo = acceptObjs[iObj].getCutResult( "EFCalo");
380  bool isPassedHLT = acceptObjs[iObj].getCutResult( "HLT" );
381 
382  const auto *eg = pairObj.first;
383  if(xAOD::EgammaHelpers::isElectron(eg)){
384  ATH_MSG_DEBUG("Offline Electron with pidword " << pidword);
385  const xAOD::Electron* el =static_cast<const xAOD::Electron*> (eg);
386  pid=pidwordDec(*el);
387  ATH_MSG_DEBUG("Electron pid " << pid);
388  }
389  if(pid){
390  if(!isPassedL1Calo){
391  fillLabel(monGroup, "InefficiencyCounts", "L1Calo" );
392  }if(!isPassedL2Calo){
393  fillLabel(monGroup, "InefficiencyCounts", "L2Calo" );
394  }if(!isPassedL2){
395  fillLabel(monGroup, "InefficiencyCounts", "L2" );
396  }if(!isPassedEFCalo){
397  fillLabel(monGroup, "InefficiencyCounts", "EFCalo" );
398  }if(!isPassedHLT){
399  fillLabel(monGroup, "InefficiencyCounts", "HLT" );
400  }
401  }
402  iObj++;
403  }
404 }

◆ fillL1Calo()

void TrigEgammaMonitorAnalysisAlgorithm::fillL1Calo ( const std::string &  trigger,
const std::vector< const xAOD::EmTauRoI * > &  l1_vec 
) const
privateinherited

Definition at line 562 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

563 {
564  auto monGroup = getGroup(trigger+"_Distributions_L1Calo");
565 
566  std::vector<float> eta_vec, phi_vec, energy_vec, roi_et_vec, emIso_vec, hadCore_vec;
567 
568  auto eta_col = Monitored::Collection( "eta" , eta_vec );
569  auto phi_col = Monitored::Collection( "phi" , phi_vec );
570  auto energy_col = Monitored::Collection( "energy" , energy_vec );
571  auto roi_et_col = Monitored::Collection( "roi_et" , roi_et_vec );
572  auto emIso_col = Monitored::Collection( "emIso" , emIso_vec );
573  auto hadCore_col = Monitored::Collection( "hadCore" , hadCore_vec );
574 
575  for( const auto *l1 : l1_vec )
576  {
577  if(!l1) continue;
578  eta_vec.push_back( l1->eta() );
579  phi_vec.push_back( l1->phi() );
580  energy_vec.push_back( l1->emClus()/Gaudi::Units::GeV );
581  roi_et_vec.push_back( l1->eT()/Gaudi::Units::GeV );
582  emIso_vec.push_back( l1->emIsol()/Gaudi::Units::GeV );
583  hadCore_vec.push_back( l1->hadCore()/Gaudi::Units::GeV );
584  }
585 
586  fill( monGroup, eta_col, phi_col, energy_col, roi_et_col, emIso_col, hadCore_col );
587 
588 }

◆ fillL1CaloAbsResolution()

void TrigEgammaMonitorAnalysisAlgorithm::fillL1CaloAbsResolution ( const std::string &  trigger,
const std::vector< std::pair< const xAOD::Egamma *, const xAOD::EmTauRoI * >> &  pairObjs 
) const
privateinherited

Definition at line 1000 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

1002 {
1003  auto monGroup = getGroup( trigger + "_AbsResolutions_L1Calo" );
1004 
1005  std::vector<float> eta_vec, res_et_vec;
1006 
1007  auto eta_col = Monitored::Collection( "eta" , eta_vec );
1008  auto res_et_col = Monitored::Collection( "res_et" , res_et_vec );
1009 
1010 
1011  for (const auto & pairObj : pairObjs){
1012  const auto *off = pairObj.first;
1013  const auto *l1 = pairObj.second;
1014  ATH_MSG_DEBUG("Fill L1CaloAbsResolution");
1015  if(off->type()==xAOD::Type::Electron){
1016  const xAOD::Electron* eloff =static_cast<const xAOD::Electron*> (off);
1017  eta_vec.push_back( l1->eta() );
1018  res_et_vec.push_back( (l1->emClus()-getEt(eloff))/Gaudi::Units::GeV ) ;
1019  }
1020  }
1021 
1022  fill( monGroup, eta_col, res_et_col );
1023 }

◆ fillL1CaloResolution()

void TrigEgammaMonitorAnalysisAlgorithm::fillL1CaloResolution ( const std::string &  trigger,
const std::vector< std::pair< const xAOD::Egamma *, const xAOD::EmTauRoI * >> &  pairObjs 
) const
privateinherited

Definition at line 970 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

972 {
973  auto monGroup = getGroup( trigger + "_Resolutions_L1Calo" );
974 
975  std::vector<float> eta_vec, res_et_vec;
976 
977  auto eta_col = Monitored::Collection( "eta" , eta_vec );
978  auto res_et_col = Monitored::Collection( "res_et" , res_et_vec );
979 
980 
981  for (const auto & pairObj : pairObjs){
982  const auto *off = pairObj.first;
983  const auto *l1 = pairObj.second;
984  ATH_MSG_DEBUG("Fill L1CaloResolution");
985  if(off->type()==xAOD::Type::Electron){
986  const xAOD::Electron* eloff =static_cast<const xAOD::Electron*> (off);
987  eta_vec.push_back( l1->eta() );
988  res_et_vec.push_back( (l1->emClus()-getEt(eloff))/getEt(eloff) ) ;
989 
990  }
991  }
992 
993  fill( monGroup, eta_col, res_et_col );
994 }

◆ fillL1eEM()

void TrigEgammaMonitorAnalysisAlgorithm::fillL1eEM ( const std::string &  trigger,
const std::vector< const xAOD::eFexEMRoI * > &  l1_vec 
) const
privateinherited

Definition at line 591 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

592 {
593  auto monGroup = getGroup(trigger+"_Distributions_L1Calo");
594 
595  std::vector<float> eta_vec, phi_vec, et_vec, wstot_vec, reta_vec, rhad_vec;
596 
597  auto et_col = Monitored::Collection( "et" , et_vec );
598  auto eta_col = Monitored::Collection( "eta" , eta_vec );
599  auto phi_col = Monitored::Collection( "phi" , phi_vec );
600  auto wstot_col = Monitored::Collection( "Wstot" , wstot_vec );
601  auto reta_col = Monitored::Collection( "Reta" , reta_vec );
602  auto rhad_col = Monitored::Collection( "Rhad" , rhad_vec );
603 
604 
605  for( const auto *l1 : l1_vec )
606  {
607  if(!l1) continue;
608  et_vec.push_back( l1->et()/Gaudi::Units::GeV );
609  eta_vec.push_back( l1->eta() );
610  phi_vec.push_back( l1->phi() );
611  wstot_vec.push_back( l1->Wstot() );
612  reta_vec.push_back( l1->Reta() );
613  rhad_vec.push_back( l1->Rhad() );
614  }
615 
616  fill( monGroup, eta_col, phi_col, et_col, wstot_col, reta_col, rhad_col );
617 
618 }

◆ fillL2Calo()

void TrigEgammaMonitorAnalysisAlgorithm::fillL2Calo ( const std::string &  trigger,
const std::vector< const xAOD::TrigEMCluster * > &  emCluster_vec 
) const
privateinherited

Definition at line 622 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

623 {
624  auto monGroup = getGroup(trigger+"_Distributions_L2Calo");
625 
626  std::vector<float> et_vec,highet_vec, eta_vec, phi_vec;
627 
628  auto et_col = Monitored::Collection("et" , et_vec );
629  auto highet_col = Monitored::Collection("highet" , highet_vec );
630  auto eta_col = Monitored::Collection("eta", eta_vec );
631  auto phi_col = Monitored::Collection("phi", phi_vec );
632 
633  for ( const auto *emCluster : emCluster_vec )
634  {
635  if(!emCluster) continue;
636  et_vec.push_back( emCluster->et()/Gaudi::Units::GeV );
637  highet_vec.push_back( emCluster->et()/Gaudi::Units::GeV );
638  eta_vec.push_back( emCluster->eta() );
639  phi_vec.push_back( emCluster->phi() );
640  }
641 
642  fill( monGroup, et_col, eta_col, phi_col, highet_col );
643 
644 
645 }

◆ fillL2CaloResolution()

void TrigEgammaMonitorAnalysisAlgorithm::fillL2CaloResolution ( const std::string &  trigger,
const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &  pairObjs 
) const
privateinherited

Definition at line 1711 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

1714 {
1715  ATH_MSG_DEBUG("Fill L2Calo Resolution");
1716 
1717  auto monGroup = getGroup( trigger + "_Resolutions_L2Calo" );
1718 
1719  std::vector<float> res_et_vec, res_phi_vec, res_eta_vec, res_Rhad_vec, res_Rhad1_vec, res_Reta_vec, res_ethad_vec, res_ethad1_vec,
1720  res_Rphi_vec, res_weta2_vec, res_f1_vec, res_f3_vec, res_eratio_vec, et_vec, eta_vec;
1721 
1722 
1723  auto et_col = Monitored::Collection( "et" , et_vec );
1724  auto eta_col = Monitored::Collection( "eta" , eta_vec );
1725  auto res_et_col = Monitored::Collection( "res_et" , res_et_vec );
1726  auto res_eta_col = Monitored::Collection( "res_eta" , res_eta_vec );
1727  auto res_phi_col = Monitored::Collection( "res_phi" , res_phi_vec );
1728  auto res_ethad_col = Monitored::Collection( "res_ethad" , res_ethad_vec );
1729  auto res_ethad1_col = Monitored::Collection( "res_ethad1" , res_ethad1_vec );
1730  auto res_Rhad_col = Monitored::Collection( "res_Rhad" , res_Rhad_vec );
1731  auto res_Rhad1_col = Monitored::Collection( "res_Rhad1" , res_Rhad1_vec );
1732  auto res_Reta_col = Monitored::Collection( "res_Reta" , res_Reta_vec );
1733  auto res_Rphi_col = Monitored::Collection( "res_Rphi" , res_Rphi_vec );
1734  auto res_weta2_col = Monitored::Collection( "res_weta2" , res_weta2_vec );
1735  auto res_f1_col = Monitored::Collection( "res_f1" , res_f1_vec );
1736  auto res_f3_col = Monitored::Collection( "res_f3" , res_f3_vec );
1737  auto res_eratio_col = Monitored::Collection( "res_eratio" , res_eratio_vec );
1738 
1739 
1740  for ( const auto & pairObj : pairObjs ){
1741 
1742 
1743  const xAOD::Egamma *off = pairObj.first;
1744  const xAOD::TrigEMCluster *onl=nullptr;
1745 
1746  { // Get the closest electron object from the trigger starting with deltaR = 0.15
1747  float maxDeltaR=0.05;
1748  auto vec = tdt()->features<xAOD::TrigEMClusterContainer>(trigger,TrigDefs::Physics ,match()->key("FastCalo") );
1749  for(auto &featLinkInfo : vec ){
1750  if(! featLinkInfo.isValid() ) continue;
1751  const auto *feat = *(featLinkInfo.link);
1752  if(!feat) continue;
1753  float deltaR = dR( off->eta(), off->phi(), feat->eta(), feat->phi() );
1754  if( deltaR < maxDeltaR){
1755  maxDeltaR=deltaR;
1756  onl=feat;
1757  }
1758  }
1759  }
1760 
1761  if(!onl) continue;
1762 
1763  et_vec.push_back(onl->et()*1e-3);
1764  eta_vec.push_back(onl->eta());
1765  const float dummy=-999;
1766 
1767  float val_off=0.;
1768 
1769  val_off=off->caloCluster()->et();
1770  if(val_off!=0.){
1771  res_et_vec.push_back(((onl->et())-val_off)/val_off);
1772  }else{
1773  res_et_vec.push_back(dummy);
1774  }
1775 
1776  val_off=off->caloCluster()->eta();
1777  if(val_off!=0.){
1778  res_eta_vec.push_back((onl->eta()-val_off)/val_off);
1779  }else{
1780  res_eta_vec.push_back(dummy);
1781  }
1782 
1783  val_off=off->caloCluster()->phi();
1784  if(val_off!=0.){
1785  res_phi_vec.push_back((onl->phi()-val_off)/val_off);
1786  }else{
1787  res_phi_vec.push_back((onl->phi()-val_off)/val_off);
1788  }
1789 
1790  float elonl_ethad = onl->energy( CaloSampling::HEC0 ); elonl_ethad += onl->energy( CaloSampling::HEC1 );
1791  elonl_ethad += onl->energy( CaloSampling::HEC2 ); elonl_ethad += onl->energy( CaloSampling::HEC3 );
1792  elonl_ethad += onl->energy( CaloSampling::TileBar0 ); elonl_ethad += onl->energy( CaloSampling::TileExt0 );
1793  elonl_ethad += onl->energy( CaloSampling::TileBar1 ); elonl_ethad += onl->energy( CaloSampling::TileExt1 );
1794  elonl_ethad += onl->energy( CaloSampling::TileBar2 ); elonl_ethad += onl->energy( CaloSampling::TileExt2 );
1795  elonl_ethad /= TMath::CosH(onl->eta() );
1796  val_off=getShowerShape_ethad(off);
1797  if(val_off!=0.){
1798  res_ethad_vec.push_back((elonl_ethad-val_off)/val_off);
1799  }else{
1800  res_ethad_vec.push_back(dummy);
1801  }
1802 
1803  val_off=getShowerShape_ethad1(off);
1804  if(val_off!=0.){
1805  res_ethad1_vec.push_back(( (onl->ehad1()/TMath::Abs(onl->eta()) )-val_off)/val_off);
1806  }else{
1807  res_ethad1_vec.push_back(dummy);
1808  }
1809 
1810  float elonl_Rhad = elonl_ethad / onl->energy() ;
1811  val_off=getShowerShape_Rhad(off);
1812  if(val_off!=0.){
1813  res_Rhad_vec.push_back(( elonl_Rhad-val_off)/val_off);
1814  }else{
1815  res_Rhad_vec.push_back(dummy);
1816  }
1817 
1818  float elonl_Rhad1 = onl->ehad1() / onl->energy() ;
1819  val_off=getShowerShape_Rhad1(off);
1820  if(val_off!=0.){
1821  res_Rhad1_vec.push_back(( elonl_Rhad1-val_off)/val_off);
1822  }else{
1823  res_Rhad1_vec.push_back(dummy);
1824  }
1825 
1826  float onl_reta= 999.0;
1827  if ( fabsf ( onl->e277() ) > 0.01 ) onl_reta = onl->e237() / onl->e277();
1828  val_off=getShowerShape_Reta(off);
1829  if(val_off!=0.){
1830  res_Reta_vec.push_back( (onl_reta -val_off)/val_off);
1831  }else{
1832  res_Reta_vec.push_back(dummy);
1833  }
1834 
1835  val_off=getShowerShape_weta2(off);
1836  if(val_off!=0.){
1837  res_weta2_vec.push_back(( (onl->weta2())-val_off)/val_off);
1838  }else{
1839  res_weta2_vec.push_back(dummy);
1840  }
1841 
1842  float onl_f1 = onl->energy(CaloSampling::EMB1)+onl->energy(CaloSampling::EME1);
1843  onl_f1 /= onl->energy();
1844  val_off=getShowerShape_f1(off);
1845  if(val_off!=0.){
1846  res_f1_vec.push_back(( (onl_f1)-val_off)/val_off);
1847  }else{
1848  res_f1_vec.push_back(dummy);
1849  }
1850 
1851  float onl_f3 = onl->energy(CaloSampling::EMB3)+onl->energy(CaloSampling::EME3);
1852  onl_f3 /= onl->energy();
1853  val_off=getShowerShape_f3(off);
1854  if(val_off!=0.){
1855  res_f3_vec.push_back(( (onl_f3)-val_off)/val_off);
1856  }else{
1857  res_f3_vec.push_back(dummy);
1858  }
1859 
1860  float onl_eratio = 999.0;
1861  if ( fabsf(onl->emaxs1() + onl->e2tsts1()) > 0.01 )
1862  onl_eratio = (onl->emaxs1() - onl->e2tsts1()) / (onl->emaxs1() + onl->e2tsts1());
1863  val_off=getShowerShape_Eratio(off);
1864  if(val_off!=0.){
1865  res_eratio_vec.push_back(( (onl_eratio)-val_off)/val_off);
1866  }else{
1867  res_eratio_vec.push_back(dummy);
1868  }
1869 
1870  }// Loop over all pair objects
1871 
1872 
1873  // Fill everything
1874  fill( monGroup ,
1875  et_col ,
1876  eta_col ,
1877  res_et_col ,
1878  res_eta_col ,
1879  res_phi_col ,
1880  res_ethad_col ,
1881  res_ethad1_col ,
1882  res_Rhad_col ,
1883  res_Rhad1_col ,
1884  res_Reta_col ,
1885  res_Rphi_col ,
1886  res_weta2_col ,
1887  res_f1_col ,
1888  res_f3_col ,
1889  res_eratio_col
1890  );
1891 
1892 }

◆ fillL2Electron()

void TrigEgammaMonitorAnalysisAlgorithm::fillL2Electron ( const std::string &  trigger,
const std::vector< const xAOD::TrigElectron * > &  el_vec 
) const
privateinherited

Definition at line 649 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

650 {
651 
652  auto monGroup = getGroup(trigger+"_Distributions_L2Electron");
653 
654  std::vector<float> et_vec, eta_vec, phi_vec, highet_vec;
655 
656  auto et_col = Monitored::Collection("et" , et_vec );
657  auto highet_col = Monitored::Collection("highet" , highet_vec );
658  auto eta_col = Monitored::Collection("eta", eta_vec );
659  auto phi_col = Monitored::Collection("phi", phi_vec );
660 
661  for ( const auto *el : el_vec )
662  {
663  if(!el) continue;
664  et_vec.push_back( el->pt()/Gaudi::Units::GeV );
665  highet_vec.push_back( el->pt()/Gaudi::Units::GeV );
666  eta_vec.push_back( el->eta() );
667  phi_vec.push_back( el->phi() );
668  }
669 
670  fill( monGroup, et_col, eta_col, phi_col, highet_col );
671 }

◆ fillL2Photon()

void TrigEgammaMonitorAnalysisAlgorithm::fillL2Photon ( const std::string &  trigger,
const std::vector< const xAOD::TrigPhoton * > &  eg_vec 
) const
privateinherited

Definition at line 673 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

674 {
675 
676  auto monGroup = getGroup(trigger+"_Distributions_L2Photon");
677 
678  std::vector<float> et_vec, eta_vec, phi_vec, highet_vec;
679 
680  auto et_col = Monitored::Collection("et" , et_vec );
681  auto highet_col = Monitored::Collection("highet" , highet_vec );
682  auto eta_col = Monitored::Collection("eta", eta_vec );
683  auto phi_col = Monitored::Collection("phi", phi_vec );
684 
685  for ( const auto *ph : ph_vec )
686  {
687  if(!ph) continue;
688  et_vec.push_back( ph->pt()/Gaudi::Units::GeV );
689  highet_vec.push_back( ph->pt()/Gaudi::Units::GeV );
690  eta_vec.push_back( ph->eta() );
691  phi_vec.push_back( ph->phi() );
692  }
693 
694  fill( monGroup, et_col, eta_col, phi_col, highet_col );
695 }

◆ fillLabel()

void TrigEgammaMonitorAnalysisAlgorithm::fillLabel ( const ToolHandle< GenericMonitoringTool > &  groupHandle,
const std::string &  histname,
const std::string &  label 
) const
protectedinherited

Definition at line 30 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

33 {
34  auto mon = Monitored::Scalar<std::string>( histname, label );
35  fill( groupHandle, mon );
36 }

◆ fillResolutions()

void TrigEgammaMonitorAnalysisAlgorithm::fillResolutions ( const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &  pairObjs,
const TrigInfo &  info 
) const
protectedinherited

***********************************************************************************************************************8

Definition at line 914 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

916 {
917 
918  std::vector< std::pair< const xAOD::Egamma*, const xAOD::EmTauRoI * >> pair_l1_vec;
919  std::vector< std::pair< const xAOD::Egamma*, const TrigCompositeUtils::Decision * >> pair_eg_vec;
920  const std::string trigger = info.trigger;
921 
922  SG::Decorator<bool> pidnameDec("is"+info.pidname);
923  for( auto pairObj : pairObjs ){
924 
925  const xAOD::Egamma* eg = pairObj.first;
926  const auto *feat = pairObj.second;
927 
928  if (feat){
929 
930  //
931  // Get only off and l1 where the offline object passed by the offline pid selector
932  //
933  const auto *l1 = match()->getL1Feature( feat );
934  if(eg->type()==xAOD::Type::Electron){
935  const xAOD::Electron* el = static_cast<const xAOD::Electron*>(eg);
936  float et = getEt(el)/Gaudi::Units::GeV;
937  if( et < info.etthr-5.0) continue;
938  if(!pidnameDec(*eg)) continue;
939  pair_eg_vec.emplace_back(el,feat);
940  if(l1) pair_l1_vec.emplace_back(eg,l1 );
941  }
942  else if(eg->type()==xAOD::Type::Photon){
943  float et = getCluster_et(eg)/Gaudi::Units::GeV;
944  if( et < info.etthr-5.0) continue;
945  pair_eg_vec.emplace_back(eg,feat);
946  if(l1) pair_l1_vec.emplace_back(eg,l1 );
947  }
948  }
949 
950  }
951 
952  // Fill L1Calo for all level 1 objects found
953  fillL1CaloResolution( trigger, pair_l1_vec );
954  fillL1CaloAbsResolution( trigger, pair_l1_vec );
955  fillL2CaloResolution( trigger, pair_eg_vec );
956 
957  // Fill HLT electron for all onl objects found
958  if ( info.signature=="Electron"){
959  fillHLTElectronResolution( trigger, pair_eg_vec, info );
960  }
961  else if ( info.signature=="Photon"){
962  fillHLTPhotonResolution( trigger, pair_eg_vec, info );
963  }
964 
965 }

◆ fillShowerShapes()

void TrigEgammaMonitorAnalysisAlgorithm::fillShowerShapes ( const std::string &  trigger,
const std::vector< const xAOD::Egamma * > &  eg_vec,
bool  online 
) const
privateinherited

Definition at line 753 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

754 {
755 
756  ATH_MSG_DEBUG("Fill SS distributions: " << trigger);
757  auto monGroup = getGroup( trigger + ( online ? "_Distributions_HLT" : "_Distributions_Offline") );
758 
759  std::vector<float> Rhad_vec, Rhad1_vec, Reta_vec, Rphi_vec, weta1_vec, weta2_vec,
760  f1_vec, f3_vec, eratio_vec, et_vec, highet_vec , eta_vec, phi_vec, topoetcone20_vec, topoetcone40_shift_vec,
761  topoetcone20_rel_vec, topoetcone40_shift_rel_vec;
762 
763 
764  auto Rhad_col = Monitored::Collection("Rhad" , Rhad_vec );
765  auto Rhad1_col = Monitored::Collection("Rhad1" , Rhad1_vec );
766  auto Reta_col = Monitored::Collection("Reta" , Reta_vec );
767  auto Rphi_col = Monitored::Collection("Rphi" , Rphi_vec );
768  auto weta1_col = Monitored::Collection("weta1" , weta1_vec );
769  auto weta2_col = Monitored::Collection("weta2" , weta2_vec );
770  auto f1_col = Monitored::Collection("f1" , f1_vec );
771  auto f3_col = Monitored::Collection("f3" , f3_vec );
772  auto eratio_col = Monitored::Collection("eratio" , eratio_vec );
773  auto et_col = Monitored::Collection("et" , et_vec );
774  auto highet_col = Monitored::Collection("highet" , highet_vec );
775  auto eta_col = Monitored::Collection("eta" , eta_vec );
776  auto phi_col = Monitored::Collection("phi" , phi_vec );
777  auto topoetcone20_col = Monitored::Collection("topoetcone20", topoetcone20_vec);
778  auto topoetcone40_shift_col = Monitored::Collection("topoetcone40_shift", topoetcone40_shift_vec );
779  auto topoetcone20_rel_col = Monitored::Collection("topoetcone20_rel", topoetcone20_rel_vec);
780  auto topoetcone40_shift_rel_col = Monitored::Collection("topoetcone40_shift_rel", topoetcone40_shift_rel_vec );
781 
782  for ( const auto *eg : eg_vec ){
783 
784  if(!eg) continue;
785 
786  Rhad_vec.push_back( getShowerShape_Rhad(eg));
787  Rhad1_vec.push_back( getShowerShape_Rhad(eg));
788  Reta_vec.push_back( getShowerShape_Reta(eg));
789  Rphi_vec.push_back( getShowerShape_Rphi(eg));
790  weta1_vec.push_back( getShowerShape_weta1(eg));
791  weta2_vec.push_back( getShowerShape_weta2(eg));
792  f1_vec.push_back( getShowerShape_f1(eg));
793  f3_vec.push_back( getShowerShape_f3(eg));
794  eratio_vec.push_back( getShowerShape_Eratio(eg));
795  et_vec.push_back( eg->pt()/Gaudi::Units::GeV);
796  highet_vec.push_back( eg->pt()/Gaudi::Units::GeV);
797  eta_vec.push_back( eg->eta());
798  phi_vec.push_back( eg->phi());
799  topoetcone20_vec.push_back( getIsolation_topoetcone20(eg)/Gaudi::Units::GeV);
800  topoetcone40_shift_vec.push_back( (getIsolation_topoetcone40(eg)-2450)/Gaudi::Units::GeV );
801 
802  if (eg->pt() > 0) {
803  topoetcone20_rel_vec.push_back( getIsolation_topoetcone20(eg)/eg->pt());
804  topoetcone40_shift_rel_vec.push_back( (getIsolation_topoetcone40(eg)-2450)/eg->pt() );
805  }
806 
807  }// Loop over egamma objects
808 
809  fill( monGroup, Rhad_col, Rhad1_col, Reta_col, Rphi_col, weta1_col, weta2_col,
810  f1_col, f3_col, eratio_col, et_col, highet_col , eta_col, phi_col, topoetcone20_col, topoetcone40_shift_col,
811  topoetcone20_rel_col, topoetcone40_shift_rel_col );
812 
813 }

◆ fillTracking()

void TrigEgammaMonitorAnalysisAlgorithm::fillTracking ( const std::string &  trigger,
const std::vector< const xAOD::Electron * > &  eg_vec,
bool  online 
) const
privateinherited

Definition at line 816 of file TrigEgammaMonitorAnalysisAlgorithm.cxx.

817 {
818 
819  ATH_MSG_DEBUG("Fill tracking");
820 
821  auto monGroup = getGroup( trigger + ( online ? "_Distributions_HLT" : "_Distributions_Offline") );
822 
823  std::vector<float> deta1_vec, deta1_EMECA_vec, deta1_EMECC_vec, deta1_EMEBA_vec, deta1_EMEBC_vec, deta2_vec, dphi2_vec,
824  dphiresc_vec, eprobht_vec, npixhits_vec, nscthits_vec, charge_vec, ptcone20_vec, ptvarcone20_vec, z0_vec, d0_vec, d0sig_vec,
825  pt_vec,pt_trk_vec, ptcone20_rel_vec, ptvarcone20_rel_vec, eta_vec, mu_vec;
826 
827  auto deta1_col = Monitored::Collection( "deta1" , deta1_vec );
828  auto deta1_EMECA_col = Monitored::Collection( "deta1_EMECA" , deta1_EMECA_vec );
829  auto deta1_EMECC_col = Monitored::Collection( "deta1_EMECC" , deta1_EMECC_vec );
830  auto deta1_EMEBA_col = Monitored::Collection( "deta1_EMEBA" , deta1_EMEBA_vec );
831  auto deta1_EMEBC_col = Monitored::Collection( "deta1_EMEBC" , deta1_EMEBC_vec );
832  auto deta2_col = Monitored::Collection( "deta2" , deta2_vec );
833  auto dphi2_col = Monitored::Collection( "dphi2" , dphi2_vec );
834  auto dphiresc_col = Monitored::Collection( "dphiresc" , dphiresc_vec );
835  auto eprobht_col = Monitored::Collection( "eprobht" , eprobht_vec );
836  auto npixhits_col = Monitored::Collection( "npixhits" , npixhits_vec );
837  auto nscthits_col = Monitored::Collection( "nscthits" , nscthits_vec );
838  auto charge_col = Monitored::Collection( "charge" , charge_vec );
839  auto ptcone20_col = Monitored::Collection( "ptcone20" , ptcone20_vec );
840  auto ptvarcone20_col = Monitored::Collection( "ptvarcone20" , ptvarcone20_vec );
841  auto z0_col = Monitored::Collection( "z0" , z0_vec );
842  auto d0_col = Monitored::Collection( "d0" , d0_vec );
843  auto d0sig_col = Monitored::Collection( "d0sig" , d0sig_vec );
844  auto pt_col = Monitored::Collection( "pt" , pt_vec );
845  auto pt_trk_col = Monitored::Collection( "pt_track" , pt_trk_vec );
846  auto ptcone20_rel_col = Monitored::Collection( "ptcone20_rel", ptcone20_rel_vec );
847  auto ptvarcone20_rel_col = Monitored::Collection( "ptvarcone20_rel" , ptvarcone20_rel_vec );
848 
849  auto eta_col = Monitored::Collection( "eta" , eta_vec );
850  auto mu_col = Monitored::Collection( "mu" , mu_vec );
851 
852 
853  for ( const auto *eg : eg_vec ){
854 
855  if(!eg) continue;
856 
857  float cleta = 0.;
858  if(eg->caloCluster()) cleta=eg->caloCluster()->eta();
859  else cleta=eg->eta();
860 
861  deta1_vec.push_back( getCaloTrackMatch_deltaEta1(eg));
862 
863  if(cleta > 1.375 && cleta < 3.2)
864  deta1_EMECA_vec.push_back( getCaloTrackMatch_deltaEta1(eg));
865  if(cleta < -1.375 && cleta > -3.2)
866  deta1_EMECC_vec.push_back( getCaloTrackMatch_deltaEta1(eg));
867  if(cleta > 0 && cleta < 1.375)
868  deta1_EMEBA_vec.push_back( getCaloTrackMatch_deltaEta1(eg));
869  if(cleta < 0 && cleta > -1.375)
870  deta1_EMEBC_vec.push_back( getCaloTrackMatch_deltaEta1(eg));
871 
872  deta2_vec.push_back( getCaloTrackMatch_deltaEta2(eg));
873  dphi2_vec.push_back( getCaloTrackMatch_deltaPhi2(eg));
874  dphiresc_vec.push_back( getCaloTrackMatch_deltaPhiRescaled2(eg));
875  eprobht_vec.push_back( getTrackSummaryFloat_eProbabilityHT(eg));
876  npixhits_vec.push_back( getTrackSummary_numberOfPixelHits(eg));
877  nscthits_vec.push_back( getTrackSummary_numberOfSCTHits(eg));
878  charge_vec.push_back( eg->charge());
879  ptcone20_vec.push_back( getIsolation_ptcone20(eg)/Gaudi::Units::GeV);
880  ptvarcone20_vec.push_back( getIsolation_ptvarcone20(eg)/Gaudi::Units::GeV);
881 
882  // Quantities directly from tracks
883  ATH_MSG_DEBUG("Get track Quantities");
884  z0_vec.push_back( getTrack_z0(eg));
885  d0_vec.push_back( getTrack_d0(eg));
886  d0sig_vec.push_back(getD0sig(eg));
887  pt_vec.push_back( eg->pt()/Gaudi::Units::GeV);
888  pt_trk_vec.push_back( getTrack_pt(eg)/Gaudi::Units::GeV);
889  eta_vec.push_back(eg->eta());
890 
891  const float avgmu=lbAverageInteractionsPerCrossing( Gaudi::Hive::currentContext() );
892  mu_vec.push_back( avgmu );
893 
894  if (eg->pt() > 0) {
895  ptcone20_rel_vec.push_back( getIsolation_ptcone20(eg)/eg->pt());
896  ptvarcone20_rel_vec.push_back( getIsolation_ptvarcone20(eg)/eg->pt());
897  }
898 
899  }
900 
901 
902  fill( monGroup, deta1_col, deta1_EMECA_col, deta1_EMECC_col, deta1_EMEBA_col, deta1_EMEBC_col, deta2_col, dphi2_col,
903  dphiresc_col, eprobht_col, npixhits_col, nscthits_col, charge_col, ptcone20_col, ptvarcone20_col, z0_col, d0_col, d0sig_col,
904  pt_col, ptcone20_rel_col, ptvarcone20_rel_col, eta_col, mu_col,pt_trk_col);
905 }

◆ filterPassed()

virtual bool AthReentrantAlgorithm::filterPassed ( const EventContext &  ctx) const
inlinevirtualinherited

Definition at line 135 of file AthReentrantAlgorithm.h.

135  {
136  return execState( ctx ).filterPassed();
137  }

◆ getD0sig()

float TrigEgammaMonitorBaseAlgorithm::getD0sig ( const xAOD::Electron *  eg) const
protectedinherited

Definition at line 362 of file TrigEgammaMonitorBaseAlgorithm.cxx.

362  {
363  const xAOD::TrackParticle* t = eg->trackParticle();
364  float d0sigma=0.;
365  if (t)
366  {
367  float vard0 = t->definingParametersCovMatrix()(0,0);
368  if (vard0 > 0) {
369  d0sigma=sqrtf(vard0);
370  }
371  else return -99.;
372 
373  if (fabs(d0sigma) < 1e-6) return -99.;
374  return t->d0()/d0sigma;
375  }
376  else return -99.;
377 }

◆ getDEmaxs1()

float TrigEgammaMonitorBaseAlgorithm::getDEmaxs1 ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 312 of file TrigEgammaMonitorBaseAlgorithm.cxx.

312  {
313  if(eg){
314  float emax2=0.;
315  eg->showerShapeValue(emax2, xAOD::EgammaParameters::e2tsts1);
316  float emax=0.;
317  eg->showerShapeValue(emax, xAOD::EgammaParameters::emaxs1);
318  float den = emax+emax2;
319 
320  if (fabs(den) < 1e-6) return -99.;
321 
322  float val = (emax-emax2)/(den);
323  return val;
324  }
325  else return -99.;
326 }

◆ getE0Eaccordion()

float TrigEgammaMonitorBaseAlgorithm::getE0Eaccordion ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 428 of file TrigEgammaMonitorBaseAlgorithm.cxx.

428  {
429  if(eg && (eg->caloCluster())){
430  const xAOD::CaloCluster* cluster = eg->caloCluster();
431  float ebe0 = cluster->energyBE(0);
432  float ebe1 = cluster->energyBE(1);
433  float ebe2 = cluster->energyBE(2);
434  float ebe3 = cluster->energyBE(3);
435  float eacc = ebe1+ebe2+ebe3;
436  if(eacc==0.) return 0.;
437  return (ebe0/eacc);
438  }
439  else return 0.;
440 }

◆ getEaccordion()

float TrigEgammaMonitorBaseAlgorithm::getEaccordion ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 416 of file TrigEgammaMonitorBaseAlgorithm.cxx.

416  {
417  if(eg && (eg->caloCluster())){
418  const xAOD::CaloCluster* cluster = eg->caloCluster();
419  float ebe1 = cluster->energyBE(1);
420  float ebe2 = cluster->energyBE(2);
421  float ebe3 = cluster->energyBE(3);
422  return (ebe1+ebe2+ebe3);
423  }
424  else return 0.;
425 }

◆ getEnergyBE0()

float TrigEgammaMonitorBaseAlgorithm::getEnergyBE0 ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 380 of file TrigEgammaMonitorBaseAlgorithm.cxx.

380  {
381  if(eg && (eg->caloCluster())){
382  const xAOD::CaloCluster* cluster = eg->caloCluster();
383  return cluster->energyBE(0);
384  }
385  else return 0;
386 }

◆ getEnergyBE1()

float TrigEgammaMonitorBaseAlgorithm::getEnergyBE1 ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 389 of file TrigEgammaMonitorBaseAlgorithm.cxx.

389  {
390  if(eg && (eg->caloCluster())){
391  const xAOD::CaloCluster* cluster = eg->caloCluster();
392  return cluster->energyBE(1);
393  }
394  else return 0.;
395 }

◆ getEnergyBE2()

float TrigEgammaMonitorBaseAlgorithm::getEnergyBE2 ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 398 of file TrigEgammaMonitorBaseAlgorithm.cxx.

398  {
399  if(eg && (eg->caloCluster())){
400  const xAOD::CaloCluster* cluster = eg->caloCluster();
401  return cluster->energyBE(2);
402  }
403  else return 0.;
404 }

◆ getEnergyBE3()

float TrigEgammaMonitorBaseAlgorithm::getEnergyBE3 ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 407 of file TrigEgammaMonitorBaseAlgorithm.cxx.

407  {
408  if(eg && (eg->caloCluster())){
409  const xAOD::CaloCluster* cluster = eg->caloCluster();
410  return cluster->energyBE(3);
411  }
412  else return 0.;
413 }

◆ getEt()

float TrigEgammaMonitorBaseAlgorithm::getEt ( const xAOD::Electron *  eg) const
protectedinherited

Definition at line 291 of file TrigEgammaMonitorBaseAlgorithm.cxx.

291  {
292  if(eg && (eg->caloCluster()) && (eg->trackParticle())){
293  const xAOD::TrackParticle *trk=eg->trackParticle();
294  const xAOD::CaloCluster *clus=eg->caloCluster();
295  float eta = fabs(trk->eta());
296  return clus->e()/cosh(eta);
297  }
298  else return -99.;
299 }

◆ getEta2()

float TrigEgammaMonitorBaseAlgorithm::getEta2 ( const xAOD::Egamma *  eg) const
protectedinherited

Features helper.

Helper functions now part of base class

Definition at line 282 of file TrigEgammaMonitorBaseAlgorithm.cxx.

282  {
283  if(eg && (eg->caloCluster())){
284  const xAOD::CaloCluster* cluster = eg->caloCluster();
285  return fabs(cluster->etaBE(2));
286  }
287  else return -99.;
288 }

◆ getEtCluster37()

float TrigEgammaMonitorBaseAlgorithm::getEtCluster37 ( const xAOD::Egamma *  eg) const
protectedinherited

Definition at line 302 of file TrigEgammaMonitorBaseAlgorithm.cxx.

302  {
303  if(eg && (eg->caloCluster())){
304  const xAOD::CaloCluster* cluster = eg->caloCluster();
305  float eta2 = fabs(cluster->etaBE(2));
306  return cluster->e()/cosh(eta2);
307  }
308  else return -99.;
309 }

◆ 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 107 of file AthMonitorAlgorithm.cxx.

107  {
108  return SG::ReadHandle<xAOD::EventInfo>(m_EventInfoKey, ctx);
109 }

◆ 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 164 of file AthMonitorAlgorithm.cxx.

164  {
165  // get the pointer to the tool, and check that it exists
166  auto idx = m_toolLookupMap.find(name);
167  if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
168  return m_tools[idx->second];
169  }
170  else {
171  if (!isInitialized()) {
172  ATH_MSG_FATAL(
173  "It seems that the AthMonitorAlgorithm::initialize was not called "
174  "in derived class initialize method");
175  } else {
176  std::string available = std::accumulate(
177  m_toolLookupMap.begin(), m_toolLookupMap.end(), std::string(""),
178  [](const std::string& s, auto h) { return s + "," + h.first; });
179  ATH_MSG_FATAL("The tool " << name << " could not be found in the tool array of the "
180  << "monitoring algorithm " << m_name << ". This probably reflects a discrepancy between "
181  << "your python configuration and c++ filling code. Note: your available groups are {"
182  << available << "}.");
183  }
184  }
185  return m_dummy;
186 }

◆ getL1Item()

std::string TrigEgammaMonitorBaseAlgorithm::getL1Item ( const std::string &  trigger) const
protectedinherited

Creates static map to return L1 item from trigger name.

Definition at line 724 of file TrigEgammaMonitorBaseAlgorithm.cxx.

724  {
725  std::vector<std::string> parts;
726  boost::split(parts,trigger,boost::is_any_of("_"));
727  // L1EMXX
728  std::string l1seed = parts.back();
729  return l1seed;
730 }

◆ getSigmaD0()

float TrigEgammaMonitorBaseAlgorithm::getSigmaD0 ( const xAOD::Electron *  eg) const
protectedinherited

Definition at line 346 of file TrigEgammaMonitorBaseAlgorithm.cxx.

346  {
347  const xAOD::TrackParticle* t = eg->trackParticle();
348  float d0sigma=0.;
349  if (t)
350  {
351  float vard0 = t->definingParametersCovMatrix()(0,0);
352  if (vard0 > 0) {
353  d0sigma=sqrtf(vard0);
354  }
355  else return -99.;
356  return d0sigma;
357  }
358  else return -99.;
359 }

◆ 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 189 of file AthMonitorAlgorithm.cxx.

189  {
190  return m_trigDecTool;
191 }

◆ getTrigInfo()

TrigInfo TrigEgammaMonitorBaseAlgorithm::getTrigInfo ( const std::string &  trigger) const
protectedinherited

Get the trigger info parsed from the chain name (only single lepton triggers)

Definition at line 583 of file TrigEgammaMonitorBaseAlgorithm.cxx.

583  {
584  return m_trigInfo.at(trigger);
585 }

◆ getTrigInfoMap()

const std::map<std::string,TrigInfo>& TrigEgammaMonitorBaseAlgorithm::getTrigInfoMap ( )
inlineprotectedinherited

Helper methods.

Get the trig info map

Definition at line 127 of file TrigEgammaMonitorBaseAlgorithm.h.

127 { return m_trigInfo; }

◆ initialize()

StatusCode TrigEgammaMonitorElectronAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from TrigEgammaMonitorAnalysisAlgorithm.

Definition at line 22 of file TrigEgammaMonitorElectronAlgorithm.cxx.

23 {
24  ATH_CHECK(TrigEgammaMonitorAnalysisAlgorithm::initialize());
25 
26  ATH_CHECK(m_offElectronKey.initialize());
27  ATH_CHECK(m_electronIsolationKeyPtCone20.initialize());
28  ATH_CHECK( m_eventInfoDecorKey.initialize() );
29 
30  for(auto& trigName : m_trigInputList)
31  {
32  if(getTrigInfoMap().count(trigName) != 0){
33  ATH_MSG_DEBUG("Trigger already booked, removing from trigger list " << trigName);
34  }else {
35  m_trigList.push_back(trigName);
36  setTrigInfo(trigName);
37  }
38  }
39 
40 
41  return StatusCode::SUCCESS;
42 }

◆ inputHandles()

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

Return this algorithm's input handles.

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

◆ isClonable()

bool AthReentrantAlgorithm::isClonable ( ) const
overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

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

Definition at line 44 of file AthReentrantAlgorithm.cxx.

45 {
46  // Reentrant algorithms are clonable.
47  return true;
48 }

◆ isHLTTruncated()

bool TrigEgammaMonitorBaseAlgorithm::isHLTTruncated ( ) const
protectedinherited

Definition at line 733 of file TrigEgammaMonitorBaseAlgorithm.cxx.

733  {
734  return m_trigdec->ExperimentalAndExpertMethods().isHLTTruncated();
735 }

◆ isIsolated()

bool TrigEgammaMonitorBaseAlgorithm::isIsolated ( const xAOD::Electron *  eg,
const std::string &  isolation 
) const
protectedinherited

Check if electron fulfils isolation criteria.

Definition at line 112 of file TrigEgammaMonitorBaseAlgorithm.cxx.

112  {
113  ATH_MSG_DEBUG("Apply Isolation " << isolation);
114  float ptcone20=0;
115  bool isoStat=eg->isolationValue(ptcone20, xAOD::Iso::ptcone20);
116  if (!isoStat) {
117  ATH_MSG_DEBUG("Electron doesn't provide isolation for ptcone20");
118  return false;
119  }
120  if (!(fabs(eg->pt()) > 0)) {
121  ATH_MSG_DEBUG("Electron pt is zero, can't calculate relative isolation");
122  return false;
123  }
124  ATH_MSG_DEBUG("ptcone20 " << ptcone20);
125  float ptcone20_rel = ptcone20/eg->pt();
126  ATH_MSG_DEBUG("Relative isolation value " << ptcone20_rel);
127  if (isolation == "loose"){
128  if (ptcone20_rel > 0.1) {
129  ATH_MSG_DEBUG("Probe failing isolation");
130  return false;
131  } else {
132  ATH_MSG_DEBUG("Probe passing isolation");
133  return true;
134  }
135  }
136  else {
137  ATH_MSG_DEBUG("No valid working point defined for " << isolation << " continue without isolation");
138  }
139  return false;
140 }

◆ isPrescaled()

bool TrigEgammaMonitorBaseAlgorithm::isPrescaled ( const std::string &  trigger) const
protectedinherited

Check if the event is prescaled.

Definition at line 144 of file TrigEgammaMonitorBaseAlgorithm.cxx.

144  {
145 
146  bool efprescale=false;
147  bool l1prescale=false;
148  bool prescale=false;
149  bool rerun=true; //assume rerun for l1
150  std::string l1item="";
151 
152  if(trigger.starts_with( "L1" ))
153  l1item=trigger;
154  if(trigger.starts_with("HLT")){
155  l1item = getL1Item(trigger);
156  const unsigned int bit=tdt()->isPassedBits(trigger);
157  efprescale=bit & TrigDefs::EF_prescaled;
158  rerun=bit&TrigDefs::EF_resurrected; //Rerun, only check for HLT
159  }
160 
161 
162  ATH_MSG_DEBUG("Checking prescale for " << trigger << " " << l1item);
163  const unsigned int l1bit=tdt()->isPassedBits(l1item);
164  bool l1_afterpre=l1bit&TrigDefs::L1_isPassedAfterPrescale;
165  bool l1_beforepre=l1bit&TrigDefs::L1_isPassedBeforePrescale;
166  l1prescale=l1_beforepre && !l1_afterpre;
167  prescale=efprescale || l1prescale;
168  ATH_MSG_DEBUG("L1 prescale " << l1item << " " << l1prescale << " before " << l1_beforepre << " after " << l1_afterpre);
169  ATH_MSG_DEBUG("EF prescale " << trigger << " " << efprescale << " Prescale " << prescale);
170  if(rerun) return false; // Rerun use the event
171  if(prescale) return true; // Prescaled, reject event
172  return false; // Not prescaled, use event
173 }

◆ match()

const ToolHandle<TrigEgammaMatchingToolMT>& TrigEgammaMonitorBaseAlgorithm::match ( ) const
inlineprotectedinherited

Get the e/g match tool.

Definition at line 135 of file TrigEgammaMonitorBaseAlgorithm.h.

135 {return m_matchTool;}

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ 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 336 of file AthMonitorAlgorithm.cxx.

336  {
337  std::string item;
338  std::stringstream ss(line);
339 
340  ATH_MSG_DEBUG( "AthMonitorAlgorithm::parseList()" );
341 
342  while ( std::getline(ss, item, ',') ) {
343  std::stringstream iss(item); // remove whitespace
344  iss >> item;
345  result.push_back(item);
346  }
347 
348  return StatusCode::SUCCESS;
349 }

◆ 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();
383  PBASE::renounce (h);
384  }

◆ 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  {
365  handlesArray.renounce();
366  }

◆ rTRT()

float TrigEgammaMonitorBaseAlgorithm::rTRT ( const xAOD::Electron *  eg) const
protectedinherited

Definition at line 329 of file TrigEgammaMonitorBaseAlgorithm.cxx.

329  {
330  if(eg && eg->trackParticle()){
331  uint8_t trtHits = 0;
332  eg->trackParticleSummaryValue(trtHits,xAOD::numberOfTRTHits);
333  uint8_t trtHTHits = 0;
334  eg->trackParticleSummaryValue(trtHTHits,xAOD::numberOfTRTHighThresholdHits);
335  if (fabs(trtHits) < 1e-6) {
336  return -99.;
337  }
338  else{
339  return ( (double)trtHTHits / (double)trtHits );
340  }
341  }
342  else return -99.;
343 }

◆ setAccept()

asg::AcceptData TrigEgammaMonitorBaseAlgorithm::setAccept ( const TrigCompositeUtils::Decision *  dec,
const TrigInfo &  info 
) const
protectedinherited

Set the accept object for all trigger levels.

Definition at line 177 of file TrigEgammaMonitorBaseAlgorithm.cxx.

177  {
178 
179  ATH_MSG_DEBUG("setAccept");
180 
181  unsigned int condition=TrigDefs::includeFailedDecisions;
182 
183  asg::AcceptData acceptData (&m_accept);
184 
185  bool passedL1Calo=false;
186  bool passedL2Calo=false;
187  bool passedEFCalo=false;
188  bool passedL2=false;
189  bool passedEFTrk=false;
190  bool passedEF=false;
191 
192  if (dec) {
193 
194  auto trigger = info.trigger;
195  // Step 1
196  passedL1Calo = match()->ancestorPassed<TrigRoiDescriptorCollection>( dec , trigger , "initialRois", condition);
197 
198  if( passedL1Calo ){ // HLT item get full decision
199  // Step 2
200  passedL2Calo = match()->ancestorPassed<xAOD::TrigEMClusterContainer>(dec, trigger, match()->key("FastCalo"), condition);
201 
202  if(passedL2Calo){
203 
204  // Step 3
205  if(info.signature == "Electron"){
206  std::string key = match()->key("FastElectrons");
207  if(info.lrt) key = match()->key("FastElectrons_LRT");
208  passedL2 = match()->ancestorPassed<xAOD::TrigElectronContainer>(dec, trigger, key, condition);
209  }else if(info.signature == "Photon"){
210  passedL2 = match()->ancestorPassed<xAOD::TrigPhotonContainer>(dec, trigger, match()->key("FastPhotons"), condition);
211  }
212 
213  if(passedL2){
214 
215  // Step 4
216  std::string key = match()->key("PrecisionCalo_Electron");
217  if(info.signature == "Photon") key = match()->key("PrecisionCalo_Photon");
218  if(info.lrt) key = match()->key("PrecisionCalo_LRT");
219  if(info.ion) key = match()->key("PrecisionCalo_HI");
220 
221  passedEFCalo = match()->ancestorPassed<xAOD::CaloClusterContainer>(dec, trigger, key, condition);
222 
223  if(passedEFCalo){
224 
225  // Step 5
226  passedEFTrk=true;// Assume true for photons
227 
228  // Step 6
229  if(info.signature == "Electron"){
230  if( info.etcut || info.idperf){// etcut or idperf
231  passedEF = true; // since we dont run the preciseElectron step
232  }else{
233  std::string key = match()->key("Electrons_GSF");
234  if(info.lrt) key = match()->key("Electrons_LRT");
235  if(info.nogsf) key = match()->key("Electrons");
236  passedEF = match()->ancestorPassed<xAOD::ElectronContainer>(dec, trigger, key, condition);
237  }
238 
239  }else if(info.signature == "Photon"){
240  if (info.etcut){
241  passedEF = true; // since we dont run the precisePhoton step
242  }else{
243  passedEF = match()->ancestorPassed<xAOD::PhotonContainer>(dec, trigger, match()->key("Photons"), condition);
244  }
245  }
246  } // EFCalo
247  }// L2
248  }// L2Calo
249  }// L2Calo
250 
251  }
252 
253  acceptData.setCutResult("L1Calo",passedL1Calo);
254  acceptData.setCutResult("L2Calo",passedL2Calo);
255  acceptData.setCutResult("L2",passedL2);
256  acceptData.setCutResult("EFCalo",passedEFCalo);
257  acceptData.setCutResult("EFTrack",passedEFTrk);
258  acceptData.setCutResult("HLT",passedEF);
259  ATH_MSG_DEBUG("Accept results:");
260  ATH_MSG_DEBUG("L1: "<< passedL1Calo);
261  ATH_MSG_DEBUG("L2Calo: " << passedL2Calo);
262  ATH_MSG_DEBUG("L2: "<< passedL2);
263  ATH_MSG_DEBUG("EFCalo: "<< passedEFCalo);
264  ATH_MSG_DEBUG("HLT: "<<passedEF);
265 
266  return acceptData;
267 }

◆ setFilterPassed()

virtual void AthReentrantAlgorithm::setFilterPassed ( bool  state,
const EventContext &  ctx 
) const
inlinevirtualinherited

Definition at line 139 of file AthReentrantAlgorithm.h.

139  {
140  execState( ctx ).setFilterPassed( state );
141  }

◆ setTrigInfo()

void TrigEgammaMonitorBaseAlgorithm::setTrigInfo ( const std::string &  trigger)
protectedinherited

Set the trigger info parsed from the chain name.

Definition at line 591 of file TrigEgammaMonitorBaseAlgorithm.cxx.

591  {
592 
593  /********************************************
594  // Trigger Information struct
595  typedef struct _triginfo
596  {
597  // L1 information
598  bool L1Legacy;
599  std::string L1Threshold; //EM22VHI
600  // HLT information
601  std::string trigger; //Trigger Name
602  std::string signature; //Electron or Photon
603  float etthr; // HLT Et threshold
604  // if trigger is etcut OR idperf, pidname should be default (usually lhloose)
605  std::string pidname; // Offline loose, medium, tight, etc...
606  // extra HLT information
607  bool idperf; // Performance chain
608  bool etcut; // Et cut only chain
609  bool nogsf; // chain without gsf reconstruction
610  bool lrt; // LRT chain
611  bool ion; // Heavy Ion chain
612  std::string isolation;
613  bool isolated;
614  } TrigInfo;
615  *******************************************/
616 
617  std::map<std::string, std::string> pidMap = { {"vloose" , "loose" },
618  {"loose" , "loose" },
619  {"medium" , "medium" },
620  {"tight" , "tight" },
621  {"loose1" , "loose" },
622  {"medium1" , "medium" },
623  {"tight1" , "tight" },
624  {"lhvloose" , "lhvloose" },
625  {"lhloose" , "lhloose" },
626  {"lhmedium" , "lhmedium" },
627  {"lhtight" , "lhtight" },
628  {"dnnloose" , "dnnloose" },
629  {"dnnmedium", "dnnmedium"},
630  {"dnntight" , "dnntight" } };
631 
632  std::vector<std::string> isoNames = {"ivarloose","ivarmedium","ivartight","icaloloose","icalomedium","icalotight"};
633 
634  bool nogsf = false;
635  bool lrt = false;
636  bool ion = false;
637  bool etcut = false;
638  bool idperf = false;
639  bool isolated = false;
640 
641  std::string isolation="";
642  bool l1legacy=true;
643 
644  std::string hltinfo=trigger;
645  std::string signature = "";
646  float threshold = 0;
647  // HLT_e/gXX_(pidname/etcut/idperf)_*_L1EMXX to e/gXX_(pidname/etcut/idperf)_*_L1EMXX
648  if(boost::contains(hltinfo,"HLT")) hltinfo.erase(0,4);
649 
650 
651  std::vector<std::string> parts;
652  boost::split(parts,hltinfo,boost::is_any_of("_"));
653  std::string pidname;
654 
655  // e/gXX_(pidname/etcut/idperf)_*_L1EMXX
656  if(boost::contains(parts.at(0),"e")) {
657  signature = "Electron";
658  pidname = m_defaultProbePidElectron;
659  }else if(boost::contains(parts.at(0),"g")) {
660  signature = "Photon";
661  pidname = m_defaultProbePidPhoton;
662  }else {
663  ATH_MSG_ERROR("Cannot set trigger type from name");
664  }
665 
666  ATH_MSG_DEBUG(parts.at(1));
667  if(parts.at(1) == "idperf"){
668  ATH_MSG_DEBUG("This is idperf");
669  idperf=true;
670  }
671  else if( parts.at(1)== "etcut"){
672  ATH_MSG_DEBUG("This is etcut");
673  etcut=true;
674  }
675  else { // remap online pidname to offline pidname
676  ATH_MSG_DEBUG("This is nominal");
677  pidname = pidMap.at(parts.at(1));
678  }
679 
680 
681  // extra information
682  nogsf = boost::contains(trigger,"nogsf");
683  lrt = boost::contains(trigger,"lrt");
684  ion = boost::contains(trigger,"ion");
685 
686  for(auto& iso : isoNames){
687  if(boost::contains(trigger, iso)){
688  isolation=iso; isolated=true; break;
689  }
690  }
691 
692  // Get the threshold
693  std::string str_thr = parts.at(0);
694  str_thr.erase( 0, 1);
695  threshold = atof(str_thr.c_str());
696 
697  // L1EMXX
698  std::string l1seed = getL1Item(trigger);
699  l1legacy = !boost::contains(l1seed, "eEM");
700 
701 
702  ATH_MSG_DEBUG("=================== Chain Parser =======================");
703  ATH_MSG_DEBUG( "trigger : " << trigger );
704  ATH_MSG_DEBUG( "threshold : " << threshold);
705  ATH_MSG_DEBUG( "Pidname : " << pidname );
706  ATH_MSG_DEBUG( "signature : " << signature);
707  ATH_MSG_DEBUG( "etcut : " << (etcut?"Yes":"No"));
708  ATH_MSG_DEBUG( "idperf : " << (idperf?"Yes":"No"));
709  ATH_MSG_DEBUG( "nogsf : " << (nogsf?"Yes":"No"));
710  ATH_MSG_DEBUG( "lrt : " << (lrt?"Yes":"No"));
711  ATH_MSG_DEBUG( "HeavyIon : " << (ion?"Yes":"No"));
712  ATH_MSG_DEBUG( "Isolation : " << isolation);
713  ATH_MSG_DEBUG( "Isolated : " << (isolated?"Yes":"No"));
714  ATH_MSG_DEBUG( "L1Seed : " << l1seed << " (Is Legacy? " << (l1legacy?"Yes":"No") << ")");
715  ATH_MSG_DEBUG("========================================================");
716 
717  TrigInfo info{l1legacy,l1seed,trigger,signature,threshold,pidname,idperf,etcut,nogsf,lrt,ion,isolation,isolated};
718  m_trigInfo[trigger] = info;
719 
720 }

◆ sysExecute()

StatusCode AthReentrantAlgorithm::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 67 of file AthReentrantAlgorithm.cxx.

68 {
69  return Gaudi::Algorithm::sysExecute (ctx);
70 }

◆ sysInitialize()

StatusCode AthReentrantAlgorithm::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 InputMakerBase, and HypoBase.

Definition at line 96 of file AthReentrantAlgorithm.cxx.

96  {
97  StatusCode sc=AthCommonDataStore<AthCommonMsg<Gaudi::Algorithm>>::sysInitialize();
98 
99  if (sc.isFailure()) {
100  return sc;
101  }
102 
103  ServiceHandle<ICondSvc> cs("CondSvc",name());
104  for (auto h : outputHandles()) {
105  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
106  // do this inside the loop so we don't create the CondSvc until needed
107  if ( cs.retrieve().isFailure() ) {
108  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
109  return StatusCode::SUCCESS;
110  }
111  if (cs->regHandle(this,*h).isFailure()) {
112  sc = StatusCode::FAILURE;
113  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
114  << " with CondSvc");
115  }
116  }
117  }
118  return sc;
119 }

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

◆ tdt()

const ToolHandle<Trig::TrigDecisionTool>& TrigEgammaMonitorBaseAlgorithm::tdt ( ) const
inlineprotectedinherited

Get the TDT

Definition at line 133 of file TrigEgammaMonitorBaseAlgorithm.h.

133 {return m_trigdec;};

◆ 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 194 of file AthMonitorAlgorithm.cxx.

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

◆ 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) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_accept

asg::AcceptInfo TrigEgammaMonitorBaseAlgorithm::m_accept
protectedinherited

Definition at line 122 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_anatype

Gaudi::Property<std::string> TrigEgammaMonitorElectronAlgorithm::m_anatype {this, "Analysis", "Electron"}
private

Directory name for each algorithm.

Definition at line 43 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_dataType

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::m_dataType
protectedinherited

Instance of the DataType_t enum.

Definition at line 351 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 353 of file AthMonitorAlgorithm.h.

◆ m_defaultLBDuration

Gaudi::Property<float> AthMonitorAlgorithm::m_defaultLBDuration {this,"DefaultLBDuration",60.}
protectedinherited

Default duration of one lumi block.

Definition at line 360 of file AthMonitorAlgorithm.h.

◆ m_defaultProbePidElectron

Gaudi::Property<std::string> TrigEgammaMonitorBaseAlgorithm::m_defaultProbePidElectron {this, "DefaultProbeSelectionElectron", "lhloose"}
protectedinherited

default probe pid for electron trigitems that don't have pid in their name

Definition at line 109 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_defaultProbePidPhoton

Gaudi::Property<std::string> TrigEgammaMonitorBaseAlgorithm::m_defaultProbePidPhoton {this, "DefaultProbeSelectionPhoton", "loose"}
protectedinherited

default probe pid for photon trigitems that don't have pid in their name

Definition at line 111 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_detailedHists

Gaudi::Property<bool> TrigEgammaMonitorBaseAlgorithm::m_detailedHists {this, "DetailedHistograms", false}
protectedinherited

Include more detailed histograms.

Definition at line 119 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_detailLevel

Gaudi::Property<int> AthMonitorAlgorithm::m_detailLevel {this,"DetailLevel",0}
protectedinherited

Sets the level of detail used in the monitoring.

Definition at line 361 of file AthMonitorAlgorithm.h.

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

Gaudi::Property<std::vector<std::string> > TrigEgammaMonitorBaseAlgorithm::m_dnnname {this, "DNNResultNames", {}, }
protectedinherited

dnn names

Definition at line 117 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_doEmulation

Gaudi::Property<bool> TrigEgammaMonitorBaseAlgorithm::m_doEmulation {this, "DoEmulation", false }
protectedinherited

Do emulation.

Definition at line 104 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_DQFilterTools

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

Array of Data Quality filter tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

◆ m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy
privateinherited

Definition at line 369 of file AthMonitorAlgorithm.h.

◆ m_electronDNNTool

ToolHandleArray<IAsgElectronLikelihoodTool> TrigEgammaMonitorBaseAlgorithm::m_electronDNNTool { this, "ElectronDNNSelectorTool", {},"DNN tools" }
protectedinherited

Offline DNN Selectors.

Definition at line 98 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_electronIsEMTool

ToolHandleArray<IAsgElectronIsEMSelector> TrigEgammaMonitorBaseAlgorithm::m_electronIsEMTool {this,"ElectronIsEMSelector",{}}
protectedinherited

Offline isEM Selectors.

Definition at line 94 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_electronIsolationKeyPtCone20

SG::ReadDecorHandleKey<xAOD::ElectronContainer> TrigEgammaMonitorElectronAlgorithm::m_electronIsolationKeyPtCone20 { this, "PtCone20Decoration", m_offElectronKey, "ptcone20", "Decoration key for the ptcone20 isolation decoration" }
private

Ensure offline electron isolation decoration is retrieved after being created.

Definition at line 60 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_electronLHTool

ToolHandleArray<IAsgElectronLikelihoodTool> TrigEgammaMonitorBaseAlgorithm::m_electronLHTool {this,"ElectronLikelihoodTool",{}}
protectedinherited

Offline LH Selectors.

Definition at line 96 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_emulatorTool

ToolHandle<Trig::TrigEgammaEmulationToolMT> TrigEgammaMonitorBaseAlgorithm::m_emulatorTool
protectedinherited

Definition at line 92 of file TrigEgammaMonitorBaseAlgorithm.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 372 of file AthMonitorAlgorithm.h.

◆ m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment
protectedinherited

Instance of the Environment_t enum.

Definition at line 350 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 352 of file AthMonitorAlgorithm.h.

◆ m_eventInfoDecorKey

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

Definition at line 58 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_EventInfoKey

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

Key for retrieving EventInfo from StoreGate.

Definition at line 362 of file AthMonitorAlgorithm.h.

◆ 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 AthReentrantAlgorithm::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 153 of file AthReentrantAlgorithm.h.

◆ m_fileKey

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

Internal Athena name for file.

Definition at line 358 of file AthMonitorAlgorithm.h.

◆ m_forceEtThr

Gaudi::Property<bool> TrigEgammaMonitorElectronAlgorithm::m_forceEtThr { this, "ForceEtThreshold", true}
private

Definition at line 51 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_forcePidSelection

Gaudi::Property<bool> TrigEgammaMonitorElectronAlgorithm::m_forcePidSelection { this, "ForcePidSelection", true}
private

force pid and crack selection into electron navigation

Definition at line 49 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_forceProbeIsolation

Gaudi::Property<bool> TrigEgammaMonitorElectronAlgorithm::m_forceProbeIsolation {this, "ForceProbeIsolation", false}
private

force probe isolation

Definition at line 47 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_forceVetoVeryLoose

Gaudi::Property<bool> TrigEgammaMonitorElectronAlgorithm::m_forceVetoVeryLoose { this, "ForceVetoVeryLoose", false}
private

Definition at line 53 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_isemname

Gaudi::Property<std::vector<std::string> > TrigEgammaMonitorBaseAlgorithm::m_isemname {this, "isEMResultNames", {} }
protectedinherited

isem names

Definition at line 113 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_lbDurationDataKey

SG::ReadCondHandleKey<LBDurationCondData> AthMonitorAlgorithm::m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
protectedinherited

Definition at line 345 of file AthMonitorAlgorithm.h.

◆ m_lhname

Gaudi::Property<std::vector<std::string> > TrigEgammaMonitorBaseAlgorithm::m_lhname {this, "LHResultNames", {} }
protectedinherited

lh names

Definition at line 115 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_lumiDataKey

SG::ReadCondHandleKey<LuminosityCondData> AthMonitorAlgorithm::m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
protectedinherited

Definition at line 343 of file AthMonitorAlgorithm.h.

◆ m_matchTool

ToolHandle<TrigEgammaMatchingToolMT> TrigEgammaMonitorBaseAlgorithm::m_matchTool
protectedinherited

Definition at line 90 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_name

std::string AthMonitorAlgorithm::m_name
privateinherited

Definition at line 366 of file AthMonitorAlgorithm.h.

◆ m_offElectronKey

SG::ReadHandleKey<xAOD::ElectronContainer> TrigEgammaMonitorElectronAlgorithm::m_offElectronKey { this, "ElectronKey", "Electrons", ""}
private

Event Wise offline ElectronContainer Access and end iterator.

Definition at line 57 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_offProbeIsolation

Gaudi::Property<std::string> TrigEgammaMonitorElectronAlgorithm::m_offProbeIsolation { this, "OfflineProbeIsolation", "Loose"}
private

Define isolation working point for Probe electron.

Definition at line 45 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_photonIsEMTool

ToolHandleArray<IAsgPhotonIsEMSelector> TrigEgammaMonitorBaseAlgorithm::m_photonIsEMTool {this,"PhotonIsEMSelector",{}}
protectedinherited

Offline isEM Photon Selectors.

Definition at line 100 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_rmCrack

Gaudi::Property<bool> TrigEgammaMonitorElectronAlgorithm::m_rmCrack { this, "RemoveCrack", true}
private

Remove crack region for Probe default True.

Definition at line 55 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_toolLookupMap

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

Definition at line 367 of file AthMonitorAlgorithm.h.

◆ m_tools

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

Array of Generic Monitoring Tools.

Definition at line 338 of file AthMonitorAlgorithm.h.

◆ m_tp

Gaudi::Property<bool> TrigEgammaMonitorBaseAlgorithm::m_tp {this, "TPTrigger", false }
protectedinherited

TP Trigger Analysis.

Definition at line 107 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_trigdec

ToolHandle<Trig::TrigDecisionTool> TrigEgammaMonitorBaseAlgorithm::m_trigdec
privateinherited

Trigger decision tool.

Definition at line 82 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_trigDecTool

PublicToolHandle<Trig::TrigDecisionTool> AthMonitorAlgorithm::m_trigDecTool {this, "TrigDecisionTool",""}
protectedinherited

Tool to tell whether a specific trigger is passed.

Definition at line 340 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 355 of file AthMonitorAlgorithm.h.

◆ m_trigInfo

std::map<std::string,TrigInfo> TrigEgammaMonitorBaseAlgorithm::m_trigInfo
privateinherited

creates map of trigger name and TrigInfo struct

Definition at line 84 of file TrigEgammaMonitorBaseAlgorithm.h.

◆ m_trigInputList

Gaudi::Property<std::vector<std::string> > TrigEgammaMonitorElectronAlgorithm::m_trigInputList {this, "TriggerList", {}}
private

List of triggers from menu.

Definition at line 41 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_trigList

std::vector<std::string> TrigEgammaMonitorElectronAlgorithm::m_trigList
private

List of triggers to study.

Definition at line 39 of file TrigEgammaMonitorElectronAlgorithm.h.

◆ m_trigLiveFractionDataKey

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

Definition at line 347 of file AthMonitorAlgorithm.h.

◆ m_useLumi

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

Allows use of various luminosity functions.

Definition at line 359 of file AthMonitorAlgorithm.h.

◆ 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 356 of file AthMonitorAlgorithm.h.

The documentation for this class was generated from the following files:
grepfile.info
info
Definition: grepfile.py:38
TrigEgammaMonitorElectronAlgorithm::m_offProbeIsolation
Gaudi::Property< std::string > m_offProbeIsolation
Define isolation working point for Probe electron.
Definition: TrigEgammaMonitorElectronAlgorithm.h:45
TrigEgammaMonitorAnalysisAlgorithm::fillEfficiencies
void fillEfficiencies(const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * > > &, const TrigInfo &) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:42
xAOD::CaloCluster_v1::phi
virtual double phi() const
The azimuthal angle ( ) of the particle.
Definition: CaloCluster_v1.cxx:256
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
TrigEgammaMonitorBaseAlgorithm::m_photonIsEMTool
ToolHandleArray< IAsgPhotonIsEMSelector > m_photonIsEMTool
Offline isEM Photon Selectors.
Definition: TrigEgammaMonitorBaseAlgorithm.h:100
python.TIDAMonTool.monGroup
def monGroup(analysis_chain)
Definition: TIDAMonTool.py:297
TrigEgammaMonitorBaseAlgorithm::isIsolated
bool isIsolated(const xAOD::Electron *, const std::string &) const
Check if electron fulfils isolation criteria.
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:112
TrigEgammaMonitorBaseAlgorithm::dR
float dR(const float, const float, const float, const float) const
Get delta R.
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:272
et
Extra patterns decribing particle interation process.
AthMonitorAlgorithm::Environment_t::tier0Raw
@ tier0Raw
GeV
#define GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:17
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
checkFileSG.line
line
Definition: checkFileSG.py:75
python.SystemOfUnits.s
int s
Definition: SystemOfUnits.py:131
CaloCell_ID_FCS::TileExt2
@ TileExt2
Definition: FastCaloSim_CaloCell_ID.h:39
xAOD::Electron
Electron_v1 Electron
Definition of the current "egamma version".
Definition: Event/xAOD/xAODEgamma/xAODEgamma/Electron.h:17
TrigEgammaMonitorBaseAlgorithm::getTrigInfoMap
const std::map< std::string, TrigInfo > & getTrigInfoMap()
Helper methods.
Definition: TrigEgammaMonitorBaseAlgorithm.h:127
get_generator_info.result
result
Definition: get_generator_info.py:21
TrigEgammaMonitorElectronAlgorithm::m_forceVetoVeryLoose
Gaudi::Property< bool > m_forceVetoVeryLoose
Definition: TrigEgammaMonitorElectronAlgorithm.h:53
constants.EMB1
int EMB1
Definition: Calorimeter/CaloClusterCorrection/python/constants.py:53
ParticleGun_SamplingFraction.eta2
eta2
Definition: ParticleGun_SamplingFraction.py:96
PowhegControl_ttHplus_NLO.ss
ss
Definition: PowhegControl_ttHplus_NLO.py:83
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:67
TrigEgammaMonitorBaseAlgorithm::m_trigdec
ToolHandle< Trig::TrigDecisionTool > m_trigdec
Trigger decision tool.
Definition: TrigEgammaMonitorBaseAlgorithm.h:77
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
xAOD::uint8_t
uint8_t
Definition: Muon_v1.cxx:557
AthMonitorAlgorithm::m_trigDecTool
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool
Tool to tell whether a specific trigger is passed.
Definition: AthMonitorAlgorithm.h:340
ParticleTest.eg
eg
Definition: ParticleTest.py:29
_triginfo
Definition: TrigEgammaMonitorBaseAlgorithm.h:42
AthMonitorAlgorithm::Environment_t::tier0
@ tier0
AthMonitorAlgorithm::Environment_t::AOD
@ AOD
TrigEgammaMonitorBaseAlgorithm::getEt
float getEt(const xAOD::Electron *eg) const
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:291
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:83
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
TrigEgammaMonitorElectronAlgorithm::m_rmCrack
Gaudi::Property< bool > m_rmCrack
Remove crack region for Probe default True.
Definition: TrigEgammaMonitorElectronAlgorithm.h:55
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
PlotCalibFromCool.label
label
Definition: PlotCalibFromCool.py:78
CaloCell_ID_FCS::TileExt0
@ TileExt0
Definition: FastCaloSim_CaloCell_ID.h:37
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: TrackParticle_v1.cxx:77
CaloCell_ID_FCS::TileBar1
@ TileBar1
Definition: FastCaloSim_CaloCell_ID.h:32
TrigEgammaMonitorElectronAlgorithm::m_eventInfoDecorKey
SG::ReadDecorHandleKey< xAOD::EventInfo > m_eventInfoDecorKey
Definition: TrigEgammaMonitorElectronAlgorithm.h:58
accumulate
bool accumulate(AccumulateMap &map, std::vector< module_t > const &modules, FPGATrackSimMatrixAccumulator const &acc)
Accumulates an accumulator (e.g.
Definition: FPGATrackSimMatrixAccumulator.cxx:22
AthMonitorAlgorithm::m_vTrigChainNames
std::vector< std::string > m_vTrigChainNames
Vector of trigger chain names parsed from trigger chain string.
Definition: AthMonitorAlgorithm.h:356
TrigEgammaMonitorAnalysisAlgorithm::fillDistributions
void fillDistributions(const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * > > &, const TrigInfo &) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:407
AthCommonDataStore::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 > renounce(T &h)
Definition: AthCommonDataStore.h:380
AthMonitorAlgorithm::m_EventInfoKey
SG::ReadHandleKey< xAOD::EventInfo > m_EventInfoKey
Key for retrieving EventInfo from StoreGate.
Definition: AthMonitorAlgorithm.h:362
dirname
std::string dirname(std::string name)
Definition: utils.cxx:200
test_pyathena.pt
pt
Definition: test_pyathena.py:11
xAOD::passBits
passBits
Definition: TrigPassBits_v1.cxx:115
xAOD::Iso::ptcone20
@ ptcone20
Track isolation.
Definition: IsolationType.h:40
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
xAOD::CaloCluster_v1::et
double et() const
Definition: CaloCluster_v1.h:856
AthMonitorAlgorithm::m_toolLookupMap
std::unordered_map< std::string, size_t > m_toolLookupMap
Definition: AthMonitorAlgorithm.h:367
AthMonitorAlgorithm::m_environment
AthMonitorAlgorithm::Environment_t m_environment
Instance of the Environment_t enum.
Definition: AthMonitorAlgorithm.h:350
xAOD::eta1
setEt setPhi setE277 setWeta2 eta1
Definition: TrigEMCluster_v1.cxx:41
LArG4GenerateShowerLib.condition
condition
Definition: LArG4GenerateShowerLib.py:19
xAOD::numberOfTRTHits
@ numberOfTRTHits
number of TRT hits [unit8_t].
Definition: TrackingPrimitives.h:275
xAOD::Egamma_v1
Definition: Egamma_v1.h:56
vec
std::vector< size_t > vec
Definition: CombinationsGeneratorTest.cxx:12
CaloCell_ID_FCS::HEC2
@ HEC2
Definition: FastCaloSim_CaloCell_ID.h:29
TrigEgammaMonitorBaseAlgorithm::tdt
const ToolHandle< Trig::TrigDecisionTool > & tdt() const
Get the TDT
Definition: TrigEgammaMonitorBaseAlgorithm.h:133
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
TrigEgammaMonitorAnalysisAlgorithm::fillL1Calo
void fillL1Calo(const std::string &trigger, const std::vector< const xAOD::EmTauRoI * > &l1_vec) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:562
TrigEgammaMonitorAnalysisAlgorithm::fillShowerShapes
void fillShowerShapes(const std::string &trigger, const std::vector< const xAOD::Egamma * > &eg_vec, bool online) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:753
xAOD::EventInfo_v1::LAr
@ LAr
The LAr calorimeter.
Definition: EventInfo_v1.h:335
TrigEgammaMonitorElectronAlgorithm::m_trigList
std::vector< std::string > m_trigList
List of triggers to study.
Definition: TrigEgammaMonitorElectronAlgorithm.h:39
AthMonitorAlgorithm::trigChainsArePassed
bool trigChainsArePassed(const std::vector< std::string > &vTrigNames) const
Check whether triggers are passed.
Definition: AthMonitorAlgorithm.cxx:194
XMLtoHeader.count
count
Definition: XMLtoHeader.py:85
AthMonitorAlgorithm::Environment_t::user
@ user
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
xAOD::numberOfTRTHighThresholdHits
@ numberOfTRTHighThresholdHits
number of TRT hits which pass the high threshold (only xenon counted) [unit8_t].
Definition: TrackingPrimitives.h:278
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
TrigEgammaMonitorElectronAlgorithm::m_offElectronKey
SG::ReadHandleKey< xAOD::ElectronContainer > m_offElectronKey
Event Wise offline ElectronContainer Access and end iterator.
Definition: TrigEgammaMonitorElectronAlgorithm.h:57
Monitored::Collection
ValuesCollection< T > Collection(std::string name, const T &collection)
Declare a monitored (double-convertible) collection.
Definition: MonitoredCollection.h:38
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
AthMonitorAlgorithm::m_dummy
const ToolHandle< GenericMonitoringTool > m_dummy
Definition: AthMonitorAlgorithm.h:369
TrigEgammaMonitorBaseAlgorithm::setAccept
asg::AcceptData setAccept(const TrigCompositeUtils::Decision *, const TrigInfo &) const
Set the accept object for all trigger levels.
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:177
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
TrigEgammaMonitorAnalysisAlgorithm::fillResolutions
void fillResolutions(const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const TrigInfo &info) const
*****************************************************************************************************...
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:914
xAOD::CaloCluster_v1::etaBE
float etaBE(const unsigned layer) const
Get the eta in one layer of the EM Calo.
Definition: CaloCluster_v1.cxx:644
TrigEgammaMonitorAnalysisAlgorithm::fillL2Electron
void fillL2Electron(const std::string &trigger, const std::vector< const xAOD::TrigElectron * > &el_vec) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:649
PyPoolBrowser.item
item
Definition: PyPoolBrowser.py:129
TrigEgammaMonitorBaseAlgorithm::isHLTTruncated
bool isHLTTruncated() const
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:733
calibdata.valid
list valid
Definition: calibdata.py:45
LArEventBitInfo::NOISEBURSTVETO
@ NOISEBURSTVETO
Definition: LArEventBitInfo.h:13
TileSynchronizeBch.online
online
Definition: TileSynchronizeBch.py:88
TrigEgammaMonitorBaseAlgorithm::getTrigInfo
TrigInfo getTrigInfo(const std::string &) const
Get the trigger info parsed from the chain name (only single lepton triggers)
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:583
TrigEgammaMonitorBaseAlgorithm::m_accept
asg::AcceptInfo m_accept
Definition: TrigEgammaMonitorBaseAlgorithm.h:122
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:59
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
xAOD::EgammaHelpers::isConvertedPhoton
bool isConvertedPhoton(const xAOD::Egamma *eg, bool excludeTRT=false)
is the object a converted photon
Definition: EgammaxAODHelpers.cxx:25
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthMonitorAlgorithm::m_DQFilterTools
ToolHandleArray< IDQFilterTool > m_DQFilterTools
Array of Data Quality filter tools.
Definition: AthMonitorAlgorithm.h:341
CaloCell_ID_FCS::HEC1
@ HEC1
Definition: FastCaloSim_CaloCell_ID.h:28
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
PixelAthClusterMonAlgCfg.histname
histname
Definition: PixelAthClusterMonAlgCfg.py:106
TrigEgammaMonitorAnalysisAlgorithm::fillL1CaloResolution
void fillL1CaloResolution(const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const xAOD::EmTauRoI * >> &pairObjs) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:970
CheckAppliedSFs.e3
e3
Definition: CheckAppliedSFs.py:264
TrigEgammaMonitorBaseAlgorithm::m_trigInfo
std::map< std::string, TrigInfo > m_trigInfo
creates map of trigger name and TrigInfo struct
Definition: TrigEgammaMonitorBaseAlgorithm.h:84
TrigEgammaMonitorBaseAlgorithm::m_defaultProbePidElectron
Gaudi::Property< std::string > m_defaultProbePidElectron
default probe pid for electron trigitems that don't have pid in their name
Definition: TrigEgammaMonitorBaseAlgorithm.h:109
SG::Decorator< bool >
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
xAOD::CaloCluster_v1::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition: CaloCluster_v1.cxx:251
TrigEgammaMonitorAnalysisAlgorithm::TrigEgammaMonitorAnalysisAlgorithm
TrigEgammaMonitorAnalysisAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:8
TrigEgammaMonitorBaseAlgorithm::setTrigInfo
void setTrigInfo(const std::string &)
Set the trigger info parsed from the chain name.
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:591
CaloCell_ID_FCS::TileBar0
@ TileBar0
Definition: FastCaloSim_CaloCell_ID.h:31
TrigEgammaMonitorAnalysisAlgorithm::fillL2CaloResolution
void fillL2CaloResolution(const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:1711
AthMonitorAlgorithm::DataType_t::heavyIonCollisions
@ heavyIonCollisions
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
xAOD::EgammaHelpers::isElectron
bool isElectron(const xAOD::Egamma *eg)
is the object an electron (not Fwd)
Definition: EgammaxAODHelpers.cxx:12
TrigEgammaMonitorAnalysisAlgorithm::fillHLTPhotonResolution
void fillHLTPhotonResolution(const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const TrigInfo &info) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:1406
tolower
void tolower(std::string &s)
Definition: AthenaSummarySvc.cxx:111
AthMonitorAlgorithm::fillHistograms
virtual StatusCode fillHistograms(const EventContext &ctx) const =0
adds event to the monitoring histograms
TrigEgammaMonitorBaseAlgorithm::m_emulatorTool
ToolHandle< Trig::TrigEgammaEmulationToolMT > m_emulatorTool
Definition: TrigEgammaMonitorBaseAlgorithm.h:92
contains
bool contains(const std::string &s, const std::string &regx)
does a string contain the substring
Definition: hcg.cxx:111
xAOD::Egamma_v1::caloCluster
const xAOD::CaloCluster * caloCluster(size_t index=0) const
Pointer to the xAOD::CaloCluster/s that define the electron candidate.
Definition: Egamma_v1.cxx:388
Amg::transform
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Definition: GeoPrimitivesHelpers.h:156
TRT::Track::d0
@ d0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:62
ParticleGun_EoverP_Config.pid
pid
Definition: ParticleGun_EoverP_Config.py:62
plotIsoValidation.el
el
Definition: plotIsoValidation.py:197
constants.EME1
int EME1
Definition: Calorimeter/CaloClusterCorrection/python/constants.py:55
python.xAODType.dummy
dummy
Definition: xAODType.py:4
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
TrigEgammaMonitorBaseAlgorithm::m_doEmulation
Gaudi::Property< bool > m_doEmulation
Do emulation.
Definition: TrigEgammaMonitorBaseAlgorithm.h:104
xAOD::Egamma_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition: Egamma_v1.cxx:75
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
AthMonitorAlgorithm::fill
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.
python.BuildSignatureFlags.cosmics
AthConfigFlags cosmics(AthConfigFlags flags, str instanceName, str recoMode)
Definition: BuildSignatureFlags.py:524
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
TrigEgammaMonitorBaseAlgorithm::m_electronIsEMTool
ToolHandleArray< IAsgElectronIsEMSelector > m_electronIsEMTool
Offline isEM Selectors.
Definition: TrigEgammaMonitorBaseAlgorithm.h:94
TRT::Track::z0
@ z0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:63
TrigEgammaMonitorAnalysisAlgorithm::fillL1CaloAbsResolution
void fillL1CaloAbsResolution(const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const xAOD::EmTauRoI * >> &pairObjs) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:1000
xAOD::TrigComposite_v1
Class used to describe composite objects in the HLT.
Definition: TrigComposite_v1.h:52
TrigEgammaMonitorBaseAlgorithm::m_electronLHTool
ToolHandleArray< IAsgElectronLikelihoodTool > m_electronLHTool
Offline LH Selectors.
Definition: TrigEgammaMonitorBaseAlgorithm.h:96
AthMonitorAlgorithm::GetEventInfo
SG::ReadHandle< xAOD::EventInfo > GetEventInfo(const EventContext &) const
Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
Definition: AthMonitorAlgorithm.cxx:107
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
TrigEgammaMonitorAnalysisAlgorithm::fillL2Calo
void fillL2Calo(const std::string &trigger, const std::vector< const xAOD::TrigEMCluster * > &emCluster_vec) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:622
TrigEgammaMonitorElectronAlgorithm::executeNavigation
StatusCode executeNavigation(const EventContext &ctx, const std::string &trigItem, float, const std::string &, std::vector< std::pair< std::shared_ptr< const xAOD::Egamma >, const TrigCompositeUtils::Decision * >> &) const
navigation method called by executeNavigation
Definition: TrigEgammaMonitorElectronAlgorithm.cxx:103
AthMonitorAlgorithm::m_dataType
AthMonitorAlgorithm::DataType_t m_dataType
Instance of the DataType_t enum.
Definition: AthMonitorAlgorithm.h:351
AthReentrantAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Definition: AthReentrantAlgorithm.h:153
CxxUtils::atof
double atof(std::string_view str)
Converts a string into a double / float.
Definition: Control/CxxUtils/Root/StringUtils.cxx:91
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
AthMonitorAlgorithm::DataType_t::cosmics
@ cosmics
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
TrigEgammaMonitorBaseAlgorithm::m_isemname
Gaudi::Property< std::vector< std::string > > m_isemname
isem names
Definition: TrigEgammaMonitorBaseAlgorithm.h:113
AthMonitorAlgorithm::m_tools
ToolHandleArray< GenericMonitoringTool > m_tools
Array of Generic Monitoring Tools.
Definition: AthMonitorAlgorithm.h:338
TrigEgammaMonitorBaseAlgorithm::ApplyElectronPid
bool ApplyElectronPid(const xAOD::Electron *eg, const std::string &) const
Get offline electron decision.
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:49
AthMonitorAlgorithm::Environment_t::online
@ online
TrigEgammaMonitorElectronAlgorithm::m_trigInputList
Gaudi::Property< std::vector< std::string > > m_trigInputList
List of triggers from menu.
Definition: TrigEgammaMonitorElectronAlgorithm.h:41
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
TrigEgammaMonitorAnalysisAlgorithm::initialize
virtual StatusCode initialize() override
initialize
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:17
TrigEgammaMonitorAnalysisAlgorithm::fillEFCalo
void fillEFCalo(const std::string &trigger, const std::vector< const xAOD::CaloCluster * > &clus_vec) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:698
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
TrigEgammaMonitorAnalysisAlgorithm::fillL1eEM
void fillL1eEM(const std::string &trigger, const std::vector< const xAOD::eFexEMRoI * > &l1_vec) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:591
CaloCell_ID_FCS::TileExt1
@ TileExt1
Definition: FastCaloSim_CaloCell_ID.h:38
CaloCell_ID_FCS::EME3
@ EME3
Definition: FastCaloSim_CaloCell_ID.h:26
AthMonitorAlgorithm::Environment_t::tier0ESD
@ tier0ESD
threshold
Definition: chainparser.cxx:74
AthMonitorAlgorithm::Environment_t::altprod
@ altprod
AthReentrantAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthReentrantAlgorithm.cxx:96
AthMonitorAlgorithm::m_enforceExpressTriggers
Gaudi::Property< bool > m_enforceExpressTriggers
Definition: AthMonitorAlgorithm.h:372
item
Definition: ItemListSvc.h:43
WriteBchToCool.user
user
Definition: WriteBchToCool.py:76
TrigEgammaMonitorBaseAlgorithm::m_matchTool
ToolHandle< TrigEgammaMatchingToolMT > m_matchTool
Definition: TrigEgammaMonitorBaseAlgorithm.h:90
xAOD::Electron_v1
Definition: Electron_v1.h:34
xAOD::CaloCluster_v1::PHICALOFRAME
@ PHICALOFRAME
Phi in the calo frame (for egamma)
Definition: CaloCluster_v1.h:188
AthMonitorAlgorithm::lbAverageInteractionsPerCrossing
virtual float lbAverageInteractionsPerCrossing(const EventContext &ctx=Gaudi::Hive::currentContext()) const
Calculate the average mu, i.e.
Definition: AthMonitorAlgorithm.cxx:222
TrigEgammaMonitorAnalysisAlgorithm::fillTracking
void fillTracking(const std::string &trigger, const std::vector< const xAOD::Electron * > &eg_vec, bool online) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:816
TrigEgammaMonitorBaseAlgorithm::match
const ToolHandle< TrigEgammaMatchingToolMT > & match() const
Get the e/g match tool.
Definition: TrigEgammaMonitorBaseAlgorithm.h:135
xAOD::Photon
Photon_v1 Photon
Definition of the current "egamma version".
Definition: Event/xAOD/xAODEgamma/xAODEgamma/Photon.h:17
CaloCell_ID_FCS::HEC0
@ HEC0
Definition: FastCaloSim_CaloCell_ID.h:27
TrigEgammaMonitorAnalysisAlgorithm::fillLabel
void fillLabel(const ToolHandle< GenericMonitoringTool > &groupHandle, const std::string &histname, const std::string &label) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:30
ATH_LIKELY
#define ATH_LIKELY(x)
Definition: AthUnlikelyMacros.h:16
TrigEgammaMonitorBaseAlgorithm::m_lhname
Gaudi::Property< std::vector< std::string > > m_lhname
lh names
Definition: TrigEgammaMonitorBaseAlgorithm.h:115
TrigEgammaMonitorBaseAlgorithm::m_electronDNNTool
ToolHandleArray< IAsgElectronLikelihoodTool > m_electronDNNTool
Offline DNN Selectors.
Definition: TrigEgammaMonitorBaseAlgorithm.h:98
a
TList * a
Definition: liststreamerinfos.cxx:10
TrigEgammaMonitorBaseAlgorithm::getD0sig
float getD0sig(const xAOD::Electron *eg) const
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:362
CaloLCW_tf.group
group
Definition: CaloLCW_tf.py:28
h
python.CaloScaleNoiseConfig.str
str
Definition: CaloScaleNoiseConfig.py:78
xAOD::Photon_v1
Definition: Photon_v1.h:37
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
AthMonitorAlgorithm::DataType_t::collisions
@ collisions
TrigRoiDescriptorCollection
Definition: TrigRoiDescriptorCollection.h:21
Pythia8_RapidityOrderMPI.val
val
Definition: Pythia8_RapidityOrderMPI.py:14
AthMonitorAlgorithm::m_name
std::string m_name
Definition: AthMonitorAlgorithm.h:366
plotBeamSpotMon.mon
mon
Definition: plotBeamSpotMon.py:67
xAOD::CaloCluster_v1::energyBE
float energyBE(const unsigned layer) const
Get the energy in one layer of the EM Calo.
Definition: CaloCluster_v1.cxx:630
TrigEgammaMonitorElectronAlgorithm::m_forcePidSelection
Gaudi::Property< bool > m_forcePidSelection
force pid and crack selection into electron navigation
Definition: TrigEgammaMonitorElectronAlgorithm.h:49
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
TrigEgammaMonitorBaseAlgorithm::m_detailedHists
Gaudi::Property< bool > m_detailedHists
Include more detailed histograms.
Definition: TrigEgammaMonitorBaseAlgorithm.h:119
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
xAOD::Egamma_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: Egamma_v1.cxx:65
doL1CaloHVCorrections.parts
parts
Definition: doL1CaloHVCorrections.py:334
xAOD::Egamma_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: Egamma_v1.cxx:70
str
Definition: BTagTrackIpAccessor.cxx:11
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
TrigEgammaMonitorAnalysisAlgorithm::fillInefficiency
void fillInefficiency(const std::string &pidword, const TrigInfo &info, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const std::vector< asg::AcceptData > &acceptObjs) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:364
TrigEgammaMonitorElectronAlgorithm::m_forceProbeIsolation
Gaudi::Property< bool > m_forceProbeIsolation
force probe isolation
Definition: TrigEgammaMonitorElectronAlgorithm.h:47
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
veto
std::vector< std::string > veto
these patterns are anded
Definition: listroot.cxx:191
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
skel.l1
l1
Definition: skel.GENtoEVGEN.py:398
Monitored::Scalar
Declare a monitored scalar variable.
Definition: MonitoredScalar.h:34
TrigEgammaMonitorElectronAlgorithm::m_forceEtThr
Gaudi::Property< bool > m_forceEtThr
Definition: TrigEgammaMonitorElectronAlgorithm.h:51
xAOD::TrigEMCluster_v1
Description of a trigger EM cluster.
Definition: TrigEMCluster_v1.h:28
AthMonitorAlgorithm::DataType_t::userDefined
@ userDefined
xAOD::EgammaParameters::e2tsts1
@ e2tsts1
energy of the cell corresponding to second energy maximum in the first sampling
Definition: EgammaEnums.h:108
CaloCell_ID_FCS::HEC3
@ HEC3
Definition: FastCaloSim_CaloCell_ID.h:30
xAOD::CaloCluster_v1::ETACALOFRAME
@ ETACALOFRAME
Eta in the calo frame (for egamma)
Definition: CaloCluster_v1.h:187
asg::AcceptData
Definition: AcceptData.h:30
CaloCell_ID_FCS::EMB3
@ EMB3
Definition: FastCaloSim_CaloCell_ID.h:22
TrigEgammaMonitorBaseAlgorithm::m_defaultProbePidPhoton
Gaudi::Property< std::string > m_defaultProbePidPhoton
default probe pid for photon trigitems that don't have pid in their name
Definition: TrigEgammaMonitorBaseAlgorithm.h:111
makeComparison.deltaR
float deltaR
Definition: makeComparison.py:36
TrigEgammaMonitorAnalysisAlgorithm::fillL2Photon
void fillL2Photon(const std::string &trigger, const std::vector< const xAOD::TrigPhoton * > &eg_vec) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:673
CaloCell_ID_FCS::TileBar2
@ TileBar2
Definition: FastCaloSim_CaloCell_ID.h:33
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
xAOD::EgammaParameters::emaxs1
@ emaxs1
energy of strip with maximal energy deposit
Definition: EgammaEnums.h:145
WriteBchToCool.onl
onl
Definition: WriteBchToCool.py:137
xAOD::CaloCluster_v1::e
virtual double e() const
The total energy of the particle.
Definition: CaloCluster_v1.cxx:265
dumpTgcDigiThreshold.threshold
list threshold
Definition: dumpTgcDigiThreshold.py:34
Trk::split
@ split
Definition: LayerMaterialProperties.h:38
TrigEgammaMonitorBaseAlgorithm::getL1Item
std::string getL1Item(const std::string &trigger) const
Creates static map to return L1 item from trigger name.
Definition: TrigEgammaMonitorBaseAlgorithm.cxx:724
VP1PartSpect::Pi
@ Pi
Definition: VP1PartSpectFlags.h:27
AthMonitorAlgorithm::getGroup
const ToolHandle< GenericMonitoringTool > & getGroup(const std::string &name) const
Get a specific monitoring tool from the tool handle array.
Definition: AthMonitorAlgorithm.cxx:164
AthMonitorAlgorithm::DataType_t::monteCarlo
@ monteCarlo
TrigEgammaMonitorElectronAlgorithm::m_electronIsolationKeyPtCone20
SG::ReadDecorHandleKey< xAOD::ElectronContainer > m_electronIsolationKeyPtCone20
Ensure offline electron isolation decoration is retrieved after being created.
Definition: TrigEgammaMonitorElectronAlgorithm.h:60
fitman.k
k
Definition: fitman.py:528
TrigEgammaMonitorAnalysisAlgorithm::fillHLTElectronResolution
void fillHLTElectronResolution(const std::string &trigger, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const TrigInfo &info) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:1028
ServiceHandle< ICondSvc >
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
TrigEgammaMonitorAnalysisAlgorithm::fillEfficiency
void fillEfficiency(const std::string &subgroup, const std::string &level, const std::string &pidword, const TrigInfo &info, const std::vector< std::pair< const xAOD::Egamma *, const TrigCompositeUtils::Decision * >> &pairObjs, const std::vector< asg::AcceptData > &acceptObjs, const std::string &dirname) const
Definition: TrigEgammaMonitorAnalysisAlgorithm.cxx:156
collisions
Definition: collisions.py:1
Generated on Sun Dec 22 2024 21:50:15 for ATLAS Offline Software by   doxygen 1.8.18