ATLAS Offline Software
Loading...
Searching...
No Matches
TrigEDMChecker.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
3*/
4
6
7#include "TrigEDMChecker.h"
8
9#include "AthViews/ViewHelper.h"
10#include "AthViews/View.h"
11
12#include "TrigConfHLTUtils/HLTUtils.h"
13#include "TrigNavStructure/TriggerElement.h"
14#include "TrigT1Interfaces/RecEmTauRoI.h"
15
16#include "xAODEgamma/Electron.h"
17#include "xAODEgamma/ElectronContainer.h"
18#include "xAODEgamma/Photon.h"
19#include "xAODEgamma/PhotonContainer.h"
20
21#include "xAODJet/JetConstituentVector.h"
22#include "xAODJet/JetContainer.h"
23
24#include "xAODMuon/MuonContainer.h"
25
26#include "xAODTau/TauDefs.h"
27#include "xAODTau/TauJet.h"
28#include "xAODTau/TauJetAuxContainer.h"
29#include "xAODTau/TauJetContainer.h"
30
31#include "xAODTracking/TrackParticleContainer.h"
32#include "xAODTracking/VertexContainer.h"
33
34#include "xAODTrigBphys/TrigBphys.h"
35#include "xAODTrigBphys/TrigBphysContainer.h"
36
37#include "xAODTrigCalo/TrigEMCluster.h"
38#include "xAODTrigCalo/TrigEMClusterContainer.h"
39
40#include "xAODTrigEgamma/TrigElectron.h"
41#include "xAODTrigEgamma/TrigElectronContainer.h"
42#include "xAODTrigEgamma/TrigPhoton.h"
43#include "xAODTrigEgamma/TrigPhotonContainer.h"
44
45#include "xAODTrigMinBias/TrigSpacePointCounts.h"
46#include "xAODTrigMinBias/TrigSpacePointCountsContainer.h"
47#include "xAODTrigMinBias/TrigT2MbtsBits.h"
48#include "xAODTrigMinBias/TrigT2MbtsBitsContainer.h"
49#include "xAODTrigMinBias/TrigTrackCounts.h"
50#include "xAODTrigMinBias/TrigTrackCountsContainer.h"
51#include "xAODTrigMinBias/TrigVertexCounts.h"
52#include "xAODTrigMinBias/TrigVertexCountsContainer.h"
53
54#include "xAODTrigMissingET/TrigMissingETAuxContainer.h"
55#include "xAODTrigMissingET/TrigMissingETContainer.h"
56
57#include "xAODTrigger/TrigPassBits.h"
58#include "xAODTrigger/TrigPassBitsContainer.h"
59#include "xAODTrigger/TriggerMenuContainer.h"
60
61#include "TrigSteeringEvent/TrigRoiDescriptor.h"
62#include "TrigRoiConversion/RoiSerialise.h"
63#include "xAODTrigger/RoiDescriptorStore.h"
64
65#include <iostream>
66#include <fstream>
67#include <queue>
68
69
70StatusCode TrigEDMChecker::initialize() {
71
72 ATH_CHECK( m_navigationHandleKey.initialize() );
73 ATH_CHECK( m_decisionsKey.initialize() );
74 ATH_CHECK( m_navigationTool.retrieve() );
75
76 ATH_MSG_DEBUG("Initializing TrigEDMChecker");
77
78 ATH_MSG_INFO("REGTEST Initializing...");
79 ATH_MSG_INFO("REGTEST m_doDumpAll = " << m_doDumpAll );
80 ATH_MSG_INFO("REGTEST m_doDumpLVL1_ROI = " << m_doDumpLVL1_ROI);
81 ATH_MSG_INFO("REGTEST m_doDumpxAODTrigMissingET = " << m_doDumpxAODTrigMissingET );
82 ATH_MSG_INFO("REGTEST m_doDumpTrigL2BphysContainer = " << m_doDumpTrigL2BphysContainer );
83 ATH_MSG_INFO("REGTEST m_doDumpTrigEFBphysContainer = " << m_doDumpTrigEFBphysContainer );
84 ATH_MSG_INFO("REGTEST m_doDumpxAODJetContainer = " << m_doDumpxAODJetContainer );
85 ATH_MSG_INFO("REGTEST m_doDumpxAODMuonContainer = " << m_doDumpxAODMuonContainer );
86 ATH_MSG_INFO("REGTEST m_doDumpxAODTrigElectronContainer= " << m_doDumpxAODTrigElectronContainer );
87 ATH_MSG_INFO("REGTEST m_doDumpxAODTrigPhotonContainer = " << m_doDumpxAODTrigPhotonContainer );
88 ATH_MSG_INFO("REGTEST m_doDumpxAODElectronContainer = " << m_doDumpxAODElectronContainer );
89 ATH_MSG_INFO("REGTEST m_doDumpxAODPhotonContainer = " << m_doDumpxAODPhotonContainer );
90 ATH_MSG_INFO("REGTEST m_doDumpxAODTrackParticle = " << m_doDumpxAODTrackParticle );
91 ATH_MSG_INFO("REGTEST m_doDumpxAODVertex = " << m_doDumpxAODVertex );
92 ATH_MSG_INFO("REGTEST m_doDumpxAODTauJetContainer = " << m_doDumpxAODTauJetContainer );
93 ATH_MSG_INFO("REGTEST m_doDumpxAODTrigMinBias = " << m_doDumpxAODTrigMinBias );
94 ATH_MSG_INFO("REGTEST m_doDumpStoreGate = " << m_doDumpStoreGate );
95 ATH_MSG_INFO("REGTEST m_doDumpAllTrigComposite = " << m_doDumpAllTrigComposite );
96 ATH_MSG_INFO("REGTEST m_dumpTrigCompositeContainers = " << m_dumpTrigCompositeContainers );
97 ATH_MSG_INFO("REGTEST m_doDumpTrigCompsiteNavigation = " << m_doDumpTrigCompsiteNavigation );
98 ATH_MSG_INFO("REGTEST m_doTDTCheck = " << m_doTDTCheck );
99
100 if(m_doDumpxAODMuonContainer || m_doDumpAll) {
101 ATH_CHECK( m_muonPrinter.retrieve() );
102 }
103 else m_muonPrinter.disable(); // to avoid auto-retrieval
104
105 if (m_doDumpTrigCompsiteNavigation) {
106 ATH_CHECK( m_clidSvc.retrieve() );
107 }
108
109 if (m_doTDTCheck || m_doDumpTrigCompsiteNavigation) {
110 ATH_CHECK( m_trigDec.retrieve() );
111 ATH_MSG_INFO("TDT Executing with navigation format: " << m_trigDec->getNavigationFormat());
112 }
113 ATH_CHECK(m_muonTracksKey.initialize(m_doDumpAll || m_doTDTCheck));
114
115 return StatusCode::SUCCESS;
116}
117
118
119StatusCode TrigEDMChecker::execute(const EventContext& ctx) {
120
121 ATH_MSG_INFO( " ==========START of event===========" );
122
123 if(m_doDumpAll || m_doDumpLVL1_ROI ){
124 StatusCode sc = dumpLVL1_ROI();
125 if (sc.isFailure()) {
126 ATH_MSG_ERROR("The method dumpLVL1_ROI() failed");
127 }
128 }
129
130 if(m_doDumpAll || m_doDumpxAODTrigMissingET){
131 StatusCode sc = dumpxAODTrigMissingET();
132 if (sc.isFailure()) {
133 ATH_MSG_ERROR("The method dumpxAODTrigMissingET() failed");
134 }
135 }
136
137 if(m_doDumpAll || m_doDumpxAODTrigEMCluster){
138 StatusCode sc = dumpxAODTrigEMCluster();
139 if (sc.isFailure()) {
140 ATH_MSG_ERROR("The method dumpxAODTrigEMCluster() failed");
141 }
142 }
143
144 if(m_doDumpAll || m_doDumpxAODTrigEMClusterContainer){
145 StatusCode sc = dumpxAODTrigEMClusterContainer();
146 if (sc.isFailure()) {
147 ATH_MSG_ERROR("The method dumpxAODTrigEMClusterContainer() failed");
148 }
149 }
150
151 if(m_doDumpAll || m_doDumpxAODJetContainer){
152 StatusCode sc = dumpxAODJetContainer();
153 if (sc.isFailure()) {
154 ATH_MSG_ERROR("The method dumpxAODJetContainer() failed");
155 }
156 }
157
158 if(m_doDumpAll || m_doDumpTrigL2BphysContainer){
159 StatusCode sc = dumpTrigL2BphysContainer();
160 if (sc.isFailure()) {
161 ATH_MSG_ERROR("The method dumpTrigL2BphysContainer() failed");
162 }
163 }
164
165 if(m_doDumpAll || m_doDumpTrigEFBphysContainer){
166 StatusCode sc = dumpTrigEFBphysContainer();
167 if (sc.isFailure()) {
168 ATH_MSG_ERROR("The method dumpTrigEFBphysContainer() failed");
169 }
170 }
171
172 if(m_doDumpAll || m_doDumpxAODMuonContainer) {
173 StatusCode sc = dumpxAODMuonContainer();
174 if(sc.isFailure()) {
175 ATH_MSG_ERROR("The method dumpxAODMuonContainer() failed");
176 }
177 }
178
179 if(m_doDumpAll || m_doDumpxAODTrigElectronContainer){
180 StatusCode sc = dumpxAODTrigElectronContainer();
181 if (sc.isFailure()) {
182 ATH_MSG_ERROR("The method dumpxAODTrigElectronContainer() failed");
183 }
184 }
185
186 if(m_doDumpAll || m_doDumpxAODTrigPhotonContainer){
187 StatusCode sc = dumpxAODTrigPhotonContainer();
188 if (sc.isFailure()) {
189 ATH_MSG_ERROR("The method dumpxAODTrigElectronContainer() failed");
190 }
191 }
192
193 if(m_doDumpAll || m_doDumpxAODElectronContainer){
194 StatusCode sc = dumpxAODElectronContainer();
195 if (sc.isFailure()) {
196 ATH_MSG_ERROR("The method dumpxAODTrigElectronContainer() failed");
197 }
198 }
199
200 if(m_doDumpAll || m_doDumpxAODPhotonContainer){
201 StatusCode sc = dumpxAODPhotonContainer();
202 if (sc.isFailure()) {
203 ATH_MSG_ERROR("The method dumpxAODTrigElectronContainer() failed");
204 }
205 }
206
207 if(m_doDumpAll || m_doDumpxAODTauJetContainer){
208 StatusCode sc = dumpxAODTauJetContainer();
209 if (sc.isFailure()) {
210 ATH_MSG_ERROR("The method dumpxAODTauJetContainer() failed");
211 }
212 }
213
214 if(m_doDumpAll || m_doDumpxAODTrackParticle){
215 StatusCode sc = dumpxAODTrackParticle();
216 if (sc.isFailure()) {
217 ATH_MSG_ERROR("The method dumpxAODTrackParticle() failed");
218 }
219 }
220
221 if(m_doDumpAll || m_doDumpxAODVertex){
222 StatusCode sc = dumpxAODVertex();
223 if (sc.isFailure()) {
224 ATH_MSG_ERROR("The method dumpxAODVertex() failed");
225 }
226 }
227
228 if (m_doDumpAll || m_doDumpxAODTrigMinBias){
229 StatusCode sc = dumpxAODTrigMinBias();
230 if (sc.isFailure()) {
231 ATH_MSG_ERROR("The method dumpxAODTrigMinBias() failed");
232 }
233 }
234
235 if (m_doDumpTrigPassBits){
236 StatusCode sc = dumpTrigPassBits();
237 if (sc.isFailure()) {
238 ATH_MSG_ERROR("The method dumpTrigPassBits() failed");
239 }
240 }
241
242 if (m_doDumpAll || m_doDumpStoreGate) {
243 ATH_MSG_DEBUG(evtStore()->dump());
244 }
245
246 if (m_doDumpAll || m_doDumpNavigation) {
247 StatusCode sc = dumpNavigation(ctx);
248 if ( sc.isFailure() ) {
249 ATH_MSG_ERROR("The method dumpNavigation() failed");
250 }
251 }
252
253 if (m_doDumpAll || m_doTDTCheck) {
254 ATH_CHECK(dumpTDT(ctx));
255 }
256
257 if (m_doDumpAll || m_doDumpAllTrigComposite || m_dumpTrigCompositeContainers.size() > 0) {
258 ATH_CHECK( dumpTrigComposite() );
259 }
260
261 if (m_doDumpAll || m_doDumpTrigCompsiteNavigation) {
262 std::string trigCompositeSteering;
263 bool pass;
264 ATH_CHECK(TrigCompositeNavigationToDot(trigCompositeSteering, pass));
265 const std::string evtNumber = std::to_string(ctx.eventID().event_number());
266 const std::string passStr = (pass ? "Pass" : "Fail");
267 std::ofstream ofile(std::string("NavGraph_" + m_dumpNavForChain + "_Ev" + evtNumber + "_" + passStr + ".dot").c_str());
268 ofile << trigCompositeSteering;
269 }
270
271
272
273 ATH_MSG_INFO( " ==========END of event===========" );
274 return StatusCode::SUCCESS;
275
276}
277
279StatusCode TrigEDMChecker::dumpTrigPassBits(){
280 const std::string name="HLT_xAOD__TrigPassBitsContainer_passbits";
281 const xAOD::TrigPassBitsContainer *xbitscont=nullptr;
282 StatusCode sc = evtStore()->retrieve(xbitscont,name);
283 if (sc.isFailure() ){
284 ATH_MSG_INFO("Cannot retrieve TrigPassBits");
285 }
286 else {
287 ATH_MSG_INFO("Size of PassBits container : " << xbitscont->size());
288 for(const auto bits:*xbitscont){
289 if(bits==nullptr){
290 ATH_MSG_INFO("TrigPassBits point nullptr ");
291 continue;
292 }
293 ATH_MSG_DEBUG("Analyzing bits for " << bits->containerClid() << " of size " << bits->size() << " with bit size " << bits->passBits().size());
294 }
295
296 for(const xAOD::TrigPassBits* bits : *xbitscont){
297 if(bits==nullptr){
298 ATH_MSG_INFO("TrigPassBits point nullptr ");
299 continue;
300 }
301 ATH_MSG_DEBUG("Analyzing bits for " << bits->containerClid() << " of size " << bits->size() << " with bit size " << bits->passBits().size());
302 }
303 }
304 return StatusCode::SUCCESS;
305}
306
307void TrigEDMChecker::dumpTrigSpacePointCounts()
308{
309 ATH_MSG_INFO("MinBias in dumpTrigSpacePointCounts()");
310
311 std::string METTag="HLT_xAOD__TrigSpacePointCountsContainer_spacepoints";
312
313 const xAOD::TrigSpacePointCountsContainer* SpacePointCountsCont=0;
314 StatusCode sc = evtStore()->retrieve(SpacePointCountsCont,METTag);
315
316 if (sc.isFailure())
317 ATH_MSG_INFO("failed to retrieve " << METTag);
318 else {
319 ATH_MSG_INFO("Accessing " << METTag << " with " << SpacePointCountsCont->size() << " elements");
320
321 std::string s; char buff[128];
322 std::vector<float> getVec;
323 float sum;
324
325 // Loop over container content
326 for(uint i = 0; i < SpacePointCountsCont->size(); i++) {
327 getVec = SpacePointCountsCont->at(i)->contentsPixelClusEndcapC();
328 sum = 0.;
329 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
330 snprintf(buff, sizeof(buff), "REGTEST %s SUM of contentsPixelClusEndcapC() = %10.2f ", s.c_str(), sum );
331 ATH_MSG_INFO(buff);
332
333 getVec = SpacePointCountsCont->at(i)->contentsPixelClusBarrel();
334 sum = 0.;
335 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
336 snprintf(buff, sizeof(buff), "REGTEST %s SUM of contentsPixelClusBarrel() = %10.2f ", s.c_str(), sum );
337 ATH_MSG_INFO(buff);
338
339 getVec = SpacePointCountsCont->at(i)->contentsPixelClusEndcapA();
340 sum = 0.;
341 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
342 snprintf(buff, sizeof(buff), "REGTEST %s SUM of contentsPixelClusEndcapA() = %10.2f ", s.c_str(), sum );
343 ATH_MSG_INFO(buff);
344
345 snprintf(buff, sizeof(buff), "REGTEST %s pixelClusTotBins() = %u ", s.c_str(), SpacePointCountsCont->at(i)->pixelClusTotBins() );
346 ATH_MSG_INFO(buff);
347
348 snprintf(buff, sizeof(buff), "REGTEST %s pixelClusTotMin() = %10.2f ", s.c_str(), SpacePointCountsCont->at(i)->pixelClusTotMin() );
349 ATH_MSG_INFO(buff);
350
351 snprintf(buff, sizeof(buff), "REGTEST %s pixelClusTotMax() = %10.2f ", s.c_str(), SpacePointCountsCont->at(i)->pixelClusTotMax() );
352 ATH_MSG_INFO(buff);
353
354 snprintf(buff, sizeof(buff), "REGTEST %s pixelClusSizeBins() = %u ", s.c_str(), SpacePointCountsCont->at(i)->pixelClusSizeBins() );
355 ATH_MSG_INFO(buff);
356
357 snprintf(buff, sizeof(buff), "REGTEST %s pixelClusSizeMin() = %10.2f ", s.c_str(), SpacePointCountsCont->at(i)->pixelClusSizeMin() );
358 ATH_MSG_INFO(buff);
359
360 snprintf(buff, sizeof(buff), "REGTEST %s pixelClusSizeMax() = %10.2f ", s.c_str(), SpacePointCountsCont->at(i)->pixelClusSizeMax() );
361 ATH_MSG_INFO(buff);
362
363 snprintf(buff, sizeof(buff), "REGTEST %s sctSpEndcapC() = %u ", s.c_str(), SpacePointCountsCont->at(i)->sctSpEndcapC() );
364 ATH_MSG_INFO(buff);
365
366 snprintf(buff, sizeof(buff), "REGTEST %s sctSpBarrel() = %u ", s.c_str(), SpacePointCountsCont->at(i)->sctSpBarrel() );
367 ATH_MSG_INFO(buff);
368
369 snprintf(buff, sizeof(buff), "REGTEST %s sctSpEndcapA() = %u ", s.c_str(), SpacePointCountsCont->at(i)->sctSpEndcapA() );
370 ATH_MSG_INFO(buff);
371 }
372 }
373}
374
375void TrigEDMChecker::dumpTrigT2MBTSBits(){
376 ATH_MSG_INFO("MinBias in dumpTrigT2MBTSBits()");
377
378 std::string METTag="HLT_xAOD__TrigT2MbtsBitsContainer_T2Mbts";
379
380 const xAOD::TrigT2MbtsBitsContainer* T2MbtsBitsCont=0;
381 StatusCode sc = evtStore()->retrieve(T2MbtsBitsCont,METTag);
382
383 if (sc.isFailure())
384 ATH_MSG_INFO("failed to retrieve " << METTag);
385 else {
386 ATH_MSG_INFO("Accessing " << METTag << " with " << T2MbtsBitsCont->size() << " elements");
387
388 std::string s; char buff[380];
389 std::vector<float> getVec;
390 float sum;
391
392 // Loop over container content
393 for(uint i = 0; i < T2MbtsBitsCont->size(); i++) {
394 getVec = T2MbtsBitsCont->at(i)->triggerEnergies();
395 sum = 0.;
396 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
397 snprintf(buff, sizeof(buff), "REGTEST %s SUM of triggerEnergies() = %10.2f ", s.c_str(), sum );
398 ATH_MSG_INFO(buff);
399
400 getVec = T2MbtsBitsCont->at(i)->triggerTimes();
401 sum = 0.;
402 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
403 snprintf(buff, sizeof(buff), "REGTEST %s SUM of triggerTimes() = %10.2f ", s.c_str(), sum );
404 ATH_MSG_INFO(buff);
405 }
406 }
407}
408
409void TrigEDMChecker::dumpTrigVertexCounts(){
410 ATH_MSG_INFO("MinBias in dumpTrigVertexCounts()");
411
412 std::string METTag="HLT_xAOD__TrigVertexCountsContainer_vertexcounts";
413
414 const xAOD::TrigVertexCountsContainer* T2VertexCountsCont=0;
415 StatusCode sc = evtStore()->retrieve(T2VertexCountsCont,METTag);
416
417 if (sc.isFailure())
418 ATH_MSG_INFO("failed to retrieve " << METTag);
419 else {
420 ATH_MSG_INFO("Accessing " << METTag << " with " << T2VertexCountsCont->size() << " elements");
421
422 std::string s; char buff[380];
423 std::vector<float> fgetVec;
424 float fsum(0.);
425 std::vector<unsigned int> ugetVec;
426 unsigned int usum(0);
427
428 // Loop over container content
429 for(uint i = 0; i < T2VertexCountsCont->size(); i++) {
430 ugetVec = T2VertexCountsCont->at(i)->vtxNtrks();
431 for (uint j = 0; j < ugetVec.size(); ++j) usum += ugetVec[j];
432 snprintf(buff, sizeof(buff), "REGTEST %s SUM of vtxNtrks() = %u ", s.c_str(), usum );
433 ATH_MSG_INFO(buff);
434
435 fgetVec = T2VertexCountsCont->at(i)->vtxTrkPtSqSum();
436 for (uint j = 0; j < fgetVec.size(); ++j) fsum += fgetVec[j];
437 snprintf(buff, sizeof(buff), "REGTEST %s SUM of vtxTrkPtSqSum() = %10.2f ", s.c_str(), fsum );
438 ATH_MSG_INFO(buff);
439 }
440 }
441}
442
443void TrigEDMChecker::dumpTrigTrackCounts(){
444 ATH_MSG_INFO("MinBias in dumpTrigTrackCounts()");
445
446 std::string METTag="HLT_xAOD__TrigTrackCountsContainer_trackcounts";
447
448 const xAOD::TrigTrackCountsContainer* T2TrackCountsCont=0;
449 StatusCode sc = evtStore()->retrieve(T2TrackCountsCont,METTag);
450
451 if (sc.isFailure())
452 ATH_MSG_INFO("failed to retrieve " << METTag);
453 else {
454 ATH_MSG_INFO("Accessing " << METTag << " with " << T2TrackCountsCont->size() << " elements");
455
456 std::string s; char buff[380];
457 std::vector<float> getVec;
458 float sum;
459
460 // Loop over container content
461 for(uint i = 0; i < T2TrackCountsCont->size(); i++) {
462 getVec = T2TrackCountsCont->at(i)->z0_pt();
463 sum = 0.;
464 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
465 snprintf(buff, sizeof(buff), "REGTEST %s SUM of z0_pt = %10.2f ", s.c_str(), sum );
466 ATH_MSG_INFO(buff);
467
468 getVec = T2TrackCountsCont->at(i)->eta_phi();
469 sum = 0.;
470 for (uint j = 0; j < getVec.size(); ++j) sum += getVec[j];
471 snprintf(buff, sizeof(buff), "REGTEST %s SUM of eta_phi() = %10.2f ", s.c_str(), sum );
472 ATH_MSG_INFO(buff);
473
474 snprintf(buff, sizeof(buff), "REGTEST %s z0Bins() = %u ", s.c_str(), T2TrackCountsCont->at(i)->z0Bins() );
475 ATH_MSG_INFO(buff);
476
477 snprintf(buff, sizeof(buff), "REGTEST %s z0Min() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->z0Min() );
478 ATH_MSG_INFO(buff);
479
480 snprintf(buff, sizeof(buff), "REGTEST %s z0Max() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->z0Max() );
481 ATH_MSG_INFO(buff);
482
483 snprintf(buff, sizeof(buff), "REGTEST %s ptBins() = %u ", s.c_str(), T2TrackCountsCont->at(i)->ptBins() );
484 ATH_MSG_INFO(buff);
485
486 snprintf(buff, sizeof(buff), "REGTEST %s ptMin() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->ptMin() );
487 ATH_MSG_INFO(buff);
488
489 snprintf(buff, sizeof(buff), "REGTEST %s ptMax() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->ptMax() );
490 ATH_MSG_INFO(buff);
491
492 snprintf(buff, sizeof(buff), "REGTEST %s etaBins() = %u ", s.c_str(), T2TrackCountsCont->at(i)->etaBins() );
493 ATH_MSG_INFO(buff);
494
495 snprintf(buff, sizeof(buff), "REGTEST %s etaMin() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->etaMin() );
496 ATH_MSG_INFO(buff);
497
498 snprintf(buff, sizeof(buff), "REGTEST %s etaMax() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->etaMax() );
499 ATH_MSG_INFO(buff);
500
501 snprintf(buff, sizeof(buff), "REGTEST %s phiBins() = %u ", s.c_str(), T2TrackCountsCont->at(i)->phiBins() );
502 ATH_MSG_INFO(buff);
503
504 snprintf(buff, sizeof(buff), "REGTEST %s phiMin() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->phiMin() );
505 ATH_MSG_INFO(buff);
506
507 snprintf(buff, sizeof(buff), "REGTEST %s phiMax() = %10.2f ", s.c_str(), T2TrackCountsCont->at(i)->phiMax() );
508 ATH_MSG_INFO(buff);
509 }
510 }
511}
512
513StatusCode TrigEDMChecker::dumpxAODTrigMinBias() {
514
515 dumpTrigSpacePointCounts();
516 dumpTrigT2MBTSBits();
517 dumpTrigVertexCounts();
518 dumpTrigTrackCounts();
519
520 return StatusCode::SUCCESS;
521}
522
524
525
526
527StatusCode TrigEDMChecker::dumpxAODTrigMissingET() {
528
529 ATH_MSG_INFO("dumpxAODTrigMissingET()");
530
531 int ntag=4;
532 std::string METTags[]={"HLT_xAOD__TrigMissingETContainer_EFJetEtSum","HLT_xAOD__TrigMissingETContainer_TrigEFMissingET", "HLT_xAOD__TrigMissingETContainer_TrigL2MissingET_FEB","HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl"};
533
534 for(int itag=0; itag <ntag; itag++) {
535
536 const xAOD::TrigMissingETContainer* MissingETCont=0;
537 StatusCode sc = evtStore()->retrieve(MissingETCont,METTags[itag]);
538 if (sc.isFailure())
539 ATH_MSG_INFO("failed to retrieve " << METTags[itag]);
540 else {
541 ATH_MSG_INFO("Accessing " << METTags[itag] << " with " << MissingETCont->size() << " elements");
542
543 // Loop over container content
544 for(uint i = 0; i < MissingETCont->size(); i++) {
545
546 std::string s; char buff[3000];
547
548 snprintf(buff, sizeof(buff), "REGTEST %s Ex = %10.2f CLHEP::MeV", s.c_str(), MissingETCont->at(i)->ex() );
549 ATH_MSG_INFO(buff);
550 snprintf(buff, sizeof(buff), "REGTEST %s Ey = %10.2f CLHEP::MeV", s.c_str(), MissingETCont->at(i)->ey() );
551 ATH_MSG_INFO(buff);
552 snprintf(buff, sizeof(buff), "REGTEST %s Ez = %10.2f CLHEP::MeV", s.c_str(), MissingETCont->at(i)->ez() );
553 ATH_MSG_INFO(buff);
554 snprintf(buff, sizeof(buff), "REGTEST %s SumET = %10.2f CLHEP::MeV", s.c_str(), MissingETCont->at(i)->sumEt() );
555 ATH_MSG_INFO(buff);
556 snprintf(buff, sizeof(buff), "REGTEST %s SumE = %10.2f CLHEP::MeV", s.c_str(), MissingETCont->at(i)->sumE() );
557 ATH_MSG_INFO(buff);
558 snprintf(buff, sizeof(buff), "REGTEST %s Flag = %d", s.c_str(), MissingETCont->at(i)->flag() );
559 ATH_MSG_INFO(buff);
560 snprintf(buff, sizeof(buff), "REGTEST %s Flag = %d", s.c_str(), MissingETCont->at(i)->roiWord() );
561 ATH_MSG_INFO(buff);
562
563 unsigned int Nc = MissingETCont->at(i)->getNumberOfComponents();
564 if (Nc > 0) {
565 s="REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
566 s+="/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
567 ATH_MSG_INFO(s);
568 }
569
570 for(uint j = 0; j < Nc; j++) {
571
572 std::string name = MissingETCont->at(i)->nameOfComponent(j);
573 const short status = MissingETCont->at(i)->statusComponent(j);
574 const unsigned short usedChan = MissingETCont->at(i)->usedChannelsComponent(j);
575 const short sumOfSigns = MissingETCont->at(i)->sumOfSignsComponent(j);
576 const float calib0 = MissingETCont->at(i)->calib0Component(j);
577 const float calib1 = MissingETCont->at(i)->calib1Component(j);
578 const float ex = MissingETCont->at(i)->exComponent(j);
579 const float ey = MissingETCont->at(i)->eyComponent(j);
580 const float ez = MissingETCont->at(i)->ezComponent(j);
581 const float sumE = MissingETCont->at(i)->sumEComponent(j);
582 const float sumEt = MissingETCont->at(i)->sumEtComponent(j);
583
584 snprintf(buff, sizeof(buff),
585 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
586 name.c_str(), status, usedChan, sumOfSigns, calib1, calib0,
587 ex, ey, ez, sumE, sumEt);
588 ATH_MSG_INFO(buff);
589 }
590 }
591 }
592 }
593
594 return StatusCode::SUCCESS;
595}
596
598
599StatusCode TrigEDMChecker::dumpLVL1_ROI() {
600
601 ATH_MSG_DEBUG("in dumpLVL1_ROI()");
602
603 ATH_MSG_INFO("REGTEST ==========START of LVL1_ROI DUMP===========");
604
605 const LVL1_ROI * lvl1ROI;
606 StatusCode sc = evtStore()->retrieve(lvl1ROI);
607 if (sc.isFailure() ) {
608 ATH_MSG_INFO("REGTEST No LVL1_ROI found");
609 return StatusCode::SUCCESS;
610 }
611
612 ATH_MSG_INFO("REGTEST LVL1_ROI retrieved");
613
614 LVL1_ROI::emtaus_type::const_iterator itEMTau =
615 (lvl1ROI->getEmTauROIs()).begin();
616 LVL1_ROI::emtaus_type::const_iterator itEMTau_e =
617 (lvl1ROI->getEmTauROIs()).end();
618 int j = 0;
619 for( ; itEMTau != itEMTau_e; ++itEMTau, ++j) {
620 ATH_MSG_INFO("REGTEST Looking at LVL1_ROI " << j);
621 ATH_MSG_INFO("REGTEST LVL1 EmTauROI Eta is " << itEMTau->getEta());
622 ATH_MSG_INFO("REGTEST LVL1 EmTauROI Phi is " << itEMTau->getPhi());
623 ATH_MSG_INFO("REGTEST LVL1 EmTauROI EMClus is " << itEMTau->getEMClus());
624 ATH_MSG_INFO("REGTEST LVL1 EmTauROI TauClus is " << itEMTau->getTauClus());
625 ATH_MSG_INFO("REGTEST LVL1 EmTauROI EMIsol is " << itEMTau->getEMIsol());
626 ATH_MSG_INFO("REGTEST LVL1 EmTauROI HadIsol is " << itEMTau->getHadIsol());
627 ATH_MSG_INFO("REGTEST LVL1 EmTauROI Core is " << itEMTau->getCore());
628 ATH_MSG_INFO("REGTEST LVL1 EmTauROI HadCore is " << itEMTau->getHadCore());
629 ATH_MSG_INFO("REGTEST LVL1 EmTauROI roiWord is " << itEMTau->getROIWord());
630 }
631
632 ATH_MSG_INFO("REGTEST ==========END of LVL1_ROI DUMP===========");
633 ATH_MSG_DEBUG("dumpLVL1_ROI() succeeded");
634 return StatusCode::SUCCESS;
635}
636
638
639StatusCode TrigEDMChecker::dumpxAODMuonContainer() {
640
641 ATH_MSG_DEBUG("In dumpxAODMuonContainer");
642
643 ATH_MSG_INFO( "REGTEST ==========START of xAOD::MuonContainer DUMP===========" );
644
645 const xAOD::MuonContainer* muonCont=0;
646 StatusCode sc = evtStore()->retrieve(muonCont,"HLT_xAOD__MuonContainer_MuonEFInfo");
647 if (sc.isFailure()) {
648 ATH_MSG_INFO("REGTEST No muon container HLT_xAOD__MuonContainer_MuonEFInfo");
649 return StatusCode::SUCCESS;
650 }
651
652 std::string output = m_muonPrinter->print( *muonCont );
653 msg(MSG::INFO) << output << endmsg;
654
655 ATH_MSG_INFO( "REGTEST ==========END of xAOD::MuonContainer DUMP===========" );
656
657 return StatusCode::SUCCESS;
658
659}
660
662
663StatusCode TrigEDMChecker::dumpxAODTrigElectronContainer() {
664
665 ATH_MSG_DEBUG("In dumpxAODElectronContainer");
666
667 ATH_MSG_INFO( "REGTEST ==========START of xAOD::TrigElectronContainer DUMP===========" );
668
669 const xAOD::TrigElectronContainer* elCont=0;
670 StatusCode sc = evtStore()->retrieve(elCont,"HLT_xAOD__TrigElectronContainer_L2ElectronFex");
671 if (sc.isFailure()) {
672 ATH_MSG_INFO("REGTEST No TrigElectron container HLT_xAOD__TrigElectronContainer_L2ElectronFex");
673 return StatusCode::SUCCESS;
674 }
675
676 for (const auto eg : *elCont){
677 ATH_MSG_INFO("REGTEST TrigElectron->Phi() returns " << eg->phi());
678 ATH_MSG_INFO("REGTEST TrigElectron->Eta() returns " << eg->eta());
679 ATH_MSG_INFO("REGTEST TrigElectron->rEta returns " << eg->rcore());
680 ATH_MSG_INFO("REGTEST TrigElectron->eratio() returns " << eg->eratio());
681 ATH_MSG_INFO("REGTEST TrigElectron->pt() returns " << eg->pt());
682 ATH_MSG_INFO("REGTEST TrigElectron->etHad() returns " << eg->etHad());
683 ATH_MSG_INFO("REGTEST TrigElectron->f1() returns " << eg->f1());
684 ATH_MSG_INFO("REGTEST TrigElectron caloEta = " << eg->caloEta());
685 ATH_MSG_INFO("REGTEST TrigElectron dPhiCalo" << eg->trkClusDphi());
686 ATH_MSG_INFO("REGTEST TrigElectron dEtaCalo" << eg->trkClusDeta());
687 ATH_MSG_INFO("REGTEST TrigElectron pTcalo" << eg->pt());
688 ATH_MSG_INFO("REGTEST TrigElectron eTOverPt" << eg->etOverPt());
689 ATH_MSG_INFO("REGTEST TrigElectron nTRTHits" << eg->nTRTHits());
690 ATH_MSG_INFO("REGTEST TrigElectron nStrawHits" << eg->nTRTHiThresholdHits());
691 ATH_MSG_INFO("REGTEST TrigElectron Check EMCluster");
692 if(eg->emCluster()){
693 ATH_MSG_INFO("REGTEST TrigElectron EMCluster retrieved");
694 ATH_MSG_INFO("REGTEST TrigElectron emCluster->energy() returns " << eg->emCluster()->energy());
695 ATH_MSG_INFO("REGTEST TrigElectron emCluster->phi() returns " << eg->emCluster()->phi());
696 ATH_MSG_INFO("REGTEST TrigElectron emCluster->eta() returns " << eg->emCluster()->eta());
697 ATH_MSG_INFO("REGTEST TrigElectron emCluster check Element Link");
698 ATH_MSG_INFO("REGTEST TrigElectron emCluster energy = " << eg->emCluster()->energy());
699 ATH_MSG_INFO("REGTEST TrigElectron ElementLink emCluster energy = " << (*eg->emClusterLink())->energy());
700 }
701 else ATH_MSG_INFO("REGTEST TrigElectron No EMCluster retrieved!");
702 ATH_MSG_INFO("REGTEST TrigElectron Check TrackParticle");
703 if(eg->trackParticle()){
704 ATH_MSG_INFO("REGTEST TrigElectron TrackParticle retrieved");
705 ATH_MSG_INFO("REGTEST TrigElectron trackParticle->pt() returns " << eg->trackParticle()->pt());
706 ATH_MSG_INFO("REGTEST TrigElectron trackParticle->phi() returns " << eg->trackParticle()->phi());
707 ATH_MSG_INFO("REGTEST TrigElectron trackParticle->eta() returns " << eg->trackParticle()->eta());
708 ATH_MSG_INFO("REGTEST TrigElectron check TrackParticle Element Link");
709 ATH_MSG_INFO("REGTEST TrigElectron TrackParticle pt = " << eg->trackParticle()->pt());
710 ATH_MSG_INFO("REGTEST TrigElectron ElementLink TrackParticle pt = " << (*eg->trackParticleLink())->pt());
711 }
712 else ATH_MSG_INFO("REGTEST TrigElectron No TrackParticle retrieved!");
713 }
714 ATH_MSG_INFO( "REGTEST ==========END of xAOD::TrigElectronContainer DUMP===========" );
715
716 return StatusCode::SUCCESS;
717
718}
719
721StatusCode TrigEDMChecker::dumpxAODTrigPhotonContainer() {
722
723 ATH_MSG_DEBUG("In dumpxAODTrigPhotonContainer");
724
725 ATH_MSG_INFO( "REGTEST ==========START of xAOD::TrigPhotonContainer DUMP===========" );
726
727 const xAOD::TrigPhotonContainer* phCont=0;
728 StatusCode sc = evtStore()->retrieve(phCont,"HLT_xAOD__TrigPhotonContainer_L2PhotonFex");
729 if (sc.isFailure()) {
730 ATH_MSG_INFO("REGTEST No TrigPhoton container HLT_xAOD__TrigPhotonContainer_L2PhotonFex");
731 return StatusCode::SUCCESS;
732 }
733
734 for (const auto eg : *phCont){
735
736 ATH_MSG_INFO("REGTEST TrigPhoton->Phi() returns " << eg->phi());
737 ATH_MSG_INFO("REGTEST TrigPhoton->Eta() returns " << eg->eta());
738 ATH_MSG_INFO("REGTEST TrigPhoton->dPhi() returns " << eg->dPhi());
739 ATH_MSG_INFO("REGTEST TrigPhoton->dEta() returns " << eg->dEta());
740 ATH_MSG_INFO("REGTEST TrigPhoton->rEta returns " << eg->rcore());
741 ATH_MSG_INFO("REGTEST TrigPhoton->eratio() returns " << eg->eratio());
742 ATH_MSG_INFO("REGTEST TrigPhoton->pt() returns " << eg->pt());
743 ATH_MSG_INFO("REGTEST TrigPhoton->etHad() returns " << eg->etHad());
744 ATH_MSG_INFO("REGTEST TrigPhoton->f1() returns " << eg->f1());
745 ATH_MSG_INFO("REGTEST TrigPhoton Check EMCluster");
746 if(eg->emCluster()){
747 ATH_MSG_INFO("REGTEST TrigPhoton EMCluster retrieved");
748 ATH_MSG_INFO("REGTEST TrigPhoton emCluster->energy() returns " << eg->emCluster()->energy());
749 ATH_MSG_INFO("REGTEST TrigPhoton emCluster->phi() returns " << eg->emCluster()->phi());
750 ATH_MSG_INFO("REGTEST TrigPhoton emCluster->eta() returns " << eg->emCluster()->eta());
751 ATH_MSG_INFO("REGTEST TrigPhoton emCluster check Element Link");
752 ATH_MSG_INFO("REGTEST TrigPhoton emCluster energy = " << eg->emCluster()->energy());
753 ATH_MSG_INFO("REGTEST TrigPhoton ElementLink emCluster energy = " << (*eg->emClusterLink())->energy());
754 }
755 else ATH_MSG_INFO("REGTEST TrigPhoton No EMCluster retrieved!");
756 }
757 ATH_MSG_INFO( "REGTEST ==========END of xAOD::TrigPhotonContainer DUMP===========" );
758
759 return StatusCode::SUCCESS;
760
761}
762
764StatusCode TrigEDMChecker::dumpxAODElectronContainer() {
765
766 ATH_MSG_DEBUG("In dumpxAODElectronContainer");
767
768 ATH_MSG_INFO( "REGTEST ==========START of xAOD::ElectronContainer DUMP===========" );
769
770 const xAOD::ElectronContainer* elCont=0;
771 StatusCode sc = evtStore()->retrieve(elCont,"HLT_xAOD__ElectronContainer_egamma_Electrons");
772 if (sc.isFailure()) {
773 ATH_MSG_INFO("REGTEST No Electron container HLT_xAOD__ElectronContainer_egamma_Electrons");
774 return StatusCode::SUCCESS;
775 }
776 float val_float=-99;
777 unsigned int isEMbit=0;
778 bool pid=false;
779 //DEBUG output for Egamma container
780 ATH_MSG_INFO(" REGTEST: xAOD Reconstruction variables: ");
781 // //Cluster and ShowerShape info
782 //
783 static const SG::AuxElement::Accessor< float > accLH("LHValue");
784 static const SG::AuxElement::Accessor< float > accLHCalo("LHCaloValue");
785 static const SG::AuxElement::Accessor<ElementLink<xAOD::CaloClusterContainer> > orig ("originalCaloCluster");
786 for (const auto eg : *elCont){
787 //REGTEST printout
788 if (eg) {
789 ATH_MSG_INFO(" REGTEST: egamma energy: " << eg->e() );
790 ATH_MSG_INFO(" REGTEST: egamma eta: " << eg->eta() );
791 ATH_MSG_INFO(" REGTEST: egamma phi: " << eg->phi() );
792 if(eg->selectionisEM(isEMbit,"isEMVLoose"))
793 ATH_MSG_INFO(" REGTEST: isEMVLoose " << std::hex << isEMbit << std::dec);
794 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
795 if(eg->selectionisEM(isEMbit,"isEMLoose"))
796 ATH_MSG_INFO(" REGTEST: isEMLoose " << std::hex << isEMbit << std::dec);
797 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
798 if(eg->selectionisEM(isEMbit,"isEMMedium"))
799 ATH_MSG_INFO(" REGTEST: isEMMedium " << std::hex << isEMbit << std::dec);
800 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
801 if(eg->selectionisEM(isEMbit,"isEMTight"))
802 ATH_MSG_INFO(" REGTEST: isEMTight " << std::hex << isEMbit << std::dec);
803 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
804 if(eg->selectionisEM(isEMbit,"isEMLHVLoose"))
805 ATH_MSG_INFO(" REGTEST: isEMLHVLoose " << std::hex << isEMbit << std::dec);
806 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
807 if(eg->selectionisEM(isEMbit,"isEMLHLoose"))
808 ATH_MSG_INFO(" REGTEST: isEMLHLoose " << std::hex << isEMbit << std::dec);
809 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
810 if(eg->selectionisEM(isEMbit,"isEMLHMedium"))
811 ATH_MSG_INFO(" REGTEST: isEMLHMedium " << std::hex << isEMbit << std::dec);
812 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
813 if(eg->selectionisEM(isEMbit,"isEMLHTight"))
814 ATH_MSG_INFO(" REGTEST: isEMLHTight " << std::hex << isEMbit << std::dec);
815 else ATH_MSG_WARNING(" REGTEST: Missing Aux info");
816 if(accLH.isAvailable(*eg))
817 ATH_MSG_INFO(" REGTEST: LHValue " << accLH(*eg));
818 else
819 ATH_MSG_WARNING(" REGTEST: Missing Aux info ");
820 if(accLHCalo.isAvailable(*eg))
821 ATH_MSG_INFO(" REGTEST: LHValue " << accLHCalo(*eg));
822 else
823 ATH_MSG_WARNING(" REGTEST: Missing Aux info ");
824 if(eg->passSelection(pid,"LHVLoose"))
825 ATH_MSG_INFO(" REGTEST: LHVLoose " << pid);
826 else
827 ATH_MSG_WARNING(" REGTEST: Missing Aux info ");
828 if(eg->passSelection(pid,"LHLoose"))
829 ATH_MSG_INFO(" REGTEST: LHLoose " << pid);
830 else
831 ATH_MSG_WARNING(" REGTEST: Missing Aux info ");
832 if(eg->passSelection(pid,"LHMedium"))
833 ATH_MSG_INFO(" REGTEST: LHMedium " << pid);
834 else
835 ATH_MSG_WARNING(" REGTEST: Missing Aux info ");
836 if(eg->passSelection(pid,"LHTight"))
837 ATH_MSG_INFO(" REGTEST: LHTight " << pid);
838 else
839 ATH_MSG_WARNING(" REGTEST: Missing Aux info ");
840 } else{
841 ATH_MSG_INFO(" REGTEST: problems with egamma pointer" );
842 return StatusCode::SUCCESS;
843 }
844 ATH_MSG_INFO(" REGTEST: caloCluster variables ");
845 if (eg->caloCluster()) {
846 ATH_MSG_INFO(" REGTEST: egamma cluster transverse energy: " << eg->caloCluster()->et() );
847 ATH_MSG_INFO(" REGTEST: egamma cluster eta: " << eg->caloCluster()->eta() );
848 ATH_MSG_INFO(" REGTEST: egamma cluster phi: " << eg->caloCluster()->phi() );
849 double tmpeta = -999.;
850 double tmpphi = -999.;
851 eg->caloCluster()->retrieveMoment(xAOD::CaloCluster::ETACALOFRAME,tmpeta);
852 eg->caloCluster()->retrieveMoment(xAOD::CaloCluster::PHICALOFRAME,tmpphi);
853 ATH_MSG_INFO(" REGTEST: egamma cluster calo-frame coords. etaCalo = " << tmpeta);
854 ATH_MSG_INFO(" REGTEST: egamma cluster calo-frame coords. phiCalo = " << tmpphi);
855 } else{
856 ATH_MSG_INFO(" REGTEST: problems with egamma cluster pointer" );
857 }
858 ATH_MSG_INFO("REGTEST: Check the original (uncalibrated)");
859 if (!orig.isAvailable(*eg->caloCluster()) || !orig(*eg->caloCluster()).isValid()){
860 ATH_MSG_INFO("Problem with original cluster link");
861 }
862 else {
863 const xAOD::CaloCluster *origClus = *orig(*eg->caloCluster());
864 ATH_MSG_INFO("REGTEST:: Compare new and old clusters");
865 ATH_MSG_INFO("REGTEST:: Original Cluster e,eta,phi" << origClus->e() << " " << origClus->eta() << " " << origClus->phi());
866 ATH_MSG_INFO("REGTEST:: MVA Cluster e,eta,phi" << eg->caloCluster()->e() << " " << eg->caloCluster()->eta() << " " << eg->caloCluster()->phi());
867 }
868 ATH_MSG_INFO(" REGTEST: trackmatch variables ");
869 if(eg->trackParticle()){
870 ATH_MSG_INFO(" REGTEST: pt= " << eg->trackParticle()->pt());
871 ATH_MSG_INFO(" REGTEST: charge= " << eg->trackParticle()->charge());
872 ATH_MSG_INFO(" REGTEST: E/p= " << eg->caloCluster()->et() / eg->trackParticle()->pt() );
873 eg->trackCaloMatchValue(val_float,xAOD::EgammaParameters::deltaEta1);
874 ATH_MSG_INFO(" REGTEST: Delta eta 1st sampling= " << val_float);
875 eg->trackCaloMatchValue(val_float,xAOD::EgammaParameters::deltaPhi2);
876 ATH_MSG_INFO(" REGTEST: Delta phi 2nd sampling= " << val_float);
877 } else{
878 ATH_MSG_INFO(" REGTEST: no electron eg->trackParticle() pointer");
879 }
880
881 //msg() << MSG::VERBOSE << " REGTEST: cluster variables " << endmsg;
882 //clus = eg->caloCluster();
883 ATH_MSG_INFO(" REGTEST: EMShower variables ");
884 eg->showerShapeValue(val_float,xAOD::EgammaParameters::ethad);
885 ATH_MSG_INFO(" REGTEST: ethad = " << val_float);
886 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e011);
887 ATH_MSG_INFO(" REGTEST: e011 = " << val_float);
888 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e132);
889 ATH_MSG_INFO(" REGTEST: e132 = " << val_float);
890 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e237);
891 ATH_MSG_INFO(" REGTEST: e237 = " << val_float);
892 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e335);
893 ATH_MSG_INFO(" REGTEST: e335 = " << val_float);
894 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e2ts1);
895 ATH_MSG_INFO(" REGTEST: e2ts1 = " << val_float);
896 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e2tsts1);
897 ATH_MSG_INFO(" REGTEST: e2tsts1 = " << val_float);
898 eg->isolationValue(val_float,xAOD::Iso::ptcone20);
899 ATH_MSG_INFO(" REGTEST: ptcone20 = " << val_float);
900 eg->isolationValue(val_float,xAOD::Iso::ptcone30);
901 ATH_MSG_INFO(" REGTEST: ptcone30 = " << val_float);
902 eg->isolationValue(val_float,xAOD::Iso::ptcone40);
903 ATH_MSG_INFO(" REGTEST: ptcone40 = " << val_float);
904 eg->isolationValue(val_float,xAOD::Iso::etcone20);
905 ATH_MSG_INFO(" REGTEST: etcone20 = " << val_float);
906 eg->isolationValue(val_float,xAOD::Iso::etcone30);
907 ATH_MSG_INFO(" REGTEST: etcone30 = " << val_float);
908 eg->isolationValue(val_float,xAOD::Iso::etcone40);
909 ATH_MSG_INFO(" REGTEST: etcone40 = " << val_float);
910 //DEBUG info for Electrons which by definition have a track match
911
912 }
913 ATH_MSG_INFO( "REGTEST ==========END of xAOD::ElectronContainer DUMP===========" );
914
915 return StatusCode::SUCCESS;
916
917}
918
920StatusCode TrigEDMChecker::dumpxAODPhotonContainer() {
921
922 ATH_MSG_DEBUG("In dumpxAODPhotonContainer");
923
924 ATH_MSG_INFO( "REGTEST ==========START of xAOD::PhotonContainer DUMP===========" );
925
926 const xAOD::PhotonContainer* phCont=0;
927 StatusCode sc = evtStore()->retrieve(phCont,"HLT_xAOD__PhotonContainer_egamma_Photons");
928 if (sc.isFailure()) {
929 ATH_MSG_INFO("REGTEST No Photon container HLT_xAOD__PhotonContainer_egamma_Photons");
930 return StatusCode::SUCCESS;
931 }
932
933 float val_float=-99;
934 unsigned int isEMbit=0;
935 //DEBUG output for xAOD::PhotonContainer
936 ATH_MSG_INFO(" REGTEST: xAOD Reconstruction variables: ");
937 //Cluster and ShowerShape info
938 static const SG::AuxElement::Accessor<ElementLink<xAOD::CaloClusterContainer> > orig ("originalCaloCluster");
939 for (const auto eg : *phCont){
940 //REGTEST printout
941 if (eg) {
942 ATH_MSG_INFO(" REGTEST: egamma energy: " << eg->e() );
943 ATH_MSG_INFO(" REGTEST: egamma eta: " << eg->eta() );
944 ATH_MSG_INFO(" REGTEST: egamma phi: " << eg->phi() );
945 ATH_MSG_INFO(" REGTEST: isEMLoose " << eg->selectionisEM(isEMbit,"isEMLoose"));
946 ATH_MSG_INFO(" REGTEST: isEMLoose bit " << std::hex << isEMbit << std::dec);
947 ATH_MSG_INFO(" REGTEST: isEMMedium " << eg->selectionisEM(isEMbit,"isEMMedium"));
948 ATH_MSG_INFO(" REGTEST: isEMMedium bit " << std::hex << isEMbit << std::dec);
949 ATH_MSG_INFO(" REGTEST: isEMTight " << eg->selectionisEM(isEMbit,"isEMTight"));
950 ATH_MSG_INFO(" REGTEST: isEMTight bit " << std::hex << isEMbit << std::dec);
951 } else{
952 ATH_MSG_INFO(" REGTEST: problems with egamma pointer" );
953 return StatusCode::SUCCESS;
954 }
955 ATH_MSG_INFO(" REGTEST: caloCluster variables ");
956 if (eg->caloCluster()) {
957 ATH_MSG_INFO(" REGTEST: egamma cluster transverse energy: " << eg->caloCluster()->et() );
958 ATH_MSG_INFO(" REGTEST: egamma cluster eta: " << eg->caloCluster()->eta() );
959 ATH_MSG_INFO(" REGTEST: egamma cluster phi: " << eg->caloCluster()->phi() );
960 double tmpeta = -999.;
961 double tmpphi = -999.;
962 eg->caloCluster()->retrieveMoment(xAOD::CaloCluster::ETACALOFRAME,tmpeta);
963 eg->caloCluster()->retrieveMoment(xAOD::CaloCluster::PHICALOFRAME,tmpphi);
964 ATH_MSG_INFO(" REGTEST: egamma cluster calo-frame coords. etaCalo = " << tmpeta);
965 ATH_MSG_INFO(" REGTEST: egamma cluster calo-frame coords. phiCalo = " << tmpphi);
966 } else{
967 ATH_MSG_INFO(" REGTEST: problems with egamma cluster pointer" );
968 }
969 ATH_MSG_INFO("REGTEST: Check the original (uncalibrated)");
970 if (!orig.isAvailable(*eg->caloCluster()) || !orig(*eg->caloCluster()).isValid()){
971 ATH_MSG_INFO("Problem with original cluster link");
972 }
973 else {
974 const xAOD::CaloCluster *origClus = *orig(*eg->caloCluster());
975 ATH_MSG_INFO("REGTEST:: Compare new and old clusters");
976 ATH_MSG_INFO("REGTEST:: Original Cluster e,eta,phi" << origClus->e() << " " << origClus->eta() << " " << origClus->phi());
977 ATH_MSG_INFO("REGTEST:: MVA Cluster e,eta,phi" << eg->caloCluster()->e() << " " << eg->caloCluster()->eta() << " " << eg->caloCluster()->phi());
978 }
979 //msg() << MSG::VERBOSE << " REGTEST: cluster variables " << endmsg;
980 //clus = eg->caloCluster();
981 ATH_MSG_INFO(" REGTEST: EMShower variables ");
982 eg->showerShapeValue(val_float,xAOD::EgammaParameters::ethad);
983 ATH_MSG_INFO(" REGTEST: ethad = " << val_float);
984 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e011);
985 ATH_MSG_INFO(" REGTEST: e011 = " << val_float);
986 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e132);
987 ATH_MSG_INFO(" REGTEST: e132 = " << val_float);
988 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e237);
989 ATH_MSG_INFO(" REGTEST: e237 = " << val_float);
990 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e335);
991 ATH_MSG_INFO(" REGTEST: e335 = " << val_float);
992 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e2ts1);
993 ATH_MSG_INFO(" REGTEST: e2ts1 = " << val_float);
994 eg->showerShapeValue(val_float,xAOD::EgammaParameters::e2tsts1);
995 ATH_MSG_INFO(" REGTEST: e2tsts1 = " << val_float);
996 eg->isolationValue(val_float,xAOD::Iso::etcone20);
997 ATH_MSG_INFO(" REGTEST: etcone20 = " << val_float);
998 eg->isolationValue(val_float,xAOD::Iso::etcone30);
999 ATH_MSG_INFO(" REGTEST: etcone30 = " << val_float);
1000 eg->isolationValue(val_float,xAOD::Iso::etcone40);
1001 ATH_MSG_INFO(" REGTEST: etcone40 = " << val_float);
1002 //DEBUG info for Electrons which by definition have a track match
1003
1004 }
1005 ATH_MSG_INFO( "REGTEST ==========END of xAOD::PhotonContainer DUMP===========" );
1006
1007 return StatusCode::SUCCESS;
1008
1009}
1010
1012
1013StatusCode TrigEDMChecker::dumpTrigEFBphysContainer() {
1014
1015 ATH_MSG_DEBUG("in dumpTrigEFBphysContainer()");
1016
1017 ATH_MSG_INFO("REGTEST ==========START of TrigEFBphysContainer DUMP===========");
1018
1019 std::string EFBphysTags[]={"HLT_xAOD__TrigBphysContainer_EFBMuMuFex",
1020 "HLT_xAOD__TrigBphysContainer_EFBMuMuXFex",
1021 "HLT_xAOD__TrigBphysContainer_EFDsPhiPiFex",
1022 "HLT_xAOD__TrigBphysContainer_EFMuPairs",
1023 "HLT_xAOD__TrigBphysContainer_EFMultiMuFex",
1024 "HLT_xAOD__TrigBphysContainer_EFTrackMass"
1025 };
1026
1027 int ntag= (int) sizeof(EFBphysTags) / sizeof(EFBphysTags[0]);
1028
1029
1030 for (int itag=0; itag<ntag; itag++){
1031 const xAOD::TrigBphysContainer* trigEFBphys;
1032 StatusCode sc = evtStore()->retrieve(trigEFBphys, EFBphysTags[itag]);
1033 if (sc.isFailure()) {
1034 ATH_MSG_INFO("REGTEST No TrigEFBphysContainer found with tag " << EFBphysTags[itag]);
1035 continue;
1036 }
1037
1038 ATH_MSG_INFO("REGTEST TrigEFBphysContainer found with tag " << EFBphysTags[itag]
1039 << " and size " << trigEFBphys->size());
1040
1041 // for (int i=0; trigEFBphys != lastTrigEFBphys; ++trigEFBphys, ++i) {
1042
1043 //mLog << MSG::INFO << "REGTEST Looking at TrigEFBphysContainer " << i << endmsg;
1044
1045 xAOD::TrigBphysContainer::const_iterator EFBphysItr = trigEFBphys->begin();
1046 xAOD::TrigBphysContainer::const_iterator EFBphysItrE = trigEFBphys->end();
1047
1048 for (int j=0; EFBphysItr != EFBphysItrE; ++EFBphysItr, ++j ) {
1049
1050 ATH_MSG_INFO("REGTEST Looking at TrigEFBphys " << j);
1051
1052 ATH_MSG_INFO("REGTEST TrigEFBphys->eta() returns " << (*EFBphysItr)->eta());
1053 ATH_MSG_INFO("REGTEST TrigEFBphys->phi() returns " << (*EFBphysItr)->phi());
1054 ATH_MSG_INFO("REGTEST TrigEFBphys->mass() returns " << (*EFBphysItr)->mass());
1055 ATH_MSG_INFO("REGTEST TrigEFBphys->fitmass() returns " << (*EFBphysItr)->fitmass());
1056 // ATH_MSG_INFO("REGTEST TrigEFBphys->isValid() returns " << (*EFBphysItr)->isValid());
1057 ATH_MSG_INFO("REGTEST TrigEFBphys->roiId() returns " << (*EFBphysItr)->roiId());
1058 ATH_MSG_INFO("REGTEST TrigEFBphys->particleType() returns " << (*EFBphysItr)->particleType());
1059
1060 if( (*EFBphysItr)->secondaryDecay() != NULL){
1061 const xAOD::TrigBphys * psecond =(*EFBphysItr)->secondaryDecay();
1062 ATH_MSG_INFO("REGTEST Secondary decay info: ");
1063 ATH_MSG_INFO("REGTEST pSecondDecay->eta() returns " << psecond->eta());
1064 ATH_MSG_INFO("REGTEST pSecondDecay->phi() returns " << psecond->phi());
1065 ATH_MSG_INFO("REGTEST pSecondDecay->mass() returns " << psecond->mass());
1066 ATH_MSG_INFO("REGTEST pSecondDecay->fitmass() returns " << psecond->fitmass());
1067 // ATH_MSG_INFO("REGTEST pSecondDecay->isValid() returns " << (*EFBphysItr)->secondaryDecayLink()->isValid());
1068 ATH_MSG_INFO("REGTEST pSecondDecay->roiId() returns " << psecond->roiId());
1069 ATH_MSG_INFO("REGTEST pSecondDecay->particleType() returns " << psecond->particleType());
1070
1071 } // end if secondary exists
1072
1073
1074
1075 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*EFBphysItr)->trackParticleLinks();
1076 if (trackVector.size() != 0) {
1077 ATH_MSG_INFO(" REGTEST got track vector size: " << trackVector.size());
1078 } else {
1079 ATH_MSG_INFO(" REGTEST no track vector!!! " );
1080 }
1081 std::vector<ElementLink<xAOD::TrackParticleContainer> >::const_iterator trkIt=trackVector.begin();
1082 for (int itrk=0 ; trkIt!= trackVector.end(); ++itrk, ++trkIt) {
1083 if (!(trkIt->isValid())) {
1084 ATH_MSG_WARNING("TrackParticleContainer::Invalid ElementLink to track ");
1085 continue;
1086 }
1087 //const Trk::Perigee* trackPerigee=(*(*trkIt))->measuredPerigee();
1088 const Trk::Perigee* trackPerigee=&((*(*trkIt))->perigeeParameters());
1089
1090 // msg() << MSG::VERBOSE << "track, iterator, pointer " << itrk << " " << *trkIt << " " << *(*trkIt) << endmsg;
1091 double phi = trackPerigee->parameters()[Trk::phi];
1092 double theta = trackPerigee->parameters()[Trk::theta];
1093 double px = trackPerigee->momentum()[Trk::px];
1094 double py = trackPerigee->momentum()[Trk::py];
1095 double pt = sqrt(px*px + py*py);
1096 double eta = -std::log(tan(theta/2));
1097
1098 ATH_MSG_INFO("track " << itrk << " pt phi eta " << pt << " " <<
1099 phi << " " << eta);
1100 }
1101
1102 }
1103 }
1104 ATH_MSG_INFO("REGTEST ==========END of TrigEFBphysContainer DUMP===========");
1105 ATH_MSG_DEBUG("dumpTrigEFBphysContainer() succeeded");
1106
1107 return StatusCode::SUCCESS;
1108}
1109
1111StatusCode TrigEDMChecker::dumpTrigL2BphysContainer() {
1112
1113 ATH_MSG_DEBUG("in dumpTrigL2BphysContainer()");
1114 ATH_MSG_INFO("REGTEST ==========START of TrigL2BphysContainer DUMP===========");
1115
1116 std::string L2BphysTags[]={"HLT_xAOD__TrigBphysContainer_L2BMuMuFex",
1117 "HLT_xAOD__TrigBphysContainer_L2BMuMuXFex",
1118 "HLT_xAOD__TrigBphysContainer_L2DsPhiPiFexDs",
1119 "HLT_xAOD__TrigBphysContainer_L2DsPhiPiFexPhi",
1120 "HLT_xAOD__TrigBphysContainer_L2JpsieeFex",
1121 "HLT_xAOD__TrigBphysContainer_L2MultiMuFex",
1122 "HLT_xAOD__TrigBphysContainer_L2TrackMass",
1123 };
1124 const int ntag = (int) sizeof(L2BphysTags) / sizeof(L2BphysTags[0]);
1125
1126
1127 for (int itag=0; itag<ntag; itag++){
1128 const xAOD::TrigBphysContainer* trigL2Bphys;
1129 StatusCode sc = evtStore()->retrieve(trigL2Bphys, L2BphysTags[itag]);
1130 if (sc.isFailure()) {
1131 ATH_MSG_INFO("REGTEST No TrigL2BphysContainer found with tag " << L2BphysTags[itag]);
1132 continue;
1133 }
1134
1135 ATH_MSG_INFO("REGTEST TrigL2BphysContainer found with tag " << L2BphysTags[itag]
1136 << " and size " << trigL2Bphys->size());
1137
1138 // for (int i=0; trigL2Bphys != lastTrigL2Bphys; ++trigL2Bphys, ++i) {
1139
1140 // mLog << MSG::INFO << "REGTEST Looking at TrigL2BphysContainer " << i << endmsg;
1141
1142 xAOD::TrigBphysContainer::const_iterator L2BphysItr = trigL2Bphys->begin();
1143 xAOD::TrigBphysContainer::const_iterator L2BphysItrE = trigL2Bphys->end();
1144
1145 for (int j=0; L2BphysItr != L2BphysItrE; ++L2BphysItr, ++j ) {
1146
1147 ATH_MSG_INFO("REGTEST Looking at TrigL2Bphys " << j);
1148
1149 ATH_MSG_INFO("REGTEST TrigL2Bphys->eta() returns " << (*L2BphysItr)->eta());
1150 ATH_MSG_INFO("REGTEST TrigL2Bphys->phi() returns " << (*L2BphysItr)->phi());
1151 ATH_MSG_INFO("REGTEST TrigL2Bphys->mass() returns " << (*L2BphysItr)->mass());
1152 ATH_MSG_INFO("REGTEST TrigL2Bphys->fitmass() returns " << (*L2BphysItr)->fitmass());
1153 // ATH_MSG_INFO("REGTEST TrigL2Bphys->isValid() returns " << (*L2BphysItr)->isValid());
1154 ATH_MSG_INFO("REGTEST TrigL2Bphys->roiId() returns " << (*L2BphysItr)->roiId());
1155 ATH_MSG_INFO("REGTEST TrigL2Bphys->particleType() returns " << (*L2BphysItr)->particleType());
1156
1157 if( (*L2BphysItr)->secondaryDecay() != NULL){
1158 const xAOD::TrigBphys * psecond =(*L2BphysItr)->secondaryDecay();
1159 ATH_MSG_INFO("REGTEST Secondary decay info: ");
1160 ATH_MSG_INFO("REGTEST pSecondDecay->eta() returns " << psecond->eta());
1161 ATH_MSG_INFO("REGTEST pSecondDecay->phi() returns " << psecond->phi());
1162 ATH_MSG_INFO("REGTEST pSecondDecay->mass() returns " << psecond->mass());
1163 ATH_MSG_INFO("REGTEST pSecondDecay->fitmass() returns " << psecond->fitmass());
1164 // ATH_MSG_INFO("REGTEST pSecondDecay->isValid() returns " << (*L2BphysItr)->secondaryDecayLink()->isValid());
1165 ATH_MSG_INFO("REGTEST pSecondDecay->roiId() returns " << psecond->roiId());
1166 ATH_MSG_INFO("REGTEST pSecondDecay->particleType() returns " << psecond->particleType());
1167 } // end if secondary exists
1168
1169 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*L2BphysItr)->trackParticleLinks();
1170 if (trackVector.size() != 0) {
1171 ATH_MSG_INFO(" REGTEST got track vector size: " << trackVector.size());
1172 } else {
1173 ATH_MSG_INFO(" REGTEST no track vector!!! " );
1174 }
1175 std::vector<ElementLink<xAOD::TrackParticleContainer> >::const_iterator trkIt=trackVector.begin();
1176 for (int itrk=0 ; trkIt!= trackVector.end(); ++itrk, ++trkIt) {
1177 if (!(trkIt->isValid())) {
1178 ATH_MSG_WARNING("TrackParticleContainer::Invalid ElementLink to track ");
1179 continue;
1180 }
1181 //const Trk::Perigee* trackPerigee=(*(*trkIt))->measuredPerigee();
1182 const Trk::Perigee* trackPerigee=&((*(*trkIt))->perigeeParameters());
1183
1184 // msg() << MSG::VERBOSE << "track, iterator, pointer " << itrk << " " << *trkIt << " " << *(*trkIt) << endmsg;
1185 double phi = trackPerigee->parameters()[Trk::phi];
1186 double theta = trackPerigee->parameters()[Trk::theta];
1187 double px = trackPerigee->momentum()[Trk::px];
1188 double py = trackPerigee->momentum()[Trk::py];
1189 double pt = sqrt(px*px + py*py);
1190 double eta = -std::log(tan(theta/2));
1191
1192 ATH_MSG_INFO("track " << itrk << " pt phi eta " << pt << " " <<
1193 phi << " " << eta);
1194 }
1195 }
1196 }
1197
1198 ATH_MSG_INFO("REGTEST ==========END of TrigL2BphysContainer DUMP===========");
1199 ATH_MSG_DEBUG("dumpTrigL2BphysContainer() succeeded");
1200
1201 return StatusCode::SUCCESS;
1202}
1203
1204
1205
1207
1208StatusCode TrigEDMChecker::dumpxAODJetContainer() {
1209
1210 ATH_MSG_DEBUG("in dumpxAODJetContainer()");
1211
1212 ATH_MSG_INFO("REGTEST ==========START of xAOD::JetContainer DUMP===========");
1213
1214 std::string containerName[30] = {"HLT_xAOD__JetContainer_a4tcemjesFS",
1215 "HLT_xAOD__JetContainer_a4tcemsubFS",
1216 "HLT_xAOD__JetContainer_a4tcemsubjesFS",
1217 "HLT_xAOD__JetContainer_a4tcemnojcalibFS",
1218 "HLT_xAOD__JetContainer_a4tcemjesPS",
1219 "HLT_xAOD__JetContainer_a4tcemnojcalibPS",
1220 "HLT_xAOD__JetContainer_a4tclcwjesFS",
1221 "HLT_xAOD__JetContainer_a4tclcwsubFS",
1222 "HLT_xAOD__JetContainer_a4tclcwsubjesFS",
1223 "HLT_xAOD__JetContainer_a4tclcwnojcalibFS",
1224 "HLT_xAOD__JetContainer_a4tclcwjesPS",
1225 "HLT_xAOD__JetContainer_a4tclcwnojcalibPS",
1226 "HLT_xAOD__JetContainer_a4TTemnojcalibFS",
1227 "HLT_xAOD__JetContainer_a4TThadnojcalibFS",
1228 "HLT_xAOD__JetContainer_a10tcemjesFS",
1229 "HLT_xAOD__JetContainer_a10tcemsubFS",
1230 "HLT_xAOD__JetContainer_a10tcemsubjesFS",
1231 "HLT_xAOD__JetContainer_a10tcemnojcalibFS",
1232 "HLT_xAOD__JetContainer_a10tcemjesPS",
1233 "HLT_xAOD__JetContainer_a10tcemnojcalibPS",
1234 "HLT_xAOD__JetContainer_a10tclcwjesFS",
1235 "HLT_xAOD__JetContainer_a10tclcwsubFS",
1236 "HLT_xAOD__JetContainer_a10tclcwsubjesFS",
1237 "HLT_xAOD__JetContainer_a10tclcwnojcalibFS",
1238 "HLT_xAOD__JetContainer_a10tclcwjesPS",
1239 "HLT_xAOD__JetContainer_a10tclcwnojcalibPS",
1240 "HLT_xAOD__JetContainer_a10TTemnojcalibFS",
1241 "HLT_xAOD__JetContainer_a10TThadnojcalibFS",
1242 "HLT_xAOD__JetContainer_a10r_tcemsubjesFS",
1243 "HLT_xAOD__JetContainer_TrigHLTJetDSSelectorCollection"};
1244 float containerSizeParameter[30] = {0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4,
1245 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 0.4};
1246 int containerInputType[30] = {1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 2,
1247 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 2, 1, 1};
1248 int containerSignalState[30] = {0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1,
1249 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0};
1250 bool onefilled = false;
1251
1252 for (int icont=0; icont<30; ++icont) {
1253 //retrieve jet collection
1254 const xAOD::JetContainer* jetCont = 0;
1255 StatusCode sc = evtStore()->retrieve(jetCont, containerName[icont]) ;
1256
1257 if (sc.isFailure()) {
1258 ATH_MSG_WARNING("REGTEST Cannot retrieve jet container");
1259 continue;
1260 }
1261
1262 int jetContsize = jetCont->size();
1263 ATH_MSG_INFO("REGTEST Got jet container " << containerName[icont] << ", size: " << jetContsize);
1264 if (jetContsize != 0) {
1265 onefilled = true;
1266 int i = 0;
1267 for(const auto thisjet : *jetCont) {
1268 ++i;
1269 ATH_MSG_INFO( "REGTEST Looking at jet " << i);
1270 if (thisjet) {
1271 //checks jet variables
1272 ATH_MSG_DEBUG("REGTEST Checking jet variables");
1273 ATH_MSG_INFO( "REGTEST pt: " << thisjet->pt() );
1274 ATH_MSG_INFO( "REGTEST eta: " << thisjet->eta() );
1275 ATH_MSG_INFO( "REGTEST phi: " << thisjet->phi() );
1276 ATH_MSG_INFO( "REGTEST m: " << thisjet->m() );
1277 ATH_MSG_INFO( "REGTEST e: " << thisjet->e() );
1278 ATH_MSG_INFO( "REGTEST rapidity: " << thisjet->rapidity() );
1279 ATH_MSG_INFO( "REGTEST px: " << thisjet->px() );
1280 ATH_MSG_INFO( "REGTEST py: " << thisjet->py() );
1281 ATH_MSG_INFO( "REGTEST pz: " << thisjet->pz() );
1282 ATH_MSG_INFO( "REGTEST type: " << thisjet->type() );
1283 ATH_MSG_INFO( "REGTEST algorithm (kt: 0, cam: 1, antikt: 2, ...): " << thisjet->getAlgorithmType() << "; should be 2");
1284 if(thisjet->getAlgorithmType() != 2) ATH_MSG_WARNING("Jet algorithm different from container");
1285 ATH_MSG_INFO( "REGTEST size parameter: " << thisjet->getSizeParameter() << "; should be " << containerSizeParameter[icont]);
1286 if(thisjet->getSizeParameter() != containerSizeParameter[icont]) ATH_MSG_WARNING("Jet size different from container");
1287 ATH_MSG_INFO( "REGTEST input (LCTopo: 0, EMTopo: 1, TopoTower: 2, ...): " << thisjet->getInputType() << "; should be " << containerInputType[icont]);
1288 if(thisjet->getInputType() != containerInputType[icont]) ATH_MSG_WARNING("Jet input different from container");
1289 ATH_MSG_INFO( "REGTEST constituents signal state (uncalibrated: 0, calibrated: 1): " << thisjet->getConstituentsSignalState() << "; should be " << containerSignalState[icont]);
1290 if(thisjet->getConstituentsSignalState() != containerSignalState[icont]) ATH_MSG_WARNING("Jet constituents' signal state different from container");
1291 ATH_MSG_INFO( "REGTEST number of constituents: " << thisjet->numConstituents() );
1292
1293
1294 //checks the constituents
1295 ATH_MSG_DEBUG("REGTEST Checking jet constituents");
1296
1297 xAOD::JetConstituentVector constitCont = thisjet->getConstituents();
1298 unsigned int constitContsize = constitCont.size();
1299 ATH_MSG_INFO("REGTEST Got constituent vector, size: " << constitContsize << "; should be " << thisjet->numConstituents());
1300 if(constitContsize != thisjet->numConstituents()) ATH_MSG_WARNING("Constituents container size different from number of constituents");
1301
1302 if (constitContsize != 0) {
1303// int j = 0;
1304 // for (const auto thisconstit : constitCont) {
1305 // ++j;
1306 // ATH_MSG_INFO( "REGTEST Looking at constituent " << j);
1307 // if(thisconstit){
1308 // ATH_MSG_INFO( "REGTEST constituent pt: " << thisconstit->pt() );
1309 // ATH_MSG_INFO( "REGTEST constituent eta: " << thisconstit->eta() );
1310 // ATH_MSG_INFO( "REGTEST constituent phi: " << thisconstit->phi() );
1311 // ATH_MSG_INFO( "REGTEST constituent m: " << thisconstit->m() );
1312 // ATH_MSG_INFO( "REGTEST constituent e: " << thisconstit->e() );
1313 // ATH_MSG_INFO( "REGTEST constituent type (CaloCluster: 1, Jet: 2, ...): " << thisconstit->type() );
1314 // }
1315 // else{
1316 // ATH_MSG_WARNING("REGTEST Problem with constituent pointer");
1317 // return StatusCode::SUCCESS;
1318 // }
1319 // }
1320 // ATH_MSG_INFO("REGTEST size of constituent vector == number of displayed constituents: "<< (constitContsize == j) );
1321 // if (constitContsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1322 }
1323 // else ATH_MSG_WARNING("REGTEST This jet has 0 constituents!");
1324
1325 //checks int attributes
1326 ATH_MSG_DEBUG("REGTEST Checking int attributes");
1327 int valueint;
1328 if( thisjet->getAttribute(xAOD::JetAttribute::GhostMuonSegmentCount, valueint)) ATH_MSG_INFO("REGTEST GhostMuonSegmentCount: " << valueint);
1329 if( thisjet->getAttribute(xAOD::JetAttribute::GhostTrackCount, valueint)) ATH_MSG_INFO("REGTEST GhostTrackCount: " << valueint);
1330 if( thisjet->getAttribute(xAOD::JetAttribute::GhostTruthParticleCount, valueint)) ATH_MSG_INFO("REGTEST GhostTruthParticleCount: " << valueint);
1331 if( thisjet->getAttribute(xAOD::JetAttribute::FracSamplingMaxIndex, valueint)) ATH_MSG_INFO( "REGTEST FracSamplingMaxIndex: " << valueint);
1332
1333
1334 //checks float attributes
1335 ATH_MSG_DEBUG("REGTEST Checking float attributes");
1336 float value;
1337 if( thisjet->getAttribute(xAOD::JetAttribute::ActiveArea, value)) ATH_MSG_INFO( "REGTEST ActiveArea: " << value );
1338 // Four-vector type if( thisjet->getAttribute(xAOD::JetAttribute::ActiveArea4vec, value)) ATH_MSG_INFO( "REGTEST ActiveArea4vec: " << value );
1339 if( thisjet->getAttribute(xAOD::JetAttribute::AverageLArQF, value)) ATH_MSG_INFO( "REGTEST AverageLArQF: " << value );
1340 if( thisjet->getAttribute(xAOD::JetAttribute::BchCorrCell, value)) ATH_MSG_INFO( "REGTEST BchCorrCell: " << value );
1341 if( thisjet->getAttribute(xAOD::JetAttribute::BchCorrDotx, value)) ATH_MSG_INFO( "REGTEST BchCorrDotx: " << value );
1342 if( thisjet->getAttribute(xAOD::JetAttribute::BchCorrJet, value)) ATH_MSG_INFO( "REGTEST BchCorrJet: " << value );
1343 if( thisjet->getAttribute(xAOD::JetAttribute::BchCorrJetForCell, value)) ATH_MSG_INFO( "REGTEST BchCorrJetForCell: " << value );
1344 if( thisjet->getAttribute(xAOD::JetAttribute::CentroidR, value)) ATH_MSG_INFO( "REGTEST CentroidR: " << value );
1345 if( thisjet->getAttribute(xAOD::JetAttribute::HECQuality, value)) ATH_MSG_INFO( "REGTEST HECQuality: " << value );
1346 if( thisjet->getAttribute(xAOD::JetAttribute::IsoKR20Par, value)) ATH_MSG_INFO( "REGTEST IsoKR20Par: " << value );
1347 if( thisjet->getAttribute(xAOD::JetAttribute::IsoKR20Perp, value)) ATH_MSG_INFO( "REGTEST IsoKR20Perp: " << value );
1348 // ElementLink<DataVector<xAOD::Vertex> > type if( thisjet->getAttribute(xAOD::JetAttribute::HighestJVFVtx, value)) ATH_MSG_INFO( "REGTEST HighestJVFVtx: " << value );
1349 // ??? type if( thisjet->getAttribute(xAOD::JetAttribute::JetLabel, value)) ATH_MSG_INFO( "REGTEST JetLabel: " << value );
1350 if( thisjet->getAttribute(xAOD::JetAttribute::KtDR, value)) ATH_MSG_INFO( "REGTEST KtDR: " << value );
1351 if( thisjet->getAttribute(xAOD::JetAttribute::LArBadHVEnergy, value)) ATH_MSG_INFO( "REGTEST LArBadHVEnergy: " << value );
1352 if( thisjet->getAttribute(xAOD::JetAttribute::LArBadHVRatio, value)) ATH_MSG_INFO( "REGTEST LArBadHVRatio: " << value );
1353 if( thisjet->getAttribute(xAOD::JetAttribute::LArQuality, value)) ATH_MSG_INFO( "REGTEST LArQuality: " << value );
1354 if( thisjet->getAttribute(xAOD::JetAttribute::NegativeE, value)) ATH_MSG_INFO( "REGTEST NegativeE: " << value );
1355 // no tools available yet if( thisjet->getAttribute(xAOD::JetAttribute::NumTowers, value)) ATH_MSG_INFO( "REGTEST NumTowers: " << value );
1356 if( thisjet->getAttribute(xAOD::JetAttribute::SamplingMax, value)) ATH_MSG_INFO( "REGTEST SamplingMax: " << value );
1357 if( thisjet->getAttribute(xAOD::JetAttribute::Timing, value)) ATH_MSG_INFO( "REGTEST Timing: " << value );
1358 if( thisjet->getAttribute(xAOD::JetAttribute::VoronoiArea, value)) ATH_MSG_INFO( "REGTEST VoronoiArea: " << value );
1359 if( thisjet->getAttribute(xAOD::JetAttribute::VoronoiAreaE, value)) ATH_MSG_INFO( "REGTEST VoronoiAreaE: " << value );
1360 if( thisjet->getAttribute(xAOD::JetAttribute::VoronoiAreaPx, value)) ATH_MSG_INFO( "REGTEST VoronoiAreaPx: " << value );
1361 if( thisjet->getAttribute(xAOD::JetAttribute::VoronoiAreaPy, value)) ATH_MSG_INFO( "REGTEST VoronoiAreaPy: " << value );
1362 if( thisjet->getAttribute(xAOD::JetAttribute::VoronoiAreaPz, value)) ATH_MSG_INFO( "REGTEST VoronoiAreaPz: " << value );
1363 if( thisjet->getAttribute(xAOD::JetAttribute::Width, value)) ATH_MSG_INFO( "REGTEST WIDTH: " << value );
1364 if( thisjet->getAttribute(xAOD::JetAttribute::FracSamplingMax, value)) ATH_MSG_INFO( "REGTEST FracSamplingMax: " << value );
1365 if( thisjet->getAttribute(xAOD::JetAttribute::EMFrac, value)) ATH_MSG_INFO( "REGTEST EMFrac: " << value );
1366 if( thisjet->getAttribute(xAOD::JetAttribute::HECFrac, value)) ATH_MSG_INFO( "REGTEST HECFrac: " << value );
1367 if( thisjet->getAttribute(xAOD::JetAttribute::isBadLoose, value)) ATH_MSG_INFO( "REGTEST isBadLoose: " << value );
1368 if( thisjet->getAttribute(xAOD::JetAttribute::isBadMedium, value)) ATH_MSG_INFO( "REGTEST isBadMedium: " << value );
1369 if( thisjet->getAttribute(xAOD::JetAttribute::isBadTight, value)) ATH_MSG_INFO( "REGTEST isBadTight: " << value );
1370 // unknown attribute if( thisjet->getAttribute(xAOD::JetAttribute::isUgly, value)) ATH_MSG_INFO( "REGTEST isUgly: " << value );
1371 if( thisjet->getAttribute(xAOD::JetAttribute::N90Constituents, value)) ATH_MSG_INFO( "REGTEST N90Constituents: " << value );
1372 if( thisjet->getAttribute(xAOD::JetAttribute::N90Cells, value)) ATH_MSG_INFO( "REGTEST N90Cells: " << value );
1373 if( thisjet->getAttribute(xAOD::JetAttribute::OotFracClusters10, value)) ATH_MSG_INFO( "REGTEST OotFracClusters10: " << value );
1374 if( thisjet->getAttribute(xAOD::JetAttribute::OotFracClusters5, value)) ATH_MSG_INFO( "REGTEST OotFracClusters5: " << value );
1375 if( thisjet->getAttribute(xAOD::JetAttribute::OotFracCells5, value)) ATH_MSG_INFO( "REGTEST OotFracCells5: " << value );
1376 if( thisjet->getAttribute(xAOD::JetAttribute::OotFracCells10, value)) ATH_MSG_INFO( "REGTEST OotFracCells10: " << value );
1377 if( thisjet->getAttribute(xAOD::JetAttribute::PtTruth, value)) ATH_MSG_INFO( "REGTEST PtTruth: " << value );
1378 if( thisjet->getAttribute(xAOD::JetAttribute::Tau1, value)) ATH_MSG_INFO( "REGTEST Tau1: " << value );
1379 if( thisjet->getAttribute(xAOD::JetAttribute::Tau2, value)) ATH_MSG_INFO( "REGTEST Tau2: " << value );
1380 if( thisjet->getAttribute(xAOD::JetAttribute::Tau3, value)) ATH_MSG_INFO( "REGTEST Tau3: " << value );
1381 // unknown attribute if( thisjet->getAttribute(xAOD::JetAttribute::Split12, value)) ATH_MSG_INFO( "REGTEST Split12: " << value );
1382 // unknown attribute if( thisjet->getAttribute(xAOD::JetAttribute::Split23, value)) ATH_MSG_INFO( "REGTEST Split23: " << value );
1383 // unknown attribute if( thisjet->getAttribute(xAOD::JetAttribute::Split34, value)) ATH_MSG_INFO( "REGTEST Split34: " << value );
1384 if( thisjet->getAttribute(xAOD::JetAttribute::Dip12, value)) ATH_MSG_INFO( "REGTEST Dip12: " << value );
1385 if( thisjet->getAttribute(xAOD::JetAttribute::Dip13, value)) ATH_MSG_INFO( "REGTEST Dip13: " << value );
1386 if( thisjet->getAttribute(xAOD::JetAttribute::Dip23, value)) ATH_MSG_INFO( "REGTEST Dip23: " << value );
1387 if( thisjet->getAttribute(xAOD::JetAttribute::DipExcl12, value)) ATH_MSG_INFO( "REGTEST DipExcl12: " << value );
1388 if( thisjet->getAttribute(xAOD::JetAttribute::ThrustMin, value)) ATH_MSG_INFO( "REGTEST ThrustMin: " << value );
1389 if( thisjet->getAttribute(xAOD::JetAttribute::ThrustMaj, value)) ATH_MSG_INFO( "REGTEST ThrustMaj: " << value );
1390 if( thisjet->getAttribute(xAOD::JetAttribute::FoxWolfram0, value)) ATH_MSG_INFO( "REGTEST FoxWolfram0: " << value );
1391 if( thisjet->getAttribute(xAOD::JetAttribute::FoxWolfram1, value)) ATH_MSG_INFO( "REGTEST FoxWolfram1: " << value );
1392 if( thisjet->getAttribute(xAOD::JetAttribute::FoxWolfram2, value)) ATH_MSG_INFO( "REGTEST FoxWolfram2: " << value );
1393 if( thisjet->getAttribute(xAOD::JetAttribute::FoxWolfram3, value)) ATH_MSG_INFO( "REGTEST FoxWolfram3: " << value );
1394 if( thisjet->getAttribute(xAOD::JetAttribute::FoxWolfram4, value)) ATH_MSG_INFO( "REGTEST FoxWolfram4: " << value );
1395 if( thisjet->getAttribute(xAOD::JetAttribute::Sphericity, value)) ATH_MSG_INFO( "REGTEST Sphericity: " << value );
1396 if( thisjet->getAttribute(xAOD::JetAttribute::Aplanarity, value)) ATH_MSG_INFO( "REGTEST Aplanarity: " << value );
1397
1398 //checks vector<int> attributes
1399 ATH_MSG_DEBUG("REGTEST Checking vector<int> attributes");
1400 std::vector<int> vecvalueint;
1401 if (thisjet->getAttribute(xAOD::JetAttribute::NumTrkPt1000, vecvalueint)) {
1402 int vecsize = vecvalueint.size();
1403 ATH_MSG_INFO("REGTEST Got NumTrkPt1000 vector, size: " << vecsize);
1404 if (vecsize != 0) {
1405 int j = 0;
1406 for(const auto & thisvalue : vecvalueint){
1407 ++j;
1408 ATH_MSG_INFO("REGTEST NumTrkPt1000 #" << j << ": " << thisvalue);
1409 }
1410 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1411 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1412 }
1413 }
1414 if (thisjet->getAttribute(xAOD::JetAttribute::NumTrkPt500, vecvalueint)) {
1415 int vecsize = vecvalueint.size();
1416 ATH_MSG_INFO("REGTEST Got Got NumTrkPt500 vector, size: " << vecsize);
1417 if (vecsize != 0) {
1418 int j = 0;
1419 for(const auto & thisvalue : vecvalueint){
1420 ++j;
1421 ATH_MSG_INFO("REGTEST NumTrkPt500 #" << j << ": " << thisvalue);
1422 }
1423 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1424 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1425 }
1426 }
1427
1428 //checks vector<float> attributes
1429 ATH_MSG_DEBUG("REGTEST Checking vector<float> attributes");
1430 std::vector<float> vecvalue;
1431 if (thisjet->getAttribute(xAOD::JetAttribute::JVF, vecvalue)) {
1432 int vecsize = vecvalue.size();
1433 ATH_MSG_INFO("REGTEST Got JVF vector, size: " << vecsize);
1434 if (vecsize != 0) {
1435 int j = 0;
1436 for(const auto & thisvalue : vecvalue){
1437 ++j;
1438 ATH_MSG_INFO("REGTEST JVF #" << j << ": " << thisvalue);
1439 }
1440 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1441 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1442 }
1443 }
1444 if (thisjet->getAttribute(xAOD::JetAttribute::SumPtTrkPt1000, vecvalue)) {
1445 int vecsize = vecvalue.size();
1446 ATH_MSG_INFO("REGTEST Got SumPtTrkPt1000 vector, size: " << vecsize);
1447 if (vecsize != 0) {
1448 int j = 0;
1449 for(const auto & thisvalue : vecvalue){
1450 ++j;
1451 ATH_MSG_INFO("REGTEST SumPtTrkPt1000 #" << j << ": " << thisvalue);
1452 }
1453 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1454 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1455 }
1456 }
1457 if (thisjet->getAttribute(xAOD::JetAttribute::SumPtTrkPt500, vecvalue)) {
1458 int vecsize = vecvalue.size();
1459 ATH_MSG_INFO("REGTEST Got SumPtTrkPt500 vector, size: " << vecsize);
1460 if (vecsize != 0) {
1461 int j = 0;
1462 for(const auto & thisvalue : vecvalue){
1463 ++j;
1464 ATH_MSG_INFO("REGTEST SumPtTrkPt500 #" << j << ": " << thisvalue);
1465 }
1466 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1467 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1468 }
1469 }
1470 if (thisjet->getAttribute(xAOD::JetAttribute::TrackWidthPt1000, vecvalue)) {
1471 int vecsize = vecvalue.size();
1472 ATH_MSG_INFO("REGTEST Got TrackWidthPt1000 vector, size: " << vecsize);
1473 if (vecsize != 0) {
1474 int j = 0;
1475 for(const auto & thisvalue : vecvalue){
1476 ++j;
1477 ATH_MSG_INFO("REGTEST TrackWidthPt1000 #" << j << ": " << thisvalue);
1478 }
1479 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1480 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1481 }
1482 }
1483 if (thisjet->getAttribute(xAOD::JetAttribute::TrackWidthPt500, vecvalue)) {
1484 int vecsize = vecvalue.size();
1485 ATH_MSG_INFO("REGTEST Got TrackWidthPt500 vector, size: " << vecsize);
1486 if (vecsize != 0) {
1487 int j = 0;
1488 for(const auto & thisvalue : vecvalue){
1489 ++j;
1490 ATH_MSG_INFO("REGTEST TrackWidthPt500 #" << j << ": " << thisvalue);
1491 }
1492 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1493 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1494 }
1495 }
1496 if (thisjet->getAttribute(xAOD::JetAttribute::EnergyPerSampling, vecvalue)) {
1497 int vecsize = vecvalue.size();
1498 ATH_MSG_INFO("REGTEST Got EnergyPerSampling vector, size: " << vecsize);
1499 if (vecsize != 0) {
1500 int j = 0;
1501 for(const auto & thisvalue : vecvalue){
1502 ++j;
1503 ATH_MSG_INFO("REGTEST EnergyPerSampling #" << j << ": " << thisvalue);
1504 }
1505 ATH_MSG_INFO("REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
1506 if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this attribute");
1507 }
1508 }
1509
1510
1511 //checks associated objects
1512 // ATH_MSG_DEBUG("REGTEST Checking associated objects");
1513 // std::vector<const xAOD::TrackParticle*> track;
1514 // if( thisjet->getAssociatedObjects(xAOD::JetAttribute::GhostTrack, track) ){
1515 // int vecsize = track.size();
1516 // ATH_MSG_INFO("REGTEST Got GhostTrack vector, size: " << vecsize);
1517 // if (vecsize != 0) {
1518 // int j = 0;
1519 // for(const auto & thistrack : track){
1520 // ++j;
1521 // //checks only one associated variable, just making sure getting the object worked
1522 // if (thistrack) ATH_MSG_INFO("REGTEST z0 for GhostTrack #" << j << ": " << thistrack->z0());
1523 // else{
1524 // ATH_MSG_WARNING("REGTEST Problem with attribute pointer");
1525 // return StatusCode::SUCCESS;
1526 // }
1527 // }
1528 // ATH_MSG_INFO("REGTEST size of associated object vector == number of displayed attributes: " << (vecsize == j) );
1529 // if (vecsize != j) ATH_MSG_WARNING("REGTEST Problem with displaying this associated object");
1530 // }
1531 // }
1532 }
1533 else{
1534 ATH_MSG_WARNING("REGTEST Problem with jet pointer");
1535 return StatusCode::SUCCESS;
1536 }
1537 }
1538
1539 if (jetContsize == i) ATH_MSG_INFO("REGTEST size of jet container == number of displayed jets: " << (jetContsize == i) );
1540 else ATH_MSG_WARNING("REGTEST Problem with displaying jets");
1541 }
1542 }
1543
1544 if (!onefilled) ATH_MSG_DEBUG("There was no filled jet containers");
1545
1546 ATH_MSG_INFO("REGTEST ==========END of xAOD::JetContainer DUMP===========");
1547
1548 ATH_MSG_DEBUG("leaving dumpxAODJetContainer()");
1549
1550 return StatusCode::SUCCESS;
1551}
1552
1554
1555StatusCode TrigEDMChecker::dumpxAODTrigEMCluster() {
1556
1557 ATH_MSG_DEBUG("in dumpxAODTrigEMCluster()");
1558
1559 ATH_MSG_INFO("REGTEST ==========START of TrigEMCluster DUMP===========");
1560
1561 SG::ConstIterator< xAOD::TrigEMCluster > EMCluster;
1562 SG::ConstIterator< xAOD::TrigEMCluster > lastEMCluster;
1563
1564 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
1565 if (sc.isFailure()) {
1566 ATH_MSG_INFO("REGTEST No xAOD::TrigEMCluster found");
1567 return StatusCode::SUCCESS;
1568 }
1569 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster retrieved");
1570
1571
1572 for (int i=0; EMCluster != lastEMCluster; ++EMCluster, ++i) {
1573
1574 const xAOD::TrigEMCluster* thisEMCluster = &(*EMCluster);
1575
1576 ATH_MSG_INFO("REGTEST Looking at xAOD::TrigEMCluster " << i);
1577
1578 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->energy() returns " << thisEMCluster->energy());
1579 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->e() returns " << thisEMCluster->energy());
1580 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->phi() returns " << thisEMCluster->phi());
1581
1582 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->eta() returns " << thisEMCluster->eta());
1583 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->e237() returns " << thisEMCluster->e237());
1584 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->e277() returns " << thisEMCluster->e277());
1585 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->fracs1() returns " << thisEMCluster->fracs1());
1586 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->weta2() returns " << thisEMCluster->weta2());
1587 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->ehad1() returns " << thisEMCluster->ehad1());
1588 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->Eta1() returns " << thisEMCluster->eta1());
1589 }
1590
1591 ATH_MSG_INFO("REGTEST ==========END of xAOD::TrigEMCluster DUMP===========");
1592 ATH_MSG_DEBUG("dumpxAODTrigEMCluster() succeeded");
1593
1594 return StatusCode::SUCCESS;
1595}
1596
1598
1599StatusCode TrigEDMChecker::dumpxAODTrigEMClusterContainer() {
1600
1601 ATH_MSG_DEBUG("in dumpxAODTrigEMClusterContainer()");
1602
1603 ATH_MSG_INFO("REGTEST ==========START of xAODTrigEMClusterContainer DUMP===========");
1604
1605 SG::ConstIterator< xAOD::TrigEMClusterContainer > EMCluster;
1606 SG::ConstIterator< xAOD::TrigEMClusterContainer > lastEMCluster;
1607
1608 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
1609 if (sc.isFailure()) {
1610 ATH_MSG_INFO("REGTEST No xAOD::TrigEMClusterContainer found");
1611 return StatusCode::SUCCESS;
1612 }
1613 ATH_MSG_INFO("REGTEST xAOD::TrigEMClusterContainer retrieved");
1614
1615
1616 for (int i=0; EMCluster != lastEMCluster; ++EMCluster, ++i) {
1617
1618 ATH_MSG_INFO("REGTEST Looking at xAOD::TrigEMClusterContainer " << i);
1619
1620 xAOD::TrigEMClusterContainer::const_iterator EMClusterItr = EMCluster->begin();
1621 xAOD::TrigEMClusterContainer::const_iterator EMClusterItrE = EMCluster->end();
1622
1623 for (int j=0; EMClusterItr != EMClusterItrE; ++EMClusterItr, ++j ) {
1624
1625 ATH_MSG_INFO("REGTEST Looking at xAOD::TrigEMCluster " << j);
1626 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->energy() returns " << (*EMClusterItr)->energy());
1627 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->et() returns " << (*EMClusterItr)->et());
1628 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->eta() returns " << (*EMClusterItr)->eta());
1629 ATH_MSG_INFO("REGTEST xAOD::TrigEMCluster->phi() returns " << (*EMClusterItr)->phi());
1630 //msg() <<MSG::INFO << "REGTEST xAOD::TrigEMCluster->print() gives" << endmsg;
1631 //int level = msg().level();
1632 // little trick to print out stuff
1633 //msg().setLevel(MSG::DEBUG);
1634 //(*EMClusterItr)->print(msg());
1635 //msg().setLevel(level);
1636 }
1637 }
1638
1639 ATH_MSG_INFO("REGTEST ==========END of TrigEMClusterContainer DUMP===========");
1640 ATH_MSG_DEBUG("dumpTrigEMClusterContainer() succeeded");
1641
1642 return StatusCode::SUCCESS;
1643}
1644
1646StatusCode TrigEDMChecker::dumpxAODTauJetContainer() {
1647
1648 ATH_MSG_DEBUG("In dumpxAODTauJetContainer");
1649 ATH_MSG_INFO( "REGTEST ==========START of xAOD::TauJetContainer DUMP===========" );
1650 const xAOD::TauJetContainer * TauJetcont = 0;
1651 StatusCode sc = evtStore() -> retrieve (TauJetcont, "HLT_xAOD__TauJetContainer_TrigTauRecMerged");
1652
1653 if (sc.isFailure()) {
1654
1655 ATH_MSG_INFO("REGTEST No Tau container HLT_xAOD__TauJetContainer_TrigTauRecMerged");
1656
1657
1658 return StatusCode::SUCCESS;
1659 }
1660
1661
1662 for(xAOD::TauJetContainer::const_iterator tauIt = TauJetcont->begin(); tauIt != TauJetcont->end();++tauIt){
1663
1664 ATH_MSG_INFO( "REGTEST (*tauIt)->eta() returns " << (*tauIt)->eta() );
1665 ATH_MSG_INFO( "REGTEST (*tauIt)->phi() returns " << (*tauIt)->phi() );
1666 ATH_MSG_INFO( "REGTEST (*tauIt)->pt() returns " << (*tauIt)->pt() );
1667
1668 // for numTracks()
1669 int EFnTracks = -1;
1670 #ifndef XAODTAU_VERSIONS_TAUJET_V3_H
1671 EFnTracks = (*tauIt)->nTracks();
1672 #else
1673 (*tauIt)->detail(xAOD::TauJetParameters::nChargedTracks, EFnTracks);
1674 #endif
1675
1676 ATH_MSG_INFO( "REGTEST (*tauIt)->nTracks() returns " << EFnTracks );
1677
1678 // for nTracksIsolation()
1679 int EFWidenTrack = -1;
1680 #ifndef XAODTAU_VERSIONS_TAUJET_V3_H
1681 EFWidenTrack = (*tauIt)->nWideTracks();
1682 #else
1683 (*tauIt)->detail(xAOD::TauJetParameters::nIsolatedTracks, EFWidenTrack);
1684 #endif
1685
1686 ATH_MSG_INFO( "REGTEST (*tauIt)->nWideTracks() returns " << EFWidenTrack );
1687
1688 //bool test = false;
1689 float trkAvgDist=0;
1690 float etOvPtLead=0;
1691 float emRadius=0;
1692 float hadRadius=0;
1693 float IsoFrac=0;
1694 float centFrac=0;
1695 float ipSigLeadTrk=0;
1696 float trFlightPathSig=0;
1697 float dRmax=0;
1698 float massTrkSys=0;
1699 float PSSFraction=0;
1700 float EMPOverTrkSysP=0;
1701 float ChPiEMEOverCaloEME=0;
1702 float EtEm=0;
1703 float EtHad=0;
1704
1705 if ( (*tauIt)->detail(xAOD::TauJetParameters::trkAvgDist,trkAvgDist))
1706 ATH_MSG_INFO( "REGTEST TauDetails->trkAvgDist() returns " << trkAvgDist);
1707
1708 if ( (*tauIt)->detail(xAOD::TauJetParameters::etOverPtLeadTrk,etOvPtLead))
1709 ATH_MSG_INFO( "REGTEST TauDetails->etOverPtLeadTrk() returns " << etOvPtLead);
1710
1711 if ( (*tauIt)->detail(xAOD::TauJetParameters::EMRadius,emRadius))
1712 ATH_MSG_INFO( "REGTEST TauDetails->EMRadius() returns " << emRadius);
1713
1714 if ( (*tauIt)->detail(xAOD::TauJetParameters::hadRadius,hadRadius))
1715 ATH_MSG_INFO( "REGTEST TauDetails->hadRadius() returns " << hadRadius);
1716
1717 if ( (*tauIt)->detail(xAOD::TauJetParameters::isolFrac,IsoFrac))
1718 ATH_MSG_INFO( "REGTEST TauDetails->isolFrac() returns " << IsoFrac);
1719
1720 if ( (*tauIt)->detail(xAOD::TauJetParameters::centFrac,centFrac))
1721 ATH_MSG_INFO( "REGTEST TauDetails->centFrac() returns " << centFrac);
1722
1723 if ( (*tauIt)->detail(xAOD::TauJetParameters::ipSigLeadTrk,ipSigLeadTrk))
1724 ATH_MSG_INFO( "REGTEST TauDetails->ipSigLeadTrk() returns " << ipSigLeadTrk);
1725
1726 if ( (*tauIt)->detail(xAOD::TauJetParameters::trFlightPathSig,trFlightPathSig))
1727 ATH_MSG_INFO( "REGTEST TauDetails->trFlightPathSig() returns " << trFlightPathSig);
1728
1729 if ( (*tauIt)->detail(xAOD::TauJetParameters::dRmax,dRmax))
1730 ATH_MSG_INFO( "REGTEST TauDetails->dRmax() returns " << dRmax);
1731
1732 if ( (*tauIt)->detail(xAOD::TauJetParameters::massTrkSys,massTrkSys)){
1733 massTrkSys /=1000;
1734 ATH_MSG_INFO( "REGTEST TauDetails->massTrkSys() returns " << massTrkSys);}
1735
1736 if ( (*tauIt)->detail(xAOD::TauJetParameters::PSSFraction,PSSFraction))
1737 ATH_MSG_INFO( "REGTEST TauDetails->PSSFraction() returns " << PSSFraction);
1738
1739 if ( (*tauIt)->detail(xAOD::TauJetParameters::EMPOverTrkSysP,EMPOverTrkSysP))
1740 ATH_MSG_INFO( "REGTEST TauDetails->EMPOverTrkSysP() returns " << EMPOverTrkSysP);
1741
1742 if ( (*tauIt)->detail(xAOD::TauJetParameters::ChPiEMEOverCaloEME,ChPiEMEOverCaloEME))
1743 ATH_MSG_INFO( "REGTEST TauDetails->ChPiEMEOverCaloEME() returns " << ChPiEMEOverCaloEME);
1744
1745 if ( (*tauIt)->detail(xAOD::TauJetParameters::etEMAtEMScale,EtEm)){
1746 EtEm /=1000;
1747 ATH_MSG_INFO( "REGTEST TauDetails->etEMAtEMScale() returns " << EtEm);}
1748
1749 if ( (*tauIt)->detail(xAOD::TauJetParameters::etHadAtEMScale,EtHad)){
1750 EtHad /=1000;
1751 ATH_MSG_INFO( "REGTEST TauDetails->etHadAtEMScale() returns " << EtHad);}
1752
1753
1754 if( !(*tauIt)->jetLink().isValid() ) {
1755 ATH_MSG_WARNING("tau does not have jet seed");
1756 return StatusCode::SUCCESS;
1757 }
1758
1759 const xAOD::Jet* pJetSeed = *((*tauIt)->jetLink());
1760
1761 xAOD::JetConstituentVector::const_iterator clusItr = pJetSeed->getConstituents().begin();
1762 xAOD::JetConstituentVector::const_iterator clusItrE = pJetSeed->getConstituents().end();
1763
1764 for (int clusCount = 0; clusItr != clusItrE; ++clusItr, ++clusCount) {
1765
1766 ATH_MSG_INFO( "REGTEST Tau Cluster " << clusCount << " pt = " << (*clusItr)->pt()
1767 << " eta = " << (*clusItr)->eta()
1768 << " phi = " << (*clusItr)->phi() );
1769
1770 }
1771
1772
1773 for (unsigned int trackNum = 0; trackNum < (*tauIt)->nTracks(); ++trackNum) {
1774
1775 const xAOD::TrackParticle *linkTrack = (*tauIt)->track(trackNum)->track();
1776 if (!linkTrack) {
1777 ATH_MSG_WARNING("can't get tau linked track");
1778 return StatusCode::SUCCESS;
1779 } else {
1780 ATH_MSG_DEBUG("Got the tau linked track");
1781 }
1782
1783 ATH_MSG_INFO( "REGTEST Tau linked track " << trackNum << " pt = " << linkTrack->pt()
1784 << " eta = " << linkTrack->eta()
1785 << " phi = " << linkTrack->phi() );
1786
1787 }
1788
1789
1790 }// end for
1791
1792 return StatusCode::SUCCESS;
1793
1794}
1795
1797
1798StatusCode TrigEDMChecker::dumpxAODTrackParticle() {
1799
1800 ATH_MSG_DEBUG("In dumpxAODTrackParticle()");
1801
1802 ATH_MSG_INFO("REGTEST ==========START of xAOD::TrackParticle DUMP===========");
1803
1804 std::vector<std::string> SGkeys;
1805 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Bjet_EFID");
1806 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Bphysics_EFID");
1807 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Electron_EFID");
1808 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_FullScan_EFID");
1809 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Muon_EFID");
1810 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Photon_EFID");
1811 SGkeys.push_back("HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Tau_EFID");
1812
1813 StatusCode returnsc = StatusCode::SUCCESS;
1814
1815 for (unsigned int SGkey = 0; SGkey < SGkeys.size(); ++SGkey) {
1816 const xAOD::TrackParticleContainer* trackParticleContainer=0;
1817 StatusCode sc = evtStore()->retrieve(trackParticleContainer,SGkeys.at(SGkey));
1818 if (sc.isFailure()) {
1819 ATH_MSG_INFO("REGTEST No track particle container found with key " << SGkeys.at(SGkey));
1820 continue;
1821 }
1822 ATH_MSG_INFO("REGTEST TrackParticleContainer retrieved with key " << SGkeys.at(SGkey)
1823 << " and size " << trackParticleContainer->size());
1824
1825 xAOD::TrackParticleContainer::const_iterator trackParticleItr = trackParticleContainer->begin();
1826 xAOD::TrackParticleContainer::const_iterator trackParticleLast = trackParticleContainer->end();
1827
1828 for (int index = 0; trackParticleItr != trackParticleLast; ++trackParticleItr, ++index) {
1829 ATH_MSG_INFO("REGTEST Looking at Track Particle " << index);
1830
1831 ATH_MSG_INFO("REGTEST IParticle functions:");
1832 ATH_MSG_INFO("REGTEST pt: " << (*trackParticleItr)->pt()
1833 << "/eta: " << (*trackParticleItr)->eta()
1834 << "/phi: " << (*trackParticleItr)->phi()
1835 << "/m: " << (*trackParticleItr)->m()
1836 << "/e: " << (*trackParticleItr)->e()
1837 << "/rapidity: " << (*trackParticleItr)->rapidity());
1838
1839 ATH_MSG_INFO("REGTEST Defining parameters functions:");
1840 ATH_MSG_INFO("REGTEST charge: " << (*trackParticleItr)->charge()
1841 << "/d0: " << (*trackParticleItr)->d0()
1842 << "/z0: " << (*trackParticleItr)->z0()
1843 << "/phi0: " << (*trackParticleItr)->phi0()
1844 << "/theta: " << (*trackParticleItr)->theta()
1845 << "/qOverP: " << (*trackParticleItr)->qOverP()
1846 << "/vx: " << (*trackParticleItr)->vx()
1847 << "/vy: " << (*trackParticleItr)->vy()
1848 << "/vz: " << (*trackParticleItr)->vz());
1849
1850 // Curvilinear functions skipped
1851
1852 ATH_MSG_INFO("REGTEST Fit quality functions:");
1853 ATH_MSG_INFO("REGTEST chiSquared: " << (*trackParticleItr)->chiSquared()
1854 << "/numberDoF: " << (*trackParticleItr)->numberDoF());
1855
1856 // TrackInfo functions skipped
1857
1858 ATH_MSG_INFO("REGTEST summaryValue variables:");
1859 msg() << MSG::INFO << "REGTEST ";
1860 uint8_t numberOfBLayerHits = 0;
1861 if ( (*trackParticleItr)->summaryValue(numberOfBLayerHits, xAOD::numberOfBLayerHits) ) {
1862 msg() << "/numberOfBLayerHits: " << static_cast<int>(numberOfBLayerHits);
1863 } else {
1864 msg() << "/numberOfBLayerHits not found";
1865 }
1866
1867 uint8_t numberOfPixelHits = 0;
1868 if ( (*trackParticleItr)->summaryValue(numberOfPixelHits, xAOD::numberOfPixelHits) ) {
1869 msg() << "/numberOfPixelHits: " << static_cast<int>(numberOfPixelHits);
1870 } else {
1871 msg() << "/numberOfPixelHits not found";
1872 }
1873
1874 uint8_t numberOfPixelHoles = 0;
1875 if ( (*trackParticleItr)->summaryValue(numberOfPixelHoles, xAOD::numberOfPixelHoles) ) {
1876 msg() << "/numberOfPixelHoles: " << static_cast<int>(numberOfPixelHoles);
1877 } else {
1878 msg() << "/numberOfPixelHoles not found";
1879 }
1880
1881 uint8_t numberOfSCTHits = 0;
1882 if ( (*trackParticleItr)->summaryValue(numberOfSCTHits, xAOD::numberOfSCTHits) ) {
1883 msg() << "/numberOfSCTHits: " << static_cast<int>(numberOfSCTHits);
1884 } else {
1885 msg() << "/numberOfSCTHits not found";
1886 }
1887
1888 uint8_t numberOfSCTHoles = 0;
1889 if ( (*trackParticleItr)->summaryValue(numberOfSCTHoles, xAOD::numberOfSCTHoles) ) {
1890 msg() << "/numberOfSCTHoles: " << static_cast<int>(numberOfSCTHoles);
1891 } else {
1892 msg() << "/numberOfSCTHoles not found";
1893 }
1894
1895 uint8_t numberOfTRTHits = 0;
1896 if ( (*trackParticleItr)->summaryValue(numberOfTRTHits, xAOD::numberOfTRTHits) ) {
1897 msg() << "/numberOfTRTHits: " << static_cast<int>(numberOfTRTHits);
1898 } else {
1899 msg() << "/numberOfTRTHits not found";
1900 }
1901
1902 uint8_t numberOfTRTHoles = 0;
1903 if ( (*trackParticleItr)->summaryValue(numberOfTRTHoles, xAOD::numberOfTRTHoles) ) {
1904 msg() << "/numberOfTRTHoles: " << static_cast<int>(numberOfTRTHoles);
1905 } else {
1906 msg() << "/numberOfTRTHoles not found";
1907 }
1908 msg() << endmsg;
1909 }
1910 }
1911
1912 ATH_MSG_INFO("REGTEST ==========END of xAOD::TrackParticle DUMP===========");
1913 ATH_MSG_DEBUG("dumpxAODTrackParticles() succeeded");
1914
1915 return returnsc;
1916}
1917
1918StatusCode TrigEDMChecker::dumpxAODVertex() {
1919
1920 ATH_MSG_DEBUG("In dumpxAODVertex()");
1921
1922 ATH_MSG_INFO("REGTEST ==========START of xAOD::Vertex DUMP===========");
1923
1924 const xAOD::VertexContainer* vertexContainer=0;
1925 StatusCode sc = evtStore()->retrieve(vertexContainer,"HLT_xAOD__VertexContainer_xPrimVx");
1926 if (sc.isFailure()) {
1927 ATH_MSG_INFO("REGTEST No vertex container");
1928 return StatusCode::FAILURE;
1929 }
1930 ATH_MSG_INFO("REGTEST VertexContainer retrieved");
1931
1932 xAOD::VertexContainer::const_iterator vertexItr = vertexContainer->begin();
1933 xAOD::VertexContainer::const_iterator vertexLast = vertexContainer->end();
1934
1935 for (int index = 0; vertexItr != vertexLast; ++vertexItr, ++index) {
1936 ATH_MSG_INFO("REGTEST Looking at Vertex " << index);
1937
1938 ATH_MSG_INFO("REGTEST Public Member Functions:");
1939 ATH_MSG_INFO("REGTEST x: " << (*vertexItr)->x()
1940 << "/y: " << (*vertexItr)->y()
1941 << "/z: " << (*vertexItr)->z());
1942
1943 ATH_MSG_INFO("REGTEST Public Member Functions:");
1944 ATH_MSG_INFO("REGTEST chiSquared: " << (*vertexItr)->chiSquared()
1945 << "/numberDoF: " << (*vertexItr)->numberDoF());
1946 }
1947
1948 ATH_MSG_INFO("REGTEST ==========END of xAOD::Vertex DUMP===========");
1949 ATH_MSG_DEBUG("dumpxAODVertex() succeeded");
1950
1951 return StatusCode::SUCCESS;
1952}
1953
1954StatusCode TrigEDMChecker::dumpTDT(const EventContext& ctx) {
1955 using namespace TrigCompositeUtils; // LinkInfo
1956 ATH_MSG_INFO( "REGTEST ==========START of TDT DUMP===========" );
1957 // Note: This minimal TDT dumper is for use during run-3 dev
1958 std::string chain = m_dumpNavForChain;
1959 if (chain.empty()) {
1960 chain = "HLT_.*";
1961 }
1962 std::vector<std::string> confChains = m_trigDec->getListOfTriggers(chain);
1963 for (const auto& item : confChains) {
1964 bool passed = m_trigDec->isPassed(item);
1965 ATH_MSG_INFO(" HLT Item " << item << " (numeric ID " << TrigConf::HLTUtils::string2hash(item, "Identifier") << ") passed raw? " << passed);
1966 if (m_trigDec->getNavigationFormat() == "TriggerElement") {
1967 ATH_MSG_DEBUG(" Skipping Run 2 features in this dumper");
1968 continue;
1969 }
1970 std::vector< LinkInfo<xAOD::IParticleContainer> > passFeatures = m_trigDec->features<xAOD::IParticleContainer>(item);
1971 if (passFeatures.size()) {
1972 ATH_MSG_INFO(" " << item << " Passed Final IParticle features size: " << passFeatures.size());
1973 for (const LinkInfo<xAOD::IParticleContainer>& li : passFeatures) {
1974 if (!li.isValid()) {
1975 ATH_MSG_WARNING(" Unable to access feature - link invalid.");
1976 } else {
1977 try {
1978 std::string state = "ACTIVE";
1979 if (li.state == ActiveState::INACTIVE) state = "INACTIVE";
1980 else if (li.state == ActiveState::UNSET) state = "UNSET";
1981 ATH_MSG_INFO(" IParticle Feature from " << li.link.dataID() << " index:" << li.link.index() << " pt:" << (*li.link)->pt() << " eta:" << (*li.link)->eta() << " phi:" << (*li.link)->phi() << " state:" << state);
1982 } catch (const std::exception& e) {
1983 ATH_MSG_WARNING(" Unable to dereference feature {" << e.what() << "}");
1984 }
1985 }
1986 }
1987 }
1988 std::vector< LinkInfo<xAOD::IParticleContainer> > passAndFailFeatures = m_trigDec->features<xAOD::IParticleContainer>(item, TrigDefs::includeFailedDecisions);
1989 if (passAndFailFeatures.size()) {
1990 ATH_MSG_INFO(" " << item << " Passed+Failed Final IParticle features size: " << passAndFailFeatures.size());
1991 for (const LinkInfo<xAOD::IParticleContainer>& li : passAndFailFeatures) {
1992 if (!li.isValid()) {
1993 ATH_MSG_WARNING(" Unable to access feature - link invalid.");
1994 } else {
1995 try {
1996 std::string state = "ACTIVE";
1997 if (li.state == ActiveState::INACTIVE) state = "INACTIVE";
1998 else if (li.state == ActiveState::UNSET) state = "UNSET";
1999 ATH_MSG_INFO(" IParticle Feature from " << li.link.dataID() << " index:" << li.link.index() << " pt:" << (*li.link)->pt() << " eta:" << (*li.link)->eta() << " phi:" << (*li.link)->phi() << " state:" << state);
2000 } catch (const std::exception& e) {
2001 ATH_MSG_WARNING(" Unable to dereference feature {" << e.what() << "}");
2002 }
2003 }
2004 }
2005 }
2006 std::vector< LinkInfo<xAOD::IParticleContainer> > allFeatures = m_trigDec->features<xAOD::IParticleContainer>(item, TrigDefs::includeFailedDecisions, "", TrigDefs::allFeaturesOfType);
2007 if (allFeatures.size()) {
2008 ATH_MSG_INFO(" " << item << " Passed+Failed ALL IParticle features size: " << allFeatures.size());
2009 for (const LinkInfo<xAOD::IParticleContainer>& li : allFeatures) {
2010 if (!li.isValid()) {
2011 ATH_MSG_WARNING(" Unable to access feature - link invalid.");
2012 } else {
2013 try {
2014 std::string state = "ACTIVE";
2015 if (li.state == ActiveState::INACTIVE) state = "INACTIVE";
2016 else if (li.state == ActiveState::UNSET) state = "UNSET";
2017 ATH_MSG_INFO(" IParticle Feature from " << li.link.dataID() << " index:" << li.link.index() << " pt:" << (*li.link)->pt() << " eta:" << (*li.link)->eta() << " phi:" << (*li.link)->phi() << " state:" << state);
2018 } catch (const std::exception& e) {
2019 ATH_MSG_WARNING(" Unable to dereference feature {" << e.what() << "}");
2020 }
2021 }
2022 }
2023 }
2024 }
2025
2026 if (m_trigDec->getNavigationFormat() == "TrigComposite") {
2027 // Check associateToEventView helper function
2028 std::vector< LinkInfo<xAOD::IParticleContainer> > muons = m_trigDec->features<xAOD::IParticleContainer>("HLT_mu24_idperf_L1MU20", TrigDefs::Physics, "HLT_MuonL2CBInfo");
2029 SG::ReadHandle<xAOD::TrackParticleContainer> muonTracksReadHandle(m_muonTracksKey, ctx);
2030 for (const LinkInfo<xAOD::IParticleContainer>& mu : muons) {
2031 // Note: auto here refers to type std::pair< xAOD::TrackParticleContainer::const_iterator, xAOD::TrackParticleContainer::const_iterator>
2032 const auto roiTrackItPair = m_trigDec->associateToEventView<xAOD::TrackParticleContainer>(muonTracksReadHandle, mu, "roi");
2033 const xAOD::TrackParticleContainer::const_iterator startIt = roiTrackItPair.first;
2034 const xAOD::TrackParticleContainer::const_iterator stopIt = roiTrackItPair.second;
2035 ATH_MSG_INFO("Muon pT: " << (*mu.link)->pt() << " is from the same ROI as tracks with index "
2036 << std::distance(muonTracksReadHandle->begin(), startIt) << "-" << std::distance(muonTracksReadHandle->begin(), stopIt)
2037 << ", which is " << std::distance(startIt, stopIt) << " tracks, out of " << muonTracksReadHandle->size() << " total tracks.");
2038 for (xAOD::TrackParticleContainer::const_iterator it = startIt; it != stopIt; ++it) {
2039 ATH_MSG_VERBOSE(" -- Track " << std::distance(startIt, it) << " in this ROI, pT: " << (*it)->pt() );
2040 }
2041 }
2042 }
2043
2044 ATH_MSG_INFO( "REGTEST ==========END of TDT DUMP===========" );
2045 return StatusCode::SUCCESS;
2046}
2047
2048StatusCode TrigEDMChecker::dumpTrigComposite() {
2049 ATH_MSG_INFO( "REGTEST ==========START of xAOD::TrigCompositeContainer DUMP===========" );
2050
2051 if (m_doDumpAllTrigComposite) {
2052 m_dumpTrigCompositeContainers.clear();
2053 const CLID TrigCompositeCLID = static_cast<CLID>( ClassID_traits< xAOD::TrigCompositeContainer >::ID() );
2054 evtStore()->keys(TrigCompositeCLID, m_dumpTrigCompositeContainers.value());
2055 std::string typeNameTC;
2056 ATH_CHECK(m_clidSvc->getTypeNameOfID(TrigCompositeCLID, typeNameTC));
2057 ATH_MSG_DEBUG("dumpTrigComposite got " << m_dumpTrigCompositeContainers.size() << " keys for " << typeNameTC);
2058 } else {
2059 ATH_MSG_DEBUG("Using supplied " << m_dumpTrigCompositeContainers.size() << " keys");
2060 }
2061
2062 for ( const std::string & key: m_dumpTrigCompositeContainers ) {
2063 // get the collection
2064 if ( not evtStore()->contains<xAOD::TrigCompositeContainer>(key) ) {
2065 ATH_MSG_WARNING("Absent TrigCompositeContainer: " << key );
2066 continue;
2067 }
2068 ATH_MSG_DEBUG( "#################### Dumping container of : " << key );
2069 const xAOD::TrigCompositeContainer* cont= nullptr;
2070 ATH_CHECK( evtStore()->retrieve( cont, key ) );
2071
2072 size_t count = 0;
2073 for ( auto tc: *cont ) {
2074 ATH_MSG_DEBUG("########## ELEMENT " << count++);
2075 ATH_MSG_DEBUG(*tc);
2076 // Get the objects we know of
2077 for (size_t i = 0; i < tc->linkColNames().size(); ++i) ATH_CHECK(checkTrigCompositeElementLink(tc, i));
2078 }
2079 }
2080 ATH_MSG_INFO( "REGTEST ==========END of xAOD::TrigCompositeContainer DUMP===========" );
2081 return StatusCode::SUCCESS;
2082}
2083
2084
2085
2086StatusCode TrigEDMChecker::checkTrigCompositeElementLink(const xAOD::TrigComposite* tc, size_t element) {
2087
2088 const std::string name = tc->linkColNames().at(element);
2089 const CLID clid = static_cast<CLID>(tc->linkColClids().at(element));
2090
2091 if (clid == ClassID_traits< TrigRoiDescriptorCollection >::ID()) {
2092
2093 const ElementLink<TrigRoiDescriptorCollection> elementLink = tc->objectLink<TrigRoiDescriptorCollection>(name);
2094 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to TrigRoiDescriptorCollection, link name:'" << name << "'");
2095 else ATH_MSG_DEBUG(" Dereferenced link '" << name << "'' to TrigRoiDescriptor:" << *elementLink);
2096
2097 } else if (clid == ClassID_traits< DataVector< LVL1::RecEmTauRoI > >::ID()) { // There could be a few ROI types....
2098 // CLASS_DEF( DataVector< LVL1::RecEmTauRoI >, 6256, 1 )
2099
2100 const ElementLink<DataVector< LVL1::RecEmTauRoI >> elementLink = tc->objectLink<DataVector< LVL1::RecEmTauRoI >>(name);
2101 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to LVL1::RecEmTauRoI, link name:'" << name << "'");
2102 else ATH_MSG_DEBUG(" Dereferenced link '" << name << "' to LVL1::RecEmTauRoI:" << *elementLink);
2103
2104 } else if (clid == ClassID_traits< xAOD::TrigCompositeContainer >::ID()) {
2105
2106 const ElementLink<xAOD::TrigCompositeContainer> elementLink = tc->objectLink<xAOD::TrigCompositeContainer>(name);
2107 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to TrigComposite, link name:'" << name << "'");
2108 else ATH_MSG_DEBUG(" Dereferenced link '" << name << "' to TrigComposite, TC name:'" << (*elementLink)->name() << "'");
2109
2110 } else if (clid == ClassID_traits< ViewContainer >::ID()) {
2111
2112 const ElementLink<ViewContainer> elementLink = tc->objectLink<ViewContainer>(name);
2113 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to View, link name:'" << name << "'");
2114 else ATH_MSG_DEBUG(" Dereferenced link '" << name << "' to View:'" << *elementLink);
2115
2116 } else if (name == "feature") {
2117
2118 if (clid == ClassID_traits< xAOD::TrigEMClusterContainer >::ID()) {
2119
2120 const ElementLink<xAOD::TrigEMClusterContainer> elementLink = tc->objectLink<xAOD::TrigEMClusterContainer>(name);
2121 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to xAOD::TrigEMClusterContainer 'feature'");
2122 else ATH_MSG_DEBUG(" Dereferenced xAOD::TrigEMClusterContainer link 'feature', Energy:" << (*elementLink)->energy());
2123
2124 } else if (clid == ClassID_traits< xAOD::TrigMissingETContainer >::ID()) {
2125
2126 const ElementLink<xAOD::TrigMissingETContainer> elementLink = tc->objectLink<xAOD::TrigMissingETContainer>(name);
2127 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to xAOD::TrigMissingETContainer 'feature'");
2128 else ATH_MSG_DEBUG(" Dereferenced xAOD::TrigMissingETContainer link 'feature', ex:" << (*elementLink)->ex() << " ey:" << (*elementLink)->ey());
2129
2130 } else {
2131
2132 try {
2133 const ElementLink<xAOD::IParticleContainer> elementLink = tc->objectLink<xAOD::IParticleContainer>(name);
2134 if (!elementLink.isValid()) ATH_MSG_WARNING(" Invalid element link to 'feature'");
2135 else ATH_MSG_DEBUG(" Dereferenced IParticle link 'feature', pt:" << (*elementLink)->pt() << " eta:" << (*elementLink)->eta() << " phi:" << (*elementLink)->phi());
2136 } catch(std::runtime_error& e) {
2137 ATH_MSG_WARNING(" Cannot dereference 'feature' as IParticle: '" << e.what() << "'");
2138 }
2139
2140 }
2141
2142 } else {
2143 ATH_MSG_DEBUG(" Ignoring link to '" << name << "' with link CLID " << clid);
2144 }
2145
2146 return StatusCode::SUCCESS;
2147
2148}
2149
2150
2151StatusCode TrigEDMChecker::TrigCompositeNavigationToDot(std::string& returnValue, bool& pass) {
2152
2153 using namespace TrigCompositeUtils;
2154
2155 // This constexpr is evaluated at compile time
2156 const CLID TrigCompositeCLID = static_cast<CLID>( ClassID_traits< xAOD::TrigCompositeContainer >::ID() );
2157 std::vector<std::string> keys;
2158 if ( m_dumpTrigCompositeContainers.size() == 0 ) {
2159 evtStore()->keys(TrigCompositeCLID, keys);
2160 }
2161 else {
2162 keys = m_dumpTrigCompositeContainers;
2163 }
2164 std::string typeNameTC;
2165 ATH_CHECK(m_clidSvc->getTypeNameOfID(TrigCompositeCLID, typeNameTC));
2166 ATH_MSG_DEBUG("Got " << keys.size() << " keys for " << typeNameTC);
2167
2168 HLT::Identifier chainID(m_dumpNavForChain);
2169 DecisionIDContainer chainIDs;
2170 chainIDs.insert( chainID.numeric() );
2171
2172 std::set<int> converted;
2173
2174 const Trig::ChainGroup* cg = m_trigDec->getChainGroup(m_dumpNavForChain);
2175 pass = cg->isPassed(TrigDefs::requireDecision);
2176 std::vector<std::string> chains = cg->getListOfTriggers();
2177 for (const std::string& chain : chains) {
2178 const TrigConf::HLTChain* hltChain = m_trigDec->ExperimentalAndExpertMethods().getChainConfigurationDetails(chain);
2179 const HLT::Identifier chainID_tmp( hltChain->chain_name() );
2180 chainIDs.insert( chainID_tmp.numeric() );
2181 const std::vector<size_t> legMultiplicites = hltChain->leg_multiplicities();
2182 if (legMultiplicites.size() == 0) {
2183 ATH_MSG_ERROR("chain " << chainID_tmp << " has invalid configuration, no multiplicity data.");
2184 } else if (legMultiplicites.size() > 1) {
2185 // For multi-leg chains, the DecisionIDs are handled per leg.
2186 // We don't care here exactly how many objects are required per leg, just that there are two-or-more legs
2187 for (size_t legNumeral = 0; legNumeral < legMultiplicites.size(); ++legNumeral) {
2188 const HLT::Identifier legID = TrigCompositeUtils::createLegName(chainID_tmp, legNumeral);
2189 chainIDs.insert( legID.numeric() );
2190 }
2191 }
2192 }
2193
2194 // First retrieve them all (this should not be needed in future)
2195 const DecisionContainer* container = nullptr;
2196 for (const std::string& key : keys) ATH_CHECK( evtStore()->retrieve( container, key ) );
2197
2198 std::stringstream ss;
2199 ss << "digraph {" << std::endl;
2200 ss << " node [shape=rectangle]" << std::endl;
2201 ss << " rankdir = BT" << std::endl;
2202
2203 // Now process them
2204 for (const std::string& key : keys) {
2205 if ( not m_doDumpAllTrigComposite ) {
2206 if ( not key.starts_with( "HLTNav_") ) { // Nav containers should always start with HLTNav_
2207 continue;
2208 }
2209 }
2210 ATH_CHECK( evtStore()->retrieve( container, key ) );
2211 ATH_MSG_DEBUG("Processing collection " << key << " to be added to the navigation graph");
2212 // ss << " rank=same" << std::endl; // dot cannot handle this is seems
2213 bool writtenHeader = false;
2214 for (const Decision* tc : *container ) {
2215 // Output my ID in the graph.
2216 const DecisionContainer* container = dynamic_cast<const DecisionContainer*>( tc->container() );
2217 const ElementLink<DecisionContainer> selfEL = ElementLink<DecisionContainer>(*container, tc->index());
2218 std::vector<ElementLink<DecisionContainer>> seedELs = tc->objectCollectionLinks<DecisionContainer>("seed");
2219 const bool isHypoAlgNode = tc->name() == "H";
2220 const bool isComboHypoAlgNode = tc->name() == "CH";
2221 const std::vector<DecisionID>& decisions = tc->decisions();
2222 const uint32_t selfKey = selfEL.key();
2223 const uint32_t selfIndex = selfEL.index();
2224 if (m_dumpNavForChain != "") {
2225 bool doDump = false;
2226 // Check me
2227 for (DecisionID id : decisions) {
2228 if (chainIDs.count(id) == 1) {
2229 doDump = true;
2230 break;
2231 }
2232 }
2233 // Check my seeds
2234 if (!doDump and (isHypoAlgNode or isComboHypoAlgNode) and not m_excludeFailedHypoNodes) {
2235 for (const ElementLink<DecisionContainer>& s : seedELs) {
2236 const std::vector<DecisionID>& seedDecisions = (*s)->decisions();
2237 for (DecisionID id : seedDecisions) {
2238 if (chainIDs.count(id) == 1) {
2239 doDump = true;
2240 break;
2241 }
2242 }
2243 }
2244 }
2245 if (!doDump) {
2246 continue;
2247 }
2248 }
2249 if (!writtenHeader) {
2250 writtenHeader = true;
2251 ss << " subgraph " << key << " {" << std::endl;
2252 ss << " label=\"" << key << "\"" << std::endl;
2253 }
2254 const std::string scheme = "rdpu9";
2255 std::string color = "1";
2256 if (tc->name() == "L1") { color = "1"; }
2257 else if (tc->name() == "F") { color = "2"; }
2258 else if (tc->name() == "IM") { color = "3"; }
2259 else if (tc->name() == "H") { color = "4"; }
2260 else if (tc->name() == "CH") { color = "5"; }
2261 else if (tc->name() == "SF") { color = "6"; }
2262 else if (tc->name() == "HLTPassRaw") { color = "7"; }
2263 ss << " \"" << selfKey << "_" << selfIndex << "\" [colorscheme="<<scheme<<",style=filled,fillcolor="<<color<<",label=<<B>Container</B>=" << typeNameTC;
2264 if (tc->name() != "") ss << " <B>Name</B>=" << tc->name();
2265 ss << "<BR/><B>Key</B>=" << key << "<BR/><B>Index</B>=" << selfIndex;
2266 const bool isRemapped = tc->isRemapped();
2267 if (isHypoAlgNode) ss << " <B>linksRemapped</B>=" << (isRemapped ? "Y" : "N");
2268 if (decisions.size() > 0) {
2269 ss << "<BR/><B>Pass</B>=";
2270 size_t c = 0;
2271 for (unsigned decisionID : decisions) {
2272 HLT::Identifier dID(decisionID);
2273 std::string highlight = (dID.numeric() == chainID.numeric() ? "<B>[CHAIN:" : "");
2274 if (highlight == "" and chainIDs.count(dID.numeric()) == 1 and TrigCompositeUtils::isLegId(dID)) {
2275 highlight = "<B>[LEG" + std::to_string(TrigCompositeUtils::getIndexFromLeg(dID)) + ":";
2276 }
2277 ss << std::hex << highlight << decisionID << (!highlight.empty() ? "]</B>" : "") << std::dec << ",";
2278 if (c++ == 5) {
2279 ss << "<BR/>";
2280 c = 0;
2281 }
2282 }
2283 }
2284 ss << ">]" << std::endl;
2285 // Output all the things I link to
2286 size_t seedCount = 0;
2287 for (size_t i = 0; i < tc->linkColNames().size(); ++i) {
2288 const std::string link = tc->linkColNames().at(i);
2289 if (link == "seed" || link == "seed__COLL") {
2290 ElementLink<DecisionContainer> seedEL = seedELs.at(seedCount++);
2291 const uint32_t seedKey = tc->linkColKeys().at(i);
2292 const uint32_t seedIndex = tc->linkColIndices().at(i);
2293 ATH_CHECK( seedKey == seedEL.key() );
2294 ATH_CHECK( seedIndex == seedEL.index() );
2295 if (m_dumpNavForChain != "") { // Only print "seed" link to nodes we include in our search
2296 const std::vector<DecisionID> seedDecisions = (*seedEL)->decisions();
2297 bool doSeedLink = false;
2298 for (DecisionID id : seedDecisions) {
2299 if (chainIDs.count(id) == 1) {
2300 doSeedLink = true;
2301 break;
2302 }
2303 }
2304 if (!doSeedLink) {
2305 continue;
2306 }
2307 }
2308 ss << " \"" << selfKey << "_" << selfIndex << "\" -> \"" << seedKey << "_" << seedIndex << "\" [colorscheme="<<scheme<<",color=9,fontcolor=8,label=\"seed\"]" << std::endl;
2309 } else {
2310 // Start with my class ID
2311 std::string linkColour = "12";
2312 std::string linkBackground = "11";
2313 const std::string extScheme = "paired12";
2314 if (link == "roi") { linkColour="2"; linkBackground="1"; }
2315 else if (link == "initialRoI") { linkColour="2"; linkBackground="1"; }
2316 else if (link == "initialRecRoI") { linkColour="8"; linkBackground="7"; }
2317 else if (link == "feature") { linkColour="4"; linkBackground="3"; }
2318 else if (link == "view") { linkColour="10"; linkBackground="9"; }
2319 const CLID linkCLID = static_cast<CLID>( tc->linkColClids().at(i) );
2320 // Use it to get my class name
2321 std::string tname;
2322 ATH_CHECK(m_clidSvc->getTypeNameOfID(linkCLID, tname));
2323 // Now get the sgkey I'm linking to & the index
2324 const SG::sgkey_t key = (isRemapped ? static_cast<SG::sgkey_t>( tc->linkColKeysRemap().at(i) ) : static_cast<SG::sgkey_t>( tc->linkColKeys().at(i) ));
2325 const unsigned index = (isRemapped ? tc->linkColIndicesRemap().at(i) : tc->linkColIndices().at(i));
2326 // Look it up
2327 CLID checkCLID;
2328 const std::string* keyStr = evtStore()->keyToString(key, checkCLID); // I don't own this str
2329 if (keyStr != nullptr && checkCLID != linkCLID) {
2330 std::string tnameOfCheck;
2331 m_clidSvc->getTypeNameOfID(checkCLID, tnameOfCheck).ignore(); // Might be invalid. But we don't care.
2332 ATH_MSG_ERROR("Inconsistent CLID " << checkCLID << " [" << tnameOfCheck << "] stored in storegate for key " << key
2333 << ". We were expecting " << linkCLID << " [" << tname << "]");
2334 }
2335
2336 std::string tnameEscape;
2337 for (std::string::const_iterator i = tname.begin(); i != tname.end(); ++i) {
2338 unsigned char c = *i;
2339 if (c == '<') {
2340 tnameEscape += "&lt;";
2341 } else if (c == '>') {
2342 tnameEscape += "&gt;";
2343 } else {
2344 tnameEscape += c;
2345 }
2346 }
2347
2348 // Print
2349 ss << " \"" << selfKey << "_" << selfIndex << "\" -> \"" << key << "_" << index << "\" ";
2350 ss << "[colorscheme="<<extScheme<<",color="<<linkColour<<",fontcolor="<<linkColour<<",arrowhead=empty,label=\"" << link << "\"]" << std::endl;
2351
2352 // Check if we are linking to self (e.g. a dummy-feature), don't output a new box for this
2353 const bool linkToSelf = (selfKey == key and selfIndex == index);
2354
2355 if (converted.count(key + index) == 0 and not linkToSelf) {
2356 ss << " \"" << key << "_" << index << "\" [colorscheme="<<extScheme<<",style=filled,fillcolor="<<linkBackground<<",label=<<B>Container</B>=" << tnameEscape << "<BR/><B>Key</B>=";
2357 if (keyStr != nullptr) ss << *keyStr;
2358 else ss << "[<I>KEY "<< key <<" NOT IN STORE</I>] ";
2359 ss << "<BR/><B>Index</B>=" << index << ">]";
2360 }
2361
2362 converted.insert(key + index);
2363 }
2364 }
2365 }
2366 if (writtenHeader) {
2367 ss << " }" << std::endl;
2368 }
2369 }
2370
2371 ss << "}" << std::endl;
2372
2373 returnValue.assign( ss.str() );
2374 return StatusCode::SUCCESS;
2375}
2376
2377StatusCode TrigEDMChecker::dumpNavigation(const EventContext& ctx)
2378{
2379 // Get object from store
2380 const xAOD::TrigNavigation * navigationHandle = nullptr;
2381 ATH_CHECK( evtStore()->retrieve( navigationHandle, m_navigationHandleKey.key() ) );
2382 // Proper version doesn't work - conversion issue?
2383 //SG::ReadHandle< xAOD::TrigNavigation > navigationHandle = SG::ReadHandle< xAOD::TrigNavigation >( m_navigationHandleKey );
2384 //if ( !navigationHandle.isValid() ) ATH_MSG_FATAL( "Could not retrieve navigation" );
2385
2386 // Get serialised navigation info
2387 const std::vector< unsigned int > serialisedNavigation = navigationHandle->serialized();
2388 ATH_MSG_INFO( "Serialised navigation size: " << serialisedNavigation.size() );
2389
2390 // Convert the input
2391 HLT::Navigation* testNav = m_navigationTool.get();
2392 testNav->deserialize( serialisedNavigation );
2393
2394 // Make a map of TE name hashes
2395 const xAOD::TriggerMenuContainer * testMenu = nullptr;
2396 ATH_CHECK( inputMetaStore()->retrieve( testMenu, "TriggerMenu" ) );
2397 std::map< int, std::string > hash2string;
2398 for ( auto const& sequence : testMenu->front()->sequenceInputTEs() ) {
2399 for ( auto const& name : sequence ) {
2400 int hash = TrigConf::HLTUtils::string2hash( name );
2401 hash2string[ hash ] = name;
2402 }
2403 }
2404
2405 // Map TE names to chain names
2406 unsigned int chainCounter = 0;
2407 std::map< int, std::string > hash2chain;
2408 for ( auto const& chain : testMenu->front()->chainSignatureOutputTEs() ) {
2409
2410 // Find the chain name
2411 std::string chainName = testMenu->front()->chainNames()[ chainCounter ];
2412 ++chainCounter;
2413
2414 // Find all associated TEs
2415 for ( auto const& signature : chain ) {
2416 for ( auto const& name : signature ) {
2417 int hash = TrigConf::HLTUtils::string2hash( name );
2418 hash2string[ hash ] = name; // for decoding
2419 hash2chain[ hash ] = chainName;
2420 }
2421 }
2422 }
2423
2424 // Define a map of TE features, to the TEs that use them. Needs a custom sort lambda
2425 auto cmpLambda = []( const HLT::TriggerElement::FeatureAccessHelper &lhs, const HLT::TriggerElement::FeatureAccessHelper &rhs) {
2426
2427 // Compare indices if CLID matches
2428 if ( lhs.getCLID() == rhs.getCLID() ) return ( lhs.getIndex() < rhs.getIndex() );
2429
2430 // Compare CLIDs
2431 else return ( lhs.getCLID() < rhs.getCLID() );
2432 };
2433 std::map< HLT::TriggerElement::FeatureAccessHelper, std::vector< HLT::TriggerElement* >, decltype(cmpLambda) > feature2element(cmpLambda);
2434
2435 // Retrieve all TE features and add them to the map
2436 std::vector< HLT::TriggerElement* > allTEs;
2437 testNav->getAll( allTEs, false );
2438 for ( auto element : allTEs ) {
2439
2440 // Add TE features to the big map
2441 for ( auto helper : element->getFeatureAccessHelpers() ) {
2442 feature2element[ helper ].push_back( element );
2443 }
2444 }
2445
2446 // Debug - output all TEs and their ancestors
2447 // No duplication - only print terminal nodes
2448 for ( auto element : allTEs ) {
2449 if ( testNav->isTerminalNode( element ) ) {
2450 ATH_MSG_INFO( "+++++++++++ " << hash2string[ element->getId() ] << " is terminal node" );
2451 ATH_MSG_INFO( "ptr: " << element );
2452 std::queue< HLT::TriggerElement* > allAncestors;
2453 allAncestors.push( element );
2454 while ( allAncestors.size() ) {
2455
2456 HLT::TriggerElement * thisElement = allAncestors.front();
2457 allAncestors.pop();
2458 auto theseAncestors = thisElement->getRelated( HLT::TriggerElement::Relation::seededByRelation );
2459
2460 // Dump TE
2461 ATH_MSG_INFO( "te: " << thisElement->getId() << " " << hash2string[ thisElement->getId() ] );
2462 ATH_MSG_INFO( " chain: " << hash2chain[ thisElement->getId() ] );
2463 for ( const auto& helper : thisElement->getFeatureAccessHelpers() ) {
2464 ATH_MSG_INFO( " feat: " << helper );
2465 }
2466 ATH_MSG_INFO( theseAncestors.size() << " ancestors" );
2467
2468 // Examine ancestors
2469 for ( auto ancestor : theseAncestors ) {
2470 allAncestors.push( ancestor );
2471 }
2472 }
2473 }
2474 }
2475
2476 // Make the decision container
2477 SG::WriteHandle< TrigCompositeUtils::DecisionContainer > outputNavigation = TrigCompositeUtils::createAndStore( m_decisionsKey, ctx );
2478 auto decisionOutput = outputNavigation.ptr();
2479
2480 // Find unique chains associated with a feature
2481 std::map< HLT::TriggerElement const*, std::vector< int > > element2decisions;
2482 for ( const auto& pair : feature2element ) {
2483
2484 // Get the feature info
2485 std::string featureName = testNav->label( pair.first.getCLID(), pair.first.getIndex().subTypeIndex() );
2486 auto sgKey = evtStore()->stringToKey( featureName, pair.first.getCLID() );
2487
2488 // Store RoIs with appropriate label ?
2489 std::string storeFeatureName = "feature";
2490/* if ( pair.first.getCLID() == ClassID_traits< TrigRoiDescriptor >::ID() ) {
2491 storeFeatureName = "roi";
2492 }*/
2493
2494 // Make a decision object for the feature
2495 auto decision = TrigCompositeUtils::newDecisionIn( decisionOutput );
2496 decision->typelessSetObjectLink( storeFeatureName, sgKey, pair.first.getCLID(), pair.first.getIndex().objectsBegin(), pair.first.getIndex().objectsEnd() );
2497
2498 // Examine associated TEs, look for chains
2499 std::set< std::string > passedChains;
2500 for ( HLT::TriggerElement const* element : pair.second ) {
2501
2502 // TODO - find out what chains actually passed!
2503 passedChains.insert( hash2chain[ element->getId() ] );
2504
2505 // Index the TE
2506 int decisionNumber = decisionOutput->size() - 1;
2507 element2decisions[ element ].push_back( decisionNumber );
2508 }
2509
2510 // Store unique chains in the decision
2511 for ( auto& chain : passedChains ) {
2512 TrigCompositeUtils::addDecisionID( TrigConf::HLTUtils::string2hash( chain ), decision );
2513 }
2514 }
2515
2516 // Store decision ancestry (had to go through once before to ensure indices populated)
2517 unsigned int decisionCounter = 0;
2518 for ( const auto& pair : feature2element ) {
2519
2520 // Get current decision
2521 auto decision = decisionOutput->at( decisionCounter );
2522 ++decisionCounter;
2523
2524 // Find ancestor TEs
2525 for ( auto element : pair.second ) {
2526 auto theseAncestors = element->getRelated( HLT::TriggerElement::Relation::seededByRelation );
2527 for ( auto ancestor : theseAncestors ) {
2528 for ( int decisionIndex : element2decisions[ ancestor ] ) {
2529 TrigCompositeUtils::linkToPrevious( decision, m_decisionsKey.key(), decisionIndex );
2530 }
2531 }
2532 }
2533 }
2534
2535 return StatusCode::SUCCESS;
2536}
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
theta
Scalar theta() const
theta method
Definition AmgMatrixBasePlugin.h:75
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:54
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition AthMsgStreamMacros.h:33
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition AthMsgStreamMacros.h:28
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition AthMsgStreamMacros.h:32
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
CLID
uint32_t CLID
The Class ID type.
Definition Event/xAOD/xAODCore/xAODCore/ClassID_traits.h:47
linkTrack
ElementLink< xAOD::TrackParticleContainer > linkTrack(const xAOD::TrackParticle *trk)
Definition ITrigBphysState.h:62
JetConstituentVector.h
This file defines helper classes to deal with jet constituents.
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
uint
unsigned int uint
Definition LArOFPhaseFill.cxx:19
ss
static Double_t ss
Definition LArPhysWaveHECTool.cxx:37
tc
static Double_t tc
Definition LArPhysWaveHECTool.cxx:38
sc
static Double_t sc
Definition LArPhysWaveHECTool.cxx:37
DecisionID
unsigned int DecisionID
Definition TrigComposite_v1.h:27
Decision
xAOD::TrigComposite Decision
Definition TrigComposite.h:20
DecisionIDContainer
std::set< DecisionID > DecisionIDContainer
Definition TrigComposite_v1.h:28
DecisionContainer
xAOD::TrigCompositeContainer DecisionContainer
Definition TrigCompositeContainer.h:21
size
size_t size() const
Number of registered mappings.
ViewContainer
DataVector< SG::View > ViewContainer
View container for recording in StoreGate.
Definition View.h:290
AthAnalysisAlgorithm::inputMetaStore
const ServiceHandle< StoreGateSvc > & inputMetaStore() const
Const accessor for the input metadata store.
Definition AthAnalysisAlgorithm.h:49
AthCommonAlgorithm< Gaudi::Algorithm >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
Definition AthCommonDataStore.h:85
AthHistogramming::hash
hash_t hash(const std::string &histName) const
Method to calculate a 32-bit hash from a string.
Definition AthHistogramming.h:430
DataVector
Derived DataVector<T>.
Definition DataVector.h:795
DataVector< xAOD::TrigBphys_v1 >::const_iterator
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
DataVector::at
const T * at(size_type n) const
Access an element, as an rvalue.
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.
HLT::Identifier::numeric
TrigCompositeUtils::DecisionID numeric() const
numeric ID
Definition TrigCompositeUtils/TrigCompositeUtils/HLTIdentifier.h:41
HLT::NavigationCore::deserialize
bool deserialize(const std::vector< uint32_t > &input)
Definition NavigationCore.cxx:151
HLT::Navigation
The Navigation class, organizes TriggerElements into the tree structure.
Definition Navigation.h:100
HLT::TrigNavStructure::getAll
void getAll(std::vector< TriggerElement * > &output, const bool activeOnly=true) const
The query returning a collection of all TriggerElements.
Definition TrigNavStructure.cxx:361
HLT::TrigNavStructure::isTerminalNode
static bool isTerminalNode(const TriggerElement *te)
queries if node is terminal (no more TriggerElement are seeded by it)
Definition TrigNavStructure.cxx:515
HLT::TrigNavStructure::label
std::string label(class_id_type clid, const index_or_label_type &sti_or_label) const
Definition TrigNavStructure.cxx:773
HLT::TriggerElement::FeatureAccessHelper
the FeatureAccessHelper is a class used to keep track of features attached to this TE.
Definition TrigNavStructure/TrigNavStructure/TriggerElement.h:192
HLT::TriggerElement::FeatureAccessHelper::getIndex
const ObjectIndex & getIndex() const
index in the external ojects array
Definition TrigNavStructure/TrigNavStructure/TriggerElement.h:209
HLT::TriggerElement
TriggerElement is the basic ingreedient of the interface between HLT algorithms and the navigation It...
Definition TrigNavStructure/TrigNavStructure/TriggerElement.h:27
HLT::TriggerElement::getRelated
const std::vector< TriggerElement * > & getRelated(Relation rel) const
returns reference to the likns to other TriggerElements related by relation r
Definition TrigNavStructure/TrigNavStructure/TriggerElement.h:127
HLT::TriggerElement::getFeatureAccessHelpers
const std::vector< FeatureAccessHelper > & getFeatureAccessHelpers() const
returns all features which ara attached to this TE
Definition TrigNavStructure/TrigNavStructure/TriggerElement.h:238
LVL1_ROI
Top level AOD object storing LVL1 RoIs.
Definition LVL1_ROI.h:43
LVL1_ROI::getEmTauROIs
const emtaus_type & getEmTauROIs() const
Get all the em/tau RoIs in the event.
Definition LVL1_ROI.h:65
SG::ConstIterator
a const_iterator facade to DataHandle.
Definition SGIterator.h:164
SG::WriteHandle::ptr
pointer_type ptr()
Dereference the pointer.
TrigConf::HLTChain
HLT chain configuration information.
Definition TrigConfHLTData/TrigConfHLTData/HLTChain.h:35
TrigConf::HLTChain::leg_multiplicities
const std::vector< size_t > & leg_multiplicities() const
Definition TrigConfHLTData/TrigConfHLTData/HLTChain.h:83
TrigConf::HLTUtils::string2hash
static HLTHash string2hash(const std::string &, const std::string &category="TE")
hash function translating TE names into identifiers
TrigEDMChecker::m_decisionsKey
SG::WriteHandleKey< TrigCompositeUtils::DecisionContainer > m_decisionsKey
Definition TrigEDMChecker.h:135
TrigEDMChecker::m_doDumpxAODElectronContainer
Gaudi::Property< bool > m_doDumpxAODElectronContainer
Definition TrigEDMChecker.h:67
TrigEDMChecker::dumpxAODTrigEMClusterContainer
StatusCode dumpxAODTrigEMClusterContainer()
Definition TrigEDMChecker.cxx:1599
TrigEDMChecker::m_doTDTCheck
Gaudi::Property< bool > m_doTDTCheck
Definition TrigEDMChecker.h:85
TrigEDMChecker::dumpTrigPassBits
StatusCode dumpTrigPassBits()
Definition TrigEDMChecker.cxx:279
TrigEDMChecker::execute
virtual StatusCode execute(const EventContext &ctx) override
Execute method.
Definition TrigEDMChecker.cxx:119
TrigEDMChecker::dumpTrigEFBphysContainer
StatusCode dumpTrigEFBphysContainer()
Definition TrigEDMChecker.cxx:1013
TrigEDMChecker::dumpxAODTrackParticle
StatusCode dumpxAODTrackParticle()
Definition TrigEDMChecker.cxx:1798
TrigEDMChecker::dumpNavigation
StatusCode dumpNavigation(const EventContext &ctx)
Definition TrigEDMChecker.cxx:2377
TrigEDMChecker::m_muonPrinter
ToolHandle< Rec::IMuonPrintingTool > m_muonPrinter
Definition TrigEDMChecker.h:130
TrigEDMChecker::m_trigDec
PublicToolHandle< Trig::TrigDecisionTool > m_trigDec
Definition TrigEDMChecker.h:137
TrigEDMChecker::m_doDumpxAODPhotonContainer
Gaudi::Property< bool > m_doDumpxAODPhotonContainer
Definition TrigEDMChecker.h:70
TrigEDMChecker::m_doDumpxAODMuonContainer
Gaudi::Property< bool > m_doDumpxAODMuonContainer
Definition TrigEDMChecker.h:58
TrigEDMChecker::dumpTrigComposite
StatusCode dumpTrigComposite()
Dump information on TrigComposite collections.
Definition TrigEDMChecker.cxx:2048
TrigEDMChecker::m_doDumpxAODJetContainer
Gaudi::Property< bool > m_doDumpxAODJetContainer
Definition TrigEDMChecker.h:43
TrigEDMChecker::dumpxAODTrigMinBias
StatusCode dumpxAODTrigMinBias()
Definition TrigEDMChecker.cxx:513
TrigEDMChecker::dumpTrigL2BphysContainer
StatusCode dumpTrigL2BphysContainer()
Definition TrigEDMChecker.cxx:1111
TrigEDMChecker::m_dumpTrigCompositeContainers
Gaudi::Property< std::vector< std::string > > m_dumpTrigCompositeContainers
Definition TrigEDMChecker.h:96
TrigEDMChecker::m_doDumpLVL1_ROI
Gaudi::Property< bool > m_doDumpLVL1_ROI
Definition TrigEDMChecker.h:37
TrigEDMChecker::TrigCompositeNavigationToDot
StatusCode TrigCompositeNavigationToDot(std::string &returnValue, bool &pass)
Construct graph of HLT navigation in Run-3.
Definition TrigEDMChecker.cxx:2151
TrigEDMChecker::checkTrigCompositeElementLink
StatusCode checkTrigCompositeElementLink(const xAOD::TrigComposite *tc, size_t element)
Dump details on element links within TrigComposites.
Definition TrigEDMChecker.cxx:2086
TrigEDMChecker::dumpxAODTauJetContainer
StatusCode dumpxAODTauJetContainer()
Definition TrigEDMChecker.cxx:1646
TrigEDMChecker::m_doDumpTrigPassBits
Gaudi::Property< bool > m_doDumpTrigPassBits
Definition TrigEDMChecker.h:34
TrigEDMChecker::m_doDumpTrigEFBphysContainer
Gaudi::Property< bool > m_doDumpTrigEFBphysContainer
Definition TrigEDMChecker.h:49
TrigEDMChecker::m_navigationHandleKey
SG::ReadHandleKey< xAOD::TrigNavigation > m_navigationHandleKey
Definition TrigEDMChecker.h:134
TrigEDMChecker::dumpxAODTrigEMCluster
StatusCode dumpxAODTrigEMCluster()
Definition TrigEDMChecker.cxx:1555
TrigEDMChecker::dumpLVL1_ROI
StatusCode dumpLVL1_ROI()
Definition TrigEDMChecker.cxx:599
TrigEDMChecker::dumpxAODTrigMissingET
StatusCode dumpxAODTrigMissingET()
Definition TrigEDMChecker.cxx:527
TrigEDMChecker::m_navigationTool
ToolHandle< HLT::Navigation > m_navigationTool
Definition TrigEDMChecker.h:136
TrigEDMChecker::m_doDumpxAODVertex
Gaudi::Property< bool > m_doDumpxAODVertex
Definition TrigEDMChecker.h:79
TrigEDMChecker::m_doDumpxAODTrackParticle
Gaudi::Property< bool > m_doDumpxAODTrackParticle
Definition TrigEDMChecker.h:76
TrigEDMChecker::dumpxAODElectronContainer
StatusCode dumpxAODElectronContainer()
Definition TrigEDMChecker.cxx:764
TrigEDMChecker::dumpxAODPhotonContainer
StatusCode dumpxAODPhotonContainer()
Definition TrigEDMChecker.cxx:920
TrigEDMChecker::m_doDumpTrigL2BphysContainer
Gaudi::Property< bool > m_doDumpTrigL2BphysContainer
Definition TrigEDMChecker.h:46
TrigEDMChecker::m_doDumpxAODTrigElectronContainer
Gaudi::Property< bool > m_doDumpxAODTrigElectronContainer
Definition TrigEDMChecker.h:61
TrigEDMChecker::m_clidSvc
ServiceHandle< IClassIDSvc > m_clidSvc
Definition TrigEDMChecker.h:131
TrigEDMChecker::dumpxAODVertex
StatusCode dumpxAODVertex()
Definition TrigEDMChecker.cxx:1918
TrigEDMChecker::m_doDumpAllTrigComposite
Gaudi::Property< bool > m_doDumpAllTrigComposite
Definition TrigEDMChecker.h:95
TrigEDMChecker::dumpxAODTrigPhotonContainer
StatusCode dumpxAODTrigPhotonContainer()
Definition TrigEDMChecker.cxx:721
TrigEDMChecker::dumpTrigSpacePointCounts
void dumpTrigSpacePointCounts()
Definition TrigEDMChecker.cxx:307
TrigEDMChecker::m_doDumpAll
Gaudi::Property< bool > m_doDumpAll
Definition TrigEDMChecker.h:32
TrigEDMChecker::m_doDumpStoreGate
Gaudi::Property< bool > m_doDumpStoreGate
Definition TrigEDMChecker.h:82
TrigEDMChecker::m_doDumpTrigCompsiteNavigation
Gaudi::Property< bool > m_doDumpTrigCompsiteNavigation
Definition TrigEDMChecker.h:119
TrigEDMChecker::dumpTDT
StatusCode dumpTDT(const EventContext &ctx)
Definition TrigEDMChecker.cxx:1954
TrigEDMChecker::m_excludeFailedHypoNodes
Gaudi::Property< bool > m_excludeFailedHypoNodes
Definition TrigEDMChecker.h:102
TrigEDMChecker::m_doDumpxAODTrigEMClusterContainer
Gaudi::Property< bool > m_doDumpxAODTrigEMClusterContainer
Definition TrigEDMChecker.h:55
TrigEDMChecker::dumpxAODMuonContainer
StatusCode dumpxAODMuonContainer()
Definition TrigEDMChecker.cxx:639
TrigEDMChecker::m_doDumpxAODTrigMinBias
Gaudi::Property< bool > m_doDumpxAODTrigMinBias
Definition TrigEDMChecker.h:88
TrigEDMChecker::m_doDumpxAODTauJetContainer
Gaudi::Property< bool > m_doDumpxAODTauJetContainer
Definition TrigEDMChecker.h:73
TrigEDMChecker::dumpxAODTrigElectronContainer
StatusCode dumpxAODTrigElectronContainer()
Definition TrigEDMChecker.cxx:663
TrigEDMChecker::m_muonTracksKey
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_muonTracksKey
Definition TrigEDMChecker.h:133
TrigEDMChecker::m_dumpNavForChain
Gaudi::Property< std::string > m_dumpNavForChain
Definition TrigEDMChecker.h:101
TrigEDMChecker::dumpxAODJetContainer
StatusCode dumpxAODJetContainer()
Definition TrigEDMChecker.cxx:1208
TrigEDMChecker::m_doDumpNavigation
Gaudi::Property< bool > m_doDumpNavigation
Definition TrigEDMChecker.h:98
TrigEDMChecker::m_doDumpxAODTrigPhotonContainer
Gaudi::Property< bool > m_doDumpxAODTrigPhotonContainer
Definition TrigEDMChecker.h:64
TrigEDMChecker::m_doDumpxAODTrigMissingET
Gaudi::Property< bool > m_doDumpxAODTrigMissingET
Definition TrigEDMChecker.h:40
TrigEDMChecker::m_doDumpxAODTrigEMCluster
Gaudi::Property< bool > m_doDumpxAODTrigEMCluster
Definition TrigEDMChecker.h:52
TrigEDMChecker::initialize
virtual StatusCode initialize() override
Adapted from code by A.Hamilton to check trigger EDM; R.Goncalo 21/11/07.
Definition TrigEDMChecker.cxx:70
Trig::ChainGroup::isPassed
bool isPassed(unsigned int condition=TrigDefs::Physics) const
tells if chain group passed
Definition ChainGroup.cxx:208
Trig::ChainGroup::getListOfTriggers
std::vector< std::string > getListOfTriggers() const
Definition ChainGroup.cxx:467
Trk::ParametersBase::momentum
const Amg::Vector3D & momentum() const
Access method for the momentum.
pair
STL class.
xAOD::CaloCluster_v1::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition CaloCluster_v1.cxx:251
xAOD::CaloCluster_v1::e
virtual double e() const
The total energy of the particle.
Definition CaloCluster_v1.cxx:265
xAOD::CaloCluster_v1::phi
virtual double phi() const
The azimuthal angle ( ) of the particle.
Definition CaloCluster_v1.cxx:256
xAOD::CaloCluster_v1::ETACALOFRAME
@ ETACALOFRAME
Eta in the calo frame (for egamma).
Definition CaloCluster_v1.h:190
xAOD::CaloCluster_v1::PHICALOFRAME
@ PHICALOFRAME
Phi in the calo frame (for egamma).
Definition CaloCluster_v1.h:191
xAOD::JetConstituentVector
A vector of jet constituents at the scale used during jet finding.
Definition JetConstituentVector.h:117
xAOD::JetConstituentVector::begin
iterator begin() const
iterator on the first constituent
Definition JetConstituentVector.cxx:103
xAOD::JetConstituentVector::size
size_t size() const
number of constituents
Definition JetConstituentVector.cxx:102
xAOD::JetConstituentVector::end
iterator end() const
iterator after the last constituent
Definition JetConstituentVector.cxx:104
xAOD::Jet_v1::getConstituents
JetConstituentVector getConstituents() const
Return a vector of consituents. The object behaves like vector<const IParticle*>. See JetConstituentV...
Definition Jet_v1.cxx:149
xAOD::TrigBphys_v1::roiId
uint32_t roiId() const
accessor method: ID of L1 RoI
xAOD::TrigBphys_v1::eta
float eta() const
accessor method: eta
xAOD::TrigBphys_v1::particleType
pType particleType() const
accessor method: particle Type
xAOD::TrigBphys_v1::phi
float phi() const
accessor method: phi
xAOD::TrigBphys_v1::fitmass
float fitmass() const
accessor method: mass from vertex fit
xAOD::TrigBphys_v1::mass
float mass() const
accessor method: mass
xAOD::TrigEMCluster_v1::eta
float eta() const
get Eta (calibrated)
xAOD::TrigEMCluster_v1::phi
float phi() const
get Phi (calibrated)
xAOD::TrigEMCluster_v1::e237
float e237() const
get Energy in a 3x7 cluster (no calibration) around hottest cell
xAOD::TrigEMCluster_v1::eta1
float eta1() const
get Eta sampling 1 (strip layer)
xAOD::TrigEMCluster_v1::e277
float e277() const
get Energy in a 7x7 cluster (no calibration) around hottest cell
xAOD::TrigEMCluster_v1::ehad1
float ehad1() const
get hadronic Energy (first hadronic layer)
xAOD::TrigEMCluster_v1::weta2
float weta2() const
get cluster width (based on a 3x5 cluster - 2nd layer)
xAOD::TrigEMCluster_v1::energy
float energy() const
get Energy (calibrated)
xAOD::TrigEMCluster_v1::fracs1
float fracs1() const
get Energy in a 7 strips (around hottest strip) minus energy in 3 strips divided by energy in 3 strip...
xAOD::TrigMissingET_v1::flag
int flag() const
Get the status flag.
xAOD::TrigMissingET_v1::ex
float ex() const
Get the x component of the missing energy.
xAOD::TrigMissingET_v1::exComponent
float exComponent(unsigned int index) const
Get the x component fo the energy of a component.
Definition TrigMissingET_v1.cxx:119
xAOD::TrigMissingET_v1::nameOfComponent
const std::string & nameOfComponent(unsigned int index) const
Get the name of a component.
Definition TrigMissingET_v1.cxx:108
xAOD::TrigMissingET_v1::eyComponent
float eyComponent(unsigned int index) const
Get the y component fo the energy of a component.
Definition TrigMissingET_v1.cxx:130
xAOD::TrigMissingET_v1::roiWord
uint32_t roiWord() const
Get the ROI word.
xAOD::TrigMissingET_v1::sumE
float sumE() const
Get the sum of the E of the missing energy.
xAOD::TrigMissingET_v1::calib0Component
float calib0Component(unsigned int index) const
Get Calib0.
Definition TrigMissingET_v1.cxx:174
xAOD::TrigMissingET_v1::calib1Component
float calib1Component(unsigned int index) const
Get Calib1.
Definition TrigMissingET_v1.cxx:185
xAOD::TrigMissingET_v1::getNumberOfComponents
unsigned int getNumberOfComponents() const
Get the number of components.
Definition TrigMissingET_v1.cxx:103
xAOD::TrigMissingET_v1::ezComponent
float ezComponent(unsigned int index) const
Get the z component fo the energy of a component.
Definition TrigMissingET_v1.cxx:141
xAOD::TrigMissingET_v1::sumEt
float sumEt() const
Get the sum of the ET of the missing energy.
xAOD::TrigMissingET_v1::statusComponent
short statusComponent(unsigned int index) const
Get Status.
Definition TrigMissingET_v1.cxx:207
xAOD::TrigMissingET_v1::ez
float ez() const
Get the z component of the missing energy.
xAOD::TrigMissingET_v1::sumOfSignsComponent
float sumOfSignsComponent(unsigned int index) const
Get SumOfSigns.
Definition TrigMissingET_v1.cxx:196
xAOD::TrigMissingET_v1::sumEComponent
float sumEComponent(unsigned int index) const
Get SumE.
Definition TrigMissingET_v1.cxx:163
xAOD::TrigMissingET_v1::ey
float ey() const
Get the y component of the missing energy.
xAOD::TrigMissingET_v1::sumEtComponent
float sumEtComponent(unsigned int index) const
Get SumEt.
Definition TrigMissingET_v1.cxx:152
xAOD::TrigMissingET_v1::usedChannelsComponent
unsigned short usedChannelsComponent(unsigned int index) const
Get used Channels.
Definition TrigMissingET_v1.cxx:218
xAOD::TrigNavigation_v1::serialized
const std::vector< unsigned int > & serialized() const
expose the navigation information (in serialized form)
xAOD::TrigSpacePointCounts_v1::sctSpEndcapA
unsigned int sctSpEndcapA() const
xAOD::TrigSpacePointCounts_v1::pixelClusSizeMin
float pixelClusSizeMin() const
xAOD::TrigSpacePointCounts_v1::pixelClusTotBins
unsigned int pixelClusTotBins() const
xAOD::TrigSpacePointCounts_v1::sctSpBarrel
unsigned int sctSpBarrel() const
xAOD::TrigSpacePointCounts_v1::pixelClusSizeBins
unsigned int pixelClusSizeBins() const
xAOD::TrigSpacePointCounts_v1::contentsPixelClusEndcapC
const std::vector< float > & contentsPixelClusEndcapC() const
getters
xAOD::TrigSpacePointCounts_v1::contentsPixelClusEndcapA
const std::vector< float > & contentsPixelClusEndcapA() const
xAOD::TrigSpacePointCounts_v1::sctSpEndcapC
unsigned int sctSpEndcapC() const
xAOD::TrigSpacePointCounts_v1::pixelClusSizeMax
float pixelClusSizeMax() const
xAOD::TrigSpacePointCounts_v1::pixelClusTotMax
float pixelClusTotMax() const
xAOD::TrigSpacePointCounts_v1::pixelClusTotMin
float pixelClusTotMin() const
xAOD::TrigSpacePointCounts_v1::contentsPixelClusBarrel
const std::vector< float > & contentsPixelClusBarrel() const
xAOD::TrigT2MbtsBits_v1::triggerEnergies
const std::vector< float > & triggerEnergies() const
Return the trigger energies of each counter.
xAOD::TrigT2MbtsBits_v1::triggerTimes
const std::vector< float > & triggerTimes() const
Return the relative times of the triggers.
xAOD::TrigTrackCounts_v1::phiMin
float phiMin() const
xAOD::TrigTrackCounts_v1::eta_phi
const std::vector< float > & eta_phi() const
xAOD::TrigTrackCounts_v1::phiBins
unsigned int phiBins() const
xAOD::TrigTrackCounts_v1::etaMin
float etaMin() const
xAOD::TrigTrackCounts_v1::etaBins
unsigned int etaBins() const
xAOD::TrigTrackCounts_v1::z0_pt
const std::vector< float > & z0_pt() const
getters
xAOD::TrigTrackCounts_v1::ptMin
float ptMin() const
xAOD::TrigTrackCounts_v1::phiMax
float phiMax() const
xAOD::TrigTrackCounts_v1::etaMax
float etaMax() const
xAOD::TrigTrackCounts_v1::z0Max
float z0Max() const
xAOD::TrigTrackCounts_v1::z0Min
float z0Min() const
xAOD::TrigTrackCounts_v1::ptMax
float ptMax() const
xAOD::TrigTrackCounts_v1::z0Bins
unsigned int z0Bins() const
xAOD::TrigTrackCounts_v1::ptBins
unsigned int ptBins() const
xAOD::TrigVertexCounts_v1::vtxTrkPtSqSum
const std::vector< float > & vtxTrkPtSqSum() const
xAOD::TrigVertexCounts_v1::vtxNtrks
const std::vector< unsigned int > & vtxNtrks() const
getters
contains
bool contains(const std::string &s, const std::string &regx)
does a string contain the substring
Definition hcg.cxx:116
count
int count(std::string s, const std::string &regx)
count how many occurances of a regx are in a string
Definition hcg.cxx:148
SG::sgkey_t
uint32_t sgkey_t
Type used for hashed StoreGate key+CLID pairs.
Definition sgkey_t.h:32
TrigCompositeUtils::createLegName
HLT::Identifier createLegName(const HLT::Identifier &chainIdentifier, size_t counter)
Generate the HLT::Identifier which corresponds to a specific leg of a given chain.
Definition TrigCompositeUtilsRoot.cxx:166
TrigCompositeUtils::UNSET
@ UNSET
Default property of state.
Definition ActiveState.h:19
TrigCompositeUtils::INACTIVE
@ INACTIVE
The link was inactive for all of the HLT Chains requested in the TDT.
Definition ActiveState.h:21
TrigCompositeUtils::newDecisionIn
Decision * newDecisionIn(DecisionContainer *dc, const std::string &name)
Helper method to create a Decision object, place it in the container and return a pointer to it.
Definition TrigCompositeUtilsRoot.cxx:44
TrigCompositeUtils::getIndexFromLeg
int32_t getIndexFromLeg(const HLT::Identifier &legIdentifier)
Extract the numeric index of a leg identifier.
Definition TrigCompositeUtilsRoot.cxx:191
TrigCompositeUtils::createAndStore
SG::WriteHandle< DecisionContainer > createAndStore(const SG::WriteHandleKey< DecisionContainer > &key, const EventContext &ctx)
Creates and right away records the DecisionContainer with the key.
Definition TrigCompositeUtilsRoot.cxx:28
TrigCompositeUtils::linkToPrevious
void linkToPrevious(Decision *d, const std::string &previousCollectionKey, size_t previousIndex)
Links to the previous object, location of previous 'seed' decision supplied by hand.
Definition TrigCompositeUtilsRoot.cxx:140
TrigCompositeUtils::addDecisionID
void addDecisionID(DecisionID id, Decision *d)
Appends the decision (given as ID) to the decision object.
Definition TrigCompositeUtilsRoot.cxx:59
TrigCompositeUtils::isLegId
bool isLegId(const HLT::Identifier &legIdentifier)
Recognise whether the chain ID is a leg ID.
Definition TrigCompositeUtilsRoot.cxx:219
TrigDefs::includeFailedDecisions
static const unsigned int includeFailedDecisions
Run3 synonym of alsoDeactivateTEs.
Definition TrigEvent/TrigDecisionInterface/TrigDecisionInterface/Conditions.h:40
TrigDefs::allFeaturesOfType
static const unsigned int allFeaturesOfType
Run 3 "enum". Return all features along legs (still with type and container checks).
Definition TrigEvent/TrigDecisionInterface/TrigDecisionInterface/Conditions.h:70
Trk::Perigee
ParametersT< TrackParametersDim, Charged, PerigeeSurface > Perigee
Definition Tracking/TrkEvent/TrkParameters/TrkParameters/TrackParameters.h:33
Trk::theta
@ theta
Definition ParamDefs.h:66
Trk::phi
@ phi
Definition ParamDefs.h:75
Trk::px
@ px
Definition ParamDefs.h:59
Trk::py
@ py
Definition ParamDefs.h:60
dump
-event-from-file
index
Definition index.py:1
xAOD::EgammaParameters::deltaPhi2
@ deltaPhi2
difference between the cluster phi (second sampling) and the phi of the track extrapolated to the sec...
Definition EgammaEnums.h:205
xAOD::EgammaParameters::deltaEta1
@ deltaEta1
difference between the cluster eta (first sampling) and the eta of the track extrapolated to the firs...
Definition EgammaEnums.h:185
xAOD::EgammaParameters::e237
@ e237
uncalibrated energy (sum of cells) of the middle sampling in a rectangle of size 3x7
Definition EgammaEnums.h:78
xAOD::EgammaParameters::ethad
@ ethad
ET leakage into hadronic calorimeter with exclusion of energy in CaloSampling::TileGap3.
Definition EgammaEnums.h:46
xAOD::EgammaParameters::e2ts1
@ e2ts1
2nd max in strips calc by summing 3 strips
Definition EgammaEnums.h:106
xAOD::EgammaParameters::e335
@ e335
uncalibrated energy (sum of cells) of the third sampling in a rectangle of size 3x5
Definition EgammaEnums.h:87
xAOD::EgammaParameters::e132
@ e132
uncalibrated energy (sum of cells) in strips in a 3x2 window in cells in eta X phi
Definition EgammaEnums.h:37
xAOD::EgammaParameters::e011
@ e011
uncalibrated energy (sum of cells) in presampler in a 1x1 window in cells in eta X phi
Definition EgammaEnums.h:31
xAOD::EgammaParameters::e2tsts1
@ e2tsts1
energy of the cell corresponding to second energy maximum in the first sampling
Definition EgammaEnums.h:109
xAOD::Iso::etcone20
@ etcone20
Calorimeter isolation.
Definition IsolationType.h:32
xAOD::Iso::ptcone20
@ ptcone20
Track isolation.
Definition IsolationType.h:40
xAOD::TauJetParameters::hadRadius
@ hadRadius
Get hadron calorimeter radius.
Definition TauDefs.h:192
xAOD::TauJetParameters::EMRadius
@ EMRadius
Get E_T radius.
Definition TauDefs.h:190
xAOD::TauJetParameters::etHadAtEMScale
@ etHadAtEMScale
Get Hadronic energy at EM scale.
Definition TauDefs.h:196
xAOD::TauJetParameters::isolFrac
@ isolFrac
Get isolation fraction.
Definition TauDefs.h:198
xAOD::TauJetParameters::trkAvgDist
@ trkAvgDist
Get calibrated EM transverse energy (DEPRECATED since r19).
Definition TauDefs.h:214
xAOD::TauJetParameters::etEMAtEMScale
@ etEMAtEMScale
Get EM energy at EM scale.
Definition TauDefs.h:194
xAOD::TauJetParameters::centFrac
@ centFrac
Get centrality fraction.
Definition TauDefs.h:200
xAOD::TauJetParameters::dRmax
@ dRmax
Get maximal dR of tracks associated to calo-seeded tau.
Definition TauDefs.h:226
xAOD::Jet
Jet_v1 Jet
Definition of the current "jet version".
Definition Event/xAOD/xAODJet/xAODJet/Jet.h:17
xAOD::TrigVertexCountsContainer
TrigVertexCountsContainer_v1 TrigVertexCountsContainer
Definition TrigVertexCountsContainer.h:11
xAOD::PhotonContainer
PhotonContainer_v1 PhotonContainer
Definition of the current "photon container version".
Definition Event/xAOD/xAODEgamma/xAODEgamma/PhotonContainer.h:17
xAOD::TrigPassBitsContainer
TrigPassBitsContainer_v1 TrigPassBitsContainer
Define the latest version of the trig pass bits container class.
Definition TrigPassBitsContainer.h:17
xAOD::TrigMissingETContainer
TrigMissingETContainer_v1 TrigMissingETContainer
DataVector of TrigMissingET - the current version.
Definition Event/xAOD/xAODTrigMissingET/xAODTrigMissingET/TrigMissingETContainer.h:13
xAOD::TrigCompositeContainer
TrigCompositeContainer_v1 TrigCompositeContainer
Declare the latest version of the container.
Definition TrigCompositeContainer.h:17
xAOD::TrigElectronContainer
TrigElectronContainer_v1 TrigElectronContainer
Declare the latest version of the container.
Definition Event/xAOD/xAODTrigEgamma/xAODTrigEgamma/TrigElectronContainer.h:17
xAOD::ElectronContainer
ElectronContainer_v1 ElectronContainer
Definition of the current "electron container version".
Definition Event/xAOD/xAODEgamma/xAODEgamma/ElectronContainer.h:17
xAOD::TrigBphys
TrigBphys_v1 TrigBphys
Definition TrigBphys.h:18
xAOD::TriggerMenuContainer
TriggerMenuContainer_v1 TriggerMenuContainer
Define the latest version of the trigger menu container class.
Definition TriggerMenuContainer.h:17
xAOD::CaloCluster
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Definition Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloCluster.h:19
xAOD::TrigComposite
TrigComposite_v1 TrigComposite
Declare the latest version of the class.
Definition TrigComposite.h:16
xAOD::TrigSpacePointCountsContainer
TrigSpacePointCountsContainer_v1 TrigSpacePointCountsContainer
Definition TrigSpacePointCountsContainer.h:11
xAOD::TrigEMClusterContainer
TrigEMClusterContainer_v1 TrigEMClusterContainer
Define the latest version of the trigger EM cluster container.
Definition Event/xAOD/xAODTrigCalo/xAODTrigCalo/TrigEMClusterContainer.h:17
xAOD::TrigTrackCountsContainer
TrigTrackCountsContainer_v1 TrigTrackCountsContainer
Definition TrigTrackCountsContainer.h:11
xAOD::TrigPhotonContainer
TrigPhotonContainer_v1 TrigPhotonContainer
Declare the latest version of the container.
Definition Event/xAOD/xAODTrigEgamma/xAODTrigEgamma/TrigPhotonContainer.h:17
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::VertexContainer
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".
Definition VertexContainer.h:14
xAOD::TrigPassBits
TrigPassBits_v1 TrigPassBits
Define the latest version of the trigger pass bits class.
Definition Event/xAOD/xAODTrigger/xAODTrigger/TrigPassBits.h:16
xAOD::TrackParticleContainer
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticleContainer.h:14
xAOD::TrigEMCluster
TrigEMCluster_v1 TrigEMCluster
Define the latest version of the trigger EM cluster class.
Definition Event/xAOD/xAODTrigCalo/xAODTrigCalo/TrigEMCluster.h:17
xAOD::TrigNavigation
TrigNavigation_v1 TrigNavigation
Define the latest version of the trigger navigation class.
Definition TrigNavigation.h:16
xAOD::JetContainer
JetContainer_v1 JetContainer
Definition of the current "jet container version".
Definition JetContainer.h:17
xAOD::TauJetContainer
TauJetContainer_v3 TauJetContainer
Definition of the current "taujet container version".
Definition TauJetContainer.h:17
xAOD::MuonContainer
MuonContainer_v1 MuonContainer
Definition of the current "Muon container version".
Definition Event/xAOD/xAODMuon/xAODMuon/MuonContainer.h:14
xAOD::TrigBphysContainer
TrigBphysContainer_v1 TrigBphysContainer
Definition TrigBphysContainer.h:12
xAOD::numberOfPixelHoles
@ numberOfPixelHoles
number of pixel layers on track with absence of hits [unit8_t].
Definition TrackingPrimitives.h:262
xAOD::numberOfTRTHits
@ numberOfTRTHits
number of TRT hits [unit8_t].
Definition TrackingPrimitives.h:276
xAOD::numberOfTRTHoles
@ numberOfTRTHoles
number of TRT holes [unit8_t].
Definition TrackingPrimitives.h:278
xAOD::numberOfBLayerHits
@ numberOfBLayerHits
these are the hits in the first pixel layer, i.e.
Definition TrackingPrimitives.h:232
xAOD::numberOfSCTHits
@ numberOfSCTHits
number of hits in SCT [unit8_t].
Definition TrackingPrimitives.h:269
xAOD::numberOfPixelHits
@ numberOfPixelHits
these are the pixel hits, including the b-layer [unit8_t].
Definition TrackingPrimitives.h:260
xAOD::numberOfSCTHoles
@ numberOfSCTHoles
number of SCT holes [unit8_t].
Definition TrackingPrimitives.h:271
xAOD::IParticleContainer
DataVector< IParticle > IParticleContainer
Simple convenience declaration of IParticleContainer.
Definition xAOD/xAODBase/xAODBase/IParticleContainer.h:32
ClassID_traits
Default, invalid implementation of ClassID_traits.
Definition Control/AthenaKernel/AthenaKernel/ClassID_traits.h:37
ClassID_traits::ID
static constexpr CLID ID()
Definition Control/AthenaKernel/AthenaKernel/ClassID_traits.h:44
LinkInfo
Helper to keep a Decision object, ElementLink and ActiveState (with respect to some requested ChainGr...
Definition LinkInfo.h:22
msg
MsgStream & msg
Definition testRead.cxx:32
Generated on for ATLAS Offline Software by doxygen 1.17.0