Merge branch 'HLTResultMT-decoder' into 'master'

HLTResultMT ByteStream decoder algorithm

See merge request atlas/athena!16284

HLT/Trigger/TrigControl/TrigExamples/TrigExPartialEB/share/MTCalibPeb.py

@@ -85,7 +85,8 @@ from TrigOutputHandling.TrigOutputHandlingConf import TriggerEDMSerialiserTool,
 
 # Tool serialising EDM objects to fill the HLT result
 serialiser = TriggerEDMSerialiserTool()
-serialiser.CollectionsToSerialize = [] # nothing at the moment
+serialiser.CollectionsToSerialize = ["xAOD::TrigCompositeContainer_v1#MTCalibPebDecisions",
+                                     "xAOD::TrigCompositeAuxContainer_v1#MTCalibPebDecisionsAux."]
 
 # Tool adding stream tags to HLT result
 stmaker = StreamTagMakerTool()

Trigger/TrigEvent/TrigSteeringEvent/TrigSteeringEvent/HLTResultMT.h

@@ -125,11 +125,11 @@ namespace HLT {
     void setStatus(const std::vector<uint32_t>& status);
 
     /// Replace error codes with the given codes
-    void setStatus(const std::vector<uint32_t>& errorCodes,
-                   const eformat::helper::Status firstStatusWord = {
-                     eformat::GenericStatus::DATA_CORRUPTION,
-                     eformat::FullEventStatus::PSC_PROBLEM
-                   });
+    void setErrorCodes(const std::vector<uint32_t>& errorCodes,
+                       const eformat::helper::Status firstStatusWord = {
+                         eformat::GenericStatus::DATA_CORRUPTION,
+                         eformat::FullEventStatus::PSC_PROBLEM
+                       });
 
     /** @brief Append an error code
      *
@@ -166,6 +166,9 @@ namespace HLT {
   };
 } // namespace HLT
 
+/// operator<< overload for printing to std::ostream
+std::ostream& operator<<(std::ostream& str, const HLT::HLTResultMT& hltResult);
+
 CLASS_DEF(HLT::HLTResultMT, 172156324, 1)
 
 #endif // TRIGSTEERINGEVENT_HLTResultMT_H

Trigger/TrigEvent/TrigSteeringEvent/src/HLTResultMT.cxx

@@ -189,8 +189,8 @@ void HLT::HLTResultMT::setStatus(const std::vector<uint32_t>& status) {
 }
 
 // -----------------------------------------------------------------------------
-void HLT::HLTResultMT::setStatus(const std::vector<uint32_t>& errorCodes,
-                                 const eformat::helper::Status firstStatusWord) {
+void HLT::HLTResultMT::setErrorCodes(const std::vector<uint32_t>& errorCodes,
+                                     const eformat::helper::Status firstStatusWord) {
   m_status.clear();
   m_status.push_back(firstStatusWord.code());
   m_status.insert(m_status.end(),errorCodes.cbegin(),errorCodes.cend());
@@ -203,3 +203,52 @@ void HLT::HLTResultMT::addErrorCode(const uint32_t& errorCode,
   else m_status[0] |= firstStatusWord.code();
   m_status.push_back(errorCode);
 }
+
+// =============================================================================
+std::ostream& operator<<(std::ostream& str, const HLT::HLTResultMT& hltResult) {
+  auto printWord = [&str](const uint32_t word, const size_t width=8){
+    str << "0x" << std::hex << std::setw(width) << std::setfill('0') << word << " " << std::dec;
+  };
+  str << "Printing HLTResultMT:" << std::endl;
+
+  // Status
+  str << "--> Status words = ";
+  for (const uint32_t word : hltResult.getStatus()) {
+    printWord(word);
+  }
+  str << std::endl;
+
+  // Stream tags
+  str << "--> Stream tags  = ";
+  bool first = true;
+  for (const eformat::helper::StreamTag& st : hltResult.getStreamTags()) {
+    if (first) first=false;
+    else str << "                   ";
+    str << "{" << st.name << ", " << st.type << ", obeysLB=" << st.obeys_lumiblock << ", robs=[";
+    for (const auto& robid : st.robs) printWord(robid);
+    str << "], dets = [";
+    for (const auto& detid : st.dets) printWord(detid,2);
+    str << "]}" << std::endl;
+  }
+
+  // HLT bits
+  std::vector<uint32_t> hltBitWords;
+  hltBitWords.resize(hltResult.getHltBits().num_blocks());
+  boost::to_block_range(hltResult.getHltBits(),hltBitWords.begin());
+  str << "--> HLT bits     = ";
+  for (const uint32_t word : hltBitWords) {
+    printWord(word);
+  }
+  str << std::endl;
+
+  // Payload size
+  str << "--> Payload size = ";
+  first = true;
+  for (const auto& p : hltResult.getSerialisedData()) {
+    if (first) first=false;
+    else str << "                   ";
+    str << "{module " << p.first << ": " << p.second.size() << " words}" << std::endl;
+  }
+
+  return str;
+}

Trigger/TrigSteer/TrigHLTResultByteStream/src/HLTResultMTByteStreamCnv.cxx

@@ -61,10 +61,9 @@ StatusCode HLT::HLTResultMTByteStreamCnv::finalize() {
 // Implementation of Converter::createObj
 // =============================================================================
 StatusCode HLT::HLTResultMTByteStreamCnv::createObj(IOpaqueAddress* /*pAddr*/, DataObject*& /*pObj*/) {
-  ATH_MSG_VERBOSE("start of " << __FUNCTION__);
-  // Not yet implemented
-  ATH_MSG_VERBOSE("end of " << __FUNCTION__);
-  return StatusCode::SUCCESS;
+  ATH_REPORT_ERROR(StatusCode::FAILURE) << "Using BS converter to decode HLTResultMT is not supported!"
+                                        << " Use HLTResultMTByteStreamDecoderAlg instead";
+  return StatusCode::FAILURE;
 }
 
 // =============================================================================
@@ -124,8 +123,8 @@ StatusCode HLT::HLTResultMTByteStreamCnv::createRep(DataObject* pObj, IOpaqueAdd
   // Loop over all module IDs and fill the ROBFragments
   ATH_MSG_DEBUG("Iterating over modules to assemble output data");
   for (const auto& p : hltResult->getSerialisedData()) {
-    uint16_t moduleId = p.first;
-    std::vector<uint32_t> data = p.second;
+    const uint16_t moduleId = p.first;
+    const std::vector<uint32_t>& data = p.second;
 
     // Get a pointer to ROD data vector to be filled
     eformat::helper::SourceIdentifier sid(eformat::TDAQ_HLT,moduleId);
@@ -137,7 +136,8 @@ StatusCode HLT::HLTResultMTByteStreamCnv::createRep(DataObject* pObj, IOpaqueAdd
 
     // Fill the ROD data vector
     rod->assign(data.cbegin(), data.cend());
-    ATH_MSG_DEBUG("Assembled data for module 0x" << MSG::hex << sid.code() << MSG::dec);
+    ATH_MSG_DEBUG("Assembled data for module 0x" << MSG::hex << sid.code() << MSG::dec << " with "
+                  << data.size() << " words of serialised payload");
   }
 
   ATH_MSG_DEBUG("Appending data to RawEventWrite");

Trigger/TrigSteer/TrigHLTResultByteStream/src/HLTResultMTByteStreamDecoderAlg.cxx

+/*
+  Copyright (C) 2002-2018 CERN for the benefit of the ATLAS collaboration
+*/
+
+#include "HLTResultMTByteStreamDecoderAlg.h"
+
+// =============================================================================
+// Standard constructor
+// =============================================================================
+HLTResultMTByteStreamDecoderAlg::HLTResultMTByteStreamDecoderAlg(const std::string& name, ISvcLocator* svcLoc)
+: AthReentrantAlgorithm(name, svcLoc),
+  m_robDataProviderSvc("ROBDataProviderSvc", name) {}
+
+// =============================================================================
+// Implementation of AthReentrantAlgorithm::initialize
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderAlg::initialize() {
+  ATH_MSG_DEBUG("Initialising " << name());
+  ATH_CHECK(m_hltResultWHKey.initialize());
+  ATH_CHECK(m_robDataProviderSvc.retrieve());
+
+  // Convert module IDs to HLT result ROB IDs
+  for (const uint16_t moduleId : m_moduleIdsToDecode) {
+    eformat::helper::SourceIdentifier sid(eformat::TDAQ_HLT, moduleId);
+    m_robIdsToDecode.push_back(sid.code());
+  }
+
+  return StatusCode::SUCCESS;
+}
+
+// =============================================================================
+// Implementation of AthReentrantAlgorithm::finalize
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderAlg::finalize() {
+  ATH_MSG_DEBUG("Finalising " << name());
+  ATH_CHECK(m_robDataProviderSvc.release());
+  return StatusCode::SUCCESS;
+}
+
+// =============================================================================
+// Implementation of AthReentrantAlgorithm::execute_r
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderAlg::execute_r(const EventContext& eventContext) const {
+  ATH_MSG_DEBUG("Executing " << name());
+
+  // Create and record the HLTResultMT object
+  auto hltResult = SG::makeHandle(m_hltResultWHKey,eventContext);
+  ATH_CHECK( hltResult.record(std::make_unique<HLT::HLTResultMT>()) );
+  ATH_MSG_DEBUG("Recorded HLTResultMT with key " << m_hltResultWHKey.key());
+
+  // Retrieve the RawEvent pointer
+  const RawEvent* re = m_robDataProviderSvc->getEvent(eventContext); // Owned by ROBDataProvider or InputSvc
+  if (!re) {
+    ATH_MSG_ERROR("Failed to retrieve the RawEvent pointer from ROBDataProvider");
+    return StatusCode::FAILURE;
+  }
+
+  // Fill the result object from ByteStream event header
+  ATH_CHECK(m_decoderTool->decodeHeader(re, *hltResult));
+
+  // Read the HLT result payload
+  IROBDataProviderSvc::VROBFRAG vrobf; // vector of ROBFragments to be filled
+  m_robDataProviderSvc->getROBData(eventContext, m_robIdsToDecode, vrobf, name());
+  ATH_CHECK(m_decoderTool->decodePayload(vrobf, *hltResult));
+
+  // Print the result
+  ATH_MSG_DEBUG(*hltResult);
+
+  return StatusCode::SUCCESS;
+}

Trigger/TrigSteer/TrigHLTResultByteStream/src/HLTResultMTByteStreamDecoderAlg.h

+/*
+  Copyright (C) 2002-2018 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef TRIGHLTRESULTBYTESTREAM_HLTResultMTByteStreamDecoderAlg_H
+#define TRIGHLTRESULTBYTESTREAM_HLTResultMTByteStreamDecoderAlg_H
+
+// Trigger includes
+#include "HLTResultMTByteStreamDecoderTool.h"
+#include "TrigSteeringEvent/HLTResultMT.h"
+
+// Athena includes
+#include "AthenaBaseComps/AthReentrantAlgorithm.h"
+#include "ByteStreamCnvSvcBase/IROBDataProviderSvc.h"
+#include "StoreGate/WriteHandle.h"
+
+/** @class HLTResultMTByteStreamDecoderAlg
+ *  @brief Algorithm creating HLTResultMT from ByteStream representation
+ **/
+class HLTResultMTByteStreamDecoderAlg : public AthReentrantAlgorithm {
+public:
+  /// Standard constructor
+  HLTResultMTByteStreamDecoderAlg(const std::string& name, ISvcLocator* svcLoc);
+
+  // ------------------------- AthReentrantAlgorithm methods -------------------
+  virtual StatusCode initialize() override;
+  virtual StatusCode finalize() override;
+  virtual StatusCode execute_r(const EventContext& eventContext) const override;
+
+private:
+  // ------------------------- Properties --------------------------------------
+  /// Module IDs to decode
+  Gaudi::Property<std::vector<uint16_t>> m_moduleIdsToDecode {
+    this, "ModuleIdsToDecode", {0,1,2,3,4,5,6},
+    "List of module IDs from which payload will be decoded"
+  };
+  /// StoreGate key for the output HLTResultMT
+  SG::WriteHandleKey<HLT::HLTResultMT> m_hltResultWHKey {this, "HLTResultWHKey", "HLTResultMT",
+                                                         "Key of the output HLTResultMT object"};
+
+  // ------------------------- Tool/Service handles ----------------------------
+  /// Tool performing the decoding work
+  ToolHandle<HLTResultMTByteStreamDecoderTool> m_decoderTool {this, "DecoderTool", "HLTResultMTByteStreamDecoderTool",
+                                                              "Tool performing the decoding work"};
+  /// ROBDataProvider service handle
+  ServiceHandle<IROBDataProviderSvc> m_robDataProviderSvc;
+
+  // ------------------------- Other private members ---------------------------
+  /// Full ROB IDs constructed from the ModuleIdsToDecode property
+  std::vector<uint32_t> m_robIdsToDecode;
+};
+
+#endif // TRIGHLTRESULTBYTESTREAM_HLTResultMTByteStreamDecoderAlg_H

Trigger/TrigSteer/TrigHLTResultByteStream/src/HLTResultMTByteStreamDecoderTool.cxx

+/*
+  Copyright (C) 2002-2018 CERN for the benefit of the ATLAS collaboration
+*/
+
+// Trigger includes
+#include "HLTResultMTByteStreamDecoderTool.h"
+
+// TDAQ includes
+#include "eformat/StreamTag.h"
+
+// =============================================================================
+// Standard constructor
+// =============================================================================
+HLTResultMTByteStreamDecoderTool::HLTResultMTByteStreamDecoderTool(const std::string& type,
+                                                                   const std::string& name,
+                                                                   const IInterface* parent)
+: AthAlgTool(type,name,parent) {}
+
+// =============================================================================
+// Implementation of AthAlgTool::initialize
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderTool::initialize() {
+  ATH_MSG_DEBUG("Initialising " << name());
+  return StatusCode::SUCCESS;
+}
+
+// =============================================================================
+// Implementation of AthAlgTool::finalize
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderTool::finalize() {
+  ATH_MSG_DEBUG("Finalising " << name());
+  return StatusCode::SUCCESS;
+}
+
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderTool::decodeHeader(const RawEvent* rawEvent, HLT::HLTResultMT& resultToFill) const {
+
+  if (!rawEvent) {
+    ATH_MSG_ERROR("Decoding requested with nullptr RawEvent");
+    return StatusCode::FAILURE;
+  }
+  ATH_MSG_DEBUG("Decoding HLTResultMT from ByteStream event " << rawEvent->global_id());
+
+  // ---------------------------------------------------------------------------
+  // Read the status words (error codes from online event processing)
+  // ---------------------------------------------------------------------------
+  std::vector<uint32_t> statusWords;
+  // try-catch around eformat calls and raw-pointer operations
+  try {
+    const uint32_t nStatus = rawEvent->nstatus();
+    const uint32_t* rawStatus = rawEvent->status(); // owned by rawEvent
+    for (uint32_t i=0; i<nStatus; ++i)
+      statusWords.push_back(rawStatus[i]);
+  }
+  catch (const std::exception& ex) {
+    ATH_MSG_ERROR("std::exception caught when reading status words: " << ex.what());
+    return StatusCode::FAILURE;
+  }
+  catch (...) {
+    ATH_MSG_ERROR("Unknown exception caught when reading status words");
+    return StatusCode::FAILURE;
+  }
+  resultToFill.setStatus(statusWords);
+  ATH_MSG_DEBUG("Successfully read " << statusWords.size() << " status words");
+
+  // ---------------------------------------------------------------------------
+  // Read the stream tags
+  // ---------------------------------------------------------------------------
+  std::vector<eformat::helper::StreamTag> streamTags;
+  // try-catch around eformat calls and raw-pointer operations
+  try {
+    eformat::helper::decode(rawEvent->nstream_tag(),rawEvent->stream_tag(),streamTags);
+  }
+  catch (const std::exception& ex) {
+    ATH_MSG_ERROR("std::exception caught when reading stream tags: " << ex.what());
+    return StatusCode::FAILURE;
+  }
+  catch (...) {
+    ATH_MSG_ERROR("Unknown exception caught when reading stream tags");
+    return StatusCode::FAILURE;
+  }
+  resultToFill.setStreamTags(streamTags);
+  ATH_MSG_DEBUG("Successfully read " << streamTags.size() << " stream tags");
+
+  // ---------------------------------------------------------------------------
+  // Read the HLT bits
+  // ---------------------------------------------------------------------------
+  std::vector<uint32_t> hltBitWords;
+  // try-catch around eformat calls and raw-pointer operations
+  try {
+    const uint32_t nHltInfo = rawEvent->nhlt_info();
+    const uint32_t* hltInfo = rawEvent->hlt_info(); // owned by rawEvent
+    for (uint32_t i=0; i<nHltInfo; ++i)
+      hltBitWords.push_back(hltInfo[i]);
+  }
+  catch (const std::exception& ex) {
+    ATH_MSG_ERROR("std::exception caught when reading HLT bits: " << ex.what());
+    return StatusCode::FAILURE;
+  }
+  catch (...) {
+    ATH_MSG_ERROR("Unknown exception caught when reading HLT bits");
+    return StatusCode::FAILURE;
+  }
+  resultToFill.setHltBits( {hltBitWords.begin(), hltBitWords.end()} );
+  ATH_MSG_DEBUG("Successfully read " << hltBitWords.size() << " HLT bit words");
+
+  return StatusCode::SUCCESS;
+}
+
+// =============================================================================
+StatusCode HLTResultMTByteStreamDecoderTool::decodePayload(const std::vector<const OFFLINE_FRAGMENTS_NAMESPACE::ROBFragment*>& vrobf,
+                                                           HLT::HLTResultMT& resultToFill) const {
+  for (const OFFLINE_FRAGMENTS_NAMESPACE::ROBFragment* robf : vrobf) {
+    eformat::helper::SourceIdentifier sid(robf->rob_source_id());
+    ATH_MSG_DEBUG("Reading ROBFragment " << sid.human());
+    std::vector<uint32_t> data;
+    // try-catch around eformat calls and raw-pointer operations
+    try {
+      const uint32_t nRodData = robf->rod_ndata();
+      const uint32_t* rodData = robf->rod_data(); // owned by robf
+      for (uint32_t i=0; i<nRodData; ++i)
+        data.push_back(rodData[i]);
+    }
+    catch (const std::exception& ex) {
+      ATH_MSG_ERROR("std::exception caught when reading HLT result payload: " << ex.what());
+      return StatusCode::FAILURE;
+    }
+    catch (...) {
+      ATH_MSG_ERROR("Unknown exception caught when reading HLT result payload");
+      return StatusCode::FAILURE;
+    }
+    resultToFill.addSerialisedDataWithCheck(sid.module_id(), data);
+    ATH_MSG_DEBUG("Successfully read " << data.size() << " words of HLT result payload for module ID "
+                  << sid.module_id());
+  }
+  return StatusCode::SUCCESS;
+}

Trigger/TrigSteer/TrigHLTResultByteStream/src/HLTResultMTByteStreamDecoderTool.h

+/*
+  Copyright (C) 2002-2018 CERN for the benefit of the ATLAS collaboration
+*/
+
+#ifndef TRIGHLTRESULTBYTESTREAM_HLTResultMTByteStreamDecoderTool_H
+#define TRIGHLTRESULTBYTESTREAM_HLTResultMTByteStreamDecoderTool_H
+
+// Trigger includes
+#include "TrigSteeringEvent/HLTResultMT.h"
+
+// Athena includes
+#include "AthenaBaseComps/AthAlgTool.h"
+#include "ByteStreamData/RawEvent.h"
+
+/** @class HLTResultMTByteStreamDecoderTool
+ *  @brief Tool filling an HLTResultMT from ByteStream event data
+ **/
+class HLTResultMTByteStreamDecoderTool : virtual public AthAlgTool {
+public:
+  /// Standard constructor
+  HLTResultMTByteStreamDecoderTool(const std::string& type, const std::string& name, const IInterface* parent);
+
+  // ------------------------- AthAlgTool methods ------------------------------
+  virtual StatusCode initialize() override;
+  virtual StatusCode finalize() override;
+
+  // ------------------------- Specific methods of this class ------------------
+  /// Fills the HLTResultMT object from the ByteStream event header
+  StatusCode decodeHeader(const RawEvent* rawEvent, HLT::HLTResultMT& resultToFill) const;
+  /// Fills the HLTResultMT object from vector of ROBFragments
+  StatusCode decodePayload(const std::vector<const OFFLINE_FRAGMENTS_NAMESPACE::ROBFragment*>& vrobf,
+                           HLT::HLTResultMT& resultToFill) const;
+};
+
+#endif // TRIGHLTRESULTBYTESTREAM_HLTResultMTByteStreamDecoderTool_H

Trigger/TrigSteer/TrigHLTResultByteStream/src/components/TrigHLTResultByteStream_entries.cxx

@@ -2,8 +2,12 @@
 #include "TrigHLTResultByteStream/HLTResultByteStreamCnv.h"
 #include "TrigHLTResultByteStream/HLTResultMTByteStreamCnv.h"
 #include "TrigHLTResultByteStream/HLTResultByteStreamTool.h"
+#include "../HLTResultMTByteStreamDecoderAlg.h"
+#include "../HLTResultMTByteStreamDecoderTool.h"
 
 DECLARE_CONVERTER( HLT::HLTResultByteStreamCnv )
 DECLARE_CONVERTER( HLT::HLTResultMTByteStreamCnv )
 DECLARE_COMPONENT( HLT::HLTResultByteStreamTool )
+DECLARE_COMPONENT( HLTResultMTByteStreamDecoderAlg )
+DECLARE_COMPONENT( HLTResultMTByteStreamDecoderTool )
 

Trigger/TrigSteer/TrigOutputHandling/src/TriggerEDMDeserialiserAlg.cxx

@@ -51,7 +51,10 @@ StatusCode TriggerEDMDeserialiserAlg::execute_r(const EventContext& context) con
   auto resultHandle = SG::makeHandle( m_resultKey, context );
   const Payload* dataptr = nullptr;
   // TODO: revise if there are use cases where result may be not available in some events
-  ATH_CHECK( resultHandle->getSerialisedData( m_moduleID, dataptr ) );
+  if ( resultHandle->getSerialisedData( m_moduleID, dataptr ).isFailure() ) {
+    ATH_MSG_DEBUG("No payload available with moduleId " << m_moduleID << " in this event");
+    return StatusCode::SUCCESS;
+  }
   PayloadIterator start = dataptr->begin();
   
   while ( start != dataptr->end() )  {
@@ -68,7 +71,7 @@ StatusCode TriggerEDMDeserialiserAlg::execute_r(const EventContext& context) con
     RootType classDesc = RootType::ByName( transientType+"_v1" ); // TODO remove this dirty hack, needsdiscussion how to find the real type
     size_t usedBytes{ bsize };
     void* obj = m_serializerSvc->deserialize( buff.get(), usedBytes, classDesc );
-    ATH_MSG_DEBUG( "Obtained object " << obj << " which used " << usedBytes << " from available " << bsize  );
+    ATH_MSG_DEBUG( "Obtained object " << obj << " which used " << usedBytes << " bytes from available " << bsize  );
     // for the moment I do not know what do with the raw prt
 
     if ( obj ) {