ATLAS Offline Software
Loading...
Searching...
No Matches
T_AnalysisConfig_Tier0.h
Go to the documentation of this file.
1/* emacs: this is -*- c++ -*- */
22
23#ifndef TrigInDetAnalysisExample_T_AnalysisConfig_Tier0_H
24#define TrigInDetAnalysisExample_T_AnalysisConfig_Tier0_H
25
26
27#include "TrigInDetAnalysis/TIDAEvent.h"
28#include "TrigInDetAnalysis/TIDAVertex.h"
29#include "TrigInDetAnalysisUtils/T_AnalysisConfig.h"
30
31#include "TrigInDetAnalysisExample/Analysis_Tier0.h"
32#include "TrigInDetAnalysisExample/VtxAnalysis.h"
33#include "TrigInDetAnalysisExample/ChainString.h"
34#include "TrigInDetAnalysisExample/TIDATools.h"
35
36#include "TTree.h"
37#include "TFile.h"
38
39
40// McParticleEvent includes
41#include "McParticleEvent/TruthParticleContainer.h"
42
43#include "GeneratorObjects/McEventCollection.h"
44#include "AtlasHepMC/GenEvent.h"
45#include "AtlasHepMC/GenVertex.h"
46#include "AtlasHepMC/GenParticle.h"
47
48#include "EventInfo/EventInfo.h"
49#include "EventInfo/EventID.h"
50#include "xAODEventInfo/EventInfo.h"
51
52
53
54#include "TrigInDetAnalysis/TIDDirectory.h"
55
56#include "TrigInDetAnalysis/Filter_AcceptAll.h"
57
58#include "TrigInDetAnalysisUtils/TIDARoiDescriptorBuilder.h"
59#include "TrigInDetAnalysisUtils/Filter_etaPT.h"
60#include "TrigInDetAnalysisUtils/Filter_RoiSelector.h"
61#include "TrigInDetAnalysisUtils/Associator_BestMatch.h"
62#include "TrigInDetAnalysisUtils/Filters.h"
63
64
65#include "VxVertex/VxContainer.h"
66
67#include "muonEvent/MuonContainer.h"
68
69#include "egammaEvent/ElectronContainer.h"
70
71#include "tauEvent/TauJetContainer.h"
72
73
74#include "TrigSteeringEvent/HLTResult.h"
75#include "TrigDecisionTool/ExpertMethods.h"
76
77#include "TrigSteeringEvent/TrigRoiDescriptorCollection.h"
78
79#include "xAODTracking/TrackParticle.h"
80#include "xAODTracking/TrackParticleContainer.h"
81
82
83
84
85
86template<typename T>
87class T_AnalysisConfig_Tier0 : public T_AnalysisConfig<T> {
88
89public:
90
91 // Full constructor: test/reference/selection
92 // - analysisInstanceName: the name of the analysis chain being created
93 // - xxxChainName: the name of the chain to be used as test/reference/selection; must be "StoreGate" in case of direct access to SG containers
94 // - xxxType: the type of tracks to be retrieved from the test/reference/selection chain or container
95 // - xxxKey: the key for tracks to be retrieved from the test/reference/selection chain or container
96 // - all standard operations are performed in loops over 0=test 1=reference 2=selection
97 T_AnalysisConfig_Tier0(const std::string& analysisInstanceName,
98 const std::string& testChainName, const std::string& testType, const std::string& testKey,
99 const std::string& referenceChainName, const std::string& referenceType, const std::string& referenceKey,
100 TrackFilter* testFilter, TrackFilter* referenceFilter,
101 TrackAssociator* associator,
102 TrackAnalysis* analysis) :
103 T_AnalysisConfig<T>( analysisInstanceName,
104 testChainName, testType, testKey,
105 referenceChainName, referenceType, referenceKey,
106 testFilter, referenceFilter,
107 associator,
108 analysis ),
109 _analysis(0),
110 m_useBeamCondSvc(false),
111 m_doOffline(true),
112 m_doMuons(false),
113 m_doElectrons(false),
114 m_doTaus(false),
115 m_doBjets(false),
116 m_hasTruthMap(false),
117 m_doTauThreeProng(false),
118 m_tauEtCutOffline(false),
119 m_NRois(0),
120 m_NRefTracks(0),
121 m_NTestTracks(0),
122 m_runPurity(false),
123 m_shifter(false),
124 m_pTthreshold(0),
125 m_first(true),
126 m_containTracks(false)
127 {
128 m_event = new TIDA::Event();
129 m_chainNames.push_back(testChainName);
130
131#if 0
132 ChainString& chain = m_chainNames.back();
133
134 std::cout << "\nT_AnalysisConfig_Tier0::name: " << name() << "\t" << this << std::endl;
135 std::cout << "T_AnalysisConfig_Tier0::chain specification: " << testChainName << " -> " << chain << "\t" << chain.raw() << std::endl;
136 std::cout << "\tchain: " << chain.head() << std::endl;
137 std::cout << "\tkey: " << chain.tail() << std::endl;
138 std::cout << "\troi: " << chain.roi() << std::endl;
139 std::cout << "\tvtx: " << chain.vtx() << std::endl;
140 std::cout << "\tte: " << chain.element() << std::endl;
141
142 std::cout << "\tpost: " << chain.post() << std::endl;
143 std::cout << "\tpt: " << chain.postvalue("pt") << std::endl;
144
145 std::cout << "\tcontainTracks: " << m_containTracks << std::endl;
146#endif
147
148 m_testType = testType;
149 }
150
151 virtual ~T_AnalysisConfig_Tier0() { delete m_event; }
152
153 void setRunPurity( bool b ) { m_runPurity=b; }
154
155 void setShifter( bool b ) { m_shifter=b; }
156
157 void useBeamCondSvc( bool b ) { m_useBeamCondSvc = b; }
158
159 void containTracks( bool b ) { m_containTracks = b; }
160
161public:
162
163 Analysis_Tier0* _analysis;
164
166 // void verbose( std::ostream& s ) {
167 // if( m_provider->msg().level() <= MSG::VERBOSE ) m_provider->msg(MSG::VERBOSE) << s << endmsg;
168 // }
169
170protected:
171
172 using T_AnalysisConfig<T>::m_provider;
173 using T_AnalysisConfig<T>::m_tdt;
174 using T_AnalysisConfig<T>::m_mcTruth;
175
176 using T_AnalysisConfig<T>::name;
177 using T_AnalysisConfig<T>::m_analysis;
178
179 using T_AnalysisConfig<T>::m_selectorTest;
180 using T_AnalysisConfig<T>::m_selectorRef;
181 using T_AnalysisConfig<T>::m_associator;
182 using T_AnalysisConfig<T>::m_filters;
183
184 // using T_AnalysisConfig<T>::selectTracks<TrigInDetTrackCollection>;
185
186 // using T_AnalysisConfig<T>::selectTracks;
187
188 virtual void loop() {
189
190 if( m_provider->msg().level() <= MSG::VERBOSE) {
191 m_provider->msg(MSG::VERBOSE) << "AnalysisConfig_Tier0::loop() for " << T_AnalysisConfig<T>::m_analysisInstanceName << endmsg;
192 }
193
194 // get (offline) beam position
195 double xbeam = 0;
196 double ybeam = 0;
197
198 if ( m_first ) {
199
200 m_first = false;
201
202 m_provider->msg(MSG::VERBOSE) << " using beam position\tx=" << xbeam << "\ty=" << ybeam << endmsg;
203
204 if (m_provider->msg().level() <= MSG::VERBOSE) {
205
206 std::vector<std::string> configuredChains = (*(m_tdt))->getListOfTriggers("L2_.*, EF_.*, HLT_.*");
207
208 for ( unsigned i=0 ; i<configuredChains.size() ; i++ ) {
209 m_provider->msg(MSG::VERBOSE) << "Chain " << configuredChains[i] << endmsg;
210 }
211
212 }
213
214
215 std::vector<ChainString>::iterator chainitr = m_chainNames.begin();
216
217 std::vector<ChainString> chains;
218
222 while ( chainitr!=m_chainNames.end() ) {
223
225 ChainString& chainName = (*chainitr);
226
227 m_provider->msg(MSG::INFO) << "process chain " << chainName << endmsg;
228
229 if ( chainName.head() == "" ) {
230
231 std::string selectChain = chainName.raw();
232
233 chains.push_back( ChainString(selectChain) );
234
235 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
236 m_provider->msg(MSG::VERBOSE) << "Matching chain " << selectChain << " (" << chainName.head() << ")" << endmsg;
237 }
238
239 }
240 else {
241
243 std::vector<std::string> selectChains = (*(m_tdt))->getListOfTriggers( chainName.head() );
244
245
246 for ( unsigned iselected=0 ; iselected<selectChains.size() ; iselected++ ) {
247
248 selectChains[iselected] = chainName.subs( selectChains[iselected] );
249
250#if 0
251 std::cout << "sorting:: chain specification: " << chainName << "\traw:" << chainName.raw() << std::endl;
252 std::cout << "\tchain: " << chainName.head() << std::endl;
253 std::cout << "\tkey: " << chainName.tail() << std::endl;
254 std::cout << "\troi: " << chainName.roi() << std::endl;
255 std::cout << "\tvtx: " << chainName.vtx() << std::endl;
256 std::cout << "\tte: " << chainName.element() << std::endl;
257 std::cout << "\tind: " << chainName.extra() << std::endl;
258#endif
259
261 chains.push_back( ChainString(selectChains[iselected]) );
262
263 if(m_provider->msg().level() <= MSG::VERBOSE) {
264 m_provider->msg(MSG::VERBOSE) << "Matching chain " << selectChains[iselected] << " (" << chainName.head() << ")" << endmsg;
265 }
266
267 }
268 }
269
270 ++chainitr;
271 }
272
273 // m_chainNames.insert( m_chainNames.end(), chains.begin(), chains.end() );
274 m_chainNames = chains;
275
276 for ( unsigned ic=0 ; ic<m_chainNames.size() ; ic++ ) m_provider->msg(MSG::VERBOSE) << "Analyse chain " << m_chainNames[ic] << endmsg;
277
278 }
279
280
282 Filter_True filter;
283
284 Filter_etaPT filter_etaPT(5,200);
285 Filter_Combined filter_truth( &filter_etaPT, &filter_etaPT);
286
288 // Filter_Combined filterRef (&filter_offline, &filter_vertex);
289
290 int iRefFilter = 1;
291 int iTestFilter = 0;
292
293 if ( m_runPurity ) {
294 iRefFilter = 0;
295 iTestFilter = 1;
296 }
297
298 Filter_Combined filterRef( m_filters[iRefFilter][0], &filter );
299 Filter_Combined filterTest( m_filters[iTestFilter][0], &filter );
300
301 TrigTrackSelector selectorTruth( &filter_truth );
302 TrigTrackSelector selectorRef( &filterRef );
303 m_selectorRef = &selectorRef;
304 TrigTrackSelector selectorTest( &filterTest );
305 m_selectorTest = &selectorTest;
306
307 if ( xbeam!=0 || ybeam!=0 ) {
308 m_selectorRef->setBeamline( xbeam, ybeam );
309 }
310
312
313 // clear the ntuple TIDA::Event class
314 m_event->clear();
315
317#ifndef XAODTRACKING_TRACKPARTICLE_H
318 const EventInfo* pEventInfo;
319#else
320 const xAOD::EventInfo* pEventInfo;
321#endif
322 unsigned run_number = 0;
323 uint64_t event_number = 0;
324 unsigned lumi_block = 0;
325 unsigned bunch_crossing_id = 0;
326 unsigned time_stamp = 0;
327 double mu_val = 0;
328
329 if ( this->template retrieve(pEventInfo, "EventInfo").isFailure() ) {
330 m_provider->msg(MSG::WARNING) << "Failed to get EventInfo " << endmsg;
331 } else {
332
333#ifndef XAODTRACKING_TRACKPARTICLE_H
334 run_number = pEventInfo->event_ID()->run_number();
335 event_number = pEventInfo->event_ID()->event_number();
336 lumi_block = pEventInfo->event_ID()->lumi_block();
337 time_stamp = pEventInfo->event_ID()->time_stamp();
338 bunch_crossing_id = pEventInfo->event_ID()->bunch_crossing_id();
339 mu_val = pEventInfo->averageInteractionsPerCrossing();
340#else
341 run_number = pEventInfo->runNumber();
342 event_number = pEventInfo->eventNumber();
343 lumi_block = pEventInfo->lumiBlock();
344 time_stamp = pEventInfo->timeStamp();
345 bunch_crossing_id = pEventInfo->bcid();
346 mu_val = pEventInfo->averageInteractionsPerCrossing();
347#endif
348 }
349
350 if(m_provider->msg().level() <= MSG::VERBOSE){
351 m_provider->msg(MSG::VERBOSE) << "run " << run_number
352 << "\tevent " << event_number
353 << "\tlb " << lumi_block << endmsg;
354 }
355
356 // m_provider->msg(MSG::INFO) << "run " << run_number
357 // << "\tevent " << event_number
358 // << "\tlb " << lumi_block << endmsg;
359
360 // std::cout << "run " << run_number << "\tevent " << event_number << "\tlb " << lumi_block << std::endl;
361
362
363 // clear the ntuple TIDA::Event class
364 m_event->clear();
365
366 m_event->run_number(run_number);
367 m_event->event_number(event_number);
368 m_event->lumi_block(lumi_block);
369 m_event->time_stamp(time_stamp);
370 m_event->bunch_crossing_id(bunch_crossing_id);
371 m_event->mu(mu_val);
372
375
376 bool analyse = false;
377
378 // Check HLTResult
379
380 for ( unsigned ichain=0 ; ichain<m_chainNames.size() ; ichain++ ) {
381
382 const std::string& chainname = m_chainNames[ichain].head();
383
384 if ( chainname == "" ) analyse = true;
385 else {
386
387 //Only for trigger chains
388 if ( chainname.find("L2") == std::string::npos &&
389 chainname.find("EF") == std::string::npos &&
390 chainname.find("HLT") == std::string::npos ) continue;
391
392 if ( m_provider->msg().level() <= MSG::DEBUG ) {
393 m_provider->msg(MSG::DEBUG) << "Chain " << chainname
394 << "\tpass " << (*m_tdt)->isPassed(chainname)
395 << "\tpres " << (*m_tdt)->getPrescale(chainname) << endmsg;
396 }
397
398 // std::cout << "Chain " << chainname << "\tpass " << (*m_tdt)->isPassed(chainname)
399 // << "\tpres " << (*m_tdt)->getPrescale(chainname) << std::endl;
400
401 if ( (*(m_tdt))->isPassed(chainname) || (*(m_tdt))->getPrescale(chainname) ) analyse = true;
402
403 }
404 }
405
406 // Remove this code to skip on truncated HLT results
407 // Need to leave the code here for the time being however, since we will
408 // still need in the future a more robust test to achieve this same
409 // functionality
410 //
411 // if ( (*m_tdt)->ExperimentalAndExpertMethods().isHLTTruncated() ) {
412 // m_provider->msg(MSG::WARNING) << "HLTResult truncated, skipping event" << endmsg;
413 // return;
414 // }
415
416 if ( !this->m_keepAllEvents && !analyse ) {
417 // m_provider->msg(MSG::VERBOSE) << "No chains passed unprescaled - not processing this event" << endmsg;
418 if(m_provider->msg().level() <= MSG::VERBOSE)
419 m_provider->msg(MSG::VERBOSE) << "No chains passed unprescaled - not processing this event" << endmsg;
420 return;
421 }
422
424
425 selectorTruth.clear();
426
427 if(m_provider->msg().level() <= MSG::VERBOSE)
428 m_provider->msg(MSG::VERBOSE) << "MC Truth flag " << m_mcTruth << endmsg;
429
430 const TrigInDetTrackTruthMap* truthMap = 0;
431
432 if ( m_mcTruth ) {
433 if(m_provider->msg().level() <= MSG::VERBOSE ) m_provider->msg(MSG::VERBOSE) << "getting Truth" << endmsg;
434
435 if ( this->template retrieve(truthMap, "TrigInDetTrackTruthMap").isFailure()) {
436 if(m_provider->msg().level() <= MSG::VERBOSE)
437 m_provider->msg(MSG::VERBOSE) << "TrigInDetTrackTruthMap not found" << endmsg;
438 m_hasTruthMap = false;
439 }
440 else {
441 if(m_provider->msg().level() <= MSG::VERBOSE)
442 m_provider->msg(MSG::VERBOSE) << "TrigInDetTrackTruthMap found" << endmsg;
443 m_hasTruthMap = true;
444 }
445 }
446
447
449
450 std::vector<TIDA::Vertex> vertices;
451 std::vector<TIDA::Vertex> vertices_rec;
452
453 std::vector<double> refbeamspot;
454 std::vector<double> testbeamspot;
455
456
457
458#ifndef XAODTRACKING_TRACKPARTICLE_H
459
460 const VxContainer* primaryVtxCollection;
461
462 if ( m_doOffline ) {
463 if ( m_provider->evtStore()->template contains<VxContainer>("VxPrimaryCandidate") ) {
464 if ( this->template retrieve(primaryVtxCollection, "VxPrimaryCandidate").isFailure()) {
465 if (m_provider->msg().level() <= MSG::WARNING) m_provider->msg(MSG::WARNING) << "Primary vertex container not found" << endmsg;
466 }
467 else {
468 VxContainer::const_iterator vtxitr = primaryVtxCollection->begin();
469 for ( ; vtxitr != primaryVtxCollection->end(); ++vtxitr) {
470 if ( (*vtxitr)->vxTrackAtVertex()->size()>0 ) {
471 vertices.push_back( TIDA::Vertex( (*vtxitr)->recVertex().position().x(),
472 (*vtxitr)->recVertex().position().y(),
473 (*vtxitr)->recVertex().position().z(),
474 0,0,0,
475 (*vtxitr)->vxTrackAtVertex()->size() ) );
476 }
477 }
478 }
479
480 // filter_vertex.setVertex(vertices);
481 }
482 }
483
484#else
485
486 //std::vector<TIDA::Vertex> vertices;
487
488 m_provider->msg(MSG::VERBOSE) << "fetching AOD Primary vertex container" << endmsg;
489
490 const xAOD::VertexContainer* xaodVtxCollection = 0;
491
492 if ( this->template retrieve( xaodVtxCollection, "PrimaryVertices" ).isFailure()) {
493 if (m_provider->msg().level() <= MSG::WARNING) m_provider->msg(MSG::WARNING) << "xAOD Primary vertex container not found with key " << "PrimaryVertices" << endmsg;
494 }
495
496 if ( xaodVtxCollection!=0 ) {
497
498 m_provider->msg(MSG::VERBOSE) << "xAOD Primary vertex container " << xaodVtxCollection->size() << " entries" << endmsg;
499
500 xAOD::VertexContainer::const_iterator vtxitr = xaodVtxCollection->begin();
501 for ( ; vtxitr != xaodVtxCollection->end(); ++vtxitr ) {
502 if ( (*vtxitr)->nTrackParticles()>0 && (*vtxitr)->vertexType()!=0 ) {
503 vertices.push_back( TIDA::Vertex( (*vtxitr)->x(),
504 (*vtxitr)->y(),
505 (*vtxitr)->z(),
507 (*vtxitr)->covariancePosition()(Trk::x,Trk::x),
508 (*vtxitr)->covariancePosition()(Trk::y,Trk::y),
509 (*vtxitr)->covariancePosition()(Trk::z,Trk::z),
510 (*vtxitr)->nTrackParticles(),
512 (*vtxitr)->chiSquared(),
513 (*vtxitr)->numberDoF() ) );
514 }
515 }
516 }
517
518
519#endif
520
521
523
524 if ( m_mcTruth ) {
525 m_event->addChain( "Truth" );
526 m_event->back().addRoi(TIDARoiDescriptor());
527 m_event->back().back().addTracks(selectorTruth.tracks());
528 }
529
531
532 if ( m_doOffline ) {
533 for ( unsigned i=0 ; i<vertices.size() ; i++ ) {
534 if(m_provider->msg().level() <= MSG::VERBOSE)
535 m_provider->msg(MSG::VERBOSE) << "vertex " << i << " " << vertices[i] << endmsg;
536 m_event->addVertex(vertices[i]);
537 }
538 }
539
541
542 // int Noff = 0;
543 std::vector<TIDA::Track*> offline_tracks;
544 std::vector<TIDA::Track*> electron_tracks;
545 std::vector<TIDA::Track*> muon_tracks;
546
547 std::vector<TIDA::Track*> ref_tracks;
548 std::vector<TIDA::Track*> test_tracks;
549
550 offline_tracks.clear();
551 electron_tracks.clear();
552 muon_tracks.clear();
553
554 ref_tracks.clear();
555 test_tracks.clear();
556
558 for ( unsigned ichain=0 ; ichain<m_chainNames.size() ; ichain++ ) {
559
560 test_tracks.clear();
561
563
564 // std::string& chainname = chains[ichain];
565 const std::string& chainname = m_chainNames[ichain].head();
566 const std::string& key = m_chainNames[ichain].tail();
567 const std::string& vtx_name = m_chainNames[ichain].vtx();
568 //no currently used but retained in case
569 //const std::string& roi_name = m_chainNames[ichain].roi();
570 //const std::string& te_name = m_chainNames[ichain].element();
571
572 m_pTthreshold = 0;
573
574 if ( m_chainNames[ichain].postcount() ) {
575 std::string ptvalue = m_chainNames[ichain].postvalue("pt");
576 if ( ptvalue!="" ) m_pTthreshold = std::stod(ptvalue);
577 }
578
579 // std::cout << "\tchain " << m_chainNames[ichain] << "\tchainname " << chainname << "\tvtx " << vtx_name << "\troi " << roi_name << std::endl;
580
581 unsigned _decisiontype = TrigDefs::Physics;
582 unsigned decisiontype;
583
584 if ( this->requireDecision() ) _decisiontype = TrigDefs::requireDecision;
585
586
587 if ( m_chainNames[ichain].passed() ) decisiontype = _decisiontype;
588 else decisiontype = TrigDefs::alsoDeactivateTEs;
589
590 // if ( decisiontype==TrigDefs::requireDecision ) std::cout << "\tSUTT TrigDefs::requireDecision " << decisiontype << std::endl;
591 // if ( decisiontype==TrigDefs::Physics ) std::cout << "\tSUTT TrigDefs::Physics " << decisiontype << std::endl;
592
593
595 const std::string& key_index_string = m_chainNames[ichain].extra();
596 unsigned key_index = 0;
597 if ( key_index_string!="" ) key_index = std::atoi( key_index_string.c_str() );
598
599 if ( chainname!="" && m_provider->msg().level() <= MSG::VERBOSE ) {
600
601 m_provider->msg(MSG::VERBOSE) << "status for chain " << chainname
602 << "\tpass " << (*m_tdt)->isPassed(chainname)
603 << "\tprescale " << (*m_tdt)->getPrescale(chainname) << endmsg;
604
605 m_provider->msg(MSG::VERBOSE) << "fetching features for chain " << chainname << endmsg;
606
607 m_provider->msg(MSG::VERBOSE) << chainname << "\tpassed: " << (*m_tdt)->isPassed( chainname ) << endmsg;
608 }
609
610 // std::cout << "\tstatus for chain " << chainname
611 // << "\tpass " << (*m_tdt)->isPassed( chainname )
612 // << "\tpassdt " << (*m_tdt)->isPassed( chainname, decisiontype )
613 // << "\tprescale " << (*m_tdt)->getPrescale( chainname ) << std::endl;
614
615
616 // m_provider->msg(MSG::INFO) << chainname << "\tpassed: " << (*m_tdt)->isPassed( chainname ) << "\t" << m_chainNames[ichain] << "\trun " << run_number << "\tevent " << event_number << endmsg;
617
618
619 if ( chainname!="" && !this->m_keepAllEvents && !(*m_tdt)->isPassed( chainname, decisiontype ) ) continue;
620
623 // Trig::FeatureContainer f = (*m_tdt)->features( chainname, TrigDefs::alsoDeactivateTEs);
624
627
628 ChainString& chainConfig = m_chainNames[ichain];
629
630 std::string chainName = chainConfig.head();
631
632 m_event->addChain( chainConfig );
633
634 TIDA::Chain& chain = m_event->back();
635
636 if ( chainName == "" ) {
637
638 m_selectorTest->clear();
639
641
642 TIDARoiDescriptor* roiInfo = new TIDARoiDescriptor(true);
643
644 chain.addRoi( *roiInfo );
645
646
647# ifdef XAODTRACKING_TRACKPARTICLE_H
648 if ( m_provider->evtStore()->template contains<xAOD::TrackParticleContainer>(key) ) {
649 this->template selectTracks<xAOD::TrackParticleContainer>( m_selectorTest, key );
650 refbeamspot = this->template getBeamspot<xAOD::TrackParticleContainer>( key );
651 }
652# endif
653
654 const std::vector<TIDA::Track*>& testtracks = m_selectorTest->tracks();
655
656 chain.back().addTracks(testtracks);
657
658 delete roiInfo;
659
660 }
661 else {
662
663 Trig::FeatureContainer f = (*m_tdt)->features( chainname, decisiontype );
664 Trig::FeatureContainer::combination_const_iterator c(f.getCombinations().begin());
665 Trig::FeatureContainer::combination_const_iterator cEnd(f.getCombinations().end());
666
667
668 if ( c==cEnd ) {
669 if(m_provider->msg().level() <= MSG::VERBOSE){
670 m_provider->msg(MSG::VERBOSE) << "No combinations: skipping this chain " << chainname << endmsg;
671 }
672 continue;
673 }
674
675 if(m_provider->msg().level() <= MSG::VERBOSE) {
676 m_provider->msg(MSG::VERBOSE) << "combinations for chain " << chainname << " " << (cEnd-c) << endmsg;
677 }
678
679 unsigned icomb = 0;
680
681 for( ; c!=cEnd ; ++c ) {
682
683 icomb++;
684
685 // now add rois to this ntuple chain
686
687 // Get seeding RoI
688 // std::vector< Trig::Feature<TrigRoiDescriptor> > initRois = c->get<TrigRoiDescriptor>("initialRoI", TrigDefs::alsoDeactivateTEs);
689 // std::vector< Trig::Feature<TrigRoiDescriptor> > initRois = c->get<TrigRoiDescriptor>("forID", TrigDefs::alsoDeactivateTEs);
690
691 std::vector< Trig::Feature<TrigRoiDescriptor> > initRois;
692
693 std::string roi_key = m_chainNames[ichain].roi();
694
695 if ( roi_key=="SuperRoi" && icomb>1 ) continue;
696
697 if ( roi_key!="" ) {
698 initRois = c->get<TrigRoiDescriptor>(roi_key, decisiontype );
699 }
700 else {
701 initRois = c->get<TrigRoiDescriptor>("forID", decisiontype );
702 if ( initRois.empty() ) initRois = c->get<TrigRoiDescriptor>("", decisiontype );
703 if ( initRois.empty() ) initRois = c->get<TrigRoiDescriptor>("initialRoI", decisiontype );
704 }
705
706 // std::cout << "initRois.size() " << initRois.size() << std::endl;
707
708 if ( initRois.empty() ) continue;
709
710
711 // for ( unsigned ir=0 ; ir<initRois.size() ; ir++ ) {
712 // std::cout << "\t" << ir << "\t" << *initRois[ir].cptr() << std::endl;
713 // }
714
715 // Skip chains seeded by multiple RoIs: not yet implemented
716 if(initRois.size()>1 && roi_key!="SuperRoi" ) {
717 if(m_provider->msg().level() <= MSG::VERBOSE)
718 m_provider->msg(MSG::VERBOSE) << " More than one initial RoI found for seeded chain " << chainname << ": not yet supported" << endmsg;
719 continue;
720 }
721
722 TIDARoiDescriptor* roiInfo = 0;
723
724 if( !initRois.empty() ) {
725 const TrigRoiDescriptor* roid = initRois[0].cptr();
726
727 if(m_provider->msg().level() <= MSG::VERBOSE)
728 m_provider->msg(MSG::VERBOSE) << " RoI descriptor for seeded chain " << chainname << " " << *roid << endmsg;
729
730 roiInfo = new TIDARoiDescriptor(TIDARoiDescriptorBuilder(*roid));
731 // roiInfo->etaHalfWidth(m_roiInfo->etaHalfWidth());
732 // roiInfo->phiHalfWidth(m_roiInfo->phiHalfWidth());
733 // roiInfo->etaHalfWidth(roid->etaHalfWidth());
734 // roiInfo->phiHalfWidth(roid->phiHalfWidth());
735 // roiInfo->zedHalfWidth(roid->zedHalfWidth());
736
737 if(m_provider->msg().level() <= MSG::VERBOSE)
738 m_provider->msg(MSG::VERBOSE) << "using chain roi " << *roid << endmsg;
739
740 }
741 else {
742 roiInfo = new TIDARoiDescriptor();
743 // roiInfo->etaHalfWidth(5);
744 // roiInfo->phiHalfWidth(M_PI);
745 // roiInfo->zedHalfWidth(m_roiInfo->zedHalfWidth());
746
747 if(m_provider->msg().level() <= MSG::WARNING)
748 m_provider->msg(MSG::WARNING) << "roi not found" << endmsg;
749
750 }
751
752 if(m_provider->msg().level() <= MSG::VERBOSE) m_provider->msg(MSG::VERBOSE) << *roiInfo << endmsg;
753
754 m_selectorTest->clear();
755
756 m_provider->msg(MSG::VERBOSE) << "Searching for collection " << key << endmsg;
757 // std::cout << "Searching for collection " << key << std::endl;
758
760 if ( key.find("InDetTrigParticleCreation")!=std::string::npos ||
761 key.find("_IDTrig")!=std::string::npos ||
762 key.find("_EFID")!=std::string::npos ||
763 chainName.find("EF_")!=std::string::npos ||
764 chainName.find("HLT_")!=std::string::npos ) {
765# ifdef XAODTRACKING_TRACKPARTICLE_H
766 if ( this->template selectTracks<xAOD::TrackParticleContainer>( m_selectorTest, c, key ) ) testbeamspot = this->template getBeamspot<xAOD::TrackParticleContainer>( c, key );
767 else if ( this->template selectTracks<Rec::TrackParticleContainer>( m_selectorTest, c, key ) );
768# else
769 if ( this->template selectTracks<Rec::TrackParticleContainer>( m_selectorTest, c, key ) );
770# endif
771 else if ( this->template selectTracks<TrackCollection>( m_selectorTest, c, key ) );
772 else if ( this->template selectTracks<TrigInDetTrackCollection>( m_selectorTest, c, truthMap, key, key_index ) );
773 else {
774 //m_provider->msg(MSG::WARNING) << "No track collection " << key << " found" << endmsg;
775 }
776 }
777 else {
779 if ( chainName.find("L2_")!=std::string::npos ) {
780 if ( this->template selectTracks<TrigInDetTrackCollection>( m_selectorTest, c, truthMap, key, key_index ) );
781 else if ( this->template selectTracks<Rec::TrackParticleContainer>( m_selectorTest, c, key ) );
782 else if ( this->template selectTracks<TrackCollection>( m_selectorTest, c, key ) );
783# ifdef XAODTRACKING_TRACKPARTICLE_H
784 else if ( this->template selectTracks<xAOD::TrackParticleContainer>( m_selectorTest, c, key ) ) testbeamspot = this->template getBeamspot<xAOD::TrackParticleContainer>( c, key );
785# endif
786 else m_provider->msg(MSG::WARNING) << "No track collection " << key << " found" << endmsg;
787 }
788 }
789
790
791 const std::vector<TIDA::Track*>& testtracks = m_selectorTest->tracks();
792
793 m_provider->msg(MSG::VERBOSE) << "test tracks.size() " << testtracks.size() << endmsg;
794 // std::cout << "test tracks.size() " << testtracks.size() << std::endl;
795
796 // std::cout << "\ttest tracks.size() " << testtracks.size() << std::endl;
797
798 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
799 m_provider->msg(MSG::VERBOSE) << "test tracks.size() " << testtracks.size() << endmsg;
800 for ( int ii=testtracks.size() ; ii-- ; ) {
801 m_provider->msg(MSG::VERBOSE) << " test track " << ii << " " << *testtracks[ii] << endmsg;
802 //test_tracks.push_back(testtracks.at(ii));
803 }
804 }
805
806 chain.addRoi( *roiInfo );
807
808 chain.back().addTracks(testtracks);
809
810 delete roiInfo;
811
813
815
816 // std::cout << "vertex " << vtx_name << "\tchain " << chainName << "\tconfig " << chainConfig << std::endl;
817
818 if ( vtx_name!="" ) {
819
820 m_provider->msg(MSG::VERBOSE) << "\tFetch xAOD::VertexContainer for chain " << chainConfig << " with key " << vtx_name << endmsg;
821
822 std::vector< Trig::Feature<xAOD::VertexContainer> > xaodtrigvertices = c->get<xAOD::VertexContainer>(vtx_name);
823
824 if ( xaodtrigvertices.empty() ) {
825 if ( m_provider->msg().level() <= MSG::DEBUG ) {
826 m_provider->msg(MSG::WARNING) << "\tNo xAOD::VertexContainer for chain " << chainConfig << " for key " << vtx_name << endmsg;
827 }
828 }
829 else {
830
831 m_provider->msg(MSG::VERBOSE) << "\txAOD::VertexContainer found with size " << xaodtrigvertices.size() << "\t" << vtx_name << endmsg;
832
833 for ( unsigned iv=0 ; iv<xaodtrigvertices.size() ; iv++ ) {
834
835 const xAOD::VertexContainer* vert = xaodtrigvertices[iv].cptr();
836
837 m_provider->msg(MSG::VERBOSE) << "\t" << iv << " xAOD VxContainer for " << chainConfig << " " << vert << " key " << vtx_name << endmsg;
838
839 xAOD::VertexContainer::const_iterator vtxitr = vert->begin();
840
841 for ( ; vtxitr != vert->end(); ++vtxitr) {
844 // if ( ( (*vtxitr)->nTrackParticles()>0 && (*vtxitr)->vertexType()!=0 ) || vtx_name=="EFHistoPrmVtx" ) {
845 if ( (*vtxitr)->vertexType()!=0 || vtx_name=="EFHistoPrmVtx" ) {
846 chain.back().addVertex( TIDA::Vertex( (*vtxitr)->x(),
847 (*vtxitr)->y(),
848 (*vtxitr)->z(),
850 (*vtxitr)->covariancePosition()(Trk::x,Trk::x),
851 (*vtxitr)->covariancePosition()(Trk::y,Trk::y),
852 (*vtxitr)->covariancePosition()(Trk::z,Trk::z),
853 (*vtxitr)->nTrackParticles(),
855 (*vtxitr)->chiSquared(),
856 (*vtxitr)->numberDoF() ) );
857
858 }
859 }
860 }
861
862 }
863
864 }
865
866 }
867
868 }
869
870
871 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
872 m_provider->msg(MSG::VERBOSE) << "event: " << *m_event << endmsg;
873 }
874
876
877 for ( unsigned iroi=0 ; iroi<chain.size() ; iroi++ ) {
878
879 m_selectorRef->clear();
880
881 if ( this->filterOnRoi() ) {
882 filterRef.setRoi( &chain.rois().at(iroi).roi() );
883 filterRef.containtracks( m_containTracks );
884 }
885 else filterRef.setRoi( 0 );
886
887 test_tracks.clear();
888
889
903
904 if ( m_provider->msg().level() <= MSG::VERBOSE )
905 m_provider->msg(MSG::VERBOSE) << "MC Truth flag " << m_mcTruth << endmsg;
906
907 bool foundTruth = false;
908
909 if ( !m_doOffline && m_mcTruth ) {
910
911 if ( this->filterOnRoi() ) filter_truth.setRoi( &chain.rois().at(iroi).roi() );
912 else filter_truth.setRoi( 0 ); // don't filter on RoI unless needed
913
914 selectorTruth.clear();
915
916 if ( m_provider->msg().level() <= MSG::VERBOSE )
917 m_provider->msg(MSG::VERBOSE) << "getting Truth" << endmsg;
918
919 if ( m_provider->evtStore()->template contains<TruthParticleContainer>("INav4MomTruthEvent") ) {
920 //ESD
921 this->template selectTracks<TruthParticleContainer>( &selectorTruth, "INav4MomTruthEvent" );
922 foundTruth = true;
923 }
924 else if ( m_provider->evtStore()->template contains<TruthParticleContainer>("SpclMC") ) {
926 this->template selectTracks<TruthParticleContainer>( &selectorTruth, "SpclMC");
927 foundTruth = true;
928 }
929 else if ( m_provider->evtStore()->template contains<TruthParticleContainer>("") ) {
931 this->template selectTracks<TruthParticleContainer>( &selectorTruth, "");
932 foundTruth = true;
933 }
934 else
935 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
936 m_provider->msg(MSG::VERBOSE) << "Truth not found - none whatsoever!" << endmsg;
937 }
938 }
939
940
941 if ( !m_doOffline && m_mcTruth && !foundTruth ) {
942
943 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
944 m_provider->msg(MSG::VERBOSE) << "getting Truth" << endmsg;
945 }
946
948
949 const DataHandle<McEventCollection> mcevent;
950
952
953 std::string keys[4] = { "GEN_AOD", "TruthEvent", "", "G4Truth" };
954
955 std::string key = "";
956
957 bool foundcollection = false;
958
959 for ( int ik=0 ; ik<4 ; ik++ ) {
960
961 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
962 m_provider->msg(MSG::VERBOSE) << "Try McEventCollection: " << keys[ik] << endmsg;
963 }
964
965 if ( !m_provider->evtStore()->template contains<McEventCollection>(keys[ik]) ) {
966 if( m_provider->msg().level() <= MSG::VERBOSE )
967 m_provider->msg(MSG::VERBOSE) << "No McEventCollection: " << keys[ik] << endmsg;
968 continue;
969 }
970
971 if ( m_provider->msg().level() <= MSG::VERBOSE )
972 m_provider->msg(MSG::VERBOSE) << "evtStore()->retrieve( mcevent, " << keys[ik] << " )" << endmsg;
973
974 if ( this->template retrieve( mcevent, keys[ik] ).isFailure() ) {
975 if ( m_provider->msg().level() <= MSG::VERBOSE )
976 m_provider->msg(MSG::VERBOSE) << "Failed to get McEventCollection: " << keys[ik] << endmsg;
977 }
978 else {
980 key = keys[ik];
981 if(m_provider->msg().level() <= MSG::VERBOSE)
982 m_provider->msg(MSG::VERBOSE) << "Found McEventCollection: " << key << endmsg;
983 foundcollection = true;
984 break;
985 }
986 }
987
989 if ( !foundcollection ) {
990 if(m_provider->msg().level() <= MSG::VERBOSE)
991 m_provider->msg(MSG::WARNING) << "No MC Truth Collections of any sort, whatsoever!!!" << endmsg;
992
993 // m_tree->Fill();
994 // return StatusCode::FAILURE;
995
996 return;
997 }
998
999 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
1000 m_provider->msg(MSG::VERBOSE) << "Found McEventCollection: " << key << "\tNevents " << mcevent->size() << endmsg;
1001 }
1002
1003 McEventCollection::const_iterator evitr = mcevent->begin();
1004 McEventCollection::const_iterator evend = mcevent->end();
1005
1006 unsigned ie = 0;
1007 unsigned ip = 0;
1008
1009 unsigned ie_ip = 0;
1010
1011 while ( evitr!=evend ) {
1012
1013 int _ip = 0;
1014
1015 int pid = HepMC::signal_process_id((*evitr));
1016
1017 //The logic should be clarified here
1018 if ( pid!=0 ) {
1019
1020 for (auto pitr: *(*evitr)) {
1021
1022 selectorTruth.selectTrack( pitr );
1023
1024 ++_ip;
1025
1026 }
1027
1028 }
1029 ++ie;
1030 ++evitr;
1031
1032 if ( _ip>0 ) {
1034 // m_provider->msg(MSG::VERBOSE) << "Found " << ie << "\tpid " << pid << "\t with " << ip << " TruthParticles (GenParticles)" << endmsg;
1035 ++ie_ip;
1036 ip += _ip;
1037 }
1038 }
1039
1040 if(m_provider->msg().level() <= MSG::VERBOSE){
1041 m_provider->msg(MSG::VERBOSE) << "Found " << ip << " TruthParticles (GenParticles) in " << ie_ip << " GenEvents out of " << ie << endmsg;
1042 m_provider->msg(MSG::VERBOSE) << "selected " << selectorTruth.size() << " TruthParticles (GenParticles)" << endmsg;
1043 }
1044
1045 if(selectorTruth.size() > 0) foundTruth = true;
1046
1047 if ( !(ip>0) ) {
1048 if (m_provider->msg().level() <= MSG::VERBOSE) m_provider->msg(MSG::WARNING) << "NO TRUTH PARTICLES - returning" << endmsg;
1049 return;
1050 }
1051
1052 }
1053
1055
1056 // m_provider->msg(MSG::VERBOSE) << " Offline tracks " << endmsg;
1057
1058 if ( m_doOffline ) {
1059
1060# ifdef XAODTRACKING_TRACKPARTICLE_H
1061 if ( m_provider->evtStore()->template contains<xAOD::TrackParticleContainer>("InDetTrackParticles") ) {
1062 this->template selectTracks<xAOD::TrackParticleContainer>( m_selectorRef, "InDetTrackParticles" );
1063 refbeamspot = this->template getBeamspot<xAOD::TrackParticleContainer>( "InDetTrackParticles" );
1064 }
1065 else if (m_provider->evtStore()->template contains<Rec::TrackParticleContainer>("TrackParticleCandidate") ) {
1066 this->template selectTracks<Rec::TrackParticleContainer>( m_selectorRef, "TrackParticleCandidate" );
1067 }
1068# else
1069 if (m_provider->evtStore()->template contains<Rec::TrackParticleContainer>("TrackParticleCandidate") ) {
1070 this->template selectTracks<Rec::TrackParticleContainer>( m_selectorRef, "TrackParticleCandidate" );
1071 }
1072# endif
1073 else if ( m_provider->msg().level() <= MSG::WARNING ) {
1074 m_provider->msg(MSG::WARNING) << " Offline tracks not found " << endmsg;
1075 }
1076
1077 ref_tracks = m_selectorRef->tracks();
1078
1079 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
1080 m_provider->msg(MSG::VERBOSE) << "ref tracks.size() " << m_selectorRef->tracks().size() << endmsg;
1081 for ( int ii=m_selectorRef->tracks().size() ; ii-- ; ) {
1082 m_provider->msg(MSG::VERBOSE) << " ref track " << ii << " " << *m_selectorRef->tracks()[ii] << endmsg;
1083 }
1084 }
1085
1086 }
1087 else {
1089 if ( m_mcTruth && foundTruth ){
1090 ref_tracks=selectorTruth.tracks();
1091 }
1092 }
1093
1094
1095
1096 test_tracks.clear();
1097
1098 for ( unsigned itrk=0 ; itrk<chain.rois().at(iroi).tracks().size() ; itrk++ ) {
1099 test_tracks.push_back(&(chain.rois().at(iroi).tracks().at(itrk)));
1100 }
1101
1102
1103 // std::cout << "sutt track multiplicities: offline " << offline_tracks.size() << "\ttest " << test_tracks.size() << std::endl;
1104
1105 _analysis->setvertices( vertices.size() );
1106
1107 if ( refbeamspot.size()>0 ) _analysis->setBeamRef( refbeamspot );
1108 if ( testbeamspot.size()>0 ) _analysis->setBeamTest( testbeamspot );
1109
1112
1113 if ( m_runPurity ) {
1114
1115 if ( this->getUseHighestPT() ) HighestPTOnly( test_tracks );
1116
1117 if ( m_pTthreshold>0 ) FilterPT( test_tracks, m_pTthreshold );
1118
1120 m_NRois++;
1121 m_NRefTracks += test_tracks.size();
1122 m_NTestTracks += ref_tracks.size();
1123
1125 m_associator->match( test_tracks, ref_tracks );
1126
1127 _analysis->execute( test_tracks, ref_tracks, m_associator );
1128
1129 }
1130 else {
1131
1133
1134 if ( this->getUseHighestPT() ) HighestPTOnly( ref_tracks );
1135
1137
1138 if ( m_pTthreshold>0 ) FilterPT( ref_tracks, m_pTthreshold );
1139
1141 m_NRois++;
1142 m_NRefTracks += ref_tracks.size();
1143 m_NTestTracks += test_tracks.size();
1144
1146 m_associator->match( ref_tracks, test_tracks );
1147
1148 // std::cout << "SUTT: execute : N tracks " << ref_tracks.size() << " " << test_tracks.size() << std::endl;
1149
1150 _analysis->setroi( &chain.rois().at(iroi).roi() );
1151 _analysis->execute( ref_tracks, test_tracks, m_associator );
1152
1153 // std::cout << "chain " << m_chainNames[ichain] << " " << "\tvtx name " << vtx_name << std::endl;
1154
1155 if ( vtx_name!="" ) {
1157 std::vector<TIDA::Vertex> vr = chain.rois().at(iroi).vertices();
1158 std::vector<TIDA::Vertex*> vtx_rec;
1159 for ( unsigned iv=0 ; iv<vr.size() ; iv++ ) vtx_rec.push_back( &vr[iv] );
1160
1161 std::vector<TIDA::Vertex*> vtx;
1162 if ( this->getVtxIndex()<0 ) {
1163 for ( unsigned iv=0 ; iv<vertices.size() ; iv++ ) vtx.push_back( &vertices[iv] );
1164 }
1165 else {
1166 if ( vertices.size()>unsigned(this->getVtxIndex()) ) vtx.push_back( &vertices[this->getVtxIndex()] );
1167 }
1168
1169 _analysis->execute_vtx( vtx, vtx_rec, m_event );
1170 }
1171
1172 }
1173
1174 if ( _analysis->debug() ) {
1175 m_provider->msg(MSG::INFO) << "Missing track for " << m_chainNames[ichain]
1176 << "\trun " << run_number
1177 << "\tevent " << event_number
1178 << "\tlb " << lumi_block << endmsg;
1179 }
1180
1181 }
1182
1183 }
1184
1185 if ( m_provider->msg().level() <= MSG::VERBOSE ) {
1186 m_provider->msg(MSG::VERBOSE) << "\n\nEvent " << *m_event << endmsg;
1187 }
1188 }
1189
1190
1191
1192 virtual void book() {
1193
1194 if(m_provider->msg().level() <= MSG::VERBOSE)
1195 m_provider->msg(MSG::VERBOSE) << "AnalysisConfig_Tier0::book() " << name() << endmsg;
1196
1197 m_provider->msg(MSG::ERROR) << "AnalysisConfig_Tier0::book() should no longer ever be called: " << name() << endmsg;
1198
1199 }
1200
1201
1202
1203 virtual void finalize() {
1204
1205 if(m_provider->msg().level() <= MSG::VERBOSE)
1206 m_provider->msg(MSG::VERBOSE) << "AnalysisConfig_Tier0::finalise() " << m_provider->name() << endmsg;
1207
1208 m_analysis->finalise();
1209
1210 m_provider->msg(MSG::INFO) << m_provider->name() << " " << m_chainNames[0] << " \tNRois processed: " << m_NRois << "\tRef tracks: " << m_NRefTracks << "\tTestTracks: " << m_NTestTracks << endmsg;
1211
1212 if(m_provider->msg().level() <= MSG::VERBOSE)
1213 m_provider->msg(MSG::VERBOSE) << m_provider->name() << " finalised" << endmsg;
1214
1215 }
1216
1217
1218protected:
1219
1220 bool m_useBeamCondSvc;
1221
1222 TIDA::Event* m_event;
1223
1224 std::vector<ChainString> m_chainNames;
1225 std::vector<Analysis_Tier0*> m_analyses;
1226 std::string m_testType;
1227
1228 bool m_doOffline;
1229 bool m_doMuons;
1230 bool m_doElectrons;
1231 bool m_doTaus;
1232 bool m_doBjets;
1233 bool m_hasTruthMap;
1234 bool m_doTauThreeProng;
1235 bool m_tauEtCutOffline;
1236
1237 std::string m_outputFileName;
1238
1240 int m_NRois;
1241 int m_NRefTracks;
1242 int m_NTestTracks;
1243
1244 bool m_runPurity;
1245
1246 bool m_shifter;
1247
1248 double m_pTthreshold;
1249
1250 bool m_first;
1251
1252 bool m_containTracks;
1253
1254};
1255
1256
1257
1258#endif // TrigInDetAnalysisExample_T_AnalysisConfig_Tier0_H
1259
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:63
Analysis_Tier0.h
Associator_BestMatch.h
Erkcan's track matchers.
EventID.h
This class provides a unique identification for each event, in terms of run/event number and/or a tim...
Filter_AcceptAll.h
Filter_RoiSelector.h
Filter_etaPT.h
generic track filter to filter on eta and PT only
passed
bool passed(DecisionID id, const DecisionIDContainer &)
checks if required decision ID is in the set of IDs in the container
Definition TrigCompositeUtilsRoot.cxx:118
TIDAEvent.h
Basic event class to contain a vector of chains for trigger analysis.
TIDATools.h
useful tool for the TrigInDetAnalysis class code
HighestPTOnly
void HighestPTOnly(std::vector< T * > &tracks)
Definition TIDATools.h:20
FilterPT
void FilterPT(std::vector< T * > &tracks, double pt)
Definition TIDATools.h:41
T_AnalysisConfig.h
TrackAssociator
TIDA::Associator< TIDA::Track > TrackAssociator
Definition TrackAssociator.h:26
Analysis_Tier0
Definition Analysis_Tier0.h:30
Analysis_Tier0::execute_vtx
virtual void execute_vtx(const std::vector< TIDA::Vertex * > &vtx0, const std::vector< TIDA::Vertex * > &vtx1, const TIDA::Event *tevt=0)
Definition Analysis_Tier0.cxx:755
ChainString::extra
const std::string & extra() const
Definition ChainString.h:38
ChainString::roi
const std::string & roi() const
Definition ChainString.h:35
ChainString::subs
std::string subs(std::string s) const
Definition ChainString.cxx:127
ChainString::raw
const std::string & raw() const
Definition ChainString.h:42
ChainString::head
const std::string & head() const
Definition ChainString.h:33
ChainString::tail
const std::string & tail() const
Definition ChainString.h:34
ChainString::vtx
const std::string & vtx() const
Definition ChainString.h:36
ChainString::element
const std::string & element() const
Definition ChainString.h:37
DataHandle
an iterator over instances of a given type in StoreGateSvc.
Definition DataHandle.h:43
DataVector< Trk::VxCandidate >::const_iterator
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
DataVector::end
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
DataVector::begin
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
EventInfo::event_ID
EventID * event_ID()
the unique identification of the event.
Definition EventInfo/EventInfo/EventInfo.h:224
EventInfo::averageInteractionsPerCrossing
float averageInteractionsPerCrossing() const
average interactions per crossing for all BCIDs - for out-of-time pile-up
Definition EventInfo/EventInfo/EventInfo.h:296
Filter_Combined
Definition Filters.h:225
Filter_Combined::setRoi
void setRoi(TIDARoiDescriptor *r)
Definition Filters.h:236
Filter_Combined::containtracks
void containtracks(bool b=true)
set / unset the flag to determine whether tracks should be fully contained in the RoI or not
Definition Filters.h:242
Filter_True
default simple filter which accepts all tracks
Definition Filters.h:26
Filter_etaPT
Definition Filter_etaPT.h:27
TIDARoiDescriptor
Describes the Region of Ineterest geometry It has basically 8 parameters.
Definition TIDARoiDescriptor.h:40
T_AnalysisConfig_Tier0::m_chainNames
std::vector< ChainString > m_chainNames
Definition T_AnalysisConfig_Tier0.h:1224
T_AnalysisConfig_Tier0::setRunPurity
void setRunPurity(bool b)
Definition T_AnalysisConfig_Tier0.h:153
T_AnalysisConfig_Tier0::m_first
bool m_first
Definition T_AnalysisConfig_Tier0.h:1250
T_AnalysisConfig_Tier0::m_doOffline
bool m_doOffline
Definition T_AnalysisConfig_Tier0.h:1228
T_AnalysisConfig_Tier0::m_NTestTracks
int m_NTestTracks
Definition T_AnalysisConfig_Tier0.h:1242
T_AnalysisConfig_Tier0::m_doTauThreeProng
bool m_doTauThreeProng
Definition T_AnalysisConfig_Tier0.h:1234
T_AnalysisConfig_Tier0::containTracks
void containTracks(bool b)
Definition T_AnalysisConfig_Tier0.h:159
T_AnalysisConfig_Tier0::m_doBjets
bool m_doBjets
Definition T_AnalysisConfig_Tier0.h:1232
T_AnalysisConfig_Tier0::T_AnalysisConfig_Tier0
T_AnalysisConfig_Tier0(const std::string &analysisInstanceName, const std::string &testChainName, const std::string &testType, const std::string &testKey, const std::string &referenceChainName, const std::string &referenceType, const std::string &referenceKey, TrackFilter *testFilter, TrackFilter *referenceFilter, TrackAssociator *associator, TrackAnalysis *analysis)
Definition T_AnalysisConfig_Tier0.h:97
T_AnalysisConfig_Tier0::_analysis
Analysis_Tier0 * _analysis
Definition T_AnalysisConfig_Tier0.h:163
T_AnalysisConfig_Tier0::m_doElectrons
bool m_doElectrons
Definition T_AnalysisConfig_Tier0.h:1230
T_AnalysisConfig_Tier0::m_pTthreshold
double m_pTthreshold
Definition T_AnalysisConfig_Tier0.h:1248
T_AnalysisConfig_Tier0::m_testType
std::string m_testType
Definition T_AnalysisConfig_Tier0.h:1226
T_AnalysisConfig_Tier0::loop
virtual void loop()
Definition T_AnalysisConfig_Tier0.h:188
T_AnalysisConfig_Tier0::m_NRois
int m_NRois
output stats
Definition T_AnalysisConfig_Tier0.h:1240
T_AnalysisConfig_Tier0::finalize
virtual void finalize()
Definition T_AnalysisConfig_Tier0.h:1203
T_AnalysisConfig_Tier0::useBeamCondSvc
void useBeamCondSvc(bool b)
Definition T_AnalysisConfig_Tier0.h:157
T_AnalysisConfig_Tier0::m_useBeamCondSvc
bool m_useBeamCondSvc
Definition T_AnalysisConfig_Tier0.h:1220
T_AnalysisConfig_Tier0::m_analyses
std::vector< Analysis_Tier0 * > m_analyses
Definition T_AnalysisConfig_Tier0.h:1225
T_AnalysisConfig_Tier0::m_doTaus
bool m_doTaus
Definition T_AnalysisConfig_Tier0.h:1231
T_AnalysisConfig_Tier0::m_hasTruthMap
bool m_hasTruthMap
Definition T_AnalysisConfig_Tier0.h:1233
T_AnalysisConfig_Tier0::m_tauEtCutOffline
bool m_tauEtCutOffline
Definition T_AnalysisConfig_Tier0.h:1235
T_AnalysisConfig_Tier0::m_shifter
bool m_shifter
Definition T_AnalysisConfig_Tier0.h:1246
T_AnalysisConfig_Tier0::m_outputFileName
std::string m_outputFileName
Definition T_AnalysisConfig_Tier0.h:1237
T_AnalysisConfig_Tier0::setShifter
void setShifter(bool b)
Definition T_AnalysisConfig_Tier0.h:155
T_AnalysisConfig_Tier0::m_NRefTracks
int m_NRefTracks
Definition T_AnalysisConfig_Tier0.h:1241
T_AnalysisConfig_Tier0::m_doMuons
bool m_doMuons
Definition T_AnalysisConfig_Tier0.h:1229
T_AnalysisConfig_Tier0::book
virtual void book()
Definition T_AnalysisConfig_Tier0.h:1192
T_AnalysisConfig_Tier0::~T_AnalysisConfig_Tier0
virtual ~T_AnalysisConfig_Tier0()
Definition T_AnalysisConfig_Tier0.h:151
T_AnalysisConfig_Tier0::m_runPurity
bool m_runPurity
Definition T_AnalysisConfig_Tier0.h:1244
T_AnalysisConfig_Tier0::m_containTracks
bool m_containTracks
Definition T_AnalysisConfig_Tier0.h:1252
T_AnalysisConfig_Tier0::m_event
TIDA::Event * m_event
Definition T_AnalysisConfig_Tier0.h:1222
T_AnalysisConfig< T >::analysis
const TrackAnalysis * analysis() const
Definition T_AnalysisConfig.h:258
T_AnalysisConfig< T >::getVtxIndex
int getVtxIndex() const
Definition T_AnalysisConfig.h:282
T_AnalysisConfig< T >::m_analysis
TrackAnalysis * m_analysis
Definition T_AnalysisConfig.h:986
T_AnalysisConfig< T >::requireDecision
bool requireDecision() const
Definition T_AnalysisConfig.h:288
T_AnalysisConfig< T >::name
const std::string & name() const
Definition T_AnalysisConfig.h:235
T_AnalysisConfig< T >::m_mcTruth
bool m_mcTruth
Definition T_AnalysisConfig.h:988
T_AnalysisConfig< T >::filterOnRoi
bool filterOnRoi() const
Definition T_AnalysisConfig.h:284
T_AnalysisConfig< T >::getUseHighestPT
bool getUseHighestPT() const
Definition T_AnalysisConfig.h:279
T_AnalysisConfig< T >::m_selectorTest
TrigTrackSelector * m_selectorTest
Definition T_AnalysisConfig.h:982
T_AnalysisConfig< T >::m_provider
T * m_provider
Definition T_AnalysisConfig.h:957
T_AnalysisConfig::m_analysisInstanceName
std::string m_analysisInstanceName
Definition T_AnalysisConfig.h:965
T_AnalysisConfig< T >::selectTracks
bool selectTracks(TrigTrackSelector *selector, const ElementLink< TrigRoiDescriptorCollection > &roi_link, const std::string &key="")
Definition T_AnalysisConfig.h:441
T_AnalysisConfig< T >::m_filters
std::vector< std::vector< TrackFilter * > > m_filters
Definition T_AnalysisConfig.h:979
T_AnalysisConfig< T >::m_selectorRef
TrigTrackSelector * m_selectorRef
Definition T_AnalysisConfig.h:981
T_AnalysisConfig< T >::retrieve
StatusCode retrieve(Collection const *&collection, const std::string &key="")
Definition T_AnalysisConfig.h:573
T_AnalysisConfig< T >::getBeamspot
std::vector< double > getBeamspot(const std::string &key)
Definition T_AnalysisConfig.h:634
T_AnalysisConfig< T >::m_tdt
ToolHandle< Trig::TrigDecisionTool > * m_tdt
Definition T_AnalysisConfig.h:959
T_AnalysisConfig< T >::T_AnalysisConfig
T_AnalysisConfig(const std::string &analysisInstanceName, const std::string &testChainName, const std::string &testType, const std::string &testKey, const std::string &referenceChainName, const std::string &referenceType, const std::string &referenceKey, const std::string &selectionChainName, const std::string &selectionType, const std::string &selectionKey, TrackFilter *testFilter, TrackFilter *referenceFilter, TrackFilter *selectionFilter, TrackAssociator *associator, TrackAnalysis *analysis)
Definition T_AnalysisConfig.h:69
T_AnalysisConfig< T >::m_associator
TrackAssociator * m_associator
Definition T_AnalysisConfig.h:985
T_AnalysisConfig< T >::m_keepAllEvents
bool m_keepAllEvents
Definition T_AnalysisConfig.h:995
TrigRoiDescriptor
nope - should be used for standalone also, perhaps need to protect the class def bits ifndef XAOD_ANA...
Definition TrigRoiDescriptor.h:56
TrigTrackSelector::selectTrack
bool selectTrack(const TrigInDetTrack *track, const TrigInDetTrackTruthMap *truthMap=0)
neater code to make use of vector function also for a single ancestor pdgid, instead of the full code...
Definition TrigTrackSelector.cxx:116
TrigTrackSelector::clear
virtual void clear() override
Definition TrigTrackSelector.h:102
Trig::FeatureContainer::combination_const_iterator
std::vector< Combination >::const_iterator combination_const_iterator
Definition FeatureContainer.h:45
contains
bool contains(const std::string &s, const std::string &regx)
does a string contain the substring
Definition hcg.cxx:114
HepMC::signal_process_id
int signal_process_id(const GenEvent &e)
Definition GenEvent.h:636
Trk::x
@ x
Definition ParamDefs.h:55
Trk::z
@ z
global position (cartesian)
Definition ParamDefs.h:57
Trk::y
@ y
Definition ParamDefs.h:56
xAOD::EventInfo
EventInfo_v1 EventInfo
Definition of the latest event info version.
Definition IEventInfoCnvTool.h:18
xAOD::VertexContainer
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".
Definition VertexContainer.h:14
Generated on for ATLAS Offline Software by doxygen 1.14.0