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

#include <MonitorTnPAlgorithm.h>

Inheritance diagram for MonitorTnPAlgorithm:
Collaboration diagram for MonitorTnPAlgorithm:

Public Types

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

Public Member Functions

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

Protected Member Functions

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

Protected Attributes

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

Private Types

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

Private Member Functions

Regions GetRegion (Float_t eta) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadHandleKey< xAOD::ElectronContainerm_ParticleContainerKey {this, "ParticleContainerName", "Electrons", "Name of electron container" }
SG::ReadDecorHandleKeyArray< xAOD::ElectronContainerm_ElectronIsolationKey {this, "ElectronIsolationKey", {"Electrons.ptcone20", "Electrons.topoetcone40"} }
Gaudi::Property< std::string > m_ParticlePrefix {this,"ParticlePrefix","electron","Name of the particle prefix to be used to define hists"}
Gaudi::Property< std::string > m_RecoName {this,"RecoName","LHLoose","Name of particle flavor in egamma reco"}
Gaudi::Property< std::string > m_TnPType {this,"TnPType","Z","Variable to describe what is the TnP resonance name"}
Gaudi::Property< Float_t > m_MassPeak {this, "MassPeak", ParticleConstants::ZMassInMeV, "Resonance peak position"}
Gaudi::Property< Float_t > m_ElectronEtCut {this, "ElectronEtCut", 15000, "Et cut for electron"}
Gaudi::Property< Float_t > m_MassLowerCut {this, "MassLowerCut", 70000, "Lower Mass cut"}
Gaudi::Property< Float_t > m_MassUpperCut {this, "MassUpperCut", 110000, "Upper Mass Cut"}
std::string m_name
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringToolm_dummy
Gaudi::Property< bool > m_enforceExpressTriggers
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 23 of file MonitorTnPAlgorithm.h.

Member Typedef Documentation

◆ MonVarVec_t

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

Definition at line 370 of file AthMonitorAlgorithm.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ DataType_t

enum class AthMonitorAlgorithm::DataType_t
stronginherited

Specifies what type of input data is being monitored.

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

Enumerator
userDefined 
monteCarlo 
collisions 
cosmics 
heavyIonCollisions 

Definition at line 194 of file AthMonitorAlgorithm.h.

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

◆ Environment_t

enum class AthMonitorAlgorithm::Environment_t
stronginherited

Specifies the processing environment.

The running environment may be used to select which histograms are produced, and where they are located in the output. For example, the output paths of the histograms are different for the "user", "online" and the various offline flags. Strings of the same names may be given as jobOptions.

Enumerator
user 
online 
tier0 
tier0Raw 
tier0ESD 
AOD 
altprod 

Definition at line 175 of file AthMonitorAlgorithm.h.

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

◆ Regions

Constructor & Destructor Documentation

◆ MonitorTnPAlgorithm()

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

Definition at line 10 of file MonitorTnPAlgorithm.cxx.

11:AthMonitorAlgorithm(name,pSvcLocator) {
12}
AthMonitorAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.

◆ ~MonitorTnPAlgorithm()

virtual MonitorTnPAlgorithm::~MonitorTnPAlgorithm ( )
virtualdefault

Member Function Documentation

◆ cardinality()

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

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

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ dataType()

DataType_t AthMonitorAlgorithm::dataType ( ) const
inlineinherited

Accessor functions for the data type.

Returns
the current value of the class's DataType_t instance.

Definition at line 224 of file AthMonitorAlgorithm.h.

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

◆ dataTypeStringToEnum()

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

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

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

Definition at line 144 of file AthMonitorAlgorithm.cxx.

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

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

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

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ environment()

Environment_t AthMonitorAlgorithm::environment ( ) const
inlineinherited

Accessor functions for the environment.

Returns
the current value of the class's Environment_t instance.

Definition at line 208 of file AthMonitorAlgorithm.h.

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

◆ envStringToEnum()

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

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

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

Definition at line 116 of file AthMonitorAlgorithm.cxx.

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

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

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

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters
ctxevent context for reentrant Athena call
Returns
StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

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

◆ extraDeps_update_handler()

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

Add StoreName to extra input/output deps as needed.

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

◆ extraOutputDeps()

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

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

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

◆ fillElectronProbe()

StatusCode MonitorTnPAlgorithm::fillElectronProbe ( const xAOD::Electron * el,
const bool isleading,
const EventContext & ctx ) const

Definition at line 154 of file MonitorTnPAlgorithm.cxx.

154 {
155 using namespace Monitored;
156
157 auto lbNProbeCandidates = Monitored::Scalar<u_int16_t>("LBEvoN",0);
158
159 lbNProbeCandidates = GetEventInfo(ctx)->lumiBlock();
160
161 if (isleading) {
162 auto EtaZ = Monitored::Scalar<Float_t>("EtaZ",0.0);
163 auto PhiZ = Monitored::Scalar<Float_t>("PhiZ",0.0);
164 EtaZ = el->eta();
165 PhiZ = el->phi();
166 fill("MonitorTnP",EtaZ,PhiZ);
167 }
168
169 auto is_TightID = Monitored::Scalar<bool>("is_Tight",false);
170
171 bool lead_isLHTight = false;
172 if ( !el->passSelection(lead_isLHTight,"LHTight") ) return StatusCode::FAILURE;
173 is_TightID = lead_isLHTight;
174
175 auto is_Iso = Monitored::Scalar<bool>("is_Iso",false);
176 auto is_IsoandTight = Monitored::Scalar<bool>("is_IsoandTight",false);
177
178 Float_t myet = el->pt(); // in MeV (/Gaudi::Units::GeV; // in GeV)
179 Float_t myptcone20 = -999.;
180 Float_t mytopoetcone40 = -999.;
181
182 bool isIsolated = false;
183 if ( ! el->isolationValue(mytopoetcone40,xAOD::Iso::topoetcone40) ) return StatusCode::FAILURE;
184
185 if (myet>0) {
186 if ( ! el->isolationValue(myptcone20,xAOD::Iso::ptcone20) ) return StatusCode::FAILURE;
187 if ( myptcone20/myet < 0.15 ) isIsolated = true;
188 }
189
190 is_Iso = isIsolated;
191 is_IsoandTight = isIsolated && lead_isLHTight;
192
193 auto et = Monitored::Scalar<Float_t>("Etprobe",0.0);
194 auto eta = Monitored::Scalar<Float_t>("Etaprobe",0.0);
195 auto phi = Monitored::Scalar<Float_t>("Phiprobe",0.0);
196 auto time = Monitored::Scalar<Float_t>("Timeprobe",0.0);
197 auto topoetcone40 = Monitored::Scalar<Float_t>("TopoEtCone40probe",0.0);
198 auto ptcone20 = Monitored::Scalar<Float_t>("PtCone20probe",0.0);
199
200 // Particle variables per Region
201
202 // BARREL
203
204 auto et_barrel = Monitored::Scalar<Float_t>("EtProbeinBARREL",0.0);
205 auto eta_barrel = Monitored::Scalar<Float_t>("EtaProbeinBARREL",0.0);
206 auto phi_barrel = Monitored::Scalar<Float_t>("PhiProbeinBARREL",0.0);
207 auto time_barrel = Monitored::Scalar<Float_t>("TimeProbeinBARREL",0.0);
208 auto ehad1_barrel = Monitored::Scalar<Float_t>("Ehad1ProbeinBARREL",0.0);
209 auto eoverp_barrel = Monitored::Scalar<Float_t>("EoverPProbeinBARREL",0.0);
210 auto coreem_barrel = Monitored::Scalar<Float_t>("CoreEMProbeinBARREL",0.0);
211 auto f0_barrel = Monitored::Scalar<Float_t>("F0ProbeinBARREL",0.0);
212 auto f1_barrel = Monitored::Scalar<Float_t>("F1ProbeinBARREL",0.0);
213 auto f2_barrel = Monitored::Scalar<Float_t>("F2ProbeinBARREL",0.0);
214 auto f3_barrel = Monitored::Scalar<Float_t>("F3ProbeinBARREL",0.0);
215 auto re233e237_barrel = Monitored::Scalar<Float_t>("Re233e237ProbeinBARREL",0.0);
216 auto re237e277_barrel = Monitored::Scalar<Float_t>("Re237e277ProbeinBARREL",0.0);
217 auto nofblayerhits_barrel = Monitored::Scalar<u_int8_t>("NOfBLayerHitsProbeinBARREL",0);
218 auto nofpixelhits_barrel = Monitored::Scalar<u_int8_t>("NOfPixelHitsProbeinBARREL",0);
219 auto nofscthits_barrel = Monitored::Scalar<u_int8_t>("NOfSCTHitsProbeinBARREL",0);
220 auto noftrthits_barrel = Monitored::Scalar<u_int8_t>("NOfTRTHitsProbeinBARREL",0);
221 auto noftrthighthresholdhits_barrel = Monitored::Scalar<u_int8_t>("NOfTRTHighThresholdHitsProbeinBARREL",0);
222 auto deltaeta1_barrel = Monitored::Scalar<Float_t>("DeltaEta1ProbeinBARREL",0.0);
223 auto deltaphi2_barrel = Monitored::Scalar<Float_t>("DeltaPhi2ProbeinBARREL",0.0);
224 auto trackd0_barrel = Monitored::Scalar<Float_t>("Trackd0ProbeinBARREL",0.0);
225
226 // ENDCAP
227
228 auto et_endcap = Monitored::Scalar<Float_t>("EtProbeinENDCAP",0.0);
229 auto eta_endcap = Monitored::Scalar<Float_t>("EtaProbeinENDCAP",0.0);
230 auto phi_endcap = Monitored::Scalar<Float_t>("PhiProbeinENDCAP",0.0);
231 auto time_endcap = Monitored::Scalar<Float_t>("TimeProbeinENDCAP",0.0);
232 auto ehad1_endcap = Monitored::Scalar<Float_t>("Ehad1ProbeinENDCAP",0.0);
233 auto eoverp_endcap = Monitored::Scalar<Float_t>("EoverPProbeinENDCAP",0.0);
234 auto coreem_endcap = Monitored::Scalar<Float_t>("CoreEMProbeinENDCAP",0.0);
235 auto f0_endcap = Monitored::Scalar<Float_t>("F0ProbeinENDCAP",0.0);
236 auto f1_endcap = Monitored::Scalar<Float_t>("F1ProbeinENDCAP",0.0);
237 auto f2_endcap = Monitored::Scalar<Float_t>("F2ProbeinENDCAP",0.0);
238 auto f3_endcap = Monitored::Scalar<Float_t>("F3ProbeinENDCAP",0.0);
239 auto re233e237_endcap = Monitored::Scalar<Float_t>("Re233e237ProbeinENDCAP",0.0);
240 auto re237e277_endcap = Monitored::Scalar<Float_t>("Re237e277ProbeinENDCAP",0.0);
241 auto nofblayerhits_endcap = Monitored::Scalar<u_int8_t>("NOfBLayerHitsProbeinENDCAP",0);
242 auto nofpixelhits_endcap = Monitored::Scalar<u_int8_t>("NOfPixelHitsProbeinENDCAP",0);
243 auto nofscthits_endcap = Monitored::Scalar<u_int8_t>("NOfSCTHitsProbeinENDCAP",0);
244 auto noftrthits_endcap = Monitored::Scalar<u_int8_t>("NOfTRTHitsProbeinENDCAP",0);
245 auto noftrthighthresholdhits_endcap = Monitored::Scalar<u_int8_t>("NOfTRTHighThresholdHitsProbeinENDCAP",0);
246 auto deltaeta1_endcap = Monitored::Scalar<Float_t>("DeltaEta1ProbeinENDCAP",0.0);
247 auto deltaphi2_endcap = Monitored::Scalar<Float_t>("DeltaPhi2ProbeinENDCAP",0.0);
248 auto trackd0_endcap = Monitored::Scalar<Float_t>("Trackd0ProbeinENDCAP",0.0);
249
250 // CRACK
251
252 auto et_crack = Monitored::Scalar<Float_t>("EtProbeinCRACK",0.0);
253 auto eta_crack = Monitored::Scalar<Float_t>("EtaProbeinCRACK",0.0);
254 auto phi_crack = Monitored::Scalar<Float_t>("PhiProbeinCRACK",0.0);
255 auto time_crack = Monitored::Scalar<Float_t>("TimeProbeinCRACK",0.0);
256 auto ehad1_crack = Monitored::Scalar<Float_t>("Ehad1ProbeinCRACK",0.0);
257 auto eoverp_crack = Monitored::Scalar<Float_t>("EoverPProbeinCRACK",0.0);
258 auto coreem_crack = Monitored::Scalar<Float_t>("CoreEMProbeinCRACK",0.0);
259 auto f0_crack = Monitored::Scalar<Float_t>("F0ProbeinCRACK",0.0);
260 auto f1_crack = Monitored::Scalar<Float_t>("F1ProbeinCRACK",0.0);
261 auto f2_crack = Monitored::Scalar<Float_t>("F2ProbeinCRACK",0.0);
262 auto f3_crack = Monitored::Scalar<Float_t>("F3ProbeinCRACK",0.0);
263 auto re233e237_crack = Monitored::Scalar<Float_t>("Re233e237ProbeinCRACK",0.0);
264 auto re237e277_crack = Monitored::Scalar<Float_t>("Re237e277ProbeinCRACK",0.0);
265 auto nofblayerhits_crack = Monitored::Scalar<u_int8_t>("NOfBLayerHitsProbeinCRACK",0);
266 auto nofpixelhits_crack = Monitored::Scalar<u_int8_t>("NOfPixelHitsProbeinCRACK",0);
267 auto nofscthits_crack = Monitored::Scalar<u_int8_t>("NOfSCTHitsProbeinCRACK",0);
268 auto noftrthits_crack = Monitored::Scalar<u_int8_t>("NOfTRTHitsProbeinCRACK",0);
269 auto noftrthighthresholdhits_crack = Monitored::Scalar<u_int8_t>("NOfTRTHighThresholdHitsProbeinCRACK",0);
270 auto deltaeta1_crack = Monitored::Scalar<Float_t>("DeltaEta1ProbeinCRACK",0.0);
271 auto deltaphi2_crack = Monitored::Scalar<Float_t>("DeltaPhi2ProbeinCRACK",0.0);
272 auto trackd0_crack = Monitored::Scalar<Float_t>("Trackd0ProbeinCRACK",0.0);
273
274
275 Float_t myeta = el->eta();
276 Float_t myphi = el->phi();
277
278 auto is_pt_gt_4gev = Monitored::Scalar<bool>("is_pt_gt_4gev",false);
279 is_pt_gt_4gev = myet > 4000. ;
280
281 auto is_pt_gt_4gev_barrel = Monitored::Scalar<bool>("is_pt_gt_4gevBARREL",false);
282 is_pt_gt_4gev_barrel = myet > 4000. ;
283
284 auto is_pt_gt_4gev_endcap = Monitored::Scalar<bool>("is_pt_gt_4gevENDCAP",false);
285 is_pt_gt_4gev_endcap = myet > 4000. ;
286
287 auto is_pt_gt_4gev_crack = Monitored::Scalar<bool>("is_pt_gt_4gevCRACK",false);
288 is_pt_gt_4gev_crack = myet > 4000. ;
289
290 // Isolation Energy
291
292 el->isolationValue(myptcone20,xAOD::Iso::ptcone20);
293 ptcone20 = myptcone20;
294
295 // Shower shape variable details
296 Float_t myehad1 = 0.0;
297 Float_t myecore = 0.0;
298 Float_t myf0 = 0.0;
299 Float_t myf1 = 0.0;
300 Float_t myf2 = 0.0;
301 Float_t myf3 = 0.0;
302 Float_t e233 = 0.0;
303 Float_t e237 = 0.0;
304 Float_t e277 = 0.0;
305 Float_t myre233e237 = 0.0;
306 Float_t myre237e277 = 0.0;
307
308 el->showerShapeValue(myehad1, xAOD::EgammaParameters::ehad1);
309 el->showerShapeValue(myecore, xAOD::EgammaParameters::ecore);
310
311 el->showerShapeValue(e237, xAOD::EgammaParameters::e237);
312 el->showerShapeValue(e233, xAOD::EgammaParameters::e233);
313 el->showerShapeValue(e277, xAOD::EgammaParameters::e277);
314
315 if (e237!=0) myre233e237 = e233 / e237;
316 if (e277!=0) myre237e277 = e237 / e277;
317
318 // Cluster track match details
319 Float_t mydeltaeta1 = -999.0;
320 el->trackCaloMatchValue(mydeltaeta1, xAOD::EgammaParameters::deltaEta1);
321 Float_t mydeltaphi2 = -999.0;
322 el->trackCaloMatchValue(mydeltaphi2, xAOD::EgammaParameters::deltaPhi2);
323
324 Float_t myeoverp = -999.0;
325 Float_t mytrackd0 = -999.0;
326
327 // associated track details
328 const xAOD::TrackParticle *t = el->trackParticle();
329 double trackp = 0;
330 u_int8_t mynofblayerhits=-1;
331 u_int8_t mynofpixelhits=-1;
332 u_int8_t mynofscthits=-1;
333 u_int8_t mynoftrthits=-1;
334 u_int8_t mynoftrthighthresholdhits=-1;
335
336 if (t) {
337 trackp = t->pt()*cosh(t->eta());
338 // retrieve track summary information
339 t->summaryValue(mynofblayerhits,xAOD::numberOfInnermostPixelLayerHits);
340 t->summaryValue(mynofpixelhits,xAOD::numberOfPixelHits);
341 t->summaryValue(mynofscthits,xAOD::numberOfSCTHits);
342 t->summaryValue(mynoftrthits,xAOD::numberOfTRTHits);
343 t->summaryValue(mynoftrthighthresholdhits,xAOD::numberOfTRTHighThresholdHits);
344 mytrackd0 = t->d0();
345 }
346 else
347 ATH_MSG_DEBUG("Error attempting to retrieve associated track");
348
349 // Associated cluster details
350
351 Float_t mytime=0.;
352 const xAOD::CaloCluster *aCluster = el->caloCluster();
353 if (aCluster) {
354 mytime = aCluster->time();
355 time = mytime ;
356 // Shower shape variable details
357 double ec = aCluster->et()*cosh(aCluster->eta());
358 if (ec!=0) myf0 = aCluster->energyBE(0)/ec;
359 if (ec!=0) myf1 = aCluster->energyBE(1)/ec;
360 if (ec!=0) myf2 = aCluster->energyBE(2)/ec;
361 if (ec!=0) myf3 = aCluster->energyBE(3)/ec;
362
363 if(trackp !=0) myeoverp = aCluster->e()/trackp;
364
365 //ATH_MSG_WARNING("Test eoverp : " << myeoverp);
366 }
367
368 // Fill per region histograms
369 auto region = GetRegion(myeta);
370 switch(region){
371 case BARREL :
372
373 nofblayerhits_barrel = mynofblayerhits ; nofpixelhits_barrel = mynofpixelhits ; nofscthits_barrel = mynofscthits ;
374 noftrthits_barrel = mynoftrthits ; noftrthighthresholdhits_barrel = mynoftrthighthresholdhits;
375 deltaeta1_barrel = mydeltaeta1; deltaphi2_barrel = mydeltaphi2; trackd0_barrel = mytrackd0;
376 eoverp_barrel = myeoverp;
377 et_barrel = myet ; eta_barrel = myeta ; phi_barrel = myphi ;
378 time_barrel = mytime; ehad1_barrel = myehad1; coreem_barrel = myecore;
379 f0_barrel = myf0; f1_barrel = myf1; f2_barrel = myf2; f3_barrel = myf3; re233e237_barrel = myre233e237; re237e277_barrel = myre237e277;
380 fill("MonitorTnP",nofblayerhits_barrel, nofpixelhits_barrel, nofscthits_barrel, noftrthits_barrel, noftrthighthresholdhits_barrel,
381 deltaeta1_barrel, deltaphi2_barrel,trackd0_barrel,eoverp_barrel,
382 et_barrel,eta_barrel,phi_barrel, time_barrel, ehad1_barrel,coreem_barrel, f0_barrel,f1_barrel,f2_barrel, f3_barrel,
383 re233e237_barrel,re237e277_barrel,is_pt_gt_4gev_barrel);
384 break;
385
386 case ENDCAP :
387
388 nofblayerhits_endcap = mynofblayerhits; nofpixelhits_endcap = mynofpixelhits ; nofscthits_endcap = mynofscthits ;
389 noftrthits_endcap = mynoftrthits ; noftrthighthresholdhits_endcap = mynoftrthighthresholdhits;
390 deltaeta1_endcap = mydeltaeta1; deltaphi2_endcap = mydeltaphi2; trackd0_endcap = mytrackd0;
391 eoverp_endcap = myeoverp;
392 et_endcap = myet ; eta_endcap = myeta ; phi_endcap = myphi ;
393 time_endcap = mytime; ehad1_endcap = myehad1; coreem_endcap = myecore;
394 f0_endcap = myf0; f1_endcap = myf1; f2_endcap = myf2; f3_endcap = myf3; re233e237_endcap = myre233e237; re237e277_endcap = myre237e277;
395 fill("MonitorTnP",nofblayerhits_endcap,nofpixelhits_endcap,nofscthits_endcap,noftrthits_endcap,noftrthighthresholdhits_endcap,
396 deltaeta1_endcap, deltaphi2_endcap,trackd0_endcap,eoverp_endcap,et_endcap,eta_endcap,phi_endcap,
397 time_endcap, ehad1_endcap,coreem_endcap,
398 f0_endcap,f1_endcap,f2_endcap,f3_endcap,re233e237_endcap,re237e277_endcap,is_pt_gt_4gev_endcap);
399 break;
400
401 case CRACK :
402 nofblayerhits_crack = mynofblayerhits; nofpixelhits_crack = mynofpixelhits ; nofscthits_crack = mynofscthits ;
403 noftrthits_crack = mynoftrthits ; noftrthighthresholdhits_crack = mynoftrthighthresholdhits;
404 deltaeta1_crack = mydeltaeta1; deltaphi2_crack = mydeltaphi2; trackd0_crack = mytrackd0;
405 eoverp_crack = myeoverp;
406 et_crack = myet ; eta_crack = myeta ; phi_crack = myphi ;
407 time_crack = mytime; ehad1_crack = myehad1; coreem_crack = myecore;
408 f0_crack = myf0; f1_crack = myf1; f2_crack = myf2; f3_crack = myf3; re233e237_crack = myre233e237; re237e277_crack = myre237e277;
409 fill("MonitorTnP",nofblayerhits_crack, nofpixelhits_crack,nofscthits_crack,noftrthits_crack,noftrthighthresholdhits_crack,
410 deltaeta1_crack, deltaphi2_crack,trackd0_crack,eoverp_crack,et_crack,eta_crack,phi_crack,
411 time_crack, ehad1_crack,coreem_crack,f0_crack,f1_crack,f2_crack,f3_crack,re233e237_crack,re237e277_crack,is_pt_gt_4gev_crack);
412 break;
413 default :
414 //ATH_MSG_WARNING("found an electron outside the |eta| > 2.47 acceptance");
415 break;
416 }
417
418 et = myet ; eta = myeta ; phi = myphi ; time = mytime ;
419 ptcone20 = myptcone20; topoetcone40 = mytopoetcone40 ;
420 fill("MonitorTnP",et,eta,phi,time,is_pt_gt_4gev,
421 ptcone20,topoetcone40,lbNProbeCandidates,is_TightID,is_Iso,is_IsoandTight);
422
423 return StatusCode::SUCCESS;
424}
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
float et(const xAOD::jFexSRJetRoI *j)
SG::ReadHandle< xAOD::EventInfo > GetEventInfo(const EventContext &) const
Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
Regions GetRegion(Float_t eta) const
flt_t time() const
Access cluster time.
virtual double eta() const
The pseudorapidity ( ) of the particle.
virtual double e() const
The total energy of the particle.
float energyBE(const unsigned layer) const
Get the energy in one layer of the EM Calo.
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.
time(flags, cells_name, *args, **kw)
@ deltaPhi2
difference between the cluster phi (second sampling) and the phi of the track extrapolated to the sec...
@ deltaEta1
difference between the cluster eta (first sampling) and the eta of the track extrapolated to the firs...
@ e277
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 7x7
Definition EgammaEnums.h:81
@ e237
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x7
Definition EgammaEnums.h:78
@ ecore
core energy in em calo E(core) = E0(3x3) + E1(15x2) + E2(5x5) + E3(3x5)
@ ehad1
E leakage into 1st sampling of had calo (CaloSampling::HEC0 + CaloSampling::TileBar0 + CaloSampling::...
Definition EgammaEnums.h:49
@ e233
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x3 (in cell units e...
Definition EgammaEnums.h:69
@ ptcone20
Track isolation.
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
TrackParticle_v1 TrackParticle
Reference the current persistent version:
@ numberOfTRTHits
number of TRT hits [unit8_t].
@ numberOfSCTHits
number of hits in SCT [unit8_t].
@ numberOfInnermostPixelLayerHits
these are the hits in the 0th pixel barrel layer
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
@ numberOfTRTHighThresholdHits
number of TRT hits which pass the high threshold (only xenon counted) [unit8_t].

◆ fillHistograms()

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

adds event to the monitoring histograms

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

Parameters
ctxforwarded from execute
Returns
StatusCode

Implements AthMonitorAlgorithm.

Definition at line 22 of file MonitorTnPAlgorithm.cxx.

22 {
23 using namespace Monitored;
24
25 // Only monitor good LAr Events :
26
28 if (error_state==xAOD::EventInfo::Error) {
29 ATH_MSG_DEBUG("LAr event data integrity error");
30 return StatusCode::SUCCESS;
31 }
32
33 //
34 // now, fill the specific Electron information
35 //
36 // get the Electron container
37
38 SG::ReadHandle<xAOD::ElectronContainer> electrons(m_ParticleContainerKey, ctx);
39 ATH_CHECK(electrons.isValid());
40
41 // Useful variables to find the leading pair of electrons
42
43 int mynp = 0;
44
45 float lead_et=-999;
46 float subl_et=-999;
47 const xAOD::Electron *lead_el = nullptr;
48 const xAOD::Electron *subl_el = nullptr;
49
50 for (const auto *const e_iter : *electrons) {
51 // Check that the electron meets our requirements
52 bool isGood;
53 if (! e_iter->passSelection(isGood,m_RecoName)) {
54 ATH_MSG_WARNING("Misconfiguration: " << m_RecoName << " is not a valid working point for electrons");
55 break; // no point in continuing
56 }
57 if(isGood) {mynp++;} else continue;
58
59 // Look for two highest pt electrons
60 float etloc = e_iter->pt();
61 if(etloc>lead_et) {
62 subl_et = lead_et;
63 subl_el = lead_el;
64 lead_et = etloc;
65 lead_el = e_iter;
66 } else if(etloc>subl_et) {
67 subl_et = etloc;
68 subl_el = e_iter;
69 }
70 }
71
72 // Check that there are at least a good pair of electrons
73 if (mynp<2) return StatusCode::SUCCESS;
74
75 // Check that this pair meet the T&P selection and is in the correct mass range
76
77 // Et cuts
78 if(lead_et<m_ElectronEtCut || lead_et<0) return StatusCode::SUCCESS;
79 if(subl_et<m_ElectronEtCut || subl_et<0) return StatusCode::SUCCESS;
80
81 // Basic kinematics
82 float lead_eta = lead_el->eta();
83 float lead_phi = lead_el->phi();
84 float subl_eta = subl_el->eta();
85 float subl_phi = subl_el->phi();
86 ATH_MSG_DEBUG("Leading electron (eta,phi,et,q): (" << lead_eta << ", " << lead_phi << ", " << lead_et << ", " << lead_el->charge() << ")");
87 ATH_MSG_DEBUG("Subleading electron (eta,phi,et,q): (" << subl_eta << ", " << subl_phi << ", " << subl_et << ", " << subl_el->charge() << ")");
88
89 // Eta cuts
90 if(fabs(lead_eta)>2.47 || fabs(subl_eta)>2.47) return StatusCode::SUCCESS;
91 // Check charges
92 if (lead_el->charge()*subl_el->charge()>=0) return StatusCode::SUCCESS;
93
94 // Mass window
95 Float_t mass = (lead_el->p4()+subl_el->p4()).M();
96
97 // Monitor the electron pair mass before cutting on the mass range
98
99 auto TnPMass = Monitored::Scalar<Float_t>("MassZ",0.);
100 auto TnPMass_barrel = Monitored::Scalar<Float_t>("MassZ_BARREL",0.);
101 auto TnPMass_endcap = Monitored::Scalar<Float_t>("MassZ_ENDCAP",0.);
102 auto TnPMass_crack = Monitored::Scalar<Float_t>("MassZ_CRACK",0.);
103 TnPMass = mass;
104 fill("MonitorTnP",TnPMass);
105
106 auto region = GetRegion(lead_eta);
107 switch(region){
108 case BARREL :
109 TnPMass_barrel = mass;
110 fill("MonitorTnP", TnPMass_barrel);
111 break;
112 case ENDCAP :
113 TnPMass_endcap = mass;
114 fill("MonitorTnP", TnPMass_endcap);
115 break;
116 case CRACK :
117 TnPMass_crack = mass;
118 fill("MonitorTnP", TnPMass_crack);
119 break;
120 default :
121 //ATH_MSG_WARNING("found an electron outside the |eta| > 2.47 acceptance");
122 break;
123 }
124
125 ATH_MSG_DEBUG("ee mass & cuts: (" << mass << ", " << m_MassLowerCut << ", " << m_MassUpperCut << ")");
126 if(mass<m_MassLowerCut || mass>m_MassUpperCut) return StatusCode::SUCCESS;
127
128 // Register this good Tag & Probe candidate
129 // count the number of good Z candidate per LB
130
131 auto lbNZ = Monitored::Scalar<u_int16_t>("LB",0);
132 lbNZ = GetEventInfo(ctx)->lumiBlock();
133 fill("MonitorTnP",lbNZ);
134
135 bool lead_isLHTight = false;
136 if ( !lead_el->passSelection(lead_isLHTight,"LHTight") ) return StatusCode::FAILURE;
137
138 bool subl_isLHTight = false;
139 if ( !subl_el->passSelection(subl_isLHTight,"LHTight") ) return StatusCode::FAILURE;
140
141 // If leading electron is LHTight use subleading as probe
142 if(lead_isLHTight) {
143 ATH_CHECK(fillElectronProbe(subl_el, true, ctx));
144 }
145 // If subleading electron is LHTight use leading as probe
146 if(subl_isLHTight) {
147 ATH_CHECK(fillElectronProbe(lead_el, false, ctx));
148 }
149
150return StatusCode::SUCCESS;
151
152}
#define ATH_CHECK
Evaluate an expression and check for errors.
StatusCode fillElectronProbe(const xAOD::Electron *el, const bool isleading, const EventContext &ctx) const
SG::ReadHandleKey< xAOD::ElectronContainer > m_ParticleContainerKey
Gaudi::Property< Float_t > m_MassUpperCut
Gaudi::Property< Float_t > m_MassLowerCut
Gaudi::Property< std::string > m_RecoName
Gaudi::Property< Float_t > m_ElectronEtCut
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition Egamma_v1.cxx:71
bool passSelection(bool &value, const std::string &menu) const
Check if the egamma object pass a selection menu (using the name) If the menu decision is stored in t...
virtual FourMom_t p4() const override final
The full 4-momentum of the particle as a TLoretzVector.
Definition Egamma_v1.cxx:94
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition Egamma_v1.cxx:76
float charge() const
Obtain the charge of the object.
@ LAr
The LAr calorimeter.
EventFlagErrorState
States that a given sub-detector could be in.
@ Error
The sub-detector issued an error.
Electron_v1 Electron
Definition of the current "egamma version".

◆ filterPassed()

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

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

◆ GetEventInfo()

SG::ReadHandle< xAOD::EventInfo > AthMonitorAlgorithm::GetEventInfo ( const EventContext & ctx) const
inherited

Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)

Parameters
ctxEventContext for the event
Returns
a SG::ReadHandle<xAOD::EventInfo>

Definition at line 111 of file AthMonitorAlgorithm.cxx.

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

◆ getGroup()

const ToolHandle< GenericMonitoringTool > & AthMonitorAlgorithm::getGroup ( const std::string & name) const
inherited

Get a specific monitoring tool from the tool handle array.

Finds a specific GenericMonitoringTool instance from the list of monitoring tools (a ToolHandleArray). Throws a FATAL warning if the object found is empty.

Parameters
namestring name of the desired tool
Returns
reference to the desired monitoring tool

Definition at line 168 of file AthMonitorAlgorithm.cxx.

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

◆ GetRegion()

Regions MonitorTnPAlgorithm::GetRegion ( Float_t eta) const
inlineprivate

Definition at line 45 of file MonitorTnPAlgorithm.h.

45 {
46 Float_t aeta = fabs(eta);
47 // check if object is in barrel
48 if ( aeta < 1.37 ) return MonitorTnPAlgorithm::BARREL;
49 // check if object is in end-caps
50 if ( aeta > 1.52 && aeta < 2.47 ) return MonitorTnPAlgorithm::ENDCAP;
51 // check if object is in crack region
52 if ( aeta > 1.37 && aeta < 1.52 ) return MonitorTnPAlgorithm::CRACK;
53 return MonitorTnPAlgorithm::NREGIONS; // out of acceptance
54 }

◆ getTrigDecisionTool()

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

Get the trigger decision tool member.

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

Returns
m_trigDecTool

Definition at line 198 of file AthMonitorAlgorithm.cxx.

198 {
199 return m_trigDecTool;
200}

◆ initialize()

StatusCode MonitorTnPAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 14 of file MonitorTnPAlgorithm.cxx.

14 {
15 using namespace Monitored;
19 return StatusCode::SUCCESS;
20}
virtual StatusCode initialize() override
initialize
SG::ReadDecorHandleKeyArray< xAOD::ElectronContainer > m_ElectronIsolationKey

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

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ parseList()

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

Parse a string into a vector.

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

Parameters
lineThe input string.
resultThe parsed output vector of strings.
Returns
StatusCode

Definition at line 345 of file AthMonitorAlgorithm.cxx.

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

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

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

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

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

◆ sysExecute()

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

Execute an algorithm.

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

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

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

◆ sysInitialize()

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

Override sysInitialize.

Override sysInitialize from the base class.

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

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

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

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

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

◆ sysStart()

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

Handle START transition.

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

◆ trigChainsArePassed()

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

Check whether triggers are passed.

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

Parameters
vTrigNamesList of trigger names.
Returns
If empty input, default to true. If at least one trigger is specified, returns whether at least one trigger was passed.

Definition at line 203 of file AthMonitorAlgorithm.cxx.

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

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_dataType

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::m_dataType
protectedinherited

Instance of the DataType_t enum.

Definition at line 356 of file AthMonitorAlgorithm.h.

◆ m_dataTypeStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_dataTypeStr {this,"DataType","userDefined"}
protectedinherited

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

Definition at line 358 of file AthMonitorAlgorithm.h.

358{this,"DataType","userDefined"};

◆ m_defaultLBDuration

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

Default duration of one lumi block.

Definition at line 365 of file AthMonitorAlgorithm.h.

365{this,"DefaultLBDuration",60.};

◆ m_detailLevel

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

Sets the level of detail used in the monitoring.

Definition at line 366 of file AthMonitorAlgorithm.h.

366{this,"DetailLevel",0};

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_DQFilterTools

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

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

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

◆ m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy
privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

◆ m_ElectronEtCut

Gaudi::Property<Float_t> MonitorTnPAlgorithm::m_ElectronEtCut {this, "ElectronEtCut", 15000, "Et cut for electron"}
private

Definition at line 41 of file MonitorTnPAlgorithm.h.

41{this, "ElectronEtCut", 15000, "Et cut for electron"};

◆ m_ElectronIsolationKey

SG::ReadDecorHandleKeyArray<xAOD::ElectronContainer> MonitorTnPAlgorithm::m_ElectronIsolationKey {this, "ElectronIsolationKey", {"Electrons.ptcone20", "Electrons.topoetcone40"} }
private

Definition at line 35 of file MonitorTnPAlgorithm.h.

35{this, "ElectronIsolationKey", {"Electrons.ptcone20", "Electrons.topoetcone40"} };

◆ m_enforceExpressTriggers

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

Definition at line 377 of file AthMonitorAlgorithm.h.

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

◆ m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment
protectedinherited

Instance of the Environment_t enum.

Definition at line 355 of file AthMonitorAlgorithm.h.

◆ m_environmentStr

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

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

Definition at line 357 of file AthMonitorAlgorithm.h.

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

◆ m_EventInfoKey

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

Key for retrieving EventInfo from StoreGate.

Definition at line 367 of file AthMonitorAlgorithm.h.

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

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

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

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_fileKey

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

Internal Athena name for file.

Definition at line 363 of file AthMonitorAlgorithm.h.

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

◆ m_lbDurationDataKey

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

Definition at line 350 of file AthMonitorAlgorithm.h.

351{this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"};

◆ m_lumiDataKey

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

Definition at line 348 of file AthMonitorAlgorithm.h.

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

◆ m_MassLowerCut

Gaudi::Property<Float_t> MonitorTnPAlgorithm::m_MassLowerCut {this, "MassLowerCut", 70000, "Lower Mass cut"}
private

Definition at line 42 of file MonitorTnPAlgorithm.h.

42{this, "MassLowerCut", 70000, "Lower Mass cut"};

◆ m_MassPeak

Gaudi::Property<Float_t> MonitorTnPAlgorithm::m_MassPeak {this, "MassPeak", ParticleConstants::ZMassInMeV, "Resonance peak position"}
private

Definition at line 40 of file MonitorTnPAlgorithm.h.

40{this, "MassPeak", ParticleConstants::ZMassInMeV, "Resonance peak position"};
constexpr double ZMassInMeV
the mass of the Z0 boson (in MeV)

◆ m_MassUpperCut

Gaudi::Property<Float_t> MonitorTnPAlgorithm::m_MassUpperCut {this, "MassUpperCut", 110000, "Upper Mass Cut"}
private

Definition at line 43 of file MonitorTnPAlgorithm.h.

43{this, "MassUpperCut", 110000, "Upper Mass Cut"};

◆ m_name

std::string AthMonitorAlgorithm::m_name
privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

◆ m_ParticleContainerKey

SG::ReadHandleKey<xAOD::ElectronContainer> MonitorTnPAlgorithm::m_ParticleContainerKey {this, "ParticleContainerName", "Electrons", "Name of electron container" }
private

Definition at line 34 of file MonitorTnPAlgorithm.h.

34{this, "ParticleContainerName", "Electrons", "Name of electron container" };

◆ m_ParticlePrefix

Gaudi::Property<std::string> MonitorTnPAlgorithm::m_ParticlePrefix {this,"ParticlePrefix","electron","Name of the particle prefix to be used to define hists"}
private

Definition at line 36 of file MonitorTnPAlgorithm.h.

36{this,"ParticlePrefix","electron","Name of the particle prefix to be used to define hists"};

◆ m_RecoName

Gaudi::Property<std::string> MonitorTnPAlgorithm::m_RecoName {this,"RecoName","LHLoose","Name of particle flavor in egamma reco"}
private

Definition at line 37 of file MonitorTnPAlgorithm.h.

37{this,"RecoName","LHLoose","Name of particle flavor in egamma reco"};

◆ m_TnPType

Gaudi::Property<std::string> MonitorTnPAlgorithm::m_TnPType {this,"TnPType","Z","Variable to describe what is the TnP resonance name"}
private

Definition at line 38 of file MonitorTnPAlgorithm.h.

38{this,"TnPType","Z","Variable to describe what is the TnP resonance name"};

◆ m_toolLookupMap

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

Definition at line 372 of file AthMonitorAlgorithm.h.

◆ m_tools

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

Array of Generic Monitoring Tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

341{this,"GMTools",{}};

◆ m_trigDecTool

PublicToolHandle<Trig::TrigDecisionTool> AthMonitorAlgorithm::m_trigDecTool
protectedinherited

Tool to tell whether a specific trigger is passed.

Definition at line 345 of file AthMonitorAlgorithm.h.

◆ m_triggerChainString

Gaudi::Property<std::string> AthMonitorAlgorithm::m_triggerChainString {this,"TriggerChain",""}
protectedinherited

Trigger chain string pulled from the job option and parsed into a vector.

Definition at line 360 of file AthMonitorAlgorithm.h.

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

◆ m_trigLiveFractionDataKey

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

Definition at line 352 of file AthMonitorAlgorithm.h.

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

◆ m_useLumi

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

Allows use of various luminosity functions.

Definition at line 364 of file AthMonitorAlgorithm.h.

364{this,"EnableLumi",false};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_vTrigChainNames

std::vector<std::string> AthMonitorAlgorithm::m_vTrigChainNames
protectedinherited

Vector of trigger chain names parsed from trigger chain string.

Definition at line 361 of file AthMonitorAlgorithm.h.


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