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

A tool implementing a specialized tau-jet overlap removal. More...

#include <TauAntiTauJetOverlapTool.h>

Inheritance diagram for ORUtils::TauAntiTauJetOverlapTool:
Collaboration diagram for ORUtils::TauAntiTauJetOverlapTool:

Public Member Functions

 TauAntiTauJetOverlapTool (const std::string &name)
 Create proper constructor for Athena. More...
 
virtual StatusCode findOverlaps (const xAOD::IParticleContainer &cont1, const xAOD::IParticleContainer &cont2) const override
 Identify overlapping taus and jets. More...
 
virtual StatusCode findOverlaps (const xAOD::JetContainer &jets, const xAOD::TauJetContainer &taus) const
 Identify overlapping taus and jets. More...
 
virtual void print () const =0
 Print the state of the tool. More...
 
StatusCode initialize () override final
 Initialize base class functionality. 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 > &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
 

Protected Member Functions

virtual StatusCode handleOverlap (const xAOD::IParticle *testParticle, const xAOD::IParticle *refParticle) const
 Common helper method to handle an overlap result. More...
 
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...
 
Helper methods
virtual StatusCode initializeDerived () override
 Initialize the tool. More...
 
bool isBJet (const xAOD::Jet &jet) const
 Is this jet a b-jet? Returns false if bjet ID not configured. More...
 
bool isSurvivingTau (const xAOD::TauJet &tau) const
 Is this an ID tau? Returns false if tau ID not configured. More...
 
bool isSurvivingAntiTau (const xAOD::TauJet &tau) const
 Is this an anti-tau? Returns false if anti-tau ID not configured. More...
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

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
 
Configurable properties
std::string m_bJetLabel
 Input jet decoration which labels a bjet. More...
 
std::string m_tauLabel
 Decoration labelling an IDed tau. More...
 
std::string m_antiTauLabel
 Decoration labelling an anti-tau. More...
 
float m_dR
 Flat delta-R cone for matching objects. More...
 
bool m_useRapidity
 Calculate deltaR using rapidity. More...
 
Utilities
std::unique_ptr< BJetHelperm_bJetHelper
 BJet helper. More...
 
std::unique_ptr< IParticleAssociatorm_dRMatcher
 Delta-R matcher. More...
 
std::unique_ptr< OverlapDecorationHelperm_tauDecHelper
 Decoration helper for the IDed taus. More...
 
std::unique_ptr< OverlapDecorationHelperm_antiTauDecHelper
 Decoration helper for the anti-taus. More...
 
std::string m_antiTauEventCategoryDecorName
 
SG::ReadHandleKey< xAOD::EventInfom_evtKey {this, "EventInfoKey", "EventInfo", "xAOD::EventInfo ReadHandleKey"}
 

Detailed Description

A tool implementing a specialized tau-jet overlap removal.

This tool was implemented by request for the h->bbtautau analysis which uses a loose-ish btag working point and also considers anti-tau ID.

Object precedence: ID-tau > bjet > anti-tau > light-jet

The procedure works as follows.

  1. Remove bjets overlapping with ID taus.
  2. Remove anti-taus from remaining bjets.
  3. Remove light jets from remaining ID-taus and anti-taus.

To fully utilize the above procedure you have to set the following properties:

Author
Steve Farrell Steve.nosp@m.n.Fa.nosp@m.rrell.nosp@m.@cer.nosp@m.n.ch

Definition at line 54 of file TauAntiTauJetOverlapTool.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

◆ TauAntiTauJetOverlapTool()

ORUtils::TauAntiTauJetOverlapTool::TauAntiTauJetOverlapTool ( const std::string &  name)

Create proper constructor for Athena.

Standalone constructor

Definition at line 23 of file TauAntiTauJetOverlapTool.cxx.

25  m_bJetLabel(""),
26  m_tauLabel(""),
27  m_antiTauLabel(""),
28  m_dR(0.2),
29  m_useRapidity(true),
30  m_bJetHelper(nullptr),
31  m_dRMatcher(nullptr),
32  m_antiTauDecHelper(nullptr)
33  {
34  declareProperty("BJetLabel", m_bJetLabel,
35  "Input b-jet flag. Disabled by default.");
36  declareProperty("TauLabel", m_tauLabel,
37  "Decoration which labels ID-ed taus");
38  declareProperty("AntiTauLabel", m_antiTauLabel,
39  "Decoration which labels anti-taus");
40  declareProperty("antiTauEventCategory", m_antiTauEventCategoryDecorName="antiTauEventCategory",
41  "Decoration which labels event type");
42  declareProperty("DR", m_dR, "Maximum dR for overlap match");
43  declareProperty("UseRapidity", m_useRapidity,
44  "Calculate delta-R using rapidity");
45  }

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  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 > &  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 > &  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 > &  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 > &  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

◆ findOverlaps() [1/2]

StatusCode ORUtils::TauAntiTauJetOverlapTool::findOverlaps ( const xAOD::IParticleContainer cont1,
const xAOD::IParticleContainer cont2 
) const
overridevirtual

Identify overlapping taus and jets.

Implements ORUtils::IOverlapTool.

Definition at line 83 of file TauAntiTauJetOverlapTool.cxx.

86  {
87  // Check the container types
88  if(typeid(cont1) != typeid(xAOD::JetContainer) &&
89  typeid(cont1) != typeid(ConstDataVector<xAOD::JetContainer>)) {
90  ATH_MSG_ERROR("First container arg is not of type JetContainer!");
91  return StatusCode::FAILURE;
92  }
93  if(typeid(cont2) != typeid(xAOD::TauJetContainer) &&
94  typeid(cont2) != typeid(ConstDataVector<xAOD::TauJetContainer>)) {
95  ATH_MSG_ERROR("Second container arg is not of type TauJetContainer!");
96  return StatusCode::FAILURE;
97  }
98  ATH_CHECK( findOverlaps(static_cast<const xAOD::JetContainer&>(cont1),
99  static_cast<const xAOD::TauJetContainer&>(cont2)) );
100  return StatusCode::SUCCESS;
101  }

◆ findOverlaps() [2/2]

StatusCode ORUtils::TauAntiTauJetOverlapTool::findOverlaps ( const xAOD::JetContainer jets,
const xAOD::TauJetContainer taus 
) const
virtual

Identify overlapping taus and jets.

The above method calls this one.

Definition at line 106 of file TauAntiTauJetOverlapTool.cxx.

109  {
110  ATH_MSG_DEBUG("Removing overlapping taus and jets");
111  const EventContext& ctx = Gaudi::Hive::currentContext();
112 
113  // Initialize output decorations if necessary
114  m_decHelper->initializeDecorations(taus);
115  m_decHelper->initializeDecorations(jets);
116 
117  // Start by discarding all taus which are not ID or anti-ID
118  for(const auto tau : taus) {
119  if(!m_decHelper->isSurvivingObject(*tau)) continue;
120  if(isSurvivingTau(*tau)) continue;
121  if(isSurvivingAntiTau(*tau)) continue;
122  // remove it with trivial overlap with itself
123  ATH_CHECK( handleOverlap(tau, tau) );
124  }
125 
126  // Remove bjets overlapping with ID taus
127  int ntaus = 0;
128  for(const auto tau : taus) {
129  if(!m_decHelper->isSurvivingObject(*tau)) continue;
130  // Only consider ID taus
131  if(!isSurvivingTau(*tau)) continue;
132  ntaus++;
133  for(const auto jet : jets) {
134  if(!m_decHelper->isSurvivingObject(*jet)) continue;
135  if(m_dRMatcher->objectsMatch(*tau, *jet)){
136  ATH_CHECK( handleOverlap(jet, tau) );
137  }
138  }
139  }
140 
141  // Remove anti-taus from remaining bjets
142  for(const auto jet : jets) {
143  if(!m_decHelper->isSurvivingObject(*jet)) continue;
144  if(!isBJet(*jet)) continue;
145  for(const auto tau : taus) {
146  if(!m_decHelper->isSurvivingObject(*jet)) continue;
147  if(!isSurvivingAntiTau(*tau)) continue;
148  if(m_dRMatcher->objectsMatch(*tau, *jet)) {
149  ATH_CHECK( handleOverlap(tau, jet) );
150  }
151  }
152  }
153 
154  int nantitaus = 0;
155  int antiTauCategory = 0;
157  for(const auto tau : taus) {
158  if(!m_decHelper->isSurvivingObject(*tau)) continue;
159  if(!isSurvivingAntiTau(*tau)) continue;
160  nantitaus++;
161  antiTauCategory = categoryAcc(*tau);
162  }
163 
164  int nAntiTauMax = int(ntaus<antiTauCategory);
165 
166  // AntiTauCategory = 1 for lephad event, 2 for hadhad events
167  // nAntiTauMax = 1 if we didn't get enough ID taus, 0 otherwise
169  auto eventIndex = eventInfo->eventNumber(); // pseudo-random selection of anti-taus // pseudo-random selection of anti-taus
170  if (nantitaus > 0) {
171  int selIndex = eventIndex%nantitaus;
172 
173  int nSelectedAntitaus = 0;
174  int idx = 0;
175  for(const auto tau : taus) {
176  if(!m_decHelper->isSurvivingObject(*tau)) continue;
177  if(!isSurvivingAntiTau(*tau) ) continue;
178  if (idx == selIndex && nSelectedAntitaus < nAntiTauMax) nSelectedAntitaus++;
179  else {
180  // remove excess anti-taus by applying OR fail (it trivially overlaps with itself)
181  ATH_CHECK( handleOverlap(tau, tau) );
182  }
183  idx++;
184  }
185  }
186 
187  // Remove light jets from remaining anti-taus.
188  for(const auto tau : taus) {
189  if(!m_decHelper->isSurvivingObject(*tau)) continue;
190  if(!isSurvivingAntiTau(*tau)) continue;
191  // if isSurviving
192  for(const auto jet : jets) {
193  if(!m_decHelper->isSurvivingObject(*jet)) continue;
194  // We don't need to check the bjet label, but it might save CPU.
195  if(isBJet(*jet)) continue;
196  if(m_dRMatcher->objectsMatch(*tau, *jet)){
197  ATH_CHECK( handleOverlap(jet, tau) );
198  }
199  }
200  }
201 
202  return StatusCode::SUCCESS;
203  }

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

◆ handleOverlap()

StatusCode ORUtils::BaseOverlapTool::handleOverlap ( const xAOD::IParticle testParticle,
const xAOD::IParticle refParticle 
) const
protectedvirtualinherited

Common helper method to handle an overlap result.

This method for now applies the user-defined object priorities set as the input decoration, so if the testParticle is higher prio than the refParticle, nothing happens. Otherwise, this method prints a debug message, decorates the testParticle as overlap, and optionally sets the object element link.

Definition at line 63 of file BaseOverlapTool.cxx.

66  {
67  // Apply user-priority override
68  if(!m_enableUserPrio ||
69  m_decHelper->getObjectPriority(*testParticle) <=
70  m_decHelper->getObjectPriority(*refParticle))
71  {
72  ATH_MSG_DEBUG(" Found overlap " << testParticle->type() <<
73  " pt " << testParticle->pt()*invGeV);
74  m_decHelper->setObjectFail(*testParticle);
75  if(m_objLinkHelper) {
76  ATH_CHECK( m_objLinkHelper->addObjectLink(*testParticle, *refParticle) );
77  }
78  }
79  return StatusCode::SUCCESS;
80  }

◆ initialize()

StatusCode ORUtils::BaseOverlapTool::initialize ( )
finaloverridevirtualinherited

Initialize base class functionality.

Concrete tool specific initialization can be implemented by overriding the initializeDerived() method.

Reimplemented from asg::AsgTool.

Definition at line 38 of file BaseOverlapTool.cxx.

39  {
40  ATH_MSG_DEBUG("Initializing " << name());
41  ATH_MSG_DEBUG("Base config options: InputLabel " << m_inputLabel <<
42  " OutputLabel " << m_outputLabel <<
43  " OutputPassValue " << m_outputPassValue <<
44  " UserPrio " << m_enableUserPrio);
45 
46  // Initialize the decoration helper
47  m_decHelper = std::make_unique<OverlapDecorationHelper>
49 
50  // Initialize the obj-link helper
52  m_objLinkHelper = std::make_unique<OverlapLinkHelper>("overlapObject");
53 
54  // Initialize the derived tool
56 
57  return StatusCode::SUCCESS;
58  }

◆ initializeDerived()

StatusCode ORUtils::TauAntiTauJetOverlapTool::initializeDerived ( )
overrideprotectedvirtual

Initialize the tool.

Reimplemented from ORUtils::BaseOverlapTool.

Definition at line 50 of file TauAntiTauJetOverlapTool.cxx.

51  {
52  // Initialize the b-jet helper
53  if(!m_bJetLabel.empty()) {
54  ATH_MSG_DEBUG("Configuring btag-aware OR with btag label: " << m_bJetLabel);
55  m_bJetHelper = std::make_unique<BJetHelper>(m_bJetLabel);
56  }
57 
58  // Initialize the dR matcher
60  m_dRMatcher = std::make_unique<DeltaRMatcher>(m_dR, m_useRapidity);
61 
62  // Initialize the IDed-tau decoration helper
63  if(!m_tauLabel.empty()) {
64  ATH_MSG_DEBUG("Configuring tau OR with label: " << m_tauLabel);
66  std::make_unique<OverlapDecorationHelper>
68  }
69 
70  // Initialize the anti-tau decoration helper
71  if(!m_antiTauLabel.empty()) {
72  ATH_MSG_DEBUG("Configuring anti-tau OR with label: " << m_antiTauLabel);
74  std::make_unique<OverlapDecorationHelper>
76  }
77  return StatusCode::SUCCESS;
78  }

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

◆ isBJet()

bool ORUtils::TauAntiTauJetOverlapTool::isBJet ( const xAOD::Jet jet) const
protected

Is this jet a b-jet? Returns false if bjet ID not configured.

Does not check if the jet is "surviving" OR.

Definition at line 208 of file TauAntiTauJetOverlapTool.cxx.

209  {
210  if(m_bJetHelper && m_bJetHelper->isBJet(jet)) return true;
211  return false;
212  }

◆ isSurvivingAntiTau()

bool ORUtils::TauAntiTauJetOverlapTool::isSurvivingAntiTau ( const xAOD::TauJet tau) const
protected

Is this an anti-tau? Returns false if anti-tau ID not configured.

This one does check if the tau is "surviving" OR.

Definition at line 226 of file TauAntiTauJetOverlapTool.cxx.

227  {
228  if(m_antiTauDecHelper && m_antiTauDecHelper->isSurvivingObject(tau))
229  return true;
230  return false;
231  }

◆ isSurvivingTau()

bool ORUtils::TauAntiTauJetOverlapTool::isSurvivingTau ( const xAOD::TauJet tau) const
protected

Is this an ID tau? Returns false if tau ID not configured.

This one does check if the tau is "surviving" OR.

Definition at line 217 of file TauAntiTauJetOverlapTool.cxx.

218  {
219  if(m_tauDecHelper && m_tauDecHelper->isSurvivingObject(tau)) return true;
220  return false;
221  }

◆ 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

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

std::unique_ptr<OverlapDecorationHelper> ORUtils::TauAntiTauJetOverlapTool::m_antiTauDecHelper
private

Decoration helper for the anti-taus.

Definition at line 132 of file TauAntiTauJetOverlapTool.h.

◆ m_antiTauEventCategoryDecorName

std::string ORUtils::TauAntiTauJetOverlapTool::m_antiTauEventCategoryDecorName
private

Definition at line 134 of file TauAntiTauJetOverlapTool.h.

◆ m_antiTauLabel

std::string ORUtils::TauAntiTauJetOverlapTool::m_antiTauLabel
private

Decoration labelling an anti-tau.

Definition at line 109 of file TauAntiTauJetOverlapTool.h.

◆ m_bJetHelper

std::unique_ptr<BJetHelper> ORUtils::TauAntiTauJetOverlapTool::m_bJetHelper
private

BJet helper.

Definition at line 123 of file TauAntiTauJetOverlapTool.h.

◆ m_bJetLabel

std::string ORUtils::TauAntiTauJetOverlapTool::m_bJetLabel
private

Input jet decoration which labels a bjet.

Definition at line 103 of file TauAntiTauJetOverlapTool.h.

◆ m_decHelper

std::unique_ptr<OverlapDecorationHelper> ORUtils::BaseOverlapTool::m_decHelper
protectedinherited

Helper for handling input/output decorations.

Definition at line 95 of file BaseOverlapTool.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_dR

float ORUtils::TauAntiTauJetOverlapTool::m_dR
private

Flat delta-R cone for matching objects.

Definition at line 112 of file TauAntiTauJetOverlapTool.h.

◆ m_dRMatcher

std::unique_ptr<IParticleAssociator> ORUtils::TauAntiTauJetOverlapTool::m_dRMatcher
private

Delta-R matcher.

Definition at line 126 of file TauAntiTauJetOverlapTool.h.

◆ m_enableUserPrio

bool ORUtils::BaseOverlapTool::m_enableUserPrio
protectedinherited

Enable user-priority scoring.

Definition at line 86 of file BaseOverlapTool.h.

◆ m_evtKey

SG::ReadHandleKey<xAOD::EventInfo> ORUtils::TauAntiTauJetOverlapTool::m_evtKey {this, "EventInfoKey", "EventInfo", "xAOD::EventInfo ReadHandleKey"}
private

Definition at line 135 of file TauAntiTauJetOverlapTool.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_inputLabel

std::string ORUtils::BaseOverlapTool::m_inputLabel
protectedinherited

Input object decoration which specifies which objects to look at.

Definition at line 74 of file BaseOverlapTool.h.

◆ m_linkOverlapObjects

bool ORUtils::BaseOverlapTool::m_linkOverlapObjects
protectedinherited

Flag to toggle overlap object links.

Definition at line 83 of file BaseOverlapTool.h.

◆ m_objLinkHelper

std::unique_ptr<OverlapLinkHelper> ORUtils::BaseOverlapTool::m_objLinkHelper
protectedinherited

Helper for linking overlap objects.

Definition at line 98 of file BaseOverlapTool.h.

◆ m_outputLabel

std::string ORUtils::BaseOverlapTool::m_outputLabel
protectedinherited

Output object decoration which specifies overlapping objects.

Definition at line 76 of file BaseOverlapTool.h.

◆ m_outputPassValue

bool ORUtils::BaseOverlapTool::m_outputPassValue
protectedinherited

Toggle the output flag logic.

If true, then non-overlapping objects will be assigned "true".

Definition at line 80 of file BaseOverlapTool.h.

◆ m_tauDecHelper

std::unique_ptr<OverlapDecorationHelper> ORUtils::TauAntiTauJetOverlapTool::m_tauDecHelper
private

Decoration helper for the IDed taus.

Definition at line 129 of file TauAntiTauJetOverlapTool.h.

◆ m_tauLabel

std::string ORUtils::TauAntiTauJetOverlapTool::m_tauLabel
private

Decoration labelling an IDed tau.

Definition at line 106 of file TauAntiTauJetOverlapTool.h.

◆ m_useRapidity

bool ORUtils::TauAntiTauJetOverlapTool::m_useRapidity
private

Calculate deltaR using rapidity.

Definition at line 114 of file TauAntiTauJetOverlapTool.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:
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:392
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
ORUtils::TauAntiTauJetOverlapTool::findOverlaps
virtual StatusCode findOverlaps(const xAOD::IParticleContainer &cont1, const xAOD::IParticleContainer &cont2) const override
Identify overlapping taus and jets.
Definition: TauAntiTauJetOverlapTool.cxx:84
ORUtils::TauAntiTauJetOverlapTool::isSurvivingAntiTau
bool isSurvivingAntiTau(const xAOD::TauJet &tau) const
Is this an anti-tau? Returns false if anti-tau ID not configured.
Definition: TauAntiTauJetOverlapTool.cxx:226
ORUtils::TauAntiTauJetOverlapTool::m_tauDecHelper
std::unique_ptr< OverlapDecorationHelper > m_tauDecHelper
Decoration helper for the IDed taus.
Definition: TauAntiTauJetOverlapTool.h:129
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
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
xAOD::IParticle::type
virtual Type::ObjectType type() const =0
The type of the object as a simple enumeration.
ORUtils::TauAntiTauJetOverlapTool::m_useRapidity
bool m_useRapidity
Calculate deltaR using rapidity.
Definition: TauAntiTauJetOverlapTool.h:114
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:54
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ORUtils::BaseOverlapTool::m_inputLabel
std::string m_inputLabel
Input object decoration which specifies which objects to look at.
Definition: BaseOverlapTool.h:74
ORUtils::BaseOverlapTool::m_enableUserPrio
bool m_enableUserPrio
Enable user-priority scoring.
Definition: BaseOverlapTool.h:86
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
ORUtils::BaseOverlapTool::m_objLinkHelper
std::unique_ptr< OverlapLinkHelper > m_objLinkHelper
Helper for linking overlap objects.
Definition: BaseOverlapTool.h:98
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
ORUtils::TauAntiTauJetOverlapTool::m_tauLabel
std::string m_tauLabel
Decoration labelling an IDed tau.
Definition: TauAntiTauJetOverlapTool.h:106
ORUtils::BaseOverlapTool::initializeDerived
virtual StatusCode initializeDerived()
Initialization for derived tools.
Definition: BaseOverlapTool.h:57
jet
Definition: JetCalibTools_PlotJESFactors.cxx:23
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
ORUtils::TauAntiTauJetOverlapTool::m_bJetLabel
std::string m_bJetLabel
Input jet decoration which labels a bjet.
Definition: TauAntiTauJetOverlapTool.h:103
ORUtils::BaseOverlapTool::m_outputPassValue
bool m_outputPassValue
Toggle the output flag logic.
Definition: BaseOverlapTool.h:80
ORUtils::BaseOverlapTool::m_decHelper
std::unique_ptr< OverlapDecorationHelper > m_decHelper
Helper for handling input/output decorations.
Definition: BaseOverlapTool.h:95
python.xAODType.dummy
dummy
Definition: xAODType.py:4
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
ORUtils::BaseOverlapTool::m_linkOverlapObjects
bool m_linkOverlapObjects
Flag to toggle overlap object links.
Definition: BaseOverlapTool.h:83
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
ORUtils::TauAntiTauJetOverlapTool::m_dRMatcher
std::unique_ptr< IParticleAssociator > m_dRMatcher
Delta-R matcher.
Definition: TauAntiTauJetOverlapTool.h:126
ORUtils::TauAntiTauJetOverlapTool::m_bJetHelper
std::unique_ptr< BJetHelper > m_bJetHelper
BJet helper.
Definition: TauAntiTauJetOverlapTool.h:123
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
DataVector
Derived DataVector<T>.
Definition: DataVector.h:581
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
xAOD::IParticle::pt
virtual double pt() const =0
The transverse momentum ( ) of the particle.
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
ORUtils::TauAntiTauJetOverlapTool::m_evtKey
SG::ReadHandleKey< xAOD::EventInfo > m_evtKey
Definition: TauAntiTauJetOverlapTool.h:135
ORUtils::TauAntiTauJetOverlapTool::isSurvivingTau
bool isSurvivingTau(const xAOD::TauJet &tau) const
Is this an ID tau? Returns false if tau ID not configured.
Definition: TauAntiTauJetOverlapTool.cxx:217
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
ORUtils::BaseOverlapTool::BaseOverlapTool
BaseOverlapTool(const std::string &name)
Create proper constructor for Athena.
Definition: BaseOverlapTool.cxx:20
ORUtils::BaseOverlapTool::handleOverlap
virtual StatusCode handleOverlap(const xAOD::IParticle *testParticle, const xAOD::IParticle *refParticle) const
Common helper method to handle an overlap result.
Definition: BaseOverlapTool.cxx:64
ORUtils::BaseOverlapTool::m_outputLabel
std::string m_outputLabel
Output object decoration which specifies overlapping objects.
Definition: BaseOverlapTool.h:76
ORUtils::TauAntiTauJetOverlapTool::isBJet
bool isBJet(const xAOD::Jet &jet) const
Is this jet a b-jet? Returns false if bjet ID not configured.
Definition: TauAntiTauJetOverlapTool.cxx:208
a
TList * a
Definition: liststreamerinfos.cxx:10
ORUtils::TauAntiTauJetOverlapTool::m_antiTauLabel
std::string m_antiTauLabel
Decoration labelling an anti-tau.
Definition: TauAntiTauJetOverlapTool.h:109
h
ConstDataVector
DataVector adapter that acts like it holds const pointers.
Definition: ConstDataVector.h:76
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
defineDB.jets
list jets
Definition: JetTagCalibration/share/defineDB.py:24
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
ORUtils::TauAntiTauJetOverlapTool::m_antiTauDecHelper
std::unique_ptr< OverlapDecorationHelper > m_antiTauDecHelper
Decoration helper for the anti-taus.
Definition: TauAntiTauJetOverlapTool.h:132
SG::DataProxy
Definition: DataProxy.h:44
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
ORUtils::TauAntiTauJetOverlapTool::m_dR
float m_dR
Flat delta-R cone for matching objects.
Definition: TauAntiTauJetOverlapTool.h:112
fitman.k
k
Definition: fitman.py:528
invGeV
constexpr float invGeV
Definition: PFTrackSelector.cxx:10
ORUtils::TauAntiTauJetOverlapTool::m_antiTauEventCategoryDecorName
std::string m_antiTauEventCategoryDecorName
Definition: TauAntiTauJetOverlapTool.h:134
EventInfoCnvParams::eventIndex
thread_local event_number_t eventIndex
Definition: IEvtIdModifierSvc.h:34