ATLAS Offline Software
Loading...
Searching...
No Matches
HLT::MET::TrkMHTFex Class Reference

#include <TrkMHTFex.h>

Inheritance diagram for HLT::MET::TrkMHTFex:
Collaboration diagram for HLT::MET::TrkMHTFex:

Public Member Functions

 TrkMHTFex (const std::string &name, ISvcLocator *pSvcLocator)
 Constructor.
virtual StatusCode initialize () override
 Initialise the fex.
virtual StatusCode execute (const EventContext &context) const override
 Run the algorithm.
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

StatusCode initializeBase (const std::vector< std::string > &componentNames)
 Initialize the base class.
virtual StatusCode monitor (const xAOD::TrigMissingET &met, MonGroupBuilder &monitors) const
 Add monitor variables from an xAOD::TrigMissingET object.
virtual StatusCode flagMET (xAOD::TrigMissingET &met) const
 Flag suspicious values in the output MET.
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

std::vector< std::string > m_componentNames
 The names of the output MET components.

Private Types

template<typename T>
using deferred_t = std::optional<T>
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

virtual StatusCode fillMET (xAOD::TrigMissingET &met, const EventContext &context, MonGroupBuilder &monitors) const override
 Calculate and fill the output MET value.
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::JetContainerm_jetKey
 Input jets and jvt.
SG::ReadHandleKey< xAOD::TrackParticleContainerm_trackKey
 Input tracks.
SG::ReadHandleKey< xAOD::VertexContainerm_vertexKey
 Input vertices.
SG::ReadHandleKey< jet::TrackVertexAssociationm_tvaKey
 The track-vertex assocation name.
SG::ReadDecorHandleKey< xAOD::JetContainerm_jvtKey
 The Jvt selection name.
Gaudi::Property< std::string > m_trackGAName
 The ghost-association aux element name.
Gaudi::Property< float > m_forwardJetPt
 Pt selection on forward jets.
Gaudi::Property< float > m_minJvtJetPt
 Minimum pt selection for JVT on central jets.
Gaudi::Property< float > m_maxJvtJetPt
 Maximum pt selection for JVT on central jets.
Gaudi::Property< float > m_jvtSelection
 The JVT selection.
ToolHandle< InDet::IInDetTrackSelectionToolm_trackSelTool
 The track selection tool.
Gaudi::Property< float > m_tstPtCeil
 The maximum pT for tracks in the soft term.
deferred_t< SG::AuxElement::ConstAccessor< std::vector< ElementLink< xAOD::IParticleContainer > > > > m_trackGA
SG::WriteHandleKey< xAOD::TrigMissingETContainerm_metContainerKey
 The output MET object.
ToolHandle< GenericMonitoringToolm_monTool
 The monitoring tool.
Gaudi::Property< float > m_maxComponentMETSumEtRatio
Gaudi::Property< float > m_maxGlobalMETSumEtRatio
bool m_baseInitialised {false}
 Whether or not this class has been correctly initialised.
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 36 of file TrkMHTFex.h.

Member Typedef Documentation

◆ deferred_t

template<typename T>
using HLT::MET::TrkMHTFex::deferred_t = std::optional<T>
private

Definition at line 40 of file TrkMHTFex.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TrkMHTFex()

HLT::MET::TrkMHTFex::TrkMHTFex ( const std::string & name,
ISvcLocator * pSvcLocator )

Constructor.

Definition at line 29 of file TrkMHTFex.cxx.

29 :
30 FexBase(name, pSvcLocator)
31 {}
FexBase(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.
Definition FexBase.cxx:43

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}

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

◆ 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 HLT::MET::FexBase::execute ( const EventContext & context) const
overridevirtualinherited

Run the algorithm.

Definition at line 60 of file FexBase.cxx.

61 {
62 if (!m_baseInitialised) {
63 ATH_MSG_ERROR("Base class was not initialised! This means that the "
64 << "derived class was not correctly written!");
65 return StatusCode::FAILURE;
66 }
67 ATH_MSG_DEBUG("Executing " << name() << "...");
68 Monitored::Timer totalTimer("TIME_Total");
69 // Create the output
70 auto metCont = std::make_unique<xAOD::TrigMissingETContainer>();
71 auto metContAux = std::make_unique<xAOD::TrigMissingETAuxContainer>();
72 metCont->setStore(metContAux.get() );
73
74 // Create the actual output object
75 metCont->push_back(std::make_unique<xAOD::TrigMissingET>());
76 xAOD::TrigMissingET* met = metCont->back();
77
78 // Initialise the components
80 // We also need to initialise all of the values. This is not done by the EDM
81 // class so you can easily get missing aux element errors if you don't do
82 // this
83 met->setEx(0);
84 met->setEy(0);
85 met->setEz(0);
86 met->setSumEt(0);
87 met->setSumE(0);
88 met->setFlag(0);
89 met->setRoiWord(0);
90
91 MonGroupBuilder monitors;
92 CHECK( fillMET(*met, context, monitors) );
93 // Add flags
94 CHECK( flagMET(*met) );
95 // Add extra monitors
96 CHECK( monitor(*met, monitors) );
97 // Create the actual group and trigger the monitoring
98 monitors.build(m_monTool);
99
100 // Output REGTEST information
101 if (msgLvl(MSG::DEBUG) ) {
102 ATH_MSG_DEBUG( "REGTEST " << METComponent(*met) );
103 ATH_MSG_DEBUG( "REGTEST flag = " << met->flag() );
104 ATH_MSG_DEBUG( "REGTEST Name, status, values: ");
105 for (std::size_t ii = 0; ii < met->getNumberOfComponents(); ++ii)
106 ATH_MSG_DEBUG( "REGTEST "
107 << met->nameOfComponent(ii) << ", "
108 << met->statusComponent(ii) << ", "
109 << METComponent(ii, *met) );
110 }
111
112 // Push this output to the store
113 auto handle = SG::makeHandle(m_metContainerKey, context);
114 CHECK( handle.record( std::move(metCont), std::move(metContAux) ) );
115 return StatusCode::SUCCESS;
116 }
#define ATH_MSG_ERROR(x)
#define ATH_MSG_DEBUG(x)
#define CHECK(...)
Evaluate an expression and check for errors.
bool msgLvl(const MSG::Level lvl) const
virtual StatusCode flagMET(xAOD::TrigMissingET &met) const
Flag suspicious values in the output MET.
Definition FexBase.cxx:118
std::vector< std::string > m_componentNames
The names of the output MET components.
Definition FexBase.h:48
bool m_baseInitialised
Whether or not this class has been correctly initialised.
Definition FexBase.h:119
ToolHandle< GenericMonitoringTool > m_monTool
The monitoring tool.
Definition FexBase.h:106
virtual StatusCode fillMET(xAOD::TrigMissingET &met, const EventContext &context, MonGroupBuilder &monitors) const =0
Calculate and fill the output MET value.
SG::WriteHandleKey< xAOD::TrigMissingETContainer > m_metContainerKey
The output MET object.
Definition FexBase.h:103
virtual StatusCode monitor(const xAOD::TrigMissingET &met, MonGroupBuilder &monitors) const
Add monitor variables from an xAOD::TrigMissingET object.
Definition FexBase.cxx:141
int flag() const
Get the status flag.
void setEy(float s)
Set the y component of the missing energy.
void setRoiWord(uint32_t s)
Set the ROI word.
void setEx(float s)
Set the x component of the missing energy.
const std::string & nameOfComponent(unsigned int index) const
Get the name of a component.
void setSumE(float s)
Set the sum of the ET of the missing energy.
void setSumEt(float s)
Set the sum of the ET of the missing energy.
unsigned int getNumberOfComponents() const
Get the number of components.
short statusComponent(unsigned int index) const
Get Status.
void defineComponents(const std::vector< std::string > &componentNames)
Define each of the components by name.
void setFlag(int s)
Set the status flag.
void setEz(float s)
Set the z component of the missing energy.
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
TrigMissingET_v1 TrigMissingET
Define the most recent version of the TrigMissingET class.

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

◆ fillMET()

StatusCode HLT::MET::TrkMHTFex::fillMET ( xAOD::TrigMissingET & met,
const EventContext & context,
MonGroupBuilder & monitors ) const
overrideprivatevirtual

Calculate and fill the output MET value.

Parameters
metThe object to fill
contextThe event context
monitors[out]Any extra monitors to fill

Implements HLT::MET::FexBase.

Definition at line 56 of file TrkMHTFex.cxx.

60 {
61 // Retrieve the inputs
62 auto jets = SG::makeHandle(m_jetKey, context);
63 if (!jets.isValid())
64 {
65 ATH_MSG_ERROR("Failed to retrieve " << m_jetKey);
66 return StatusCode::FAILURE;
67 }
68 auto tracks = SG::makeHandle(m_trackKey, context);
69 if (!tracks.isValid())
70 {
71 ATH_MSG_ERROR("Failed to retrieve " << m_trackKey);
72 return StatusCode::FAILURE;
73 }
74 auto vertices = SG::makeHandle(m_vertexKey, context);
75 if (!vertices.isValid())
76 {
77 ATH_MSG_ERROR("Failed to retrieve " << m_vertexKey);
78 return StatusCode::FAILURE;
79 }
80 auto tva = SG::makeHandle(m_tvaKey, context);
81 auto jvtAcc = SG::makeHandle<float>(m_jvtKey, context);
82 auto& trackLinksAcc = *m_trackGA;
83
84 // Work out which is the primary vertex
85 // I do not know if there is a way to see which vertex was counted as
86 // primary by the JVT calculation but I see no reason to suspect that it
87 // would be anything other than the one found by this method
88 const xAOD::Vertex* priVtx = nullptr;
89 for (const xAOD::Vertex* vtx : *vertices)
90 if (vtx->vertexType() == xAOD::VxType::PriVtx) {
91 priVtx = vtx;
92 break;
93 }
94
95 // Prepare the output values
96 std::array<METComponent, 4> mhtSums;
97 METComponent tstSum;
98
99 // If no primary vertex is identified, the algorithm will return 0
100 if (priVtx) {
101 // Prepare a list of primary vertex tracks
102 std::vector<const xAOD::TrackParticle*> tstTracks;
103 tstTracks.reserve(tracks->size() );
104 for (const xAOD::TrackParticle* itrk : *tracks)
105 if (tva->associatedVertex(itrk) == priVtx &&
106 itrk->pt() < m_tstPtCeil &&
107 m_trackSelTool->accept(*itrk, priVtx) )
108 tstTracks.push_back(itrk);
109 // Keep indices of tracks associated to hard-scatter jets
110 std::unordered_set<std::size_t> jetTrackIndices;
111 // Iterate over the jets
112 for (const xAOD::Jet* ijet : *jets) {
113 if (isCentral(ijet) ) {
114 if (ijet->pt() < m_minJvtJetPt)
115 // Skip central jets below the JVT pT threshold
116 continue;
117 if (ijet->pt() < m_maxJvtJetPt && jvtAcc(*ijet) < m_jvtSelection)
118 // Skip central jets below the maximum JVT pT threshold and with JVT
119 // below threshold
120 continue;
121 } // end if isCentral(ijet)
122 else {
123 if (ijet->pt() < m_forwardJetPt)
124 // Skip forward jets below the forward jet pT threshold
125 continue;
126 }
127 // The index of the mht METComponent
128 std::size_t componentIdx = 0;
129 if (ijet->eta() < 0)
130 ++componentIdx;
131 if (!isCentral(ijet) )
132 componentIdx += 2;
133 // Add the jet into the correct component
134 mhtSums.at(componentIdx) += ijet->p4();
135 // Iterate over the ghost-associated tracks
136 for (const auto& link : trackLinksAcc(*ijet) ) {
137 if (link.getDataPtr() != tracks.ptr() ) {
138 ATH_MSG_ERROR("Tracks linked with key '" << m_trackGAName
139 << "' do not match tracks retrieved with '"
140 << m_trackKey << "'!");
141 return StatusCode::FAILURE;
142 }
143 jetTrackIndices.insert(link.index() );
144 }
145 }
146 // Build up the track soft term
147 for (const xAOD::TrackParticle* itrk : tstTracks)
148 if (jetTrackIndices.count(itrk->index() ) == 0)
149 // Only include tracks that weren't associated to a hard-scatter jet
150 tstSum += itrk->p4();
151 } //> end if priVtx
152 // Save the sum over components
153 std::accumulate(mhtSums.begin(), mhtSums.end(), tstSum).fillMET(met);
154 // Save each component
155 for (std::size_t ii = 0; ii < 4; ++ii)
156 mhtSums.at(ii).fillMETComponent(ii, met);
157 // Save the tst
158 tstSum.fillMETComponent(4, met);
159
160 // Debugging information
161 ATH_MSG_DEBUG("Jets: "
162 << std::accumulate(mhtSums.begin(), mhtSums.end(), METComponent{}) );
163
164 return StatusCode::SUCCESS;
165 }
ToolHandle< InDet::IInDetTrackSelectionTool > m_trackSelTool
The track selection tool.
Definition TrkMHTFex.h:87
Gaudi::Property< float > m_tstPtCeil
The maximum pT for tracks in the soft term.
Definition TrkMHTFex.h:92
SG::ReadDecorHandleKey< xAOD::JetContainer > m_jvtKey
The Jvt selection name.
Definition TrkMHTFex.h:67
SG::ReadHandleKey< jet::TrackVertexAssociation > m_tvaKey
The track-vertex assocation name.
Definition TrkMHTFex.h:63
SG::ReadHandleKey< xAOD::JetContainer > m_jetKey
Input jets and jvt.
Definition TrkMHTFex.h:53
deferred_t< SG::AuxElement::ConstAccessor< std::vector< ElementLink< xAOD::IParticleContainer > > > > m_trackGA
Definition TrkMHTFex.h:113
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackKey
Input tracks.
Definition TrkMHTFex.h:56
SG::ReadHandleKey< xAOD::VertexContainer > m_vertexKey
Input vertices.
Definition TrkMHTFex.h:60
Gaudi::Property< std::string > m_trackGAName
The ghost-association aux element name.
Definition TrkMHTFex.h:70
Gaudi::Property< float > m_minJvtJetPt
Minimum pt selection for JVT on central jets.
Definition TrkMHTFex.h:76
Gaudi::Property< float > m_forwardJetPt
Pt selection on forward jets.
Definition TrkMHTFex.h:73
Gaudi::Property< float > m_maxJvtJetPt
Maximum pt selection for JVT on central jets.
Definition TrkMHTFex.h:80
Gaudi::Property< float > m_jvtSelection
The JVT selection.
Definition TrkMHTFex.h:84
@ PriVtx
Primary vertex.
Jet_v1 Jet
Definition of the current "jet version".
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Vertex_v1 Vertex
Define the latest version of the vertex class.

◆ filterPassed()

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

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

◆ flagMET()

StatusCode HLT::MET::FexBase::flagMET ( xAOD::TrigMissingET & met) const
protectedvirtualinherited

Flag suspicious values in the output MET.

Parameters
metThe xAOD::TrigMissingET object to flag

Definition at line 118 of file FexBase.cxx.

119 {
120 // Start with the components
121 // Keep a flag to OR into the main value
122 int overall = 0;
123 for (std::size_t ii = 0; ii < met.getNumberOfComponents(); ++ii) {
124 METComponent component(ii, met);
125 if (component.sumEt > 0 &&
126 component.met() / component.sumEt > m_maxComponentMETSumEtRatio) {
127 component.status |= StatusFlag::ComponentBigMEtSEtRatio;
129 }
130 met.setStatusComponent(ii, component.status);
131 }
132 METComponent total(met);
133 total.status |= overall;
134 if (total.sumEt > 0 &&
135 total.met() / total.sumEt > m_maxGlobalMETSumEtRatio)
136 total.status |= StatusFlag::GlobalBigMEtSEtRatio;
137 met.setFlag(total.status);
138 return StatusCode::SUCCESS;
139 }
Gaudi::Property< float > m_maxGlobalMETSumEtRatio
Definition FexBase.h:111
Gaudi::Property< float > m_maxComponentMETSumEtRatio
Definition FexBase.h:108
void setStatusComponent(unsigned int index, float value)
Set Status.

◆ initialize()

StatusCode HLT::MET::TrkMHTFex::initialize ( )
overridevirtual

Initialise the fex.

Definition at line 33 of file TrkMHTFex.cxx.

34 {
35 CHECK( m_jetKey.initialize() );
36 CHECK( m_trackKey.initialize() );
37 CHECK( m_vertexKey.initialize() );
38 CHECK( m_tvaKey.initialize() );
39 CHECK( m_trackSelTool.retrieve() );
40
41 // Update the decor keys if necessary
42 if (m_jvtKey.key().find(".") == std::string::npos)
43 m_jvtKey = m_jetKey.key() + "." + m_jvtKey.key();
44 else if (SG::contKeyFromKey(m_jvtKey.key()) != m_jetKey.key())
45 {
46 ATH_MSG_ERROR("JVT container key does not match jet key!");
47 return StatusCode::FAILURE;
48 }
49 CHECK(m_jvtKey.initialize());
50 m_trackGA.emplace(m_trackGAName);
51
52 return initializeBase(
53 {"JETB1", "JETB2", "JETE1", "JETE2", "TrackSoftTerm"});
54 }
StatusCode initializeBase(const std::vector< std::string > &componentNames)
Initialize the base class.
Definition FexBase.cxx:48
std::string contKeyFromKey(const std::string &key)
Extract the container part of key.

◆ initializeBase()

StatusCode HLT::MET::FexBase::initializeBase ( const std::vector< std::string > & componentNames)
protectedinherited

Initialize the base class.

Parameters
componentNamesThe names of any components in the output MET object.

The xAOD::TrigMissingET writes out a main 'met' value but can also write out extra MET components - either for debugging, monitoring or further use offline. The object has to be initialized with this names, they cannot be extended after this, so this information must be available for the base class to use.

Definition at line 48 of file FexBase.cxx.

50 {
51 ATH_MSG_DEBUG("Initialising FexBase base class");
52 m_baseInitialised = true;
53 m_componentNames = componentNames;
54 CHECK( m_metContainerKey.initialize() );
55 if (!m_monTool.empty())
56 CHECK( m_monTool.retrieve() );
57 return StatusCode::SUCCESS;
58 }

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

◆ monitor()

StatusCode HLT::MET::FexBase::monitor ( const xAOD::TrigMissingET & met,
MonGroupBuilder & monitors ) const
protectedvirtualinherited

Add monitor variables from an xAOD::TrigMissingET object.

Parameters
metThe object to fill from
monitorsExtra variables to be monitored

Definition at line 141 of file FexBase.cxx.

144 {
145 METComponent metComponent(met);
146 // Increase the capacity of the vector
147 monitors.increaseCapacity(16, true);
148 // Add the standard variables
149 monitors.add(Monitored::Scalar("EF_MEx_log", toLog10GeV(metComponent.mpx) ) );
150 monitors.add(Monitored::Scalar("EF_MEy_log", toLog10GeV(metComponent.mpy) ) );
151 monitors.add(Monitored::Scalar("EF_MEz_log", toLog10GeV(metComponent.mpz) ) );
152 monitors.add(Monitored::Scalar("EF_MET_log", toLog10GeV(metComponent.met() ) ) );
153 monitors.add(Monitored::Scalar("EF_ME_log", toLog10GeV(metComponent.magnitude() ) ) );
154 monitors.add(Monitored::Scalar("EF_SumEt_log", toLog10GeV(metComponent.sumEt) ) );
155 monitors.add(Monitored::Scalar("EF_SumE_log", toLog10GeV(metComponent.sumE, -9e9) ) );
156 monitors.add(Monitored::Scalar("EF_MEx_lin", toLinGeV(metComponent.mpx) ) );
157 monitors.add(Monitored::Scalar("EF_MEy_lin", toLinGeV(metComponent.mpy) ) );
158 monitors.add(Monitored::Scalar("EF_MEz_lin", toLinGeV(metComponent.mpz) ) );
159 monitors.add(Monitored::Scalar("EF_MET_lin", toLinGeV(metComponent.met() ) ) );
160 monitors.add(Monitored::Scalar("EF_ME_lin", toLinGeV(metComponent.magnitude() ) ) );
161 monitors.add(Monitored::Scalar("EF_SumEt_lin", toLinGeV(metComponent.sumEt) ) );
162 monitors.add(Monitored::Scalar("EF_SumE_lin", toLinGeV(metComponent.sumE) ) );
163 monitors.add(Monitored::Scalar(
164 "EF_XS", toLinGeV(metComponent.met() ) / toLinGeV(metComponent.sumEt, 1) ) );
165 monitors.add(Monitored::Scalar("EF_MET_phi", metComponent.phi() ) );
166 return StatusCode::SUCCESS;
167 }

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

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

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

bool HLT::MET::FexBase::m_baseInitialised {false}
privateinherited

Whether or not this class has been correctly initialised.

Definition at line 119 of file FexBase.h.

119{false};

◆ m_componentNames

std::vector<std::string> HLT::MET::FexBase::m_componentNames
protectedinherited

The names of the output MET components.

Definition at line 48 of file FexBase.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_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_forwardJetPt

Gaudi::Property<float> HLT::MET::TrkMHTFex::m_forwardJetPt
private
Initial value:
{
this, "ForwardJetPt", 0., "The pT cut applied to forward jets"}

Pt selection on forward jets.

Definition at line 73 of file TrkMHTFex.h.

73 {
74 this, "ForwardJetPt", 0., "The pT cut applied to forward jets"};

◆ m_jetKey

SG::ReadHandleKey<xAOD::JetContainer> HLT::MET::TrkMHTFex::m_jetKey
private
Initial value:
{
this, "JetName", "jets", "Input jet collection"}

Input jets and jvt.

Definition at line 53 of file TrkMHTFex.h.

53 {
54 this, "JetName", "jets", "Input jet collection"};

◆ m_jvtKey

SG::ReadDecorHandleKey<xAOD::JetContainer> HLT::MET::TrkMHTFex::m_jvtKey
private
Initial value:
{
this, "JvtName", "Jvt", "The name of the JVT decorator"}

The Jvt selection name.

Definition at line 67 of file TrkMHTFex.h.

67 {
68 this, "JvtName", "Jvt", "The name of the JVT decorator"};

◆ m_jvtSelection

Gaudi::Property<float> HLT::MET::TrkMHTFex::m_jvtSelection
private
Initial value:
{
this, "JVTCut", 0.9, "The JVT selection"}

The JVT selection.

Definition at line 84 of file TrkMHTFex.h.

84 {
85 this, "JVTCut", 0.9, "The JVT selection"};

◆ m_maxComponentMETSumEtRatio

Gaudi::Property<float> HLT::MET::FexBase::m_maxComponentMETSumEtRatio
privateinherited
Initial value:
{
"MaxComponentMETSumEtRatio", 1.,
"The maximum MET/SumEt ratio per component"}

Definition at line 108 of file FexBase.h.

108 {
109 "MaxComponentMETSumEtRatio", 1.,
110 "The maximum MET/SumEt ratio per component"};

◆ m_maxGlobalMETSumEtRatio

Gaudi::Property<float> HLT::MET::FexBase::m_maxGlobalMETSumEtRatio
privateinherited
Initial value:
{
"MaxGlobalMETSumEtRatio", 1.,
"The maximum MET/SumEt ratio for the total value"}

Definition at line 111 of file FexBase.h.

111 {
112 "MaxGlobalMETSumEtRatio", 1.,
113 "The maximum MET/SumEt ratio for the total value"};

◆ m_maxJvtJetPt

Gaudi::Property<float> HLT::MET::TrkMHTFex::m_maxJvtJetPt
private
Initial value:
{
this, "MaximumJvtJetPt", 50. * Gaudi::Units::GeV,
"The maximum pT for central jets for JVT selections to be used"}

Maximum pt selection for JVT on central jets.

Definition at line 80 of file TrkMHTFex.h.

80 {
81 this, "MaximumJvtJetPt", 50. * Gaudi::Units::GeV,
82 "The maximum pT for central jets for JVT selections to be used"};

◆ m_metContainerKey

SG::WriteHandleKey<xAOD::TrigMissingETContainer> HLT::MET::FexBase::m_metContainerKey
privateinherited
Initial value:
{
this, "METContainerKey", "HLT_MET", "Name of the output MET object"}

The output MET object.

Definition at line 103 of file FexBase.h.

103 {
104 this, "METContainerKey", "HLT_MET", "Name of the output MET object"};

◆ m_minJvtJetPt

Gaudi::Property<float> HLT::MET::TrkMHTFex::m_minJvtJetPt
private
Initial value:
{
this, "MinimumJvtJetPt", 20. * Gaudi::Units::GeV,
"The minimum pT for central jets for JVT selections to be used"}

Minimum pt selection for JVT on central jets.

Definition at line 76 of file TrkMHTFex.h.

76 {
77 this, "MinimumJvtJetPt", 20. * Gaudi::Units::GeV,
78 "The minimum pT for central jets for JVT selections to be used"};

◆ m_monTool

ToolHandle<GenericMonitoringTool> HLT::MET::FexBase::m_monTool
privateinherited
Initial value:
{
this, "MonTool", "", "Monitoring tool"}

The monitoring tool.

Definition at line 106 of file FexBase.h.

106 {
107 this, "MonTool", "", "Monitoring tool"};

◆ m_trackGA

deferred_t<SG::AuxElement::ConstAccessor<std::vector<ElementLink<xAOD::IParticleContainer> > > > HLT::MET::TrkMHTFex::m_trackGA
private

Definition at line 113 of file TrkMHTFex.h.

◆ m_trackGAName

Gaudi::Property<std::string> HLT::MET::TrkMHTFex::m_trackGAName
private
Initial value:
{
this, "TrackLinkName", "tracks", "The name of the jet track links"}

The ghost-association aux element name.

Definition at line 70 of file TrkMHTFex.h.

70 {
71 this, "TrackLinkName", "tracks", "The name of the jet track links"};

◆ m_trackKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> HLT::MET::TrkMHTFex::m_trackKey
private
Initial value:
{
this, "TrackName", "tracks",
"Input track collection. Should be ghost-associated to the jets"}

Input tracks.

Definition at line 56 of file TrkMHTFex.h.

56 {
57 this, "TrackName", "tracks",
58 "Input track collection. Should be ghost-associated to the jets"};

◆ m_trackSelTool

ToolHandle<InDet::IInDetTrackSelectionTool> HLT::MET::TrkMHTFex::m_trackSelTool
private
Initial value:
{
this, "TrackSelTool", "InDet::InDetTrackSelectionTool/TrackSelTool",
"The track selection tool to use for tracks entering the track soft "
"term."}

The track selection tool.

Definition at line 87 of file TrkMHTFex.h.

87 {
88 this, "TrackSelTool", "InDet::InDetTrackSelectionTool/TrackSelTool",
89 "The track selection tool to use for tracks entering the track soft "
90 "term."};

◆ m_tstPtCeil

Gaudi::Property<float> HLT::MET::TrkMHTFex::m_tstPtCeil
private
Initial value:
{
this, "TrackSoftTermPtCeiling", 20. * Gaudi::Units::GeV,
"The maximum pT for tracks entering the track soft term."}

The maximum pT for tracks in the soft term.

Definition at line 92 of file TrkMHTFex.h.

92 {
93 this, "TrackSoftTermPtCeiling", 20. * Gaudi::Units::GeV,
94 "The maximum pT for tracks entering the track soft term."};

◆ m_tvaKey

SG::ReadHandleKey<jet::TrackVertexAssociation> HLT::MET::TrkMHTFex::m_tvaKey
private
Initial value:
{
this, "TVAName", "trackVtxAssoc",
"The name of the track-vertex association"}

The track-vertex assocation name.

Definition at line 63 of file TrkMHTFex.h.

63 {
64 this, "TVAName", "trackVtxAssoc",
65 "The name of the track-vertex association"};

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexKey

SG::ReadHandleKey<xAOD::VertexContainer> HLT::MET::TrkMHTFex::m_vertexKey
private
Initial value:
{
this, "VertexName", "vertices", "Input vertex collection"}

Input vertices.

Definition at line 60 of file TrkMHTFex.h.

60 {
61 this, "VertexName", "vertices", "Input vertex collection"};

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.


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