21#include "Monitoring/Tags.h"
22#include "Monitoring/Metric.h"
23#include "Monitoring/Monitoring.h"
25#include "rapidjson/document.h"
26#include "rapidjson/prettywriter.h"
27#include "rapidjson/filereadstream.h"
30#include "TObjString.h"
34#include <arrow/dataset/file_base.h>
35#include <arrow/dataset/dataset.h>
39#if __has_include(<TJAlienFile.h>)
40#include <TJAlienFile.h>
50using namespace rapidjson;
55 fileNameHolder->fileName = fileName;
57 return fileNameHolder;
61 : mAlienSupport(alienSupport),
65 std::vector<char const*> capabilitiesSpecs = {
66 "O2Framework:RNTupleObjectReadingCapability",
67 "O2Framework:TTreeObjectReadingCapability",
70 std::vector<LoadablePlugin> plugins;
71 for (
auto spec : capabilitiesSpecs) {
73 for (
auto& extra : morePlugins) {
74 plugins.push_back(extra);
78 PluginManager::loadFromPlugin<RootObjectReadingCapability, RootObjectReadingCapabilityPlugin>(plugins, mFactory.
capabilities);
83 LOGP(info,
"DataInputDescriptor");
84 LOGP(info,
" Table name : {}",
tablename);
85 LOGP(info,
" Tree name : {}",
treename);
88 LOGP(info,
" Input files : {}", mfilenames.size());
89 for (
auto fn : mfilenames) {
90 LOGP(info,
" {} {}", fn->fileName, fn->numberOfTimeFrames);
97 return (minputfilesFile.empty() && minputfilesFilePtr) ? (std::string)*minputfilesFilePtr : minputfilesFile;
102 return (mFilenameRegex.empty() && mFilenameRegexPtr) ? (std::string)*mFilenameRegexPtr : mFilenameRegex;
113 if (fn->
fileName.rfind(
"file://", 0) == 0) {
115 }
else if (!mAlienSupport && fn->
fileName.rfind(
"alien://", 0) == 0 && !gGrid) {
116 LOGP(
debug,
"AliEn file requested. Enabling support.");
117 TGrid::Connect(
"alien://");
118 mAlienSupport =
true;
122 mfilenames.emplace_back(fn);
136 if (!
origin.starts_with(
"AOD")) {
141 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
143 if (rootFS->GetFile()->GetName() ==
filename) {
149 mCurrentFilesystem = std::make_shared<TFileFileSystem>(TFile::Open(
filename.c_str()), 50 * 1024 * 1024, mFactory);
150 if (!mCurrentFilesystem.get()) {
151 throw std::runtime_error(fmt::format(
"Couldn't open file \"{}\"!",
filename));
153 rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
157 mParentFileMap = (TMap*)rootFS->GetFile()->Get(
"parentFiles");
160 if (
pos == std::string::npos) {
161 throw std::runtime_error(fmt::format(
"Invalid syntax in aod-parent-base-path-replacement: \"{}\"", mContext.
parentFileReplacement.c_str()));
166 auto it = mParentFileMap->MakeIterator();
167 while (
auto obj = it->Next()) {
168 auto objString = (TObjString*)mParentFileMap->GetValue(obj);
169 objString->String().ReplaceAll(from.c_str(), to.c_str());
175 if (mfilenames[
counter]->numberOfTimeFrames <= 0) {
176 const std::regex TFRegex = std::regex(
"/?DF_([0-9]+)(|-.*)$");
177 TList* keyList = rootFS->GetFile()->GetListOfKeys();
178 std::vector<std::string> finalList;
184 std::unordered_set<size_t> seen;
185 for (
auto key : *keyList) {
186 std::smatch matchResult;
187 std::string keyName = ((TObjString*)
key)->GetString().Data();
188 bool match = std::regex_match(keyName, matchResult, TFRegex);
190 auto folderNumber = std::stoul(matchResult[1].
str());
191 if (seen.find(folderNumber) == seen.end()) {
192 seen.insert(folderNumber);
193 mfilenames[
counter]->listOfTimeFrameNumbers.emplace_back(folderNumber);
198 if (mParentFileMap !=
nullptr) {
200 std::sort(mfilenames[
counter]->listOfTimeFrameNumbers.begin(), mfilenames[
counter]->listOfTimeFrameNumbers.end(),
201 [
this](
long const& l1,
long const& l2) ->
bool {
202 auto p1 = (TObjString*)this->mParentFileMap->GetValue((
"DF_" + std::to_string(l1)).c_str());
203 auto p2 = (TObjString*)this->mParentFileMap->GetValue((
"DF_" + std::to_string(l2)).c_str());
204 return p1->GetString().CompareTo(p2->GetString()) < 0;
207 std::sort(mfilenames[
counter]->listOfTimeFrameNumbers.begin(), mfilenames[
counter]->listOfTimeFrameNumbers.end());
210 mfilenames[
counter]->alreadyRead.resize(mfilenames[
counter]->alreadyRead.size() + mfilenames[
counter]->listOfTimeFrameNumbers.size(),
false);
211 mfilenames[
counter]->numberOfTimeFrames = mfilenames[
counter]->listOfTimeFrameNumbers.size();
215 mCurrentFileStartedAt = uv_hrtime();
230 if (mfilenames[
counter]->numberOfTimeFrames > 0 && numTF >= mfilenames[
counter]->numberOfTimeFrames) {
234 return (mfilenames[
counter]->listOfTimeFrameNumbers)[numTF];
245 if (mfilenames[
counter]->numberOfTimeFrames > 0 && numTF >= mfilenames[
counter]->numberOfTimeFrames) {
249 mfilenames[
counter]->alreadyRead[numTF] =
true;
251 return {fmt::format(
"DF_{}", mfilenames[
counter]->listOfTimeFrameNumbers[numTF]), mCurrentFilesystem};
256 if (!mParentFileMap) {
260 auto folderName = fmt::format(
"DF_{}", mfilenames[
counter]->listOfTimeFrameNumbers[numTF]);
261 auto parentFileName = (TObjString*)mParentFileMap->GetValue(folderName.c_str());
263 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
264 if (!parentFileName) {
265 throw std::runtime_error(fmt::format(R
"(parent file map exists but does not contain the current DF "{}" in file "{}")", folderName.c_str(), rootFS->GetFile()->GetName()));
271 auto parentRootFS = std::dynamic_pointer_cast<TFileFileSystem>(mParentFile->mCurrentFilesystem);
272 if (parentFileName->GetString().CompareTo(parentRootFS->GetFile()->GetName()) == 0) {
277 mParentFile =
nullptr;
282 throw std::runtime_error(fmt::format(R
"(while looking for tree "{}", the parent file was requested but we are already at level {} of maximal allowed level {} for DF "{}" in file "{}")", treename.c_str(), mLevel, mContext.allowedParentLevel, folderName.c_str(),
283 rootFS->GetFile()->GetName()));
286 LOGP(info, "Opening parent file {} for DF {}", parentFileName->GetString().Data(), folderName.c_str());
288 mParentFile->mdefaultFilenamesPtr =
new std::vector<FileNameHolder*>;
289 mParentFile->mdefaultFilenamesPtr->emplace_back(
makeFileNameHolder(parentFileName->GetString().Data()));
297 return mfilenames.at(
counter)->numberOfTimeFrames;
303 return std::count(
list.begin(),
list.end(),
true);
308 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
309 auto f =
dynamic_cast<TFile*
>(rootFS->GetFile());
310 std::string monitoringInfo(fmt::format(
"lfn={},size={}",
f->GetName(),
f->GetSize()));
311#if __has_include(<TJAlienFile.h>)
312 auto alienFile =
dynamic_cast<TJAlienFile*
>(
f);
314 monitoringInfo += fmt::format(
",se={},open_time={:.1f}", alienFile->GetSE(), alienFile->GetElapsed());
318 mContext.
monitoring->send(o2::monitoring::Metric{monitoringInfo,
"aod-file-open-info"}.addTag(o2::monitoring::tags::Key::Subsystem, o2::monitoring::tags::Value::DPL));
320 LOGP(info,
"Opening file: {}", monitoringInfo);
325 int64_t wait_time = (int64_t)uv_hrtime() - (int64_t)mCurrentFileStartedAt - (int64_t)mIOTime;
329 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
330 auto f =
dynamic_cast<TFile*
>(rootFS->GetFile());
331 std::string monitoringInfo(fmt::format(
"lfn={},size={},total_df={},read_df={},read_bytes={},read_calls={},io_time={:.1f},wait_time={:.1f},level={}",
f->GetName(),
333 ((
float)mIOTime / 1e9), ((
float)wait_time / 1e9), mLevel));
334#if __has_include(<TJAlienFile.h>)
335 auto alienFile =
dynamic_cast<TJAlienFile*
>(
f);
337 monitoringInfo += fmt::format(
",se={},open_time={:.1f}", alienFile->GetSE(), alienFile->GetElapsed());
341 mContext.
monitoring->send(o2::monitoring::Metric{monitoringInfo,
"aod-file-read-info"}.addTag(o2::monitoring::tags::Key::Subsystem, o2::monitoring::tags::Value::DPL));
343 LOGP(info,
"Read info: {}", monitoringInfo);
348 if (mCurrentFilesystem.get()) {
352 mParentFile =
nullptr;
355 delete mParentFileMap;
356 mParentFileMap =
nullptr;
359 mCurrentFilesystem.reset();
371 if (!fileName.empty()) {
374 std::ifstream filelist(fileName);
375 if (!filelist.is_open()) {
376 throw std::runtime_error(fmt::format(R
"(Couldn't open file "{}")", fileName));
378 while (std::getline(filelist, fileName)) {
380 fileName.erase(std::remove_if(fileName.begin(), fileName.end(), ::isspace), fileName.end());
392 if (mdefaultFilenamesPtr) {
393 for (
auto fileNameHolder : *mdefaultFilenamesPtr) {
407 auto dfList = mfilenames[
file]->listOfTimeFrameNumbers;
408 auto it = std::find_if(dfList.begin(), dfList.end(), [dfName](
size_t i) { return fmt::format(
"DF_{}", i) == dfName; });
409 if (it == dfList.end()) {
412 return it - dfList.begin();
427 void (*dump_)(
const char*);
428 if (
void* sym = dlsym(
nullptr,
"igprof_dump_now")) {
429 dump_ = __extension__(
void (*)(
const char*)) sym;
431 std::string
filename = fmt::format(
"reader-memory-dump-{}.gz", uv_hrtime());
455 if (!folder.filesystem()) {
460 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(folder.filesystem());
464 throw std::runtime_error(fmt::format(R
"(Not a TFile filesystem!)"));
467 std::shared_ptr<arrow::dataset::FileFormat>
format;
470 auto fullpath = arrow::dataset::FileSource{folder.path() +
"/" +
treename, folder.filesystem()};
473 auto objectPath = capability.lfn2objectPath(fullpath.path());
474 void* handle = capability.getHandle(rootFS, objectPath);
476 format = capability.factory().format();
477 creator = capability.factory().deferredOutputStreamer;
486 LOGP(
debug,
"Could not find tree {}. Trying in parent file.", fullpath.path());
488 if (parentFile !=
nullptr) {
489 int parentNumTF = parentFile->findDFNumber(0, folder.path());
490 if (parentNumTF == -1) {
491 auto parentRootFS = std::dynamic_pointer_cast<TFileFileSystem>(parentFile->mCurrentFilesystem);
492 throw std::runtime_error(fmt::format(R
"(DF {} listed in parent file map but not found in the corresponding file "{}")", folder.path(), parentRootFS->GetFile()->GetName()));
495 return parentFile->readTree(outputs, dh, 0, parentNumTF,
treename, totalSizeCompressed, totalSizeUncompressed);
497 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
498 throw std::runtime_error(fmt::format(R
"(Couldn't get TTree "{}" from "{}". Please check https://aliceo2group.github.io/analysis-framework/docs/troubleshooting/#tree-not-found for more information.)", fullpath.path(), rootFS->GetFile()->GetName()));
501 auto schemaOpt =
format->Inspect(fullpath);
502 auto physicalSchema = schemaOpt;
503 std::vector<std::shared_ptr<arrow::Field>>
fields;
504 for (
auto& original : (*schemaOpt)->fields()) {
505 if (original->name().ends_with(
"_size")) {
508 fields.push_back(original);
510 auto datasetSchema = std::make_shared<arrow::Schema>(
fields);
512 auto fragment =
format->MakeFragment(fullpath, {}, *physicalSchema);
524 f2b->fill(datasetSchema,
format);
532 if (inputFiles.size() == 1 && !inputFiles[0].empty() && inputFiles[0][0] ==
'@') {
535 for (
auto inputFile : inputFiles) {
545 for (
auto fn : mdefaultInputFiles) {
548 mdefaultInputFiles.clear();
549 mdefaultDataInputDescriptor =
nullptr;
551 for (
auto fn : mdataInputDescriptors) {
554 mdataInputDescriptors.clear();
559 mdataInputDescriptors.clear();
560 mdefaultInputFiles.clear();
561 mFilenameRegex = std::string(
"");
566 if (mdefaultDataInputDescriptor) {
567 delete mdefaultDataInputDescriptor;
574 mdefaultDataInputDescriptor->
tablename =
"any";
575 mdefaultDataInputDescriptor->
treename =
"any";
584 FILE*
f = fopen(fnjson.c_str(),
"r");
586 LOGP(error,
"Could not open JSON file \"{}\"!", fnjson);
591 char readBuffer[65536];
592 FileReadStream inputStream(
f, readBuffer,
sizeof(readBuffer));
596 jsonDoc.ParseStream(inputStream);
597 auto status = readJsonDocument(&jsonDoc);
605bool DataInputDirector::readJsonDocument(Document* jsonDoc)
608 std::string fileName(
"");
609 const char* itemName;
612 if (jsonDoc->HasParseError()) {
613 LOGP(error,
"Check the JSON document! There is a problem with the format!");
618 itemName =
"InputDirector";
619 const Value& didirItem = (*jsonDoc)[itemName];
620 if (!didirItem.IsObject()) {
621 LOGP(info,
"No \"{}\" object found in the JSON document!", itemName);
626 itemName =
"debugmode";
627 if (didirItem.HasMember(itemName)) {
628 if (didirItem[itemName].IsBool()) {
629 mDebugMode = (didirItem[itemName].GetBool());
631 LOGP(error,
"Check the JSON document! Item \"{}\" must be a boolean!", itemName);
641 PrettyWriter<StringBuffer> writer(
buffer);
642 didirItem.Accept(writer);
643 LOGP(info,
"InputDirector object: {}", std::string(
buffer.GetString()));
646 itemName =
"fileregex";
647 if (didirItem.HasMember(itemName)) {
648 if (didirItem[itemName].IsString()) {
651 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
656 itemName =
"resfiles";
657 if (didirItem.HasMember(itemName)) {
658 if (didirItem[itemName].IsString()) {
659 fileName = didirItem[itemName].GetString();
660 if (fileName.size() && fileName[0] ==
'@') {
661 fileName.erase(0, 1);
667 }
else if (didirItem[itemName].IsArray()) {
669 auto fns = didirItem[itemName].GetArray();
670 for (
auto& fn : fns) {
674 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string or an array!", itemName);
679 itemName =
"InputDescriptors";
680 if (didirItem.HasMember(itemName)) {
681 if (!didirItem[itemName].IsArray()) {
682 LOGP(error,
"Check the JSON document! Item \"{}\" must be an array!", itemName);
687 for (
auto& didescItem : didirItem[itemName].GetArray()) {
688 if (!didescItem.IsObject()) {
689 LOGP(error,
"Check the JSON document! \"{}\" must be objects!", itemName);
693 auto didesc =
new DataInputDescriptor(mAlienSupport, 0, mContext);
694 didesc->setDefaultInputfiles(&mdefaultInputFiles);
697 if (didescItem.HasMember(itemName)) {
698 if (didescItem[itemName].IsString()) {
699 didesc->tablename = didescItem[itemName].GetString();
702 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
706 LOGP(error,
"Check the JSON document! Item \"{}\" is missing!", itemName);
710 itemName =
"treename";
711 if (didescItem.HasMember(itemName)) {
712 if (didescItem[itemName].IsString()) {
713 didesc->treename = didescItem[itemName].GetString();
715 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
720 didesc->treename =
m[2];
723 itemName =
"fileregex";
724 if (didescItem.HasMember(itemName)) {
725 if (didescItem[itemName].IsString()) {
726 if (didesc->getNumberInputfiles() == 0) {
727 didesc->setFilenamesRegex(didescItem[itemName].GetString());
730 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
734 if (didesc->getNumberInputfiles() == 0) {
735 didesc->setFilenamesRegex(mFilenameRegexPtr);
739 itemName =
"resfiles";
740 if (didescItem.HasMember(itemName)) {
741 if (didescItem[itemName].IsString()) {
742 fileName = didescItem[itemName].GetString();
743 if (fileName.size() && fileName[0] ==
'@') {
744 didesc->setInputfilesFile(fileName.erase(0, 1));
746 if (didesc->getFilenamesRegexString().empty() ||
747 std::regex_match(fileName, didesc->getFilenamesRegex())) {
751 }
else if (didescItem[itemName].IsArray()) {
752 auto fns = didescItem[itemName].GetArray();
753 for (
auto& fn : fns) {
754 if (didesc->getFilenamesRegexString().empty() ||
755 std::regex_match(fn.GetString(), didesc->getFilenamesRegex())) {
760 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string or an array!", itemName);
764 didesc->setInputfilesFile(minputfilesFilePtr);
768 if (didesc->fillInputfiles() > 0) {
769 mdataInputDescriptors.emplace_back(didesc);
772 LOGP(info,
"This DataInputDescriptor is ignored because its file list is empty!");
774 mAlienSupport &= didesc->isAlienSupportOn();
802 for (
auto didesc : mdataInputDescriptors) {
803 if (didesc->matcher->match(dh, context)) {
817 didesc = mdefaultDataInputDescriptor;
829 didesc = mdefaultDataInputDescriptor;
840 didesc = mdefaultDataInputDescriptor;
849 std::string treename;
854 treename = didesc->treename;
859 didesc = mdefaultDataInputDescriptor;
864 auto result = didesc->readTree(outputs, dh,
counter, numTF, treename, totalSizeCompressed, totalSizeUncompressed);
871 for (
auto didesc : mdataInputDescriptors) {
872 didesc->closeInputFile();
876bool DataInputDirector::isValid()
880 for (
auto didesc : mdataInputDescriptors) {
881 status &= didesc->getNumberInputfiles() == numberFiles;
890 for (
auto didesc : mdataInputDescriptors) {
891 status &= (didesc->getNumberInputfiles() <=
counter);
899 LOGP(info,
"DataInputDirector");
900 LOGP(info,
" Default input files file : {}", minputfilesFile);
901 LOGP(info,
" Default file name regex : {}", mFilenameRegex);
902 LOGP(info,
" Default file names : {}", mdefaultInputFiles.size());
903 for (
auto const& fn : mdefaultInputFiles) {
904 LOGP(info,
" {} {}", fn->fileName, fn->numberOfTimeFrames);
906 LOGP(info,
" Default DataInputDescriptor:");
907 mdefaultDataInputDescriptor->
printOut();
909 for (
auto const& didesc : mdataInputDescriptors) {
header::DataOrigin origin
o2::monitoring::tags::Value Value
std::vector< std::shared_ptr< arrow::Field > > fields
#define O2_DECLARE_DYNAMIC_LOG(name)
#define O2_SIGNPOST_ACTION(log, callback)
decltype(auto) make(const Output &spec, Args... args)
std::function< std::shared_ptr< arrow::io::OutputStream >(std::shared_ptr< arrow::dataset::FileFragment >, const std::shared_ptr< arrow::ResizableBuffer > &buffer)> StreamerCreator
void setLabel(const char *label)
bool match(const std::vector< std::string > &queries, const char *pattern)
typedef void(APIENTRYP PFNGLCULLFACEPROC)(GLenum mode)
GLint GLint GLsizei GLint GLenum format
std::string getTreeName(header::DataHeader dh)
Defining PrimaryVertex explicitly as messageable.
FileNameHolder * makeFileNameHolder(std::string fileName)
CalculateDelta(uint64_t &target)
static std::unique_ptr< data_matcher::DataDescriptorMatcher > buildNode(std::string const &nodeString)
static std::vector< std::string > getTokens(std::string const &nodeString)
std::string parentFileReplacement
o2::monitoring::Monitoring * monitoring
static std::vector< LoadablePlugin > parsePluginSpecString(char const *str)
Parse a comma separated list of <library>:<plugin-name> plugin declarations.
std::vector< RootObjectReadingCapability > capabilities