CbmMcbm2018RawConverterGdpb.cxx

Go to the documentation of this file.

/* Copyright (C) 2019 Facility for Antiproton and Ion Research in Europe, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Pierre-Alain Loizeau [committer] */
 
// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                     CbmMcbm2018RawConverterGdpb                   -----
// -----               Created 01.07.2019 by P.-A. Loizeau                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------
 
#include "CbmMcbm2018RawConverterGdpb.h"
 
#include "CbmFormatMsHeaderPrintout.h"
#include "CbmMcbm2018TofPar.h"
#include "CbmStsDigi.h"
 
#include "FairParGenericSet.h"
#include "FairRootManager.h"
#include "FairRun.h"
#include "FairRunOnline.h"
#include "FairRuntimeDb.h"
#include <Logger.h>
 
#include <iomanip>
 
CbmMcbm2018RawConverterGdpb::CbmMcbm2018RawConverterGdpb(UInt_t /*uNbGdpb*/)
  : CbmMcbmUnpack()
  , fParCList(nullptr)
  , fUnpackPar(nullptr)
  , fuNrOfDpbs(0)
  , fDpbIdIndexMap()
  , fvbMaskedComponents()
  , fvMsComponentsList()
  , fuNbCoreMsPerTs(0)
  , fuNbOverMsPerTs(0)
  , fuNbMsLoop(0)
  , fbIgnoreOverlapMs(kFALSE)
  , fdMsSizeInNs(-1.0)
  , fdTsCoreSizeInNs(-1.0)
  , fulCurrentTsIdx(0)
  , fulCurrentMsIdx(0)
  , fdTsStartTime(0.0)
  , fdTsStopTimeCore(0.0)
  , fdMsTime(0.0)
  , fuMsIndex(0)
  , fuCurrentEquipmentId(0)
  , fuCurrDpbId(0)
  , fuCurrDpbIdx(0)
  , fvGdpbMessages()
{
}
 
CbmMcbm2018RawConverterGdpb::~CbmMcbm2018RawConverterGdpb() {}
 
Bool_t CbmMcbm2018RawConverterGdpb::Init()
{
  LOG(info) << "CbmMcbm2018RawConverterGdpb::Init";
  LOG(info) << "Initializing mCBM gDPB 2018 Raw Messages Converter";
 
  return kTRUE;
}
 
void CbmMcbm2018RawConverterGdpb::SetParContainers()
{
  LOG(info) << "Setting parameter containers for " << GetName();
 
  if (nullptr == fParCList) fParCList = new TList();
  fUnpackPar = new CbmMcbm2018TofPar("CbmMcbm2018TofPar");
  fParCList->Add(fUnpackPar);
 
  for (Int_t iparC = 0; iparC < fParCList->GetEntries(); ++iparC) {
    FairParGenericSet* tempObj = (FairParGenericSet*) (fParCList->At(iparC));
    fParCList->Remove(tempObj);
 
    std::string sParamName {tempObj->GetName()};
    FairParGenericSet* newObj =
      dynamic_cast<FairParGenericSet*>(FairRun::Instance()->GetRuntimeDb()->getContainer(sParamName.data()));
 
    if (nullptr == newObj) {
      LOG(error) << "Failed to obtain parameter container " << sParamName << ", for parameter index " << iparC;
      return;
    }  // if( nullptr == newObj )
 
    fParCList->AddAt(newObj, iparC);
    //      delete tempObj;
  }  // for( Int_t iparC = 0; iparC < fParCList->GetEntries(); ++iparC )
}
 
Bool_t CbmMcbm2018RawConverterGdpb::InitContainers()
{
  LOG(info) << "Init parameter containers for " << GetName();
 
  fUnpackPar = dynamic_cast<CbmMcbm2018TofPar*>(FairRun::Instance()->GetRuntimeDb()->getContainer("CbmMcbm2018TofPar"));
  if (nullptr == fUnpackPar) {
    LOG(error) << "Failed to obtain parameter container CbmMcbm2018TofPar";
    return kFALSE;
  }  // if( nullptr == fUnpackPar )
 
  fuNrOfDpbs = fUnpackPar->GetNrOfGdpbs();
  LOG(info) << "Nr. of STS DPBs:       " << fuNrOfDpbs;
 
  fDpbIdIndexMap.clear();
  for (UInt_t uDpb = 0; uDpb < fuNrOfDpbs; ++uDpb) {
    fDpbIdIndexMap[fUnpackPar->GetGdpbId(uDpb)] = uDpb;
    LOG(info) << "Eq. ID for DPB #" << std::setw(2) << uDpb << " = 0x" << std::setw(4) << std::hex
              << fUnpackPar->GetGdpbId(uDpb) << std::dec << " => " << fDpbIdIndexMap[fUnpackPar->GetGdpbId(uDpb)];
  }  // for( UInt_t uDpb = 0; uDpb < fuNrOfDpbs; ++uDpb )
 
  FairRootManager* ioman = FairRootManager::Instance();
  if (NULL == ioman) { LOG(fatal) << "No FairRootManager instance"; }
 
  fvGdpbMessages.resize(fuNrOfDpbs);
  for (UInt_t uDpb = 0; uDpb < fuNrOfDpbs; ++uDpb) {
    fvGdpbMessages[uDpb] = new std::vector<gdpbv100::Message>;
    if (NULL == fvGdpbMessages[uDpb]) {
      LOG(fatal) << "Failed creating the gDPB messages vector ";
    }  // if( NULL == fvGdpbMessages[ uDpb ] )
    ioman->RegisterAny(Form("CbmGdpbMsg%03u", uDpb), fvGdpbMessages[uDpb], kTRUE);
  }  // for( UInt_t uDpb = 0; uDpb < fuNrOfDpbs; ++uDpb )
 
  return kTRUE;
}
 
Bool_t CbmMcbm2018RawConverterGdpb::ReInitContainers()
{
  LOG(info) << "ReInit parameter containers for " << GetName();
 
  return kTRUE;
}
 
void CbmMcbm2018RawConverterGdpb::AddMsComponentToList(size_t component, UShort_t usDetectorId)
{
  for (UInt_t uCompIdx = 0; uCompIdx < fvMsComponentsList.size(); ++uCompIdx)
    if (component == fvMsComponentsList[uCompIdx]) return;
 
  fvMsComponentsList.push_back(component);
 
  LOG(info) << "CbmMcbm2018RawConverterGdpb::AddMsComponentToList => Component " << component << " with detector ID 0x"
            << std::hex << usDetectorId << std::dec << " added to list";
}
 
Bool_t CbmMcbm2018RawConverterGdpb::DoUnpack(const fles::Timeslice& ts, size_t /*component*/)
{
  fulCurrentTsIdx = ts.index();
  fdTsStartTime   = static_cast<Double_t>(ts.descriptor(0, 0).idx);
 
  if (0 == fulCurrentTsIdx) { return kTRUE; }  // if( 0 == fulCurrentTsIdx )
 
  if (-1.0 == fdTsCoreSizeInNs) {
    fuNbCoreMsPerTs  = ts.num_core_microslices();
    fuNbOverMsPerTs  = ts.num_microslices(0) - ts.num_core_microslices();
    fdTsCoreSizeInNs = fdMsSizeInNs * (fuNbCoreMsPerTs);
    fdTsFullSizeInNs = fdMsSizeInNs * (fuNbCoreMsPerTs + fuNbOverMsPerTs);
    LOG(info) << "Timeslice parameters: each TS has " << fuNbCoreMsPerTs << " Core MS and " << fuNbOverMsPerTs
              << " Overlap MS, for a core duration of " << fdTsCoreSizeInNs << " ns and a full duration of "
              << fdTsFullSizeInNs << " ns";
 
    fuNbMsLoop = fuNbCoreMsPerTs;
    if (kFALSE == fbIgnoreOverlapMs) fuNbMsLoop += fuNbOverMsPerTs;
    LOG(info) << "In each TS " << fuNbMsLoop << " MS will be looped over";
  }  // if( -1.0 == fdTsCoreSizeInNs )
 
  fdTsStopTimeCore = fdTsStartTime + fdTsCoreSizeInNs;
 
  for (UInt_t uMsCompIdx = 0; uMsCompIdx < fvMsComponentsList.size(); ++uMsCompIdx) {
    for (fuMsIndex = 0; fuMsIndex < fuNbMsLoop; fuMsIndex++) {
      UInt_t uMsComp = fvMsComponentsList[uMsCompIdx];
 
      auto msDescriptor        = ts.descriptor(uMsComp, fuMsIndex);
      fuCurrentEquipmentId     = msDescriptor.eq_id;
      const uint8_t* msContent = reinterpret_cast<const uint8_t*>(ts.content(uMsComp, fuMsIndex));
 
      uint32_t uSize  = msDescriptor.size;
      fulCurrentMsIdx = msDescriptor.idx;
      //         Double_t dMsTime = (1e-9) * static_cast<double>(fulCurrentMsIdx);
      LOG(debug) << "Microslice: " << fulCurrentMsIdx << " from EqId " << std::hex << fuCurrentEquipmentId << std::dec
                 << " has size: " << uSize;
 
      if (0 == fvbMaskedComponents.size()) fvbMaskedComponents.resize(ts.num_components(), kFALSE);
 
      fuCurrDpbId = static_cast<uint32_t>(fuCurrentEquipmentId & 0xFFFF);
 
      auto it = fDpbIdIndexMap.find(fuCurrDpbId);
      if (it == fDpbIdIndexMap.end()) {
        if (kFALSE == fvbMaskedComponents[uMsCompIdx]) {
          LOG(info) << "-------------------------------------------------------"
                       "--------";
          /*
                LOG(info) << "hi hv eqid flag si sv idx/start        crc      size     offset";
                LOG(info) << Form( "%02x %02x %04x %04x %02x %02x %016llx %08x %08x %016llx",
                                  static_cast<unsigned int>(msDescriptor.hdr_id),
                                  static_cast<unsigned int>(msDescriptor.hdr_ver), msDescriptor.eq_id, msDescriptor.flags,
                                  static_cast<unsigned int>(msDescriptor.sys_id),
                                  static_cast<unsigned int>(msDescriptor.sys_ver), msDescriptor.idx, msDescriptor.crc,
                                  msDescriptor.size, msDescriptor.offset );
*/
          LOG(info) << FormatMsHeaderPrintout(msDescriptor);
          LOG(warning) << "Could not find the gDPB index for AFCK id 0x" << std::hex << fuCurrDpbId << std::dec
                       << " in timeslice " << fulCurrentTsIdx << " in microslice " << fuMsIndex << " component "
                       << uMsCompIdx << "\n"
                       << "If valid this index has to be added in the TOF "
                          "parameter file in the DbpIdArray field";
          fvbMaskedComponents[uMsCompIdx] = kTRUE;
        }  // if( kFALSE == fvbMaskedComponents[ uMsComp ] )
        else
          return kTRUE;
 
 
        return kFALSE;
      }  // if( it == fDpbIdIndexMap.end() )
      else
        fuCurrDpbIdx = fDpbIdIndexMap[fuCurrDpbId];
 
      if (0 != (uSize % kuBytesPerMessage))
        LOG(error) << "The input microslice buffer does NOT "
                   << "contain only complete gDPB messages!";
 
      uint32_t uNbMessages = (uSize - (uSize % kuBytesPerMessage)) / kuBytesPerMessage;
 
      const uint64_t* pInBuff = reinterpret_cast<const uint64_t*>(msContent);
      for (uint32_t uIdx = 0; uIdx < uNbMessages; uIdx++) {
        uint64_t ulData = static_cast<uint64_t>(pInBuff[uIdx]);
 
        if (0 == uIdx) { continue; }  // if( 0 == uIdx )
 
        gdpbv100::Message mess(ulData);
 
        fvGdpbMessages[fuCurrDpbIdx]->push_back(mess);
      }  // for (uint32_t uIdx = 0; uIdx < uNbMessages; uIdx ++)
    }    // for( fuMsIndex = 0; fuMsIndex < uNbMsLoop; fuMsIndex ++ )
  }      // for( UInt_t uMsCompIdx = 0; uMsCompIdx < fvMsComponentsList.size(); ++uMsCompIdx )
 
  if (0 == fulCurrentTsIdx % 10000) LOG(info) << "Processed TS " << fulCurrentTsIdx;
 
  return kTRUE;
}
 
void CbmMcbm2018RawConverterGdpb::Reset()
{
  for (UInt_t uDpb = 0; uDpb < fuNrOfDpbs; ++uDpb)
    fvGdpbMessages[uDpb]->clear();
}
 
void CbmMcbm2018RawConverterGdpb::Finish() {}
 
ClassImp(CbmMcbm2018RawConverterGdpb)