EventInfoCnvTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
// Gaudi/Athena include(s):
#include "AthenaKernel/errorcheck.h"
 
// EDM include(s):
#include "EventInfo/EventInfo.h"
#include "EventInfo/EventID.h"
#include "EventInfo/EventType.h"
#include "EventInfo/TriggerInfo.h"
#include "EventInfo/PileUpEventInfo.h"
#include "AthenaPoolUtilities/AthenaAttributeList.h"
#include "AthenaPoolUtilities/CondAttrListCollection.h"
#include "xAODEventInfo/EventInfo.h"
#include "xAODEventInfo/EventInfoContainer.h"
#include "StoreGate/ReadCondHandle.h"
 
// Local include(s):
#include "EventInfoCnvTool.h"
 
namespace xAODMaker {
 
  // maps to copy the flags
  static const std::map< xAOD::EventInfo::EventFlagSubDet,
                         EventInfo::EventFlagSubDet > subDetMap = 
    { {xAOD::EventInfo::Pixel,      EventInfo::Pixel},
      {xAOD::EventInfo::SCT,        EventInfo::SCT},
      {xAOD::EventInfo::TRT,        EventInfo::TRT},
      {xAOD::EventInfo::LAr,        EventInfo::LAr},
      {xAOD::EventInfo::Tile,       EventInfo::Tile},
      {xAOD::EventInfo::Muon,       EventInfo::Muon},
      {xAOD::EventInfo::ForwardDet, EventInfo::ForwardDet},
      {xAOD::EventInfo::Core,       EventInfo::Core},
      {xAOD::EventInfo::Background, EventInfo::Background},
      {xAOD::EventInfo::Lumi,       EventInfo::Lumi} };
 
  static const std::map< EventInfo::EventFlagErrorState,
                         xAOD::EventInfo::EventFlagErrorState > errorStateMap =
    { { EventInfo::NotSet,  xAOD::EventInfo::NotSet },
      { EventInfo::Warning, xAOD::EventInfo::Warning },
      { EventInfo::Error,   xAOD::EventInfo::Error}
    };
 
#if !defined(XAOD_ANALYSIS) && !defined(GENERATIONBASE)
   static const std::string INDET_BEAMPOS = "/Indet/Beampos";
#endif
 
 
   EventInfoCnvTool::EventInfoCnvTool( const std::string& type,
                                       const std::string& name,
                                       const IInterface* parent )
      : AthAlgTool( type, name, parent )
   {
      // Declare the interface(s) provided by the tool:
      declareInterface< IEventInfoCnvTool >( this );
   }
 
   StatusCode EventInfoCnvTool::initialize()
   {
 
#if !defined(XAOD_ANALYSIS) && !defined(GENERATIONBASE)
      // Check if the beam position will be available or not:
      if( detStore()->contains< AthenaAttributeList >( INDET_BEAMPOS ) ) {
         m_beamSpotInformationAvailable = true;
      } else {
         // This is normal in some cases, but tell the user anyway
         ATH_MSG_INFO( "Beam spot information not available" );
         ATH_MSG_INFO( "Will not fill beam spot information into "
                       "xAOD::EventInfo" );
         m_beamSpotInformationAvailable = false;
      }
      if(m_disableBeamSpot){
         ATH_MSG_WARNING( "Beam spot information manually disabled on EventInfo object" );
         m_beamSpotInformationAvailable = false;
      }
 
      ATH_CHECK(m_beamSpotKey.initialize(m_beamSpotInformationAvailable));
#endif // not XAOD_ANALYSIS or GENERATIONBASE
 
#if !defined(XAOD_ANALYSIS) && !defined(SIMULATIONBASE) && !defined(GENERATIONBASE)
      CHECK( m_lumiDataKey.initialize (SG::AllowEmpty) );
#endif // not XAOD_ANALYSIS or GENERATIONBASE
 
      // Return gracefully:
      return StatusCode::SUCCESS;
   }
 
   StatusCode EventInfoCnvTool::convert( const EventInfo* aod,
                                         xAOD::EventInfo* xaod,
                                         bool pileUpInfo,
                                         bool copyPileUpLinks,
                                         [[maybe_unused]] const EventContext& ctx /*= Gaudi::Hive::currentContext()*/) const
   {
 
      if( ! aod ) {
         ATH_MSG_WARNING( "Null pointer received for input!" );
         return StatusCode::SUCCESS;
      }
 
      // Copy the basic properties into the xAOD object:
      if( aod->event_ID() ) {
         xaod->setRunNumber( aod->event_ID()->run_number() );
         xaod->setEventNumber( aod->event_ID()->event_number() );
         xaod->setLumiBlock( aod->event_ID()->lumi_block() );
         xaod->setTimeStamp( aod->event_ID()->time_stamp() );
         xaod->setTimeStampNSOffset( aod->event_ID()->time_stamp_ns_offset() );
         xaod->setBCID( aod->event_ID()->bunch_crossing_id() );
         xaod->setDetectorMask( aod->event_ID()->detector_mask0(),
                                aod->event_ID()->detector_mask1() );
         xaod->setDetectorMaskExt( aod->event_ID()->detector_mask2(),
                                 aod->event_ID()->detector_mask3() );
      }
 
      // Copy the event type properties:
      std::vector< float > eventWeights;
      if( aod->event_type() ) {
         EventType::NameTagPairVec detDescrTags;
         aod->event_type()->get_detdescr_tags( detDescrTags );
         xaod->setDetDescrTags( detDescrTags );
         uint32_t eventTypeBitmask = 0;
         if( aod->event_type()->test( EventType::IS_SIMULATION ) ) {
            eventTypeBitmask |= xAOD::EventInfo::IS_SIMULATION;
         }
         if( aod->event_type()->test( EventType::IS_TESTBEAM ) ) {
            eventTypeBitmask |= xAOD::EventInfo::IS_TESTBEAM;
         }
         if( aod->event_type()->test( EventType::IS_CALIBRATION ) ) {
            eventTypeBitmask |= xAOD::EventInfo::IS_CALIBRATION;
         }
         xaod->setEventTypeBitmask( eventTypeBitmask );
         // Only add MC information for simulation files:
         if( xaod->eventType( xAOD::EventInfo::IS_SIMULATION ) ) {
            xaod->setMCChannelNumber( aod->event_type()->mc_channel_number() );
            xaod->setMCEventNumber( aod->event_type()->mc_event_number() );
            eventWeights.resize( aod->event_type()->n_mc_event_weights(), 0.0 );
            for( unsigned int i = 0; i < aod->event_type()->n_mc_event_weights();
                 ++i ) {
               eventWeights[ i ] = aod->event_type()->mc_event_weight( i );
            }
         }
      }
 
      // Copy the trigger properties into the xAOD object:
      if( aod->trigger_info() && !pileUpInfo  ) {
         xaod->setStatusElement( aod->trigger_info()->statusElement() );
         xaod->setExtendedLevel1ID( aod->trigger_info()->extendedLevel1ID() );
         xaod->setLevel1TriggerType( aod->trigger_info()->level1TriggerType() );
         std::vector< xAOD::EventInfo::StreamTag > streamTags;
         std::vector< TriggerInfo::StreamTag >::const_iterator st_itr =
            aod->trigger_info()->streamTags().begin();
         std::vector< TriggerInfo::StreamTag >::const_iterator st_end =
            aod->trigger_info()->streamTags().end();
         for( ; st_itr != st_end; ++st_itr ) {
            streamTags.emplace_back( st_itr->name(), st_itr->type(),
                                            st_itr->obeysLumiblock(),
                                            st_itr->robs(), st_itr->dets() );
         }
         xaod->setStreamTags( streamTags );
      }
 
      // Copy/calculate the pileup information:
      if( ! pileUpInfo ) {
         bool haveLumi = false;
#if !defined(XAOD_ANALYSIS) && !defined(SIMULATIONBASE) && !defined(GENERATIONBASE)
         if (!m_lumiDataKey.empty()) {
           SG::ReadCondHandle<LuminosityCondData> lumiData (m_lumiDataKey, ctx);
           if (lumiData->lbAverageLuminosity() != 0 ||
               lumiData->lbAverageInteractionsPerCrossing() != 0)
           {
             float actualMu = 0.0;
             const float muToLumi = lumiData->muToLumi();
             if( std::abs( muToLumi ) > 0.00001 ) {
               unsigned int bcid = ctx.eventID().bunch_crossing_id();
               actualMu = lumiData->lbLuminosityPerBCIDVector().at(bcid) / muToLumi;
             }
             xaod->setActualInteractionsPerCrossing( actualMu );
             xaod->setAverageInteractionsPerCrossing(
               lumiData->lbAverageInteractionsPerCrossing() );
           }
           haveLumi = true;
         }
#endif // not XAOD_ANALYSIS or SIMULATIONBASE or GENERATIONBASE
         if (!haveLumi) {
            xaod->setActualInteractionsPerCrossing(
               aod->actualInteractionsPerCrossing() );
            xaod->setAverageInteractionsPerCrossing(
               aod->averageInteractionsPerCrossing() );
         }
      }
 
      // Copy the sub-detector flags:
      std::map< xAOD::EventInfo::EventFlagSubDet,
                EventInfo::EventFlagSubDet >::const_iterator sd_itr =
         subDetMap.begin();
      std::map< xAOD::EventInfo::EventFlagSubDet,
                EventInfo::EventFlagSubDet >::const_iterator sd_end =
         subDetMap.end();
      for( ; sd_itr != sd_end; ++sd_itr ) {
 
         // Lumi does not store event flags, or an error state:
        if( (int)sd_itr->first == (int)EventInfo::Lumi ) {
            continue;
         }
 
         // Set the event flags for this sub-detector:
         xaod->setEventFlags( sd_itr->first, aod->eventFlags( sd_itr->second ) );
 
         // Look up the error state of this sub-detector:
         std::map< EventInfo::EventFlagErrorState,
                   xAOD::EventInfo::EventFlagErrorState >::const_iterator state =
            errorStateMap.find( aod->errorState( sd_itr->second ) );
         if( state == errorStateMap.end() ) {
            REPORT_MESSAGE( MSG::FATAL )
               << "Unknown error state found for sub-detector "
               << sd_itr->second << ": " << aod->errorState( sd_itr->second );
            return StatusCode::FAILURE;
         }
 
         // Set the error state for the sub-detector:
         xaod->setErrorState( sd_itr->first, state->second );
      }
 
      // Check if it is a PileUpEventInfo object:
      const PileUpEventInfo* puei =
         dynamic_cast< const PileUpEventInfo* >( aod );
      if( puei && copyPileUpLinks ) {
         // Create the sub-event objects to fill into the output object:
         std::vector< xAOD::EventInfo::SubEvent > subEvents;
         PileUpEventInfo::SubEvent::const_iterator itr = puei->beginSubEvt();
         PileUpEventInfo::SubEvent::const_iterator end = puei->endSubEvt();
         for( ; itr != end; ++itr ) {
            // Look up the sub-event type:
           
            xAOD::EventInfo::PileUpType type = xAOD::EventInfo::Unknown;
            switch (itr->type()) {
            case PileUpTimeEventIndex::Signal:
              type = xAOD::EventInfo::Signal;
              break;
            case PileUpTimeEventIndex::MinimumBias:
              type = xAOD::EventInfo::MinimumBias;
              break;
            case PileUpTimeEventIndex::Cavern:
              type = xAOD::EventInfo::Cavern;
              break;
            case PileUpTimeEventIndex::HaloGas:
              type = xAOD::EventInfo::HaloGas;
              break;
            case PileUpTimeEventIndex::HighPtMinimumBias:
              type = xAOD::EventInfo::HighPtMinimumBias;
              break;
            case PileUpTimeEventIndex::ZeroBias:
              type = xAOD::EventInfo::ZeroBias;
              break;
            case PileUpTimeEventIndex::Unknown:
              break;
            default:
              ATH_MSG_WARNING( "Unknown sub-event type ("
                               << itr->type() << ") encountered" );
              break;
            }
 
            // Construct the link to the pile-up EventInfo object:
            ElementLink< xAOD::EventInfoContainer > link;
            link.resetWithKeyAndIndex( "PileUpEventInfo", itr->index() );
 
            // Add the new object
            subEvents.emplace_back( itr->time(),
                                                            itr->index(),
                                                            type, link );
         }
 
         // Add the vector to the EventInfo:
         xaod->setSubEvents( subEvents );
      }
 
#if !defined(XAOD_ANALYSIS) && !defined(GENERATIONBASE)
      // Fill the beam spot variables if the necessary service is available:
      if( m_beamSpotInformationAvailable && ( ! pileUpInfo ) ) {
         SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey, ctx };
         xaod->setBeamPos( beamSpotHandle->beamPos()[ Amg::x ],
                           beamSpotHandle->beamPos()[ Amg::y ],
                           beamSpotHandle->beamPos()[ Amg::z ] );
         xaod->setBeamPosSigma( beamSpotHandle->beamSigma( 0 ),
                                beamSpotHandle->beamSigma( 1 ),
                                beamSpotHandle->beamSigma( 2 ) );
         xaod->setBeamPosSigmaXY( beamSpotHandle->beamSigmaXY() );
         xaod->setBeamTiltXZ( beamSpotHandle->beamTilt( 0 ) );
         xaod->setBeamTiltYZ( beamSpotHandle->beamTilt( 1 ) );
         xaod->setBeamStatus( beamSpotHandle->beamStatus() );
      }
#endif // not XAOD_ANALYSIS or GENERATIONBASE
 
      if (!eventWeights.empty()) {
        if (xaod->getStore()) {
          // Ensure that this gets added as a decoration.
          xaod->getStore()->lock();
        }
        static const xAOD::EventInfo::Decorator<std::vector<float> > ew ("mcEventWeights");
        ew (*xaod) = std::move (eventWeights);
      }
 
      // Finish with some printout:
      ATH_MSG_DEBUG( "Finished conversion EventInfo="<<aod<<" xAOD::EventInfo="<<xaod<<" content=" << *xaod );
 
      // Return gracefully:
      return StatusCode::SUCCESS;
   }
 
} // namespace xAODMaker