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

#include <DiTauTruthMatchingTool.h>

Inheritance diagram for TauAnalysisTools::DiTauTruthMatchingTool:
Collaboration diagram for TauAnalysisTools::DiTauTruthMatchingTool:

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...
 
virtual StatusCode sysInitialize ()
 Function initialising the tool in the correct way in Athena. 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 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
 

Protected Member Functions

StatusCode retrieveTruthTaus (TruthTausEvent &truthTausEvent) 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

TruthTausEvent m_truthTausEvent
 
bool m_bTruthTauAvailable {}
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

virtual StatusCode initialize ()
 
virtual void getTruth (const xAOD::DiTauJet &xDiTau)
 get pointer to the truth matched particle. More...
 
virtual void getTruth (const std::vector< const xAOD::DiTauJet * > &vDiTaus)
 
StatusCode findTruthTau (const xAOD::DiTauJet &xDiTau)
 
StatusCode checkTruthMatch (const xAOD::DiTauJet &xDiTau, const xAOD::TruthParticleContainer &xTauContainer) const
 
StatusCode truthMatch (const TLorentzVector &vSubjetTLV, const xAOD::TruthParticleContainer &xTruthTauContainer, const xAOD::TruthParticle *&xTruthMatch, const xAOD::Jet *&xTruthJetMatch, TruthMatchedParticleType &eTruthMatchedParticleType) const
 
virtual void setTruthMatchingMode () override
 
virtual StatusCode retrieveTruthTaus () override
 
virtual StatusCode retrieveTruthTaus (ITruthTausEvent &truthTausEvent) const override
 
StatusCode buildTruthTausFromTruthParticles (TruthTausEvent &truthTausEvent) const
 
StatusCode examineTruthTau (const xAOD::TruthParticle &xTruthParticle) const
 
StatusCode examineTruthTauDecay (const xAOD::TruthParticle &xTruthParticle, TauTruthInfo &truthInfo) const
 
void printDecay (const xAOD::TruthParticle &xTruthParticle, int depth=0) const
 
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< double > m_dMaxDeltaR {this, "MaxDeltaR", 0.2}
 
CxxUtils::CachedValue< bool > m_bIsTruthMatchedAvailable
 
CxxUtils::CachedValue< bool > m_bIsTruthParticleLinkAvailable
 
SG::ConstAccessor< double > m_accPtVis
 
SG::ConstAccessor< double > m_accEtaVis
 
SG::ConstAccessor< double > m_accPhiVis
 
SG::ConstAccessor< double > m_accMVis
 
Gaudi::Property< bool > m_bWriteInvisibleFourMomentum { this, "WriteInvisibleFourMomentum", false}
 
Gaudi::Property< bool > m_bWriteVisibleChargedFourMomentum { this, "WriteVisibleChargedFourMomentum", false}
 
Gaudi::Property< bool > m_bWriteVisibleNeutralFourMomentum { this, "WriteVisibleNeutralFourMomentum", false}
 
Gaudi::Property< bool > m_bWriteDecayModeVector { this, "WriteDecayModeVector", true}
 
Gaudi::Property< bool > m_bWriteVertices { this, "WriteVertices", true}
 
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthTauInputContainer { this, "TruthTauContainerName", "TruthTaus", "Truth tau input container name (truth matching mode)" }
 
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthParticleContainer { this, "TruthParticleContainerName", "TruthParticles", "Truth particles input container name" }
 
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthElectronContainer { this, "TruthElectronContainerName", "TruthElectrons", "Truth electrons input container name" }
 
SG::ReadHandleKey< xAOD::TruthParticleContainerm_truthMuonContainer { this, "TruthMuonContainerName", "TruthMuons", "Truth muons input container name" }
 
SG::ReadHandleKey< xAOD::JetContainerm_truthJetContainer { this, "TruthJetContainerName", "AntiKt4TruthDressedWZJets", "Truth jets input container name" }
 
SG::WriteHandleKey< xAOD::TruthParticleContainerm_truthTauOutputContainer { this, "NewTruthTauContainerName", "TruthTaus", "Truth tau output container name" }
 
bool m_truthMatchingMode = false
 
asg::AnaToolHandle< IMCTruthClassifierm_tMCTruthClassifier
 
MetaStore_t m_inputMetaStore
 Object accessing the input metadata store. More...
 
MetaStore_t m_outputMetaStore
 Object accessing the output metadata store. More...
 
bool m_beginInputFileCalled
 Flag helping to discover when the tool misses the opening of the first input file. More...
 
bool m_useIncidents
 
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
 

Callback functions helping in metadata reading/writing

void setUseIncidents (const bool flag)
 
virtual void handle (const Incident &inc)
 Function receiving incidents from IncidentSvc/TEvent. More...
 
virtual StatusCode beginInputFile ()
 Function called when a new input file is opened. More...
 
virtual StatusCode endInputFile ()
 Function called when the currently open input file got completely processed. More...
 
virtual StatusCode beginEvent ()
 Function called when a new events is loaded. More...
 
virtual StatusCode metaDataStop ()
 Function called when the tool should write out its metadata. More...
 

Detailed Description

Definition at line 27 of file DiTauTruthMatchingTool.h.

Member Typedef Documentation

◆ MetaStore_t

Type of the metadata store object in Athena.

Definition at line 66 of file AsgMetadataTool.h.

◆ MetaStorePtr_t

Type of the metadata store pointer in standalone mode.

Definition at line 68 of file AsgMetadataTool.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Function Documentation

◆ beginEvent()

StatusCode asg::AsgMetadataTool::beginEvent ( )
protectedvirtualinherited

Function called when a new events is loaded.

Dummy implementation that can be overridden by the derived tool.

Reimplemented in AsgElectronEfficiencyCorrectionTool, TrigConf::xAODConfigTool, TauAnalysisTools::TauSelectionTool, TauAnalysisTools::DiTauSelectionTool, TauAnalysisTools::CommonSmearingTool, TauAnalysisTools::DiTauEfficiencyCorrectionsTool, Trig::TrigDecisionTool, TauAnalysisTools::TauEfficiencyCorrectionsTool, Trig::TrigConfBunchCrossingTool, xAODMaker::TriggerMenuMetaDataTool, and Trig::xAODBunchCrossingTool.

Definition at line 201 of file AsgMetadataTool.cxx.

201  {
202 
203  // Return gracefully:
204  return StatusCode::SUCCESS;
205  }

◆ beginInputFile()

StatusCode asg::AsgMetadataTool::beginInputFile ( )
protectedvirtualinherited

Function called when a new input file is opened.

Dummy implementation that can be overridden by the derived tool.

Reimplemented in AsgElectronEfficiencyCorrectionTool, TrigConf::xAODConfigTool, PMGTools::PMGTruthWeightTool, BookkeeperTool, BookkeeperDumperTool, Trig::TrigDecisionTool, xAODMaker::TriggerMenuMetaDataTool, xAODMaker::TruthMetaDataTool, Trig::xAODBunchCrossingTool, TauAnalysisTools::TauEfficiencyCorrectionsTool, and TauAnalysisTools::TauSmearingTool.

Definition at line 185 of file AsgMetadataTool.cxx.

185  {
186 
187  // Return gracefully:
188  return StatusCode::SUCCESS;
189  }

◆ buildTruthTausFromTruthParticles()

StatusCode BuildTruthTaus::buildTruthTausFromTruthParticles ( TruthTausEvent &  truthTausEvent) const
privateinherited

Definition at line 161 of file BuildTruthTaus.cxx.

162 {
163  static const SG::Accessor<char> dressedPhotonAcc ("dressedPhoton");
164 
165  bool copyDressedPhotons = false;
166 
167  for (auto xTruthParticle : *truthTausEvent.m_xTruthParticleContainer)
168  {
169  if ( xTruthParticle->isTau() )
170  {
171  auto xTruthTau = std::make_unique<xAOD::TruthParticle>();
172  xTruthTau->makePrivateStore( *xTruthParticle );
173 
174  if ( examineTruthTau(*xTruthTau).isFailure() )
175  {
176  continue;
177  }
178 
179  // The dressedPhoton decoration will likely be unlocked and thus will
180  // not be copied by the above. Copy it explicitly in that case.
181  // See ATLASRECTS-8008.
182  // First time through the loop we check to see if the decoration
183  // needs to be copied and remember for subsequent iterations
184  // (during which the destination decoration will have already
185  // been created).
186  if (truthTausEvent.m_xTruthTauContainer->empty() &&
187  dressedPhotonAcc.isAvailable ( *xTruthParticle ) &&
188  !dressedPhotonAcc.isAvailable ( *xTruthTau ))
189  {
190  copyDressedPhotons = true;
191  }
192  if (copyDressedPhotons)
193  {
194  dressedPhotonAcc( *xTruthTau ) = dressedPhotonAcc( *xTruthParticle );
195  }
196 
197  // Run classification
198  auto pClassification = m_tMCTruthClassifier->particleTruthClassifier(xTruthTau.get());
199  static const SG::Accessor<unsigned int> decClassifierParticleType("classifierParticleType");
200  static const SG::Accessor<unsigned int> decClassifierParticleOrigin("classifierParticleOrigin");
201  decClassifierParticleType(*xTruthTau) = pClassification.first;
202  decClassifierParticleOrigin(*xTruthTau) = pClassification.second;
203 
204  // create link to the original TruthParticle
205  ElementLink < xAOD::TruthParticleContainer > lTruthParticleLink(xTruthParticle, *truthTausEvent.m_xTruthParticleContainer);
206  static const SG::Accessor<ElementLink< xAOD::TruthParticleContainer > > accOriginalTruthParticle("originalTruthParticle");
207  accOriginalTruthParticle(*xTruthTau) = lTruthParticleLink;
208 
209  truthTausEvent.m_xTruthTauContainer->push_back(std::move(xTruthTau));
210  }
211  }
212  return StatusCode::SUCCESS;
213 }

◆ checkTruthMatch()

StatusCode DiTauTruthMatchingTool::checkTruthMatch ( const xAOD::DiTauJet xDiTau,
const xAOD::TruthParticleContainer xTauContainer 
) const
private

Definition at line 99 of file DiTauTruthMatchingTool.cxx.

100 {
101  std::vector<const xAOD::TruthParticle*> vTruthMatch;
102  std::vector<TruthMatchedParticleType> vTruthMatchedParticleType;
103  std::vector<const xAOD::Jet*> vTruthJetMatch;
104 
105 
106  xAOD::TruthParticleContainer xRemainingTruthTaus = xTruthTauContainer;
107 
108  static const SG::Decorator<char> decIsTruthMatched("IsTruthMatched");
109  static const SG::Decorator<char> decIsTruthHadronic("IsTruthHadronic");
110  static const SG::ConstAccessor<int> accNSubjets("n_subjets");
111  static const SG::ConstAccessor<char> accIsTruthHadronic("IsTruthHadronic");
112 
113  int n_subjets = 0;
114  if(!(accNSubjets.isAvailable(xDiTau))){
115  // n_subjets decoration is not available, recalculation on the fly
116  while (xDiTau.subjetPt(n_subjets) > 0. )
117  {
118  n_subjets++;
119  }
120  } else {
121  n_subjets = accNSubjets(xDiTau);
122  }
123 
124  // set default values for each subjet
125  for (int i = 0; i < n_subjets; ++i)
126  {
127  const xAOD::TruthParticle* xTruthMatch = nullptr;
128  TruthMatchedParticleType eTruthMatchedParticleType = Unknown;
129 
130  vTruthMatch.push_back(xTruthMatch);
131  vTruthMatchedParticleType.push_back(eTruthMatchedParticleType);
132 
133  const xAOD::Jet* xTruthJetMatch = nullptr;
134  vTruthJetMatch.push_back(xTruthJetMatch);
135  }
136 
137  // truthmatching for subjets:
138  for (int i = 0; i < n_subjets; ++i)
139  {
140  TLorentzVector vSubjetTLV;
141  vSubjetTLV.SetPtEtaPhiE(xDiTau.subjetPt(i),
142  xDiTau.subjetEta(i),
143  xDiTau.subjetPhi(i),
144  xDiTau.subjetE(i));
145  if ( truthMatch(vSubjetTLV,
146  xRemainingTruthTaus,
147  vTruthMatch.at(i),
148  vTruthJetMatch.at(i),
149  vTruthMatchedParticleType.at(i)).isFailure() )
150  {
151  ATH_MSG_WARNING("There was a failure in matching truth taus with subjet " << i);
152  return StatusCode::FAILURE;
153  }
154  if (vTruthMatch.at(i) && vTruthMatchedParticleType.at(i) == TruthHadronicTau)
155  {
156  xRemainingTruthTaus.erase( std::find(xRemainingTruthTaus.begin(),
157  xRemainingTruthTaus.end(),
158  vTruthMatch.at(i)) );
159  }
160  }
161 
162  // create links for jets
163  std::vector< ElementLink < xAOD::JetContainer > > vTruthJetLinks;
164  for (int i = 0; i < n_subjets; ++i)
165  {
166  const xAOD::Jet* xTruthJetMatch = vTruthJetMatch.at(i);
167  if(xTruthJetMatch){
168  ElementLink < xAOD::JetContainer > lTruthParticleLink(xTruthJetMatch, *m_truthTausEvent.m_xTruthJetContainerConst);
169  vTruthJetLinks.push_back(lTruthParticleLink);
170  }
171  else
172  {
173  ElementLink < xAOD::JetContainer > lTruthParticleLink;
174  vTruthJetLinks.push_back(lTruthParticleLink);
175  }
176  }
178  decTruthJetLinks ("truthJetLinks");
179  decTruthJetLinks(xDiTau) = std::move(vTruthJetLinks);
180 
181 
182  bool bTruthMatched = true;
183 
184  // create link to the original TruthParticle
185  std::vector< ElementLink < xAOD::TruthParticleContainer > > vTruthLinks;
186  for (int i = 0; i < n_subjets; ++i)
187  {
188  const xAOD::TruthParticle* xTruthMatch = vTruthMatch.at(i);
189  TruthMatchedParticleType eTruthMatchedParticleType = vTruthMatchedParticleType.at(i);
190  if (xTruthMatch)
191  {
192  if (eTruthMatchedParticleType == TruthHadronicTau)
193  {
194  ElementLink < xAOD::TruthParticleContainer > lTruthParticleLink(xTruthMatch, xTruthTauContainer);
195  vTruthLinks.push_back(lTruthParticleLink);
196  }
197  else if (eTruthMatchedParticleType == TruthMuon)
198  {
199  ElementLink <xAOD::TruthParticleContainer> lTruthParticleLink(xTruthMatch, *m_truthTausEvent.m_xTruthMuonContainerConst);
200  vTruthLinks.push_back(lTruthParticleLink);
201  }
202  else if (eTruthMatchedParticleType == TruthElectron)
203  {
204  ElementLink <xAOD::TruthParticleContainer> lTruthParticleLink(xTruthMatch, *m_truthTausEvent.m_xTruthElectronContainerConst);
205  vTruthLinks.push_back(lTruthParticleLink);
206  }
207  }
208  else
209  {
211  vTruthLinks.push_back(lTruthParticleLink);
212 
213  // ditau is not truth matched if one of the two leading subjets is not truth matched
214  if (i == 0 || i == 1) bTruthMatched = false;
215  }
216  }
217 
219  decTruthParticleLinks ("truthParticleLinks");
220  decTruthParticleLinks(xDiTau) = vTruthLinks;
222  {
224  decTruthTaus ("TruthTaus");
225  decTruthTaus(xDiTau) = vTruthLinks;
226  }
227 
228  static const SG::Decorator<float> decTruthLeadPt("TruthVisLeadPt");
229  static const SG::Decorator<float> decTruthLeadEta("TruthVisLeadEta");
230  static const SG::Decorator<float> decTruthLeadPhi("TruthVisLeadPhi");
231  static const SG::Decorator<float> decTruthLeadM("TruthVisLeadM");
232  static const SG::Decorator<float> decTruthLeadPdgID("TruthLeadPdgID");
233  static const SG::Decorator<float> decTruthSubleadPt("TruthVisSubleadPt");
234  static const SG::Decorator<float> decTruthSubleadEta("TruthVisSubleadEta");
235  static const SG::Decorator<float> decTruthSubleadPhi("TruthVisSubleadPhi");
236  static const SG::Decorator<float> decTruthSubleadM("TruthVisSubleadM");
237  static const SG::Decorator<float> decTruthSubleadPdgID("TruthSubleadPdgID");
238  static const SG::Decorator<float> decTruthDeltaR("TruthVisDeltaR");
239  static const SG::Decorator<float> decTruthMass("TruthVisMass");
240 
241  // the ditau candidate should have at least 2 subjets to be truth matched
242  if ( n_subjets < 2) {
243  decIsTruthMatched(xDiTau) = (char)false;
244  decIsTruthHadronic(xDiTau) = (char)false;
245  decTruthLeadPt(xDiTau) = -1234.;
246  decTruthLeadEta(xDiTau) = -1234.;
247  decTruthLeadPhi(xDiTau) = -1234.;
248  decTruthLeadM(xDiTau) = -1234.;
249  decTruthLeadPdgID(xDiTau) = -1234.;
250  decTruthSubleadPt(xDiTau) = -1234.;
251  decTruthSubleadEta(xDiTau) = -1234.;
252  decTruthSubleadPhi(xDiTau) = -1234.;
253  decTruthSubleadM(xDiTau) = -1234.;
254  decTruthSubleadPdgID(xDiTau) = -1234.;
255  decTruthDeltaR(xDiTau) = -1234.;
256  decTruthMass(xDiTau) = -1234.;
257  return StatusCode::SUCCESS;
258  }
259 
260  decIsTruthMatched(xDiTau) = (char)bTruthMatched;
261  if (bTruthMatched)
262  {
263  // ditau is hadronic if two leading subjets are truth matched with hadronic decay
264  decIsTruthHadronic(xDiTau) = (char)(vTruthMatchedParticleType[0]==TruthHadronicTau
265  && vTruthMatchedParticleType[1]==TruthHadronicTau );
266  }
267  else
268  decIsTruthHadronic(xDiTau) = (char)false;
269 
270  if (accIsTruthHadronic(xDiTau))
271  {
272  TLorentzVector tlvTruthTau1;
273  TLorentzVector tlvTruthTau2;
274  tlvTruthTau1.SetPtEtaPhiM(m_accPtVis(*(*vTruthLinks.at(0))),
275  m_accEtaVis(*(*vTruthLinks.at(0))),
276  m_accPhiVis(*(*vTruthLinks.at(0))),
277  m_accMVis(*(*vTruthLinks.at(0))));
278  tlvTruthTau2.SetPtEtaPhiM(m_accPtVis(*(*vTruthLinks.at(1))),
279  m_accEtaVis(*(*vTruthLinks.at(1))),
280  m_accPhiVis(*(*vTruthLinks.at(1))),
281  m_accMVis(*(*vTruthLinks.at(1))));
282 
283  decTruthLeadPt(xDiTau) = std::max(tlvTruthTau1.Pt(), tlvTruthTau2.Pt());
284  decTruthLeadEta(xDiTau) = (tlvTruthTau1.Pt() > tlvTruthTau2.Pt()) ? tlvTruthTau1.Eta() : tlvTruthTau2.Eta();
285  decTruthLeadPhi(xDiTau) = (tlvTruthTau1.Pt() > tlvTruthTau2.Pt()) ? tlvTruthTau1.Phi() : tlvTruthTau2.Phi();
286  decTruthLeadM(xDiTau) = (tlvTruthTau1.Pt() > tlvTruthTau2.Pt()) ? tlvTruthTau1.M() : tlvTruthTau2.M();
287  decTruthLeadPdgID(xDiTau) = (*vTruthLinks.at(0))->pdgId();
288  decTruthSubleadPt(xDiTau) = std::min(tlvTruthTau1.Pt(), tlvTruthTau2.Pt());
289  decTruthSubleadEta(xDiTau) = (tlvTruthTau1.Pt() > tlvTruthTau2.Pt()) ? tlvTruthTau2.Eta() : tlvTruthTau1.Eta();
290  decTruthSubleadPhi(xDiTau) = (tlvTruthTau1.Pt() > tlvTruthTau2.Pt()) ? tlvTruthTau2.Phi() : tlvTruthTau1.Phi();
291  decTruthSubleadM(xDiTau) = (tlvTruthTau1.Pt() > tlvTruthTau2.Pt()) ? tlvTruthTau2.M() : tlvTruthTau1.M();
292  decTruthSubleadPdgID(xDiTau) = (*vTruthLinks.at(1))->pdgId();
293  decTruthDeltaR(xDiTau) = tlvTruthTau1.DeltaR(tlvTruthTau2);
294  decTruthMass(xDiTau) = (tlvTruthTau1 + tlvTruthTau2).M();
295  }
296  else {
297  // set to a default value
298  decTruthLeadPt(xDiTau) = -1234.;
299  decTruthLeadEta(xDiTau) = -1234.;
300  decTruthLeadPhi(xDiTau) = -1234.;
301  decTruthLeadM(xDiTau) = -1234.;
302  decTruthLeadPdgID(xDiTau) = -1234.;
303  decTruthSubleadPt(xDiTau) = -1234.;
304  decTruthSubleadEta(xDiTau) = -1234.;
305  decTruthSubleadPhi(xDiTau) = -1234.;
306  decTruthSubleadM(xDiTau) = -1234.;
307  decTruthSubleadPdgID(xDiTau) = -1234.;
308  decTruthDeltaR(xDiTau) = -1234.;
309  decTruthMass(xDiTau) = -1234.;
310  }
311 
312  return StatusCode::SUCCESS;
313 }

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

◆ endInputFile()

StatusCode asg::AsgMetadataTool::endInputFile ( )
protectedvirtualinherited

Function called when the currently open input file got completely processed.

Dummy implementation that can be overridden by the derived tool.

Reimplemented in BookkeeperTool, BookkeeperDumperTool, xAODMaker::TriggerMenuMetaDataTool, and xAODMaker::TruthMetaDataTool.

Definition at line 193 of file AsgMetadataTool.cxx.

193  {
194 
195  // Return gracefully:
196  return StatusCode::SUCCESS;
197  }

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

◆ examineTruthTau()

StatusCode BuildTruthTaus::examineTruthTau ( const xAOD::TruthParticle xTruthParticle) const
privateinherited

Definition at line 216 of file BuildTruthTaus.cxx.

217 {
218  // skip this tau if it has no decay vertex, should not happen
219  if ( !xTruthParticle.hasDecayVtx() )
220  return StatusCode::FAILURE;
221 
222  ATH_MSG_VERBOSE("looking for charged daughters of a truth tau");
223 
224  TauTruthInfo truthInfo;
225 
226  const xAOD::TruthVertex* xDecayVertex = xTruthParticle.decayVtx();
227  if (xDecayVertex == nullptr)
228  return StatusCode::FAILURE;
229  for ( size_t iOutgoingParticle = 0; iOutgoingParticle < xDecayVertex->nOutgoingParticles(); ++iOutgoingParticle )
230  {
231  const xAOD::TruthParticle* xTruthDaughter = xDecayVertex->outgoingParticle(iOutgoingParticle);
232  if (xTruthDaughter == nullptr)
233  {
234  ATH_MSG_ERROR("Truth daughter of tau decay was not found in "<< m_truthParticleContainer.key() <<" container. Please ensure that this container has the full tau decay information or produce the TruthTaus container in AtlasDerivation.");
235  return StatusCode::FAILURE;
236  }
237 
238  // if tau decays into tau this is not a proper tau decay
239  if ( xTruthDaughter->isTau() )
240  {
241  ATH_MSG_VERBOSE("Tau decays into a tau itself -> skip this decay");
242  return StatusCode::FAILURE;
243  }
244  }
245 
246  examineTruthTauDecay(xTruthParticle, truthInfo).ignore();
247 
248  if (truthInfo.m_bIsHadronicTau)
249  ATH_MSG_VERBOSE(truthInfo.m_iNChargedDaughters << " prong hadronic truth tau was found with uniqueID "<<HepMC::uniqueID(xTruthParticle));
250  else
251  ATH_MSG_VERBOSE(truthInfo.m_iNChargedDaughters << " prong leptonic truth tau was found with uniqueID "<<HepMC::uniqueID(xTruthParticle));
252  if ( truthInfo.m_iNChargedDaughters%2 == 0 )
253  {
254  ATH_MSG_WARNING("found tau with even multiplicity: " << truthInfo.m_iNChargedDaughters);
255  printDecay(xTruthParticle);
256  }
257 
258  static const SG::Decorator<double> decPtVis("pt_vis");
259  static const SG::Decorator<double> decEtaVis("eta_vis");
260  static const SG::Decorator<double> decPhiVis("phi_vis");
261  static const SG::Decorator<double> decMVis("m_vis");
262 
263  static const SG::Decorator<size_t> decNumCharged("numCharged");
264  static const SG::Decorator<size_t> decNumChargedPion("numChargedPion");
265  static const SG::Decorator<size_t> decNumNeutral("numNeutral");
266  static const SG::Decorator<size_t> decNumNeutralPion("numNeutralPion");
267 
268  decPtVis(xTruthParticle) = truthInfo.m_vTruthVisTLV.Pt();
269  decEtaVis(xTruthParticle) = truthInfo.m_vTruthVisTLV.Eta();
270  decPhiVis(xTruthParticle) = truthInfo.m_vTruthVisTLV.Phi();
271  decMVis(xTruthParticle) = truthInfo.m_vTruthVisTLV.M();
272 
273  decNumCharged(xTruthParticle) = truthInfo.m_iNChargedDaughters;
274  decNumChargedPion(xTruthParticle) = truthInfo.m_iNChargedPions;
275  decNumNeutral(xTruthParticle) = truthInfo.m_iNNeutralPions+truthInfo.m_iNNeutralOthers;
276  decNumNeutralPion(xTruthParticle) = truthInfo.m_iNNeutralPions;
277 
278  static const SG::Decorator<char> decIsHadronicTau("IsHadronicTau");
279  decIsHadronicTau(xTruthParticle) = (char)truthInfo.m_bIsHadronicTau;
280 
282  {
283  TLorentzVector vTruthInvisTLV = xTruthParticle.p4() - truthInfo.m_vTruthVisTLV;
284  static const SG::Decorator<double> decPtInvis("pt_invis");
285  static const SG::Decorator<double> decEtaInvis("eta_invis");
286  static const SG::Decorator<double> decPhiInvis("phi_invis");
287  static const SG::Decorator<double> decMInvis("m_invis");
288  decPtInvis(xTruthParticle) = vTruthInvisTLV.Pt();
289  decEtaInvis(xTruthParticle) = vTruthInvisTLV.Eta();
290  decPhiInvis(xTruthParticle) = vTruthInvisTLV.Phi();
291  decMInvis(xTruthParticle) = vTruthInvisTLV.M();
292  }
293 
295  {
296  static const SG::Decorator<double> decPtVisCharged("pt_vis_charged");
297  static const SG::Decorator<double> decEtaVisCharged("eta_vis_charged");
298  static const SG::Decorator<double> decPhiVisCharged("phi_vis_charged");
299  static const SG::Decorator<double> decMVisCharged("m_vis_charged");
300  decPtVisCharged(xTruthParticle) = truthInfo.m_vTruthVisTLVCharged.Pt();
301  decEtaVisCharged(xTruthParticle) = truthInfo.m_vTruthVisTLVCharged.Eta();
302  decPhiVisCharged(xTruthParticle) = truthInfo.m_vTruthVisTLVCharged.Phi();
303  decMVisCharged(xTruthParticle) = truthInfo.m_vTruthVisTLVCharged.M();
304  }
305 
307  {
308  static const SG::Decorator<double> decPtVisNeutral("pt_vis_neutral");
309  static const SG::Decorator<double> decEtaVisNeutral("eta_vis_neutral");
310  static const SG::Decorator<double> decPhiVisNeutral("phi_vis_neutral");
311  static const SG::Decorator<double> decMVisNeutral("m_vis_neutral");
312  decPtVisNeutral(xTruthParticle) = truthInfo.m_vTruthVisTLVNeutral.Pt();
313  decEtaVisNeutral(xTruthParticle) = truthInfo.m_vTruthVisTLVNeutral.Eta();
314  decPhiVisNeutral(xTruthParticle) = truthInfo.m_vTruthVisTLVNeutral.Phi();
315  decMVisNeutral(xTruthParticle) = truthInfo.m_vTruthVisTLVNeutral.M();
316  }
317 
319  {
320  static const SG::Decorator<std::vector<int> > decDecayModeVector("DecayModeVector");
321  decDecayModeVector(xTruthParticle) = truthInfo.m_vDecayMode;
322  }
323 
324  if ( m_bWriteVertices )
325  {
326  // tau decay vertex
327  static const SG::Decorator<float> decDecayVertexX("decay_vertex_x");
328  static const SG::Decorator<float> decDecayVertexY("decay_vertex_y");
329  static const SG::Decorator<float> decDecayVertexZ("decay_vertex_z");
330 
331  decDecayVertexX(xTruthParticle) = truthInfo.m_vDecayVertex.X();
332  decDecayVertexY(xTruthParticle) = truthInfo.m_vDecayVertex.Y();
333  decDecayVertexZ(xTruthParticle) = truthInfo.m_vDecayVertex.Z();
334 
335  // tau production vertex
336  static const SG::Decorator<float> decProdVertexX("prod_vertex_x");
337  static const SG::Decorator<float> decProdVertexY("prod_vertex_y");
338  static const SG::Decorator<float> decProdVertexZ("prod_vertex_z");
339 
340  decProdVertexX(xTruthParticle) = truthInfo.m_vProdVertex.X();
341  decProdVertexY(xTruthParticle) = truthInfo.m_vProdVertex.Y();
342  decProdVertexZ(xTruthParticle) = truthInfo.m_vProdVertex.Z();
343  }
344 
345  return StatusCode::SUCCESS;
346 }

◆ examineTruthTauDecay()

StatusCode BuildTruthTaus::examineTruthTauDecay ( const xAOD::TruthParticle xTruthParticle,
TauTruthInfo truthInfo 
) const
privateinherited

Definition at line 349 of file BuildTruthTaus.cxx.

351 {
352  // get vertex and check if it exists
353  const xAOD::TruthVertex* xDecayVertex = xTruthParticle.decayVtx();
354  if (!xDecayVertex)
355  return StatusCode::SUCCESS;
356 
357  truthInfo.m_vDecayVertex.SetXYZ(xDecayVertex->x(),xDecayVertex->y(),xDecayVertex->z());
358 
359  if (xTruthParticle.hasProdVtx() ) {
360  const xAOD::TruthVertex* xProdVertex = xTruthParticle.prodVtx();
361  truthInfo.m_vProdVertex.SetXYZ(xProdVertex->x(),xProdVertex->y(),xProdVertex->z());
362  } else {
363  truthInfo.m_vProdVertex.SetXYZ(-1234,-1234,-1234);
364  }
365 
366  for ( size_t iOutgoingParticle = 0; iOutgoingParticle < xDecayVertex->nOutgoingParticles(); ++iOutgoingParticle )
367  {
368  const xAOD::TruthParticle* xTruthDaughter = xDecayVertex->outgoingParticle(iOutgoingParticle);
369  if (xTruthDaughter == nullptr)
370  {
371  ATH_MSG_ERROR("Truth daughter of tau decay was not found in "<< m_truthParticleContainer.key() <<" container. Please ensure that this container has the full tau decay information or produce the TruthTaus container in AtlasDerivation.");
372  return StatusCode::FAILURE;
373  }
374 
375  int iAbsPdgId = xTruthDaughter->absPdgId();
376  int iPdgId = xTruthDaughter->pdgId();
377 
378  // look at decay of unstable particles
379  if (MC::isDecayed(xTruthDaughter) || !MC::isPhysical(xTruthDaughter))
380  {
381  if ( iAbsPdgId != 111 && iAbsPdgId != 311 && iAbsPdgId != 310 && iAbsPdgId != 130 )
382  {
383  examineTruthTauDecay(*xTruthDaughter, truthInfo).ignore();
384  continue;
385  }
386  }
387 
388  // only process stable particles
389  if (!MC::isStable(xTruthDaughter) && !MC::isDecayed(xTruthDaughter))
390  continue;
391 
392  // add pdgID to vector for decay mode classification
393  truthInfo.m_vDecayMode.push_back(iPdgId);
394 
395  // if tau decays leptonically, indicated by an electron/muon neutrino then
396  // it is not a hadronic decay
397  if ( xTruthDaughter->isHadron() )
398  truthInfo.m_bIsHadronicTau = true;
399 
400  // ignore neutrinos for further progress
401  if ( xTruthDaughter->isNeutrino() )
402  {
403  ATH_MSG_VERBOSE("found neutrino decay particle with PdgId "<<iPdgId);
404  continue;
405  }
406 
407  // add momentum of non-neutrino particle to visible momentum
408  truthInfo.m_vTruthVisTLV += xTruthDaughter->p4();
410  if ( xTruthDaughter->isCharged() )
411  truthInfo.m_vTruthVisTLVCharged += xTruthDaughter->p4();
413  if ( xTruthDaughter->isNeutral() )
414  truthInfo.m_vTruthVisTLVNeutral += xTruthDaughter->p4();
415 
416  // only count charged decay particles
417  if ( xTruthDaughter->isCharged() )
418  {
419  ATH_MSG_VERBOSE("found charged decay particle with PdgId "<<iPdgId);
420  truthInfo.m_iNChargedDaughters++;
421  // count charged pions
422  if (iAbsPdgId==211) truthInfo.m_iNChargedPions++;
423  else truthInfo.m_iNChargedOthers++;
424  }
425  else
426  {
427  ATH_MSG_VERBOSE("found neutral decay particle with PdgId "<<iPdgId);
428  // count neutral pions
429  if (iAbsPdgId==111) truthInfo.m_iNNeutralPions++;
430  else truthInfo.m_iNNeutralOthers++;
431  }
432  }
433  return StatusCode::SUCCESS;
434 }

◆ 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

◆ findTruthTau()

StatusCode DiTauTruthMatchingTool::findTruthTau ( const xAOD::DiTauJet xDiTau)
private

Definition at line 82 of file DiTauTruthMatchingTool.cxx.

83 {
84  // check if decorations were already added to the first passed tau
86  static const SG::ConstAccessor<char> accIsTruthMatched("IsTruthMatched");
87  m_bIsTruthMatchedAvailable.set (accIsTruthMatched.isAvailable(xDiTau));
88  }
90  return StatusCode::SUCCESS;
91 
93  return checkTruthMatch(xDiTau, *m_truthTausEvent.m_xTruthTauContainerConst);
94  else
95  return checkTruthMatch(xDiTau, *m_truthTausEvent.m_xTruthTauContainer);
96 }

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

◆ getTruth() [1/2]

void DiTauTruthMatchingTool::getTruth ( const std::vector< const xAOD::DiTauJet * > &  vDiTaus)
privatevirtual

Implements TauAnalysisTools::IDiTauTruthMatchingTool.

Definition at line 69 of file DiTauTruthMatchingTool.cxx.

70 {
71  for (auto xDiTau : vDiTaus)
72  getTruth(*xDiTau);
73  return;
74 }

◆ getTruth() [2/2]

void DiTauTruthMatchingTool::getTruth ( const xAOD::DiTauJet xDiTau)
privatevirtual

get pointer to the truth matched particle.

if no truth particle was found a null pointer is returned

Implements TauAnalysisTools::IDiTauTruthMatchingTool.

Definition at line 57 of file DiTauTruthMatchingTool.cxx.

58 {
59  if (retrieveTruthTaus().isFailure())
60  return;
61 
62  if (findTruthTau(xDiTau).isFailure())
63  ATH_MSG_WARNING("There was a failure in finding the matched truth tau");
64 
65  return;
66 }

◆ handle()

void asg::AsgMetadataTool::handle ( const Incident &  inc)
protectedvirtualinherited

Function receiving incidents from IncidentSvc/TEvent.

Reimplemented in Trig::TrigDecisionTool.

Definition at line 135 of file AsgMetadataTool.cxx.

135  {
136 
137  // Tell the user what's happening:
138  ATH_MSG_VERBOSE( "Callback received with incident: " << inc.type() );
139 
140  // Call the appropriate member function:
141  if( inc.type() == IncidentType::BeginInputFile ) {
142  m_beginInputFileCalled = true;
143  if( beginInputFile().isFailure() ) {
144  ATH_MSG_FATAL( "Failed to call beginInputFile()" );
145  throw std::runtime_error( "Couldn't call beginInputFile()" );
146  }
147  } else if( inc.type() == IncidentType::EndInputFile ) {
148  if( endInputFile().isFailure() ) {
149  ATH_MSG_FATAL( "Failed to call endInputFile()" );
150  throw std::runtime_error( "Couldn't call endInputFile()" );
151  }
152  } else if( inc.type() == IncidentType::BeginEvent ) {
153  // If the tool didn't catch the begin input file incident for the
154  // first input file of the job, then call the appropriate function
155  // now.
156  if( ! m_beginInputFileCalled ) {
157  m_beginInputFileCalled = true;
158  if( beginInputFile().isFailure() ) {
159  ATH_MSG_FATAL( "Failed to call beginInputFile()" );
160  throw std::runtime_error( "Couldn't call beginInputFile()" );
161  }
162  }
163  if( beginEvent().isFailure() ) {
164  ATH_MSG_FATAL( "Failed to call beginEvent()" );
165  throw std::runtime_error( "Couldn't call beginEvent()" );
166  }
167 
168  #ifdef XAOD_STANDALONE
169  } else if( inc.type() == IncidentType::MetaDataStop ) {
170  if( metaDataStop().isFailure() ) {
171  ATH_MSG_FATAL( "Failed to call metaDataStop()" );
172  throw std::runtime_error( "Couldn't call metaDataStop()" );
173  }
174 
175  #endif // XAOD_STANDALONE
176  } else {
177  ATH_MSG_WARNING( "Unknown incident type received in AsgMetaDataTool: " << inc.type() );
178  }
179 
180  return;
181  }

◆ initialize()

StatusCode DiTauTruthMatchingTool::initialize ( )
privatevirtual

Reimplemented from TauAnalysisTools::BuildTruthTaus.

Definition at line 41 of file DiTauTruthMatchingTool.cxx.

42 {
43  ATH_MSG_INFO( "Initializing DiTauTruthMatchingTool" );
44 
45  // configure BuildTruthTaus in truth matching mode, not truth tau building mode
46  DiTauTruthMatchingTool::BuildTruthTaus::setTruthMatchingMode();
47 
49  {
50  ATH_MSG_FATAL("Failed initializing BuildTruthTaus");
51  return StatusCode::FAILURE;
52  }
53  return StatusCode::SUCCESS;
54 }

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

◆ inputMetaStore()

AsgMetadataTool::MetaStorePtr_t asg::AsgMetadataTool::inputMetaStore ( ) const
inherited

Accessor for the input metadata store.

Definition at line 93 of file AsgMetadataTool.cxx.

93  {
94 
95 #ifdef XAOD_STANDALONE
96  return &m_inputMetaStore;
97 #else // XAOD_STANDALONE
98  return m_inputMetaStore;
99 #endif // XAOD_STANDALONE
100  }

◆ metaDataStop()

StatusCode asg::AsgMetadataTool::metaDataStop ( )
protectedvirtualinherited

Function called when the tool should write out its metadata.

Dummy implementation that can be overridden by the derived tool.

Reimplemented in BookkeeperTool, BookkeeperDumperTool, xAODMaker::TriggerMenuMetaDataTool, and xAODMaker::TruthMetaDataTool.

Definition at line 209 of file AsgMetadataTool.cxx.

209  {
210 
211  // Return gracefully:
212  return StatusCode::SUCCESS;
213  }

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

◆ outputMetaStore()

AsgMetadataTool::MetaStorePtr_t asg::AsgMetadataTool::outputMetaStore ( ) const
inherited

Accessor for the output metadata store.

Definition at line 102 of file AsgMetadataTool.cxx.

102  {
103 
104 #ifdef XAOD_STANDALONE
105  return &m_outputMetaStore;
106 #else // XAOD_STANDALONE
107  return m_outputMetaStore;
108 #endif // XAOD_STANDALONE
109  }

◆ print() [1/2]

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

◆ print() [2/2]

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

◆ printDecay()

void BuildTruthTaus::printDecay ( const xAOD::TruthParticle xTruthParticle,
int  depth = 0 
) const
privateinherited

Definition at line 436 of file BuildTruthTaus.cxx.

437 {
438  // loop over all decay particles, print their kinematic and other properties
439 
440  const xAOD::TruthVertex* xDecayVertex = xTruthParticle.decayVtx();
441  if (xDecayVertex == nullptr)
442  return;
443 
444  for ( size_t iOutgoingParticle = 0; iOutgoingParticle < xDecayVertex->nOutgoingParticles(); ++iOutgoingParticle )
445  {
446  const xAOD::TruthParticle* xTruthDaughter = xDecayVertex->outgoingParticle(iOutgoingParticle);
447  if (xTruthDaughter == nullptr)
448  {
449  ATH_MSG_WARNING("Truth daughter of tau decay was not found in "<< m_truthParticleContainer.key() <<" container. Please ensure that this container has the full tau decay information or produce the TruthTaus container in AtlasDerivation.");
450  return;
451  }
452  ATH_MSG_WARNING("depth "<<depth
453  <<" e "<<xTruthDaughter->e()
454  <<" eta "<<xTruthDaughter->p4().Eta()
455  <<" phi "<<xTruthDaughter->p4().Phi()
456  <<" pdgid "<<xTruthDaughter->pdgId()
457  <<" status "<<xTruthDaughter->status()
458  <<" uniqueID "<<HepMC::uniqueID(xTruthDaughter));
459  printDecay(*xTruthDaughter, depth+1);
460  }
461 }

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

◆ retrieveTruthTaus() [1/3]

StatusCode BuildTruthTaus::retrieveTruthTaus ( )
overrideprivatevirtualinherited

Implements TauAnalysisTools::IBuildTruthTaus.

Definition at line 82 of file BuildTruthTaus.cxx.

83 {
85 }

◆ retrieveTruthTaus() [2/3]

StatusCode BuildTruthTaus::retrieveTruthTaus ( ITruthTausEvent truthTausEvent) const
overrideprivatevirtualinherited

Implements TauAnalysisTools::IBuildTruthTaus.

Definition at line 88 of file BuildTruthTaus.cxx.

89 {
90  return retrieveTruthTaus (dynamic_cast<TruthTausEvent&> (truthTausEvent));
91 }

◆ retrieveTruthTaus() [3/3]

StatusCode BuildTruthTaus::retrieveTruthTaus ( TruthTausEvent &  truthTausEvent) const
protectedinherited

Definition at line 94 of file BuildTruthTaus.cxx.

95 {
96  const EventContext& ctx = Gaudi::Hive::currentContext();
97 
98  // truth matching mode
99  if (m_truthMatchingMode) {
102  if (!truthElectronsHandle.isValid()) {
103  ATH_MSG_ERROR ("Could not retrieve " << truthElectronsHandle.key());
104  return StatusCode::FAILURE;
105  }
106  truthTausEvent.m_xTruthElectronContainerConst = truthElectronsHandle.cptr();
107  }
108 
109  if (!m_truthMuonContainer.empty()) {
111  if (!truthMuonsHandle.isValid()) {
112  ATH_MSG_ERROR ("Could not retrieve " << truthMuonsHandle.key());
113  return StatusCode::FAILURE;
114  }
115  truthTausEvent.m_xTruthMuonContainerConst = truthMuonsHandle.cptr();
116  }
117 
118  if (!m_truthJetContainer.empty()) {
120  if (!truthJetsHandle.isValid()) {
121  ATH_MSG_ERROR ("Could not retrieve " << truthJetsHandle.key());
122  return StatusCode::FAILURE;
123  }
124  truthTausEvent.m_xTruthJetContainerConst = truthJetsHandle.cptr();
125  }
126 
128  if (!truthTausHandle.isValid()) {
129  ATH_MSG_ERROR ("Could not retrieve " << truthTausHandle.key());
130  return StatusCode::FAILURE;
131  }
132  truthTausEvent.m_xTruthTauContainerConst = truthTausHandle.cptr();
133  }
134  // truth tau building mode
135  else {
137  if (!truthParticlesHandle.isValid()) {
138  ATH_MSG_ERROR ("Could not retrieve " << truthParticlesHandle.key());
139  return StatusCode::FAILURE;
140  }
141  truthTausEvent.m_xTruthParticleContainer = truthParticlesHandle.cptr();
142 
143  auto truthTausOutput = std::make_unique<xAOD::TruthParticleContainer>();
144  auto truthTausOutputAux = std::make_unique<xAOD::TruthParticleAuxContainer>();
145  truthTausOutput->setStore(truthTausOutputAux.get());
146  truthTausEvent.m_xTruthTauContainer = truthTausOutput.get();
147 
148  ATH_CHECK( buildTruthTausFromTruthParticles(truthTausEvent) );
149 
150  auto writeHandle = SG::makeHandle(m_truthTauOutputContainer, ctx);
151  ATH_CHECK(writeHandle.record(std::move(truthTausOutput), std::move(truthTausOutputAux)));
152  }
153 
154  return StatusCode::SUCCESS;
155 }

◆ setTruthMatchingMode()

virtual void TauAnalysisTools::BuildTruthTaus::setTruthMatchingMode ( )
inlineoverrideprivatevirtualinherited

Implements TauAnalysisTools::IBuildTruthTaus.

Definition at line 61 of file BuildTruthTaus.h.

61  {
62  m_truthMatchingMode = true;
63  };

◆ setUseIncidents()

void asg::AsgMetadataTool::setUseIncidents ( const bool  flag)
inlineprotectedinherited

Definition at line 132 of file AsgMetadataTool.h.

133  {
135  }

◆ sysInitialize()

StatusCode asg::AsgMetadataTool::sysInitialize ( )
virtualinherited

Function initialising the tool in the correct way in Athena.

This function is used to set up the callbacks from IncidentSvc in Athena at the right time during initialisation, without the user having to do anything special in his/her code.

Reimplemented from AthCommonDataStore< AthCommonMsg< AlgTool > >.

Definition at line 115 of file AsgMetadataTool.cxx.

115  {
116 
117 #ifndef XAOD_STANDALONE
118  if (m_useIncidents) {
119  // Connect to the IncidentSvc:
120  ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
121  ATH_CHECK( incSvc.retrieve() );
122 
123  // Set up the right callbacks: don't rethrow exceptions, any failure and we should end
124  incSvc->addListener( this, IncidentType::BeginEvent, 0, false );
125  }
126  // Let the base class do its thing:
127  ATH_CHECK( AlgTool::sysInitialize() );
128 
129 #endif // not XAOD_STANDALONE
130 
131  // Return gracefully:
132  return StatusCode::SUCCESS;
133  }

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

◆ truthMatch()

StatusCode DiTauTruthMatchingTool::truthMatch ( const TLorentzVector &  vSubjetTLV,
const xAOD::TruthParticleContainer xTruthTauContainer,
const xAOD::TruthParticle *&  xTruthMatch,
const xAOD::Jet *&  xTruthJetMatch,
TruthMatchedParticleType eTruthMatchedParticleType 
) const
private

Definition at line 315 of file DiTauTruthMatchingTool.cxx.

320 {
321  for (auto xTruthTauIt : xTruthTauContainer)
322  {
323  TLorentzVector vTruthVisTLV;
324  vTruthVisTLV.SetPtEtaPhiM(m_accPtVis(*xTruthTauIt),
325  m_accEtaVis(*xTruthTauIt),
326  m_accPhiVis(*xTruthTauIt),
327  m_accMVis(*xTruthTauIt));
328  if (vSubjetTLV.DeltaR(vTruthVisTLV) <= m_dMaxDeltaR)
329  {
330  static const SG::ConstAccessor<char> accIsHadronicTau("IsHadronicTau");
331  if (static_cast<bool>(accIsHadronicTau(*xTruthTauIt)))
332  eTruthMatchedParticleType = TruthHadronicTau;
333  else
334  continue; // don't let leptonic taus steal truthmatch just by chance
335 
336  xTruthMatch = xTruthTauIt;
337  break;
338  }
339  }
340 
341  if (!xTruthMatch and m_truthTausEvent.m_xTruthMuonContainerConst)
342  {
343  double dPtMax = 0.;
344  for (auto xTruthMuonIt : *m_truthTausEvent.m_xTruthMuonContainerConst)
345  {
346  if (vSubjetTLV.DeltaR(xTruthMuonIt->p4()) <= m_dMaxDeltaR)
347  {
348  if (xTruthMuonIt->pt()<dPtMax)
349  continue;
350  eTruthMatchedParticleType = TruthMuon;
351 
352  xTruthMatch = xTruthMuonIt;
353  dPtMax = xTruthMuonIt->pt();
354  }
355  }
356  }
357 
358  if (!xTruthMatch and m_truthTausEvent.m_xTruthElectronContainerConst)
359  {
360  double dPtMax = 0.;
361  for (auto xTruthElectronIt : *m_truthTausEvent.m_xTruthElectronContainerConst)
362  {
363  if (vSubjetTLV.DeltaR(xTruthElectronIt->p4()) <= m_dMaxDeltaR)
364  {
365  if (xTruthElectronIt->pt()<dPtMax)
366  continue;
367  eTruthMatchedParticleType = TruthElectron;
368  xTruthMatch = xTruthElectronIt;
369  dPtMax = xTruthElectronIt->pt();
370  }
371  }
372  }
373 
374  if (m_truthTausEvent.m_xTruthJetContainerConst)
375  {
376  double dPtMax = 0.;
377  for (auto xTruthJetIt : *m_truthTausEvent.m_xTruthJetContainerConst)
378  {
379  if (vSubjetTLV.DeltaR(xTruthJetIt->p4()) <= m_dMaxDeltaR)
380  {
381  if (xTruthJetIt->pt()<dPtMax)
382  continue;
383  xTruthJetMatch = xTruthJetIt;
384  dPtMax = xTruthJetIt->pt();
385  }
386  }
387  }
388 
389  return StatusCode::SUCCESS;
390 }

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

SG::ConstAccessor<double> TauAnalysisTools::DiTauTruthMatchingTool::m_accEtaVis
private

Definition at line 71 of file DiTauTruthMatchingTool.h.

◆ m_accMVis

SG::ConstAccessor<double> TauAnalysisTools::DiTauTruthMatchingTool::m_accMVis
private

Definition at line 73 of file DiTauTruthMatchingTool.h.

◆ m_accPhiVis

SG::ConstAccessor<double> TauAnalysisTools::DiTauTruthMatchingTool::m_accPhiVis
private

Definition at line 72 of file DiTauTruthMatchingTool.h.

◆ m_accPtVis

SG::ConstAccessor<double> TauAnalysisTools::DiTauTruthMatchingTool::m_accPtVis
private

Definition at line 70 of file DiTauTruthMatchingTool.h.

◆ m_beginInputFileCalled

bool asg::AsgMetadataTool::m_beginInputFileCalled
privateinherited

Flag helping to discover when the tool misses the opening of the first input file.

Definition at line 126 of file AsgMetadataTool.h.

◆ m_bIsTruthMatchedAvailable

CxxUtils::CachedValue<bool> TauAnalysisTools::DiTauTruthMatchingTool::m_bIsTruthMatchedAvailable
private

Definition at line 65 of file DiTauTruthMatchingTool.h.

◆ m_bIsTruthParticleLinkAvailable

CxxUtils::CachedValue<bool> TauAnalysisTools::DiTauTruthMatchingTool::m_bIsTruthParticleLinkAvailable
private

Definition at line 66 of file DiTauTruthMatchingTool.h.

◆ m_bTruthTauAvailable

bool TauAnalysisTools::BuildTruthTaus::m_bTruthTauAvailable {}
protectedinherited

Definition at line 109 of file BuildTruthTaus.h.

◆ m_bWriteDecayModeVector

Gaudi::Property<bool> TauAnalysisTools::BuildTruthTaus::m_bWriteDecayModeVector { this, "WriteDecayModeVector", true}
privateinherited

Definition at line 117 of file BuildTruthTaus.h.

◆ m_bWriteInvisibleFourMomentum

Gaudi::Property<bool> TauAnalysisTools::BuildTruthTaus::m_bWriteInvisibleFourMomentum { this, "WriteInvisibleFourMomentum", false}
privateinherited

Definition at line 114 of file BuildTruthTaus.h.

◆ m_bWriteVertices

Gaudi::Property<bool> TauAnalysisTools::BuildTruthTaus::m_bWriteVertices { this, "WriteVertices", true}
privateinherited

Definition at line 118 of file BuildTruthTaus.h.

◆ m_bWriteVisibleChargedFourMomentum

Gaudi::Property<bool> TauAnalysisTools::BuildTruthTaus::m_bWriteVisibleChargedFourMomentum { this, "WriteVisibleChargedFourMomentum", false}
privateinherited

Definition at line 115 of file BuildTruthTaus.h.

◆ m_bWriteVisibleNeutralFourMomentum

Gaudi::Property<bool> TauAnalysisTools::BuildTruthTaus::m_bWriteVisibleNeutralFourMomentum { this, "WriteVisibleNeutralFourMomentum", false}
privateinherited

Definition at line 116 of file BuildTruthTaus.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_dMaxDeltaR

Gaudi::Property<double> TauAnalysisTools::DiTauTruthMatchingTool::m_dMaxDeltaR {this, "MaxDeltaR", 0.2}
private

Definition at line 63 of file DiTauTruthMatchingTool.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_inputMetaStore

MetaStore_t asg::AsgMetadataTool::m_inputMetaStore
privateinherited

Object accessing the input metadata store.

Definition at line 119 of file AsgMetadataTool.h.

◆ m_outputMetaStore

MetaStore_t asg::AsgMetadataTool::m_outputMetaStore
privateinherited

Object accessing the output metadata store.

Definition at line 121 of file AsgMetadataTool.h.

◆ m_tMCTruthClassifier

asg::AnaToolHandle<IMCTruthClassifier> TauAnalysisTools::BuildTruthTaus::m_tMCTruthClassifier
privateinherited

Definition at line 131 of file BuildTruthTaus.h.

◆ m_truthElectronContainer

SG::ReadHandleKey<xAOD::TruthParticleContainer> TauAnalysisTools::BuildTruthTaus::m_truthElectronContainer { this, "TruthElectronContainerName", "TruthElectrons", "Truth electrons input container name" }
privateinherited

Definition at line 123 of file BuildTruthTaus.h.

◆ m_truthJetContainer

SG::ReadHandleKey<xAOD::JetContainer> TauAnalysisTools::BuildTruthTaus::m_truthJetContainer { this, "TruthJetContainerName", "AntiKt4TruthDressedWZJets", "Truth jets input container name" }
privateinherited

Definition at line 125 of file BuildTruthTaus.h.

◆ m_truthMatchingMode

bool TauAnalysisTools::BuildTruthTaus::m_truthMatchingMode = false
privateinherited

Definition at line 129 of file BuildTruthTaus.h.

◆ m_truthMuonContainer

SG::ReadHandleKey<xAOD::TruthParticleContainer> TauAnalysisTools::BuildTruthTaus::m_truthMuonContainer { this, "TruthMuonContainerName", "TruthMuons", "Truth muons input container name" }
privateinherited

Definition at line 124 of file BuildTruthTaus.h.

◆ m_truthParticleContainer

SG::ReadHandleKey<xAOD::TruthParticleContainer> TauAnalysisTools::BuildTruthTaus::m_truthParticleContainer { this, "TruthParticleContainerName", "TruthParticles", "Truth particles input container name" }
privateinherited

Definition at line 122 of file BuildTruthTaus.h.

◆ m_truthTauInputContainer

SG::ReadHandleKey<xAOD::TruthParticleContainer> TauAnalysisTools::BuildTruthTaus::m_truthTauInputContainer { this, "TruthTauContainerName", "TruthTaus", "Truth tau input container name (truth matching mode)" }
privateinherited

Definition at line 121 of file BuildTruthTaus.h.

◆ m_truthTauOutputContainer

SG::WriteHandleKey<xAOD::TruthParticleContainer> TauAnalysisTools::BuildTruthTaus::m_truthTauOutputContainer { this, "NewTruthTauContainerName", "TruthTaus", "Truth tau output container name" }
privateinherited

Definition at line 127 of file BuildTruthTaus.h.

◆ m_truthTausEvent

TruthTausEvent TauAnalysisTools::BuildTruthTaus::m_truthTausEvent
protectedinherited

Definition at line 106 of file BuildTruthTaus.h.

◆ m_useIncidents

bool asg::AsgMetadataTool::m_useIncidents
privateinherited

Definition at line 128 of file AsgMetadataTool.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:
TauAnalysisTools::BuildTruthTaus::examineTruthTau
StatusCode examineTruthTau(const xAOD::TruthParticle &xTruthParticle) const
Definition: BuildTruthTaus.cxx:216
TauAnalysisTools::BuildTruthTaus::m_bWriteInvisibleFourMomentum
Gaudi::Property< bool > m_bWriteInvisibleFourMomentum
Definition: BuildTruthTaus.h:114
xAOD::TruthVertex_v1::nOutgoingParticles
size_t nOutgoingParticles() const
Get the number of outgoing particles.
TauAnalysisTools::DiTauTruthMatchingTool::getTruth
virtual void getTruth(const xAOD::DiTauJet &xDiTau)
get pointer to the truth matched particle.
Definition: DiTauTruthMatchingTool.cxx:57
TauAnalysisTools::TruthElectron
@ TruthElectron
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:102
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
egammaParameters::depth
@ depth
pointing depth of the shower as calculated in egammaqgcld
Definition: egammaParamDefs.h:276
xAOD::TruthParticle_v1::absPdgId
int absPdgId() const
Absolute PDG ID code (often useful)
CxxUtils::CachedValue::ptr
const T * ptr() const
Return a pointer to the cached value.
asg::AsgMetadataTool::beginEvent
virtual StatusCode beginEvent()
Function called when a new events is loaded.
Definition: AsgMetadataTool.cxx:201
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:407
TauAnalysisTools::BuildTruthTaus::m_truthElectronContainer
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthElectronContainer
Definition: BuildTruthTaus.h:123
TauAnalysisTools::BuildTruthTaus::m_tMCTruthClassifier
asg::AnaToolHandle< IMCTruthClassifier > m_tMCTruthClassifier
Definition: BuildTruthTaus.h:131
TauAnalysisTools::DiTauTruthMatchingTool::m_accMVis
SG::ConstAccessor< double > m_accMVis
Definition: DiTauTruthMatchingTool.h:73
TauAnalysisTools::BuildTruthTaus::m_bWriteVertices
Gaudi::Property< bool > m_bWriteVertices
Definition: BuildTruthTaus.h:118
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
TauAnalysisTools::DiTauTruthMatchingTool::m_bIsTruthMatchedAvailable
CxxUtils::CachedValue< bool > m_bIsTruthMatchedAvailable
Definition: DiTauTruthMatchingTool.h:65
TauAnalysisTools::BuildTruthTaus::printDecay
void printDecay(const xAOD::TruthParticle &xTruthParticle, int depth=0) const
Definition: BuildTruthTaus.cxx:436
SG::Accessor< char >
TauAnalysisTools::DiTauTruthMatchingTool::m_accPtVis
SG::ConstAccessor< double > m_accPtVis
Definition: DiTauTruthMatchingTool.h:70
TauAnalysisTools::TruthHadronicTau
@ TruthHadronicTau
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:99
TauAnalysisTools::DiTauTruthMatchingTool::checkTruthMatch
StatusCode checkTruthMatch(const xAOD::DiTauJet &xDiTau, const xAOD::TruthParticleContainer &xTauContainer) const
Definition: DiTauTruthMatchingTool.cxx:99
SG::ReadHandle< xAOD::TruthParticleContainer >
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
CxxUtils::CachedValue::isValid
bool isValid() const
Test to see if the value is valid.
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
initialize
void initialize()
Definition: run_EoverP.cxx:894
TauAnalysisTools::BuildTruthTaus::retrieveTruthTaus
virtual StatusCode retrieveTruthTaus() override
Definition: BuildTruthTaus.cxx:82
xAOD::char
char
Definition: TrigDecision_v1.cxx:38
xAOD::DiTauJet_v1::subjetPhi
float subjetPhi(unsigned int numSubjet) const
Definition: DiTauJet_v1.cxx:111
TauAnalysisTools::BuildTruthTaus::m_truthTauOutputContainer
SG::WriteHandleKey< xAOD::TruthParticleContainer > m_truthTauOutputContainer
Definition: BuildTruthTaus.h:127
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
TauAnalysisTools::BuildTruthTaus::m_truthMatchingMode
bool m_truthMatchingMode
Definition: BuildTruthTaus.h:129
asg::AsgMetadataTool::m_beginInputFileCalled
bool m_beginInputFileCalled
Flag helping to discover when the tool misses the opening of the first input file.
Definition: AsgMetadataTool.h:126
SG::ConstAccessor< int >
TauAnalysisTools::BuildTruthTaus::m_bWriteDecayModeVector
Gaudi::Property< bool > m_bWriteDecayModeVector
Definition: BuildTruthTaus.h:117
TauAnalysisTools::BuildTruthTaus::m_bWriteVisibleChargedFourMomentum
Gaudi::Property< bool > m_bWriteVisibleChargedFourMomentum
Definition: BuildTruthTaus.h:115
asg::AsgMetadataTool::endInputFile
virtual StatusCode endInputFile()
Function called when the currently open input file got completely processed.
Definition: AsgMetadataTool.cxx:193
asg::AsgMetadataTool::m_outputMetaStore
MetaStore_t m_outputMetaStore
Object accessing the output metadata store.
Definition: AsgMetadataTool.h:121
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
TauAnalysisTools::Unknown
@ Unknown
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:98
SG::VarHandleKey::empty
bool empty() const
Test if the key is blank.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:150
TauAnalysisTools::DiTauTruthMatchingTool::m_accEtaVis
SG::ConstAccessor< double > m_accEtaVis
Definition: DiTauTruthMatchingTool.h:71
xAOD::TruthVertex_v1::y
float y() const
Vertex y displacement.
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
xAOD::DiTauJet_v1::subjetE
float subjetE(unsigned int numSubjet) const
Definition: DiTauJet_v1.cxx:121
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
SG::makeHandle
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
Definition: ReadCondHandle.h:274
MC::isPhysical
bool isPhysical(const T &p)
Identify if the particle is physical, i.e. is stable or decayed.
Definition: HepMCHelpers.h:51
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
xAOD::TruthParticle_v1::hasDecayVtx
bool hasDecayVtx() const
Check for a decay vertex on this particle.
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
TauAnalysisTools::DiTauTruthMatchingTool::truthMatch
StatusCode truthMatch(const TLorentzVector &vSubjetTLV, const xAOD::TruthParticleContainer &xTruthTauContainer, const xAOD::TruthParticle *&xTruthMatch, const xAOD::Jet *&xTruthJetMatch, TruthMatchedParticleType &eTruthMatchedParticleType) const
Definition: DiTauTruthMatchingTool.cxx:315
xAOD::TruthParticle_v1::e
virtual double e() const override final
The total energy of the particle.
SG::Decorator< char >
lumiFormat.i
int i
Definition: lumiFormat.py:85
TauAnalysisTools::BuildTruthTaus::m_bTruthTauAvailable
bool m_bTruthTauAvailable
Definition: BuildTruthTaus.h:109
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
xAOD::TruthParticle_v1
Class describing a truth particle in the MC record.
Definition: TruthParticle_v1.h:37
asg::AsgMetadataTool::beginInputFile
virtual StatusCode beginInputFile()
Function called when a new input file is opened.
Definition: AsgMetadataTool.cxx:185
master.flag
bool flag
Definition: master.py:29
xAOD::TruthParticle_v1::hasProdVtx
bool hasProdVtx() const
Check for a production vertex on this particle.
Definition: TruthParticle_v1.cxx:69
CalibDbCompareRT.dummy
dummy
Definition: CalibDbCompareRT.py:59
HepMC::uniqueID
int uniqueID(const T &p)
Definition: MagicNumbers.h:116
asg::AsgMetadataTool::metaDataStop
virtual StatusCode metaDataStop()
Function called when the tool should write out its metadata.
Definition: AsgMetadataTool.cxx:209
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
xAOD::TruthParticle_v1::isTau
bool isTau() const
Whether the particle is a tau (or antitau)
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
TauAnalysisTools::BuildTruthTaus::buildTruthTausFromTruthParticles
StatusCode buildTruthTausFromTruthParticles(TruthTausEvent &truthTausEvent) const
Definition: BuildTruthTaus.cxx:161
xAOD::TruthParticle_v1::isCharged
bool isCharged() const
Whether the particle is electrically charged.
TauAnalysisTools::TruthMuon
@ TruthMuon
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:101
xAOD::DiTauJet_v1::subjetEta
float subjetEta(unsigned int numSubjet) const
Definition: DiTauJet_v1.cxx:101
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
xAOD::TruthParticle_v1::decayVtx
const TruthVertex_v1 * decayVtx() const
The decay vertex of this particle.
xAOD::TruthParticle_v1::prodVtx
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Definition: TruthParticle_v1.cxx:75
xAOD::TruthVertex_v1
Class describing a truth vertex in the MC record.
Definition: TruthVertex_v1.h:37
TauAnalysisTools::DiTauTruthMatchingTool::m_accPhiVis
SG::ConstAccessor< double > m_accPhiVis
Definition: DiTauTruthMatchingTool.h:72
TauAnalysisTools::BuildTruthTaus::m_truthTausEvent
TruthTausEvent m_truthTausEvent
Definition: BuildTruthTaus.h:106
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
asg::AsgMetadataTool::m_inputMetaStore
MetaStore_t m_inputMetaStore
Object accessing the input metadata store.
Definition: AsgMetadataTool.h:119
TauAnalysisTools::TruthMatchedParticleType
TruthMatchedParticleType
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:97
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
TauAnalysisTools::BuildTruthTaus::m_truthTauInputContainer
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthTauInputContainer
Definition: BuildTruthTaus.h:121
runIDPVM.pdgId
pdgId
Definition: runIDPVM.py:91
CxxUtils::CachedValue::set
void set(const T &val) const
Set the value, assuming it is currently invalid.
xAOD::TruthVertex_v1::x
float x() const
Vertex x displacement.
DataVector::end
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
TauAnalysisTools::DiTauTruthMatchingTool::m_dMaxDeltaR
Gaudi::Property< double > m_dMaxDeltaR
Definition: DiTauTruthMatchingTool.h:63
xAOD::Jet_v1
Class describing a jet.
Definition: Jet_v1.h:57
TauAnalysisTools::BuildTruthTaus::m_truthJetContainer
SG::ReadHandleKey< xAOD::JetContainer > m_truthJetContainer
Definition: BuildTruthTaus.h:125
xAOD::TruthParticle_v1::status
int status() const
Status code.
MC::isStable
bool isStable(const T &p)
Identify if the particle is stable, i.e. has not decayed.
Definition: HepMCHelpers.h:45
IMCTruthClassifier::particleTruthClassifier
virtual std::pair< MCTruthPartClassifier::ParticleType, MCTruthPartClassifier::ParticleOrigin > particleTruthClassifier(const xAOD::TruthParticle *, MCTruthPartClassifier::Info *info=nullptr) const =0
a
TList * a
Definition: liststreamerinfos.cxx:10
h
TauAnalysisTools::BuildTruthTaus::m_truthMuonContainer
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthMuonContainer
Definition: BuildTruthTaus.h:124
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
TauAnalysisTools::DiTauTruthMatchingTool::findTruthTau
StatusCode findTruthTau(const xAOD::DiTauJet &xDiTau)
Definition: DiTauTruthMatchingTool.cxx:82
MC::isDecayed
bool isDecayed(const T &p)
Identify if the particle decayed.
Definition: HepMCHelpers.h:42
xAOD::TruthVertex_v1::z
float z() const
Vertex longitudinal distance along the beam line form the origin.
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
if
if(febId1==febId2)
Definition: LArRodBlockPhysicsV0.cxx:567
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
xAOD::TruthParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: TruthParticle_v1.cxx:161
TauAnalysisTools::BuildTruthTaus::m_bWriteVisibleNeutralFourMomentum
Gaudi::Property< bool > m_bWriteVisibleNeutralFourMomentum
Definition: BuildTruthTaus.h:116
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
xAOD::TruthParticle_v1::isNeutrino
bool isNeutrino() const
Whether the particle is a neutrino (or antineutrino)
TauAnalysisTools::BuildTruthTaus::m_truthParticleContainer
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthParticleContainer
Definition: BuildTruthTaus.h:122
xAOD::DiTauJet_v1::subjetPt
float subjetPt(unsigned int numSubjet) const
Definition: DiTauJet_v1.cxx:91
xAOD::TruthParticle_v1::p4
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
Definition: TruthParticle_v1.cxx:191
DataVector::erase
iterator erase(iterator position)
Remove element at a given position.
xAOD::TruthParticle_v1::pdgId
int pdgId() const
PDG ID code.
SG::DataProxy
Definition: DataProxy.h:45
xAOD::Jet_v1::pt
virtual double pt() const
The transverse momentum ( ) of the particle.
Definition: Jet_v1.cxx:44
TauAnalysisTools::BuildTruthTaus::examineTruthTauDecay
StatusCode examineTruthTauDecay(const xAOD::TruthParticle &xTruthParticle, TauTruthInfo &truthInfo) const
Definition: BuildTruthTaus.cxx:349
xAOD::TruthParticle_v1::isNeutral
bool isNeutral() const
Whether the particle is electrically neutral.
xAOD::TruthVertex_v1::outgoingParticle
const TruthParticle_v1 * outgoingParticle(size_t index) const
Get one of the outgoing particles.
Definition: TruthVertex_v1.cxx:120
DataVector::begin
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
fitman.k
k
Definition: fitman.py:528
xAOD::TruthParticle_v1::isHadron
bool isHadron() const
Whether the particle is a hadron.
ServiceHandle< IIncidentSvc >
asg::AsgMetadataTool::m_useIncidents
bool m_useIncidents
Definition: AsgMetadataTool.h:128