ATLAS Offline Software
Public Member Functions | Public Attributes | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
MuonPerformanceAlg Class Reference

#include <MuonPerformanceAlg.h>

Inheritance diagram for MuonPerformanceAlg:
Collaboration diagram for MuonPerformanceAlg:

Public Member Functions

 MuonPerformanceAlg (const std::string &name, ISvcLocator *pSvcLocator)
 
virtual ~MuonPerformanceAlg ()=default
 
virtual StatusCode initialize () override
 
virtual StatusCode execute () override
 
virtual StatusCode finalize () override
 
virtual unsigned int cardinality () const override final
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. 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 > &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
 

Public Attributes

std::vector< std::string > m_muonLocationList
 

Protected Member Functions

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

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

bool passID (const xAOD::TrackParticle *, bool debug) const
 
void print (const std::string &txt, const xAOD::TruthParticle *muon)
 
void print (const std::string &txt, const xAOD::Muon *muon)
 
bool selectPdg (int pdg) const
 
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

bool m_writeToFile
 name of external file to write statistics More...
 
std::string m_fileName
 
std::ofstream m_fileOutput
 output file More...
 
SG::ReadHandleKey< xAOD::MuonContainer > m_muonsNameKey {this, "MuonContainerName", "Muons", "muon container"}
 
SG::ReadDecorHandleKey< xAOD::MuonContainer > m_muonTruthParticleKey {this,"muonTruthParticleLinkKey", m_muonsNameKey, "truthParticleLink"}
 
unsigned int m_nevents
 
std::vector< std::string > m_hitCutString
 
std::vector< int > m_ntruth
 
std::vector< int > m_ntruth5
 
std::vector< int > m_ntruth10
 
std::vector< int > m_nfound
 
std::vector< int > m_nfound5
 
std::vector< int > m_nfound10
 
std::vector< int > m_nfoundr
 
std::vector< int > m_nfoundr5
 
std::vector< int > m_nfoundr10
 
std::vector< int > m_nreco
 
std::vector< int > m_nreco5
 
std::vector< int > m_nreco10
 
SG::ReadHandleKey< xAOD::EventInfo > m_eventInfo {this, "EventInfo", "EventInfo", "event info"}
 
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthMuons {this, "TruthMuons", "MuonTruthParticles", "truth muons"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthMuonRecoMuonLinkKey {this,"truthMuonRecoMuonLinkKey",m_truthMuons, "recoMuonLink"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthMuonTruthType {this,"truthMuonTruthType",m_truthMuons, "truthType"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthMuonTruthOrigin {this,"truthMuonTruthOrigin",m_truthMuons, "truthOrigin"}
 
int m_runNumber
 
int m_eventNumber
 
IntegerArrayProperty m_pdgsToBeConsidered
 
std::set< int > m_selectedPdgs
 
DataObjIDColl m_extendedExtraObjects
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 23 of file MuonPerformanceAlg.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ MuonPerformanceAlg()

MuonPerformanceAlg::MuonPerformanceAlg ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Definition at line 20 of file MuonPerformanceAlg.cxx.

20  :
21  AthAlgorithm(name, pSvcLocator), m_writeToFile(false), m_nevents(0), m_runNumber(0), m_eventNumber(0) {
22  declareProperty("writeToFile", m_writeToFile = false);
23  declareProperty("FileName", m_fileName = "MuonPerformanceAlg.txt");
24  declareProperty("ConsideredPDGs", m_pdgsToBeConsidered);
25 }

◆ ~MuonPerformanceAlg()

virtual MuonPerformanceAlg::~MuonPerformanceAlg ( )
virtualdefault

Member Function Documentation

◆ cardinality()

virtual unsigned int MuonPerformanceAlg::cardinality ( ) const
inlinefinaloverridevirtual

Definition at line 33 of file MuonPerformanceAlg.h.

33 { return 1;}

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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; }

◆ execute()

StatusCode MuonPerformanceAlg::execute ( )
overridevirtual

Definition at line 126 of file MuonPerformanceAlg.cxx.

126  {
127  SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInfo);
128 
129  m_runNumber = eventInfo->runNumber();
130  m_eventNumber = eventInfo->eventNumber();
131 
132  SG::ReadHandle<xAOD::TruthParticleContainer> TruthMuons(m_truthMuons);
133  if (!TruthMuons.isPresent()) {
134  ATH_MSG_DEBUG("no truth muon collection present");
135  return StatusCode::SUCCESS;
136  }
137  if (!TruthMuons.isValid()) {
138  ATH_MSG_WARNING(m_truthMuons.key() << " not valid");
139  return StatusCode::FAILURE;
140  }
141  ATH_MSG_VERBOSE("Retrieved truth muons " << TruthMuons->size());
142  typedef ElementLink<xAOD::MuonContainer> MuonLink;
143 
144  static const SG::ConstAccessor<int> truthTypeAcc("truthType");
145  static const SG::ConstAccessor<int> truthOriginAcc("truthOrigin");
146  static const SG::ConstAccessor<MuonLink> recoMuonLinkAcc("recoMuonLink");
147 
148  m_nevents += 1;
149 
150  static const SG::ConstAccessor<uint8_t> nprecLayersAcc("nprecLayers");
151 
152  for (const auto truthMu : *TruthMuons) {
153  MuonLink link;
154  const int theType = truthTypeAcc(*truthMu);
155  const int theOrigin = truthOriginAcc(*truthMu);
156  ATH_MSG_VERBOSE("Truth muon: pt " << truthMu->pt() << " eta " << truthMu->eta());
157  ATH_MSG_VERBOSE("first loop: type " << theType << " origin " << theOrigin);
158  if (truthMu->pt() < 2000. || std::abs(truthMu->eta()) > 2.8) continue;
159  if (std::abs(truthMu->eta()) > 2.5)
160  ATH_MSG_VERBOSE(" SA |eta| > 2.5 muon with truth prec layers " << (int)nprecLayersAcc(*truthMu));
161 
162  if (std::abs(truthMu->eta()) > 2.5 && (int)nprecLayersAcc(*truthMu) < 2) continue;
163  if (theType != 6 && theType != 7) continue;
164  if (theOrigin == 0 || theOrigin > 17) continue;
165  bool insideID = false;
166  if (std::abs(truthMu->eta()) < 2.0) insideID = true;
167  ATH_MSG_VERBOSE("Accepted Truth muon: pt " << truthMu->pt() << " eta " << truthMu->eta());
168 
169  if (!insideID) m_ntruth[0] += 1;
170  if (!insideID) m_ntruth[9] += 1;
171  if (!insideID) m_ntruth[10] += 1;
172  if (insideID)
173  for (int n = 1; n < 6; n++) m_ntruth[n] += 1;
174  if (insideID)
175  for (int n = 7; n < 11; n++) m_ntruth[n] += 1;
176 
177  if (truthMu->pt() > 5000.) {
178  if (!insideID) m_ntruth5[0] += 1;
179  if (!insideID) m_ntruth5[9] += 1;
180  if (!insideID) m_ntruth5[10] += 1;
181  if (insideID)
182  for (int n = 1; n < 6; n++) m_ntruth5[n] += 1;
183  if (insideID)
184  for (int n = 7; n < 11; n++) m_ntruth5[n] += 1;
185  }
186 
187  if (truthMu->pt() > 10000.) {
188  if (!insideID) m_ntruth10[0] += 1;
189  if (!insideID) m_ntruth10[9] += 1;
190  if (!insideID) m_ntruth10[10] += 1;
191  if (insideID)
192  for (int n = 1; n < 6; n++) m_ntruth10[n] += 1;
193  if (insideID)
194  for (int n = 7; n < 11; n++) m_ntruth10[n] += 1;
195  }
196 
197  if (recoMuonLinkAcc.isAvailable(*truthMu)) {
198  link = recoMuonLinkAcc(*truthMu);
199  ATH_MSG_VERBOSE(" link " << link.isValid());
200  if (link.isValid()) {
201  bool loose = false;
202  bool medium = false;
203  bool tight = false;
204  if (((*link)->quality() <= xAOD::Muon_v1::Loose)) loose = true;
205  if (((*link)->quality() <= xAOD::Muon_v1::Medium)) medium = true;
206  if (((*link)->quality() <= xAOD::Muon_v1::Tight)) tight = true;
207  if (!loose) ATH_MSG_WARNING("Muon did not pass Loose WP");
208  if (insideID) {
209  const xAOD::TrackParticle* tp = (*link)->primaryTrackParticle();
210  if (tp) {
211  } else
212  ATH_MSG_VERBOSE("No Track particle found on recoMuonLink");
213  // ID efficiency can only be measured on muons in the container
214 
215  m_ntruth[6] += 1;
216  if (truthMu->pt() > 5000.) m_ntruth5[6] += 1;
217  if (truthMu->pt() > 10000.) m_ntruth10[6] += 1;
218 
219  bool passesIDcuts = passID(tp, true);
220  ATH_MSG_VERBOSE(" all authors " << (*link)->allAuthors() << " passesIDcuts " << passesIDcuts);
221  if (passesIDcuts != (*link)->passesIDCuts()) {
222  ATH_MSG_DEBUG(" PROBLEM passedIDcuts from xAOD muon: " << (*link)->passesIDCuts()
223  << " BUT MuonPerformancAlg code gives " << passesIDcuts);
224  }
225  // Use xAOD Muon ID hit selection cuts
226  passesIDcuts = (*link)->passesIDCuts();
227  if (passesIDcuts) {
228  if (nprecLayersAcc(*truthMu) > 0) {
229  m_ntruth[11] += 1;
230  if (truthMu->pt() > 5000.) m_ntruth5[11] += 1;
231  if (truthMu->pt() > 10000.) m_ntruth10[11] += 1;
232  }
233  if (loose) {
234  m_nfound[10] += 1;
235  } else
236  print(" Muon not found by Loose ", truthMu);
237  if (medium) {
238  m_nfound[9] += 1;
239  } else if (loose)
240  print(" Muon not found by Medium ", truthMu);
241  if (tight) {
242  m_nfound[8] += 1;
243  } else if (medium)
244  print(" Muon not found by Tight ", truthMu);
245  m_nfound[6] += 1;
246  if (truthMu->pt() > 5000.) {
247  if (loose) m_nfound5[10] += 1;
248  if (medium) m_nfound5[9] += 1;
249  if (tight) m_nfound5[8] += 1;
250  m_nfound5[6] += 1;
251  }
252  if (truthMu->pt() > 10000.) {
253  if (loose) m_nfound10[10] += 1;
254  if (medium) m_nfound10[9] += 1;
255  if (tight) m_nfound10[8] += 1;
256  m_nfound10[6] += 1;
257  }
258  // CaloTag or Calolikelihood
259  if (((*link)->allAuthors() & 64 * 4) || ((*link)->allAuthors() & 64 * 8)) {
260  m_nfound[5] += 1;
261  if (truthMu->pt() > 5000.) m_nfound5[5] += 1;
262  if (truthMu->pt() > 10000.) m_nfound10[5] += 1;
263  } else
264  print(" Muon not found by CaloTag and Calolikelihood ", truthMu);
265  // MuidSA
266  if (((*link)->allAuthors() & 32)) {
267  m_nfound[7] += 1;
268  if (truthMu->pt() > 5000.) m_nfound5[7] += 1;
269  if (truthMu->pt() > 10000.) m_nfound10[7] += 1;
270  } else if (!((*link)->allAuthors() & 2))
271  print(" Muon not found by MuidSA ", truthMu);
272  // MuTag or MuTagIMO
273  if (((*link)->allAuthors() & 8) || ((*link)->allAuthors() & 16)) {
274  m_nfound[4] += 1;
275  if (truthMu->pt() > 5000.) m_nfound5[4] += 1;
276  if (truthMu->pt() > 10000.) m_nfound10[4] += 1;
277  } else
278  print(" Muon not found by MuTagIMO ", truthMu);
279  // MuGirl
280  if (((*link)->allAuthors() & 64)) {
281  m_nfound[3] += 1;
282  if (truthMu->pt() > 5000.) m_nfound5[3] += 1;
283  if (truthMu->pt() > 10000.) m_nfound10[3] += 1;
284  } else
285  print(" Muon not found by MuGirl ", truthMu);
286  // MuidCo
287  if (((*link)->allAuthors() & 2)) {
288  m_nfound[2] += 1;
289  if (truthMu->pt() > 5000.) m_nfound5[2] += 1;
290  if (truthMu->pt() > 10000.) m_nfound10[2] += 1;
291  } else
292  print(" Muon not found by MuidCo ", truthMu);
293  // MuidCo || STACO || Combined
294  if (((*link)->allAuthors() & 2) || ((*link)->allAuthors() & 4) ||
295  (*link)->muonType() == xAOD::Muon::MuonType::Combined) {
296  m_nfound[1] += 1;
297  if (truthMu->pt() > 5000.) m_nfound5[1] += 1;
298  if (truthMu->pt() > 10000.) m_nfound10[1] += 1;
299  if (nprecLayersAcc(*truthMu) > 0) {
300  m_nfound[11] += 1;
301  if (truthMu->pt() > 5000.) m_nfound5[11] += 1;
302  if (truthMu->pt() > 10000.) m_nfound10[11] += 1;
303  }
304  } else
305  print(" Combined Muon not found ", truthMu);
306  } else {
307  // SA or CB not passing ID cuts
308  if (((*link)->allAuthors() & 32) || (*link)->muonType() == xAOD::Muon::MuonType::Combined) {
309  m_nfound[7] += 1;
310  if (truthMu->pt() > 5000.) m_nfound5[7] += 1;
311  if (truthMu->pt() > 10000.) m_nfound10[7] += 1;
312  } else
313  print(" SA (CB) Muon not found ", truthMu);
314  }
315  }
316  // MuidSA
317  else if (!insideID) {
318  if ((((*link)->allAuthors() & 32)) || ((*link)->allAuthors() & 2) || ((*link)->allAuthors() & 4)) {
319  m_nfound[0] += 1;
320  if (truthMu->pt() > 5000.) m_nfound5[0] += 1;
321  if (truthMu->pt() > 10000.) m_nfound10[0] += 1;
322  } else
323  print(" Muon not found by MuidSA endcap ", truthMu);
324  if (loose) {
325  m_nfound[10] += 1;
326  if (truthMu->pt() > 5000.) m_nfound5[10] += 1;
327  if (truthMu->pt() > 10000.) m_nfound10[10] += 1;
328  } else
329  print(" Muon not found by Loose endcap ", truthMu);
330  if (medium) {
331  m_nfound[9] += 1;
332  if (truthMu->pt() > 5000.) m_nfound5[9] += 1;
333  if (truthMu->pt() > 10000.) m_nfound10[9] += 1;
334  } else if (loose)
335  print(" Muon not found by Medium endcap ", truthMu);
336  }
337  } else {
338  print(" No link Muon not found by CaloTag and Calolikelihood ", truthMu);
339  print(" No link Muon not found by MuidSA ", truthMu);
340  print(" No link Muon not found by MuTagIMO ", truthMu);
341  print(" No link Muon not found by MuGirl ", truthMu);
342  print(" No link Muon not found by MuidCo ", truthMu);
343  print(" No link Combined Muon not found ", truthMu);
344  }
345  } // end valid link
346  }
347 
348  for (const auto truthMu : *TruthMuons) {
349  MuonLink link;
350  if (truthMu->pt() < 2000. || std::abs(truthMu->eta()) > 3.) continue;
351  const int theType = truthTypeAcc(*truthMu);
352  const int theOrigin = truthOriginAcc(*truthMu);
353  if (theType != 6 && theType != 7) continue;
354  if (theOrigin == 0 || theOrigin > 17) continue;
355  if (std::abs(truthMu->eta()) > 2.5 && (int)nprecLayersAcc(*truthMu) < 2) continue;
356  bool insideID = false;
357  if (recoMuonLinkAcc.isAvailable(*truthMu)) {
358  link = recoMuonLinkAcc(*truthMu);
359  if (link.isValid()) {
360  const xAOD::TrackParticle* tp = (*link)->primaryTrackParticle();
361  if (tp) {
362  if ((*link)->pt() < 2000. || std::abs((*link)->eta()) > 2.8) continue;
363  if (std::abs((*link)->eta()) < 2.0) insideID = true;
364  bool loose = false;
365  bool medium = false;
366  bool tight = false;
367  if (((*link)->quality() == xAOD::Muon_v1::Loose)) loose = true;
368  if (((*link)->quality() == xAOD::Muon_v1::Medium)) loose = true;
369  if (((*link)->quality() == xAOD::Muon_v1::Tight)) loose = true;
370  if (((*link)->quality() == xAOD::Muon_v1::Medium)) medium = true;
371  if (((*link)->quality() == xAOD::Muon_v1::Tight)) medium = true;
372  if (((*link)->quality() == xAOD::Muon_v1::Tight)) tight = true;
373  if (insideID) {
374  bool passesIDcuts = passID(tp, false);
375  // Use xAOD Muon ID hit selection cuts
376  passesIDcuts = (*link)->passesIDCuts();
377  if (passesIDcuts) {
378  if (loose) m_nfoundr[10] += 1;
379  if (medium) m_nfoundr[9] += 1;
380  if (tight) m_nfoundr[8] += 1;
381  m_nfoundr[6] += 1;
382  if (((*link)->allAuthors() & 64 * 4) || ((*link)->allAuthors() & 64 * 8)) m_nfoundr[5] += 1;
383  if (((*link)->allAuthors() & 32)) m_nfoundr[7] += 1;
384  if (((*link)->allAuthors() & 8) || ((*link)->allAuthors() & 16)) m_nfoundr[4] += 1;
385  if (((*link)->allAuthors() & 64)) m_nfoundr[3] += 1;
386  if (((*link)->allAuthors() & 2)) m_nfoundr[2] += 1;
387  if (((*link)->allAuthors() & 2) || ((*link)->allAuthors() & 4) ||
388  (*link)->muonType() == xAOD::Muon::MuonType::Combined)
389  m_nfoundr[1] += 1;
390  if ((*link)->pt() > 5000.) {
391  if (loose) m_nfoundr5[10] += 1;
392  if (medium) m_nfoundr5[9] += 1;
393  if (tight) m_nfoundr5[8] += 1;
394  m_nfoundr5[6] += 1;
395  if (((*link)->allAuthors() & 64 * 4) || ((*link)->allAuthors() & 64 * 8)) m_nfoundr5[5] += 1;
396  if (((*link)->allAuthors() & 32)) m_nfoundr5[7] += 1;
397  if (((*link)->allAuthors() & 8) || ((*link)->allAuthors() & 16)) m_nfoundr5[4] += 1;
398  if (((*link)->allAuthors() & 64)) m_nfoundr5[3] += 1;
399  if (((*link)->allAuthors() & 2)) m_nfoundr5[2] += 1;
400  if (((*link)->allAuthors() & 2) || ((*link)->allAuthors() & 4) ||
401  (*link)->muonType() == xAOD::Muon::MuonType::Combined)
402  m_nfoundr5[1] += 1;
403  }
404  if ((*link)->pt() > 10000.) {
405  if (loose) m_nfoundr10[10] += 1;
406  if (medium) m_nfoundr10[9] += 1;
407  if (tight) m_nfoundr10[8] += 1;
408  m_nfoundr10[6] += 1;
409  if (((*link)->allAuthors() & 64 * 4) || ((*link)->allAuthors() & 64 * 8)) m_nfoundr10[5] += 1;
410  if (((*link)->allAuthors() & 32)) m_nfoundr10[7] += 1;
411  if (((*link)->allAuthors() & 8) || ((*link)->allAuthors() & 16)) m_nfoundr10[4] += 1;
412  if (((*link)->allAuthors() & 64)) m_nfoundr10[3] += 1;
413  if (((*link)->allAuthors() & 2)) m_nfoundr10[2] += 1;
414  if (((*link)->allAuthors() & 2) || ((*link)->allAuthors() & 4) ||
415  (*link)->muonType() == xAOD::Muon::MuonType::Combined)
416  m_nfoundr10[1] += 1;
417  }
418  } else {
419  if (((*link)->allAuthors() & 32) || (*link)->muonType() == xAOD::Muon::MuonType::Combined) {
420  m_nfoundr[7] += 1;
421  if ((*link)->pt() > 5000.) m_nfoundr5[7] += 1;
422  if ((*link)->pt() > 10000.) m_nfoundr10[7] += 1;
423  }
424  }
425  } else if (!insideID) {
426  if ((((*link)->allAuthors() & 32)) || ((*link)->allAuthors() & 2) || ((*link)->allAuthors() & 4)) {
427  m_nfoundr[0] += 1;
428  if ((*link)->pt() > 5000.) m_nfoundr5[0] += 1;
429  if ((*link)->pt() > 10000.) m_nfoundr10[0] += 1;
430  }
431  if (loose) m_nfoundr[10] += 1;
432  if (medium) m_nfoundr[9] += 1;
433  if ((*link)->pt() > 5000.) {
434  if (loose) m_nfoundr5[10] += 1;
435  if (medium) m_nfoundr5[9] += 1;
436  }
437  if ((*link)->pt() > 10000.) {
438  if (loose) m_nfoundr10[10] += 1;
439  if (medium) m_nfoundr10[9] += 1;
440  }
441  }
442  } else
443  ATH_MSG_VERBOSE("No Track particle found on recoMuonLink");
444  } // not valid muonLink
445  } // no muonLink
446  }
447 
448  SG::ReadHandle<xAOD::MuonContainer> Muons(m_muonsNameKey);
449 
450  if (!Muons.isValid()) {
451  ATH_MSG_ERROR("Couldn't retrieve Muons container with key: " << m_muonsNameKey.key());
452  return StatusCode::FAILURE;
453  }
454  ATH_MSG_VERBOSE("Retrieved muons " << Muons->size());
455 
456  for (const auto mu : *Muons) {
457  if (mu->pt() < 2000. || std::abs(mu->eta()) > 2.8) continue;
458  const xAOD::TrackParticle* tp = mu->primaryTrackParticle();
459 
460  if (tp) {
461  ElementLink<xAOD::TruthParticleContainer> truthLink;
462  bool passesIDcuts = passID(tp, false);
463  // Use xAOD Muon ID hit selection cuts
464  passesIDcuts = mu->passesIDCuts();
465 
466  bool insideID = false;
467  if (std::abs(mu->eta()) < 2.0) insideID = true;
468  static const SG::ConstAccessor<ElementLink<xAOD::TruthParticleContainer> >
469  truthParticleLinkAcc ("truthParticleLink");
470  if (truthParticleLinkAcc.isAvailable(*tp))
471  truthLink = truthParticleLinkAcc(*tp);
472  bool fake = true;
473  if (truthLink.isValid()) {
474  if (selectPdg((*truthLink)->pdgId())) fake = false;
475  }
476  bool loose = false;
477  bool medium = false;
478  bool tight = false;
479  if (mu->quality() == xAOD::Muon_v1::Loose) loose = true;
480  if (mu->quality() == xAOD::Muon_v1::Medium) loose = true;
481  if (mu->quality() == xAOD::Muon_v1::Tight) loose = true;
482  if (mu->quality() == xAOD::Muon_v1::Medium) medium = true;
483  if (mu->quality() == xAOD::Muon_v1::Tight) medium = true;
484  if (mu->quality() == xAOD::Muon_v1::Tight) tight = true;
485  if (insideID) {
486  if (passesIDcuts) {
487  if (loose) m_nreco[10] += 1;
488  if (medium) m_nreco[9] += 1;
489  if (tight) m_nreco[8] += 1;
490  m_nreco[6] += 1;
491  if ((mu->allAuthors() & 64 * 4) || (mu->allAuthors() & 64 * 8)) m_nreco[5] += 1;
492  if ((mu->allAuthors() & 32)) m_nreco[7] += 1;
493  if ((mu->allAuthors() & 8) || (mu->allAuthors() & 16)) m_nreco[4] += 1;
494  if ((mu->allAuthors() & 64)) m_nreco[3] += 1;
495  if ((mu->allAuthors() & 2)) m_nreco[2] += 1;
496  if ((mu->allAuthors() & 2) || (mu->allAuthors() & 4) || mu->muonType() == xAOD::Muon::MuonType::Combined)
497  m_nreco[1] += 1;
498  if (mu->pt() > 5000.) {
499  if (loose) m_nreco5[10] += 1;
500  if (medium) m_nreco5[9] += 1;
501  if (tight) m_nreco5[8] += 1;
502  m_nreco5[6] += 1;
503  if ((mu->allAuthors() & 64 * 4) || (mu->allAuthors() & 64 * 8)) m_nreco5[5] += 1;
504  if ((mu->allAuthors() & 32)) m_nreco5[7] += 1;
505  if ((mu->allAuthors() & 8) || (mu->allAuthors() & 16)) m_nreco5[4] += 1;
506  if ((mu->allAuthors() & 64)) m_nreco5[3] += 1;
507  if ((mu->allAuthors() & 2)) m_nreco5[2] += 1;
508  if ((mu->allAuthors() & 2) || (mu->allAuthors() & 4) || mu->muonType() == xAOD::Muon::MuonType::Combined)
509  m_nreco5[1] += 1;
510  }
511  if (mu->pt() > 10000.) {
512  if (loose) m_nreco10[10] += 1;
513  if (medium) m_nreco10[9] += 1;
514  if (tight) m_nreco10[8] += 1;
515  m_nreco10[6] += 1;
516  if ((mu->allAuthors() & 64 * 4) || (mu->allAuthors() & 64 * 8)) m_nreco10[5] += 1;
517  if ((mu->allAuthors() & 32)) m_nreco10[7] += 1;
518  if ((mu->allAuthors() & 8) || (mu->allAuthors() & 16)) m_nreco10[4] += 1;
519  if ((mu->allAuthors() & 64)) m_nreco10[3] += 1;
520  if ((mu->allAuthors() & 2)) m_nreco10[2] += 1;
521  if ((mu->allAuthors() & 2) || (mu->allAuthors() & 4) || mu->muonType() == xAOD::Muon::MuonType::Combined)
522  m_nreco10[1] += 1;
523  }
524  if (fake) {
525  if ((mu->quality() == xAOD::Muon_v1::Loose)) print(" Fake muon found by Loose ", mu);
526  if ((mu->quality() == xAOD::Muon_v1::Medium)) print(" Fake muon found by Medium ", mu);
527  if ((mu->quality() == xAOD::Muon_v1::Tight)) print(" Fake muon found by Tight ", mu);
528  if ((mu->allAuthors() & 64 * 4) || (mu->allAuthors() & 64 * 8))
529  print(" Fake muon found by CaloTag and Calolikelihood ", mu);
530  if ((mu->allAuthors() & 32)) print(" Fake muon found by MuidSA ", mu);
531  if ((mu->allAuthors() & 8) || (mu->allAuthors() & 16)) print(" Fake muon found by MuTagIMO ", mu);
532  if ((mu->allAuthors() & 64)) print(" Fake muon found by MuGirl ", mu);
533  if ((mu->allAuthors() & 2)) print(" Fake muon found by MuidCo ", mu);
534  if ((mu->allAuthors() & 2) || (mu->allAuthors() & 4) || mu->muonType() == xAOD::Muon::MuonType::Combined)
535  print(" Fake Combined muon ", mu);
536  }
537  } else {
538  if ((mu->allAuthors() & 32) || mu->muonType() == xAOD::Muon::MuonType::Combined) {
539  m_nreco[7] += 1;
540  if (fake) print(" Fake muon found by MuidSA with no ID", mu);
541  if (mu->pt() > 5000.) m_nreco5[7] += 1;
542  if (mu->pt() > 10000.) m_nreco10[7] += 1;
543  }
544  }
545  } else if (!insideID) {
546  if (loose) m_nreco[10] += 1;
547  if (medium) m_nreco[9] += 1;
548  if (loose && fake) print(" Fake muon found by Loose Endcap ", mu);
549  if (medium && fake && !loose) print(" Fake muon found by Medium Endcap ", mu);
550  if (tight && !medium && fake) print(" Fake muon found by Tight Endcap ", mu);
551  if (mu->pt() > 5000.) {
552  if (loose) m_nreco5[10] += 1;
553  if (medium) m_nreco5[9] += 1;
554  }
555  if (mu->pt() > 10000.) {
556  if (loose) m_nreco10[10] += 1;
557  if (medium) m_nreco10[9] += 1;
558  }
559  if (((mu->allAuthors() & 32)) || (mu->allAuthors() & 2) || (mu->allAuthors() & 4)) {
560  m_nreco[0] += 1;
561  if (fake) print(" Fake muon found by MuidSA Endcap ", mu);
562  if (mu->pt() > 5000.) m_nreco5[0] += 1;
563  if (mu->pt() > 10000.) m_nreco10[0] += 1;
564  }
565  }
566  }
567  }
568 
569  return StatusCode::SUCCESS;
570 }

◆ extraDeps_update_handler()

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

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

51 {
52  // If we didn't find any symlinks to add, just return the collection
53  // from the base class. Otherwise, return the extended collection.
54  if (!m_extendedExtraObjects.empty()) {
55  return m_extendedExtraObjects;
56  }
57  return Algorithm::extraOutputDeps();
58 }

◆ finalize()

StatusCode MuonPerformanceAlg::finalize ( )
overridevirtual

Definition at line 596 of file MuonPerformanceAlg.cxx.

596  {
597  std::ofstream fileOutput;
598  std::string outfile = "muonPerformance_xAOD.txt";
599  fileOutput.open(outfile.c_str(), std::ios::trunc);
600  std::ostringstream sout;
601  sout.precision(4);
602 
603  unsigned int width = 9;
604  unsigned int precision = 3;
605 
606  sout << std::endl;
607  sout << " Summary of the xAOD Muon performance " << std::endl;
608  sout << " Muon type #Truth muons Efficiency #Truth muons Efficiency #Truth muons Efficiency" << std::endl;
609  sout << " (pt>2) (pt>5) (pt>10 GeV/c) " << std::endl;
610  for (unsigned int i = 0; i < 12; ++i) {
611  if (m_ntruth[i] == 0) {
612  sout << " " << std::endl;
613  } else {
614  sout << " " << m_hitCutString[i];
615  sout << std::setw(width) << std::setprecision(precision);
616  sout << static_cast<double>(m_ntruth[i]);
617  sout << " ";
618  sout << std::setw(width) << std::setprecision(precision);
619  sout << static_cast<double>(m_nfound[i]) / static_cast<double>(m_ntruth[i]);
620  sout << std::setw(width) << std::setprecision(precision);
621  sout << static_cast<double>(m_ntruth5[i]);
622  sout << " ";
623  sout << std::setw(width) << std::setprecision(precision);
624  if (m_ntruth5[i] != 0) {
625  sout << static_cast<double>(m_nfound5[i]) / static_cast<double>(m_ntruth5[i]);
626  } else
627  sout << 0.;
628 
629  sout << std::setw(width) << std::setprecision(precision);
630  sout << static_cast<double>(m_ntruth10[i]);
631  sout << " ";
632  sout << std::setw(width) << std::setprecision(precision);
633  if (m_ntruth10[i] != 0) {
634  sout << static_cast<double>(m_nfound10[i]) / static_cast<double>(m_ntruth10[i]) << std::endl;
635  } else
636  sout << 0. << std::endl;
637  }
638  }
639  sout << " The efficiency of the ID is calculated inside |eta| < 2 for the MCP ID hit selection cuts (it uses identified - Combined, "
640  "Tagged or CaloTagged - muons with an ID track)"
641  << std::endl;
642  sout << " The efficiencies for CB all, MuidCB, MuGirl, Tag and Calo include the MCP ID cuts" << std::endl;
643  sout << " The Tight, Medium and Loose efficiencies include MCP ID cuts for muons |eta| < 2" << std::endl;
644  sout << " The SA 2.0 for |eta| >2 and SA (no ID) for |eta| < 2 doesnot include MCP ID cuts" << std::endl;
645  sout << " The Combined efficiency is defined ID tracks after MCP cuts and requiring at least 1 muon station at truth level"
646  << std::endl;
647  sout << " Fakes are calculated with the selections listed above" << std::endl;
648  sout << std::endl;
649 
650  sout << " Muon type #Fake muons rate #Fake muons rate #Fake muons rate " << std::endl;
651  sout << " (all pt>0,2) (pt>5) (pt>10 GeV/c) " << std::endl;
652  for (unsigned int i = 0; i < 11; ++i) {
653  if (i == 6) continue;
654  if (m_nevents == 0) {
655  sout << " " << std::endl;
656  } else {
657  sout << " " << m_hitCutString[i];
658  sout << std::setw(width) << std::setprecision(precision);
659  sout << static_cast<double>(m_nreco[i] - m_nfoundr[i]);
660  sout << std::setw(width) << std::setprecision(precision);
661  sout << static_cast<double>(m_nreco[i] - m_nfoundr[i]) / static_cast<double>(m_nevents);
662  sout << std::setw(width) << std::setprecision(precision);
663  sout << static_cast<double>(m_nreco5[i] - m_nfoundr5[i]);
664  sout << std::setw(width) << std::setprecision(precision);
665  sout << static_cast<double>(m_nreco5[i] - m_nfoundr5[i]) / static_cast<double>(m_nevents);
666  sout << std::setw(width) << std::setprecision(precision);
667  sout << static_cast<double>(m_nreco10[i] - m_nfoundr10[i]);
668  sout << std::setw(width) << std::setprecision(precision);
669  sout << static_cast<double>(m_nreco10[i] - m_nfoundr10[i]) / static_cast<double>(m_nevents) << std::endl;
670  }
671  }
672  sout << std::endl;
673  fileOutput << sout.str() << std::endl;
674  fileOutput.close();
675  return StatusCode::SUCCESS;
676 }

◆ initialize()

StatusCode MuonPerformanceAlg::initialize ( )
overridevirtual

Definition at line 27 of file MuonPerformanceAlg.cxx.

27  {
28  // initialize cuts, please make sure the number of bins and the sizes of the cuts + string are always the same
29  unsigned int nbins = 12;
30  m_nevents = 0;
31  m_ntruth.resize(nbins);
32  m_ntruth5.resize(nbins);
33  m_ntruth10.resize(nbins);
34  m_nfound.resize(nbins);
35  m_nfound5.resize(nbins);
36  m_nfound10.resize(nbins);
37  m_nfoundr.resize(nbins);
38  m_nfoundr5.resize(nbins);
39  m_nfoundr10.resize(nbins);
40  m_nreco.resize(nbins);
41  m_nreco5.resize(nbins);
42  m_nreco10.resize(nbins);
43  m_hitCutString = {"SA 2.0 ", "CB all ", "MuidCB ", "MuGirl ", "Tag ", "Calo ",
44  "ID ", "SA(no ID)", "Tight ", "Medium ", "Loose ", "Combined "};
45 
46  // add muons
47  if (m_pdgsToBeConsidered.value().empty()) {
48  m_selectedPdgs.insert(13);
49  m_selectedPdgs.insert(-13);
50  } else {
51  // add pdgs
52  for (auto pdg : m_pdgsToBeConsidered.value()) m_selectedPdgs.insert(pdg);
53  }
54  if (!m_muonsNameKey.key().empty()) ATH_CHECK(m_muonsNameKey.initialize());
55 
56  ATH_CHECK(m_eventInfo.initialize());
57  ATH_CHECK(m_truthMuons.initialize());
58  ATH_CHECK(m_muonTruthParticleKey.initialize());
59  ATH_CHECK(m_truthMuonRecoMuonLinkKey.initialize());
60  ATH_CHECK(m_truthMuonTruthType.initialize());
61  ATH_CHECK(m_truthMuonTruthOrigin.initialize());
62  return StatusCode::SUCCESS;
63 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

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

◆ msgLvl()

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

◆ passID()

bool MuonPerformanceAlg::passID ( const xAOD::TrackParticle *  tp,
bool  debug 
) const
private

Definition at line 67 of file MuonPerformanceAlg.cxx.

67  {
68  uint8_t nblh = 0x0;
69  uint8_t eblh = 0x0;
70  uint8_t nblo = 0x0;
71  uint8_t nphi = 0x0;
72  uint8_t npds = 0x0;
73  uint8_t nsctds = 0x0;
74  uint8_t nscthi = 0x0;
75  uint8_t nsctho = 0x0;
76  uint8_t npho = 0x0;
77  uint8_t ntrthi = 0x0;
78  uint8_t ntrtol = 0x0;
79 
80  int tightness = 0;
81 
82  bool passesIDcuts = false;
83 
84  if (tp) {
85  // check blayer
86  tp->summaryValue(eblh, xAOD::expectInnermostPixelLayerHit);
87  tp->summaryValue(nblh, xAOD::numberOfInnermostPixelLayerHits);
88  tp->summaryValue(nblo, xAOD::numberOfInnermostPixelLayerOutliers);
89  if (static_cast<int>(eblh) == 0 || static_cast<int>(nblh) + static_cast<int>(nblo) > 0) tightness += 1;
90 
91  // pixel hit counts
92  tp->summaryValue(nphi, xAOD::numberOfPixelHits);
93  tp->summaryValue(npds, xAOD::numberOfPixelDeadSensors);
94  if (static_cast<int>(nphi) + static_cast<int>(npds) >= 2) tightness += 2;
95 
96  // sct hit counts
97  tp->summaryValue(nscthi, xAOD::numberOfSCTHits);
98  tp->summaryValue(nsctds, xAOD::numberOfSCTDeadSensors);
99  if (static_cast<int>(nscthi) + static_cast<int>(nsctds) > 5) tightness += 4;
100 
101  // hole cuts
102  tp->summaryValue(npho, xAOD::numberOfPixelHoles);
103  tp->summaryValue(nsctho, xAOD::numberOfSCTHoles);
104  if (static_cast<int>(npho) + static_cast<int>(nsctho) < 2) tightness += 8;
105 
106  // trt cuts
107  tp->summaryValue(ntrthi, xAOD::numberOfTRTHits);
108  tp->summaryValue(ntrtol, xAOD::numberOfTRTOutliers);
109  int ntrt = static_cast<int>(ntrthi) + static_cast<int>(ntrtol);
110  bool pass = true;
111  if (ntrt > 5 && static_cast<int>(ntrtol) / (double)ntrt > 0.9) pass = false;
112  if ((std::abs(tp->eta()) > 0.1 && std::abs(tp->eta()) < 1.9) && ntrt <= 5) pass = false;
113  if (pass) tightness += 16;
114  //
115  // B layer is not a requirement anymore
116  //
117  if (pass && tightness > 29) passesIDcuts = true;
118 
119  if (passesIDcuts && debug) ATH_MSG_VERBOSE("ID trackparticle PASSED IDcuts with tightness " << tightness);
120  }
121  if (!passesIDcuts && debug) ATH_MSG_DEBUG("ID trackparticle FAILED IDcuts with tightness " << tightness);
122 
123  return passesIDcuts;
124 }

◆ print() [1/2]

void MuonPerformanceAlg::print ( const std::string &  txt,
const xAOD::Muon *  muon 
)
private

Definition at line 583 of file MuonPerformanceAlg.cxx.

583  {
584  int nprec = 0;
585  uint8_t nPrecision = 0;
586  if (muon->primaryTrackParticleLink().isValid()) {
587  muon->primaryTrackParticle()->summaryValue(nPrecision, xAOD::numberOfPrecisionLayers);
588  nprec = static_cast<int>(nPrecision);
589  }
590  double p = sqrt(muon->e() * muon->e() - muon->m() * muon->m());
591  ATH_MSG_DEBUG(txt << " run " << m_runNumber << " event " << m_eventNumber << std::endl
592  << " eta " << muon->eta() << " phi " << muon->phi() << " q*p (GeV) " << muon->charge() * p / 1000. << " pt (GeV) "
593  << muon->pt() / 1000. << " precisionLayers " << nprec << " nr segments " << muon->nMuonSegments());
594 }

◆ print() [2/2]

void MuonPerformanceAlg::print ( const std::string &  txt,
const xAOD::TruthParticle *  muon 
)
private

Definition at line 571 of file MuonPerformanceAlg.cxx.

571  {
572  // Truth: r 4251 z 3147 theta 0.91843 phi -0.781 q*p(GeV) 1.081e+02 pt(Gev) 8.589e+01
573  // uniqueID 1489 mother 23 production vertex: r 1.06363 z -102.414
574  int q = 1;
575  if (muon->pdgId() > 0) q = -1;
576  double p = sqrt(muon->e() * muon->e() - muon->m() * muon->m());
577  static const SG::ConstAccessor<uint8_t> nprecLayersAcc("nprecLayers");
578  ATH_MSG_DEBUG(txt << " run " << m_runNumber << " event " << m_eventNumber << std::endl
579  << " Truth: pdgId " << muon->pdgId() << " uniqueID " << HepMC::uniqueID(muon) << " eta " << muon->eta() << " phi "
580  << muon->phi() << " q*p (GeV) " << q * p / 1000. << " pt (GeV) " << muon->pt() / 1000. << " precisionLayers "
581  << static_cast<int>(nprecLayersAcc(*muon)));
582 }

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

◆ selectPdg()

bool MuonPerformanceAlg::selectPdg ( int  pdg) const
inlineprivate

Definition at line 80 of file MuonPerformanceAlg.h.

80 { return m_selectedPdgs.count(pdg); }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

66  {
67  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
68 
69  if (sc.isFailure()) {
70  return sc;
71  }
72  ServiceHandle<ICondSvc> cs("CondSvc",name());
73  for (auto h : outputHandles()) {
74  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75  // do this inside the loop so we don't create the CondSvc until needed
76  if ( cs.retrieve().isFailure() ) {
77  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78  return StatusCode::SUCCESS;
79  }
80  if (cs->regHandle(this,*h).isFailure()) {
81  sc = StatusCode::FAILURE;
82  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83  << " with CondSvc");
84  }
85  }
86  }
87  return sc;
88 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_eventInfo

SG::ReadHandleKey<xAOD::EventInfo> MuonPerformanceAlg::m_eventInfo {this, "EventInfo", "EventInfo", "event info"}
private

Definition at line 66 of file MuonPerformanceAlg.h.

◆ m_eventNumber

int MuonPerformanceAlg::m_eventNumber
private

Definition at line 73 of file MuonPerformanceAlg.h.

◆ m_evtStore

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

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_fileName

std::string MuonPerformanceAlg::m_fileName
private

Definition at line 43 of file MuonPerformanceAlg.h.

◆ m_fileOutput

std::ofstream MuonPerformanceAlg::m_fileOutput
private

output file

Definition at line 46 of file MuonPerformanceAlg.h.

◆ m_hitCutString

std::vector<std::string> MuonPerformanceAlg::m_hitCutString
private

Definition at line 52 of file MuonPerformanceAlg.h.

◆ m_muonLocationList

std::vector<std::string> MuonPerformanceAlg::m_muonLocationList

Definition at line 36 of file MuonPerformanceAlg.h.

◆ m_muonsNameKey

SG::ReadHandleKey<xAOD::MuonContainer> MuonPerformanceAlg::m_muonsNameKey {this, "MuonContainerName", "Muons", "muon container"}
private

Definition at line 49 of file MuonPerformanceAlg.h.

◆ m_muonTruthParticleKey

SG::ReadDecorHandleKey<xAOD::MuonContainer> MuonPerformanceAlg::m_muonTruthParticleKey {this,"muonTruthParticleLinkKey", m_muonsNameKey, "truthParticleLink"}
private

Definition at line 50 of file MuonPerformanceAlg.h.

◆ m_nevents

unsigned int MuonPerformanceAlg::m_nevents
private

Definition at line 51 of file MuonPerformanceAlg.h.

◆ m_nfound

std::vector<int> MuonPerformanceAlg::m_nfound
private

Definition at line 56 of file MuonPerformanceAlg.h.

◆ m_nfound10

std::vector<int> MuonPerformanceAlg::m_nfound10
private

Definition at line 58 of file MuonPerformanceAlg.h.

◆ m_nfound5

std::vector<int> MuonPerformanceAlg::m_nfound5
private

Definition at line 57 of file MuonPerformanceAlg.h.

◆ m_nfoundr

std::vector<int> MuonPerformanceAlg::m_nfoundr
private

Definition at line 59 of file MuonPerformanceAlg.h.

◆ m_nfoundr10

std::vector<int> MuonPerformanceAlg::m_nfoundr10
private

Definition at line 61 of file MuonPerformanceAlg.h.

◆ m_nfoundr5

std::vector<int> MuonPerformanceAlg::m_nfoundr5
private

Definition at line 60 of file MuonPerformanceAlg.h.

◆ m_nreco

std::vector<int> MuonPerformanceAlg::m_nreco
private

Definition at line 62 of file MuonPerformanceAlg.h.

◆ m_nreco10

std::vector<int> MuonPerformanceAlg::m_nreco10
private

Definition at line 64 of file MuonPerformanceAlg.h.

◆ m_nreco5

std::vector<int> MuonPerformanceAlg::m_nreco5
private

Definition at line 63 of file MuonPerformanceAlg.h.

◆ m_ntruth

std::vector<int> MuonPerformanceAlg::m_ntruth
private

Definition at line 53 of file MuonPerformanceAlg.h.

◆ m_ntruth10

std::vector<int> MuonPerformanceAlg::m_ntruth10
private

Definition at line 55 of file MuonPerformanceAlg.h.

◆ m_ntruth5

std::vector<int> MuonPerformanceAlg::m_ntruth5
private

Definition at line 54 of file MuonPerformanceAlg.h.

◆ m_pdgsToBeConsidered

IntegerArrayProperty MuonPerformanceAlg::m_pdgsToBeConsidered
private

Definition at line 78 of file MuonPerformanceAlg.h.

◆ m_runNumber

int MuonPerformanceAlg::m_runNumber
private

Definition at line 72 of file MuonPerformanceAlg.h.

◆ m_selectedPdgs

std::set<int> MuonPerformanceAlg::m_selectedPdgs
private

Definition at line 79 of file MuonPerformanceAlg.h.

◆ m_truthMuonRecoMuonLinkKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> MuonPerformanceAlg::m_truthMuonRecoMuonLinkKey {this,"truthMuonRecoMuonLinkKey",m_truthMuons, "recoMuonLink"}
private

Definition at line 68 of file MuonPerformanceAlg.h.

◆ m_truthMuons

SG::ReadHandleKey<xAOD::TruthParticleContainer> MuonPerformanceAlg::m_truthMuons {this, "TruthMuons", "MuonTruthParticles", "truth muons"}
private

Definition at line 67 of file MuonPerformanceAlg.h.

◆ m_truthMuonTruthOrigin

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> MuonPerformanceAlg::m_truthMuonTruthOrigin {this,"truthMuonTruthOrigin",m_truthMuons, "truthOrigin"}
private

Definition at line 70 of file MuonPerformanceAlg.h.

◆ m_truthMuonTruthType

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> MuonPerformanceAlg::m_truthMuonTruthType {this,"truthMuonTruthType",m_truthMuons, "truthType"}
private

Definition at line 69 of file MuonPerformanceAlg.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_writeToFile

bool MuonPerformanceAlg::m_writeToFile
private

name of external file to write statistics

Definition at line 42 of file MuonPerformanceAlg.h.

The documentation for this class was generated from the following files:
MuonPerformanceAlg::m_nevents
unsigned int m_nevents
Definition: MuonPerformanceAlg.h:51
MuonPerformanceAlg::print
void print(const std::string &txt, const xAOD::TruthParticle *muon)
Definition: MuonPerformanceAlg.cxx:571
MuonPerformanceAlg::m_fileName
std::string m_fileName
Definition: MuonPerformanceAlg.h:43
xAOD::numberOfPixelHoles
@ numberOfPixelHoles
number of pixel layers on track with absence of hits [unit8_t].
Definition: TrackingPrimitives.h:261
LikeEnum::Loose
@ Loose
Definition: LikelihoodEnums.h:12
xAOD::muon
@ muon
Definition: TrackingPrimitives.h:195
MuonLink
ElementLink< xAOD::MuonContainer > MuonLink
Definition: BPhysHelper.cxx:21
MuonPerformanceAlg::m_muonsNameKey
SG::ReadHandleKey< xAOD::MuonContainer > m_muonsNameKey
Definition: MuonPerformanceAlg.h:49
xAOD::uint8_t
uint8_t
Definition: Muon_v1.cxx:557
MuonPerformanceAlg::selectPdg
bool selectPdg(int pdg) const
Definition: MuonPerformanceAlg.h:80
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
MuonPerformanceAlg::m_nfoundr10
std::vector< int > m_nfoundr10
Definition: MuonPerformanceAlg.h:61
MuonPerformanceAlg::m_nreco
std::vector< int > m_nreco
Definition: MuonPerformanceAlg.h:62
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
ParticleTest.tp
tp
Definition: ParticleTest.py:25
xAOD::numberOfPixelHits
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
Definition: TrackingPrimitives.h:259
xAOD::expectInnermostPixelLayerHit
@ expectInnermostPixelLayerHit
Do we expect a 0th-layer barrel hit for this track?
Definition: TrackingPrimitives.h:236
xAOD::numberOfTRTHits
@ numberOfTRTHits
number of TRT hits [unit8_t].
Definition: TrackingPrimitives.h:275
SG::ConstAccessor< int >
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
MuonPerformanceAlg::m_writeToFile
bool m_writeToFile
name of external file to write statistics
Definition: MuonPerformanceAlg.h:42
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
MuonPerformanceAlg::m_nfound
std::vector< int > m_nfound
Definition: MuonPerformanceAlg.h:56
SCT_CalibAlgs::nbins
@ nbins
Definition: SCT_CalibNumbers.h:10
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
MuonPerformanceAlg::m_selectedPdgs
std::set< int > m_selectedPdgs
Definition: MuonPerformanceAlg.h:79
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
MuonPerformanceAlg::m_runNumber
int m_runNumber
Definition: MuonPerformanceAlg.h:72
MuonPerformanceAlg::m_muonTruthParticleKey
SG::ReadDecorHandleKey< xAOD::MuonContainer > m_muonTruthParticleKey
Definition: MuonPerformanceAlg.h:50
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:27
MuonPerformanceAlg::m_nfoundr5
std::vector< int > m_nfoundr5
Definition: MuonPerformanceAlg.h:60
xAOD::numberOfInnermostPixelLayerOutliers
@ numberOfInnermostPixelLayerOutliers
number of 0th layer barrel outliers
Definition: TrackingPrimitives.h:238
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
lumiFormat.i
int i
Definition: lumiFormat.py:85
SG::ReadDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
beamspotman.n
n
Definition: beamspotman.py:731
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
HepMC::uniqueID
int uniqueID(const T &p)
Definition: MagicNumbers.h:109
jobOptions.fileOutput
fileOutput
Definition: jobOptions.py:39
MuonPerformanceAlg::m_eventNumber
int m_eventNumber
Definition: MuonPerformanceAlg.h:73
MuonPerformanceAlg::m_nreco10
std::vector< int > m_nreco10
Definition: MuonPerformanceAlg.h:64
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
MuonPerformanceAlg::m_eventInfo
SG::ReadHandleKey< xAOD::EventInfo > m_eventInfo
Definition: MuonPerformanceAlg.h:66
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
xAOD::double
double
Definition: CompositeParticle_v1.cxx:159
xAOD::numberOfSCTHoles
@ numberOfSCTHoles
number of SCT holes [unit8_t].
Definition: TrackingPrimitives.h:270
MuonPerformanceAlg::m_truthMuons
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthMuons
Definition: MuonPerformanceAlg.h:67
Trk::Combined
@ Combined
Definition: TrackSummaryTool.h:32
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
MuonPerformanceAlg::m_ntruth
std::vector< int > m_ntruth
Definition: MuonPerformanceAlg.h:53
LikeEnum::Tight
@ Tight
Definition: LikelihoodEnums.h:15
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
MuonPerformanceAlg::m_nfound10
std::vector< int > m_nfound10
Definition: MuonPerformanceAlg.h:58
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
debug
const bool debug
Definition: MakeUncertaintyPlots.cxx:53
MuonPerformanceAlg::m_truthMuonRecoMuonLinkKey
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthMuonRecoMuonLinkKey
Definition: MuonPerformanceAlg.h:68
MuonPerformanceAlg::m_ntruth5
std::vector< int > m_ntruth5
Definition: MuonPerformanceAlg.h:54
MuonPerformanceAlg::m_ntruth10
std::vector< int > m_ntruth10
Definition: MuonPerformanceAlg.h:55
xAOD::numberOfTRTOutliers
@ numberOfTRTOutliers
number of TRT outliers [unit8_t].
Definition: TrackingPrimitives.h:276
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
a
TList * a
Definition: liststreamerinfos.cxx:10
h
Base_Fragment.width
width
Definition: Sherpa_i/share/common/Base_Fragment.py:59
Muons
Definition: Muons.py:1
LikeEnum::Medium
@ Medium
Definition: LikelihoodEnums.h:14
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
xAOD::numberOfSCTDeadSensors
@ numberOfSCTDeadSensors
number of dead SCT sensors crossed [unit8_t].
Definition: TrackingPrimitives.h:273
MuonPerformanceAlg::m_nfoundr
std::vector< int > m_nfoundr
Definition: MuonPerformanceAlg.h:59
MuonPerformanceAlg::m_truthMuonTruthOrigin
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthMuonTruthOrigin
Definition: MuonPerformanceAlg.h:70
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
extractSporadic.q
list q
Definition: extractSporadic.py:98
xAOD::numberOfSCTHits
@ numberOfSCTHits
number of hits in SCT [unit8_t].
Definition: TrackingPrimitives.h:268
AthAlgorithm::AthAlgorithm
AthAlgorithm()
Default constructor:
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
xAOD::numberOfPixelDeadSensors
@ numberOfPixelDeadSensors
number of dead pixel sensors crossed [unit8_t].
Definition: TrackingPrimitives.h:266
MuonPerformanceAlg::m_hitCutString
std::vector< std::string > m_hitCutString
Definition: MuonPerformanceAlg.h:52
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
CaloNoise_fillDB.mu
mu
Definition: CaloNoise_fillDB.py:53
PrepareReferenceFile.outfile
outfile
Definition: PrepareReferenceFile.py:42
MuonPerformanceAlg::m_truthMuonTruthType
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_truthMuonTruthType
Definition: MuonPerformanceAlg.h:69
MuonPerformanceAlg::m_nreco5
std::vector< int > m_nreco5
Definition: MuonPerformanceAlg.h:63
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
MuonPerformanceAlg::m_pdgsToBeConsidered
IntegerArrayProperty m_pdgsToBeConsidered
Definition: MuonPerformanceAlg.h:78
MuonPerformanceAlg::m_nfound5
std::vector< int > m_nfound5
Definition: MuonPerformanceAlg.h:57
fitman.k
k
Definition: fitman.py:528
MuonPerformanceAlg::passID
bool passID(const xAOD::TrackParticle *, bool debug) const
Definition: MuonPerformanceAlg.cxx:67
xAOD::numberOfInnermostPixelLayerHits
@ numberOfInnermostPixelLayerHits
these are the hits in the 0th pixel barrel layer
Definition: TrackingPrimitives.h:237
ServiceHandle< ICondSvc >
xAOD::numberOfPrecisionLayers
@ numberOfPrecisionLayers
layers with at least 3 hits [unit8_t].
Definition: TrackingPrimitives.h:288
Generated on Thu Nov 7 2024 21:52:54 for ATLAS Offline Software by   doxygen 1.8.18