ATLAS Offline Software
Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | Private Types | Private Member Functions | Private Attributes | List of all members
InDet::InDetTrackSmearingTool Class Referenceabstract

#include <InDetTrackSmearingTool.h>

Inheritance diagram for InDet::InDetTrackSmearingTool:
Collaboration diagram for InDet::InDetTrackSmearingTool:

Public Types

typedef xAOD::TrackParticleContainer xAODContainerType
 description: the type of xAOD container used More...
 
typedef xAODContainerType::base_value_type xAODObjectType
 description: the type of xAOD object used More...
 

Public Member Functions

virtual void print () const =0
 Print the state of the tool. More...
 
virtual void print () const
 Print the state of the tool. More...
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 
virtual CorrectionCode applyCorrection (xAODObjectType &inputObject)=0
 effects: apply the correction and store the output in the object returns: success More...
 
CorrectionCode correctedCopy (const xAODObjectType &inputObject, xAODObjectType *&outputObject)
 effects: make a clone of the input object, apply the correction to it and store the output in the object returns: success rationale: the interface class provides a default implementation in terms of applyCorrection, that the derived classes can override for efficiency More...
 
CorrectionCode applyContainerCorrection (xAODContainerType &inputContainer)
 effects: apply the correction to all the members of the container returns: overall success rationale: this will apply the correction to every single member (even if some members fail) rationale: the interface class provides a default implementation in terms of applyCorrection, that the derived classes can override for efficiency More...
 

Protected Member Functions

std::unique_ptr< TFile > getFile (const std::string &) const
 open and return a file with the given name. More...
 
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 More...
 
bool isActive (TrackSystematic) const
 
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Protected Attributes

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

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

virtual CP::CorrectionCode applyCorrection (xAOD::TrackParticle &ID) override
 Smearing method. More...
 
virtual CP::CorrectionCode correctedCopy (const xAOD::TrackParticle &in, xAOD::TrackParticle *&out) override
 
virtual CP::CorrectionCode applyContainerCorrection (xAOD::TrackParticleContainer &cont) override
 
virtual bool isAffectedBySystematic (const CP::SystematicVariation &) const override
 returns: whether the tool is affected by the systematic More...
 
virtual CP::SystematicSet affectingSystematics () const override
 returns: list of systematics this tool can be affected by More...
 
virtual CP::SystematicSet recommendedSystematics () const override
 returns: list of recommended systematics to use with this tool More...
 
virtual StatusCode applySystematicVariation (const CP::SystematicSet &) override
 configure the tool to apply a given list of systematic variations More...
 
float GetSmearD0Sigma (const xAOD::TrackParticle &) const
 Get smearing widths to add to IPs. More...
 
float GetSmearZ0Sigma (const xAOD::TrackParticle &) const
 
float readHistogram (const TH2 *, float pt, float eta) const
 
virtual StatusCode initialize () override
 Dummy implementation of the initialisation function. More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

Gaudi::Property< int > m_seed {this, "Seed", 0}
 
std::unique_ptr< TRandom3 > m_rnd
 
std::unique_ptr< TH2 > m_smearD0 = nullptr
 
std::unique_ptr< TH2 > m_smearZ0 = nullptr
 
Gaudi::Property< std::string > m_calibFileIP_CTIDE {this, "calibFileIP_CTIDE", ""}
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 29 of file InDetTrackSmearingTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

◆ xAODContainerType

description: the type of xAOD container used

Definition at line 46 of file CorrectionTool.h.

◆ xAODObjectType

description: the type of xAOD object used

Definition at line 48 of file CorrectionTool.h.

Member Function Documentation

◆ affectingSystematics()

CP::SystematicSet InDet::InDetTrackSmearingTool::affectingSystematics ( ) const
overrideprivatevirtual

returns: list of systematics this tool can be affected by

Implements InDet::IInDetTrackSmearingTool.

Definition at line 158 of file InDetTrackSmearingTool.cxx.

159  {
160  return RecommendedSystematics;
161  }

◆ applyContainerCorrection() [1/2]

CP::CorrectionCode InDet::InDetTrackSmearingTool::applyContainerCorrection ( xAOD::TrackParticleContainer cont)
overrideprivatevirtual

Implements InDet::IInDetTrackSmearingTool.

Definition at line 148 of file InDetTrackSmearingTool.cxx.

149  {
150  return TrackCorrTool_t::applyContainerCorrection(cont);
151  }

◆ applyContainerCorrection() [2/2]

CorrectionCode CP::CorrectionTool< xAOD::TrackParticleContainer >::applyContainerCorrection ( xAODContainerType inputContainer)
inherited

effects: apply the correction to all the members of the container returns: overall success rationale: this will apply the correction to every single member (even if some members fail) rationale: the interface class provides a default implementation in terms of applyCorrection, that the derived classes can override for efficiency

◆ applyCorrection() [1/2]

CP::CorrectionCode InDet::InDetTrackSmearingTool::applyCorrection ( xAOD::TrackParticle ID)
overrideprivatevirtual

Smearing method.

Implements InDet::IInDetTrackSmearingTool.

Definition at line 127 of file InDetTrackSmearingTool.cxx.

127  {
128 
129  float sigmaD0 = GetSmearD0Sigma( track );
130  float sigmaZ0 = GetSmearZ0Sigma( track );
131 
132  static const SG::AuxElement::Accessor< float > accD0( "d0" );
133  static const SG::AuxElement::Accessor< float > accZ0( "z0" );
134 
135  //NB: only call the RNG if the widths are greater than 0
136  if ( sigmaD0 > 0. ) accD0( track ) = m_rnd->Gaus( track.d0(), sigmaD0 );
137  if ( sigmaZ0 > 0. ) accZ0( track ) = m_rnd->Gaus( track.z0(), sigmaZ0 );
138 
139  return CP::CorrectionCode::Ok;
140  }

◆ applyCorrection() [2/2]

virtual CorrectionCode CP::CorrectionTool< xAOD::TrackParticleContainer >::applyCorrection ( xAODObjectType inputObject)
pure virtualinherited

effects: apply the correction and store the output in the object returns: success

◆ applySystematicVariation()

StatusCode InDet::InDetTrackSmearingTool::applySystematicVariation ( const CP::SystematicSet systs)
overrideprivatevirtual

configure the tool to apply a given list of systematic variations

Implements InDet::IInDetTrackSmearingTool.

Definition at line 169 of file InDetTrackSmearingTool.cxx.

170  {
171  // should do check in here and return StatusCode::FAILURE if LOOSE and TIGHT versions are set simultaneously
173  }

◆ correctedCopy() [1/2]

CP::CorrectionCode InDet::InDetTrackSmearingTool::correctedCopy ( const xAOD::TrackParticle in,
xAOD::TrackParticle *&  out 
)
overrideprivatevirtual

Implements InDet::IInDetTrackSmearingTool.

Definition at line 142 of file InDetTrackSmearingTool.cxx.

144  {
145  return TrackCorrTool_t::correctedCopy(in, out);
146  }

◆ correctedCopy() [2/2]

CorrectionCode CP::CorrectionTool< xAOD::TrackParticleContainer >::correctedCopy ( const xAODObjectType inputObject,
xAODObjectType *&  outputObject 
)
inherited

effects: make a clone of the input object, apply the correction to it and store the output in the object returns: success rationale: the interface class provides a default implementation in terms of applyCorrection, that the derived classes can override for efficiency

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< 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  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

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

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

Definition at line 199 of file AthCommonDataStore.h.

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

◆ declareProperty() [1/6]

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

Declare a new Gaudi property.

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

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

Definition at line 245 of file AthCommonDataStore.h.

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

◆ declareProperty() [2/6]

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

Declare a new Gaudi property.

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

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

Definition at line 221 of file AthCommonDataStore.h.

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

◆ declareProperty() [3/6]

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

Definition at line 259 of file AthCommonDataStore.h.

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

◆ declareProperty() [4/6]

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

Declare a new Gaudi property.

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

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

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

Declare a new Gaudi property.

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

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

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

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

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

95 { return m_detStore; }

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

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

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  }

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

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

◆ GetSmearD0Sigma()

float InDet::InDetTrackSmearingTool::GetSmearD0Sigma ( const xAOD::TrackParticle track) const
private

Get smearing widths to add to IPs.

Definition at line 68 of file InDetTrackSmearingTool.cxx.

68  {
69  float pt = 1.e-3*track.pt(); // need to convert pt to GeV
70  float eta = track.eta();
71  float sigma_D0 = 0.f;
72 
73  bool isActiveD0Meas = isActive(TRK_RES_D0_MEAS);
74 
75  if (isActiveD0Meas == 0) {
76  // pass-through D0Meas, return sigma_D0 early
77  return std::sqrt(sigma_D0);
78  }
79 
80  TH2* d0hist = m_smearD0.get();
81 
82  if (d0hist == nullptr) ATH_MSG_ERROR( "d0 histogram is null!" );
83 
84  float d0Smear = readHistogram(d0hist, pt, eta);
85 
86  // Apply a minimum smearing in the case where there is an empty bin
87  // TODO: Remove for next round of recommendations
88  if (d0Smear == 0) {
89  d0Smear = 0.002;
90  }
91 
92  sigma_D0 += d0Smear*d0Smear;
93 
94  return std::sqrt(sigma_D0);
95  }

◆ GetSmearZ0Sigma()

float InDet::InDetTrackSmearingTool::GetSmearZ0Sigma ( const xAOD::TrackParticle track) const
private

Definition at line 98 of file InDetTrackSmearingTool.cxx.

98  {
99  float pt = 1.e-3*track.pt(); // need to convert pt to GeV
100  float eta = track.eta();
101  float sigma_Z0 = 0.f;
102 
103  bool isActiveZ0Meas = isActive(TRK_RES_Z0_MEAS);
104 
105  if (isActiveZ0Meas == 0) {
106  // pass-through Z0Meas, return sigma_Z0 early
107  return std::sqrt(sigma_Z0);
108  }
109 
110  TH2* z0hist = m_smearZ0.get();
111 
112  if (z0hist == nullptr) ATH_MSG_ERROR( "z0 histogram is null!" );
113 
114  float z0Smear = readHistogram(z0hist, pt, eta);
115 
116  // Apply a minimum smearing in the case where there is an empty bin
117  // TODO: Remove for next round of recommendations
118  if (z0Smear == 0) {
119  z0Smear = 0.03;
120  }
121 
122  sigma_Z0 += z0Smear*z0Smear;
123 
124  return std::sqrt(sigma_Z0);
125  }

◆ initialize()

StatusCode InDet::InDetTrackSystematicsTool::initialize ( )
overrideprivatevirtualinherited

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, and InDet::InDetTrackBiasingTool.

Definition at line 17 of file InDetTrackSystematicsTool.cxx.

18  {
19 
20  // add the affecting systematics to the registry
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  }

◆ initObject()

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

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 }

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

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  }

◆ isAffectedBySystematic()

bool InDet::InDetTrackSmearingTool::isAffectedBySystematic ( const CP::SystematicVariation syst) const
overrideprivatevirtual

returns: whether the tool is affected by the systematic

Implements InDet::IInDetTrackSmearingTool.

Definition at line 153 of file InDetTrackSmearingTool.cxx.

154  {
156  }

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

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

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

◆ 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() [1/2]

void asg::AsgTool::print ( ) const
virtualinherited

◆ print() [2/2]

virtual void asg::IAsgTool::print ( ) const
pure virtualinherited

◆ readHistogram()

float InDet::InDetTrackSmearingTool::readHistogram ( const TH2 *  histogram,
float  pt,
float  eta 
) const
private

Definition at line 175 of file InDetTrackSmearingTool.cxx.

175  {
176  // safety measure:
177  if( eta>2.499 ) eta= 2.499;
178  if( eta<-2.499 ) eta=-2.499;
179  if ( pt < histogram->GetXaxis()->GetXmin() ) {
180  pt = histogram->GetXaxis()->GetXmin() + 0.01;
181  }
182  if ( pt >= histogram->GetXaxis()->GetXmax() ) {
183  pt = histogram->GetXaxis()->GetXmax() - 0.01;
184  }
185 
186  int bin = histogram->FindFixBin(pt, eta);
187  float val = histogram->GetBinContent(bin);
188  return val;
189  }

◆ recommendedSystematics()

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

returns: list of recommended systematics to use with this tool

Implements InDet::IInDetTrackSmearingTool.

Definition at line 163 of file InDetTrackSmearingTool.cxx.

164  {
165  // return InDetTrackSystematicsTool::recommendedSystematics();
166  return RecommendedSystematics;
167  }

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

◆ 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  {
365  handlesArray.renounce();
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 DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

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

Member Data Documentation

◆ m_activeSysts

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

Definition at line 53 of file InDetTrackSystematicsTool.h.

◆ m_calibFileIP_CTIDE

Gaudi::Property<std::string> InDet::InDetTrackSmearingTool::m_calibFileIP_CTIDE {this, "calibFileIP_CTIDE", ""}
private

Definition at line 84 of file InDetTrackSmearingTool.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_rnd

std::unique_ptr<TRandom3> InDet::InDetTrackSmearingTool::m_rnd
private

Definition at line 71 of file InDetTrackSmearingTool.h.

◆ m_seed

Gaudi::Property<int> InDet::InDetTrackSmearingTool::m_seed {this, "Seed", 0}
private

Definition at line 70 of file InDetTrackSmearingTool.h.

◆ m_smearD0

std::unique_ptr<TH2> InDet::InDetTrackSmearingTool::m_smearD0 = nullptr
private

Definition at line 80 of file InDetTrackSmearingTool.h.

◆ m_smearZ0

std::unique_ptr<TH2> InDet::InDetTrackSmearingTool::m_smearZ0 = nullptr
private

Definition at line 81 of file InDetTrackSmearingTool.h.

◆ m_sysFilterMap

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

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:
python.Dso.registry
registry
Definition: Control/AthenaServices/python/Dso.py:158
InDet::InDetTrackSmearingTool::m_rnd
std::unique_ptr< TRandom3 > m_rnd
Definition: InDetTrackSmearingTool.h:71
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:407
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
SG::Accessor
Helper class to provide type-safe access to aux data.
Definition: Control/AthContainers/AthContainers/Accessor.h:68
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:83
InDet::InDetTrackSystematicsTool::isActive
bool isActive(TrackSystematic) const
Definition: InDetTrackSystematicsTool.cxx:89
test_pyathena.pt
pt
Definition: test_pyathena.py:11
python.AthDsoLogger.out
out
Definition: AthDsoLogger.py:70
bin
Definition: BinsDiffFromStripMedian.h:43
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
InDet::InDetTrackSystematicsTool::isAffectedBySystematic
virtual bool isAffectedBySystematic(const CP::SystematicVariation &) const override
returns: whether the tool is affected by the systematic
Definition: InDetTrackSystematicsTool.cxx:30
InDet::InDetTrackSystematicsTool::m_activeSysts
const CP::SystematicSet * m_activeSysts
Definition: InDetTrackSystematicsTool.h:53
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
InDet::TrackSystematicMap
static const std::unordered_map< InDet::TrackSystematic, CP::SystematicVariation, std::hash< int > > TrackSystematicMap
Definition: InDetTrackSystematics.h:50
InDet::InDetTrackSmearingTool::readHistogram
float readHistogram(const TH2 *, float pt, float eta) const
Definition: InDetTrackSmearingTool.cxx:175
InDet::InDetTrackSmearingTool::GetSmearZ0Sigma
float GetSmearZ0Sigma(const xAOD::TrackParticle &) const
Definition: InDetTrackSmearingTool.cxx:98
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
CalibDbCompareRT.dummy
dummy
Definition: CalibDbCompareRT.py:59
CP::SystematicSet::end
const_iterator end() const
description: const iterator to the end of the set
Definition: SystematicSet.h:59
InDet::InDetTrackSystematicsTool::getFile
std::unique_ptr< TFile > getFile(const std::string &) const
open and return a file with the given name.
Definition: InDetTrackSystematicsTool.cxx:99
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
InDet::InDetTrackSmearingTool::m_smearZ0
std::unique_ptr< TH2 > m_smearZ0
Definition: InDetTrackSmearingTool.h:81
InDet::InDetTrackSmearingTool::m_smearD0
std::unique_ptr< TH2 > m_smearD0
Definition: InDetTrackSmearingTool.h:80
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
InDet::TRK_RES_D0_MEAS
@ TRK_RES_D0_MEAS
Definition: InDetTrackSystematics.h:16
InDet::RecommendedSystematics
static const CP::SystematicSet RecommendedSystematics
Definition: InDetTrackSmearingTool.cxx:28
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:283
CP::SystematicSet::find
iterator find(const SystematicVariation &sys) const
description: find an element in the set
Definition: SystematicSet.h:63
CP::CorrectionCode::Ok
@ Ok
The correction was done successfully.
Definition: CorrectionCode.h:38
a
TList * a
Definition: liststreamerinfos.cxx:10
h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
Pythia8_RapidityOrderMPI.val
val
Definition: Pythia8_RapidityOrderMPI.py:14
F
#define F(x, y, z)
Definition: MD5.cxx:112
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
InDet::TRK_RES_Z0_MEAS
@ TRK_RES_Z0_MEAS
Definition: InDetTrackSystematics.h:17
CaloCellTimeCorrFiller.filename
filename
Definition: CaloCellTimeCorrFiller.py:23
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
xAOD::track
@ track
Definition: TrackingPrimitives.h:513
python.PyAthena.obj
obj
Definition: PyAthena.py:132
SG::DataProxy
Definition: DataProxy.h:45
InDet::InDetTrackSystematicsTool::applySystematicVariation
virtual StatusCode applySystematicVariation(const CP::SystematicSet &) override
configure the tool to apply a given list of systematic variations
Definition: InDetTrackSystematicsTool.cxx:42
InDet::InDetTrackSmearingTool::GetSmearD0Sigma
float GetSmearD0Sigma(const xAOD::TrackParticle &) const
Get smearing widths to add to IPs.
Definition: InDetTrackSmearingTool.cxx:68
TSU::T
unsigned long long T
Definition: L1TopoDataTypes.h:35
histogram
std::string histogram
Definition: chains.cxx:52
fitman.k
k
Definition: fitman.py:528
CP::SystematicRegistry::getInstance
static SystematicRegistry & getInstance()
Get the singleton instance of the registry for the curren thread.
Definition: SystematicRegistry.cxx:25