ATLAS Offline Software
Loading...
Searching...
No Matches
met::METRecoTool Class Reference

#include <METRecoTool.h>

Inheritance diagram for met::METRecoTool:
Collaboration diagram for met::METRecoTool:

Public Member Functions

 METRecoTool (const std::string &name)
virtual StatusCode initialize ()
 Dummy implementation of the initialisation function.
virtual StatusCode execute () const
virtual StatusCode finalize ()
virtual void print () const
 Print the state of the tool.
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 sysInitialize () override
 Perform system initialization for an algorithm.
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
Additional helper functions, not directly mimicking Athena
template<class T>
const T * getProperty (const std::string &name) const
 Get one of the tool's properties.
const std::string & msg_level_name () const __attribute__((deprecated))
 A deprecated function for getting the message level's name.
const std::string & getName (const void *ptr) const
 Get the name of an object that is / should be in the event store.
SG::sgkey_t getKey (const void *ptr) const
 Get the (hashed) key of an object that is in the event store.

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.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

 METRecoTool ()
 Default constructor:
StatusCode buildMET (xAOD::MissingETContainer *metCont, xAOD::MissingETComponentMap *metMap) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

bool m_doMetSum
bool m_warnOfDupes
SG::WriteHandleKey< xAOD::MissingETContainerm_contname
SG::WriteHandleKey< xAOD::MissingETComponentMapm_mapname
std::string m_metfinalname
ToolHandleArray< IMETToolBasem_metbuilders
ToolHandleArray< IMETToolBasem_metrefiners
std::atomic< unsigned intm_nevt
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 58 of file METRecoTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ METRecoTool() [1/2]

met::METRecoTool::METRecoTool ( const std::string & name)

Definition at line 48 of file METRecoTool.cxx.

48 :
50 m_doMetSum(false),
51 m_metbuilders(this),
52 m_metrefiners(this),
53 m_nevt(0)
54 {
55 declareProperty( "METBuilders", m_metbuilders );
56 declareProperty( "METRefiners", m_metrefiners );
57 declareProperty( "METContainer", m_contname = "MET" );
58 declareProperty( "METComponentMap", m_mapname = "METMap" );
59 declareProperty( "WarnIfDuplicate", m_warnOfDupes = true );
60 declareProperty( "METFinalName", m_metfinalname = "" );
61 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
AsgTool(const std::string &name)
Constructor specifying the tool instance's name.
Definition AsgTool.cxx:58
ToolHandleArray< IMETToolBase > m_metbuilders
Definition METRecoTool.h:98
ToolHandleArray< IMETToolBase > m_metrefiners
Definition METRecoTool.h:99
SG::WriteHandleKey< xAOD::MissingETComponentMap > m_mapname
Definition METRecoTool.h:96
std::atomic< unsigned int > m_nevt
SG::WriteHandleKey< xAOD::MissingETContainer > m_contname
Definition METRecoTool.h:95
std::string m_metfinalname
Definition METRecoTool.h:97

◆ METRecoTool() [2/2]

met::METRecoTool::METRecoTool ( )
private

Default constructor:

Member Function Documentation

◆ buildMET()

StatusCode met::METRecoTool::buildMET ( xAOD::MissingETContainer * metCont,
xAOD::MissingETComponentMap * metMap ) const
private

Definition at line 137 of file METRecoTool.cxx.

138 {
139
140
141 MissingET* metFinal = nullptr;
142 if( m_doMetSum ) {
143 ATH_MSG_DEBUG("Building final MET sum: " << m_metfinalname);
144 metFinal = new MissingET(0.,0.,0., m_metfinalname, MissingETBase::Source::total());
145 }
146
147 // Run the MET reconstruction tools in sequence
148 for(auto tool : m_metbuilders) {
149 ATH_MSG_VERBOSE("Building new MET term with: " << tool->name() );
150 MissingET* metTerm = new MissingET(0.,0.,0.);
151 ATH_MSG_VERBOSE("Insert MET object into container");
152 metCont->push_back(metTerm);
153 ATH_MSG_VERBOSE("Insert MET object into ComponentMap");
154 MissingETComposition::add(metMap,metTerm);
155 ATH_MSG_VERBOSE("Execute tool");
156 if( tool->execute(metTerm, metMap).isFailure() ) {
157 ATH_MSG_WARNING("Failed to execute tool: " << tool->name());
158 // return StatusCode::SUCCESS;
159 }
161 // FIXME: Make a genuine decision about whether
162 // to include terms in the sum here.
165 ATH_MSG_DEBUG("Adding constructed term: " << metTerm->name() << " to sum" );
166 (*metFinal) += (*metTerm);
167 }
168 }
169
170 // Run the MET reconstruction tools in sequence
171 for(auto tool : m_metrefiners) {
172 ATH_MSG_VERBOSE("Refining MET with: " << tool->name() );
173 MissingET* metTerm = new MissingET(0.,0.,0.);
174 ATH_MSG_VERBOSE("Insert MET object into container");
175 metCont->push_back(metTerm);
176 ATH_MSG_VERBOSE("Insert MET object into ComponentMap");
177 MissingETComposition::add(metMap,metTerm);
178 ATH_MSG_VERBOSE("Execute tool");
179 if( tool->execute(metTerm, metMap).isFailure() ) {
180 ATH_MSG_WARNING("Failed to execute tool: " << tool->name());
181 // return StatusCode::SUCCESS;
182 }
183 }
184
185 if( m_doMetSum ) {
187 metFinal->setSource(source);
188 metCont->push_back(metFinal);
189 }
190 ++m_nevt;
191 return StatusCode::SUCCESS;
192 }
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
value_type push_back(value_type pElem)
Add an element to the end of the collection.
uint64_t bitmask_t
Type for status word bit mask.
MissingET_v1 MissingET
Version control by type defintion.
static Types::bitmask_t total(Region reg=Region::FullAcceptance)
Standard full reconstructed MET.
@ Refined
Indicator for genuine reconstructed MET terms.
static bool hasCategory(Types::bitmask_t bits, Category cat)
Check if bit pattern includes a given category.
static bool add(MissingETComponentMap *pMap, const MissingET *pMET, MissingETBase::Types::bitmask_t sw=MissingETBase::Status::clearedStatus())
Adding a MissingET object to the map.

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::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 }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::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< AlgTool > >::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 met::METRecoTool::execute ( ) const
virtual

Implements IMETRecoTool.

Definition at line 99 of file METRecoTool.cxx.

100 {
101 ATH_MSG_DEBUG ("In execute: " << name() << "...");
102
103 // Create a MissingETContainer with its aux store
104 auto metHandle= SG::makeHandle (m_contname);
105 //note that the method below automatically creates the MET container and its corresponding aux store (which will be named "m_contname+Aux.")
106 ATH_CHECK( metHandle.record (std::make_unique<xAOD::MissingETContainer>(), std::make_unique<xAOD::MissingETAuxContainer>()) );
107 xAOD::MissingETContainer* metCont=metHandle.ptr();
108
109
110 // Create a MissingETComponentMap with its aux store
111
112 auto metMapHandle= SG::makeHandle (m_mapname);
113 //note that the method below automatically creates the MET container and its corresponding aux store (which will be named "m_contname+Aux.")
114 ATH_CHECK( metMapHandle.record (std::make_unique<xAOD::MissingETComponentMap>(), std::make_unique<xAOD::MissingETAuxComponentMap>()) );
115 xAOD::MissingETComponentMap* metMap=metMapHandle.ptr();
116
117
118 if( buildMET(metCont, metMap).isFailure() ) {
119 ATH_MSG_WARNING("Failed in MissingET reconstruction");
120 return StatusCode::SUCCESS;
121 }
122
123 return StatusCode::SUCCESS;
124 }
#define ATH_CHECK
Evaluate an expression and check for errors.
StatusCode buildMET(xAOD::MissingETContainer *metCont, xAOD::MissingETComponentMap *metMap) const
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
MissingETComponentMap_v1 MissingETComponentMap
Version control by type definition.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::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

◆ finalize()

StatusCode met::METRecoTool::finalize ( )
virtual

Definition at line 126 of file METRecoTool.cxx.

127 {
128 ATH_MSG_INFO ("Finalizing " << name() << "...");
129 return StatusCode::SUCCESS;
130 }
#define ATH_MSG_INFO(x)

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void * ptr) const
inherited

Get the (hashed) key of an object that is in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the SG::sgkey_t key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getName
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The hashed key of the object in the store. If not found, an invalid (zero) key.

Definition at line 119 of file AsgTool.cxx.

119 {
120
121#ifdef XAOD_STANDALONE
122 // In case we use @c xAOD::TEvent, we have a direct function call
123 // for this.
124 return evtStore()->event()->getKey( ptr );
125#else
126 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
127 return ( proxy == nullptr ? 0 : proxy->sgkey() );
128#endif // XAOD_STANDALONE
129 }
ServiceHandle< StoreGateSvc > & evtStore()

◆ getName()

const std::string & asg::AsgTool::getName ( const void * ptr) const
inherited

Get the name of an object that is / should be in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the std::string name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getKey
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The string name of the object in the store. If not found, an empty string.

Definition at line 106 of file AsgTool.cxx.

106 {
107
108#ifdef XAOD_STANDALONE
109 // In case we use @c xAOD::TEvent, we have a direct function call
110 // for this.
111 return evtStore()->event()->getName( ptr );
112#else
113 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
114 static const std::string dummy = "";
115 return ( proxy == nullptr ? dummy : proxy->name() );
116#endif // XAOD_STANDALONE
117 }

◆ getProperty()

template<class T>
const T * asg::AsgTool::getProperty ( const std::string & name) const
inherited

Get one of the tool's properties.

◆ initialize()

StatusCode met::METRecoTool::initialize ( void )
virtual

Dummy implementation of the initialisation function.

It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...

Reimplemented from asg::AsgTool.

Definition at line 65 of file METRecoTool.cxx.

66 {
67 ATH_MSG_INFO ("Initializing " << name() << "...");
68
69 if( m_contname.key().empty() ) {
70 ATH_MSG_FATAL("Output MissingETContainer name must be provided.");
71 return StatusCode::FAILURE;
72 }
73
74 if( m_mapname.key().empty() ) {
75 ATH_MSG_FATAL("Output MissingETComponentMap name must be provided.");
76 return StatusCode::FAILURE;
77 }
78 ATH_CHECK( m_contname.initialize() );
79 ATH_CHECK( m_mapname.initialize() );
80
81
82 ATH_MSG_INFO ("Reconstructing MET container: " << m_contname
83 << " with composition map: " << m_mapname );
84
85 // Do we need a flag to toggle the summation?
86 // Or will we require the summation to be made
87 if( !m_metfinalname.empty() ) {
88 m_doMetSum = true;
89 ATH_MSG_INFO ("Will produce final MET sum \"" << m_metfinalname << "\"");
90 } else {
91 ATH_MSG_INFO ("Will not sum MET in this container.");
92 }
93
94 ATH_CHECK( m_metbuilders.retrieve() );
95 ATH_CHECK( m_metrefiners.retrieve() );
96 return StatusCode::SUCCESS;
97 }
#define ATH_MSG_FATAL(x)

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

◆ msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const
inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns
The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

101 {
102
103 return MSG::name( msg().level() );
104 }
MsgStream & msg() const
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition MsgLevel.cxx:19

◆ msgLvl()

bool AthCommonMsg< AlgTool >::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< AlgTool > >::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.

◆ print()

◆ 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< AlgTool > >::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< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::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< AlgTool > >::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 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ m_contname

SG::WriteHandleKey<xAOD::MissingETContainer> met::METRecoTool::m_contname
private

Definition at line 95 of file METRecoTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_doMetSum

bool met::METRecoTool::m_doMetSum
private

Definition at line 93 of file METRecoTool.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_mapname

SG::WriteHandleKey<xAOD::MissingETComponentMap> met::METRecoTool::m_mapname
private

Definition at line 96 of file METRecoTool.h.

◆ m_metbuilders

ToolHandleArray<IMETToolBase> met::METRecoTool::m_metbuilders
private

Definition at line 98 of file METRecoTool.h.

◆ m_metfinalname

std::string met::METRecoTool::m_metfinalname
private

Definition at line 97 of file METRecoTool.h.

◆ m_metrefiners

ToolHandleArray<IMETToolBase> met::METRecoTool::m_metrefiners
private

Definition at line 99 of file METRecoTool.h.

◆ m_nevt

std::atomic<unsigned int> met::METRecoTool::m_nevt
mutableprivate

Definition at line 101 of file METRecoTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_warnOfDupes

bool met::METRecoTool::m_warnOfDupes
private

Definition at line 94 of file METRecoTool.h.


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