ATLAS Offline Software
Loading...
Searching...
No Matches
InDet::InDetTrackSystematicsTool Class Referenceabstract

#include <InDetTrackSystematicsTool.h>

Inheritance diagram for InDet::InDetTrackSystematicsTool:
Collaboration diagram for InDet::InDetTrackSystematicsTool:

Public Member Functions

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

std::unique_ptr< TFile > getFile (const std::string &) const
 open and return a file with the given name.
template<class T>
StatusCode initObject (std::unique_ptr< T > &obj, const std::string &rootFileName, const std::string &objName) const
 a function to initialize an object from a root file
bool isActive (TrackSystematic) const
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.

Protected Attributes

std::unordered_map< CP::SystematicSet, CP::SystematicSetm_sysFilterMap
const CP::SystematicSetm_activeSysts = nullptr

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

virtual ASG_TOOL_CLASS2(InDetTrackSystematicsTool, CP::ISystematicsTool, CP::IReentrantSystematicsTool) public ~InDetTrackSystematicsTool ()=default
virtual StatusCode initialize () override
 Dummy implementation of the initialisation function.
virtual bool isAffectedBySystematic (const CP::SystematicVariation &) const override
 returns: whether the tool is affected by the systematic
virtual CP::SystematicSet affectingSystematics () const override=0
 returns: list of systematics this tool can be affected by
virtual CP::SystematicSet recommendedSystematics () const override
 returns: list of recommended systematics to use with this tool
virtual StatusCode applySystematicVariation (const CP::SystematicSet &) override
 configure the tool to apply a given list of systematic variations
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

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 22 of file InDetTrackSystematicsTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Function Documentation

◆ affectingSystematics()

virtual CP::SystematicSet InDet::InDetTrackSystematicsTool::affectingSystematics ( ) const
overrideprivatepure virtual

◆ applySystematicVariation()

StatusCode InDet::InDetTrackSystematicsTool::applySystematicVariation ( const CP::SystematicSet & systematics)
overrideprivatevirtual

configure the tool to apply a given list of systematic variations

Implements CP::ISystematicsTool.

Reimplemented in InDet::InclusiveTrackFilterTool, InDet::InDetTrackBiasingTool, InDet::InDetTrackSmearingTool, InDet::InDetTrackTruthFilterTool, and InDet::JetTrackFilterTool.

Definition at line 42 of file InDetTrackSystematicsTool.cxx.

43 {
44
45 // first check if we already know this systematic configuration
46 // look for it in our filter map
47 auto sysMapItr = m_sysFilterMap.find(systematics);
48
49 if (sysMapItr == m_sysFilterMap.end()) {
50
51 // filter the input systematic with the affecting systematics
52 const CP::SystematicSet affectingSysts = affectingSystematics();
53 CP::SystematicSet filteredSysts;
55 (systematics, affectingSysts, filteredSysts) != StatusCode::SUCCESS ) {
56 ATH_MSG_ERROR( "Recieved unsupported systematics: " << systematics.name() );
57 return StatusCode::FAILURE;
58 }
59
60 // check in here for incompatible systematics
61 auto isIn = [&](InDet::TrackSystematic syst) {
62 // helper function to check if a systematic exists in the filtered set
63 return filteredSysts.find( InDet::TrackSystematicMap.at( syst ) ) != filteredSysts.end();
64 };
65 bool isSetTrkLoose = isIn( TRK_FAKE_RATE_LOOSE )
66 || isIn( TRK_EFF_LOOSE_GLOBAL )
67 || isIn( TRK_EFF_LOOSE_IBL )
68 || isIn( TRK_EFF_LOOSE_PP0 )
69 || isIn( TRK_EFF_LOOSE_PHYSMODEL );
70 bool isSetTrkTight = isIn( TRK_FAKE_RATE_TIGHT )
71 || isIn( TRK_EFF_TIGHT_GLOBAL )
72 || isIn( TRK_EFF_TIGHT_IBL )
73 || isIn( TRK_EFF_TIGHT_PP0 )
74 || isIn( TRK_EFF_TIGHT_PHYSMODEL );
75 if ( isSetTrkLoose && isSetTrkTight) {
76 ATH_MSG_ERROR( "Both Loose and TightPrimary versions of systematics are simultaneously active." );
77 return StatusCode::FAILURE;
78 }
79
80 // if the systematics set looks supported, then insert it into the map for future use
81 sysMapItr = m_sysFilterMap.insert(std::make_pair(systematics, filteredSysts)).first;
82 }
83
84 m_activeSysts = &sysMapItr->second;
85
86 return StatusCode::SUCCESS;
87 }
#define ATH_MSG_ERROR(x)
static const std::vector< std::string > systematics
const_iterator end() const
description: const iterator to the end of the set
iterator find(const SystematicVariation &sys) const
description: find an element in the set
static StatusCode filterForAffectingSystematics(const SystematicSet &systConfig, const SystematicSet &affectingSystematics, SystematicSet &filteredSystematics)
description: filter the systematics for the affected systematics returns: success guarantee: strong f...
virtual CP::SystematicSet affectingSystematics() const override=0
returns: list of systematics this tool can be affected by
std::unordered_map< CP::SystematicSet, CP::SystematicSet > m_sysFilterMap

◆ 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 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

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

◆ 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

◆ getFile()

std::unique_ptr< TFile > InDet::InDetTrackSystematicsTool::getFile ( const std::string & filename) const
protected

open and return a file with the given name.

Definition at line 99 of file InDetTrackSystematicsTool.cxx.

100 {
101 // now the files are stored in the calibration area
102 // filename is configurable in each tool and should include the path starting from the package name
103 // but defaults to recommendation for current release
104 std::string filenameWithFullPath = PathResolverFindCalibFile(filename);
105
106 return std::unique_ptr<TFile>(TFile::Open(filenameWithFullPath.data(), "READ"));;
107 }
std::string PathResolverFindCalibFile(const std::string &logical_file_name)

◆ 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 InDet::InDetTrackSystematicsTool::initialize ( void )
overrideprivatevirtual

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.

Reimplemented in InDet::InclusiveTrackFilterTool, InDet::InDetTrackBiasingTool, InDet::InDetTrackSmearingTool, InDet::InDetTrackTruthFilterTool, and InDet::JetTrackFilterTool.

Definition at line 17 of file InDetTrackSystematicsTool.cxx.

18 {
19
20 // add the affecting systematics to the registry
21 auto& registry = CP::SystematicRegistry::getInstance();
22 if (registry.registerSystematics(*this) != StatusCode::SUCCESS) {
23 ATH_MSG_ERROR( "Unable to register systematics" );
24 return StatusCode::FAILURE;
25 }
26
27 return StatusCode::SUCCESS;
28 }
static SystematicRegistry & getInstance()
Get the singleton instance of the registry for the curren thread.

◆ initObject()

template<class T>
StatusCode InDet::InDetTrackSystematicsTool::initObject ( std::unique_ptr< T > & obj,
const std::string & rootFileName,
const std::string & objName ) const
protected

a function to initialize an object from a root file

Definition at line 63 of file InDetTrackSystematicsTool.h.

64{
65 if (obj != nullptr) ATH_MSG_WARNING( obj->GetName() << " is not null, yet we are now attempting to initialize from " << rootFileName );
66 std::unique_ptr<TFile> F = getFile(rootFileName);
67 if(!F || F->IsZombie()) {
68 ATH_MSG_ERROR( "Could not open file " << rootFileName );
69 return StatusCode::FAILURE;
70 }
71 T* tempObj = nullptr;
72 F->GetObject(objName.data(), tempObj);
73 if(tempObj==nullptr) {
74 ATH_MSG_ERROR( "Could not retrieve " << objName << " from file " << rootFileName );
75 return StatusCode::FAILURE;
76 }
77 obj = std::unique_ptr<T>(static_cast<T*>(tempObj->Clone()));
78 obj->SetDirectory(0);
79 F->Clear();
80 F->Close();
81 return StatusCode::SUCCESS;
82}
#define ATH_MSG_WARNING(x)
#define F(x, y, z)
Definition MD5.cxx:112
std::unique_ptr< TFile > getFile(const std::string &) const
open and return a file with the given name.
unsigned long long T

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

◆ isActive()

bool InDet::InDetTrackSystematicsTool::isActive ( TrackSystematic syst) const
protected

Definition at line 89 of file InDetTrackSystematicsTool.cxx.

90 {
91 if (m_activeSysts == nullptr) {
92 ATH_MSG_DEBUG( "applySystematicsVariation() has not been successfully called." );
93 return false;
94 }
95 const auto it_syst = m_activeSysts->find(InDet::TrackSystematicMap.at(syst));
96 return it_syst != m_activeSysts->end();
97 }
#define ATH_MSG_DEBUG(x)

◆ isAffectedBySystematic()

bool InDet::InDetTrackSystematicsTool::isAffectedBySystematic ( const CP::SystematicVariation & systematic) const
overrideprivatevirtual

returns: whether the tool is affected by the systematic

Implements CP::ISystematicsTool.

Reimplemented in InDet::InclusiveTrackFilterTool, InDet::InDetTrackBiasingTool, InDet::InDetTrackSmearingTool, InDet::InDetTrackTruthFilterTool, and InDet::JetTrackFilterTool.

Definition at line 30 of file InDetTrackSystematicsTool.cxx.

31 {
32 auto affSysts = affectingSystematics();
33 return ( affSysts.find(systematic) != affSysts.end() );
34 }

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

◆ recommendedSystematics()

CP::SystematicSet InDet::InDetTrackSystematicsTool::recommendedSystematics ( ) const
overrideprivatevirtual

returns: list of recommended systematics to use with this tool

Implements CP::IReentrantSystematicsTool.

Reimplemented in InDet::InclusiveTrackFilterTool, InDet::InDetTrackBiasingTool, InDet::InDetTrackSmearingTool, InDet::InDetTrackTruthFilterTool, and InDet::JetTrackFilterTool.

Definition at line 36 of file InDetTrackSystematicsTool.cxx.

37 {
38 // it's probably a good idea to change this behavior to return a more intelligently thought-out list
39 return affectingSystematics();
40 }

◆ 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

◆ ~InDetTrackSystematicsTool()

virtual ASG_TOOL_CLASS2(InDetTrackSystematicsTool, CP::ISystematicsTool, CP::IReentrantSystematicsTool) public InDet::InDetTrackSystematicsTool::~InDetTrackSystematicsTool ( )
privatevirtualdefault

Member Data Documentation

◆ m_activeSysts

const CP::SystematicSet* InDet::InDetTrackSystematicsTool::m_activeSysts = nullptr
protected

Definition at line 53 of file InDetTrackSystematicsTool.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_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_sysFilterMap

std::unordered_map< CP::SystematicSet, CP::SystematicSet > InDet::InDetTrackSystematicsTool::m_sysFilterMap
protected

Definition at line 51 of file InDetTrackSystematicsTool.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.


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