ctf-reader-workflow.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#include <string>
#include <vector>
#include "Framework/Logger.h"
#include "Framework/ControlService.h"
#include "Framework/ConfigParamRegistry.h"
#include "Framework/ChannelSpecHelpers.h"
#include "Framework/InputSpec.h"
#include "CommonUtils/NameConf.h"
#include "CTFWorkflow/CTFReaderSpec.h"
#include "DetectorsCommonDataFormats/DetID.h"
#include "CommonUtils/ConfigurableParam.h"
#include "Algorithm/RangeTokenizer.h"
#include "DetectorsBase/DPLWorkflowUtils.h"
 
// Specific detectors specs
#include "ITSMFTWorkflow/EntropyDecoderSpec.h"
#include "TPCWorkflow/EntropyDecoderSpec.h"
#include "TRDWorkflow/EntropyDecoderSpec.h"
#include "HMPIDWorkflow/EntropyDecoderSpec.h"
#include "FT0Workflow/EntropyDecoderSpec.h"
#include "FV0Workflow/EntropyDecoderSpec.h"
#include "FDDWorkflow/EntropyDecoderSpec.h"
#include "TOFWorkflowUtils/EntropyDecoderSpec.h"
#include "MIDWorkflow/EntropyDecoderSpec.h"
#include "MCHCTF/EntropyDecoderSpec.h"
#include "EMCALWorkflow/EntropyDecoderSpec.h"
#include "PHOSWorkflow/EntropyDecoderSpec.h"
#include "CPVWorkflow/EntropyDecoderSpec.h"
#include "ZDCWorkflow/EntropyDecoderSpec.h"
#include "CTPWorkflow/EntropyDecoderSpec.h"
#ifdef WITH_OPENMP
#include <omp.h>
#endif
 
using namespace o2::framework;
using DetID = o2::detectors::DetID;
 
// we need to add workflow options before including Framework/runDataProcessing
void customize(std::vector<o2::framework::ConfigParamSpec>& workflowOptions)
{
  // option allowing to set parameters
  std::vector<o2::framework::ConfigParamSpec> options;
  options.push_back(ConfigParamSpec{"ctf-input", VariantType::String, "none", {"comma-separated list CTF input files"}});
  options.push_back(ConfigParamSpec{"onlyDet", VariantType::String, std::string{DetID::ALL}, {"comma-separated list of detectors to accept. Overrides skipDet"}});
  options.push_back(ConfigParamSpec{"skipDet", VariantType::String, std::string{DetID::NONE}, {"comma-separate list of detectors to skip"}});
  options.push_back(ConfigParamSpec{"loop", VariantType::Int, 0, {"loop N times (infinite for N<0)"}});
  options.push_back(ConfigParamSpec{"delay", VariantType::Float, 0.f, {"delay in seconds between consecutive TFs sending"}});
  options.push_back(ConfigParamSpec{"copy-cmd", VariantType::String, "alien_cp ?src file://?dst", {"copy command for remote files or no-copy to avoid copying"}}); // Use "XrdSecPROTOCOL=sss,unix xrdcp -N root://eosaliceo2.cern.ch/?src ?dst" for direct EOS access
  options.push_back(ConfigParamSpec{"ctf-file-regex", VariantType::String, ".*o2_ctf_run.+\\.root$", {"regex string to identify CTF files"}});
  options.push_back(ConfigParamSpec{"remote-regex", VariantType::String, "^(alien://|)/alice/data/.+", {"regex string to identify remote files"}}); // Use "^/eos/aliceo2/.+" for direct EOS access
  options.push_back(ConfigParamSpec{"max-cached-files", VariantType::Int, 3, {"max CTF files queued (copied for remote source)"}});
  options.push_back(ConfigParamSpec{"allow-missing-detectors", VariantType::Bool, false, {"send empty message if detector is missing in the CTF (otherwise throw)"}});
  options.push_back(ConfigParamSpec{"send-diststf-0xccdb", VariantType::Bool, false, {"send explicit FLP/DISTSUBTIMEFRAME/0xccdb output"}});
  options.push_back(ConfigParamSpec{"ctf-reader-verbosity", VariantType::Int, 0, {"verbosity level (0: summary per detector, 1: summary per block"}});
  options.push_back(ConfigParamSpec{"ctf-data-subspec", VariantType::Int, 0, {"subspec to use for decoded CTF messages (use non-0 if CTF writer will be attached downstream)"}});
  options.push_back(ConfigParamSpec{"configKeyValues", VariantType::String, "", {"Semicolon separated key=value strings"}});
  options.push_back(ConfigParamSpec{"ir-frames-files", VariantType::String, "", {"If non empty, inject selected IRFrames from this file"}});
  options.push_back(ConfigParamSpec{"run-time-span-file", VariantType::String, "", {"If non empty, inject selected IRFrames from this text file (run, min/max orbit or unix time)"}});
  options.push_back(ConfigParamSpec{"skip-skimmed-out-tf", VariantType::Bool, false, {"Do not process TFs with empty IR-Frame coverage"}});
  options.push_back(ConfigParamSpec{"invert-irframe-selection", VariantType::Bool, false, {"Select only frames mentioned in ir-frames-file (skip-skimmed-out-tf applied to TF not selected!)"}});
  //
  options.push_back(ConfigParamSpec{"its-digits", VariantType::Bool, false, {"convert ITS clusters to digits"}});
  options.push_back(ConfigParamSpec{"mft-digits", VariantType::Bool, false, {"convert MFT clusters to digits"}});
  //
  options.push_back(ConfigParamSpec{"emcal-decoded-subspec", VariantType::Int, 0, {"subspec to use for decoded EMCAL data"}});
  //
  options.push_back(ConfigParamSpec{"timeframes-shm-limit", VariantType::String, "0", {"Minimum amount of SHM required in order to publish data"}});
  options.push_back(ConfigParamSpec{"metric-feedback-channel-format", VariantType::String, "name=metric-feedback,type=pull,method=connect,address=ipc://{}metric-feedback-{},transport=shmem,rateLogging=0", {"format for the metric-feedback channel for TF rate limiting"}});
  options.push_back(ConfigParamSpec{"combine-devices", VariantType::Bool, false, {"combine multiple DPL devices (entropy decoders)"}});
  std::swap(workflowOptions, options);
}
 
// ------------------------------------------------------------------
#include "Framework/runDataProcessing.h"
 
WorkflowSpec defineDataProcessing(ConfigContext const& configcontext)
{
  o2::conf::ConfigurableParam::updateFromString(configcontext.options().get<std::string>("configKeyValues"));
  o2::ctf::CTFReaderInp ctfInput;
 
  WorkflowSpec specs;
  std::string allowedDetectors = "ITS,TPC,TRD,TOF,PHS,CPV,EMC,HMP,MFT,MCH,MID,ZDC,FT0,FV0,FDD,CTP"; // FIXME: explicit list to avoid problem with upgrade detectors
  auto mskOnly = DetID::getMask(configcontext.options().get<std::string>("onlyDet"));
  auto mskSkip = DetID::getMask(configcontext.options().get<std::string>("skipDet"));
  if (mskOnly.any()) {
    ctfInput.detMask &= mskOnly;
  } else {
    ctfInput.detMask ^= mskSkip;
  }
  ctfInput.detMask &= DetID::getMask(allowedDetectors);
  ctfInput.inpdata = configcontext.options().get<std::string>("ctf-input");
  ctfInput.subspec = (unsigned int)configcontext.options().get<int>("ctf-data-subspec");
  ctfInput.decSSpecEMC = (unsigned int)configcontext.options().get<int>("emcal-decoded-subspec");
  if (ctfInput.inpdata.empty() || ctfInput.inpdata == "none") {
    if (!configcontext.helpOnCommandLine()) {
      throw std::runtime_error("--ctf-input <file,...> is not provided");
    }
    ctfInput.inpdata = "";
  }
 
  ctfInput.maxLoops = configcontext.options().get<int>("loop");
  if (ctfInput.maxLoops < 0) {
    ctfInput.maxLoops = 0x7fffffff;
  }
  ctfInput.delay_us = int32_t(1e6 * configcontext.options().get<float>("delay")); // delay in microseconds
  if (ctfInput.delay_us < 0) {
    ctfInput.delay_us = 0;
  }
 
  ctfInput.maxFileCache = std::max(1, configcontext.options().get<int>("max-cached-files"));
 
  ctfInput.copyCmd = configcontext.options().get<std::string>("copy-cmd");
  ctfInput.tffileRegex = configcontext.options().get<std::string>("ctf-file-regex");
  ctfInput.remoteRegex = configcontext.options().get<std::string>("remote-regex");
  ctfInput.allowMissingDetectors = configcontext.options().get<bool>("allow-missing-detectors");
  ctfInput.sup0xccdb = !configcontext.options().get<bool>("send-diststf-0xccdb");
  ctfInput.minSHM = std::stoul(configcontext.options().get<std::string>("timeframes-shm-limit"));
  ctfInput.fileIRFrames = configcontext.options().get<std::string>("ir-frames-files");
  ctfInput.fileRunTimeSpans = configcontext.options().get<std::string>("run-time-span-file");
  ctfInput.skipSkimmedOutTF = configcontext.options().get<bool>("skip-skimmed-out-tf");
  ctfInput.invertIRFramesSelection = configcontext.options().get<bool>("invert-irframe-selection");
  int verbosity = configcontext.options().get<int>("ctf-reader-verbosity");
 
  int rateLimitingIPCID = std::stoi(configcontext.options().get<std::string>("timeframes-rate-limit-ipcid"));
  std::string chanFmt = configcontext.options().get<std::string>("metric-feedback-channel-format");
  if (rateLimitingIPCID > -1 && !chanFmt.empty()) {
    ctfInput.metricChannel = fmt::format(fmt::runtime(chanFmt), o2::framework::ChannelSpecHelpers::defaultIPCFolder(), rateLimitingIPCID);
  }
  if (!ctfInput.fileRunTimeSpans.empty()) {
    ctfInput.skipSkimmedOutTF = true;
  }
  if (!ctfInput.fileIRFrames.empty() && !ctfInput.fileRunTimeSpans.empty()) {
    LOGP(fatal, "One cannot provide --ir-frames-files and --run-time-span-file options simultaneously");
  }
 
  specs.push_back(o2::ctf::getCTFReaderSpec(ctfInput));
 
  auto pipes = configcontext.options().get<std::string>("pipeline");
  std::unordered_map<std::string, int> plines;
  auto ptokens = o2::utils::Str::tokenize(pipes, ',');
  for (auto& token : ptokens) {
    auto split = token.find(":");
    if (split == std::string::npos) {
      throw std::runtime_error("bad pipeline definition. Syntax <processor>:<pipeline>");
    }
    auto key = token.substr(0, split);
    token.erase(0, split + 1);
    size_t error;
    auto value = std::stoll(token, &error, 10);
    if (token[error] != '\0') {
      throw std::runtime_error("Bad pipeline definition. Expecting integer");
    }
    if (value > 1) {
      plines[key] = value;
    }
  }
 
  std::vector<WorkflowSpec> decSpecsV;
 
  auto addSpecs = [&decSpecsV, &plines](DataProcessorSpec&& s) {
    auto entry = plines.find(s.name);
    size_t mult = (entry == plines.end() || entry->second < 2) ? 1 : entry->second;
    if (mult > decSpecsV.size()) {
      decSpecsV.resize(mult);
    }
    decSpecsV[mult - 1].push_back(s);
  };
 
  // add decoders for all allowed detectors.
  if (ctfInput.detMask[DetID::ITS]) {
    addSpecs(o2::itsmft::getEntropyDecoderSpec(DetID::getDataOrigin(DetID::ITS), verbosity, configcontext.options().get<bool>("its-digits"), ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::MFT]) {
    addSpecs(o2::itsmft::getEntropyDecoderSpec(DetID::getDataOrigin(DetID::MFT), verbosity, configcontext.options().get<bool>("mft-digits"), ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::TPC]) {
    addSpecs(o2::tpc::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::TRD]) {
    addSpecs(o2::trd::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::TOF]) {
    addSpecs(o2::tof::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::FT0]) {
    addSpecs(o2::ft0::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::FV0]) {
    addSpecs(o2::fv0::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::FDD]) {
    addSpecs(o2::fdd::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::MID]) {
    addSpecs(o2::mid::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::MCH]) {
    addSpecs(o2::mch::getEntropyDecoderSpec(verbosity, "mch-entropy-decoder", ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::EMC]) {
    addSpecs(o2::emcal::getEntropyDecoderSpec(verbosity, ctfInput.subspec, ctfInput.decSSpecEMC));
  }
  if (ctfInput.detMask[DetID::PHS]) {
    addSpecs(o2::phos::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::CPV]) {
    addSpecs(o2::cpv::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::ZDC]) {
    addSpecs(o2::zdc::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::HMP]) {
    addSpecs(o2::hmpid::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
  if (ctfInput.detMask[DetID::CTP]) {
    addSpecs(o2::ctp::getEntropyDecoderSpec(verbosity, ctfInput.subspec));
  }
 
  bool combine = configcontext.options().get<bool>("combine-devices");
  if (!combine) {
    for (auto& decSpecs : decSpecsV) {
      for (auto& s : decSpecs) {
        specs.push_back(s);
      }
    }
  } else {
    std::vector<DataProcessorSpec> remaining;
    if (decSpecsV.size() && decSpecsV[0].size()) {
      specs.push_back(specCombiner("EntropyDecoders", decSpecsV[0], remaining)); // processing w/o pipelining
    }
    bool updatePipelines = false;
    for (size_t i = 1; i < decSpecsV.size(); i++) { // add pipelined processes separately, consider combining them to separate groups (need to have modify argument of pipeline option)
      if (decSpecsV[i].size() > 1) {
        specs.push_back(specCombiner(fmt::format("EntropyDecodersP{}", i + 1), decSpecsV[i], remaining)); // processing pipelining multiplicity i+1
        updatePipelines = true;
        pipes += fmt::format(",EntropyDecodersP{}:{}", i + 1, i + 1);
      } else {
        for (auto& s : decSpecsV[i]) {
          specs.push_back(s);
        }
      }
    }
    for (auto& s : remaining) {
      specs.push_back(s);
    }
    if (updatePipelines) {
      configcontext.options().override("pipeline", pipes);
    }
  }
 
  return std::move(specs);
}