ATLAS Offline Software
Loading...
Searching...
No Matches
Trig::TrigMuonMatching Class Reference

#include <TrigMuonMatching.h>

Inheritance diagram for Trig::TrigMuonMatching:
Collaboration diagram for Trig::TrigMuonMatching:

Classes

struct  EFmuon

Public Types

Definition of the StoreGate-like object's definition
typedef ServiceHandle< StoreGateSvcMetaStore_t
 Type of the metadata store object in Athena.
typedef const ServiceHandle< StoreGateSvc > & MetaStorePtr_t
 Type of the metadata store pointer in standalone mode.

Public Member Functions

 TrigMuonMatching (const std::string &name)
virtual ~TrigMuonMatching ()
virtual StatusCode initialize (void) override
 Dummy implementation of the initialisation function.
virtual Bool_t match (const xAOD::Muon *mu, std::string_view chain, const double mindelR=0.1) const override
virtual Bool_t matchL1 (const xAOD::Muon *mu, std::string_view l1item, const double DelR=0.2) const override
virtual Bool_t matchL2SA (const xAOD::Muon *mu, std::string_view l1item, std::string_view chain, const double DelR=0.2) const override
virtual Bool_t matchL2CB (const xAOD::Muon *mu, std::string_view chain, const double DelR=0.2) const override
virtual Double_t minDelR (const xAOD::Muon *mu, std::string_view chain, const double mindelR=0.1) const override
virtual Double_t minDelRL1 (const xAOD::Muon *mu, std::string_view l1item, const double DelR=0.2) const override
virtual Bool_t matchDimuon (const xAOD::Muon *mu1, const xAOD::Muon *mu2, const std::string &chain, std::pair< Bool_t, Bool_t > &result1, std::pair< Bool_t, Bool_t > &result2, const Double_t &mindelR=0.1) override
virtual Bool_t match (const double eta, const double phi, std::string_view chain, const double mindelR=0.1) const override
virtual Bool_t matchL1 (const double eta, const double phi, std::string_view l1item, const double DelR=0.2) const override
virtual Bool_t matchDimuon (const TLorentzVector &muon1, const TLorentzVector &muon2, const std::string &chain, std::pair< Bool_t, Bool_t > &result1, std::pair< Bool_t, Bool_t > &result2, const Double_t &mindelR=0.1) override
virtual Bool_t isPassedRerun (const std::string &trigger) const override
virtual StatusCode sysInitialize ()
 Function initialising the tool in the correct way in Athena.
virtual void print () const
 Print the state of the tool.
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const
Functions providing access to the input/output metadata
MetaStorePtr_t inputMetaStore () const
 Accessor for the input metadata store.
MetaStorePtr_t outputMetaStore () const
 Accessor for the output metadata store.
Additional helper functions, not directly mimicking Athena
template<class T>
const T * getProperty (const std::string &name) const
 Get one of the tool's properties.
const std::string & msg_level_name () const __attribute__((deprecated))
 A deprecated function for getting the message level's name.
const std::string & getName (const void *ptr) const
 Get the name of an object that is / should be in the event store.
SG::sgkey_t getKey (const void *ptr) const
 Get the (hashed) key of an object that is in the event store.

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.
Callback functions helping in metadata reading/writing
void setUseIncidents (const bool flag)
virtual void handle (const Incident &inc)
 Function receiving incidents from IncidentSvc/Event.
virtual StatusCode beginInputFile ()
 Function called when a new input file is opened.
virtual StatusCode endInputFile ()
 Function called when the currently open input file got completely processed.
virtual StatusCode beginEvent ()
 Function called when a new events is loaded.
virtual StatusCode metaDataStop ()
 Function called when the tool should write out its metadata.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

std::pair< bool, bool > matchDimuon (const TLorentzVector &muon1, const TLorentzVector &muon2, const DimuonChainInfo &chainInfo, const double mindelR)
double dR (const double eta1, const double phi1, const double eta2, const double phi2) const
int getL1pt (std::string_view l1item) const
Double_t matchedTrackDetail (EFmuon &efMuonId, const EFmuon &usedEFMuonId, const double eta, const double phi, const double mindelR, std::string_view chainEventTrigger) const
bool decodeDimuonChain (DimuonChainInfo &chainInfo)
bool isEqual (const double x, const double y) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

ToolHandle< Trig::TrigDecisionToolm_trigDecTool
std::map< std::string, DimuonChainInfom_DimuonChainMap
MetaStore_t m_inputMetaStore
 Object accessing the input metadata store.
MetaStore_t m_outputMetaStore
 Object accessing the output metadata store.
bool m_beginInputFileCalled
 Flag helping to discover when the tool misses the opening of the first input file.
bool m_useIncidents
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 19 of file TrigMuonMatching.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.

Constructor & Destructor Documentation

◆ TrigMuonMatching()

Trig::TrigMuonMatching::TrigMuonMatching ( const std::string & name)

Definition at line 26 of file TrigMuonMatching.cxx.

27 : asg::AsgMetadataTool( name ),
28 m_trigDecTool("Trig::TrigDecisionTool/TrigDecisionTool")
29 {
30 declareProperty( "TriggerTool", m_trigDecTool);
31 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
ToolHandle< Trig::TrigDecisionTool > m_trigDecTool

◆ ~TrigMuonMatching()

Trig::TrigMuonMatching::~TrigMuonMatching ( )
virtual

Definition at line 33 of file TrigMuonMatching.cxx.

33 {
34 }

Member Function Documentation

◆ beginEvent()

StatusCode asg::AsgMetadataTool::beginEvent ( )
protectedvirtualinherited

◆ 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, BookkeeperDumperTool, BookkeeperTool, PMGTools::PMGTruthWeightTool, TauAnalysisTools::TauEfficiencyCorrectionsTool, TauAnalysisTools::TauSmearingTool, Trig::TrigDecisionTool, Trig::xAODBunchCrossingTool, TrigConf::xAODConfigTool, xAODMaker::TriggerMenuMetaDataTool, and xAODMaker::TruthMetaDataTool.

Definition at line 185 of file AsgMetadataTool.cxx.

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

◆ declareGaudiProperty()

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

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

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ decodeDimuonChain()

bool Trig::TrigMuonMatching::decodeDimuonChain ( DimuonChainInfo & chainInfo)
private

Definition at line 392 of file TrigMuonMatching.cxx.

393 {
394 chainInfo.isValid = false;
395 std::map<std::string, DimuonChainInfo >::const_iterator p = m_DimuonChainMap.find(chainInfo.chain);
396 if (p != m_DimuonChainMap.end()) {
397 chainInfo = p->second;
398 return chainInfo.isValid;
399 }
400 std::vector<std::string> tokens = tokenize(chainInfo.chain, '_');
401 if (tokens.size() < 2) return false;
402 if (tokens[0] != "HLT") return false;
403 chainInfo.isSymmetric = (tokens[1].substr(0, 3) == "2mu");
404
405 if(chainInfo.isSymmetric) {
406 std::string threshold = std::string("HLT_" + tokens[1].substr(1));
407 chainInfo.thresholds.first = threshold;
408 chainInfo.thresholds.second = std::move(threshold);
409 chainInfo.isValid = true;
410 }
411 else {
412 if(tokens.size() != 3) return false;
413
414 std::string high = std::string("HLT_" + tokens[1]);
415 chainInfo.thresholds.first = std::move(high);
416 chainInfo.thresholds.second = chainInfo.chain;
417 chainInfo.isValid = true;
418 return chainInfo.isValid;
419
420 }
421 m_DimuonChainMap[chainInfo.chain] = chainInfo;
422 return chainInfo.isValid;
423 }
std::vector< std::string > tokenize(std::string_view the_str, std::string_view delimiters)
Splits the string into smaller substrings.
std::map< std::string, DimuonChainInfo > m_DimuonChainMap
std::pair< std::string, std::string > thresholds

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

◆ dR()

double Trig::TrigMuonMatching::dR ( const double eta1,
const double phi1,
const double eta2,
const double phi2 ) const
private

Definition at line 282 of file TrigMuonMatching.cxx.

286 {
287 double deta = fabs(eta1 - eta2);
288 double dphi = fabs(phi1 - phi2) < TMath::Pi() ? fabs(phi1 - phi2) : 2*TMath:: \
289 Pi() - fabs(phi1 - phi2);
290 return sqrt(deta*deta + dphi*dphi);
291 }
static const double Pi

◆ 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 BookkeeperDumperTool, BookkeeperTool, TrigConf::xAODConfigTool, xAODMaker::TriggerMenuMetaDataTool, and xAODMaker::TruthMetaDataTool.

Definition at line 193 of file AsgMetadataTool.cxx.

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

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ 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::Event 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::Event, we have a direct function call
123 // for this.
124 return evtStore()->event()->getKey( ptr );
125#else
126 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
127 return ( proxy == nullptr ? 0 : proxy->sgkey() );
128#endif // XAOD_STANDALONE
129 }
ServiceHandle< StoreGateSvc > & evtStore()

◆ getL1pt()

int Trig::TrigMuonMatching::getL1pt ( std::string_view l1item) const
private

Definition at line 294 of file TrigMuonMatching.cxx.

295 {
296 int rc = -1;
297
298 std::vector<std::string> tokens = tokenize(l1item, '_');
299 std::string pt;
300 if (tokens.size() == 1) {
301 pt = tokens.front();
302
303 } else if ((tokens.size() == 2 ) and (tokens.front() == "L1")) {
304 pt = tokens[1];
305
306 } else {
307 ATH_MSG_ERROR("TrigMuonMatching::getL1pt : cannot parse " << l1item);
308 return rc;
309 }
310
311 const size_t index = pt.find("MU");
312 if (index != std::string::npos) {
313 std::istringstream iss(pt.substr(index + 2));
314 iss >> rc;
315 } else {
316 ATH_MSG_ERROR("TrigMuonMatching::getL1pt : cannot parse " << l1item);
317 }
318
319 return rc;
320 }
#define ATH_MSG_ERROR(x,...)
static Double_t rc
str index
Definition DeMoScan.py:362

◆ 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::Event 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::Event, 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.

◆ handle()

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

Function receiving incidents from IncidentSvc/Event.

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 ) {
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 ) {
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 }
#define ATH_MSG_WARNING(x,...)
#define ATH_MSG_VERBOSE(x,...)
#define ATH_MSG_FATAL(x,...)
virtual StatusCode beginInputFile()
Function called when a new input file is opened.
virtual StatusCode beginEvent()
Function called when a new events is loaded.
bool m_beginInputFileCalled
Flag helping to discover when the tool misses the opening of the first input file.
virtual StatusCode endInputFile()
Function called when the currently open input file got completely processed.
virtual StatusCode metaDataStop()
Function called when the tool should write out its metadata.

◆ initialize()

StatusCode Trig::TrigMuonMatching::initialize ( void )
overridevirtual

Dummy implementation of the initialisation function.

It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...

Reimplemented from asg::AsgTool.

Definition at line 36 of file TrigMuonMatching.cxx.

36 {
37 // Greet the user:
38 ATH_MSG_INFO( "Initialising... " );
39
40 ATH_CHECK(m_trigDecTool.retrieve());
41
42 return StatusCode::SUCCESS;
43 }
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_INFO(x,...)

◆ 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 }
MetaStore_t m_inputMetaStore
Object accessing the input metadata store.

◆ isEqual()

bool Trig::TrigMuonMatching::isEqual ( const double x,
const double y ) const
private

Definition at line 425 of file TrigMuonMatching.cxx.

426 {
427 if (std::fabs(x - y) < std::numeric_limits<float>::epsilon()) return true;
428 return false;
429 }
#define y
#define x

◆ isPassedRerun()

Bool_t Trig::TrigMuonMatching::isPassedRerun ( const std::string & trigger) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 243 of file TrigMuonMatching.cxx.

244 {
245 const unsigned int bits = m_trigDecTool->isPassedBits(trigger);
246 if( (bits & TrigDefs::EF_passedRaw) && ! ( bits & TrigDefs::EF_passThrough) && ( bits & TrigDefs::EF_resurrected)){
247 return true;
248 }
249 return false;
250 }

◆ match() [1/2]

Bool_t Trig::TrigMuonMatching::match ( const double eta,
const double phi,
std::string_view chain,
const double mindelR = 0.1 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 107 of file TrigMuonMatching.cxx.

111 {
112 if(!m_trigDecTool->isPassed(chain)){
113 return false;
114 }
115 Double_t delmin = mindelR;
116 EFmuon efmuon,dummy;
117
118 matchedTrackDetail(efmuon, dummy, eta, phi, delmin, chain);
119 return efmuon.valid;
120 }
Scalar eta() const
pseudorapidity method
Scalar phi() const
phi method
Double_t matchedTrackDetail(EFmuon &efMuonId, const EFmuon &usedEFMuonId, const double eta, const double phi, const double mindelR, std::string_view chainEventTrigger) const

◆ match() [2/2]

Bool_t Trig::TrigMuonMatching::match ( const xAOD::Muon * mu,
std::string_view chain,
const double mindelR = 0.1 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 45 of file TrigMuonMatching.cxx.

48 {
49 return match(mu->eta(),mu->phi(),chain,mindelR);
50 }
virtual Bool_t match(const xAOD::Muon *mu, std::string_view chain, const double mindelR=0.1) const override

◆ matchDimuon() [1/3]

std::pair< Bool_t, Bool_t > Trig::TrigMuonMatching::matchDimuon ( const TLorentzVector & muon1,
const TLorentzVector & muon2,
const DimuonChainInfo & chainInfo,
const double mindelR )
private

Definition at line 252 of file TrigMuonMatching.cxx.

256 {
257
258 EFmuon trkId1, trkId2, dummy;
259 const Double_t dr1 = matchedTrackDetail(trkId1, dummy, muon1.Eta(), muon1.Phi(), mindelR, chainInfo.thresholds.first);
260
261 // for full scan trigger.
262
263 const Double_t dr2 = matchedTrackDetail(trkId2, dummy, muon2.Eta(), muon2.Phi(), mindelR, chainInfo.thresholds.second);
264
265 if(trkId1.valid && trkId2.valid &&
266 isEqual(trkId1.pt, trkId2.pt) &&
267 isEqual(trkId1.eta, trkId2.eta) &&
268 isEqual(trkId1.phi, trkId2.phi)){
269 if(dr1 > dr2){
270 matchedTrackDetail(trkId1, trkId2, muon1.Eta(), muon1.Phi(), mindelR, chainInfo.thresholds.first);
271 }
272 else{
273 matchedTrackDetail(trkId2, trkId1, muon2.Eta(), muon2.Phi(), mindelR, chainInfo.thresholds.second);
274 }
275 }
276
277 return std::pair<Bool_t, Bool_t> (trkId1.valid, trkId2.valid);
278 }
bool isEqual(const double x, const double y) const

◆ matchDimuon() [2/3]

Bool_t Trig::TrigMuonMatching::matchDimuon ( const TLorentzVector & muon1,
const TLorentzVector & muon2,
const std::string & chain,
std::pair< Bool_t, Bool_t > & result1,
std::pair< Bool_t, Bool_t > & result2,
const Double_t & mindelR = 0.1 )
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 214 of file TrigMuonMatching.cxx.

220 {
221 Double_t delmin = mindelR;
222 DimuonChainInfo chainInfo(chain);
223 if (not decodeDimuonChain(chainInfo)) {
224 ATH_MSG_ERROR("TrigMuonMatching : Failed to decode chain " << chain << " matchDimuon can accept only chains named HLT_2muXX and HLT_muXX_mu8noL1.");
225 return false;
226 }
227
228 std::pair<Bool_t, Bool_t> rc12, rc21;
229 rc12 = matchDimuon(muon1, muon2, chainInfo, delmin);
230 if(chainInfo.isSymmetric){
231 rc21.first = rc12.second; rc21.second = rc12.first;
232 }
233
234 else{
235 rc21 = matchDimuon(muon2, muon1, chainInfo, delmin);
236 }
237
238 result1.first = rc12.first; result1.second = rc21.second;
239 result2.first = rc21.first; result2.second = rc12.second;
240 return true;
241 }
bool decodeDimuonChain(DimuonChainInfo &chainInfo)
virtual Bool_t matchDimuon(const xAOD::Muon *mu1, const xAOD::Muon *mu2, const std::string &chain, std::pair< Bool_t, Bool_t > &result1, std::pair< Bool_t, Bool_t > &result2, const Double_t &mindelR=0.1) override

◆ matchDimuon() [3/3]

Bool_t Trig::TrigMuonMatching::matchDimuon ( const xAOD::Muon * mu1,
const xAOD::Muon * mu2,
const std::string & chain,
std::pair< Bool_t, Bool_t > & result1,
std::pair< Bool_t, Bool_t > & result2,
const Double_t & mindelR = 0.1 )
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 93 of file TrigMuonMatching.cxx.

99 {
100 TLorentzVector mu_1, mu_2;
101 mu_1.SetPtEtaPhiM(mu1->pt(),mu1->eta(),mu1->phi(),MUONMASS);
102 mu_2.SetPtEtaPhiM(mu2->pt(),mu2->eta(),mu2->phi(),MUONMASS);
103 Bool_t valid = matchDimuon(mu_1,mu_2,chain,result1,result2,mindelR);
104 return valid;
105 }
#define MUONMASS
virtual double pt() const override
The transverse momentum ( ) of the particle.
virtual double eta() const override
The pseudorapidity ( ) of the particle.
virtual double phi() const override
The azimuthal angle ( ) of the particle.
list valid
Definition calibdata.py:44

◆ matchedTrackDetail()

Double_t Trig::TrigMuonMatching::matchedTrackDetail ( EFmuon & efMuonId,
const EFmuon & usedEFMuonId,
const double eta,
const double phi,
const double mindelR,
std::string_view chainEventTrigger ) const
private

Definition at line 326 of file TrigMuonMatching.cxx.

332 {
333 efMuonId.valid = false;
334 Double_t drmin = mindelR;
335
336 auto cg = m_trigDecTool->getChainGroup(chainForEventTrigger);
337 auto fc = cg->features();
338
339#if defined(XAOD_STANDALONE) || defined(XAOD_ANALYSIS)
340 auto MuFeatureContainers = fc.containerFeature<xAOD::MuonContainer>();
341#else
342 const std::vector< Trig::Feature<xAOD::MuonContainer> > MuFeatureContainers = fc.get<xAOD::MuonContainer>();
343#endif
344
345 for(const auto & mucont : MuFeatureContainers){
346 for(auto mu : *mucont.cptr()){
347
348 // l1 matching
349
350 if(isEqual(usedEFMuonId.pt, mu->pt()) &&
351 isEqual(usedEFMuonId.eta, mu->eta()) &&
352 isEqual(usedEFMuonId.phi, mu->phi())) continue;
353
354 Double_t dr = dR(eta,phi,mu->eta(),mu->phi());
355 if(drmin > dr){
356 drmin = dr;
357 efMuonId.pt = mu->pt();
358 efMuonId.eta = mu->eta();
359 efMuonId.phi = mu->phi();
360 efMuonId.valid = true;
361 }
362 }
363 }
364 return drmin;
365 }
const T * get(size_type n) const
Access an element, as an rvalue.
double dR(const double eta1, const double phi1, const double eta2, const double phi2) const
MuonContainer_v1 MuonContainer
Definition of the current "Muon container version".

◆ matchL1() [1/2]

Bool_t Trig::TrigMuonMatching::matchL1 ( const double eta,
const double phi,
std::string_view l1item,
const double DelR = 0.2 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 122 of file TrigMuonMatching.cxx.

126 {
127 if(!m_trigDecTool->isPassed("L1_MU.*")){
128 return false;
129 }
130 Double_t l1dr = DelR;
131 const xAOD::MuonRoIContainer* muonrois = 0;
132 StatusCode sc = evtStore()->retrieve(muonrois,"LVL1MuonRoIs");
133 if(!sc){
134 ATH_MSG_ERROR("TrigMuonMatching::matchL1 : could not retrieve LVL1MuonRoIs");
135 return false;
136 }
137 Int_t threshold = getL1pt(l1item);
138 xAOD::MuonRoIContainer::const_iterator muroi_itr = muonrois->begin();
139 xAOD::MuonRoIContainer::const_iterator muroi_end = muonrois->end();
140 for( ; muroi_itr != muroi_end; ++muroi_itr ) {
141 if((*muroi_itr)->thrValue() >= threshold*1000){
142 Double_t dR = TrigMuonMatching::dR(eta, phi, (*muroi_itr)->eta(), (*muroi_itr)->phi());
143 if(dR < l1dr){
144 return true;
145 }
146 }
147 }
148 return false;
149 }
static Double_t sc
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:861
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
int getL1pt(std::string_view l1item) const
::StatusCode StatusCode
StatusCode definition for legacy code.
MuonRoIContainer_v1 MuonRoIContainer

◆ matchL1() [2/2]

Bool_t Trig::TrigMuonMatching::matchL1 ( const xAOD::Muon * mu,
std::string_view l1item,
const double DelR = 0.2 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 52 of file TrigMuonMatching.cxx.

55 {
56 return matchL1(mu->eta(),mu->phi(),l1item,DelR);
57 }
virtual Bool_t matchL1(const xAOD::Muon *mu, std::string_view l1item, const double DelR=0.2) const override

◆ matchL2CB()

Bool_t Trig::TrigMuonMatching::matchL2CB ( const xAOD::Muon * mu,
std::string_view chain,
const double DelR = 0.2 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 193 of file TrigMuonMatching.cxx.

196 {
197 auto cg = m_trigDecTool->getChainGroup(chain);
198 auto fc = cg->features(TrigDefs::alsoDeactivateTEs);
199#if defined(XAOD_STANDALONE) || defined(XAOD_ANALYSIS)
200 auto MuFeatureContainers = fc.containerFeature<xAOD::L2CombinedMuonContainer>("",TrigDefs::alsoDeactivateTEs);
201#else
202 const std::vector< Trig::Feature<xAOD::L2CombinedMuonContainer> > MuFeatureContainers = fc.get<xAOD::L2CombinedMuonContainer>("", TrigDefs::alsoDeactivateTEs);
203#endif
204
205 for(const auto & mucont : MuFeatureContainers){
206 for(auto muon : *mucont.cptr()){
207 Double_t dR = TrigMuonMatching::dR(mu->eta(), mu->phi(), muon->eta(), muon->phi());
208 if(dR < DelR) return true;
209 }
210 }
211 return false;
212 }
L2CombinedMuonContainer_v1 L2CombinedMuonContainer
Define the latest version of the muon CB container.

◆ matchL2SA()

Bool_t Trig::TrigMuonMatching::matchL2SA ( const xAOD::Muon * mu,
std::string_view l1item,
std::string_view chain,
const double DelR = 0.2 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 151 of file TrigMuonMatching.cxx.

155 {
156 if(!m_trigDecTool->isPassed("L1_MU.*")){
157 return false;
158 }
159 const xAOD::MuonRoIContainer* muonrois = 0;
160 StatusCode sc = evtStore()->retrieve(muonrois,"LVL1MuonRoIs");
161 if(!sc){
162 ATH_MSG_ERROR("TrigMuonMatching::matchL1 : could not retrieve LVL1MuonRoIs");
163 return false;
164 }
165 xAOD::MuonRoIContainer::const_iterator muroi_itr = muonrois->begin();
166 xAOD::MuonRoIContainer::const_iterator muroi_end = muonrois->end();
167 Int_t threshold = getL1pt(l1item);
168 unsigned int ROI = 0;
169 const std::string emptyStr;
170 for( ; muroi_itr != muroi_end; ++muroi_itr ) {
171 if(!((*muroi_itr)->thrValue() >= threshold*1000)) continue;
172 ROI = (*muroi_itr)->getRoI();
173 auto cg = m_trigDecTool->getChainGroup(chain);
174 auto fc = cg->features(TrigDefs::alsoDeactivateTEs);
175#if defined(XAOD_STANDALONE) || defined(XAOD_ANALYSIS)
176 auto MuFeatureContainers = fc.containerFeature<xAOD::L2StandAloneMuonContainer>(emptyStr,TrigDefs::alsoDeactivateTEs);
177#else
178 const std::vector< Trig::Feature<xAOD::L2StandAloneMuonContainer> > MuFeatureContainers = fc.get<xAOD::L2StandAloneMuonContainer>(emptyStr, TrigDefs::alsoDeactivateTEs);
179#endif
180 for(const auto & mucont : MuFeatureContainers){
181 for(auto muon : *mucont.cptr()){
182 if(muon->roiNumber() == ROI){
183 Double_t dR = TrigMuonMatching::dR(mu->eta(), mu->phi(), muon->eta(), muon->phi());
184 if(dR < DelR) return true;
185 }
186 }
187 }
188 }
189
190 return false;
191 }
L2StandAloneMuonContainer_v2 L2StandAloneMuonContainer
Define the latest version of the muon SA container.

◆ 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 BookkeeperDumperTool, BookkeeperTool, xAODMaker::TriggerMenuMetaDataTool, and xAODMaker::TruthMetaDataTool.

Definition at line 209 of file AsgMetadataTool.cxx.

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

◆ minDelR()

Double_t Trig::TrigMuonMatching::minDelR ( const xAOD::Muon * mu,
std::string_view chain,
const double mindelR = 0.1 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 59 of file TrigMuonMatching.cxx.

62 {
63 Double_t delmin = mindelR;
64 EFmuon efmuon,dummy;
65 return matchedTrackDetail(efmuon, dummy, mu->eta(), mu->phi(), delmin, chain);
66 }

◆ minDelRL1()

Double_t Trig::TrigMuonMatching::minDelRL1 ( const xAOD::Muon * mu,
std::string_view l1item,
const double DelR = 0.2 ) const
overridevirtual

Implements Trig::ITrigMuonMatching.

Definition at line 68 of file TrigMuonMatching.cxx.

71 {
72 Double_t l1dr = DelR;
73 const xAOD::MuonRoIContainer* muonrois = 0;
74 StatusCode sc = evtStore()->retrieve(muonrois,"LVL1MuonRoIs");
75 if(!sc){
76 ATH_MSG_ERROR("TrigMuonMatching::matchL1 : could not retrieve LVL1MuonRoIs");
77 return false;
78 }
79 Int_t threshold = getL1pt(l1item);
80 xAOD::MuonRoIContainer::const_iterator muroi_itr = muonrois->begin();
81 xAOD::MuonRoIContainer::const_iterator muroi_end = muonrois->end();
82 for( ; muroi_itr != muroi_end; ++muroi_itr ) {
83 if((*muroi_itr)->thrValue() >= threshold*1000){
84 Double_t dR = TrigMuonMatching::dR(mu->eta(), mu->phi(), (*muroi_itr)->eta(), (*muroi_itr)->phi());
85 if(dR < l1dr){
86 l1dr = dR;
87 }
88 }
89 }
90 return l1dr;
91 }

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const
inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns
The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

101 {
102
103 return MSG::name( msg().level() );
104 }
MsgStream & msg() const
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition MsgLevel.cxx:19

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

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

◆ 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 }
MetaStore_t m_outputMetaStore
Object accessing the output metadata store.

◆ print()

◆ renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setUseIncidents()

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

Definition at line 132 of file AsgMetadataTool.h.

133 {
135 }
bool flag
Definition master.py:29

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

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ 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_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_DimuonChainMap

std::map<std::string, DimuonChainInfo> Trig::TrigMuonMatching::m_DimuonChainMap
private

Definition at line 97 of file TrigMuonMatching.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_trigDecTool

ToolHandle<Trig::TrigDecisionTool> Trig::TrigMuonMatching::m_trigDecTool
private

Definition at line 95 of file TrigMuonMatching.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: