7 #include "GaudiKernel/IToolSvc.h"
109 static const int maxRepWarnings = 5;
114 m_clidSvc(
"ClassIDSvc",
name )
239 return StatusCode::SUCCESS;
263 if(m_doDumpTrackParticleContainer){
265 if (
sc.isFailure()) {
266 ATH_MSG_ERROR(
"The method dumpTrackParticleContainer() failed");
270 if(m_doDumpAll || m_doDumpLVL1_ROI ){
272 if (
sc.isFailure()) {
287 if(m_doDumpAll || m_doDumpxAODTrigMissingET){
289 if (
sc.isFailure()) {
294 if(m_doDumpAll || m_doDumpMuonFeature){
296 if (
sc.isFailure()) {
301 if(m_doDumpAll || m_doDumpCombinedMuonFeature){
303 if (
sc.isFailure()) {
304 ATH_MSG_ERROR(
"The method dumpCombinedMuonFeature() failed");
306 sc = dumpCombinedMuonFeatureContainer();
307 if (
sc.isFailure()) {
308 ATH_MSG_ERROR(
"The method dumpCombinedMuonFeatureContainer() failed");
309 return StatusCode::SUCCESS;
313 if(m_doDumpAll || m_doDumpTileMuFeature) {
315 if (
sc.isFailure()) {
316 ATH_MSG_ERROR(
"The method dumpTileMuFeatureContainer() failed");
320 if(m_doDumpAll || m_doDumpTileTrackMuFeature) {
322 if (
sc.isFailure()) {
323 ATH_MSG_ERROR(
"The method dumpTileTrackMuFeatureContainer() failed");
327 if(m_doDumpAll || m_doDumpTrigEMCluster){
329 if (
sc.isFailure()) {
334 if(m_doDumpAll || m_doDumpTrigEMClusterContainer){
336 if (
sc.isFailure()) {
337 ATH_MSG_ERROR(
"The method dumpTrigEMClusterContainer() failed");
341 if(m_doDumpAll || m_doDumpxAODTrigEMCluster){
343 if (
sc.isFailure()) {
348 if(m_doDumpAll || m_doDumpxAODTrigEMClusterContainer){
350 if (
sc.isFailure()) {
351 ATH_MSG_ERROR(
"The method dumpxAODTrigEMClusterContainer() failed");
355 if(m_doDumpTrigTauClusterContainer){
357 if (
sc.isFailure()) {
358 ATH_MSG_ERROR(
"The method dumpTrigTauClusterContainer() failed");
362 if(m_doDumpAll || m_doDumpTrigPhotonContainer){
364 if (
sc.isFailure()) {
365 ATH_MSG_ERROR(
"The method dumpTrigPhotonContainer() failed");
369 if(m_doDumpAll || m_doDumpxAODJetContainer){
371 if (
sc.isFailure()) {
376 if(m_doDumpAll || m_doDumpTrigL2BphysContainer){
378 if (
sc.isFailure()) {
379 ATH_MSG_ERROR(
"The method dumpTrigL2BphysContainer() failed");
383 if(m_doDumpAll || m_doDumpTrigEFBphysContainer){
385 if (
sc.isFailure()) {
386 ATH_MSG_ERROR(
"The method dumpTrigEFBphysContainer() failed");
390 if(m_doDumpAll || m_doDumpTrigEFBjetContainer){
392 if (
sc.isFailure()) {
393 ATH_MSG_ERROR(
"The method dumpTrigEFBjetContainer() failed");
397 if(m_doDumpAll || m_doDumpTrigL2BjetContainer){
399 if (
sc.isFailure()) {
400 ATH_MSG_ERROR(
"The method dumpTrigL2BjetContainer() failed");
404 if(m_doDumpAll || m_doDumpTrigMuonEFContainer){
406 if (
sc.isFailure()) {
407 ATH_MSG_ERROR(
"The method dumpTrigMuonEFContainer() failed");
411 if(m_doDumpAll || m_doDumpTrigMuonEFInfoContainer){
413 if (
sc.isFailure()) {
414 ATH_MSG_ERROR(
"The method dumpTrigMuonEFInfoContainer() failed");
418 if(m_doDumpAll || m_doDumpTrigMuonEFIsolationContainer) {
421 ATH_MSG_ERROR(
"The method dumpTrigMuonEFIsolationContainer() failed");
425 if(m_doDumpAll || m_doDumpxAODMuonContainer) {
432 if(m_doDumpAll || m_doDumpTrigElectronContainer){
434 if (
sc.isFailure()) {
435 ATH_MSG_ERROR(
"The method dumpTrigElectronContainer() failed");
439 if(m_doDumpAll || m_doDumpxAODTrigElectronContainer){
441 if (
sc.isFailure()) {
442 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
446 if(m_doDumpAll || m_doDumpxAODTrigPhotonContainer){
448 if (
sc.isFailure()) {
449 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
453 if(m_doDumpAll || m_doDumpxAODElectronContainer){
455 if (
sc.isFailure()) {
456 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
460 if(m_doDumpAll || m_doDumpxAODPhotonContainer){
462 if (
sc.isFailure()) {
463 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
467 if(m_doDumpTrigTauContainer){
469 if (
sc.isFailure()) {
474 if(m_doDumpTrigTauTracksInfo){
476 if (
sc.isFailure()) {
481 if(m_doDumpAll || m_doDumpHLTResult){
483 if (
sc.isFailure()) {
488 if(m_doDumpTrigInDetTrackCollection){
490 if (
sc.isFailure()) {
491 ATH_MSG_ERROR(
"The method dumpTrigInDetTrackCollection() failed");
495 if(m_doDumpAll || m_doDumpTrigVertexCollection){
497 if (
sc.isFailure()) {
498 ATH_MSG_ERROR(
"The method dumpTrigVertexCollection() failed");
502 if(m_doDumpAll || m_doDumpxAODTauJetContainer){
504 if (
sc.isFailure()) {
505 ATH_MSG_ERROR(
"The method dumpxAODTauJetContainer() failed");
509 if(m_doDumpTauJetContainer){
511 if (
sc.isFailure()) {
516 if(m_doDumpAll || m_doDumpxAODTrackParticle){
518 if (
sc.isFailure()) {
523 if(m_doDumpAll || m_doDumpxAODVertex){
525 if (
sc.isFailure()) {
530 if (m_doDumpAll || m_doDumpxAODTrigMinBias){
532 if (
sc.isFailure()) {
537 if (m_doDumpTrigPassBits){
539 if (
sc.isFailure()) {
544 if (m_doDumpAll || m_doDumpStoreGate) {
548 if (m_doDumpAll || m_doDumpNavigation) {
550 if (
sc.isFailure() ) {
555 if (m_doDumpAll || m_doTDTCheck) {
559 if (m_doDumpAll || m_doDumpAllTrigComposite || m_dumpTrigCompositeContainers.size() > 0) {
563 if (m_doDumpAll || m_doDumpTrigCompsiteNavigation) {
564 std::string trigCompositeSteering;
566 ATH_CHECK(TrigCompositeNavigationToDot(trigCompositeSteering, pass));
567 const EventContext& context = Gaudi::Hive::currentContext();
568 const std::string evtNumber =
std::to_string(context.eventID().event_number());
569 const std::string passStr = (pass ?
"Pass" :
"Fail");
570 std::ofstream
ofile(std::string(
"NavGraph_" + m_dumpNavForChain +
"_Ev" + evtNumber +
"_" + passStr +
".dot").c_str());
571 ofile << trigCompositeSteering;
577 return StatusCode::SUCCESS;
583 const std::string
name=
"HLT_xAOD__TrigPassBitsContainer_passbits";
586 if (
sc.isFailure() ){
591 for(
const auto bits:*xbitscont){
596 ATH_MSG_DEBUG(
"Analyzing bits for " << bits->containerClid() <<
" of size " << bits->size() <<
" with bit size " << bits->passBits().size());
601 for (
int j=0; itr != itrE; ++itr, ++j ) {
610 return StatusCode::SUCCESS;
617 std::string METTag=
"HLT_xAOD__TrigSpacePointCountsContainer_spacepoints";
625 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << SpacePointCountsCont->
size() <<
" elements");
627 std::string
s;
char buff[128];
628 std::vector<float> getVec;
632 for(
uint i = 0;
i < SpacePointCountsCont->
size();
i++) {
633 getVec = SpacePointCountsCont->
at(
i)->contentsPixelClusEndcapC();
635 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
636 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of contentsPixelClusEndcapC() = %10.2f ",
s.c_str(),
sum );
639 getVec = SpacePointCountsCont->
at(
i)->contentsPixelClusBarrel();
641 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
642 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of contentsPixelClusBarrel() = %10.2f ",
s.c_str(),
sum );
645 getVec = SpacePointCountsCont->
at(
i)->contentsPixelClusEndcapA();
647 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
648 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of contentsPixelClusEndcapA() = %10.2f ",
s.c_str(),
sum );
651 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusTotBins() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusTotBins() );
654 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusTotMin() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusTotMin() );
657 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusTotMax() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusTotMax() );
660 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusSizeBins() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusSizeBins() );
663 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusSizeMin() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusSizeMin() );
666 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusSizeMax() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusSizeMax() );
669 snprintf(
buff,
sizeof(
buff),
"REGTEST %s sctSpEndcapC() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->sctSpEndcapC() );
672 snprintf(
buff,
sizeof(
buff),
"REGTEST %s sctSpBarrel() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->sctSpBarrel() );
675 snprintf(
buff,
sizeof(
buff),
"REGTEST %s sctSpEndcapA() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->sctSpEndcapA() );
684 std::string METTag=
"HLT_xAOD__TrigT2MbtsBitsContainer_T2Mbts";
692 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << T2MbtsBitsCont->
size() <<
" elements");
694 std::string
s;
char buff[380];
695 std::vector<float> getVec;
700 getVec = T2MbtsBitsCont->
at(
i)->triggerEnergies();
702 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
703 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of triggerEnergies() = %10.2f ",
s.c_str(),
sum );
706 getVec = T2MbtsBitsCont->
at(
i)->triggerTimes();
708 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
709 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of triggerTimes() = %10.2f ",
s.c_str(),
sum );
718 std::string METTag=
"HLT_xAOD__TrigVertexCountsContainer_vertexcounts";
726 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << T2VertexCountsCont->
size() <<
" elements");
728 std::string
s;
char buff[380];
729 std::vector<float> fgetVec;
731 std::vector<unsigned int> ugetVec;
732 unsigned int usum(0);
735 for(
uint i = 0;
i < T2VertexCountsCont->
size();
i++) {
736 ugetVec = T2VertexCountsCont->
at(
i)->vtxNtrks();
737 for (
uint j = 0; j < ugetVec.size(); ++j) usum += ugetVec[j];
738 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of vtxNtrks() = %u ",
s.c_str(), usum );
741 fgetVec = T2VertexCountsCont->
at(
i)->vtxTrkPtSqSum();
742 for (
uint j = 0; j < fgetVec.size(); ++j) fsum += fgetVec[j];
743 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of vtxTrkPtSqSum() = %10.2f ",
s.c_str(), fsum );
752 std::string METTag=
"HLT_xAOD__TrigTrackCountsContainer_trackcounts";
760 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << T2TrackCountsCont->
size() <<
" elements");
762 std::string
s;
char buff[380];
763 std::vector<float> getVec;
767 for(
uint i = 0;
i < T2TrackCountsCont->
size();
i++) {
768 getVec = T2TrackCountsCont->
at(
i)->z0_pt();
770 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
771 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of z0_pt = %10.2f ",
s.c_str(),
sum );
774 getVec = T2TrackCountsCont->
at(
i)->eta_phi();
776 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
777 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of eta_phi() = %10.2f ",
s.c_str(),
sum );
780 snprintf(
buff,
sizeof(
buff),
"REGTEST %s z0Bins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->z0Bins() );
783 snprintf(
buff,
sizeof(
buff),
"REGTEST %s z0Min() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->z0Min() );
786 snprintf(
buff,
sizeof(
buff),
"REGTEST %s z0Max() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->z0Max() );
789 snprintf(
buff,
sizeof(
buff),
"REGTEST %s ptBins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->ptBins() );
792 snprintf(
buff,
sizeof(
buff),
"REGTEST %s ptMin() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->ptMin() );
795 snprintf(
buff,
sizeof(
buff),
"REGTEST %s ptMax() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->ptMax() );
798 snprintf(
buff,
sizeof(
buff),
"REGTEST %s etaBins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->etaBins() );
801 snprintf(
buff,
sizeof(
buff),
"REGTEST %s etaMin() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->etaMin() );
804 snprintf(
buff,
sizeof(
buff),
"REGTEST %s etaMax() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->etaMax() );
807 snprintf(
buff,
sizeof(
buff),
"REGTEST %s phiBins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->phiBins() );
810 snprintf(
buff,
sizeof(
buff),
"REGTEST %s phiMin() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->phiMin() );
813 snprintf(
buff,
sizeof(
buff),
"REGTEST %s phiMax() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->phiMax() );
826 return StatusCode::SUCCESS;
838 std::string METTags[]={
"HLT_xAOD__TrigMissingETContainer_EFJetEtSum",
"HLT_xAOD__TrigMissingETContainer_TrigEFMissingET",
"HLT_xAOD__TrigMissingETContainer_TrigL2MissingET_FEB",
"HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl"};
840 for(
int itag=0; itag <ntag; itag++) {
847 ATH_MSG_INFO(
"Accessing " << METTags[itag] <<
" with " << MissingETCont->
size() <<
" elements");
852 std::string
s;
char buff[3000];
854 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ex = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->ex() );
856 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ey = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->ey() );
858 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ez = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->ez() );
860 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumET = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->sumEt() );
862 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumE = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->sumE() );
864 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Flag = %d",
s.c_str(), MissingETCont->
at(
i)->flag() );
866 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Flag = %d",
s.c_str(), MissingETCont->
at(
i)->roiWord() );
869 unsigned int Nc = MissingETCont->
at(
i)->getNumberOfComponents();
871 s=
"REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
872 s+=
"/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
876 for(
uint j = 0; j < Nc; j++) {
878 std::string
name = MissingETCont->
at(
i)->nameOfComponent(j);
879 const short status = MissingETCont->
at(
i)->statusComponent(j);
880 const unsigned short usedChan = MissingETCont->
at(
i)->usedChannelsComponent(j);
881 const short sumOfSigns = MissingETCont->
at(
i)->sumOfSignsComponent(j);
882 const float calib0 = MissingETCont->
at(
i)->calib0Component(j);
883 const float calib1 = MissingETCont->
at(
i)->calib1Component(j);
884 const float ex = MissingETCont->
at(
i)->exComponent(j);
885 const float ey = MissingETCont->
at(
i)->eyComponent(j);
886 const float ez = MissingETCont->
at(
i)->ezComponent(j);
887 const float sumE = MissingETCont->
at(
i)->sumEComponent(j);
888 const float sumEt = MissingETCont->
at(
i)->sumEtComponent(j);
891 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
892 name.c_str(),
status, usedChan, sumOfSigns, calib1, calib0,
893 ex, ey, ez, sumE, sumEt);
900 return StatusCode::SUCCESS;
912 std::string METTags[]={
"HLT_TrigMissingETContainer_TrigEFMissingET",
"HLT_TrigMissingETContainer_TrigEFMissingET_FEB",
"HLT_TrigMissingETContainer_TrigEFMissingET_topocl"};
915 for (
int itag=0; itag < ntag; itag++) {
917 StatusCode sc=evtStore()->retrieve(trigMETcont , METTags[itag]);
918 if(
sc.isFailure() ){
919 ATH_MSG_INFO(
"Failed to retrieve TrigMissingETContainer with key " << METTags[itag]);
923 ATH_MSG_INFO(
"Got TrigMissingETContainer with key \"" << METTags[itag]<<
"\"");
928 for (
int j=0; trigMETfirst != trigMETlast; ++trigMETfirst++, ++j ) {
930 ATH_MSG_INFO(
"REGTEST ==========START of TrigMissingET DUMP===========");
935 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ex = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->ex() );
937 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ey = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->ey() );
939 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ez = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->ez() );
941 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Et = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->et() );
943 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumEt = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->sumEt() );
945 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumE = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->sumE() );
947 snprintf(
buff,
sizeof(
buff),
"REGTEST %s E = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->e() );
949 snprintf(
buff,
sizeof(
buff),
"REGTEST %s flag = %10d",
s.c_str(), (*trigMETfirst)->getFlag() );
951 snprintf(
buff,
sizeof(
buff),
"REGTEST %s RoIword = %10ld",
s.c_str(), (*trigMETfirst)->RoIword() );
954 unsigned int Nc = (*trigMETfirst)->getNumOfComponents();
956 s=
"REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
957 s+=
"/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
960 for (
unsigned int i=0;
i<Nc; ++
i) {
961 std::string
name = (*trigMETfirst)->getNameOfComponent(
i);
962 const short status = (*trigMETfirst)->getStatus(
i);
963 const unsigned short usedChan = (*trigMETfirst)->getUsedChannels(
i);
964 const short sumOfSigns = (*trigMETfirst)->getSumOfSigns(
i);
965 const float calib0 = (*trigMETfirst)->getComponentCalib0(
i);
966 const float calib1 = (*trigMETfirst)->getComponentCalib1(
i);
967 const float ex = (*trigMETfirst)->getExComponent(
i);
968 const float ey = (*trigMETfirst)->getEyComponent(
i);
969 const float ez = (*trigMETfirst)->getEzComponent(
i);
970 const float sumE = (*trigMETfirst)->getSumEComponent(
i);
971 const float sumEt = (*trigMETfirst)->getSumEtComponent(
i);
974 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
975 name.c_str(),
status, usedChan, sumOfSigns, calib1, calib0,
976 ex, ey, ez, sumE, sumEt);
986 ATH_MSG_INFO(
"Trying to fetch TrigMissingET objects from older releases");
989 StatusCode sc=evtStore()->retrieve(trigMETfirst ,trigMETlast);
990 if(
sc.isFailure() ){
991 ATH_MSG_INFO(
"Failed to retrieve TrigMissingET (rel. <= 14.2.0)");
994 for( ; trigMETfirst != trigMETlast ; ++trigMETfirst ){
995 const std::string&
name(trigMETfirst.
key());
1001 if(
name.find(
"TrigEFMissingET") != std::string::npos ) {
1003 }
else if(
name.find(
"T2MissingET") != std::string::npos ){
1010 snprintf(
buff,
sizeof(
buff),
"%s Ex = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->ex() );
1012 snprintf(
buff,
sizeof(
buff),
"%s Ey = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->ey() );
1014 snprintf(
buff,
sizeof(
buff),
"%s Ez = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->ez() );
1016 snprintf(
buff,
sizeof(
buff),
"%s Et = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->et() );
1018 snprintf(
buff,
sizeof(
buff),
"%s SumE = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->sumE() );
1020 snprintf(
buff,
sizeof(
buff),
"%s SumEt = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->sumEt() );
1022 snprintf(
buff,
sizeof(
buff),
"%s E = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->e() );
1024 snprintf(
buff,
sizeof(
buff),
"%s flag = %10d",
s.c_str(), trigMETfirst->getFlag() );
1026 snprintf(
buff,
sizeof(
buff),
"%s RoIword = %10ld",
s.c_str(), trigMETfirst->RoIword() );
1029 unsigned int Nc = trigMETfirst->getNumOfComponents();
1031 s=
"REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
1032 s+=
"/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
1035 for (
unsigned int i=0;
i<Nc; ++
i) {
1036 std::string
name = trigMETfirst->getNameOfComponent(
i);
1037 const short status = trigMETfirst->getStatus(
i);
1038 const unsigned short usedChan = trigMETfirst->getUsedChannels(
i);
1039 const short sumOfSigns = trigMETfirst->getSumOfSigns(
i);
1040 const float calib0 = trigMETfirst->getComponentCalib0(
i);
1041 const float calib1 = trigMETfirst->getComponentCalib1(
i);
1042 const float ex = trigMETfirst->getExComponent(
i);
1043 const float ey = trigMETfirst->getEyComponent(
i);
1044 const float ez = trigMETfirst->getEzComponent(
i);
1045 const float sumE = trigMETfirst->getSumEComponent(
i);
1046 const float sumEt = trigMETfirst->getSumEtComponent(
i);
1049 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
1050 name.c_str(),
status, usedChan, sumOfSigns, calib1, calib0,
1051 ex, ey, ez, sumE, sumEt);
1057 ATH_MSG_INFO(
"REGTEST ==========END of TrigMissingET DUMP===========");
1059 return StatusCode::SUCCESS;
1068 ATH_MSG_INFO(
"REGTEST ==========START of TrackParticleContainer DUMP===========");
1070 std::string trackPtags[]={
"HLT_InDetTrigParticleCreation_Bjet_EFID",
1071 "HLT_InDetTrigParticleCreation_Bphysics_EFID",
1072 "HLT_InDetTrigParticleCreation_Electron_EFID",
1073 "HLT_InDetTrigParticleCreation_FullScan_EFID",
1074 "HLT_InDetTrigParticleCreation_Muon_EFID",
1075 "HLT_InDetTrigParticleCreation_Photon_EFID",
1076 "HLT_InDetTrigParticleCreation_Tau_EFID"};
1082 for (
int itag=0; itag<ntag; itag++){
1085 if (
sc.isFailure()) {
1086 ATH_MSG_INFO(
"REGTEST No TrackParticleContainer found with tag " << trackPtags[itag]);
1089 ATH_MSG_INFO(
"TrackParticleContainer found with tag " << trackPtags[itag]
1090 <<
" and size " << pTrackParticleC->
size());
1094 for (
int ind=1; trackItr != trackItrE; ++trackItr,
ind++) {
1097 << trackParticle->
charge() <<
" p "
1098 << trackParticle->
p()<<
" eta " << trackParticle->
eta()
1099 <<
" phi " << trackParticle->
phi());
1107 const auto& parameterVector = perigee->parameters();
1110 " phi " <<parameterVector[
Trk::phi]);
1119 ATH_MSG_WARNING(
" Max attached track warning reached, no further warnings given");
1130 ATH_MSG_INFO(
" vertex position (" << position[0] <<
", " <<
1131 position[1] <<
", " << position[2] <<
") ");
1136 ATH_MSG_WARNING(
" Max attached vertex warning reached, no further warnings given");
1144 const auto&
parameters = perigee->parameters();
1175 ATH_MSG_INFO(
"REGTEST ==========START of LVL1_ROI DUMP===========");
1179 if (
sc.isFailure() ) {
1181 return StatusCode::SUCCESS;
1186 LVL1_ROI::emtaus_type::const_iterator itEMTau =
1188 LVL1_ROI::emtaus_type::const_iterator itEMTau_e =
1191 for( ; itEMTau != itEMTau_e; ++itEMTau, ++j) {
1193 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI Eta is " << itEMTau->getEta());
1194 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI Phi is " << itEMTau->getPhi());
1195 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI EMClus is " << itEMTau->getEMClus());
1196 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI TauClus is " << itEMTau->getTauClus());
1197 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI EMIsol is " << itEMTau->getEMIsol());
1198 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI HadIsol is " << itEMTau->getHadIsol());
1199 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI Core is " << itEMTau->getCore());
1200 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI HadCore is " << itEMTau->getHadCore());
1201 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI roiWord is " << itEMTau->getROIWord());
1204 ATH_MSG_INFO(
"REGTEST ==========END of LVL1_ROI DUMP===========");
1206 return StatusCode::SUCCESS;
1216 ATH_MSG_INFO(
"REGTEST ==========START of TrigPhotonContainer DUMP===========");
1221 StatusCode sc = evtStore()->retrieve(trigPhoton,lastTrigPhoton);
1222 if (
sc.isFailure()) {
1224 return StatusCode::SUCCESS;
1226 ATH_MSG_INFO(
"REGTEST TrigPhotonContainers retrieved");
1229 for (
int i=0; trigPhoton != lastTrigPhoton; ++trigPhoton, ++
i) {
1231 ATH_MSG_INFO(
"REGTEST Looking at TrigPhotonContainer " <<
i);
1236 for (
int j=0; PhotonItr != PhotonItrE; ++PhotonItr, ++j ) {
1240 ATH_MSG_INFO(
"REGTEST TrigPhoton->dPhi() returns " << (*PhotonItr)->dPhi());
1241 ATH_MSG_INFO(
"REGTEST TrigPhoton->dEta() returns " << (*PhotonItr)->dEta());
1242 ATH_MSG_INFO(
"REGTEST TrigPhoton->rCore() returns " << (*PhotonItr)->rCore());
1243 ATH_MSG_INFO(
"REGTEST TrigPhoton->isValid() returns " << (*PhotonItr)->isValid());
1244 ATH_MSG_INFO(
"REGTEST TrigPhoton->Et() returns " << (*PhotonItr)->Et());
1246 if ( (*PhotonItr)->cluster() != NULL ) {
1248 ATH_MSG_INFO(
"REGTEST cluster->e() returns " << (*PhotonItr)->cluster()->e());
1249 ATH_MSG_INFO(
"REGTEST cluster->e237() returns " << (*PhotonItr)->cluster()->e237());
1250 ATH_MSG_INFO(
"REGTEST cluster->e277() returns " << (*PhotonItr)->cluster()->e277());
1251 ATH_MSG_INFO(
"REGTEST cluster->fracs1() returns " << (*PhotonItr)->cluster()->fracs1());
1252 ATH_MSG_INFO(
"REGTEST cluster->weta2() returns " << (*PhotonItr)->cluster()->weta2());
1253 ATH_MSG_INFO(
"REGTEST cluster->ehad() returns " << (*PhotonItr)->cluster()->ehad1());
1254 ATH_MSG_INFO(
"REGTEST cluster->emaxs1() returns " << (*PhotonItr)->cluster()->emaxs1());
1258 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->Pt() returns " << (*PhotonItr)->pt());
1259 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->Eta() returns " << (*PhotonItr)->eta());
1260 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->Phi() returns " << (*PhotonItr)->phi());
1261 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->m() returns " << (*PhotonItr)->m());
1272 for (
int k=0; PhotonItr2 != PhotonItr2E; ++PhotonItr2, ++
k ) {
1274 if ( (**PhotonItr) == (**PhotonItr2) ) {
1275 ATH_MSG_INFO(
"REGTEST TrigPhoton nr. " << j <<
" is the same as TrigPhoton nr. " <<
k);
1278 std::map<std::string, double> v_diff;
1279 diff(*(*PhotonItr),*(*PhotonItr2), v_diff);
1281 ATH_MSG_INFO(
"TrigPhoton " <<
k <<
" different form TrigPhoton " << j <<
" :");
1282 for (
int m=0;
it !=v_diff.end();++
it, ++
m) {
1283 ATH_MSG_INFO(
"TrigPhoton Delta_" << (*it).first <<
" = " << (*it).second);
1290 ATH_MSG_INFO(
"REGTEST ==========END of TrigPhotonContainer DUMP===========");
1293 return StatusCode::SUCCESS;
1303 ATH_MSG_INFO(
"REGTEST ==========START of TrigMuonEFContainer DUMP===========");
1308 StatusCode sc = evtStore()->retrieve(trigMuon,lastTrigMuon);
1309 if (
sc.isFailure()) {
1311 return StatusCode::SUCCESS;
1313 ATH_MSG_INFO(
"REGTEST TrigMuonEFContainers retrieved");
1315 for (
int i=0; trigMuon != lastTrigMuon; ++trigMuon, ++
i) {
1317 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFContainer " <<
i);
1322 for (
int j=0; MuonItr != MuonItrE; ++MuonItr, ++j ) {
1326 ATH_MSG_INFO(
"REGTEST TrigMuonEF->muonCode() returns " << (*MuonItr)->MuonCode());
1327 ATH_MSG_INFO(
"REGTEST TrigMuonEF->roi() returns " << (*MuonItr)->RoINum());
1328 ATH_MSG_INFO(
"REGTEST TrigMuonEF->charge() returns " << (*MuonItr)->Charge());
1330 ATH_MSG_INFO(
"REGTEST Looking at P4IPtCotThPhiM " << j);
1332 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->iPt() returns " << (*MuonItr)->iPt());
1333 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->CotTh() returns " << (*MuonItr)->cotTh());
1334 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->Phi() returns " << (*MuonItr)->phi());
1335 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->m() returns " << (*MuonItr)->m());
1339 ATH_MSG_INFO(
"REGTEST ==========END of TrigMuonEFContainer DUMP===========");
1342 return StatusCode::SUCCESS;
1349 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::MuonContainer DUMP===========" );
1353 if (
sc.isFailure()) {
1354 ATH_MSG_INFO(
"REGTEST No muon container HLT_xAOD__MuonContainer_MuonEFInfo");
1355 return StatusCode::SUCCESS;
1361 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::MuonContainer DUMP===========" );
1363 return StatusCode::SUCCESS;
1373 ATH_MSG_INFO(
"REGTEST ==========START of TrigMuonEFInfoContainer DUMP===========");
1378 StatusCode sc = evtStore()->retrieve(trigMuon,lastTrigMuon);
1379 if (
sc.isFailure()) {
1380 ATH_MSG_INFO(
"REGTEST No TrigMuonEFInfoContainer found");
1381 return StatusCode::SUCCESS;
1383 ATH_MSG_INFO(
"REGTEST TrigMuonEFInfoContainers retrieved");
1386 for (
int i=0; trigMuon != lastTrigMuon; ++trigMuon, ++
i) {
1388 ATH_MSG_INFO(
"REGTEST -+-+-+-+ Looking at TrigMuonEFInfoContainer " <<
i);
1393 for (
int j=0; MuonItr != MuonItrE; ++MuonItr, ++j ) {
1395 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFInfo " << j);
1400 if (*muonInfo == *muonInfo) {
1422 printMuonTrk(muonTrack);
1424 ATH_MSG_INFO(
"REGTEST no SpectrometerTrack track found");
1427 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFTrack ExtrapolatedTrack()");
1430 printMuonTrk(muonTrack);
1432 ATH_MSG_INFO(
"REGTEST no ExtrapolatedTrack track found");
1435 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFTrack CombinedTrack()");
1438 printMuonTrk(muonTrack);
1448 return StatusCode::SUCCESS;
1450 ATH_MSG_INFO(
"REGTEST TrigMuonEFInfo->RoINum() returns " << muonInfo->RoINum());
1454 ATH_MSG_INFO(
"REGTEST ==========END of TrigMuonEFInfoContainer DUMP===========");
1457 return StatusCode::SUCCESS;
1466 ATH_MSG_INFO(
"REGTEST ==========START of TrigMuonEFIsolationContainer DUMP===========");
1471 StatusCode sc = evtStore()->retrieve(trigMuon,lastTrigMuon);
1472 if (
sc.isFailure()) {
1473 ATH_MSG_INFO(
"REGTEST No TrigMuonEFIsolationContainer found");
1474 return StatusCode::SUCCESS;
1476 ATH_MSG_INFO(
"REGTEST TrigMuonEFIsolationContainers retrieved");
1479 for (
int i=0; trigMuon != lastTrigMuon; ++trigMuon, ++
i) {
1481 ATH_MSG_INFO(
"REGTEST -+-+-+-+ Looking at TrigMuonEFIsolationContainer " <<
i);
1486 for (
int j=0; MuonItr != MuonItrE; ++MuonItr, ++j ) {
1488 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFIsolation " << j);
1493 if (*muonIsolation == *muonIsolation) {
1511 ATH_MSG_WARNING(
"REGTEST No InfoTrack attached to this EFIsolation object");
1526 ATH_MSG_INFO(
"REGTEST No Linked MuonEFInfo, expected for pre 2012 data");
1532 ATH_MSG_INFO(
"REGTEST ==========END of TrigMuonEFIsolationContainer DUMP===========");
1533 ATH_MSG_DEBUG(
"dumpTrigMuonEFIsolationContainer() succeeded");
1535 return StatusCode::SUCCESS;
1541 ATH_MSG_INFO(
"POINTER TrigMuonEFTrack: " << muonTrack);
1567 ATH_MSG_INFO(
" POINTER TrigMuonEFCbTrack: " << muonTrack);
1584 ATH_MSG_INFO(
"REGTEST ==========START of TrigElectronContainer DUMP===========");
1589 StatusCode sc = evtStore()->retrieve(trigElec,lastTrigElec);
1590 if (
sc.isFailure()) {
1591 ATH_MSG_INFO(
"REGTEST No TrigelectronContainer found");
1592 return StatusCode::SUCCESS;
1594 ATH_MSG_INFO(
"REGTEST TrigElectronContainers retrieved");
1599 for (
int i=0; trigElec != lastTrigElec; ++trigElec, ++
i) {
1601 ATH_MSG_INFO(
"REGTEST Looking at TrigElectronContainer " <<
i);
1606 for (
int j=0; elecItr != elecItrE; ++elecItr, ++j ) {
1643 for (
int k=0; elecItr2 != elecItr2E; ++elecItr2, ++
k ) {
1645 if ( (**elecItr) == (**elecItr2) ) {
1646 ATH_MSG_INFO(
"REGTEST TrigElectron nr. " << j <<
" is the same as TrigElectron nr. " <<
k);
1649 std::map<std::string, double> v_diff;
1650 diff(*(*elecItr),*(*elecItr2), v_diff);
1652 ATH_MSG_INFO(
"TrigElectron " <<
k <<
" different form TrigElectron " << j <<
" :");
1653 for (
int m=0;
it !=v_diff.end();++
it, ++
m) {
1654 ATH_MSG_INFO(
"TrigElectron Delta_" << (*it).first <<
" = " << (*it).second);
1661 ATH_MSG_INFO(
"REGTEST ==========END of TrigElectronContainer DUMP===========");
1664 return StatusCode::SUCCESS;
1672 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrigElectronContainer DUMP===========" );
1676 if (
sc.isFailure()) {
1677 ATH_MSG_INFO(
"REGTEST No TrigElectron container HLT_xAOD__TrigElectronContainer_L2ElectronFex");
1678 return StatusCode::SUCCESS;
1681 for (
const auto eg : *elCont){
1682 ATH_MSG_INFO(
"REGTEST TrigElectron->Phi() returns " <<
eg->phi());
1683 ATH_MSG_INFO(
"REGTEST TrigElectron->Eta() returns " <<
eg->eta());
1684 ATH_MSG_INFO(
"REGTEST TrigElectron->rEta returns " <<
eg->rcore());
1685 ATH_MSG_INFO(
"REGTEST TrigElectron->eratio() returns " <<
eg->eratio());
1687 ATH_MSG_INFO(
"REGTEST TrigElectron->etHad() returns " <<
eg->etHad());
1689 ATH_MSG_INFO(
"REGTEST TrigElectron caloEta = " <<
eg->caloEta());
1690 ATH_MSG_INFO(
"REGTEST TrigElectron dPhiCalo" <<
eg->trkClusDphi());
1691 ATH_MSG_INFO(
"REGTEST TrigElectron dEtaCalo" <<
eg->trkClusDeta());
1695 ATH_MSG_INFO(
"REGTEST TrigElectron nStrawHits" <<
eg->nTRTHiThresholdHits());
1697 if(
eg->emCluster()){
1698 ATH_MSG_INFO(
"REGTEST TrigElectron EMCluster retrieved");
1699 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster->energy() returns " <<
eg->emCluster()->energy());
1700 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster->phi() returns " <<
eg->emCluster()->phi());
1701 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster->eta() returns " <<
eg->emCluster()->eta());
1702 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster check Element Link");
1703 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster energy = " <<
eg->emCluster()->energy());
1704 ATH_MSG_INFO(
"REGTEST TrigElectron ElementLink emCluster energy = " << (*
eg->emClusterLink())->energy());
1706 else ATH_MSG_INFO(
"REGTEST TrigElectron No EMCluster retrieved!");
1707 ATH_MSG_INFO(
"REGTEST TrigElectron Check TrackParticle");
1708 if(
eg->trackParticle()){
1709 ATH_MSG_INFO(
"REGTEST TrigElectron TrackParticle retrieved");
1710 ATH_MSG_INFO(
"REGTEST TrigElectron trackParticle->pt() returns " <<
eg->trackParticle()->pt());
1711 ATH_MSG_INFO(
"REGTEST TrigElectron trackParticle->phi() returns " <<
eg->trackParticle()->phi());
1712 ATH_MSG_INFO(
"REGTEST TrigElectron trackParticle->eta() returns " <<
eg->trackParticle()->eta());
1713 ATH_MSG_INFO(
"REGTEST TrigElectron check TrackParticle Element Link");
1714 ATH_MSG_INFO(
"REGTEST TrigElectron TrackParticle pt = " <<
eg->trackParticle()->pt());
1715 ATH_MSG_INFO(
"REGTEST TrigElectron ElementLink TrackParticle pt = " << (*
eg->trackParticleLink())->pt());
1717 else ATH_MSG_INFO(
"REGTEST TrigElectron No TrackParticle retrieved!");
1719 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigElectronContainer DUMP===========" );
1721 return StatusCode::SUCCESS;
1730 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrigPhotonContainer DUMP===========" );
1734 if (
sc.isFailure()) {
1735 ATH_MSG_INFO(
"REGTEST No TrigPhoton container HLT_xAOD__TrigPhotonContainer_L2PhotonFex");
1736 return StatusCode::SUCCESS;
1739 for (
const auto eg : *phCont){
1743 ATH_MSG_INFO(
"REGTEST TrigPhoton->dPhi() returns " <<
eg->dPhi());
1744 ATH_MSG_INFO(
"REGTEST TrigPhoton->dEta() returns " <<
eg->dEta());
1745 ATH_MSG_INFO(
"REGTEST TrigPhoton->rEta returns " <<
eg->rcore());
1746 ATH_MSG_INFO(
"REGTEST TrigPhoton->eratio() returns " <<
eg->eratio());
1748 ATH_MSG_INFO(
"REGTEST TrigPhoton->etHad() returns " <<
eg->etHad());
1751 if(
eg->emCluster()){
1752 ATH_MSG_INFO(
"REGTEST TrigPhoton EMCluster retrieved");
1753 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster->energy() returns " <<
eg->emCluster()->energy());
1754 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster->phi() returns " <<
eg->emCluster()->phi());
1755 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster->eta() returns " <<
eg->emCluster()->eta());
1756 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster check Element Link");
1757 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster energy = " <<
eg->emCluster()->energy());
1758 ATH_MSG_INFO(
"REGTEST TrigPhoton ElementLink emCluster energy = " << (*
eg->emClusterLink())->energy());
1760 else ATH_MSG_INFO(
"REGTEST TrigPhoton No EMCluster retrieved!");
1762 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigPhotonContainer DUMP===========" );
1764 return StatusCode::SUCCESS;
1773 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::ElectronContainer DUMP===========" );
1777 if (
sc.isFailure()) {
1778 ATH_MSG_INFO(
"REGTEST No Electron container HLT_xAOD__ElectronContainer_egamma_Electrons");
1779 return StatusCode::SUCCESS;
1781 float val_float=-99;
1782 unsigned int isEMbit=0;
1785 ATH_MSG_INFO(
" REGTEST: xAOD Reconstruction variables: ");
1790 for (
const auto eg : *elCont){
1796 if(
eg->selectionisEM(isEMbit,
"isEMVLoose"))
1797 ATH_MSG_INFO(
" REGTEST: isEMVLoose " << std::hex << isEMbit << std::dec);
1799 if(
eg->selectionisEM(isEMbit,
"isEMLoose"))
1800 ATH_MSG_INFO(
" REGTEST: isEMLoose " << std::hex << isEMbit << std::dec);
1802 if(
eg->selectionisEM(isEMbit,
"isEMMedium"))
1803 ATH_MSG_INFO(
" REGTEST: isEMMedium " << std::hex << isEMbit << std::dec);
1805 if(
eg->selectionisEM(isEMbit,
"isEMTight"))
1806 ATH_MSG_INFO(
" REGTEST: isEMTight " << std::hex << isEMbit << std::dec);
1808 if(
eg->selectionisEM(isEMbit,
"isEMLHVLoose"))
1809 ATH_MSG_INFO(
" REGTEST: isEMLHVLoose " << std::hex << isEMbit << std::dec);
1811 if(
eg->selectionisEM(isEMbit,
"isEMLHLoose"))
1812 ATH_MSG_INFO(
" REGTEST: isEMLHLoose " << std::hex << isEMbit << std::dec);
1814 if(
eg->selectionisEM(isEMbit,
"isEMLHMedium"))
1815 ATH_MSG_INFO(
" REGTEST: isEMLHMedium " << std::hex << isEMbit << std::dec);
1817 if(
eg->selectionisEM(isEMbit,
"isEMLHTight"))
1818 ATH_MSG_INFO(
" REGTEST: isEMLHTight " << std::hex << isEMbit << std::dec);
1828 if(
eg->passSelection(
pid,
"LHVLoose"))
1832 if(
eg->passSelection(
pid,
"LHLoose"))
1836 if(
eg->passSelection(
pid,
"LHMedium"))
1840 if(
eg->passSelection(
pid,
"LHTight"))
1845 ATH_MSG_INFO(
" REGTEST: problems with egamma pointer" );
1846 return StatusCode::SUCCESS;
1849 if (
eg->caloCluster()) {
1850 ATH_MSG_INFO(
" REGTEST: egamma cluster transverse energy: " <<
eg->caloCluster()->et() );
1851 ATH_MSG_INFO(
" REGTEST: egamma cluster eta: " <<
eg->caloCluster()->eta() );
1852 ATH_MSG_INFO(
" REGTEST: egamma cluster phi: " <<
eg->caloCluster()->phi() );
1853 double tmpeta = -999.;
1854 double tmpphi = -999.;
1857 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. etaCalo = " << tmpeta);
1858 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. phiCalo = " << tmpphi);
1860 ATH_MSG_INFO(
" REGTEST: problems with egamma cluster pointer" );
1862 ATH_MSG_INFO(
"REGTEST: Check the original (uncalibrated)");
1864 if (!orig.
isAvailable(*
eg->caloCluster()) || !orig(*
eg->caloCluster()).isValid()){
1869 ATH_MSG_INFO(
"REGTEST:: Compare new and old clusters");
1870 ATH_MSG_INFO(
"REGTEST:: Original Cluster e,eta,phi" << origClus->
e() <<
" " << origClus->
eta() <<
" " << origClus->
phi());
1871 ATH_MSG_INFO(
"REGTEST:: MVA Cluster e,eta,phi" <<
eg->caloCluster()->e() <<
" " <<
eg->caloCluster()->eta() <<
" " <<
eg->caloCluster()->phi());
1874 if(
eg->trackParticle()){
1876 ATH_MSG_INFO(
" REGTEST: charge= " <<
eg->trackParticle()->charge());
1877 ATH_MSG_INFO(
" REGTEST: E/p= " <<
eg->caloCluster()->et() /
eg->trackParticle()->pt() );
1879 ATH_MSG_INFO(
" REGTEST: Delta eta 1st sampling= " << val_float);
1881 ATH_MSG_INFO(
" REGTEST: Delta phi 2nd sampling= " << val_float);
1883 ATH_MSG_INFO(
" REGTEST: no electron eg->trackParticle() pointer");
1918 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::ElectronContainer DUMP===========" );
1920 return StatusCode::SUCCESS;
1929 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::PhotonContainer DUMP===========" );
1933 if (
sc.isFailure()) {
1934 ATH_MSG_INFO(
"REGTEST No Photon container HLT_xAOD__PhotonContainer_egamma_Photons");
1935 return StatusCode::SUCCESS;
1938 float val_float=-99;
1939 unsigned int isEMbit=0;
1941 ATH_MSG_INFO(
" REGTEST: xAOD Reconstruction variables: ");
1943 for (
const auto eg : *phCont){
1949 ATH_MSG_INFO(
" REGTEST: isEMLoose " <<
eg->selectionisEM(isEMbit,
"isEMLoose"));
1950 ATH_MSG_INFO(
" REGTEST: isEMLoose bit " << std::hex << isEMbit << std::dec);
1951 ATH_MSG_INFO(
" REGTEST: isEMMedium " <<
eg->selectionisEM(isEMbit,
"isEMMedium"));
1952 ATH_MSG_INFO(
" REGTEST: isEMMedium bit " << std::hex << isEMbit << std::dec);
1953 ATH_MSG_INFO(
" REGTEST: isEMTight " <<
eg->selectionisEM(isEMbit,
"isEMTight"));
1954 ATH_MSG_INFO(
" REGTEST: isEMTight bit " << std::hex << isEMbit << std::dec);
1956 ATH_MSG_INFO(
" REGTEST: problems with egamma pointer" );
1957 return StatusCode::SUCCESS;
1960 if (
eg->caloCluster()) {
1961 ATH_MSG_INFO(
" REGTEST: egamma cluster transverse energy: " <<
eg->caloCluster()->et() );
1962 ATH_MSG_INFO(
" REGTEST: egamma cluster eta: " <<
eg->caloCluster()->eta() );
1963 ATH_MSG_INFO(
" REGTEST: egamma cluster phi: " <<
eg->caloCluster()->phi() );
1964 double tmpeta = -999.;
1965 double tmpphi = -999.;
1968 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. etaCalo = " << tmpeta);
1969 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. phiCalo = " << tmpphi);
1971 ATH_MSG_INFO(
" REGTEST: problems with egamma cluster pointer" );
1973 ATH_MSG_INFO(
"REGTEST: Check the original (uncalibrated)");
1975 if (!orig.
isAvailable(*
eg->caloCluster()) || !orig(*
eg->caloCluster()).isValid()){
1980 ATH_MSG_INFO(
"REGTEST:: Compare new and old clusters");
1981 ATH_MSG_INFO(
"REGTEST:: Original Cluster e,eta,phi" << origClus->
e() <<
" " << origClus->
eta() <<
" " << origClus->
phi());
1982 ATH_MSG_INFO(
"REGTEST:: MVA Cluster e,eta,phi" <<
eg->caloCluster()->e() <<
" " <<
eg->caloCluster()->eta() <<
" " <<
eg->caloCluster()->phi());
2010 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::PhotonContainer DUMP===========" );
2012 return StatusCode::SUCCESS;
2020 ATH_MSG_INFO(
"REGTEST ==========START of TrigTauContainer DUMP===========");
2025 StatusCode sc = evtStore()->retrieve(trigTau,lastTrigTau);
2026 if (
sc.isFailure()) {
2028 return StatusCode::SUCCESS;
2032 for (
int i=0; trigTau != lastTrigTau; ++trigTau, ++
i) {
2039 for (
int j=0; TrigTauItr != TrigTauItrE; ++TrigTauItr, ++j ) {
2043 ATH_MSG_INFO(
"REGTEST TrigTau->pt() returns " << (*TrigTauItr)->pt());
2044 ATH_MSG_INFO(
"REGTEST TrigTau->eta() returns " << (*TrigTauItr)->eta());
2045 ATH_MSG_INFO(
"REGTEST TrigTau->phi() returns " << (*TrigTauItr)->phi());
2047 ATH_MSG_INFO(
"REGTEST TrigTau->nMatchedTracks() returns " << (*TrigTauItr)->nMatchedTracks());
2048 ATH_MSG_INFO(
"REGTEST TrigTau->simpleEtFlow() returns " << (*TrigTauItr)->simpleEtFlow());
2049 ATH_MSG_INFO(
"REGTEST TrigTau->m() returns " << (*TrigTauItr)->m());
2053 ATH_MSG_INFO(
"REGTEST ==========END of TrigTau DUMP===========");
2056 return StatusCode::SUCCESS;
2063 ATH_MSG_INFO(
"REGTEST ==========START of TrigTauTracksInfo DUMP===========");
2068 StatusCode sc = evtStore()->retrieve(trigTau,lastTrigTau);
2069 if (
sc.isFailure()) {
2071 return StatusCode::SUCCESS;
2075 for (
int i=0; trigTau != lastTrigTau; ++trigTau, ++
i) {
2081 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->pt() returns " << thisTrigTau->
pt());
2082 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->eta() returns " << thisTrigTau->
eta());
2083 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->phi() returns " << thisTrigTau->
phi());
2089 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->charge() returns " << thisTrigTau->
charge());
2090 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->leadingTrackPt() returns " <<
2092 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->scalarPtSumCore() returns " <<
2094 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->scalarPtSumIso() returns " <<
2096 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->threeFastestTracks().pt() returns " <<
2100 ATH_MSG_INFO(
"REGTEST ==========END of TrigTauTracksInfo DUMP===========");
2103 return StatusCode::SUCCESS;
2112 ATH_MSG_INFO(
"REGTEST ==========START of HLTResult DUMP===========");
2116 if(
sc.isFailure() || !hltResult ) {
2118 return StatusCode::SUCCESS;
2120 ATH_MSG_INFO(
"REGTEST HLTResult_L2 Successfully Retrieved");
2123 ATH_MSG_INFO(
"REGTEST HLTResult_L2->size() returns " << hltResult->
size());
2127 ATH_MSG_INFO(
"REGTEST ==========END of HLTResult DUMP===========");
2131 return StatusCode::SUCCESS;
2141 std::string TrigInDetTrackTags[]={
"HLT_TrigIDSCAN_Jet",
2142 "HLT_TrigIDSCAN_Tau",
2143 "HLT_TrigIDSCAN_eGamma",
2144 "HLT_TrigIDSCAN_Muon",
2145 "HLT_TrigIDSCAN_muonIso",
2146 "HLT_TrigIDSCAN_Bphysics",
2147 "HLT_TrigIDSCAN_FullScan",
2148 "HLT_TrigIDSCAN_Cosmics",
2149 "HLT_TrigIDSCAN_eGamma_Brem",
2150 "HLT_TrigIDSCAN_Tile",
2151 "HLT_TrigSiTrack_eGamma",
2152 "HLT_TrigSiTrack_Muon",
2153 "HLT_TrigSiTrack_muonIso",
2154 "HLT_TrigSiTrack_Tau",
2155 "HLT_TrigSiTrack_Jet",
2156 "HLT_TrigSiTrack_Bphysics",
2157 "HLT_TrigSiTrack_FullScan",
2158 "HLT_TrigSiTrack_Cosmics",
2159 "HLT_TRTSegmentFinder",
2166 bool gotTruthMap =
false;
2168 if (
evtStore()->contains<TrigInDetTrackTruthMap>(
"TrigInDetTrackTruthMap")) {
2172 ATH_MSG_WARNING(
" could not retrieve TrackTruthMap with key TrigInDetTruthMap");
2178 else ATH_MSG_DEBUG(
" didn't find any TrackTruthMap objects with key TrigInDetTruthMap");
2180 ATH_MSG_INFO(
"REGTEST ==========START of TrigInDetTrackCollection DUMP===========");
2182 for (
int iTag=0; iTag < ntag; iTag++) {
2185 if (
sc.isFailure()) {
2186 ATH_MSG_DEBUG(
"REGTEST No TrigInDetTrackCollection found with key " << TrigInDetTrackTags[iTag]);
2189 ATH_MSG_INFO(
"REGTEST TrigInDetTrackCollections retrieved with key "
2190 << TrigInDetTrackTags[iTag]);
2192 ATH_MSG_INFO(
"REGTEST TrigInDetTrackCollection->size() returns " << trigInDetTrackCollection->
size());
2197 for (
int j=0; trkItr != trkItrE; ++trkItr, ++j ) {
2199 ATH_MSG_INFO(
"REGTEST Looking at TrigInDetTrack " << j);
2201 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->algorithmId() returns " << (*trkItr)->algorithmId());
2202 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->chi2() returns " << (*trkItr)->chi2());
2203 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NStrawHits() returns " << (*trkItr)->NStrawHits());
2205 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NStraw() returns " << (*trkItr)->NStraw());
2206 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NStrawTime() returns " << (*trkItr)->NStrawTime());
2207 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NTRHits() returns " << (*trkItr)->NTRHits());
2209 ATH_MSG_INFO(
"REGTEST Looking at TrigInDetFitTrack->param()");
2212 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->a0() returns " << my_param->
a0());
2213 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->z0() returns " << my_param->
z0());
2214 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->phi0() returns " << my_param->
phi0());
2215 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->eta() returns " << my_param->
eta());
2216 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->pT() returns " << my_param->
pT());
2217 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->ea0() returns " << my_param->
ea0());
2218 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->ez0() returns " << my_param->
ez0());
2219 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->ephi0() returns " << my_param->
ephi0());
2220 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->eeta() returns " << my_param->
eeta());
2221 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->epT() returns " << my_param->
epT());
2224 ATH_MSG_INFO(
"REGTEST Looking at covarience matrix: ");
2226 const std::vector<double>*
const my_cov = my_param->
cov();
2228 ATH_MSG_INFO(
"REGTEST covarience matrix NULL (not stored)");
2230 std::vector<double> tempcov = *my_cov;
2233 for(iter = tempcov.begin(),
k=0; iter != tempcov.end(); ++iter, ++
k){
2234 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->cov() element " <<
k <<
" is " << (*iter));
2238 if (!(pTruthMap->
hasTruth((*trkItr)))) {
2254 if (pSiTruth->
cptr() != 0) {
2263 if (pTRTtruth->
cptr() != 0) {
2273 ATH_MSG_INFO(
"REGTEST ==========END of TrigInDetTrackCollection DUMP===========");
2276 return StatusCode::SUCCESS;
2285 ATH_MSG_INFO(
"REGTEST ==========START of TrigVertexCollection DUMP===========");
2290 StatusCode sc = evtStore()->retrieve(trigVertex,lastTrigVertex);
2291 if (
sc.isFailure()) {
2293 return StatusCode::SUCCESS;
2295 ATH_MSG_INFO(
"REGTEST TrigVertexCollection retrieved");
2298 for (
int i=0; trigVertex != lastTrigVertex; ++trigVertex, ++
i) {
2300 ATH_MSG_INFO(
"REGTEST Looking at TrigVertexCollection " <<
i);
2305 for (
int j=0; VertexItr != VertexItrE; ++VertexItr, ++j ) {
2309 ATH_MSG_INFO(
"REGTEST TrigVertex->x() returns " << (*VertexItr)->x());
2310 ATH_MSG_INFO(
"REGTEST TrigVertex->y() returns " << (*VertexItr)->y());
2311 ATH_MSG_INFO(
"REGTEST TrigVertex->z() returns " << (*VertexItr)->z());
2312 ATH_MSG_INFO(
"REGTEST TrigVertex->energyFraction() returns " << (*VertexItr)->energyFraction());
2313 ATH_MSG_INFO(
"REGTEST TrigVertex->ndof() returns " << (*VertexItr)->ndof());
2315 if(((*VertexItr)->tracks())!=NULL ){
2317 TrackInVertexList::const_iterator trkItr = (*VertexItr)->tracks()->begin();
2318 TrackInVertexList::const_iterator trkItrE = (*VertexItr)->tracks()->end();
2320 for (
int p=0; trkItr != trkItrE; ++trkItr, ++
p ) {
2325 ATH_MSG_INFO(
"REGTEST tracks->NStrawHits() returns " << (*trkItr)->NStrawHits());
2326 ATH_MSG_INFO(
"REGTEST tracks->NStrawTime() returns " << (*trkItr)->NStrawTime());
2327 ATH_MSG_INFO(
"REGTEST tracks->NTRHits() returns " << (*trkItr)->NTRHits());
2333 ATH_MSG_INFO(
"REGTEST ==========END of TrigVertexCollection DUMP===========");
2336 return StatusCode::SUCCESS;
2347 ATH_MSG_INFO(
"REGTEST ==========START of TrigEFBphysContainer DUMP===========");
2349 std::string EFBphysTags[]={
"HLT_xAOD__TrigBphysContainer_EFBMuMuFex",
2350 "HLT_xAOD__TrigBphysContainer_EFBMuMuXFex",
2351 "HLT_xAOD__TrigBphysContainer_EFDsPhiPiFex",
2352 "HLT_xAOD__TrigBphysContainer_EFMuPairs",
2353 "HLT_xAOD__TrigBphysContainer_EFMultiMuFex",
2354 "HLT_xAOD__TrigBphysContainer_EFTrackMass"
2357 int ntag= (
int)
sizeof(EFBphysTags) /
sizeof(EFBphysTags[0]);
2360 for (
int itag=0; itag<ntag; itag++){
2363 if (
sc.isFailure()) {
2364 ATH_MSG_INFO(
"REGTEST No TrigEFBphysContainer found with tag " << EFBphysTags[itag]);
2368 ATH_MSG_INFO(
"REGTEST TrigEFBphysContainer found with tag " << EFBphysTags[itag]
2369 <<
" and size " << trigEFBphys->
size());
2378 for (
int j=0; EFBphysItr != EFBphysItrE; ++EFBphysItr, ++j ) {
2382 ATH_MSG_INFO(
"REGTEST TrigEFBphys->eta() returns " << (*EFBphysItr)->eta());
2383 ATH_MSG_INFO(
"REGTEST TrigEFBphys->phi() returns " << (*EFBphysItr)->phi());
2384 ATH_MSG_INFO(
"REGTEST TrigEFBphys->mass() returns " << (*EFBphysItr)->mass());
2385 ATH_MSG_INFO(
"REGTEST TrigEFBphys->fitmass() returns " << (*EFBphysItr)->fitmass());
2387 ATH_MSG_INFO(
"REGTEST TrigEFBphys->roiId() returns " << (*EFBphysItr)->roiId());
2388 ATH_MSG_INFO(
"REGTEST TrigEFBphys->particleType() returns " << (*EFBphysItr)->particleType());
2390 if( (*EFBphysItr)->secondaryDecay() != NULL){
2393 ATH_MSG_INFO(
"REGTEST pSecondDecay->eta() returns " << psecond->
eta());
2394 ATH_MSG_INFO(
"REGTEST pSecondDecay->phi() returns " << psecond->
phi());
2405 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*EFBphysItr)->trackParticleLinks();
2406 if (trackVector.size() != 0) {
2407 ATH_MSG_INFO(
" REGTEST got track vector size: " << trackVector.size());
2411 std::vector<ElementLink<xAOD::TrackParticleContainer> >::const_iterator trkIt=trackVector.begin();
2412 for (
int itrk=0 ; trkIt!= trackVector.end(); ++itrk, ++trkIt) {
2413 if (!(trkIt->isValid())) {
2414 ATH_MSG_WARNING(
"TrackParticleContainer::Invalid ElementLink to track ");
2421 double phi = trackPerigee->parameters()[
Trk::phi];
2423 double px = trackPerigee->momentum()[
Trk::px];
2424 double py = trackPerigee->momentum()[
Trk::py];
2434 ATH_MSG_INFO(
"REGTEST ==========END of TrigEFBphysContainer DUMP===========");
2437 return StatusCode::SUCCESS;
2444 ATH_MSG_INFO(
"REGTEST ==========START of TrigL2BphysContainer DUMP===========");
2446 std::string L2BphysTags[]={
"HLT_xAOD__TrigBphysContainer_L2BMuMuFex",
2447 "HLT_xAOD__TrigBphysContainer_L2BMuMuXFex",
2448 "HLT_xAOD__TrigBphysContainer_L2DsPhiPiFexDs",
2449 "HLT_xAOD__TrigBphysContainer_L2DsPhiPiFexPhi",
2450 "HLT_xAOD__TrigBphysContainer_L2JpsieeFex",
2451 "HLT_xAOD__TrigBphysContainer_L2MultiMuFex",
2452 "HLT_xAOD__TrigBphysContainer_L2TrackMass",
2454 const int ntag = (
int)
sizeof(L2BphysTags) /
sizeof(L2BphysTags[0]);
2457 for (
int itag=0; itag<ntag; itag++){
2460 if (
sc.isFailure()) {
2461 ATH_MSG_INFO(
"REGTEST No TrigL2BphysContainer found with tag " << L2BphysTags[itag]);
2465 ATH_MSG_INFO(
"REGTEST TrigL2BphysContainer found with tag " << L2BphysTags[itag]
2466 <<
" and size " << trigL2Bphys->
size());
2475 for (
int j=0; L2BphysItr != L2BphysItrE; ++L2BphysItr, ++j ) {
2479 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->eta() returns " << (*L2BphysItr)->eta());
2480 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->phi() returns " << (*L2BphysItr)->phi());
2481 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->mass() returns " << (*L2BphysItr)->mass());
2482 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->fitmass() returns " << (*L2BphysItr)->fitmass());
2484 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->roiId() returns " << (*L2BphysItr)->roiId());
2485 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->particleType() returns " << (*L2BphysItr)->particleType());
2487 if( (*L2BphysItr)->secondaryDecay() != NULL){
2490 ATH_MSG_INFO(
"REGTEST pSecondDecay->eta() returns " << psecond->
eta());
2491 ATH_MSG_INFO(
"REGTEST pSecondDecay->phi() returns " << psecond->
phi());
2499 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*L2BphysItr)->trackParticleLinks();
2500 if (trackVector.size() != 0) {
2501 ATH_MSG_INFO(
" REGTEST got track vector size: " << trackVector.size());
2505 std::vector<ElementLink<xAOD::TrackParticleContainer> >::const_iterator trkIt=trackVector.begin();
2506 for (
int itrk=0 ; trkIt!= trackVector.end(); ++itrk, ++trkIt) {
2507 if (!(trkIt->isValid())) {
2508 ATH_MSG_WARNING(
"TrackParticleContainer::Invalid ElementLink to track ");
2515 double phi = trackPerigee->parameters()[
Trk::phi];
2517 double px = trackPerigee->momentum()[
Trk::px];
2518 double py = trackPerigee->momentum()[
Trk::py];
2528 ATH_MSG_INFO(
"REGTEST ==========END of TrigL2BphysContainer DUMP===========");
2531 return StatusCode::SUCCESS;
2542 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::JetContainer DUMP===========");
2544 std::string containerName[30] = {
"HLT_xAOD__JetContainer_a4tcemjesFS",
2545 "HLT_xAOD__JetContainer_a4tcemsubFS",
2546 "HLT_xAOD__JetContainer_a4tcemsubjesFS",
2547 "HLT_xAOD__JetContainer_a4tcemnojcalibFS",
2548 "HLT_xAOD__JetContainer_a4tcemjesPS",
2549 "HLT_xAOD__JetContainer_a4tcemnojcalibPS",
2550 "HLT_xAOD__JetContainer_a4tclcwjesFS",
2551 "HLT_xAOD__JetContainer_a4tclcwsubFS",
2552 "HLT_xAOD__JetContainer_a4tclcwsubjesFS",
2553 "HLT_xAOD__JetContainer_a4tclcwnojcalibFS",
2554 "HLT_xAOD__JetContainer_a4tclcwjesPS",
2555 "HLT_xAOD__JetContainer_a4tclcwnojcalibPS",
2556 "HLT_xAOD__JetContainer_a4TTemnojcalibFS",
2557 "HLT_xAOD__JetContainer_a4TThadnojcalibFS",
2558 "HLT_xAOD__JetContainer_a10tcemjesFS",
2559 "HLT_xAOD__JetContainer_a10tcemsubFS",
2560 "HLT_xAOD__JetContainer_a10tcemsubjesFS",
2561 "HLT_xAOD__JetContainer_a10tcemnojcalibFS",
2562 "HLT_xAOD__JetContainer_a10tcemjesPS",
2563 "HLT_xAOD__JetContainer_a10tcemnojcalibPS",
2564 "HLT_xAOD__JetContainer_a10tclcwjesFS",
2565 "HLT_xAOD__JetContainer_a10tclcwsubFS",
2566 "HLT_xAOD__JetContainer_a10tclcwsubjesFS",
2567 "HLT_xAOD__JetContainer_a10tclcwnojcalibFS",
2568 "HLT_xAOD__JetContainer_a10tclcwjesPS",
2569 "HLT_xAOD__JetContainer_a10tclcwnojcalibPS",
2570 "HLT_xAOD__JetContainer_a10TTemnojcalibFS",
2571 "HLT_xAOD__JetContainer_a10TThadnojcalibFS",
2572 "HLT_xAOD__JetContainer_a10r_tcemsubjesFS",
2573 "HLT_xAOD__JetContainer_TrigHLTJetDSSelectorCollection"};
2574 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,
2575 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};
2576 int containerInputType[30] = {1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 2,
2577 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 2, 1, 1};
2578 int containerSignalState[30] = {0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1,
2579 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0};
2580 bool onefilled =
false;
2582 for (
int icont=0; icont<30; ++icont) {
2587 if (
sc.isFailure()) {
2592 int jetContsize = jetCont->
size();
2593 ATH_MSG_INFO(
"REGTEST Got jet container " << containerName[icont] <<
", size: " << jetContsize);
2594 if (jetContsize != 0) {
2597 for(
const auto thisjet : *jetCont) {
2608 ATH_MSG_INFO(
"REGTEST rapidity: " << thisjet->rapidity() );
2613 ATH_MSG_INFO(
"REGTEST algorithm (kt: 0, cam: 1, antikt: 2, ...): " << thisjet->getAlgorithmType() <<
"; should be 2");
2614 if(thisjet->getAlgorithmType() != 2)
ATH_MSG_WARNING(
"Jet algorithm different from container");
2615 ATH_MSG_INFO(
"REGTEST size parameter: " << thisjet->getSizeParameter() <<
"; should be " << containerSizeParameter[icont]);
2616 if(thisjet->getSizeParameter() != containerSizeParameter[icont])
ATH_MSG_WARNING(
"Jet size different from container");
2617 ATH_MSG_INFO(
"REGTEST input (LCTopo: 0, EMTopo: 1, TopoTower: 2, ...): " << thisjet->getInputType() <<
"; should be " << containerInputType[icont]);
2618 if(thisjet->getInputType() != containerInputType[icont])
ATH_MSG_WARNING(
"Jet input different from container");
2619 ATH_MSG_INFO(
"REGTEST constituents signal state (uncalibrated: 0, calibrated: 1): " << thisjet->getConstituentsSignalState() <<
"; should be " << containerSignalState[icont]);
2620 if(thisjet->getConstituentsSignalState() != containerSignalState[icont])
ATH_MSG_WARNING(
"Jet constituents' signal state different from container");
2621 ATH_MSG_INFO(
"REGTEST number of constituents: " << thisjet->numConstituents() );
2628 unsigned int constitContsize = constitCont.
size();
2629 ATH_MSG_INFO(
"REGTEST Got constituent vector, size: " << constitContsize <<
"; should be " << thisjet->numConstituents());
2630 if(constitContsize != thisjet->numConstituents())
ATH_MSG_WARNING(
"Constituents container size different from number of constituents");
2632 if (constitContsize != 0) {
2730 std::vector<int> vecvalueint;
2732 int vecsize = vecvalueint.size();
2733 ATH_MSG_INFO(
"REGTEST Got NumTrkPt1000 vector, size: " << vecsize);
2736 for(
const auto & thisvalue : vecvalueint){
2738 ATH_MSG_INFO(
"REGTEST NumTrkPt1000 #" << j <<
": " << thisvalue);
2740 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2741 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2745 int vecsize = vecvalueint.size();
2746 ATH_MSG_INFO(
"REGTEST Got Got NumTrkPt500 vector, size: " << vecsize);
2749 for(
const auto & thisvalue : vecvalueint){
2751 ATH_MSG_INFO(
"REGTEST NumTrkPt500 #" << j <<
": " << thisvalue);
2753 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2754 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2760 std::vector<float> vecvalue;
2762 int vecsize = vecvalue.size();
2763 ATH_MSG_INFO(
"REGTEST Got JVF vector, size: " << vecsize);
2766 for(
const auto & thisvalue : vecvalue){
2768 ATH_MSG_INFO(
"REGTEST JVF #" << j <<
": " << thisvalue);
2770 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2771 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2775 int vecsize = vecvalue.size();
2776 ATH_MSG_INFO(
"REGTEST Got SumPtTrkPt1000 vector, size: " << vecsize);
2779 for(
const auto & thisvalue : vecvalue){
2781 ATH_MSG_INFO(
"REGTEST SumPtTrkPt1000 #" << j <<
": " << thisvalue);
2783 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2784 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2788 int vecsize = vecvalue.size();
2789 ATH_MSG_INFO(
"REGTEST Got SumPtTrkPt500 vector, size: " << vecsize);
2792 for(
const auto & thisvalue : vecvalue){
2794 ATH_MSG_INFO(
"REGTEST SumPtTrkPt500 #" << j <<
": " << thisvalue);
2796 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2797 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2801 int vecsize = vecvalue.size();
2802 ATH_MSG_INFO(
"REGTEST Got TrackWidthPt1000 vector, size: " << vecsize);
2805 for(
const auto & thisvalue : vecvalue){
2807 ATH_MSG_INFO(
"REGTEST TrackWidthPt1000 #" << j <<
": " << thisvalue);
2809 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2810 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2814 int vecsize = vecvalue.size();
2815 ATH_MSG_INFO(
"REGTEST Got TrackWidthPt500 vector, size: " << vecsize);
2818 for(
const auto & thisvalue : vecvalue){
2820 ATH_MSG_INFO(
"REGTEST TrackWidthPt500 #" << j <<
": " << thisvalue);
2822 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2823 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2827 int vecsize = vecvalue.size();
2828 ATH_MSG_INFO(
"REGTEST Got EnergyPerSampling vector, size: " << vecsize);
2831 for(
const auto & thisvalue : vecvalue){
2833 ATH_MSG_INFO(
"REGTEST EnergyPerSampling #" << j <<
": " << thisvalue);
2835 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2836 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2865 return StatusCode::SUCCESS;
2869 if (jetContsize ==
i)
ATH_MSG_INFO(
"REGTEST size of jet container == number of displayed jets: " << (jetContsize ==
i) );
2874 if (!onefilled)
ATH_MSG_DEBUG(
"There was no filled jet containers");
2876 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::JetContainer DUMP===========");
2880 return StatusCode::SUCCESS;
2889 ATH_MSG_INFO(
"REGTEST ==========START of TrigEFBjetContainer DUMP===========");
2894 StatusCode sc = evtStore()->retrieve(trigEFBjet,lastTrigEFBjet);
2895 if (
sc.isFailure()) {
2897 return StatusCode::SUCCESS;
2899 ATH_MSG_INFO(
"REGTEST TrigEFBjetContainers retrieved");
2902 for (
int i=0; trigEFBjet != lastTrigEFBjet; ++trigEFBjet, ++
i) {
2904 ATH_MSG_INFO(
"REGTEST Looking at TrigEFBjetContainer " <<
i);
2909 for (
int j=0; EFBjetItr != EFBjetItrE; ++EFBjetItr, ++j ) {
2913 ATH_MSG_INFO(
"REGTEST TrigEFBjet->prmVtx() returns " << (*EFBjetItr)->prmVtx());
2914 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xComb() returns " << (*EFBjetItr)->xComb());
2916 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xIP1D() returns " << (*EFBjetItr)->xIP1D());
2917 ATH_MSG_INFO(
"REGTEST TrigEFBjet->isValid() returns " << (*EFBjetItr)->isValid());
2918 ATH_MSG_INFO(
"REGTEST TrigEFBjet->roiId() returns " << (*EFBjetItr)->roiId());
2920 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xIP2D() returns " << (*EFBjetItr)->xIP2D());
2922 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xIP3D() returns " << (*EFBjetItr)->xIP3D());
2924 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xCHI2() returns " << (*EFBjetItr)->xCHI2());
2926 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xSV() returns " << (*EFBjetItr)->xSV());
2928 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xMVtx() returns " << (*EFBjetItr)->xMVtx());
2929 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xEVtx() returns " << (*EFBjetItr)->xEVtx());
2930 ATH_MSG_INFO(
"REGTEST TrigEFBjet-> xNVtx() returns " << (*EFBjetItr)-> xNVtx());
2933 ATH_MSG_INFO(
"REGTEST ==========END of TrigEFBjetContainer DUMP===========");
2936 return StatusCode::SUCCESS;
2945 ATH_MSG_INFO(
"REGTEST ==========START of TrigL2BjetContainer DUMP===========");
2950 StatusCode sc = evtStore()->retrieve(trigL2Bjet,lastTrigL2Bjet);
2951 if (
sc.isFailure()) {
2953 return StatusCode::SUCCESS;
2955 ATH_MSG_INFO(
"REGTEST TrigL2BjetContainers retrieved");
2958 for (
int i=0; trigL2Bjet != lastTrigL2Bjet; ++trigL2Bjet, ++
i) {
2960 ATH_MSG_INFO(
"REGTEST Looking at TrigL2BjetContainer " <<
i);
2965 for (
int j=0; L2BjetItr != L2BjetItrE; ++L2BjetItr, ++j ) {
2968 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->prmVtx() returns " << (*L2BjetItr)->prmVtx());
2969 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xComb() returns " << (*L2BjetItr)->xComb());
2971 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xIP1D() returns " << (*L2BjetItr)->xIP1D());
2972 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->isValid() returns " << (*L2BjetItr)->isValid());
2973 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->roiId() returns " << (*L2BjetItr)->roiId());
2975 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xIP2D() returns " << (*L2BjetItr)->xIP2D());
2977 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xIP3D() returns " << (*L2BjetItr)->xIP3D());
2979 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xCHI2() returns " << (*L2BjetItr)->xCHI2());
2981 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xSV() returns " << (*L2BjetItr)->xSV());
2982 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xMVtx() returns " << (*L2BjetItr)->xMVtx());
2983 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xEVtx() returns " << (*L2BjetItr)->xEVtx());
2984 ATH_MSG_INFO(
"REGTEST TrigL2Bjet-> xNVtx() returns " << (*L2BjetItr)-> xNVtx());
2987 ATH_MSG_INFO(
"REGTEST ==========END of TrigL2BjetContainer DUMP===========");
2990 return StatusCode::SUCCESS;
2999 ATH_MSG_INFO(
"REGTEST ==========START of MuonFeature DUMP===========");
3004 StatusCode sc = evtStore()->retrieve(MuFeature,lastMuFeature);
3005 if (
sc.isFailure()) {
3007 return StatusCode::SUCCESS;
3012 for (
int i=0; MuFeature != lastMuFeature; ++MuFeature, ++
i) {
3017 ATH_MSG_INFO(
"REGTEST MuonFeature->roiId() returns " << thisMuFeature->
roiId());
3018 ATH_MSG_INFO(
"REGTEST MuonFeature->eta() returns " << thisMuFeature->
eta());
3019 ATH_MSG_INFO(
"REGTEST MuonFeature->phi() returns " << thisMuFeature->
phi());
3021 ATH_MSG_INFO(
"REGTEST MuonFeature->pt() returns " << thisMuFeature->
pt());
3024 ATH_MSG_INFO(
"REGTEST MuonFeature->zeta() returns " << thisMuFeature->
zeta());
3026 ATH_MSG_INFO(
"REGTEST MuonFeature->beta() returns " << thisMuFeature->
beta());
3029 ATH_MSG_INFO(
"REGTEST ==========END of MuonFeature DUMP===========");
3032 return StatusCode::SUCCESS;
3041 ATH_MSG_INFO(
"REGTEST ==========START of CombinedMuonFeature DUMP===========");
3046 StatusCode sc = evtStore()->retrieve(CombMuon,lastCombMuon);
3047 if (
sc.isFailure()) {
3049 return StatusCode::SUCCESS;
3051 ATH_MSG_INFO(
"REGTEST CombinedMuonFeatures retrieved");
3053 for (
int i=0; CombMuon != lastCombMuon; ++CombMuon, ++
i) {
3057 ATH_MSG_INFO(
"REGTEST Looking at CombinedMuonFeature " <<
i);
3058 ATH_MSG_INFO(
"REGTEST TrigPhoton->pt() returns " << (thisCombMuFeature)->
pt());
3059 ATH_MSG_INFO(
"REGTEST TrigPhoton->sigma_pt() returns " << (thisCombMuFeature)->sigma_pt());
3062 if ((thisCombMuFeature)->muFastTrackLink().
isValid() ) {
3063 ATH_MSG_INFO(
"REGTEST muFastTrack->eta() returns " << (thisCombMuFeature)->muFastTrack()->
eta());
3064 ATH_MSG_INFO(
"REGTEST muFastTrack->phi() returns " << (thisCombMuFeature)->muFastTrack()->
phi());
3065 ATH_MSG_INFO(
"REGTEST muFastTrack->radius() returns " << (thisCombMuFeature)->muFastTrack()->
radius());
3066 ATH_MSG_INFO(
"REGTEST muFastTrack->zeta() returns " << (thisCombMuFeature)->muFastTrack()->zeta());
3068 ATH_MSG_INFO(
"CombinedMuonFeature has no muFastTrack" );
3073 if ((thisCombMuFeature)->IDTrackLink().
isValid() ) {
3074 ATH_MSG_INFO(
"REGTEST IDTrack->algorithmId() returns " <<(thisCombMuFeature)->IDTrack()->algorithmId());
3075 ATH_MSG_INFO(
"REGTEST IDTrack->chi2() returns " << (thisCombMuFeature)->IDTrack()->
chi2());
3076 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawHits() returns " <<(thisCombMuFeature)->IDTrack()->NStrawHits());
3077 ATH_MSG_INFO(
"REGTEST IDTrack->NStraw() returns " << (thisCombMuFeature)->IDTrack()->NStraw());
3078 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawTime() returns " <<(thisCombMuFeature)->IDTrack()->NStrawTime());
3079 ATH_MSG_INFO(
"REGTEST IDTrack->NTRHits() returns " <<(thisCombMuFeature)->IDTrack()->NTRHits());
3086 ATH_MSG_INFO(
"REGTEST ==========END of CombinedMuonFeature DUMP===========");
3089 return StatusCode::SUCCESS;
3095 ATH_MSG_INFO(
"REGTEST ==========START of CombinedMuonFeatureContainer DUMP===========");
3100 StatusCode sc = evtStore()->retrieve(CombMuon,lastCombMuon);
3101 if (
sc.isFailure()) {
3102 ATH_MSG_INFO(
"REGTEST No CombinedMuonFeatureContainer found");
3103 return StatusCode::SUCCESS;
3105 ATH_MSG_INFO(
"REGTEST CombinedMuonFeaturesContainer retrieved");
3107 for (
int j=0; CombMuon != lastCombMuon; ++CombMuon, ++j) {
3108 ATH_MSG_INFO(
"REGTEST Looking at CombinedMuonFeatureContainer " << j);
3116 ATH_MSG_INFO(
"REGTEST Looking at CombinedMuonFeature " <<
i);
3118 ATH_MSG_INFO(
"REGTEST CombinedMuonFeature->pt() returns " << (thisCombMuFeature)->
pt());
3119 ATH_MSG_INFO(
"REGTEST CombinedMuonFeature->sigma_pt() returns " << (thisCombMuFeature)->sigma_pt());
3123 ATH_MSG_INFO(
"REGTEST muFastTrack->eta() returns " << (thisCombMuFeature)->muFastTrack()->
eta());
3124 ATH_MSG_INFO(
"REGTEST muFastTrack->phi() returns " << (thisCombMuFeature)->muFastTrack()->
phi());
3125 ATH_MSG_INFO(
"REGTEST muFastTrack->radius() returns " << (thisCombMuFeature)->muFastTrack()->
radius());
3126 ATH_MSG_INFO(
"REGTEST muFastTrack->zeta() returns " << (thisCombMuFeature)->muFastTrack()->zeta());
3133 ATH_MSG_INFO(
"REGTEST IDTrack->algorithmId() returns " <<(thisCombMuFeature)->IDTrack()->algorithmId());
3134 ATH_MSG_INFO(
"REGTEST IDTrack->chi2() returns " << (thisCombMuFeature)->IDTrack()->
chi2());
3135 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawHits() returns " <<(thisCombMuFeature)->IDTrack()->NStrawHits());
3137 ATH_MSG_INFO(
"REGTEST IDTrack->NStraw() returns " << (thisCombMuFeature)->IDTrack()->NStraw());
3138 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawTime() returns " <<(thisCombMuFeature)->IDTrack()->NStrawTime());
3139 ATH_MSG_INFO(
"REGTEST IDTrack->NTRHits() returns " <<(thisCombMuFeature)->IDTrack()->NTRHits());
3146 return StatusCode::SUCCESS;
3154 ATH_MSG_INFO(
"REGTEST ==========START of TrigEMCluster DUMP===========");
3159 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3160 if (
sc.isFailure()) {
3162 return StatusCode::SUCCESS;
3166 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3173 ATH_MSG_INFO(
"REGTEST TrigEMCluster->e() returns " << thisEMCluster->
e());
3174 ATH_MSG_INFO(
"REGTEST TrigEMCluster->phi() returns " << thisEMCluster->
phi());
3176 ATH_MSG_INFO(
"REGTEST TrigEMCluster->eta() returns " << thisEMCluster->
eta());
3177 ATH_MSG_INFO(
"REGTEST TrigEMCluster->e237() returns " << thisEMCluster->
e237());
3178 ATH_MSG_INFO(
"REGTEST TrigEMCluster->e277() returns " << thisEMCluster->
e277());
3180 ATH_MSG_INFO(
"REGTEST TrigEMCluster->weta2() returns " << thisEMCluster->
weta2());
3181 ATH_MSG_INFO(
"REGTEST TrigEMCluster->ehad1() returns " << thisEMCluster->
ehad1());
3182 ATH_MSG_INFO(
"REGTEST TrigEMCluster->Eta1() returns " << thisEMCluster->
Eta1());
3185 ATH_MSG_INFO(
"REGTEST ==========END of TrigEMCluster DUMP===========");
3188 return StatusCode::SUCCESS;
3195 ATH_MSG_INFO(
"REGTEST ==========START of TrigEMCluster DUMP===========");
3200 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3201 if (
sc.isFailure()) {
3203 return StatusCode::SUCCESS;
3208 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3212 ATH_MSG_INFO(
"REGTEST Looking at xAOD::TrigEMCluster " <<
i);
3214 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->energy() returns " << thisEMCluster->
energy());
3215 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->e() returns " << thisEMCluster->
energy());
3216 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->phi() returns " << thisEMCluster->
phi());
3218 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->eta() returns " << thisEMCluster->
eta());
3219 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->e237() returns " << thisEMCluster->
e237());
3220 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->e277() returns " << thisEMCluster->
e277());
3221 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->fracs1() returns " << thisEMCluster->
fracs1());
3222 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->weta2() returns " << thisEMCluster->
weta2());
3223 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->ehad1() returns " << thisEMCluster->
ehad1());
3224 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->Eta1() returns " << thisEMCluster->
eta1());
3227 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigEMCluster DUMP===========");
3230 return StatusCode::SUCCESS;
3239 ATH_MSG_INFO(
"REGTEST ==========START of TrigTauClusterContainer/TrigTauClusterDetailsContainer DUMP===========");
3244 StatusCode sc = evtStore()->retrieve(TauCluster,lastTauCluster);
3245 if (
sc.isFailure()) {
3246 ATH_MSG_INFO(
"REGTEST No TrigTauClusterContainer found");
3247 return StatusCode::SUCCESS;
3249 ATH_MSG_INFO(
"REGTEST TrigTauClusterContainer retrieved");
3253 for (
int i=0; TauCluster != lastTauCluster; ++TauCluster, ++
i) {
3255 ATH_MSG_INFO(
"REGTEST Looking at TrigTauClusterContainer " <<
i);
3260 for (
int j=0; TauClusterItr != TauClusterItrE; ++TauClusterItr, ++j ) {
3262 ATH_MSG_INFO(
"REGTEST Looking at TrigTauCluster " << j);
3264 ATH_MSG_INFO(
"REGTEST TrigTauCluster->energy() returns " << (*TauClusterItr)->energy());
3265 ATH_MSG_INFO(
"REGTEST TrigTauCluster->et() returns " << (*TauClusterItr)->et());
3266 ATH_MSG_INFO(
"REGTEST TrigTauCluster->EMCalibEnergy() returns " << (*TauClusterItr)->EMCalibEnergy());
3268 ATH_MSG_INFO(
"REGTEST TrigTauCluster->EMenergy() returns " << (*TauClusterItr)->EMenergy());
3269 ATH_MSG_INFO(
"REGTEST TrigTauCluster->HADenergy() returns " << (*TauClusterItr)->HADenergy());
3270 ATH_MSG_INFO(
"REGTEST TrigTauCluster->eta() returns " << (*TauClusterItr)->eta());
3271 ATH_MSG_INFO(
"REGTEST TrigTauCluster->phi() returns " << (*TauClusterItr)->phi());
3272 ATH_MSG_INFO(
"REGTEST TrigTauCluster->IsoFrac() returns " << (*TauClusterItr)->IsoFrac());
3273 ATH_MSG_INFO(
"REGTEST TrigTauCluster->stripWidth() returns " << (*TauClusterItr)->stripWidth());
3274 ATH_MSG_INFO(
"REGTEST TrigTauCluster->numStripCells() returns " << (*TauClusterItr)->numStripCells());
3275 ATH_MSG_INFO(
"REGTEST TrigTauCluster->CaloRadius() returns " << (*TauClusterItr)->CaloRadius());
3276 ATH_MSG_INFO(
"REGTEST TrigTauCluster->numTotCells() returns " << (*TauClusterItr)->numTotCells());
3279 if( (*TauClusterItr)->clusterDetails() == 0 )
3282 ATH_MSG_INFO(
"REGTEST TrigTauCluster->Energy in Narrow cone : EM[0/1/2/3]="
3283 <<(*TauClusterItr)->EMenergyNarrow(0) <<
" " << (*TauClusterItr)->EMenergyNarrow(1) <<
" "
3284 <<(*TauClusterItr)->EMenergyNarrow(2) <<
" ");
3285 ATH_MSG_INFO(
"REGTEST TrigTauCluster->Energy in Narrow cone : HAD[0/1/2]="
3286 <<(*TauClusterItr)->HADenergyNarrow(0) <<
" " << (*TauClusterItr)->HADenergyNarrow(1) <<
" "
3287 <<(*TauClusterItr)->HADenergyNarrow(2));
3295 sc = evtStore()->retrieve(TauDetailsCluster,lastTauDetailsCluster);
3296 if (
sc.isFailure()) {
3297 ATH_MSG_INFO(
"REGTEST No TrigTauDetailsClusterContainer found");
3298 return StatusCode::SUCCESS;
3300 ATH_MSG_INFO(
"REGTEST TrigTauDetailsClusterContainer retrieved");
3302 for (
int i=0; TauDetailsCluster != lastTauDetailsCluster; ++TauDetailsCluster, ++
i) {
3304 ATH_MSG_INFO(
"REGTEST Looking at TrigTauClusterDetailsContainer " <<
i);
3309 for (
int j=0; TauDetailsClusterItr != TauDetailsClusterItrE; ++TauDetailsClusterItr, ++j ) {
3311 ATH_MSG_INFO(
"REGTEST Looking at TrigTauClusterDetails " << j);
3313 for(
unsigned int i = 0;
i<4; ++
i ) {
3314 ATH_MSG_INFO(
"REGTEST TrigTauClusterDetails-> EM["<<
i <<
"] Radius/EnergyNar/EnergyMid/EnergyWid returns "
3315 << (*TauDetailsClusterItr)->EMRadius(
i) <<
" "
3316 << (*TauDetailsClusterItr)->EMenergyNarrow(
i) <<
" "
3317 << (*TauDetailsClusterItr)->EMenergyMedium(
i) <<
" "
3318 << (*TauDetailsClusterItr)->EMenergyWide(
i));
3320 for(
unsigned int i = 0;
i<3; ++
i ) {
3321 ATH_MSG_INFO(
"REGTEST TrigTauClusterDetails-> HAD["<<
i <<
"] Radius/EnergyNar/EnergyMid/EnergyWid returns "
3322 << (*TauDetailsClusterItr)->HADRadius(
i) <<
" "
3323 << (*TauDetailsClusterItr)->HADenergyNarrow(
i) <<
" "
3324 << (*TauDetailsClusterItr)->HADenergyMedium(
i) <<
" "
3325 << (*TauDetailsClusterItr)->HADenergyWide(
i));
3330 if( nDetails != nClusters)
3331 ATH_MSG_WARNING(
"REGTEST inconsistent number of TrigTauClusters ("<< nClusters<<
") and TrigTauClusterDetails ("
3332 << nDetails <<
")");
3334 ATH_MSG_INFO(
"REGTEST ==========END of TrigTauClusterContainer/TrigTauClusterDetailsContainer DUMP===========");
3338 return StatusCode::SUCCESS;
3347 ATH_MSG_INFO(
"REGTEST ==========START of TrigEMClusterContainer DUMP===========");
3352 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3353 if (
sc.isFailure()) {
3354 ATH_MSG_INFO(
"REGTEST No TrigEMClusterContainer found");
3355 return StatusCode::SUCCESS;
3357 ATH_MSG_INFO(
"REGTEST TrigEMClusterContainer retrieved");
3360 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3362 ATH_MSG_INFO(
"REGTEST Looking at TrigEMClusterContainer " <<
i);
3367 for (
int j=0; EMClusterItr != EMClusterItrE; ++EMClusterItr, ++j ) {
3369 ATH_MSG_INFO(
"REGTEST Looking at TrigEMCluster " << j);
3371 ATH_MSG_INFO(
"REGTEST TrigEMCluster->energy() returns " << (*EMClusterItr)->energy());
3372 ATH_MSG_INFO(
"REGTEST TrigEMCluster->et() returns " << (*EMClusterItr)->et());
3373 ATH_MSG_INFO(
"REGTEST TrigEMCluster->eta() returns " << (*EMClusterItr)->eta());
3374 ATH_MSG_INFO(
"REGTEST TrigEMCluster->phi() returns " << (*EMClusterItr)->phi());
3379 (*EMClusterItr)->print(
msg());
3384 ATH_MSG_INFO(
"REGTEST ==========END of TrigEMClusterContainer DUMP===========");
3387 return StatusCode::SUCCESS;
3394 ATH_MSG_INFO(
"REGTEST ==========START of xAODTrigEMClusterContainer DUMP===========");
3399 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3400 if (
sc.isFailure()) {
3401 ATH_MSG_INFO(
"REGTEST No xAOD::TrigEMClusterContainer found");
3402 return StatusCode::SUCCESS;
3404 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMClusterContainer retrieved");
3407 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3409 ATH_MSG_INFO(
"REGTEST Looking at xAOD::TrigEMClusterContainer " <<
i);
3414 for (
int j=0; EMClusterItr != EMClusterItrE; ++EMClusterItr, ++j ) {
3416 ATH_MSG_INFO(
"REGTEST Looking at xAOD::TrigEMCluster " << j);
3417 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->energy() returns " << (*EMClusterItr)->energy());
3418 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->et() returns " << (*EMClusterItr)->et());
3419 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->eta() returns " << (*EMClusterItr)->eta());
3420 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->phi() returns " << (*EMClusterItr)->phi());
3430 ATH_MSG_INFO(
"REGTEST ==========END of TrigEMClusterContainer DUMP===========");
3433 return StatusCode::SUCCESS;
3441 ATH_MSG_INFO(
"REGTEST ==========START of TileMuFeatureContainer DUMP===========");
3447 if (
sc.isFailure()) {
3448 ATH_MSG_INFO(
"REGTEST No TileMuFeatureContainer found");
3449 return StatusCode::SUCCESS;
3451 ATH_MSG_INFO(
"REGTEST TileMuFeatureContainer retrieved");
3454 ATH_MSG_INFO(
"REGTEST Looking at TileMuFeatureContainer " <<
i);
3458 for(
int j=0; TileMuItr != TileMuItrE; ++TileMuItr, ++j) {
3459 ATH_MSG_INFO(
"REGTEST Looking at TileMuFeature " << j);
3460 ATH_MSG_INFO(
"REGTEST TileMuFeature->eta() returns " << (*TileMuItr)->eta());
3461 ATH_MSG_INFO(
"REGTEST TileMuFeature->phi() returns " << (*TileMuItr)->phi());
3462 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(0) returns " << (*TileMuItr)->enedep().at(0));
3463 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(1) returns " << (*TileMuItr)->enedep().at(1));
3464 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(2) returns " << (*TileMuItr)->enedep().at(2));
3465 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(3) returns " << (*TileMuItr)->enedep().at(3));
3466 ATH_MSG_INFO(
"REGTEST TileMuFeature->quality() returns " << (*TileMuItr)->quality());
3471 ATH_MSG_INFO(
"REGTEST ==========END of TileMuFeatureContainer DUMP===========");
3474 return StatusCode::SUCCESS;
3482 ATH_MSG_INFO(
"REGTEST ==========START of TileTrackMuFeatureContainer DUMP===========");
3487 StatusCode sc = evtStore()->retrieve(TileTrackMu, lastTileTrackMu);
3488 if (
sc.isFailure()) {
3489 ATH_MSG_INFO(
"REGTEST No TileTrackMuFeatureContainer found");
3490 return StatusCode::SUCCESS;
3492 ATH_MSG_INFO(
"REGTEST TileTrackMuFeatureContainer retrieved");
3494 for (
int i=0; TileTrackMu!=lastTileTrackMu; ++TileTrackMu, ++
i) {
3495 ATH_MSG_INFO(
"REGTEST Looking at TileTrackMuFeatureContainer " <<
i);
3500 for (
int j=0; TileTrackMuItr != TileTrackMuItrE; ++TileTrackMuItr, ++j) {
3501 ATH_MSG_INFO(
"REGTEST Looking at TileTrackMuFeature " << j);
3502 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->PtTR_Trk() returns " << (*TileTrackMuItr)->PtTR_Trk());
3503 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->EtaTR_Trk() returns " << (*TileTrackMuItr)->EtaTR_Trk());
3504 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->PhiTR_Trk() returns " << (*TileTrackMuItr)->PhiTR_Trk());
3505 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->Typ_IDTrk() returns " << (*TileTrackMuItr)->Typ_IDTrk());
3511 ATH_MSG_INFO(
"REGTEST No TileMuFeature (Something Wrong)");
3538 ATH_MSG_INFO(
"REGTEST IDScanLink->algorithmId() returns " <<
Track->algorithmId());
3540 ATH_MSG_INFO(
"REGTEST IDScanLink->NStrawHits() returns ");
3542 ATH_MSG_INFO(
"REGTEST IDScanLink->NStrawTime() returns " <<
Track->NStrawTime());
3544 ATH_MSG_INFO(
"REGTEST IDScanLink->param()->phi0() returns " <<
Track->param()->phi0());
3545 ATH_MSG_INFO(
"REGTEST IDScanLink->param()->eta() returns " <<
Track->param()->eta());
3546 ATH_MSG_INFO(
"REGTEST IDScanLink->param()->pT() returns " <<
Track->param()->pT());
3551 ATH_MSG_INFO(
"REGTEST ==========END of TileTrackMuFeatureContainer DUMP===========");
3552 ATH_MSG_DEBUG(
"dumpTileTrackMuFeatureContainer() succeeded");
3553 return StatusCode::SUCCESS;
3560 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TauJetContainer DUMP===========" );
3562 StatusCode sc = evtStore() ->
retrieve (TauJetcont,
"HLT_xAOD__TauJetContainer_TrigTauRecMerged");
3564 if (
sc.isFailure()) {
3566 ATH_MSG_INFO(
"REGTEST No Tau container HLT_xAOD__TauJetContainer_TrigTauRecMerged");
3569 return StatusCode::SUCCESS;
3575 ATH_MSG_INFO(
"REGTEST (*tauIt)->eta() returns " << (*tauIt)->eta() );
3576 ATH_MSG_INFO(
"REGTEST (*tauIt)->phi() returns " << (*tauIt)->phi() );
3577 ATH_MSG_INFO(
"REGTEST (*tauIt)->pt() returns " << (*tauIt)->pt() );
3581 #ifndef XAODTAU_VERSIONS_TAUJET_V3_H
3582 EFnTracks = (*tauIt)->nTracks();
3587 ATH_MSG_INFO(
"REGTEST (*tauIt)->nTracks() returns " << EFnTracks );
3590 int EFWidenTrack = -1;
3591 #ifndef XAODTAU_VERSIONS_TAUJET_V3_H
3592 EFWidenTrack = (*tauIt)->nWideTracks();
3597 ATH_MSG_INFO(
"REGTEST (*tauIt)->nWideTracks() returns " << EFWidenTrack );
3620 ATH_MSG_INFO(
"REGTEST TauDetails->etOverPtLeadTrk() returns " << etOvPtLead);
3629 ATH_MSG_INFO(
"REGTEST TauDetails->isolFrac() returns " << IsoFrac);
3658 ATH_MSG_INFO(
"REGTEST TauDetails->etEMAtEMScale() returns " << EtEm);}
3662 ATH_MSG_INFO(
"REGTEST TauDetails->etHadAtEMScale() returns " << EtHad);}
3665 if( !(*tauIt)->jetLink().isValid() ) {
3667 return StatusCode::SUCCESS;
3670 const xAOD::Jet* pJetSeed = *((*tauIt)->jetLink());
3675 for (
int clusCount = 0; clusItr != clusItrE; ++clusItr, ++clusCount) {
3677 ATH_MSG_INFO(
"REGTEST Tau Cluster " << clusCount <<
" pt = " << (*clusItr)->pt()
3678 <<
" eta = " << (*clusItr)->eta()
3679 <<
" phi = " << (*clusItr)->phi() );
3684 for (
unsigned int trackNum = 0; trackNum < (*tauIt)->nTracks(); ++trackNum) {
3689 return StatusCode::SUCCESS;
3694 ATH_MSG_INFO(
"REGTEST Tau linked track " << trackNum <<
" pt = " << linkTrack->
pt()
3695 <<
" eta = " << linkTrack->
eta()
3696 <<
" phi = " << linkTrack->
phi() );
3703 return StatusCode::SUCCESS;
3712 ATH_MSG_INFO(
"REGTEST ==========START of TauJetContainer DUMP===========");
3718 std::string TauContainerTags[]={
"HLT_TrigTauRecMerged"};
3719 for (
int itag=0; itag < ntag; itag++) {
3721 sCode=evtStore()->retrieve(TauJetcont , TauContainerTags[itag]);
3722 if( sCode.isFailure() ){
3723 ATH_MSG_INFO(
"Failed to retrieve TauJetContainer with key " << TauContainerTags[itag]);
3729 bool IsTaurec =
false;
3730 bool Is1p3p =
false;
3751 ATH_MSG_INFO(
"REGTEST (*tauIt)->eta() returns " << (*tauIt)->eta());
3752 ATH_MSG_INFO(
"REGTEST (*tauIt)->phi() returns " << (*tauIt)->phi());
3753 ATH_MSG_INFO(
"REGTEST (*tauIt)->numTrack() returns " << (*tauIt)->numTrack());
3757 ATH_MSG_INFO(
"REGTEST TauDetails->ipZ0SinThetaSigLeadTrk() returns " <<
TauDetails->ipZ0SinThetaSigLeadTrk());
3770 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_stripWidth2() returns " <<
TauDetails->seedCalo_stripWidth2());
3772 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_etEMAtEMScale() returns " <<
TauDetails->seedCalo_etEMAtEMScale());
3773 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_etHADAtEMScale() returns " <<
TauDetails->seedCalo_etHadAtEMScale());
3775 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_nIsolLooseTrk() returns " <<
TauDetails->seedCalo_nIsolLooseTrk());
3784 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_etChrgHadOverSumTrkPt () returns " <<
TauDetails->seedTrk_etChrgHadOverSumTrkPt ());
3785 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_isolFracWide() returns " <<
TauDetails->seedTrk_isolFracWide());
3786 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_etHadAtEMScale() returns " <<
TauDetails->seedTrk_etHadAtEMScale());
3787 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_etEMAtEMScale() returns " <<
TauDetails->seedTrk_etEMAtEMScale());
3793 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_sumEMCellEtOverLeadTrkPt() returns " <<
TauDetails->seedTrk_sumEMCellEtOverLeadTrkPt());
3794 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_secMaxStripEt() returns " <<
TauDetails->seedTrk_secMaxStripEt());
3798 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_nOtherCoreTrk() returns " <<
TauDetails->seedTrk_nOtherCoreTrk());
3812 ATH_MSG_INFO(
"REGTEST TauDetails->looseTrk() link present ");
3814 ATH_MSG_INFO(
"REGTEST TauDetails->looseTrk() link missing ");
3851 if (
sc.isFailure()) {
3853 return StatusCode::FAILURE;
3863 return StatusCode::SUCCESS;
3870 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrackParticle DUMP===========");
3872 std::vector<std::string> SGkeys;
3873 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Bjet_EFID");
3874 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Bphysics_EFID");
3875 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Electron_EFID");
3876 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_FullScan_EFID");
3877 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Muon_EFID");
3878 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Photon_EFID");
3879 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Tau_EFID");
3886 if (
sc.isFailure()) {
3887 ATH_MSG_INFO(
"REGTEST No track particle container found with key " << SGkeys.at(
SGkey));
3890 ATH_MSG_INFO(
"REGTEST TrackParticleContainer retrieved with key " << SGkeys.at(
SGkey)
3891 <<
" and size " << trackParticleContainer->
size());
3896 for (
int index = 0; trackParticleItr != trackParticleLast; ++trackParticleItr, ++
index) {
3900 ATH_MSG_INFO(
"REGTEST pt: " << (*trackParticleItr)->pt()
3901 <<
"/eta: " << (*trackParticleItr)->eta()
3902 <<
"/phi: " << (*trackParticleItr)->phi()
3903 <<
"/m: " << (*trackParticleItr)->m()
3904 <<
"/e: " << (*trackParticleItr)->e()
3905 <<
"/rapidity: " << (*trackParticleItr)->rapidity());
3907 ATH_MSG_INFO(
"REGTEST Defining parameters functions:");
3908 ATH_MSG_INFO(
"REGTEST charge: " << (*trackParticleItr)->charge()
3909 <<
"/d0: " << (*trackParticleItr)->d0()
3910 <<
"/z0: " << (*trackParticleItr)->z0()
3911 <<
"/phi0: " << (*trackParticleItr)->phi0()
3912 <<
"/theta: " << (*trackParticleItr)->theta()
3913 <<
"/qOverP: " << (*trackParticleItr)->qOverP()
3914 <<
"/vx: " << (*trackParticleItr)->vx()
3915 <<
"/vy: " << (*trackParticleItr)->vy()
3916 <<
"/vz: " << (*trackParticleItr)->vz());
3921 ATH_MSG_INFO(
"REGTEST chiSquared: " << (*trackParticleItr)->chiSquared()
3922 <<
"/numberDoF: " << (*trackParticleItr)->numberDoF());
3927 msg() << MSG::INFO <<
"REGTEST ";
3932 msg() <<
"/numberOfBLayerHits not found";
3939 msg() <<
"/numberOfPixelHits not found";
3946 msg() <<
"/numberOfPixelHoles not found";
3953 msg() <<
"/numberOfSCTHits not found";
3960 msg() <<
"/numberOfSCTHoles not found";
3967 msg() <<
"/numberOfTRTHits not found";
3974 msg() <<
"/numberOfTRTHoles not found";
3980 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrackParticle DUMP===========");
3990 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::Vertex DUMP===========");
3994 if (
sc.isFailure()) {
3996 return StatusCode::FAILURE;
4003 for (
int index = 0; vertexItr != vertexLast; ++vertexItr, ++
index) {
4008 <<
"/y: " << (*vertexItr)->y()
4009 <<
"/z: " << (*vertexItr)->z());
4012 ATH_MSG_INFO(
"REGTEST chiSquared: " << (*vertexItr)->chiSquared()
4013 <<
"/numberDoF: " << (*vertexItr)->numberDoF());
4016 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::Vertex DUMP===========");
4019 return StatusCode::SUCCESS;
4024 ATH_MSG_INFO(
"REGTEST ==========START of TDT DUMP===========" );
4027 if (
chain.empty()) {
4030 std::vector<std::string> confChains =
m_trigDec->getListOfTriggers(
chain);
4031 for (
const auto&
item : confChains) {
4034 if (
m_trigDec->getNavigationFormat() ==
"TriggerElement") {
4039 if (passFeatures.size()) {
4040 ATH_MSG_INFO(
" " <<
item <<
" Passed Final IParticle features size: " << passFeatures.size());
4042 if (!li.isValid()) {
4046 std::string state =
"ACTIVE";
4049 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);
4057 if (passAndFailFeatures.size()) {
4058 ATH_MSG_INFO(
" " <<
item <<
" Passed+Failed Final IParticle features size: " << passAndFailFeatures.size());
4060 if (!li.isValid()) {
4064 std::string state =
"ACTIVE";
4067 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);
4074 std::vector< LinkInfo<xAOD::IParticleContainer> > allFeatures =
m_trigDec->features<
xAOD::IParticleContainer>(
item, TrigDefs::includeFailedDecisions,
"", TrigDefs::allFeaturesOfType);
4075 if (allFeatures.size()) {
4076 ATH_MSG_INFO(
" " <<
item <<
" Passed+Failed ALL IParticle features size: " << allFeatures.size());
4078 if (!li.isValid()) {
4082 std::string state =
"ACTIVE";
4085 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);
4094 if (
m_trigDec->getNavigationFormat() ==
"TrigComposite") {
4096 std::vector< LinkInfo<xAOD::IParticleContainer> > muons =
m_trigDec->features<
xAOD::IParticleContainer>(
"HLT_mu24_idperf_L1MU20", TrigDefs::Physics,
"HLT_MuonL2CBInfo");
4103 ATH_MSG_INFO(
"Muon pT: " << (*
mu.link)->pt() <<
" is from the same ROI as tracks with index "
4105 <<
", which is " <<
std::distance(startIt, stopIt) <<
" tracks, out of " << muonTracksReadHandle->
size() <<
" total tracks.");
4112 ATH_MSG_INFO(
"REGTEST ==========END of TDT DUMP===========" );
4113 return StatusCode::SUCCESS;
4117 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrigCompositeContainer DUMP===========" );
4123 std::string typeNameTC;
4132 if ( not
evtStore()->contains<xAOD::TrigCompositeContainer>(
key) ) {
4141 for (
auto tc: *cont ) {
4148 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigCompositeContainer DUMP===========" );
4149 return StatusCode::SUCCESS;
4162 if (!elementLink.
isValid())
ATH_MSG_WARNING(
" Invalid element link to TrigRoiDescriptorCollection, link name:'" <<
name <<
"'");
4163 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"'' to TrigRoiDescriptor:" << *elementLink);
4170 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"' to LVL1::RecEmTauRoI:" << *elementLink);
4176 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"' to TrigComposite, TC name:'" << (*elementLink)->name() <<
"'");
4182 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"' to View:'" << *elementLink);
4184 }
else if (
name ==
"feature") {
4189 if (!elementLink.
isValid())
ATH_MSG_WARNING(
" Invalid element link to xAOD::TrigEMClusterContainer 'feature'");
4190 else ATH_MSG_DEBUG(
" Dereferenced xAOD::TrigEMClusterContainer link 'feature', Energy:" << (*elementLink)->energy());
4195 if (!elementLink.
isValid())
ATH_MSG_WARNING(
" Invalid element link to xAOD::TrigMissingETContainer 'feature'");
4196 else ATH_MSG_DEBUG(
" Dereferenced xAOD::TrigMissingETContainer link 'feature', ex:" << (*elementLink)->ex() <<
" ey:" << (*elementLink)->ey());
4203 else ATH_MSG_DEBUG(
" Dereferenced IParticle link 'feature', pt:" << (*elementLink)->pt() <<
" eta:" << (*elementLink)->eta() <<
" phi:" << (*elementLink)->phi());
4204 }
catch(std::runtime_error&
e) {
4205 ATH_MSG_WARNING(
" Cannot dereference 'feature' as IParticle: '" <<
e.what() <<
"'");
4214 return StatusCode::SUCCESS;
4225 std::vector<std::string>
keys;
4232 std::string typeNameTC;
4238 chainIDs.insert( chainID.
numeric() );
4243 pass = cg->
isPassed(TrigDefs::requireDecision);
4248 chainIDs.insert( chainID_tmp.
numeric() );
4250 if (legMultiplicites.size() == 0) {
4251 ATH_MSG_ERROR(
"chain " << chainID_tmp <<
" has invalid configuration, no multiplicity data.");
4252 }
else if (legMultiplicites.size() > 1) {
4255 for (
size_t legNumeral = 0; legNumeral < legMultiplicites.size(); ++legNumeral) {
4257 chainIDs.insert( legID.
numeric() );
4266 std::stringstream
ss;
4267 ss <<
"digraph {" << std::endl;
4268 ss <<
" node [shape=rectangle]" << std::endl;
4269 ss <<
" rankdir = BT" << std::endl;
4272 for (
const std::string&
key :
keys) {
4274 if ( not
key.starts_with(
"HLTNav_") ) {
4279 ATH_MSG_DEBUG(
"Processing collection " <<
key <<
" to be added to the navigation graph");
4281 bool writtenHeader =
false;
4282 for (
const Decision* tc : *container ) {
4286 std::vector<ElementLink<DecisionContainer>> seedELs = tc->objectCollectionLinks<
DecisionContainer>(
"seed");
4287 const bool isHypoAlgNode = tc->name() ==
"H";
4288 const bool isComboHypoAlgNode = tc->name() ==
"CH";
4289 const std::vector<DecisionID>&
decisions = tc->decisions();
4293 bool doDump =
false;
4296 if (chainIDs.count(
id) == 1) {
4304 const std::vector<DecisionID>& seedDecisions = (*s)->decisions();
4306 if (chainIDs.count(
id) == 1) {
4317 if (!writtenHeader) {
4318 writtenHeader =
true;
4319 ss <<
" subgraph " <<
key <<
" {" << std::endl;
4320 ss <<
" label=\"" <<
key <<
"\"" << std::endl;
4322 const std::string
scheme =
"rdpu9";
4323 std::string
color =
"1";
4324 if (tc->name() ==
"L1") {
color =
"1"; }
4325 else if (tc->name() ==
"F") {
color =
"2"; }
4326 else if (tc->name() ==
"IM") {
color =
"3"; }
4327 else if (tc->name() ==
"H") {
color =
"4"; }
4328 else if (tc->name() ==
"CH") {
color =
"5"; }
4329 else if (tc->name() ==
"SF") {
color =
"6"; }
4330 else if (tc->name() ==
"HLTPassRaw") {
color =
"7"; }
4331 ss <<
" \"" << selfKey <<
"_" << selfIndex <<
"\" [colorscheme="<<
scheme<<
",style=filled,fillcolor="<<
color<<
",label=<<B>Container</B>=" << typeNameTC;
4332 if (tc->name() !=
"")
ss <<
" <B>Name</B>=" << tc->name();
4333 ss <<
"<BR/><B>Key</B>=" <<
key <<
"<BR/><B>Index</B>=" << selfIndex;
4334 const bool isRemapped = tc->isRemapped();
4335 if (isHypoAlgNode)
ss <<
" <B>linksRemapped</B>=" << (isRemapped ?
"Y" :
"N");
4337 ss <<
"<BR/><B>Pass</B>=";
4341 std::string highlight = (dID.
numeric() == chainID.
numeric() ?
"<B>[CHAIN:" :
"");
4345 ss << std::hex << highlight << decisionID << (!highlight.empty() ?
"]</B>" :
"") << std::dec <<
",";
4352 ss <<
">]" << std::endl;
4354 size_t seedCount = 0;
4355 for (
size_t i = 0;
i < tc->linkColNames().
size(); ++
i) {
4356 const std::string link = tc->linkColNames().
at(
i);
4357 if (link ==
"seed" || link ==
"seed__COLL") {
4360 const uint32_t seedIndex = tc->linkColIndices().
at(
i);
4364 const std::vector<DecisionID> seedDecisions = (*seedEL)->decisions();
4365 bool doSeedLink =
false;
4367 if (chainIDs.count(
id) == 1) {
4376 ss <<
" \"" << selfKey <<
"_" << selfIndex <<
"\" -> \"" << seedKey <<
"_" << seedIndex <<
"\" [colorscheme="<<
scheme<<
",color=9,fontcolor=8,label=\"seed\"]" << std::endl;
4379 std::string linkColour =
"12";
4380 std::string linkBackground =
"11";
4381 const std::string extScheme =
"paired12";
4382 if (link ==
"roi") { linkColour=
"2"; linkBackground=
"1"; }
4383 else if (link ==
"initialRoI") { linkColour=
"2"; linkBackground=
"1"; }
4384 else if (link ==
"initialRecRoI") { linkColour=
"8"; linkBackground=
"7"; }
4385 else if (link ==
"feature") { linkColour=
"4"; linkBackground=
"3"; }
4386 else if (link ==
"view") { linkColour=
"10"; linkBackground=
"9"; }
4387 const CLID linkCLID =
static_cast<CLID>( tc->linkColClids().
at(
i) );
4393 const unsigned index = (isRemapped ? tc->linkColIndicesRemap().
at(
i) : tc->linkColIndices().
at(
i));
4396 const std::string* keyStr =
evtStore()->keyToString(
key, checkCLID);
4397 if (keyStr !=
nullptr && checkCLID != linkCLID) {
4398 std::string tnameOfCheck;
4399 m_clidSvc->getTypeNameOfID(checkCLID, tnameOfCheck).ignore();
4400 ATH_MSG_ERROR(
"Inconsistent CLID " << checkCLID <<
" [" << tnameOfCheck <<
"] stored in storegate for key " <<
key
4401 <<
". We were expecting " << linkCLID <<
" [" << tname <<
"]");
4404 std::string tnameEscape;
4405 for (std::string::const_iterator
i = tname.begin();
i != tname.end(); ++
i) {
4406 unsigned char c = *
i;
4408 tnameEscape +=
"<";
4409 }
else if (
c ==
'>') {
4410 tnameEscape +=
">";
4417 ss <<
" \"" << selfKey <<
"_" << selfIndex <<
"\" -> \"" <<
key <<
"_" <<
index <<
"\" ";
4418 ss <<
"[colorscheme="<<extScheme<<
",color="<<linkColour<<
",fontcolor="<<linkColour<<
",arrowhead=empty,label=\"" << link <<
"\"]" << std::endl;
4421 const bool linkToSelf = (selfKey ==
key and selfIndex ==
index);
4424 ss <<
" \"" <<
key <<
"_" <<
index <<
"\" [colorscheme="<<extScheme<<
",style=filled,fillcolor="<<linkBackground<<
",label=<<B>Container</B>=" << tnameEscape <<
"<BR/><B>Key</B>=";
4425 if (keyStr !=
nullptr)
ss << *keyStr;
4426 else ss <<
"[<I>KEY "<<
key <<
" NOT IN STORE</I>] ";
4427 ss <<
"<BR/><B>Index</B>=" <<
index <<
">]";
4434 if (writtenHeader) {
4435 ss <<
" }" << std::endl;
4439 ss <<
"}" << std::endl;
4441 returnValue.assign(
ss.str() );
4442 return StatusCode::SUCCESS;
4455 const std::vector< unsigned int > serialisedNavigation = navigationHandle->
serialized();
4456 ATH_MSG_INFO(
"Serialised navigation size: " << serialisedNavigation.size() );
4465 std::map< int, std::string > hash2string;
4466 for (
auto const& sequence :
testMenu->front()->sequenceInputTEs() ) {
4467 for (
auto const&
name : sequence ) {
4474 unsigned int chainCounter = 0;
4475 std::map< int, std::string > hash2chain;
4476 for (
auto const&
chain :
testMenu->front()->chainSignatureOutputTEs() ) {
4483 for (
auto const& signature :
chain ) {
4484 for (
auto const&
name : signature ) {
4496 if ( lhs.
getCLID() == rhs.getCLID() )
return ( lhs.
getIndex() < rhs.getIndex() );
4499 else return ( lhs.
getCLID() < rhs.getCLID() );
4501 std::map< HLT::TriggerElement::FeatureAccessHelper, std::vector< HLT::TriggerElement* >, decltype(cmpLambda) > feature2element(cmpLambda);
4504 std::vector< HLT::TriggerElement* > allTEs;
4505 testNav->
getAll( allTEs,
false );
4506 for (
auto element : allTEs ) {
4509 for (
auto helper : element->getFeatureAccessHelpers() ) {
4510 feature2element[
helper ].push_back( element );
4516 for (
auto element : allTEs ) {
4518 ATH_MSG_INFO(
"+++++++++++ " << hash2string[ element->getId() ] <<
" is terminal node" );
4520 std::queue< HLT::TriggerElement* > allAncestors;
4521 allAncestors.push( element );
4522 while ( allAncestors.size() ) {
4526 auto theseAncestors = thisElement->
getRelated( HLT::TriggerElement::Relation::seededByRelation );
4537 for (
auto ancestor : theseAncestors ) {
4538 allAncestors.push( ancestor );
4546 auto decisionOutput = outputNavigation.
ptr();
4549 std::map< HLT::TriggerElement const*, std::vector< int > > element2decisions;
4550 for (
const auto& pair : feature2element ) {
4553 std::string featureName = testNav->
label( pair.first.getCLID(), pair.first.getIndex().subTypeIndex() );
4554 auto sgKey =
evtStore()->stringToKey( featureName, pair.first.getCLID() );
4557 std::string storeFeatureName =
"feature";
4564 decision->typelessSetObjectLink( storeFeatureName, sgKey, pair.first.getCLID(), pair.first.getIndex().objectsBegin(), pair.first.getIndex().objectsEnd() );
4567 std::set< std::string > passedChains;
4571 passedChains.insert( hash2chain[ element->
getId() ] );
4574 int decisionNumber = decisionOutput->size() - 1;
4575 element2decisions[ element ].push_back( decisionNumber );
4579 for (
auto&
chain : passedChains ) {
4585 unsigned int decisionCounter = 0;
4586 for (
const auto& pair : feature2element ) {
4589 auto decision = decisionOutput->at( decisionCounter );
4593 for (
auto element : pair.second ) {
4594 auto theseAncestors = element->getRelated( HLT::TriggerElement::Relation::seededByRelation );
4595 for (
auto ancestor : theseAncestors ) {
4596 for (
int decisionIndex : element2decisions[ ancestor ] ) {
4603 return StatusCode::SUCCESS;