7 #include "GaudiKernel/IToolSvc.h"
110 static const int maxRepWarnings = 5;
115 m_clidSvc(
"ClassIDSvc",
name )
240 return StatusCode::SUCCESS;
264 if(m_doDumpTrackParticleContainer){
266 if (
sc.isFailure()) {
267 ATH_MSG_ERROR(
"The method dumpTrackParticleContainer() failed");
271 if(m_doDumpAll || m_doDumpLVL1_ROI ){
273 if (
sc.isFailure()) {
288 if(m_doDumpAll || m_doDumpxAODTrigMissingET){
290 if (
sc.isFailure()) {
295 if(m_doDumpAll || m_doDumpMuonFeature){
297 if (
sc.isFailure()) {
302 if(m_doDumpAll || m_doDumpCombinedMuonFeature){
304 if (
sc.isFailure()) {
305 ATH_MSG_ERROR(
"The method dumpCombinedMuonFeature() failed");
307 sc = dumpCombinedMuonFeatureContainer();
308 if (
sc.isFailure()) {
309 ATH_MSG_ERROR(
"The method dumpCombinedMuonFeatureContainer() failed");
310 return StatusCode::SUCCESS;
314 if(m_doDumpAll || m_doDumpTileMuFeature) {
316 if (
sc.isFailure()) {
317 ATH_MSG_ERROR(
"The method dumpTileMuFeatureContainer() failed");
321 if(m_doDumpAll || m_doDumpTileTrackMuFeature) {
323 if (
sc.isFailure()) {
324 ATH_MSG_ERROR(
"The method dumpTileTrackMuFeatureContainer() failed");
328 if(m_doDumpAll || m_doDumpTrigEMCluster){
330 if (
sc.isFailure()) {
335 if(m_doDumpAll || m_doDumpTrigEMClusterContainer){
337 if (
sc.isFailure()) {
338 ATH_MSG_ERROR(
"The method dumpTrigEMClusterContainer() failed");
342 if(m_doDumpAll || m_doDumpxAODTrigEMCluster){
344 if (
sc.isFailure()) {
349 if(m_doDumpAll || m_doDumpxAODTrigEMClusterContainer){
351 if (
sc.isFailure()) {
352 ATH_MSG_ERROR(
"The method dumpxAODTrigEMClusterContainer() failed");
356 if(m_doDumpTrigTauClusterContainer){
358 if (
sc.isFailure()) {
359 ATH_MSG_ERROR(
"The method dumpTrigTauClusterContainer() failed");
363 if(m_doDumpAll || m_doDumpTrigPhotonContainer){
365 if (
sc.isFailure()) {
366 ATH_MSG_ERROR(
"The method dumpTrigPhotonContainer() failed");
370 if(m_doDumpAll || m_doDumpxAODJetContainer){
372 if (
sc.isFailure()) {
377 if(m_doDumpAll || m_doDumpTrigL2BphysContainer){
379 if (
sc.isFailure()) {
380 ATH_MSG_ERROR(
"The method dumpTrigL2BphysContainer() failed");
384 if(m_doDumpAll || m_doDumpTrigEFBphysContainer){
386 if (
sc.isFailure()) {
387 ATH_MSG_ERROR(
"The method dumpTrigEFBphysContainer() failed");
391 if(m_doDumpAll || m_doDumpTrigEFBjetContainer){
393 if (
sc.isFailure()) {
394 ATH_MSG_ERROR(
"The method dumpTrigEFBjetContainer() failed");
398 if(m_doDumpAll || m_doDumpTrigL2BjetContainer){
400 if (
sc.isFailure()) {
401 ATH_MSG_ERROR(
"The method dumpTrigL2BjetContainer() failed");
405 if(m_doDumpAll || m_doDumpTrigMuonEFContainer){
407 if (
sc.isFailure()) {
408 ATH_MSG_ERROR(
"The method dumpTrigMuonEFContainer() failed");
412 if(m_doDumpAll || m_doDumpTrigMuonEFInfoContainer){
414 if (
sc.isFailure()) {
415 ATH_MSG_ERROR(
"The method dumpTrigMuonEFInfoContainer() failed");
419 if(m_doDumpAll || m_doDumpTrigMuonEFIsolationContainer) {
422 ATH_MSG_ERROR(
"The method dumpTrigMuonEFIsolationContainer() failed");
426 if(m_doDumpAll || m_doDumpxAODMuonContainer) {
433 if(m_doDumpAll || m_doDumpTrigElectronContainer){
435 if (
sc.isFailure()) {
436 ATH_MSG_ERROR(
"The method dumpTrigElectronContainer() failed");
440 if(m_doDumpAll || m_doDumpxAODTrigElectronContainer){
442 if (
sc.isFailure()) {
443 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
447 if(m_doDumpAll || m_doDumpxAODTrigPhotonContainer){
449 if (
sc.isFailure()) {
450 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
454 if(m_doDumpAll || m_doDumpxAODElectronContainer){
456 if (
sc.isFailure()) {
457 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
461 if(m_doDumpAll || m_doDumpxAODPhotonContainer){
463 if (
sc.isFailure()) {
464 ATH_MSG_ERROR(
"The method dumpxAODTrigElectronContainer() failed");
468 if(m_doDumpTrigTauContainer){
470 if (
sc.isFailure()) {
475 if(m_doDumpTrigTauTracksInfo){
477 if (
sc.isFailure()) {
482 if(m_doDumpAll || m_doDumpHLTResult){
484 if (
sc.isFailure()) {
489 if(m_doDumpTrigInDetTrackCollection){
491 if (
sc.isFailure()) {
492 ATH_MSG_ERROR(
"The method dumpTrigInDetTrackCollection() failed");
496 if(m_doDumpAll || m_doDumpTrigVertexCollection){
498 if (
sc.isFailure()) {
499 ATH_MSG_ERROR(
"The method dumpTrigVertexCollection() failed");
503 if(m_doDumpAll || m_doDumpxAODTauJetContainer){
505 if (
sc.isFailure()) {
506 ATH_MSG_ERROR(
"The method dumpxAODTauJetContainer() failed");
510 if(m_doDumpTauJetContainer){
512 if (
sc.isFailure()) {
517 if(m_doDumpAll || m_doDumpxAODTrackParticle){
519 if (
sc.isFailure()) {
524 if(m_doDumpAll || m_doDumpxAODVertex){
526 if (
sc.isFailure()) {
531 if (m_doDumpAll || m_doDumpxAODTrigMinBias){
533 if (
sc.isFailure()) {
538 if (m_doDumpTrigPassBits){
540 if (
sc.isFailure()) {
545 if (m_doDumpAll || m_doDumpStoreGate) {
549 if (m_doDumpAll || m_doDumpNavigation) {
551 if (
sc.isFailure() ) {
556 if (m_doDumpAll || m_doTDTCheck) {
560 if (m_doDumpAll || m_doDumpAllTrigComposite || m_dumpTrigCompositeContainers.size() > 0) {
564 if (m_doDumpAll || m_doDumpTrigCompsiteNavigation) {
565 std::string trigCompositeSteering;
567 ATH_CHECK(TrigCompositeNavigationToDot(trigCompositeSteering, pass));
568 const EventContext& context = Gaudi::Hive::currentContext();
569 const std::string evtNumber =
std::to_string(context.eventID().event_number());
570 const std::string passStr = (pass ?
"Pass" :
"Fail");
571 std::ofstream
ofile(std::string(
"NavGraph_" + m_dumpNavForChain +
"_Ev" + evtNumber +
"_" + passStr +
".dot").c_str());
572 ofile << trigCompositeSteering;
578 return StatusCode::SUCCESS;
584 const std::string
name=
"HLT_xAOD__TrigPassBitsContainer_passbits";
587 if (
sc.isFailure() ){
592 for(
const auto bits:*xbitscont){
597 ATH_MSG_DEBUG(
"Analyzing bits for " << bits->containerClid() <<
" of size " << bits->size() <<
" with bit size " << bits->passBits().size());
602 for (
int j=0; itr != itrE; ++itr, ++j ) {
611 return StatusCode::SUCCESS;
618 std::string METTag=
"HLT_xAOD__TrigSpacePointCountsContainer_spacepoints";
626 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << SpacePointCountsCont->
size() <<
" elements");
628 std::string
s;
char buff[128];
629 std::vector<float> getVec;
633 for(
uint i = 0;
i < SpacePointCountsCont->
size();
i++) {
634 getVec = SpacePointCountsCont->
at(
i)->contentsPixelClusEndcapC();
636 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
637 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of contentsPixelClusEndcapC() = %10.2f ",
s.c_str(),
sum );
640 getVec = SpacePointCountsCont->
at(
i)->contentsPixelClusBarrel();
642 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
643 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of contentsPixelClusBarrel() = %10.2f ",
s.c_str(),
sum );
646 getVec = SpacePointCountsCont->
at(
i)->contentsPixelClusEndcapA();
648 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
649 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of contentsPixelClusEndcapA() = %10.2f ",
s.c_str(),
sum );
652 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusTotBins() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusTotBins() );
655 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusTotMin() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusTotMin() );
658 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusTotMax() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusTotMax() );
661 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusSizeBins() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusSizeBins() );
664 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusSizeMin() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusSizeMin() );
667 snprintf(
buff,
sizeof(
buff),
"REGTEST %s pixelClusSizeMax() = %10.2f ",
s.c_str(), SpacePointCountsCont->
at(
i)->pixelClusSizeMax() );
670 snprintf(
buff,
sizeof(
buff),
"REGTEST %s sctSpEndcapC() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->sctSpEndcapC() );
673 snprintf(
buff,
sizeof(
buff),
"REGTEST %s sctSpBarrel() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->sctSpBarrel() );
676 snprintf(
buff,
sizeof(
buff),
"REGTEST %s sctSpEndcapA() = %u ",
s.c_str(), SpacePointCountsCont->
at(
i)->sctSpEndcapA() );
685 std::string METTag=
"HLT_xAOD__TrigT2MbtsBitsContainer_T2Mbts";
693 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << T2MbtsBitsCont->
size() <<
" elements");
695 std::string
s;
char buff[380];
696 std::vector<float> getVec;
701 getVec = T2MbtsBitsCont->
at(
i)->triggerEnergies();
703 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
704 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of triggerEnergies() = %10.2f ",
s.c_str(),
sum );
707 getVec = T2MbtsBitsCont->
at(
i)->triggerTimes();
709 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
710 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of triggerTimes() = %10.2f ",
s.c_str(),
sum );
719 std::string METTag=
"HLT_xAOD__TrigVertexCountsContainer_vertexcounts";
727 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << T2VertexCountsCont->
size() <<
" elements");
729 std::string
s;
char buff[380];
730 std::vector<float> fgetVec;
732 std::vector<unsigned int> ugetVec;
733 unsigned int usum(0);
736 for(
uint i = 0;
i < T2VertexCountsCont->
size();
i++) {
737 ugetVec = T2VertexCountsCont->
at(
i)->vtxNtrks();
738 for (
uint j = 0; j < ugetVec.size(); ++j) usum += ugetVec[j];
739 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of vtxNtrks() = %u ",
s.c_str(), usum );
742 fgetVec = T2VertexCountsCont->
at(
i)->vtxTrkPtSqSum();
743 for (
uint j = 0; j < fgetVec.size(); ++j) fsum += fgetVec[j];
744 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of vtxTrkPtSqSum() = %10.2f ",
s.c_str(), fsum );
753 std::string METTag=
"HLT_xAOD__TrigTrackCountsContainer_trackcounts";
761 ATH_MSG_INFO(
"Accessing " << METTag <<
" with " << T2TrackCountsCont->
size() <<
" elements");
763 std::string
s;
char buff[380];
764 std::vector<float> getVec;
768 for(
uint i = 0;
i < T2TrackCountsCont->
size();
i++) {
769 getVec = T2TrackCountsCont->
at(
i)->z0_pt();
771 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
772 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of z0_pt = %10.2f ",
s.c_str(),
sum );
775 getVec = T2TrackCountsCont->
at(
i)->eta_phi();
777 for (
uint j = 0; j < getVec.size(); ++j)
sum += getVec[j];
778 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SUM of eta_phi() = %10.2f ",
s.c_str(),
sum );
781 snprintf(
buff,
sizeof(
buff),
"REGTEST %s z0Bins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->z0Bins() );
784 snprintf(
buff,
sizeof(
buff),
"REGTEST %s z0Min() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->z0Min() );
787 snprintf(
buff,
sizeof(
buff),
"REGTEST %s z0Max() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->z0Max() );
790 snprintf(
buff,
sizeof(
buff),
"REGTEST %s ptBins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->ptBins() );
793 snprintf(
buff,
sizeof(
buff),
"REGTEST %s ptMin() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->ptMin() );
796 snprintf(
buff,
sizeof(
buff),
"REGTEST %s ptMax() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->ptMax() );
799 snprintf(
buff,
sizeof(
buff),
"REGTEST %s etaBins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->etaBins() );
802 snprintf(
buff,
sizeof(
buff),
"REGTEST %s etaMin() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->etaMin() );
805 snprintf(
buff,
sizeof(
buff),
"REGTEST %s etaMax() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->etaMax() );
808 snprintf(
buff,
sizeof(
buff),
"REGTEST %s phiBins() = %u ",
s.c_str(), T2TrackCountsCont->
at(
i)->phiBins() );
811 snprintf(
buff,
sizeof(
buff),
"REGTEST %s phiMin() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->phiMin() );
814 snprintf(
buff,
sizeof(
buff),
"REGTEST %s phiMax() = %10.2f ",
s.c_str(), T2TrackCountsCont->
at(
i)->phiMax() );
827 return StatusCode::SUCCESS;
839 std::string METTags[]={
"HLT_xAOD__TrigMissingETContainer_EFJetEtSum",
"HLT_xAOD__TrigMissingETContainer_TrigEFMissingET",
"HLT_xAOD__TrigMissingETContainer_TrigL2MissingET_FEB",
"HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_topocl"};
841 for(
int itag=0; itag <ntag; itag++) {
848 ATH_MSG_INFO(
"Accessing " << METTags[itag] <<
" with " << MissingETCont->
size() <<
" elements");
853 std::string
s;
char buff[3000];
855 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ex = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->ex() );
857 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ey = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->ey() );
859 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ez = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->ez() );
861 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumET = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->sumEt() );
863 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumE = %10.2f CLHEP::MeV",
s.c_str(), MissingETCont->
at(
i)->sumE() );
865 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Flag = %d",
s.c_str(), MissingETCont->
at(
i)->flag() );
867 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Flag = %d",
s.c_str(), MissingETCont->
at(
i)->roiWord() );
870 unsigned int Nc = MissingETCont->
at(
i)->getNumberOfComponents();
872 s=
"REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
873 s+=
"/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
877 for(
uint j = 0; j < Nc; j++) {
879 std::string
name = MissingETCont->
at(
i)->nameOfComponent(j);
880 const short status = MissingETCont->
at(
i)->statusComponent(j);
881 const unsigned short usedChan = MissingETCont->
at(
i)->usedChannelsComponent(j);
882 const short sumOfSigns = MissingETCont->
at(
i)->sumOfSignsComponent(j);
883 const float calib0 = MissingETCont->
at(
i)->calib0Component(j);
884 const float calib1 = MissingETCont->
at(
i)->calib1Component(j);
885 const float ex = MissingETCont->
at(
i)->exComponent(j);
886 const float ey = MissingETCont->
at(
i)->eyComponent(j);
887 const float ez = MissingETCont->
at(
i)->ezComponent(j);
888 const float sumE = MissingETCont->
at(
i)->sumEComponent(j);
889 const float sumEt = MissingETCont->
at(
i)->sumEtComponent(j);
892 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
893 name.c_str(),
status, usedChan, sumOfSigns, calib1, calib0,
894 ex, ey, ez, sumE, sumEt);
901 return StatusCode::SUCCESS;
913 std::string METTags[]={
"HLT_TrigMissingETContainer_TrigEFMissingET",
"HLT_TrigMissingETContainer_TrigEFMissingET_FEB",
"HLT_TrigMissingETContainer_TrigEFMissingET_topocl"};
916 for (
int itag=0; itag < ntag; itag++) {
918 StatusCode sc=evtStore()->retrieve(trigMETcont , METTags[itag]);
919 if(
sc.isFailure() ){
920 ATH_MSG_INFO(
"Failed to retrieve TrigMissingETContainer with key " << METTags[itag]);
924 ATH_MSG_INFO(
"Got TrigMissingETContainer with key \"" << METTags[itag]<<
"\"");
929 for (
int j=0; trigMETfirst != trigMETlast; ++trigMETfirst++, ++j ) {
931 ATH_MSG_INFO(
"REGTEST ==========START of TrigMissingET DUMP===========");
936 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ex = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->ex() );
938 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ey = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->ey() );
940 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Ez = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->ez() );
942 snprintf(
buff,
sizeof(
buff),
"REGTEST %s Et = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->et() );
944 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumEt = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->sumEt() );
946 snprintf(
buff,
sizeof(
buff),
"REGTEST %s SumE = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->sumE() );
948 snprintf(
buff,
sizeof(
buff),
"REGTEST %s E = %10.2f CLHEP::MeV",
s.c_str(), (*trigMETfirst)->e() );
950 snprintf(
buff,
sizeof(
buff),
"REGTEST %s flag = %10d",
s.c_str(), (*trigMETfirst)->getFlag() );
952 snprintf(
buff,
sizeof(
buff),
"REGTEST %s RoIword = %10ld",
s.c_str(), (*trigMETfirst)->RoIword() );
955 unsigned int Nc = (*trigMETfirst)->getNumOfComponents();
957 s=
"REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
958 s+=
"/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
961 for (
unsigned int i=0;
i<Nc; ++
i) {
962 std::string
name = (*trigMETfirst)->getNameOfComponent(
i);
963 const short status = (*trigMETfirst)->getStatus(
i);
964 const unsigned short usedChan = (*trigMETfirst)->getUsedChannels(
i);
965 const short sumOfSigns = (*trigMETfirst)->getSumOfSigns(
i);
966 const float calib0 = (*trigMETfirst)->getComponentCalib0(
i);
967 const float calib1 = (*trigMETfirst)->getComponentCalib1(
i);
968 const float ex = (*trigMETfirst)->getExComponent(
i);
969 const float ey = (*trigMETfirst)->getEyComponent(
i);
970 const float ez = (*trigMETfirst)->getEzComponent(
i);
971 const float sumE = (*trigMETfirst)->getSumEComponent(
i);
972 const float sumEt = (*trigMETfirst)->getSumEtComponent(
i);
975 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
976 name.c_str(),
status, usedChan, sumOfSigns, calib1, calib0,
977 ex, ey, ez, sumE, sumEt);
987 ATH_MSG_INFO(
"Trying to fetch TrigMissingET objects from older releases");
990 StatusCode sc=evtStore()->retrieve(trigMETfirst ,trigMETlast);
991 if(
sc.isFailure() ){
992 ATH_MSG_INFO(
"Failed to retrieve TrigMissingET (rel. <= 14.2.0)");
995 for( ; trigMETfirst != trigMETlast ; ++trigMETfirst ){
996 const std::string&
name(trigMETfirst.
key());
1002 if(
name.find(
"TrigEFMissingET") != std::string::npos ) {
1004 }
else if(
name.find(
"T2MissingET") != std::string::npos ){
1011 snprintf(
buff,
sizeof(
buff),
"%s Ex = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->ex() );
1013 snprintf(
buff,
sizeof(
buff),
"%s Ey = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->ey() );
1015 snprintf(
buff,
sizeof(
buff),
"%s Ez = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->ez() );
1017 snprintf(
buff,
sizeof(
buff),
"%s Et = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->et() );
1019 snprintf(
buff,
sizeof(
buff),
"%s SumE = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->sumE() );
1021 snprintf(
buff,
sizeof(
buff),
"%s SumEt = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->sumEt() );
1023 snprintf(
buff,
sizeof(
buff),
"%s E = %10.2f CLHEP::MeV",
s.c_str(), trigMETfirst->e() );
1025 snprintf(
buff,
sizeof(
buff),
"%s flag = %10d",
s.c_str(), trigMETfirst->getFlag() );
1027 snprintf(
buff,
sizeof(
buff),
"%s RoIword = %10ld",
s.c_str(), trigMETfirst->RoIword() );
1030 unsigned int Nc = trigMETfirst->getNumOfComponents();
1032 s=
"REGTEST __name____status_usedChannels__sumOfSigns__calib1_calib0";
1033 s+=
"/MeV__ex/MeV_____ey/MeV_____ez/MeV___sumE/MeV__sumEt/CLHEP::MeV";
1036 for (
unsigned int i=0;
i<Nc; ++
i) {
1037 std::string
name = trigMETfirst->getNameOfComponent(
i);
1038 const short status = trigMETfirst->getStatus(
i);
1039 const unsigned short usedChan = trigMETfirst->getUsedChannels(
i);
1040 const short sumOfSigns = trigMETfirst->getSumOfSigns(
i);
1041 const float calib0 = trigMETfirst->getComponentCalib0(
i);
1042 const float calib1 = trigMETfirst->getComponentCalib1(
i);
1043 const float ex = trigMETfirst->getExComponent(
i);
1044 const float ey = trigMETfirst->getEyComponent(
i);
1045 const float ez = trigMETfirst->getEzComponent(
i);
1046 const float sumE = trigMETfirst->getSumEComponent(
i);
1047 const float sumEt = trigMETfirst->getSumEtComponent(
i);
1050 "REGTEST %s %6d %12d %10d %6.2f %6.3f %10.2f %10.2f %10.2f %10.2f %10.2f",
1051 name.c_str(),
status, usedChan, sumOfSigns, calib1, calib0,
1052 ex, ey, ez, sumE, sumEt);
1058 ATH_MSG_INFO(
"REGTEST ==========END of TrigMissingET DUMP===========");
1060 return StatusCode::SUCCESS;
1069 ATH_MSG_INFO(
"REGTEST ==========START of TrackParticleContainer DUMP===========");
1071 std::string trackPtags[]={
"HLT_InDetTrigParticleCreation_Bjet_EFID",
1072 "HLT_InDetTrigParticleCreation_Bphysics_EFID",
1073 "HLT_InDetTrigParticleCreation_Electron_EFID",
1074 "HLT_InDetTrigParticleCreation_FullScan_EFID",
1075 "HLT_InDetTrigParticleCreation_Muon_EFID",
1076 "HLT_InDetTrigParticleCreation_Photon_EFID",
1077 "HLT_InDetTrigParticleCreation_Tau_EFID"};
1083 for (
int itag=0; itag<ntag; itag++){
1086 if (
sc.isFailure()) {
1087 ATH_MSG_INFO(
"REGTEST No TrackParticleContainer found with tag " << trackPtags[itag]);
1090 ATH_MSG_INFO(
"TrackParticleContainer found with tag " << trackPtags[itag]
1091 <<
" and size " << pTrackParticleC->
size());
1095 for (
int ind=1; trackItr != trackItrE; ++trackItr,
ind++) {
1098 << trackParticle->
charge() <<
" p "
1099 << trackParticle->
p()<<
" eta " << trackParticle->
eta()
1100 <<
" phi " << trackParticle->
phi());
1108 const auto& parameterVector = perigee->parameters();
1111 " phi " <<parameterVector[
Trk::phi]);
1120 ATH_MSG_WARNING(
" Max attached track warning reached, no further warnings given");
1131 ATH_MSG_INFO(
" vertex position (" << position[0] <<
", " <<
1132 position[1] <<
", " << position[2] <<
") ");
1137 ATH_MSG_WARNING(
" Max attached vertex warning reached, no further warnings given");
1145 const auto&
parameters = perigee->parameters();
1176 ATH_MSG_INFO(
"REGTEST ==========START of LVL1_ROI DUMP===========");
1180 if (
sc.isFailure() ) {
1182 return StatusCode::SUCCESS;
1187 LVL1_ROI::emtaus_type::const_iterator itEMTau =
1189 LVL1_ROI::emtaus_type::const_iterator itEMTau_e =
1192 for( ; itEMTau != itEMTau_e; ++itEMTau, ++j) {
1194 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI Eta is " << itEMTau->getEta());
1195 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI Phi is " << itEMTau->getPhi());
1196 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI EMClus is " << itEMTau->getEMClus());
1197 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI TauClus is " << itEMTau->getTauClus());
1198 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI EMIsol is " << itEMTau->getEMIsol());
1199 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI HadIsol is " << itEMTau->getHadIsol());
1200 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI Core is " << itEMTau->getCore());
1201 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI HadCore is " << itEMTau->getHadCore());
1202 ATH_MSG_INFO(
"REGTEST LVL1 EmTauROI roiWord is " << itEMTau->getROIWord());
1205 ATH_MSG_INFO(
"REGTEST ==========END of LVL1_ROI DUMP===========");
1207 return StatusCode::SUCCESS;
1217 ATH_MSG_INFO(
"REGTEST ==========START of TrigPhotonContainer DUMP===========");
1222 StatusCode sc = evtStore()->retrieve(trigPhoton,lastTrigPhoton);
1223 if (
sc.isFailure()) {
1225 return StatusCode::SUCCESS;
1227 ATH_MSG_INFO(
"REGTEST TrigPhotonContainers retrieved");
1230 for (
int i=0; trigPhoton != lastTrigPhoton; ++trigPhoton, ++
i) {
1232 ATH_MSG_INFO(
"REGTEST Looking at TrigPhotonContainer " <<
i);
1237 for (
int j=0; PhotonItr != PhotonItrE; ++PhotonItr, ++j ) {
1241 ATH_MSG_INFO(
"REGTEST TrigPhoton->dPhi() returns " << (*PhotonItr)->dPhi());
1242 ATH_MSG_INFO(
"REGTEST TrigPhoton->dEta() returns " << (*PhotonItr)->dEta());
1243 ATH_MSG_INFO(
"REGTEST TrigPhoton->rCore() returns " << (*PhotonItr)->rCore());
1244 ATH_MSG_INFO(
"REGTEST TrigPhoton->isValid() returns " << (*PhotonItr)->isValid());
1245 ATH_MSG_INFO(
"REGTEST TrigPhoton->Et() returns " << (*PhotonItr)->Et());
1247 if ( (*PhotonItr)->cluster() != NULL ) {
1249 ATH_MSG_INFO(
"REGTEST cluster->e() returns " << (*PhotonItr)->cluster()->e());
1250 ATH_MSG_INFO(
"REGTEST cluster->e237() returns " << (*PhotonItr)->cluster()->e237());
1251 ATH_MSG_INFO(
"REGTEST cluster->e277() returns " << (*PhotonItr)->cluster()->e277());
1252 ATH_MSG_INFO(
"REGTEST cluster->fracs1() returns " << (*PhotonItr)->cluster()->fracs1());
1253 ATH_MSG_INFO(
"REGTEST cluster->weta2() returns " << (*PhotonItr)->cluster()->weta2());
1254 ATH_MSG_INFO(
"REGTEST cluster->ehad() returns " << (*PhotonItr)->cluster()->ehad1());
1255 ATH_MSG_INFO(
"REGTEST cluster->emaxs1() returns " << (*PhotonItr)->cluster()->emaxs1());
1259 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->Pt() returns " << (*PhotonItr)->pt());
1260 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->Eta() returns " << (*PhotonItr)->eta());
1261 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->Phi() returns " << (*PhotonItr)->phi());
1262 ATH_MSG_INFO(
"REGTEST P4PtEtaPhiM->m() returns " << (*PhotonItr)->m());
1273 for (
int k=0; PhotonItr2 != PhotonItr2E; ++PhotonItr2, ++
k ) {
1275 if ( (**PhotonItr) == (**PhotonItr2) ) {
1276 ATH_MSG_INFO(
"REGTEST TrigPhoton nr. " << j <<
" is the same as TrigPhoton nr. " <<
k);
1279 std::map<std::string, double> v_diff;
1280 diff(*(*PhotonItr),*(*PhotonItr2), v_diff);
1282 ATH_MSG_INFO(
"TrigPhoton " <<
k <<
" different form TrigPhoton " << j <<
" :");
1283 for (
int m=0;
it !=v_diff.end();++
it, ++
m) {
1284 ATH_MSG_INFO(
"TrigPhoton Delta_" << (*it).first <<
" = " << (*it).second);
1291 ATH_MSG_INFO(
"REGTEST ==========END of TrigPhotonContainer DUMP===========");
1294 return StatusCode::SUCCESS;
1304 ATH_MSG_INFO(
"REGTEST ==========START of TrigMuonEFContainer DUMP===========");
1309 StatusCode sc = evtStore()->retrieve(trigMuon,lastTrigMuon);
1310 if (
sc.isFailure()) {
1312 return StatusCode::SUCCESS;
1314 ATH_MSG_INFO(
"REGTEST TrigMuonEFContainers retrieved");
1316 for (
int i=0; trigMuon != lastTrigMuon; ++trigMuon, ++
i) {
1318 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFContainer " <<
i);
1323 for (
int j=0; MuonItr != MuonItrE; ++MuonItr, ++j ) {
1327 ATH_MSG_INFO(
"REGTEST TrigMuonEF->muonCode() returns " << (*MuonItr)->MuonCode());
1328 ATH_MSG_INFO(
"REGTEST TrigMuonEF->roi() returns " << (*MuonItr)->RoINum());
1329 ATH_MSG_INFO(
"REGTEST TrigMuonEF->charge() returns " << (*MuonItr)->Charge());
1331 ATH_MSG_INFO(
"REGTEST Looking at P4IPtCotThPhiM " << j);
1333 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->iPt() returns " << (*MuonItr)->iPt());
1334 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->CotTh() returns " << (*MuonItr)->cotTh());
1335 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->Phi() returns " << (*MuonItr)->phi());
1336 ATH_MSG_INFO(
"REGTEST P4IPtCotThPhiM->m() returns " << (*MuonItr)->m());
1340 ATH_MSG_INFO(
"REGTEST ==========END of TrigMuonEFContainer DUMP===========");
1343 return StatusCode::SUCCESS;
1350 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::MuonContainer DUMP===========" );
1354 if (
sc.isFailure()) {
1355 ATH_MSG_INFO(
"REGTEST No muon container HLT_xAOD__MuonContainer_MuonEFInfo");
1356 return StatusCode::SUCCESS;
1362 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::MuonContainer DUMP===========" );
1364 return StatusCode::SUCCESS;
1374 ATH_MSG_INFO(
"REGTEST ==========START of TrigMuonEFInfoContainer DUMP===========");
1379 StatusCode sc = evtStore()->retrieve(trigMuon,lastTrigMuon);
1380 if (
sc.isFailure()) {
1381 ATH_MSG_INFO(
"REGTEST No TrigMuonEFInfoContainer found");
1382 return StatusCode::SUCCESS;
1384 ATH_MSG_INFO(
"REGTEST TrigMuonEFInfoContainers retrieved");
1387 for (
int i=0; trigMuon != lastTrigMuon; ++trigMuon, ++
i) {
1389 ATH_MSG_INFO(
"REGTEST -+-+-+-+ Looking at TrigMuonEFInfoContainer " <<
i);
1394 for (
int j=0; MuonItr != MuonItrE; ++MuonItr, ++j ) {
1396 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFInfo " << j);
1401 if (*muonInfo == *muonInfo) {
1423 printMuonTrk(muonTrack);
1425 ATH_MSG_INFO(
"REGTEST no SpectrometerTrack track found");
1428 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFTrack ExtrapolatedTrack()");
1431 printMuonTrk(muonTrack);
1433 ATH_MSG_INFO(
"REGTEST no ExtrapolatedTrack track found");
1436 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFTrack CombinedTrack()");
1439 printMuonTrk(muonTrack);
1449 return StatusCode::SUCCESS;
1451 ATH_MSG_INFO(
"REGTEST TrigMuonEFInfo->RoINum() returns " << muonInfo->RoINum());
1455 ATH_MSG_INFO(
"REGTEST ==========END of TrigMuonEFInfoContainer DUMP===========");
1458 return StatusCode::SUCCESS;
1467 ATH_MSG_INFO(
"REGTEST ==========START of TrigMuonEFIsolationContainer DUMP===========");
1472 StatusCode sc = evtStore()->retrieve(trigMuon,lastTrigMuon);
1473 if (
sc.isFailure()) {
1474 ATH_MSG_INFO(
"REGTEST No TrigMuonEFIsolationContainer found");
1475 return StatusCode::SUCCESS;
1477 ATH_MSG_INFO(
"REGTEST TrigMuonEFIsolationContainers retrieved");
1480 for (
int i=0; trigMuon != lastTrigMuon; ++trigMuon, ++
i) {
1482 ATH_MSG_INFO(
"REGTEST -+-+-+-+ Looking at TrigMuonEFIsolationContainer " <<
i);
1487 for (
int j=0; MuonItr != MuonItrE; ++MuonItr, ++j ) {
1489 ATH_MSG_INFO(
"REGTEST Looking at TrigMuonEFIsolation " << j);
1494 if (*muonIsolation == *muonIsolation) {
1512 ATH_MSG_WARNING(
"REGTEST No InfoTrack attached to this EFIsolation object");
1527 ATH_MSG_INFO(
"REGTEST No Linked MuonEFInfo, expected for pre 2012 data");
1533 ATH_MSG_INFO(
"REGTEST ==========END of TrigMuonEFIsolationContainer DUMP===========");
1534 ATH_MSG_DEBUG(
"dumpTrigMuonEFIsolationContainer() succeeded");
1536 return StatusCode::SUCCESS;
1542 ATH_MSG_INFO(
"POINTER TrigMuonEFTrack: " << muonTrack);
1568 ATH_MSG_INFO(
" POINTER TrigMuonEFCbTrack: " << muonTrack);
1585 ATH_MSG_INFO(
"REGTEST ==========START of TrigElectronContainer DUMP===========");
1590 StatusCode sc = evtStore()->retrieve(trigElec,lastTrigElec);
1591 if (
sc.isFailure()) {
1592 ATH_MSG_INFO(
"REGTEST No TrigelectronContainer found");
1593 return StatusCode::SUCCESS;
1595 ATH_MSG_INFO(
"REGTEST TrigElectronContainers retrieved");
1600 for (
int i=0; trigElec != lastTrigElec; ++trigElec, ++
i) {
1602 ATH_MSG_INFO(
"REGTEST Looking at TrigElectronContainer " <<
i);
1607 for (
int j=0; elecItr != elecItrE; ++elecItr, ++j ) {
1644 for (
int k=0; elecItr2 != elecItr2E; ++elecItr2, ++
k ) {
1646 if ( (**elecItr) == (**elecItr2) ) {
1647 ATH_MSG_INFO(
"REGTEST TrigElectron nr. " << j <<
" is the same as TrigElectron nr. " <<
k);
1650 std::map<std::string, double> v_diff;
1651 diff(*(*elecItr),*(*elecItr2), v_diff);
1653 ATH_MSG_INFO(
"TrigElectron " <<
k <<
" different form TrigElectron " << j <<
" :");
1654 for (
int m=0;
it !=v_diff.end();++
it, ++
m) {
1655 ATH_MSG_INFO(
"TrigElectron Delta_" << (*it).first <<
" = " << (*it).second);
1662 ATH_MSG_INFO(
"REGTEST ==========END of TrigElectronContainer DUMP===========");
1665 return StatusCode::SUCCESS;
1673 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrigElectronContainer DUMP===========" );
1677 if (
sc.isFailure()) {
1678 ATH_MSG_INFO(
"REGTEST No TrigElectron container HLT_xAOD__TrigElectronContainer_L2ElectronFex");
1679 return StatusCode::SUCCESS;
1682 for (
const auto eg : *elCont){
1683 ATH_MSG_INFO(
"REGTEST TrigElectron->Phi() returns " <<
eg->phi());
1684 ATH_MSG_INFO(
"REGTEST TrigElectron->Eta() returns " <<
eg->eta());
1685 ATH_MSG_INFO(
"REGTEST TrigElectron->rEta returns " <<
eg->rcore());
1686 ATH_MSG_INFO(
"REGTEST TrigElectron->eratio() returns " <<
eg->eratio());
1688 ATH_MSG_INFO(
"REGTEST TrigElectron->etHad() returns " <<
eg->etHad());
1690 ATH_MSG_INFO(
"REGTEST TrigElectron caloEta = " <<
eg->caloEta());
1691 ATH_MSG_INFO(
"REGTEST TrigElectron dPhiCalo" <<
eg->trkClusDphi());
1692 ATH_MSG_INFO(
"REGTEST TrigElectron dEtaCalo" <<
eg->trkClusDeta());
1696 ATH_MSG_INFO(
"REGTEST TrigElectron nStrawHits" <<
eg->nTRTHiThresholdHits());
1698 if(
eg->emCluster()){
1699 ATH_MSG_INFO(
"REGTEST TrigElectron EMCluster retrieved");
1700 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster->energy() returns " <<
eg->emCluster()->energy());
1701 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster->phi() returns " <<
eg->emCluster()->phi());
1702 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster->eta() returns " <<
eg->emCluster()->eta());
1703 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster check Element Link");
1704 ATH_MSG_INFO(
"REGTEST TrigElectron emCluster energy = " <<
eg->emCluster()->energy());
1705 ATH_MSG_INFO(
"REGTEST TrigElectron ElementLink emCluster energy = " << (*
eg->emClusterLink())->energy());
1707 else ATH_MSG_INFO(
"REGTEST TrigElectron No EMCluster retrieved!");
1708 ATH_MSG_INFO(
"REGTEST TrigElectron Check TrackParticle");
1709 if(
eg->trackParticle()){
1710 ATH_MSG_INFO(
"REGTEST TrigElectron TrackParticle retrieved");
1711 ATH_MSG_INFO(
"REGTEST TrigElectron trackParticle->pt() returns " <<
eg->trackParticle()->pt());
1712 ATH_MSG_INFO(
"REGTEST TrigElectron trackParticle->phi() returns " <<
eg->trackParticle()->phi());
1713 ATH_MSG_INFO(
"REGTEST TrigElectron trackParticle->eta() returns " <<
eg->trackParticle()->eta());
1714 ATH_MSG_INFO(
"REGTEST TrigElectron check TrackParticle Element Link");
1715 ATH_MSG_INFO(
"REGTEST TrigElectron TrackParticle pt = " <<
eg->trackParticle()->pt());
1716 ATH_MSG_INFO(
"REGTEST TrigElectron ElementLink TrackParticle pt = " << (*
eg->trackParticleLink())->pt());
1718 else ATH_MSG_INFO(
"REGTEST TrigElectron No TrackParticle retrieved!");
1720 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigElectronContainer DUMP===========" );
1722 return StatusCode::SUCCESS;
1731 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrigPhotonContainer DUMP===========" );
1735 if (
sc.isFailure()) {
1736 ATH_MSG_INFO(
"REGTEST No TrigPhoton container HLT_xAOD__TrigPhotonContainer_L2PhotonFex");
1737 return StatusCode::SUCCESS;
1740 for (
const auto eg : *phCont){
1744 ATH_MSG_INFO(
"REGTEST TrigPhoton->dPhi() returns " <<
eg->dPhi());
1745 ATH_MSG_INFO(
"REGTEST TrigPhoton->dEta() returns " <<
eg->dEta());
1746 ATH_MSG_INFO(
"REGTEST TrigPhoton->rEta returns " <<
eg->rcore());
1747 ATH_MSG_INFO(
"REGTEST TrigPhoton->eratio() returns " <<
eg->eratio());
1749 ATH_MSG_INFO(
"REGTEST TrigPhoton->etHad() returns " <<
eg->etHad());
1752 if(
eg->emCluster()){
1753 ATH_MSG_INFO(
"REGTEST TrigPhoton EMCluster retrieved");
1754 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster->energy() returns " <<
eg->emCluster()->energy());
1755 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster->phi() returns " <<
eg->emCluster()->phi());
1756 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster->eta() returns " <<
eg->emCluster()->eta());
1757 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster check Element Link");
1758 ATH_MSG_INFO(
"REGTEST TrigPhoton emCluster energy = " <<
eg->emCluster()->energy());
1759 ATH_MSG_INFO(
"REGTEST TrigPhoton ElementLink emCluster energy = " << (*
eg->emClusterLink())->energy());
1761 else ATH_MSG_INFO(
"REGTEST TrigPhoton No EMCluster retrieved!");
1763 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigPhotonContainer DUMP===========" );
1765 return StatusCode::SUCCESS;
1774 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::ElectronContainer DUMP===========" );
1778 if (
sc.isFailure()) {
1779 ATH_MSG_INFO(
"REGTEST No Electron container HLT_xAOD__ElectronContainer_egamma_Electrons");
1780 return StatusCode::SUCCESS;
1782 float val_float=-99;
1783 unsigned int isEMbit=0;
1786 ATH_MSG_INFO(
" REGTEST: xAOD Reconstruction variables: ");
1791 for (
const auto eg : *elCont){
1797 if(
eg->selectionisEM(isEMbit,
"isEMVLoose"))
1798 ATH_MSG_INFO(
" REGTEST: isEMVLoose " << std::hex << isEMbit << std::dec);
1800 if(
eg->selectionisEM(isEMbit,
"isEMLoose"))
1801 ATH_MSG_INFO(
" REGTEST: isEMLoose " << std::hex << isEMbit << std::dec);
1803 if(
eg->selectionisEM(isEMbit,
"isEMMedium"))
1804 ATH_MSG_INFO(
" REGTEST: isEMMedium " << std::hex << isEMbit << std::dec);
1806 if(
eg->selectionisEM(isEMbit,
"isEMTight"))
1807 ATH_MSG_INFO(
" REGTEST: isEMTight " << std::hex << isEMbit << std::dec);
1809 if(
eg->selectionisEM(isEMbit,
"isEMLHVLoose"))
1810 ATH_MSG_INFO(
" REGTEST: isEMLHVLoose " << std::hex << isEMbit << std::dec);
1812 if(
eg->selectionisEM(isEMbit,
"isEMLHLoose"))
1813 ATH_MSG_INFO(
" REGTEST: isEMLHLoose " << std::hex << isEMbit << std::dec);
1815 if(
eg->selectionisEM(isEMbit,
"isEMLHMedium"))
1816 ATH_MSG_INFO(
" REGTEST: isEMLHMedium " << std::hex << isEMbit << std::dec);
1818 if(
eg->selectionisEM(isEMbit,
"isEMLHTight"))
1819 ATH_MSG_INFO(
" REGTEST: isEMLHTight " << std::hex << isEMbit << std::dec);
1829 if(
eg->passSelection(
pid,
"LHVLoose"))
1833 if(
eg->passSelection(
pid,
"LHLoose"))
1837 if(
eg->passSelection(
pid,
"LHMedium"))
1841 if(
eg->passSelection(
pid,
"LHTight"))
1846 ATH_MSG_INFO(
" REGTEST: problems with egamma pointer" );
1847 return StatusCode::SUCCESS;
1850 if (
eg->caloCluster()) {
1851 ATH_MSG_INFO(
" REGTEST: egamma cluster transverse energy: " <<
eg->caloCluster()->et() );
1852 ATH_MSG_INFO(
" REGTEST: egamma cluster eta: " <<
eg->caloCluster()->eta() );
1853 ATH_MSG_INFO(
" REGTEST: egamma cluster phi: " <<
eg->caloCluster()->phi() );
1854 double tmpeta = -999.;
1855 double tmpphi = -999.;
1858 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. etaCalo = " << tmpeta);
1859 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. phiCalo = " << tmpphi);
1861 ATH_MSG_INFO(
" REGTEST: problems with egamma cluster pointer" );
1863 ATH_MSG_INFO(
"REGTEST: Check the original (uncalibrated)");
1865 if (!orig.
isAvailable(*
eg->caloCluster()) || !orig(*
eg->caloCluster()).isValid()){
1870 ATH_MSG_INFO(
"REGTEST:: Compare new and old clusters");
1871 ATH_MSG_INFO(
"REGTEST:: Original Cluster e,eta,phi" << origClus->
e() <<
" " << origClus->
eta() <<
" " << origClus->
phi());
1872 ATH_MSG_INFO(
"REGTEST:: MVA Cluster e,eta,phi" <<
eg->caloCluster()->e() <<
" " <<
eg->caloCluster()->eta() <<
" " <<
eg->caloCluster()->phi());
1875 if(
eg->trackParticle()){
1877 ATH_MSG_INFO(
" REGTEST: charge= " <<
eg->trackParticle()->charge());
1878 ATH_MSG_INFO(
" REGTEST: E/p= " <<
eg->caloCluster()->et() /
eg->trackParticle()->pt() );
1880 ATH_MSG_INFO(
" REGTEST: Delta eta 1st sampling= " << val_float);
1882 ATH_MSG_INFO(
" REGTEST: Delta phi 2nd sampling= " << val_float);
1884 ATH_MSG_INFO(
" REGTEST: no electron eg->trackParticle() pointer");
1919 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::ElectronContainer DUMP===========" );
1921 return StatusCode::SUCCESS;
1930 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::PhotonContainer DUMP===========" );
1934 if (
sc.isFailure()) {
1935 ATH_MSG_INFO(
"REGTEST No Photon container HLT_xAOD__PhotonContainer_egamma_Photons");
1936 return StatusCode::SUCCESS;
1939 float val_float=-99;
1940 unsigned int isEMbit=0;
1942 ATH_MSG_INFO(
" REGTEST: xAOD Reconstruction variables: ");
1944 for (
const auto eg : *phCont){
1950 ATH_MSG_INFO(
" REGTEST: isEMLoose " <<
eg->selectionisEM(isEMbit,
"isEMLoose"));
1951 ATH_MSG_INFO(
" REGTEST: isEMLoose bit " << std::hex << isEMbit << std::dec);
1952 ATH_MSG_INFO(
" REGTEST: isEMMedium " <<
eg->selectionisEM(isEMbit,
"isEMMedium"));
1953 ATH_MSG_INFO(
" REGTEST: isEMMedium bit " << std::hex << isEMbit << std::dec);
1954 ATH_MSG_INFO(
" REGTEST: isEMTight " <<
eg->selectionisEM(isEMbit,
"isEMTight"));
1955 ATH_MSG_INFO(
" REGTEST: isEMTight bit " << std::hex << isEMbit << std::dec);
1957 ATH_MSG_INFO(
" REGTEST: problems with egamma pointer" );
1958 return StatusCode::SUCCESS;
1961 if (
eg->caloCluster()) {
1962 ATH_MSG_INFO(
" REGTEST: egamma cluster transverse energy: " <<
eg->caloCluster()->et() );
1963 ATH_MSG_INFO(
" REGTEST: egamma cluster eta: " <<
eg->caloCluster()->eta() );
1964 ATH_MSG_INFO(
" REGTEST: egamma cluster phi: " <<
eg->caloCluster()->phi() );
1965 double tmpeta = -999.;
1966 double tmpphi = -999.;
1969 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. etaCalo = " << tmpeta);
1970 ATH_MSG_INFO(
" REGTEST: egamma cluster calo-frame coords. phiCalo = " << tmpphi);
1972 ATH_MSG_INFO(
" REGTEST: problems with egamma cluster pointer" );
1974 ATH_MSG_INFO(
"REGTEST: Check the original (uncalibrated)");
1976 if (!orig.
isAvailable(*
eg->caloCluster()) || !orig(*
eg->caloCluster()).isValid()){
1981 ATH_MSG_INFO(
"REGTEST:: Compare new and old clusters");
1982 ATH_MSG_INFO(
"REGTEST:: Original Cluster e,eta,phi" << origClus->
e() <<
" " << origClus->
eta() <<
" " << origClus->
phi());
1983 ATH_MSG_INFO(
"REGTEST:: MVA Cluster e,eta,phi" <<
eg->caloCluster()->e() <<
" " <<
eg->caloCluster()->eta() <<
" " <<
eg->caloCluster()->phi());
2011 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::PhotonContainer DUMP===========" );
2013 return StatusCode::SUCCESS;
2021 ATH_MSG_INFO(
"REGTEST ==========START of TrigTauContainer DUMP===========");
2026 StatusCode sc = evtStore()->retrieve(trigTau,lastTrigTau);
2027 if (
sc.isFailure()) {
2029 return StatusCode::SUCCESS;
2033 for (
int i=0; trigTau != lastTrigTau; ++trigTau, ++
i) {
2040 for (
int j=0; TrigTauItr != TrigTauItrE; ++TrigTauItr, ++j ) {
2044 ATH_MSG_INFO(
"REGTEST TrigTau->pt() returns " << (*TrigTauItr)->pt());
2045 ATH_MSG_INFO(
"REGTEST TrigTau->eta() returns " << (*TrigTauItr)->eta());
2046 ATH_MSG_INFO(
"REGTEST TrigTau->phi() returns " << (*TrigTauItr)->phi());
2048 ATH_MSG_INFO(
"REGTEST TrigTau->nMatchedTracks() returns " << (*TrigTauItr)->nMatchedTracks());
2049 ATH_MSG_INFO(
"REGTEST TrigTau->simpleEtFlow() returns " << (*TrigTauItr)->simpleEtFlow());
2050 ATH_MSG_INFO(
"REGTEST TrigTau->m() returns " << (*TrigTauItr)->m());
2054 ATH_MSG_INFO(
"REGTEST ==========END of TrigTau DUMP===========");
2057 return StatusCode::SUCCESS;
2064 ATH_MSG_INFO(
"REGTEST ==========START of TrigTauTracksInfo DUMP===========");
2069 StatusCode sc = evtStore()->retrieve(trigTau,lastTrigTau);
2070 if (
sc.isFailure()) {
2072 return StatusCode::SUCCESS;
2076 for (
int i=0; trigTau != lastTrigTau; ++trigTau, ++
i) {
2082 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->pt() returns " << thisTrigTau->
pt());
2083 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->eta() returns " << thisTrigTau->
eta());
2084 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->phi() returns " << thisTrigTau->
phi());
2090 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->charge() returns " << thisTrigTau->
charge());
2091 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->leadingTrackPt() returns " <<
2093 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->scalarPtSumCore() returns " <<
2095 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->scalarPtSumIso() returns " <<
2097 ATH_MSG_INFO(
"REGTEST TrigTauTracksInfo->threeFastestTracks().pt() returns " <<
2101 ATH_MSG_INFO(
"REGTEST ==========END of TrigTauTracksInfo DUMP===========");
2104 return StatusCode::SUCCESS;
2113 ATH_MSG_INFO(
"REGTEST ==========START of HLTResult DUMP===========");
2117 if(
sc.isFailure() || !hltResult ) {
2119 return StatusCode::SUCCESS;
2121 ATH_MSG_INFO(
"REGTEST HLTResult_L2 Successfully Retrieved");
2124 ATH_MSG_INFO(
"REGTEST HLTResult_L2->size() returns " << hltResult->
size());
2128 ATH_MSG_INFO(
"REGTEST ==========END of HLTResult DUMP===========");
2132 return StatusCode::SUCCESS;
2142 std::string TrigInDetTrackTags[]={
"HLT_TrigIDSCAN_Jet",
2143 "HLT_TrigIDSCAN_Tau",
2144 "HLT_TrigIDSCAN_eGamma",
2145 "HLT_TrigIDSCAN_Muon",
2146 "HLT_TrigIDSCAN_muonIso",
2147 "HLT_TrigIDSCAN_Bphysics",
2148 "HLT_TrigIDSCAN_FullScan",
2149 "HLT_TrigIDSCAN_Cosmics",
2150 "HLT_TrigIDSCAN_eGamma_Brem",
2151 "HLT_TrigIDSCAN_Tile",
2152 "HLT_TrigSiTrack_eGamma",
2153 "HLT_TrigSiTrack_Muon",
2154 "HLT_TrigSiTrack_muonIso",
2155 "HLT_TrigSiTrack_Tau",
2156 "HLT_TrigSiTrack_Jet",
2157 "HLT_TrigSiTrack_Bphysics",
2158 "HLT_TrigSiTrack_FullScan",
2159 "HLT_TrigSiTrack_Cosmics",
2160 "HLT_TRTSegmentFinder",
2167 bool gotTruthMap =
false;
2169 if (
evtStore()->contains<TrigInDetTrackTruthMap>(
"TrigInDetTrackTruthMap")) {
2173 ATH_MSG_WARNING(
" could not retrieve TrackTruthMap with key TrigInDetTruthMap");
2179 else ATH_MSG_DEBUG(
" didn't find any TrackTruthMap objects with key TrigInDetTruthMap");
2181 ATH_MSG_INFO(
"REGTEST ==========START of TrigInDetTrackCollection DUMP===========");
2183 for (
int iTag=0; iTag < ntag; iTag++) {
2186 if (
sc.isFailure()) {
2187 ATH_MSG_DEBUG(
"REGTEST No TrigInDetTrackCollection found with key " << TrigInDetTrackTags[iTag]);
2190 ATH_MSG_INFO(
"REGTEST TrigInDetTrackCollections retrieved with key "
2191 << TrigInDetTrackTags[iTag]);
2193 ATH_MSG_INFO(
"REGTEST TrigInDetTrackCollection->size() returns " << trigInDetTrackCollection->
size());
2198 for (
int j=0; trkItr != trkItrE; ++trkItr, ++j ) {
2200 ATH_MSG_INFO(
"REGTEST Looking at TrigInDetTrack " << j);
2202 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->algorithmId() returns " << (*trkItr)->algorithmId());
2203 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->chi2() returns " << (*trkItr)->chi2());
2204 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NStrawHits() returns " << (*trkItr)->NStrawHits());
2206 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NStraw() returns " << (*trkItr)->NStraw());
2207 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NStrawTime() returns " << (*trkItr)->NStrawTime());
2208 ATH_MSG_INFO(
"REGTEST TrigInDetTrack->NTRHits() returns " << (*trkItr)->NTRHits());
2210 ATH_MSG_INFO(
"REGTEST Looking at TrigInDetFitTrack->param()");
2213 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->a0() returns " << my_param->
a0());
2214 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->z0() returns " << my_param->
z0());
2215 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->phi0() returns " << my_param->
phi0());
2216 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->eta() returns " << my_param->
eta());
2217 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->pT() returns " << my_param->
pT());
2218 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->ea0() returns " << my_param->
ea0());
2219 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->ez0() returns " << my_param->
ez0());
2220 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->ephi0() returns " << my_param->
ephi0());
2221 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->eeta() returns " << my_param->
eeta());
2222 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->epT() returns " << my_param->
epT());
2225 ATH_MSG_INFO(
"REGTEST Looking at covarience matrix: ");
2227 const std::vector<double>*
const my_cov = my_param->
cov();
2229 ATH_MSG_INFO(
"REGTEST covarience matrix NULL (not stored)");
2231 std::vector<double> tempcov = *my_cov;
2234 for(iter = tempcov.begin(),
k=0; iter != tempcov.end(); ++iter, ++
k){
2235 ATH_MSG_INFO(
"REGTEST TrigInDetTrackFitPar->cov() element " <<
k <<
" is " << (*iter));
2239 if (!(pTruthMap->
hasTruth((*trkItr)))) {
2255 if (pSiTruth->
cptr() != 0) {
2264 if (pTRTtruth->
cptr() != 0) {
2274 ATH_MSG_INFO(
"REGTEST ==========END of TrigInDetTrackCollection DUMP===========");
2277 return StatusCode::SUCCESS;
2286 ATH_MSG_INFO(
"REGTEST ==========START of TrigVertexCollection DUMP===========");
2291 StatusCode sc = evtStore()->retrieve(trigVertex,lastTrigVertex);
2292 if (
sc.isFailure()) {
2294 return StatusCode::SUCCESS;
2296 ATH_MSG_INFO(
"REGTEST TrigVertexCollection retrieved");
2299 for (
int i=0; trigVertex != lastTrigVertex; ++trigVertex, ++
i) {
2301 ATH_MSG_INFO(
"REGTEST Looking at TrigVertexCollection " <<
i);
2306 for (
int j=0; VertexItr != VertexItrE; ++VertexItr, ++j ) {
2310 ATH_MSG_INFO(
"REGTEST TrigVertex->x() returns " << (*VertexItr)->x());
2311 ATH_MSG_INFO(
"REGTEST TrigVertex->y() returns " << (*VertexItr)->y());
2312 ATH_MSG_INFO(
"REGTEST TrigVertex->z() returns " << (*VertexItr)->z());
2313 ATH_MSG_INFO(
"REGTEST TrigVertex->energyFraction() returns " << (*VertexItr)->energyFraction());
2314 ATH_MSG_INFO(
"REGTEST TrigVertex->ndof() returns " << (*VertexItr)->ndof());
2316 if(((*VertexItr)->tracks())!=NULL ){
2318 TrackInVertexList::const_iterator trkItr = (*VertexItr)->tracks()->begin();
2319 TrackInVertexList::const_iterator trkItrE = (*VertexItr)->tracks()->end();
2321 for (
int p=0; trkItr != trkItrE; ++trkItr, ++
p ) {
2326 ATH_MSG_INFO(
"REGTEST tracks->NStrawHits() returns " << (*trkItr)->NStrawHits());
2327 ATH_MSG_INFO(
"REGTEST tracks->NStrawTime() returns " << (*trkItr)->NStrawTime());
2328 ATH_MSG_INFO(
"REGTEST tracks->NTRHits() returns " << (*trkItr)->NTRHits());
2334 ATH_MSG_INFO(
"REGTEST ==========END of TrigVertexCollection DUMP===========");
2337 return StatusCode::SUCCESS;
2348 ATH_MSG_INFO(
"REGTEST ==========START of TrigEFBphysContainer DUMP===========");
2350 std::string EFBphysTags[]={
"HLT_xAOD__TrigBphysContainer_EFBMuMuFex",
2351 "HLT_xAOD__TrigBphysContainer_EFBMuMuXFex",
2352 "HLT_xAOD__TrigBphysContainer_EFDsPhiPiFex",
2353 "HLT_xAOD__TrigBphysContainer_EFMuPairs",
2354 "HLT_xAOD__TrigBphysContainer_EFMultiMuFex",
2355 "HLT_xAOD__TrigBphysContainer_EFTrackMass"
2358 int ntag= (
int)
sizeof(EFBphysTags) /
sizeof(EFBphysTags[0]);
2361 for (
int itag=0; itag<ntag; itag++){
2364 if (
sc.isFailure()) {
2365 ATH_MSG_INFO(
"REGTEST No TrigEFBphysContainer found with tag " << EFBphysTags[itag]);
2369 ATH_MSG_INFO(
"REGTEST TrigEFBphysContainer found with tag " << EFBphysTags[itag]
2370 <<
" and size " << trigEFBphys->
size());
2379 for (
int j=0; EFBphysItr != EFBphysItrE; ++EFBphysItr, ++j ) {
2383 ATH_MSG_INFO(
"REGTEST TrigEFBphys->eta() returns " << (*EFBphysItr)->eta());
2384 ATH_MSG_INFO(
"REGTEST TrigEFBphys->phi() returns " << (*EFBphysItr)->phi());
2385 ATH_MSG_INFO(
"REGTEST TrigEFBphys->mass() returns " << (*EFBphysItr)->mass());
2386 ATH_MSG_INFO(
"REGTEST TrigEFBphys->fitmass() returns " << (*EFBphysItr)->fitmass());
2388 ATH_MSG_INFO(
"REGTEST TrigEFBphys->roiId() returns " << (*EFBphysItr)->roiId());
2389 ATH_MSG_INFO(
"REGTEST TrigEFBphys->particleType() returns " << (*EFBphysItr)->particleType());
2391 if( (*EFBphysItr)->secondaryDecay() != NULL){
2394 ATH_MSG_INFO(
"REGTEST pSecondDecay->eta() returns " << psecond->
eta());
2395 ATH_MSG_INFO(
"REGTEST pSecondDecay->phi() returns " << psecond->
phi());
2406 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*EFBphysItr)->trackParticleLinks();
2407 if (trackVector.size() != 0) {
2408 ATH_MSG_INFO(
" REGTEST got track vector size: " << trackVector.size());
2412 std::vector<ElementLink<xAOD::TrackParticleContainer> >::const_iterator trkIt=trackVector.begin();
2413 for (
int itrk=0 ; trkIt!= trackVector.end(); ++itrk, ++trkIt) {
2414 if (!(trkIt->isValid())) {
2415 ATH_MSG_WARNING(
"TrackParticleContainer::Invalid ElementLink to track ");
2422 double phi = trackPerigee->parameters()[
Trk::phi];
2424 double px = trackPerigee->momentum()[
Trk::px];
2425 double py = trackPerigee->momentum()[
Trk::py];
2435 ATH_MSG_INFO(
"REGTEST ==========END of TrigEFBphysContainer DUMP===========");
2438 return StatusCode::SUCCESS;
2445 ATH_MSG_INFO(
"REGTEST ==========START of TrigL2BphysContainer DUMP===========");
2447 std::string L2BphysTags[]={
"HLT_xAOD__TrigBphysContainer_L2BMuMuFex",
2448 "HLT_xAOD__TrigBphysContainer_L2BMuMuXFex",
2449 "HLT_xAOD__TrigBphysContainer_L2DsPhiPiFexDs",
2450 "HLT_xAOD__TrigBphysContainer_L2DsPhiPiFexPhi",
2451 "HLT_xAOD__TrigBphysContainer_L2JpsieeFex",
2452 "HLT_xAOD__TrigBphysContainer_L2MultiMuFex",
2453 "HLT_xAOD__TrigBphysContainer_L2TrackMass",
2455 const int ntag = (
int)
sizeof(L2BphysTags) /
sizeof(L2BphysTags[0]);
2458 for (
int itag=0; itag<ntag; itag++){
2461 if (
sc.isFailure()) {
2462 ATH_MSG_INFO(
"REGTEST No TrigL2BphysContainer found with tag " << L2BphysTags[itag]);
2466 ATH_MSG_INFO(
"REGTEST TrigL2BphysContainer found with tag " << L2BphysTags[itag]
2467 <<
" and size " << trigL2Bphys->
size());
2476 for (
int j=0; L2BphysItr != L2BphysItrE; ++L2BphysItr, ++j ) {
2480 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->eta() returns " << (*L2BphysItr)->eta());
2481 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->phi() returns " << (*L2BphysItr)->phi());
2482 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->mass() returns " << (*L2BphysItr)->mass());
2483 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->fitmass() returns " << (*L2BphysItr)->fitmass());
2485 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->roiId() returns " << (*L2BphysItr)->roiId());
2486 ATH_MSG_INFO(
"REGTEST TrigL2Bphys->particleType() returns " << (*L2BphysItr)->particleType());
2488 if( (*L2BphysItr)->secondaryDecay() != NULL){
2491 ATH_MSG_INFO(
"REGTEST pSecondDecay->eta() returns " << psecond->
eta());
2492 ATH_MSG_INFO(
"REGTEST pSecondDecay->phi() returns " << psecond->
phi());
2500 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*L2BphysItr)->trackParticleLinks();
2501 if (trackVector.size() != 0) {
2502 ATH_MSG_INFO(
" REGTEST got track vector size: " << trackVector.size());
2506 std::vector<ElementLink<xAOD::TrackParticleContainer> >::const_iterator trkIt=trackVector.begin();
2507 for (
int itrk=0 ; trkIt!= trackVector.end(); ++itrk, ++trkIt) {
2508 if (!(trkIt->isValid())) {
2509 ATH_MSG_WARNING(
"TrackParticleContainer::Invalid ElementLink to track ");
2516 double phi = trackPerigee->parameters()[
Trk::phi];
2518 double px = trackPerigee->momentum()[
Trk::px];
2519 double py = trackPerigee->momentum()[
Trk::py];
2529 ATH_MSG_INFO(
"REGTEST ==========END of TrigL2BphysContainer DUMP===========");
2532 return StatusCode::SUCCESS;
2543 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::JetContainer DUMP===========");
2545 std::string containerName[30] = {
"HLT_xAOD__JetContainer_a4tcemjesFS",
2546 "HLT_xAOD__JetContainer_a4tcemsubFS",
2547 "HLT_xAOD__JetContainer_a4tcemsubjesFS",
2548 "HLT_xAOD__JetContainer_a4tcemnojcalibFS",
2549 "HLT_xAOD__JetContainer_a4tcemjesPS",
2550 "HLT_xAOD__JetContainer_a4tcemnojcalibPS",
2551 "HLT_xAOD__JetContainer_a4tclcwjesFS",
2552 "HLT_xAOD__JetContainer_a4tclcwsubFS",
2553 "HLT_xAOD__JetContainer_a4tclcwsubjesFS",
2554 "HLT_xAOD__JetContainer_a4tclcwnojcalibFS",
2555 "HLT_xAOD__JetContainer_a4tclcwjesPS",
2556 "HLT_xAOD__JetContainer_a4tclcwnojcalibPS",
2557 "HLT_xAOD__JetContainer_a4TTemnojcalibFS",
2558 "HLT_xAOD__JetContainer_a4TThadnojcalibFS",
2559 "HLT_xAOD__JetContainer_a10tcemjesFS",
2560 "HLT_xAOD__JetContainer_a10tcemsubFS",
2561 "HLT_xAOD__JetContainer_a10tcemsubjesFS",
2562 "HLT_xAOD__JetContainer_a10tcemnojcalibFS",
2563 "HLT_xAOD__JetContainer_a10tcemjesPS",
2564 "HLT_xAOD__JetContainer_a10tcemnojcalibPS",
2565 "HLT_xAOD__JetContainer_a10tclcwjesFS",
2566 "HLT_xAOD__JetContainer_a10tclcwsubFS",
2567 "HLT_xAOD__JetContainer_a10tclcwsubjesFS",
2568 "HLT_xAOD__JetContainer_a10tclcwnojcalibFS",
2569 "HLT_xAOD__JetContainer_a10tclcwjesPS",
2570 "HLT_xAOD__JetContainer_a10tclcwnojcalibPS",
2571 "HLT_xAOD__JetContainer_a10TTemnojcalibFS",
2572 "HLT_xAOD__JetContainer_a10TThadnojcalibFS",
2573 "HLT_xAOD__JetContainer_a10r_tcemsubjesFS",
2574 "HLT_xAOD__JetContainer_TrigHLTJetDSSelectorCollection"};
2575 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,
2576 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};
2577 int containerInputType[30] = {1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 2,
2578 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 2, 1, 1};
2579 int containerSignalState[30] = {0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1,
2580 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0};
2581 bool onefilled =
false;
2583 for (
int icont=0; icont<30; ++icont) {
2588 if (
sc.isFailure()) {
2593 int jetContsize = jetCont->
size();
2594 ATH_MSG_INFO(
"REGTEST Got jet container " << containerName[icont] <<
", size: " << jetContsize);
2595 if (jetContsize != 0) {
2598 for(
const auto thisjet : *jetCont) {
2609 ATH_MSG_INFO(
"REGTEST rapidity: " << thisjet->rapidity() );
2614 ATH_MSG_INFO(
"REGTEST algorithm (kt: 0, cam: 1, antikt: 2, ...): " << thisjet->getAlgorithmType() <<
"; should be 2");
2615 if(thisjet->getAlgorithmType() != 2)
ATH_MSG_WARNING(
"Jet algorithm different from container");
2616 ATH_MSG_INFO(
"REGTEST size parameter: " << thisjet->getSizeParameter() <<
"; should be " << containerSizeParameter[icont]);
2617 if(thisjet->getSizeParameter() != containerSizeParameter[icont])
ATH_MSG_WARNING(
"Jet size different from container");
2618 ATH_MSG_INFO(
"REGTEST input (LCTopo: 0, EMTopo: 1, TopoTower: 2, ...): " << thisjet->getInputType() <<
"; should be " << containerInputType[icont]);
2619 if(thisjet->getInputType() != containerInputType[icont])
ATH_MSG_WARNING(
"Jet input different from container");
2620 ATH_MSG_INFO(
"REGTEST constituents signal state (uncalibrated: 0, calibrated: 1): " << thisjet->getConstituentsSignalState() <<
"; should be " << containerSignalState[icont]);
2621 if(thisjet->getConstituentsSignalState() != containerSignalState[icont])
ATH_MSG_WARNING(
"Jet constituents' signal state different from container");
2622 ATH_MSG_INFO(
"REGTEST number of constituents: " << thisjet->numConstituents() );
2629 unsigned int constitContsize = constitCont.
size();
2630 ATH_MSG_INFO(
"REGTEST Got constituent vector, size: " << constitContsize <<
"; should be " << thisjet->numConstituents());
2631 if(constitContsize != thisjet->numConstituents())
ATH_MSG_WARNING(
"Constituents container size different from number of constituents");
2633 if (constitContsize != 0) {
2731 std::vector<int> vecvalueint;
2733 int vecsize = vecvalueint.size();
2734 ATH_MSG_INFO(
"REGTEST Got NumTrkPt1000 vector, size: " << vecsize);
2737 for(
const auto & thisvalue : vecvalueint){
2739 ATH_MSG_INFO(
"REGTEST NumTrkPt1000 #" << j <<
": " << thisvalue);
2741 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2742 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2746 int vecsize = vecvalueint.size();
2747 ATH_MSG_INFO(
"REGTEST Got Got NumTrkPt500 vector, size: " << vecsize);
2750 for(
const auto & thisvalue : vecvalueint){
2752 ATH_MSG_INFO(
"REGTEST NumTrkPt500 #" << j <<
": " << thisvalue);
2754 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2755 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2761 std::vector<float> vecvalue;
2763 int vecsize = vecvalue.size();
2764 ATH_MSG_INFO(
"REGTEST Got JVF vector, size: " << vecsize);
2767 for(
const auto & thisvalue : vecvalue){
2769 ATH_MSG_INFO(
"REGTEST JVF #" << j <<
": " << thisvalue);
2771 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2772 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2776 int vecsize = vecvalue.size();
2777 ATH_MSG_INFO(
"REGTEST Got SumPtTrkPt1000 vector, size: " << vecsize);
2780 for(
const auto & thisvalue : vecvalue){
2782 ATH_MSG_INFO(
"REGTEST SumPtTrkPt1000 #" << j <<
": " << thisvalue);
2784 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2785 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2789 int vecsize = vecvalue.size();
2790 ATH_MSG_INFO(
"REGTEST Got SumPtTrkPt500 vector, size: " << vecsize);
2793 for(
const auto & thisvalue : vecvalue){
2795 ATH_MSG_INFO(
"REGTEST SumPtTrkPt500 #" << j <<
": " << thisvalue);
2797 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2798 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2802 int vecsize = vecvalue.size();
2803 ATH_MSG_INFO(
"REGTEST Got TrackWidthPt1000 vector, size: " << vecsize);
2806 for(
const auto & thisvalue : vecvalue){
2808 ATH_MSG_INFO(
"REGTEST TrackWidthPt1000 #" << j <<
": " << thisvalue);
2810 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2811 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2815 int vecsize = vecvalue.size();
2816 ATH_MSG_INFO(
"REGTEST Got TrackWidthPt500 vector, size: " << vecsize);
2819 for(
const auto & thisvalue : vecvalue){
2821 ATH_MSG_INFO(
"REGTEST TrackWidthPt500 #" << j <<
": " << thisvalue);
2823 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2824 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2828 int vecsize = vecvalue.size();
2829 ATH_MSG_INFO(
"REGTEST Got EnergyPerSampling vector, size: " << vecsize);
2832 for(
const auto & thisvalue : vecvalue){
2834 ATH_MSG_INFO(
"REGTEST EnergyPerSampling #" << j <<
": " << thisvalue);
2836 ATH_MSG_INFO(
"REGTEST size of attribute vector == number of displayed attributes: " << (vecsize == j) );
2837 if (vecsize != j)
ATH_MSG_WARNING(
"REGTEST Problem with displaying this attribute");
2866 return StatusCode::SUCCESS;
2870 if (jetContsize ==
i)
ATH_MSG_INFO(
"REGTEST size of jet container == number of displayed jets: " << (jetContsize ==
i) );
2875 if (!onefilled)
ATH_MSG_DEBUG(
"There was no filled jet containers");
2877 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::JetContainer DUMP===========");
2881 return StatusCode::SUCCESS;
2890 ATH_MSG_INFO(
"REGTEST ==========START of TrigEFBjetContainer DUMP===========");
2895 StatusCode sc = evtStore()->retrieve(trigEFBjet,lastTrigEFBjet);
2896 if (
sc.isFailure()) {
2898 return StatusCode::SUCCESS;
2900 ATH_MSG_INFO(
"REGTEST TrigEFBjetContainers retrieved");
2903 for (
int i=0; trigEFBjet != lastTrigEFBjet; ++trigEFBjet, ++
i) {
2905 ATH_MSG_INFO(
"REGTEST Looking at TrigEFBjetContainer " <<
i);
2910 for (
int j=0; EFBjetItr != EFBjetItrE; ++EFBjetItr, ++j ) {
2914 ATH_MSG_INFO(
"REGTEST TrigEFBjet->prmVtx() returns " << (*EFBjetItr)->prmVtx());
2915 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xComb() returns " << (*EFBjetItr)->xComb());
2917 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xIP1D() returns " << (*EFBjetItr)->xIP1D());
2918 ATH_MSG_INFO(
"REGTEST TrigEFBjet->isValid() returns " << (*EFBjetItr)->isValid());
2919 ATH_MSG_INFO(
"REGTEST TrigEFBjet->roiId() returns " << (*EFBjetItr)->roiId());
2921 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xIP2D() returns " << (*EFBjetItr)->xIP2D());
2923 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xIP3D() returns " << (*EFBjetItr)->xIP3D());
2925 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xCHI2() returns " << (*EFBjetItr)->xCHI2());
2927 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xSV() returns " << (*EFBjetItr)->xSV());
2929 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xMVtx() returns " << (*EFBjetItr)->xMVtx());
2930 ATH_MSG_INFO(
"REGTEST TrigEFBjet->xEVtx() returns " << (*EFBjetItr)->xEVtx());
2931 ATH_MSG_INFO(
"REGTEST TrigEFBjet-> xNVtx() returns " << (*EFBjetItr)-> xNVtx());
2934 ATH_MSG_INFO(
"REGTEST ==========END of TrigEFBjetContainer DUMP===========");
2937 return StatusCode::SUCCESS;
2946 ATH_MSG_INFO(
"REGTEST ==========START of TrigL2BjetContainer DUMP===========");
2951 StatusCode sc = evtStore()->retrieve(trigL2Bjet,lastTrigL2Bjet);
2952 if (
sc.isFailure()) {
2954 return StatusCode::SUCCESS;
2956 ATH_MSG_INFO(
"REGTEST TrigL2BjetContainers retrieved");
2959 for (
int i=0; trigL2Bjet != lastTrigL2Bjet; ++trigL2Bjet, ++
i) {
2961 ATH_MSG_INFO(
"REGTEST Looking at TrigL2BjetContainer " <<
i);
2966 for (
int j=0; L2BjetItr != L2BjetItrE; ++L2BjetItr, ++j ) {
2969 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->prmVtx() returns " << (*L2BjetItr)->prmVtx());
2970 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xComb() returns " << (*L2BjetItr)->xComb());
2972 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xIP1D() returns " << (*L2BjetItr)->xIP1D());
2973 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->isValid() returns " << (*L2BjetItr)->isValid());
2974 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->roiId() returns " << (*L2BjetItr)->roiId());
2976 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xIP2D() returns " << (*L2BjetItr)->xIP2D());
2978 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xIP3D() returns " << (*L2BjetItr)->xIP3D());
2980 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xCHI2() returns " << (*L2BjetItr)->xCHI2());
2982 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xSV() returns " << (*L2BjetItr)->xSV());
2983 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xMVtx() returns " << (*L2BjetItr)->xMVtx());
2984 ATH_MSG_INFO(
"REGTEST TrigL2Bjet->xEVtx() returns " << (*L2BjetItr)->xEVtx());
2985 ATH_MSG_INFO(
"REGTEST TrigL2Bjet-> xNVtx() returns " << (*L2BjetItr)-> xNVtx());
2988 ATH_MSG_INFO(
"REGTEST ==========END of TrigL2BjetContainer DUMP===========");
2991 return StatusCode::SUCCESS;
3000 ATH_MSG_INFO(
"REGTEST ==========START of MuonFeature DUMP===========");
3005 StatusCode sc = evtStore()->retrieve(MuFeature,lastMuFeature);
3006 if (
sc.isFailure()) {
3008 return StatusCode::SUCCESS;
3013 for (
int i=0; MuFeature != lastMuFeature; ++MuFeature, ++
i) {
3018 ATH_MSG_INFO(
"REGTEST MuonFeature->roiId() returns " << thisMuFeature->
roiId());
3019 ATH_MSG_INFO(
"REGTEST MuonFeature->eta() returns " << thisMuFeature->
eta());
3020 ATH_MSG_INFO(
"REGTEST MuonFeature->phi() returns " << thisMuFeature->
phi());
3022 ATH_MSG_INFO(
"REGTEST MuonFeature->pt() returns " << thisMuFeature->
pt());
3025 ATH_MSG_INFO(
"REGTEST MuonFeature->zeta() returns " << thisMuFeature->
zeta());
3027 ATH_MSG_INFO(
"REGTEST MuonFeature->beta() returns " << thisMuFeature->
beta());
3030 ATH_MSG_INFO(
"REGTEST ==========END of MuonFeature DUMP===========");
3033 return StatusCode::SUCCESS;
3042 ATH_MSG_INFO(
"REGTEST ==========START of CombinedMuonFeature DUMP===========");
3047 StatusCode sc = evtStore()->retrieve(CombMuon,lastCombMuon);
3048 if (
sc.isFailure()) {
3050 return StatusCode::SUCCESS;
3052 ATH_MSG_INFO(
"REGTEST CombinedMuonFeatures retrieved");
3054 for (
int i=0; CombMuon != lastCombMuon; ++CombMuon, ++
i) {
3058 ATH_MSG_INFO(
"REGTEST Looking at CombinedMuonFeature " <<
i);
3059 ATH_MSG_INFO(
"REGTEST TrigPhoton->pt() returns " << (thisCombMuFeature)->
pt());
3060 ATH_MSG_INFO(
"REGTEST TrigPhoton->sigma_pt() returns " << (thisCombMuFeature)->sigma_pt());
3063 if ((thisCombMuFeature)->muFastTrackLink().
isValid() ) {
3064 ATH_MSG_INFO(
"REGTEST muFastTrack->eta() returns " << (thisCombMuFeature)->muFastTrack()->
eta());
3065 ATH_MSG_INFO(
"REGTEST muFastTrack->phi() returns " << (thisCombMuFeature)->muFastTrack()->
phi());
3066 ATH_MSG_INFO(
"REGTEST muFastTrack->radius() returns " << (thisCombMuFeature)->muFastTrack()->
radius());
3067 ATH_MSG_INFO(
"REGTEST muFastTrack->zeta() returns " << (thisCombMuFeature)->muFastTrack()->zeta());
3069 ATH_MSG_INFO(
"CombinedMuonFeature has no muFastTrack" );
3074 if ((thisCombMuFeature)->IDTrackLink().
isValid() ) {
3075 ATH_MSG_INFO(
"REGTEST IDTrack->algorithmId() returns " <<(thisCombMuFeature)->IDTrack()->algorithmId());
3076 ATH_MSG_INFO(
"REGTEST IDTrack->chi2() returns " << (thisCombMuFeature)->IDTrack()->
chi2());
3077 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawHits() returns " <<(thisCombMuFeature)->IDTrack()->NStrawHits());
3078 ATH_MSG_INFO(
"REGTEST IDTrack->NStraw() returns " << (thisCombMuFeature)->IDTrack()->NStraw());
3079 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawTime() returns " <<(thisCombMuFeature)->IDTrack()->NStrawTime());
3080 ATH_MSG_INFO(
"REGTEST IDTrack->NTRHits() returns " <<(thisCombMuFeature)->IDTrack()->NTRHits());
3087 ATH_MSG_INFO(
"REGTEST ==========END of CombinedMuonFeature DUMP===========");
3090 return StatusCode::SUCCESS;
3096 ATH_MSG_INFO(
"REGTEST ==========START of CombinedMuonFeatureContainer DUMP===========");
3101 StatusCode sc = evtStore()->retrieve(CombMuon,lastCombMuon);
3102 if (
sc.isFailure()) {
3103 ATH_MSG_INFO(
"REGTEST No CombinedMuonFeatureContainer found");
3104 return StatusCode::SUCCESS;
3106 ATH_MSG_INFO(
"REGTEST CombinedMuonFeaturesContainer retrieved");
3108 for (
int j=0; CombMuon != lastCombMuon; ++CombMuon, ++j) {
3109 ATH_MSG_INFO(
"REGTEST Looking at CombinedMuonFeatureContainer " << j);
3117 ATH_MSG_INFO(
"REGTEST Looking at CombinedMuonFeature " <<
i);
3119 ATH_MSG_INFO(
"REGTEST CombinedMuonFeature->pt() returns " << (thisCombMuFeature)->
pt());
3120 ATH_MSG_INFO(
"REGTEST CombinedMuonFeature->sigma_pt() returns " << (thisCombMuFeature)->sigma_pt());
3124 ATH_MSG_INFO(
"REGTEST muFastTrack->eta() returns " << (thisCombMuFeature)->muFastTrack()->
eta());
3125 ATH_MSG_INFO(
"REGTEST muFastTrack->phi() returns " << (thisCombMuFeature)->muFastTrack()->
phi());
3126 ATH_MSG_INFO(
"REGTEST muFastTrack->radius() returns " << (thisCombMuFeature)->muFastTrack()->
radius());
3127 ATH_MSG_INFO(
"REGTEST muFastTrack->zeta() returns " << (thisCombMuFeature)->muFastTrack()->zeta());
3134 ATH_MSG_INFO(
"REGTEST IDTrack->algorithmId() returns " <<(thisCombMuFeature)->IDTrack()->algorithmId());
3135 ATH_MSG_INFO(
"REGTEST IDTrack->chi2() returns " << (thisCombMuFeature)->IDTrack()->
chi2());
3136 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawHits() returns " <<(thisCombMuFeature)->IDTrack()->NStrawHits());
3138 ATH_MSG_INFO(
"REGTEST IDTrack->NStraw() returns " << (thisCombMuFeature)->IDTrack()->NStraw());
3139 ATH_MSG_INFO(
"REGTEST IDTrack->NStrawTime() returns " <<(thisCombMuFeature)->IDTrack()->NStrawTime());
3140 ATH_MSG_INFO(
"REGTEST IDTrack->NTRHits() returns " <<(thisCombMuFeature)->IDTrack()->NTRHits());
3147 return StatusCode::SUCCESS;
3155 ATH_MSG_INFO(
"REGTEST ==========START of TrigEMCluster DUMP===========");
3160 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3161 if (
sc.isFailure()) {
3163 return StatusCode::SUCCESS;
3167 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3174 ATH_MSG_INFO(
"REGTEST TrigEMCluster->e() returns " << thisEMCluster->
e());
3175 ATH_MSG_INFO(
"REGTEST TrigEMCluster->phi() returns " << thisEMCluster->
phi());
3177 ATH_MSG_INFO(
"REGTEST TrigEMCluster->eta() returns " << thisEMCluster->
eta());
3178 ATH_MSG_INFO(
"REGTEST TrigEMCluster->e237() returns " << thisEMCluster->
e237());
3179 ATH_MSG_INFO(
"REGTEST TrigEMCluster->e277() returns " << thisEMCluster->
e277());
3181 ATH_MSG_INFO(
"REGTEST TrigEMCluster->weta2() returns " << thisEMCluster->
weta2());
3182 ATH_MSG_INFO(
"REGTEST TrigEMCluster->ehad1() returns " << thisEMCluster->
ehad1());
3183 ATH_MSG_INFO(
"REGTEST TrigEMCluster->Eta1() returns " << thisEMCluster->
Eta1());
3186 ATH_MSG_INFO(
"REGTEST ==========END of TrigEMCluster DUMP===========");
3189 return StatusCode::SUCCESS;
3196 ATH_MSG_INFO(
"REGTEST ==========START of TrigEMCluster DUMP===========");
3201 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3202 if (
sc.isFailure()) {
3204 return StatusCode::SUCCESS;
3209 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3213 ATH_MSG_INFO(
"REGTEST Looking at xAOD::TrigEMCluster " <<
i);
3215 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->energy() returns " << thisEMCluster->
energy());
3216 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->e() returns " << thisEMCluster->
energy());
3217 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->phi() returns " << thisEMCluster->
phi());
3219 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->eta() returns " << thisEMCluster->
eta());
3220 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->e237() returns " << thisEMCluster->
e237());
3221 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->e277() returns " << thisEMCluster->
e277());
3222 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->fracs1() returns " << thisEMCluster->
fracs1());
3223 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->weta2() returns " << thisEMCluster->
weta2());
3224 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->ehad1() returns " << thisEMCluster->
ehad1());
3225 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->Eta1() returns " << thisEMCluster->
eta1());
3228 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigEMCluster DUMP===========");
3231 return StatusCode::SUCCESS;
3240 ATH_MSG_INFO(
"REGTEST ==========START of TrigTauClusterContainer/TrigTauClusterDetailsContainer DUMP===========");
3245 StatusCode sc = evtStore()->retrieve(TauCluster,lastTauCluster);
3246 if (
sc.isFailure()) {
3247 ATH_MSG_INFO(
"REGTEST No TrigTauClusterContainer found");
3248 return StatusCode::SUCCESS;
3250 ATH_MSG_INFO(
"REGTEST TrigTauClusterContainer retrieved");
3254 for (
int i=0; TauCluster != lastTauCluster; ++TauCluster, ++
i) {
3256 ATH_MSG_INFO(
"REGTEST Looking at TrigTauClusterContainer " <<
i);
3261 for (
int j=0; TauClusterItr != TauClusterItrE; ++TauClusterItr, ++j ) {
3263 ATH_MSG_INFO(
"REGTEST Looking at TrigTauCluster " << j);
3265 ATH_MSG_INFO(
"REGTEST TrigTauCluster->energy() returns " << (*TauClusterItr)->energy());
3266 ATH_MSG_INFO(
"REGTEST TrigTauCluster->et() returns " << (*TauClusterItr)->et());
3267 ATH_MSG_INFO(
"REGTEST TrigTauCluster->EMCalibEnergy() returns " << (*TauClusterItr)->EMCalibEnergy());
3269 ATH_MSG_INFO(
"REGTEST TrigTauCluster->EMenergy() returns " << (*TauClusterItr)->EMenergy());
3270 ATH_MSG_INFO(
"REGTEST TrigTauCluster->HADenergy() returns " << (*TauClusterItr)->HADenergy());
3271 ATH_MSG_INFO(
"REGTEST TrigTauCluster->eta() returns " << (*TauClusterItr)->eta());
3272 ATH_MSG_INFO(
"REGTEST TrigTauCluster->phi() returns " << (*TauClusterItr)->phi());
3273 ATH_MSG_INFO(
"REGTEST TrigTauCluster->IsoFrac() returns " << (*TauClusterItr)->IsoFrac());
3274 ATH_MSG_INFO(
"REGTEST TrigTauCluster->stripWidth() returns " << (*TauClusterItr)->stripWidth());
3275 ATH_MSG_INFO(
"REGTEST TrigTauCluster->numStripCells() returns " << (*TauClusterItr)->numStripCells());
3276 ATH_MSG_INFO(
"REGTEST TrigTauCluster->CaloRadius() returns " << (*TauClusterItr)->CaloRadius());
3277 ATH_MSG_INFO(
"REGTEST TrigTauCluster->numTotCells() returns " << (*TauClusterItr)->numTotCells());
3280 if( (*TauClusterItr)->clusterDetails() == 0 )
3283 ATH_MSG_INFO(
"REGTEST TrigTauCluster->Energy in Narrow cone : EM[0/1/2/3]="
3284 <<(*TauClusterItr)->EMenergyNarrow(0) <<
" " << (*TauClusterItr)->EMenergyNarrow(1) <<
" "
3285 <<(*TauClusterItr)->EMenergyNarrow(2) <<
" ");
3286 ATH_MSG_INFO(
"REGTEST TrigTauCluster->Energy in Narrow cone : HAD[0/1/2]="
3287 <<(*TauClusterItr)->HADenergyNarrow(0) <<
" " << (*TauClusterItr)->HADenergyNarrow(1) <<
" "
3288 <<(*TauClusterItr)->HADenergyNarrow(2));
3296 sc = evtStore()->retrieve(TauDetailsCluster,lastTauDetailsCluster);
3297 if (
sc.isFailure()) {
3298 ATH_MSG_INFO(
"REGTEST No TrigTauDetailsClusterContainer found");
3299 return StatusCode::SUCCESS;
3301 ATH_MSG_INFO(
"REGTEST TrigTauDetailsClusterContainer retrieved");
3303 for (
int i=0; TauDetailsCluster != lastTauDetailsCluster; ++TauDetailsCluster, ++
i) {
3305 ATH_MSG_INFO(
"REGTEST Looking at TrigTauClusterDetailsContainer " <<
i);
3310 for (
int j=0; TauDetailsClusterItr != TauDetailsClusterItrE; ++TauDetailsClusterItr, ++j ) {
3312 ATH_MSG_INFO(
"REGTEST Looking at TrigTauClusterDetails " << j);
3314 for(
unsigned int i = 0;
i<4; ++
i ) {
3315 ATH_MSG_INFO(
"REGTEST TrigTauClusterDetails-> EM["<<
i <<
"] Radius/EnergyNar/EnergyMid/EnergyWid returns "
3316 << (*TauDetailsClusterItr)->EMRadius(
i) <<
" "
3317 << (*TauDetailsClusterItr)->EMenergyNarrow(
i) <<
" "
3318 << (*TauDetailsClusterItr)->EMenergyMedium(
i) <<
" "
3319 << (*TauDetailsClusterItr)->EMenergyWide(
i));
3321 for(
unsigned int i = 0;
i<3; ++
i ) {
3322 ATH_MSG_INFO(
"REGTEST TrigTauClusterDetails-> HAD["<<
i <<
"] Radius/EnergyNar/EnergyMid/EnergyWid returns "
3323 << (*TauDetailsClusterItr)->HADRadius(
i) <<
" "
3324 << (*TauDetailsClusterItr)->HADenergyNarrow(
i) <<
" "
3325 << (*TauDetailsClusterItr)->HADenergyMedium(
i) <<
" "
3326 << (*TauDetailsClusterItr)->HADenergyWide(
i));
3331 if( nDetails != nClusters)
3332 ATH_MSG_WARNING(
"REGTEST inconsistent number of TrigTauClusters ("<< nClusters<<
") and TrigTauClusterDetails ("
3333 << nDetails <<
")");
3335 ATH_MSG_INFO(
"REGTEST ==========END of TrigTauClusterContainer/TrigTauClusterDetailsContainer DUMP===========");
3339 return StatusCode::SUCCESS;
3348 ATH_MSG_INFO(
"REGTEST ==========START of TrigEMClusterContainer DUMP===========");
3353 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3354 if (
sc.isFailure()) {
3355 ATH_MSG_INFO(
"REGTEST No TrigEMClusterContainer found");
3356 return StatusCode::SUCCESS;
3358 ATH_MSG_INFO(
"REGTEST TrigEMClusterContainer retrieved");
3361 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3363 ATH_MSG_INFO(
"REGTEST Looking at TrigEMClusterContainer " <<
i);
3368 for (
int j=0; EMClusterItr != EMClusterItrE; ++EMClusterItr, ++j ) {
3370 ATH_MSG_INFO(
"REGTEST Looking at TrigEMCluster " << j);
3372 ATH_MSG_INFO(
"REGTEST TrigEMCluster->energy() returns " << (*EMClusterItr)->energy());
3373 ATH_MSG_INFO(
"REGTEST TrigEMCluster->et() returns " << (*EMClusterItr)->et());
3374 ATH_MSG_INFO(
"REGTEST TrigEMCluster->eta() returns " << (*EMClusterItr)->eta());
3375 ATH_MSG_INFO(
"REGTEST TrigEMCluster->phi() returns " << (*EMClusterItr)->phi());
3380 (*EMClusterItr)->print(
msg());
3385 ATH_MSG_INFO(
"REGTEST ==========END of TrigEMClusterContainer DUMP===========");
3388 return StatusCode::SUCCESS;
3395 ATH_MSG_INFO(
"REGTEST ==========START of xAODTrigEMClusterContainer DUMP===========");
3400 StatusCode sc = evtStore()->retrieve(EMCluster,lastEMCluster);
3401 if (
sc.isFailure()) {
3402 ATH_MSG_INFO(
"REGTEST No xAOD::TrigEMClusterContainer found");
3403 return StatusCode::SUCCESS;
3405 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMClusterContainer retrieved");
3408 for (
int i=0; EMCluster != lastEMCluster; ++EMCluster, ++
i) {
3410 ATH_MSG_INFO(
"REGTEST Looking at xAOD::TrigEMClusterContainer " <<
i);
3415 for (
int j=0; EMClusterItr != EMClusterItrE; ++EMClusterItr, ++j ) {
3417 ATH_MSG_INFO(
"REGTEST Looking at xAOD::TrigEMCluster " << j);
3418 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->energy() returns " << (*EMClusterItr)->energy());
3419 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->et() returns " << (*EMClusterItr)->et());
3420 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->eta() returns " << (*EMClusterItr)->eta());
3421 ATH_MSG_INFO(
"REGTEST xAOD::TrigEMCluster->phi() returns " << (*EMClusterItr)->phi());
3431 ATH_MSG_INFO(
"REGTEST ==========END of TrigEMClusterContainer DUMP===========");
3434 return StatusCode::SUCCESS;
3442 ATH_MSG_INFO(
"REGTEST ==========START of TileMuFeatureContainer DUMP===========");
3448 if (
sc.isFailure()) {
3449 ATH_MSG_INFO(
"REGTEST No TileMuFeatureContainer found");
3450 return StatusCode::SUCCESS;
3452 ATH_MSG_INFO(
"REGTEST TileMuFeatureContainer retrieved");
3455 ATH_MSG_INFO(
"REGTEST Looking at TileMuFeatureContainer " <<
i);
3459 for(
int j=0; TileMuItr != TileMuItrE; ++TileMuItr, ++j) {
3460 ATH_MSG_INFO(
"REGTEST Looking at TileMuFeature " << j);
3461 ATH_MSG_INFO(
"REGTEST TileMuFeature->eta() returns " << (*TileMuItr)->eta());
3462 ATH_MSG_INFO(
"REGTEST TileMuFeature->phi() returns " << (*TileMuItr)->phi());
3463 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(0) returns " << (*TileMuItr)->enedep().at(0));
3464 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(1) returns " << (*TileMuItr)->enedep().at(1));
3465 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(2) returns " << (*TileMuItr)->enedep().at(2));
3466 ATH_MSG_INFO(
"REGTEST TileMuFeature->enedep().at(3) returns " << (*TileMuItr)->enedep().at(3));
3467 ATH_MSG_INFO(
"REGTEST TileMuFeature->quality() returns " << (*TileMuItr)->quality());
3472 ATH_MSG_INFO(
"REGTEST ==========END of TileMuFeatureContainer DUMP===========");
3475 return StatusCode::SUCCESS;
3483 ATH_MSG_INFO(
"REGTEST ==========START of TileTrackMuFeatureContainer DUMP===========");
3488 StatusCode sc = evtStore()->retrieve(TileTrackMu, lastTileTrackMu);
3489 if (
sc.isFailure()) {
3490 ATH_MSG_INFO(
"REGTEST No TileTrackMuFeatureContainer found");
3491 return StatusCode::SUCCESS;
3493 ATH_MSG_INFO(
"REGTEST TileTrackMuFeatureContainer retrieved");
3495 for (
int i=0; TileTrackMu!=lastTileTrackMu; ++TileTrackMu, ++
i) {
3496 ATH_MSG_INFO(
"REGTEST Looking at TileTrackMuFeatureContainer " <<
i);
3501 for (
int j=0; TileTrackMuItr != TileTrackMuItrE; ++TileTrackMuItr, ++j) {
3502 ATH_MSG_INFO(
"REGTEST Looking at TileTrackMuFeature " << j);
3503 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->PtTR_Trk() returns " << (*TileTrackMuItr)->PtTR_Trk());
3504 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->EtaTR_Trk() returns " << (*TileTrackMuItr)->EtaTR_Trk());
3505 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->PhiTR_Trk() returns " << (*TileTrackMuItr)->PhiTR_Trk());
3506 ATH_MSG_INFO(
"REGTEST TileTrackMuFeature->Typ_IDTrk() returns " << (*TileTrackMuItr)->Typ_IDTrk());
3512 ATH_MSG_INFO(
"REGTEST No TileMuFeature (Something Wrong)");
3539 ATH_MSG_INFO(
"REGTEST IDScanLink->algorithmId() returns " <<
Track->algorithmId());
3541 ATH_MSG_INFO(
"REGTEST IDScanLink->NStrawHits() returns ");
3543 ATH_MSG_INFO(
"REGTEST IDScanLink->NStrawTime() returns " <<
Track->NStrawTime());
3545 ATH_MSG_INFO(
"REGTEST IDScanLink->param()->phi0() returns " <<
Track->param()->phi0());
3546 ATH_MSG_INFO(
"REGTEST IDScanLink->param()->eta() returns " <<
Track->param()->eta());
3547 ATH_MSG_INFO(
"REGTEST IDScanLink->param()->pT() returns " <<
Track->param()->pT());
3552 ATH_MSG_INFO(
"REGTEST ==========END of TileTrackMuFeatureContainer DUMP===========");
3553 ATH_MSG_DEBUG(
"dumpTileTrackMuFeatureContainer() succeeded");
3554 return StatusCode::SUCCESS;
3561 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TauJetContainer DUMP===========" );
3563 StatusCode sc = evtStore() ->
retrieve (TauJetcont,
"HLT_xAOD__TauJetContainer_TrigTauRecMerged");
3565 if (
sc.isFailure()) {
3567 ATH_MSG_INFO(
"REGTEST No Tau container HLT_xAOD__TauJetContainer_TrigTauRecMerged");
3570 return StatusCode::SUCCESS;
3576 ATH_MSG_INFO(
"REGTEST (*tauIt)->eta() returns " << (*tauIt)->eta() );
3577 ATH_MSG_INFO(
"REGTEST (*tauIt)->phi() returns " << (*tauIt)->phi() );
3578 ATH_MSG_INFO(
"REGTEST (*tauIt)->pt() returns " << (*tauIt)->pt() );
3582 #ifndef XAODTAU_VERSIONS_TAUJET_V3_H
3583 EFnTracks = (*tauIt)->nTracks();
3588 ATH_MSG_INFO(
"REGTEST (*tauIt)->nTracks() returns " << EFnTracks );
3591 int EFWidenTrack = -1;
3592 #ifndef XAODTAU_VERSIONS_TAUJET_V3_H
3593 EFWidenTrack = (*tauIt)->nWideTracks();
3598 ATH_MSG_INFO(
"REGTEST (*tauIt)->nWideTracks() returns " << EFWidenTrack );
3621 ATH_MSG_INFO(
"REGTEST TauDetails->etOverPtLeadTrk() returns " << etOvPtLead);
3630 ATH_MSG_INFO(
"REGTEST TauDetails->isolFrac() returns " << IsoFrac);
3659 ATH_MSG_INFO(
"REGTEST TauDetails->etEMAtEMScale() returns " << EtEm);}
3663 ATH_MSG_INFO(
"REGTEST TauDetails->etHadAtEMScale() returns " << EtHad);}
3666 if( !(*tauIt)->jetLink().isValid() ) {
3668 return StatusCode::SUCCESS;
3671 const xAOD::Jet* pJetSeed = *((*tauIt)->jetLink());
3676 for (
int clusCount = 0; clusItr != clusItrE; ++clusItr, ++clusCount) {
3678 ATH_MSG_INFO(
"REGTEST Tau Cluster " << clusCount <<
" pt = " << (*clusItr)->pt()
3679 <<
" eta = " << (*clusItr)->eta()
3680 <<
" phi = " << (*clusItr)->phi() );
3685 for (
unsigned int trackNum = 0; trackNum < (*tauIt)->nTracks(); ++trackNum) {
3690 return StatusCode::SUCCESS;
3695 ATH_MSG_INFO(
"REGTEST Tau linked track " << trackNum <<
" pt = " << linkTrack->
pt()
3696 <<
" eta = " << linkTrack->
eta()
3697 <<
" phi = " << linkTrack->
phi() );
3704 return StatusCode::SUCCESS;
3713 ATH_MSG_INFO(
"REGTEST ==========START of TauJetContainer DUMP===========");
3719 std::string TauContainerTags[]={
"HLT_TrigTauRecMerged"};
3720 for (
int itag=0; itag < ntag; itag++) {
3722 sCode=evtStore()->retrieve(TauJetcont , TauContainerTags[itag]);
3723 if( sCode.isFailure() ){
3724 ATH_MSG_INFO(
"Failed to retrieve TauJetContainer with key " << TauContainerTags[itag]);
3730 bool IsTaurec =
false;
3731 bool Is1p3p =
false;
3752 ATH_MSG_INFO(
"REGTEST (*tauIt)->eta() returns " << (*tauIt)->eta());
3753 ATH_MSG_INFO(
"REGTEST (*tauIt)->phi() returns " << (*tauIt)->phi());
3754 ATH_MSG_INFO(
"REGTEST (*tauIt)->numTrack() returns " << (*tauIt)->numTrack());
3758 ATH_MSG_INFO(
"REGTEST TauDetails->ipZ0SinThetaSigLeadTrk() returns " <<
TauDetails->ipZ0SinThetaSigLeadTrk());
3771 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_stripWidth2() returns " <<
TauDetails->seedCalo_stripWidth2());
3773 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_etEMAtEMScale() returns " <<
TauDetails->seedCalo_etEMAtEMScale());
3774 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_etHADAtEMScale() returns " <<
TauDetails->seedCalo_etHadAtEMScale());
3776 ATH_MSG_INFO(
"REGTEST TauDetails->seedCalo_nIsolLooseTrk() returns " <<
TauDetails->seedCalo_nIsolLooseTrk());
3785 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_etChrgHadOverSumTrkPt () returns " <<
TauDetails->seedTrk_etChrgHadOverSumTrkPt ());
3786 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_isolFracWide() returns " <<
TauDetails->seedTrk_isolFracWide());
3787 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_etHadAtEMScale() returns " <<
TauDetails->seedTrk_etHadAtEMScale());
3788 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_etEMAtEMScale() returns " <<
TauDetails->seedTrk_etEMAtEMScale());
3794 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_sumEMCellEtOverLeadTrkPt() returns " <<
TauDetails->seedTrk_sumEMCellEtOverLeadTrkPt());
3795 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_secMaxStripEt() returns " <<
TauDetails->seedTrk_secMaxStripEt());
3799 ATH_MSG_INFO(
"REGTEST TauDetails->seedTrk_nOtherCoreTrk() returns " <<
TauDetails->seedTrk_nOtherCoreTrk());
3813 ATH_MSG_INFO(
"REGTEST TauDetails->looseTrk() link present ");
3815 ATH_MSG_INFO(
"REGTEST TauDetails->looseTrk() link missing ");
3852 if (
sc.isFailure()) {
3854 return StatusCode::FAILURE;
3864 return StatusCode::SUCCESS;
3871 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrackParticle DUMP===========");
3873 std::vector<std::string> SGkeys;
3874 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Bjet_EFID");
3875 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Bphysics_EFID");
3876 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Electron_EFID");
3877 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_FullScan_EFID");
3878 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Muon_EFID");
3879 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Photon_EFID");
3880 SGkeys.push_back(
"HLT_xAOD__TrackParticleContainer_InDetTrigTrackingxAODCnv_Tau_EFID");
3887 if (
sc.isFailure()) {
3888 ATH_MSG_INFO(
"REGTEST No track particle container found with key " << SGkeys.at(
SGkey));
3891 ATH_MSG_INFO(
"REGTEST TrackParticleContainer retrieved with key " << SGkeys.at(
SGkey)
3892 <<
" and size " << trackParticleContainer->
size());
3897 for (
int index = 0; trackParticleItr != trackParticleLast; ++trackParticleItr, ++
index) {
3901 ATH_MSG_INFO(
"REGTEST pt: " << (*trackParticleItr)->pt()
3902 <<
"/eta: " << (*trackParticleItr)->eta()
3903 <<
"/phi: " << (*trackParticleItr)->phi()
3904 <<
"/m: " << (*trackParticleItr)->m()
3905 <<
"/e: " << (*trackParticleItr)->e()
3906 <<
"/rapidity: " << (*trackParticleItr)->rapidity());
3908 ATH_MSG_INFO(
"REGTEST Defining parameters functions:");
3909 ATH_MSG_INFO(
"REGTEST charge: " << (*trackParticleItr)->charge()
3910 <<
"/d0: " << (*trackParticleItr)->d0()
3911 <<
"/z0: " << (*trackParticleItr)->z0()
3912 <<
"/phi0: " << (*trackParticleItr)->phi0()
3913 <<
"/theta: " << (*trackParticleItr)->theta()
3914 <<
"/qOverP: " << (*trackParticleItr)->qOverP()
3915 <<
"/vx: " << (*trackParticleItr)->vx()
3916 <<
"/vy: " << (*trackParticleItr)->vy()
3917 <<
"/vz: " << (*trackParticleItr)->vz());
3922 ATH_MSG_INFO(
"REGTEST chiSquared: " << (*trackParticleItr)->chiSquared()
3923 <<
"/numberDoF: " << (*trackParticleItr)->numberDoF());
3928 msg() << MSG::INFO <<
"REGTEST ";
3933 msg() <<
"/numberOfBLayerHits not found";
3940 msg() <<
"/numberOfPixelHits not found";
3947 msg() <<
"/numberOfPixelHoles not found";
3954 msg() <<
"/numberOfSCTHits not found";
3961 msg() <<
"/numberOfSCTHoles not found";
3968 msg() <<
"/numberOfTRTHits not found";
3975 msg() <<
"/numberOfTRTHoles not found";
3981 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrackParticle DUMP===========");
3991 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::Vertex DUMP===========");
3995 if (
sc.isFailure()) {
3997 return StatusCode::FAILURE;
4004 for (
int index = 0; vertexItr != vertexLast; ++vertexItr, ++
index) {
4009 <<
"/y: " << (*vertexItr)->y()
4010 <<
"/z: " << (*vertexItr)->z());
4013 ATH_MSG_INFO(
"REGTEST chiSquared: " << (*vertexItr)->chiSquared()
4014 <<
"/numberDoF: " << (*vertexItr)->numberDoF());
4017 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::Vertex DUMP===========");
4020 return StatusCode::SUCCESS;
4025 ATH_MSG_INFO(
"REGTEST ==========START of TDT DUMP===========" );
4028 if (
chain.empty()) {
4031 std::vector<std::string> confChains =
m_trigDec->getListOfTriggers(
chain);
4032 for (
const auto&
item : confChains) {
4035 if (
m_trigDec->getNavigationFormat() ==
"TriggerElement") {
4040 if (passFeatures.size()) {
4041 ATH_MSG_INFO(
" " <<
item <<
" Passed Final IParticle features size: " << passFeatures.size());
4043 if (!li.isValid()) {
4047 std::string state =
"ACTIVE";
4050 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);
4058 if (passAndFailFeatures.size()) {
4059 ATH_MSG_INFO(
" " <<
item <<
" Passed+Failed Final IParticle features size: " << passAndFailFeatures.size());
4061 if (!li.isValid()) {
4065 std::string state =
"ACTIVE";
4068 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);
4075 std::vector< LinkInfo<xAOD::IParticleContainer> > allFeatures =
m_trigDec->features<
xAOD::IParticleContainer>(
item, TrigDefs::includeFailedDecisions,
"", TrigDefs::allFeaturesOfType);
4076 if (allFeatures.size()) {
4077 ATH_MSG_INFO(
" " <<
item <<
" Passed+Failed ALL IParticle features size: " << allFeatures.size());
4079 if (!li.isValid()) {
4083 std::string state =
"ACTIVE";
4086 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);
4095 if (
m_trigDec->getNavigationFormat() ==
"TrigComposite") {
4097 std::vector< LinkInfo<xAOD::IParticleContainer> > muons =
m_trigDec->features<
xAOD::IParticleContainer>(
"HLT_mu24_idperf_L1MU20", TrigDefs::Physics,
"HLT_MuonL2CBInfo");
4104 ATH_MSG_INFO(
"Muon pT: " << (*
mu.link)->pt() <<
" is from the same ROI as tracks with index "
4106 <<
", which is " <<
std::distance(startIt, stopIt) <<
" tracks, out of " << muonTracksReadHandle->
size() <<
" total tracks.");
4113 ATH_MSG_INFO(
"REGTEST ==========END of TDT DUMP===========" );
4114 return StatusCode::SUCCESS;
4118 ATH_MSG_INFO(
"REGTEST ==========START of xAOD::TrigCompositeContainer DUMP===========" );
4124 std::string typeNameTC;
4133 if ( not
evtStore()->contains<xAOD::TrigCompositeContainer>(
key) ) {
4142 for (
auto tc: *cont ) {
4149 ATH_MSG_INFO(
"REGTEST ==========END of xAOD::TrigCompositeContainer DUMP===========" );
4150 return StatusCode::SUCCESS;
4163 if (!elementLink.
isValid())
ATH_MSG_WARNING(
" Invalid element link to TrigRoiDescriptorCollection, link name:'" <<
name <<
"'");
4164 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"'' to TrigRoiDescriptor:" << *elementLink);
4171 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"' to LVL1::RecEmTauRoI:" << *elementLink);
4177 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"' to TrigComposite, TC name:'" << (*elementLink)->name() <<
"'");
4183 else ATH_MSG_DEBUG(
" Dereferenced link '" <<
name <<
"' to View:'" << *elementLink);
4185 }
else if (
name ==
"feature") {
4190 if (!elementLink.
isValid())
ATH_MSG_WARNING(
" Invalid element link to xAOD::TrigEMClusterContainer 'feature'");
4191 else ATH_MSG_DEBUG(
" Dereferenced xAOD::TrigEMClusterContainer link 'feature', Energy:" << (*elementLink)->energy());
4196 if (!elementLink.
isValid())
ATH_MSG_WARNING(
" Invalid element link to xAOD::TrigMissingETContainer 'feature'");
4197 else ATH_MSG_DEBUG(
" Dereferenced xAOD::TrigMissingETContainer link 'feature', ex:" << (*elementLink)->ex() <<
" ey:" << (*elementLink)->ey());
4204 else ATH_MSG_DEBUG(
" Dereferenced IParticle link 'feature', pt:" << (*elementLink)->pt() <<
" eta:" << (*elementLink)->eta() <<
" phi:" << (*elementLink)->phi());
4205 }
catch(std::runtime_error&
e) {
4206 ATH_MSG_WARNING(
" Cannot dereference 'feature' as IParticle: '" <<
e.what() <<
"'");
4215 return StatusCode::SUCCESS;
4226 std::vector<std::string>
keys;
4233 std::string typeNameTC;
4239 chainIDs.insert( chainID.
numeric() );
4244 pass = cg->
isPassed(TrigDefs::requireDecision);
4249 chainIDs.insert( chainID_tmp.
numeric() );
4251 if (legMultiplicites.size() == 0) {
4252 ATH_MSG_ERROR(
"chain " << chainID_tmp <<
" has invalid configuration, no multiplicity data.");
4253 }
else if (legMultiplicites.size() > 1) {
4256 for (
size_t legNumeral = 0; legNumeral < legMultiplicites.size(); ++legNumeral) {
4258 chainIDs.insert( legID.
numeric() );
4267 std::stringstream
ss;
4268 ss <<
"digraph {" << std::endl;
4269 ss <<
" node [shape=rectangle]" << std::endl;
4270 ss <<
" rankdir = BT" << std::endl;
4273 for (
const std::string&
key :
keys) {
4280 ATH_MSG_DEBUG(
"Processing collection " <<
key <<
" to be added to the navigation graph");
4282 bool writtenHeader =
false;
4283 for (
const Decision* tc : *container ) {
4287 std::vector<ElementLink<DecisionContainer>> seedELs = tc->objectCollectionLinks<
DecisionContainer>(
"seed");
4288 const bool isHypoAlgNode = tc->name() ==
"H";
4289 const bool isComboHypoAlgNode = tc->name() ==
"CH";
4290 const std::vector<DecisionID>&
decisions = tc->decisions();
4294 bool doDump =
false;
4297 if (chainIDs.count(
id) == 1) {
4305 const std::vector<DecisionID>& seedDecisions = (*s)->decisions();
4307 if (chainIDs.count(
id) == 1) {
4318 if (!writtenHeader) {
4319 writtenHeader =
true;
4320 ss <<
" subgraph " <<
key <<
" {" << std::endl;
4321 ss <<
" label=\"" <<
key <<
"\"" << std::endl;
4323 const std::string
scheme =
"rdpu9";
4324 std::string
color =
"1";
4325 if (tc->name() ==
"L1") {
color =
"1"; }
4326 else if (tc->name() ==
"F") {
color =
"2"; }
4327 else if (tc->name() ==
"IM") {
color =
"3"; }
4328 else if (tc->name() ==
"H") {
color =
"4"; }
4329 else if (tc->name() ==
"CH") {
color =
"5"; }
4330 else if (tc->name() ==
"SF") {
color =
"6"; }
4331 else if (tc->name() ==
"HLTPassRaw") {
color =
"7"; }
4332 ss <<
" \"" << selfKey <<
"_" << selfIndex <<
"\" [colorscheme="<<
scheme<<
",style=filled,fillcolor="<<
color<<
",label=<<B>Container</B>=" << typeNameTC;
4333 if (tc->name() !=
"")
ss <<
" <B>Name</B>=" << tc->name();
4334 ss <<
"<BR/><B>Key</B>=" <<
key <<
"<BR/><B>Index</B>=" << selfIndex;
4335 const bool isRemapped = tc->isRemapped();
4336 if (isHypoAlgNode)
ss <<
" <B>linksRemapped</B>=" << (isRemapped ?
"Y" :
"N");
4338 ss <<
"<BR/><B>Pass</B>=";
4342 std::string highlight = (dID.
numeric() == chainID.
numeric() ?
"<B>[CHAIN:" :
"");
4346 ss << std::hex << highlight << decisionID << (!highlight.empty() ?
"]</B>" :
"") << std::dec <<
",";
4353 ss <<
">]" << std::endl;
4355 size_t seedCount = 0;
4356 for (
size_t i = 0;
i < tc->linkColNames().
size(); ++
i) {
4357 const std::string link = tc->linkColNames().
at(
i);
4358 if (link ==
"seed" || link ==
"seed__COLL") {
4361 const uint32_t seedIndex = tc->linkColIndices().
at(
i);
4365 const std::vector<DecisionID> seedDecisions = (*seedEL)->decisions();
4366 bool doSeedLink =
false;
4368 if (chainIDs.count(
id) == 1) {
4377 ss <<
" \"" << selfKey <<
"_" << selfIndex <<
"\" -> \"" << seedKey <<
"_" << seedIndex <<
"\" [colorscheme="<<
scheme<<
",color=9,fontcolor=8,label=\"seed\"]" << std::endl;
4380 std::string linkColour =
"12";
4381 std::string linkBackground =
"11";
4382 const std::string extScheme =
"paired12";
4383 if (link ==
"roi") { linkColour=
"2"; linkBackground=
"1"; }
4384 else if (link ==
"initialRoI") { linkColour=
"2"; linkBackground=
"1"; }
4385 else if (link ==
"initialRecRoI") { linkColour=
"8"; linkBackground=
"7"; }
4386 else if (link ==
"feature") { linkColour=
"4"; linkBackground=
"3"; }
4387 else if (link ==
"view") { linkColour=
"10"; linkBackground=
"9"; }
4388 const CLID linkCLID =
static_cast<CLID>( tc->linkColClids().
at(
i) );
4394 const unsigned index = (isRemapped ? tc->linkColIndicesRemap().
at(
i) : tc->linkColIndices().
at(
i));
4397 const std::string* keyStr =
evtStore()->keyToString(
key, checkCLID);
4398 if (keyStr !=
nullptr && checkCLID != linkCLID) {
4399 std::string tnameOfCheck;
4400 m_clidSvc->getTypeNameOfID(checkCLID, tnameOfCheck).ignore();
4401 ATH_MSG_ERROR(
"Inconsistent CLID " << checkCLID <<
" [" << tnameOfCheck <<
"] stored in storegate for key " <<
key
4402 <<
". We were expecting " << linkCLID <<
" [" << tname <<
"]");
4405 std::string tnameEscape;
4406 for (std::string::const_iterator
i = tname.begin();
i != tname.end(); ++
i) {
4407 unsigned char c = *
i;
4409 tnameEscape +=
"<";
4410 }
else if (
c ==
'>') {
4411 tnameEscape +=
">";
4418 ss <<
" \"" << selfKey <<
"_" << selfIndex <<
"\" -> \"" <<
key <<
"_" <<
index <<
"\" ";
4419 ss <<
"[colorscheme="<<extScheme<<
",color="<<linkColour<<
",fontcolor="<<linkColour<<
",arrowhead=empty,label=\"" << link <<
"\"]" << std::endl;
4422 const bool linkToSelf = (selfKey ==
key and selfIndex ==
index);
4425 ss <<
" \"" <<
key <<
"_" <<
index <<
"\" [colorscheme="<<extScheme<<
",style=filled,fillcolor="<<linkBackground<<
",label=<<B>Container</B>=" << tnameEscape <<
"<BR/><B>Key</B>=";
4426 if (keyStr !=
nullptr)
ss << *keyStr;
4427 else ss <<
"[<I>KEY "<<
key <<
" NOT IN STORE</I>] ";
4428 ss <<
"<BR/><B>Index</B>=" <<
index <<
">]";
4435 if (writtenHeader) {
4436 ss <<
" }" << std::endl;
4440 ss <<
"}" << std::endl;
4442 returnValue.assign(
ss.str() );
4443 return StatusCode::SUCCESS;
4456 const std::vector< unsigned int > serialisedNavigation = navigationHandle->
serialized();
4457 ATH_MSG_INFO(
"Serialised navigation size: " << serialisedNavigation.size() );
4466 std::map< int, std::string > hash2string;
4467 for (
auto const& sequence :
testMenu->front()->sequenceInputTEs() ) {
4468 for (
auto const&
name : sequence ) {
4475 unsigned int chainCounter = 0;
4476 std::map< int, std::string > hash2chain;
4477 for (
auto const&
chain :
testMenu->front()->chainSignatureOutputTEs() ) {
4484 for (
auto const& signature :
chain ) {
4485 for (
auto const&
name : signature ) {
4497 if ( lhs.
getCLID() == rhs.getCLID() )
return ( lhs.
getIndex() < rhs.getIndex() );
4500 else return ( lhs.
getCLID() < rhs.getCLID() );
4502 std::map< HLT::TriggerElement::FeatureAccessHelper, std::vector< HLT::TriggerElement* >, decltype(cmpLambda) > feature2element(cmpLambda);
4505 std::vector< HLT::TriggerElement* > allTEs;
4506 testNav->
getAll( allTEs,
false );
4507 for (
auto element : allTEs ) {
4510 for (
auto helper : element->getFeatureAccessHelpers() ) {
4511 feature2element[
helper ].push_back( element );
4517 for (
auto element : allTEs ) {
4519 ATH_MSG_INFO(
"+++++++++++ " << hash2string[ element->getId() ] <<
" is terminal node" );
4521 std::queue< HLT::TriggerElement* > allAncestors;
4522 allAncestors.push( element );
4523 while ( allAncestors.size() ) {
4527 auto theseAncestors = thisElement->
getRelated( HLT::TriggerElement::Relation::seededByRelation );
4538 for (
auto ancestor : theseAncestors ) {
4539 allAncestors.push( ancestor );
4547 auto decisionOutput = outputNavigation.
ptr();
4550 std::map< HLT::TriggerElement const*, std::vector< int > > element2decisions;
4551 for (
const auto& pair : feature2element ) {
4554 std::string featureName = testNav->
label( pair.first.getCLID(), pair.first.getIndex().subTypeIndex() );
4555 auto sgKey =
evtStore()->stringToKey( featureName, pair.first.getCLID() );
4558 std::string storeFeatureName =
"feature";
4565 decision->typelessSetObjectLink( storeFeatureName, sgKey, pair.first.getCLID(), pair.first.getIndex().objectsBegin(), pair.first.getIndex().objectsEnd() );
4568 std::set< std::string > passedChains;
4572 passedChains.insert( hash2chain[ element->
getId() ] );
4575 int decisionNumber = decisionOutput->size() - 1;
4576 element2decisions[ element ].push_back( decisionNumber );
4580 for (
auto&
chain : passedChains ) {
4586 unsigned int decisionCounter = 0;
4587 for (
const auto& pair : feature2element ) {
4590 auto decision = decisionOutput->at( decisionCounter );
4594 for (
auto element : pair.second ) {
4595 auto theseAncestors = element->getRelated( HLT::TriggerElement::Relation::seededByRelation );
4596 for (
auto ancestor : theseAncestors ) {
4597 for (
int decisionIndex : element2decisions[ ancestor ] ) {
4604 return StatusCode::SUCCESS;