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);
138 if (wantedParentLevel == -1 && !
origin.starts_with(
"AOD")) {
143 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
145 if (rootFS->GetFile()->GetName() ==
filename) {
151 TFile* tfile =
nullptr;
152 bool externalFile =
false;
160 if (tfile ==
nullptr) {
161 tfile = TFile::Open(
filename.c_str());
163 mCurrentFilesystem = std::make_shared<TFileFileSystem>(tfile, 50 * 1024 * 1024, mFactory, !externalFile);
164 if (!mCurrentFilesystem.get()) {
165 throw std::runtime_error(fmt::format(
"Couldn't open file \"{}\"!",
filename));
167 rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
171 mParentFileMap = (TMap*)rootFS->GetFile()->Get(
"parentFiles");
174 if (
pos == std::string::npos) {
175 throw std::runtime_error(fmt::format(
"Invalid syntax in aod-parent-base-path-replacement: \"{}\"", mContext.
parentFileReplacement.c_str()));
180 auto it = mParentFileMap->MakeIterator();
181 while (
auto obj = it->Next()) {
182 auto objString = (TObjString*)mParentFileMap->GetValue(obj);
183 objString->String().ReplaceAll(from.c_str(), to.c_str());
189 if (mfilenames[
counter]->numberOfTimeFrames <= 0) {
190 const std::regex TFRegex = std::regex(
"/?DF_([0-9]+)(|-.*)$");
191 TList* keyList = rootFS->GetFile()->GetListOfKeys();
192 std::vector<std::string> finalList;
198 std::unordered_set<size_t> seen;
199 for (
auto key : *keyList) {
200 std::smatch matchResult;
201 std::string keyName = ((TObjString*)
key)->GetString().Data();
202 bool match = std::regex_match(keyName, matchResult, TFRegex);
204 auto folderNumber = std::stoul(matchResult[1].
str());
205 if (seen.find(folderNumber) == seen.end()) {
206 seen.insert(folderNumber);
207 mfilenames[
counter]->listOfTimeFrameNumbers.emplace_back(folderNumber);
212 if (mParentFileMap !=
nullptr) {
214 std::sort(mfilenames[
counter]->listOfTimeFrameNumbers.begin(), mfilenames[
counter]->listOfTimeFrameNumbers.end(),
215 [
this](
long const& l1,
long const& l2) ->
bool {
216 auto p1 = (TObjString*)this->mParentFileMap->GetValue((
"DF_" + std::to_string(l1)).c_str());
217 auto p2 = (TObjString*)this->mParentFileMap->GetValue((
"DF_" + std::to_string(l2)).c_str());
218 return p1->GetString().CompareTo(p2->GetString()) < 0;
221 std::sort(mfilenames[
counter]->listOfTimeFrameNumbers.begin(), mfilenames[
counter]->listOfTimeFrameNumbers.end());
224 mfilenames[
counter]->alreadyRead.resize(mfilenames[
counter]->alreadyRead.size() + mfilenames[
counter]->listOfTimeFrameNumbers.size(),
false);
225 mfilenames[
counter]->numberOfTimeFrames = mfilenames[
counter]->listOfTimeFrameNumbers.size();
229 mCurrentFileStartedAt = uv_hrtime();
244 if (mfilenames[
counter]->numberOfTimeFrames > 0 && numTF >= mfilenames[
counter]->numberOfTimeFrames) {
248 return (mfilenames[
counter]->listOfTimeFrameNumbers)[numTF];
254 return {
nullptr, -1};
256 auto folderName = fmt::format(
"DF_{}", mfilenames[
counter]->listOfTimeFrameNumbers[numTF]);
258 if (parentFile ==
nullptr) {
259 return {
nullptr, -1};
261 return {parentFile, parentFile->findDFNumber(0, folderName)};
267 if (wantedParentLevel != -1 && mLevel < wantedParentLevel) {
269 if (parentFile ==
nullptr || parentNumTF == -1) {
272 return parentFile->getFileFolder(0, parentNumTF, wantedParentLevel, wantedOrigin);
281 if (mfilenames[
counter]->numberOfTimeFrames > 0 && numTF >= mfilenames[
counter]->numberOfTimeFrames) {
285 mfilenames[
counter]->alreadyRead[numTF] =
true;
287 return {fmt::format(
"DF_{}", mfilenames[
counter]->listOfTimeFrameNumbers[numTF]), mCurrentFilesystem};
292 if (!mParentFileMap) {
296 auto folderName = fmt::format(
"DF_{}", mfilenames[
counter]->listOfTimeFrameNumbers[numTF]);
297 auto parentFileName = (TObjString*)mParentFileMap->GetValue(folderName.c_str());
299 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
300 if (!parentFileName) {
301 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()));
307 auto parentRootFS = std::dynamic_pointer_cast<TFileFileSystem>(mParentFile->mCurrentFilesystem);
308 if (parentFileName->GetString().CompareTo(parentRootFS->GetFile()->GetName()) == 0) {
313 mParentFile =
nullptr;
318 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(),
319 rootFS->GetFile()->GetName()));
322 LOGP(info, "Opening parent file {} for DF {}", parentFileName->GetString().Data(), folderName.c_str());
324 mParentFile->mdefaultFilenamesPtr =
new std::vector<FileNameHolder*>;
325 mParentFile->mdefaultFilenamesPtr->emplace_back(
makeFileNameHolder(parentFileName->GetString().Data()));
327 mParentFile->
setFile(0, wantedParentLevel, wantedOrigin);
333 return mfilenames.at(
counter)->numberOfTimeFrames;
338 auto& list = mfilenames.at(
counter)->alreadyRead;
339 return std::count(list.begin(), list.end(),
true);
344 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
345 auto f =
dynamic_cast<TFile*
>(rootFS->GetFile());
346 std::string monitoringInfo(fmt::format(
"lfn={},size={}",
f->GetName(),
f->GetSize()));
347#if __has_include(<TJAlienFile.h>)
348 auto alienFile =
dynamic_cast<TJAlienFile*
>(
f);
350 monitoringInfo += fmt::format(
",se={},open_time={:.1f}", alienFile->GetSE(), alienFile->GetElapsed());
354 mContext.
monitoring->send(o2::monitoring::Metric{monitoringInfo,
"aod-file-open-info"}.addTag(o2::monitoring::tags::Key::Subsystem, o2::monitoring::tags::Value::DPL));
356 LOGP(info,
"Opening file: {}", monitoringInfo);
361 int64_t wait_time = (int64_t)uv_hrtime() - (int64_t)mCurrentFileStartedAt - (int64_t)mIOTime;
365 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
366 auto f =
dynamic_cast<TFile*
>(rootFS->GetFile());
367 std::string monitoringInfo(fmt::format(
"lfn={},size={},total_df={},read_df={},read_bytes={},read_calls={},io_time={:.1f},wait_time={:.1f},level={}",
f->GetName(),
369 ((
float)mIOTime / 1e9), ((
float)wait_time / 1e9), mLevel));
370#if __has_include(<TJAlienFile.h>)
371 auto alienFile =
dynamic_cast<TJAlienFile*
>(
f);
373 monitoringInfo += fmt::format(
",se={},open_time={:.1f}", alienFile->GetSE(), alienFile->GetElapsed());
377 mContext.
monitoring->send(o2::monitoring::Metric{monitoringInfo,
"aod-file-read-info"}.addTag(o2::monitoring::tags::Key::Subsystem, o2::monitoring::tags::Value::DPL));
379 LOGP(info,
"Read info: {}", monitoringInfo);
384 if (mCurrentFilesystem.get()) {
388 mParentFile =
nullptr;
391 delete mParentFileMap;
392 mParentFileMap =
nullptr;
395 mCurrentFilesystem.reset();
407 if (!fileName.empty()) {
410 std::ifstream filelist(fileName);
411 if (!filelist.is_open()) {
412 throw std::runtime_error(fmt::format(R
"(Couldn't open file "{}")", fileName));
414 while (std::getline(filelist, fileName)) {
416 fileName.erase(std::remove_if(fileName.begin(), fileName.end(), ::isspace), fileName.end());
428 if (mdefaultFilenamesPtr) {
429 for (
auto fileNameHolder : *mdefaultFilenamesPtr) {
443 auto dfList = mfilenames[file]->listOfTimeFrameNumbers;
444 auto it = std::find_if(dfList.begin(), dfList.end(), [dfName](
size_t i) { return fmt::format(
"DF_{}", i) == dfName; });
445 if (it == dfList.end()) {
448 return it - dfList.begin();
463 void (*dump_)(
const char*);
464 if (
void* sym = dlsym(
nullptr,
"igprof_dump_now")) {
465 dump_ = __extension__(
void (*)(
const char*)) sym;
467 std::string
filename = fmt::format(
"reader-memory-dump-{}.gz", uv_hrtime());
489 std::string wantedOrigin = dh.
dataOrigin.
as<std::string>();
494 if (wantedLevel != -1 && mLevel < wantedLevel) {
496 if (parentFile ==
nullptr) {
497 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
498 throw std::runtime_error(fmt::format(R
"(No parent file found for "{}" while looking for level {} in "{}")", treename, wantedLevel, rootFS->GetFile()->GetName()));
500 if (parentNumTF == -1) {
501 auto parentRootFS = std::dynamic_pointer_cast<TFileFileSystem>(parentFile->mCurrentFilesystem);
502 throw std::runtime_error(fmt::format(R
"(DF not found in parent file "{}")", parentRootFS->GetFile()->GetName()));
505 return parentFile->readTree(outputs, dh, 0, parentNumTF,
treename, totalSizeCompressed, totalSizeUncompressed);
509 if (!folder.filesystem()) {
514 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(folder.filesystem());
518 throw std::runtime_error(fmt::format(R
"(Not a TFile filesystem!)"));
521 std::shared_ptr<arrow::dataset::FileFormat>
format;
524 auto fullpath = arrow::dataset::FileSource{folder.path() +
"/" +
treename, folder.filesystem()};
527 auto objectPath = capability.lfn2objectPath(fullpath.path());
528 void* handle = capability.getHandle(rootFS, objectPath);
530 format = capability.factory().format();
531 creator = capability.factory().deferredOutputStreamer;
540 LOGP(
debug,
"Could not find tree {}. Trying in parent file.", fullpath.path());
542 if (parentFile !=
nullptr) {
543 int parentNumTF = parentFile->findDFNumber(0, folder.path());
544 if (parentNumTF == -1) {
545 auto parentRootFS = std::dynamic_pointer_cast<TFileFileSystem>(parentFile->mCurrentFilesystem);
546 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()));
549 return parentFile->readTree(outputs, dh, 0, parentNumTF,
treename, totalSizeCompressed, totalSizeUncompressed);
551 auto rootFS = std::dynamic_pointer_cast<TFileFileSystem>(mCurrentFilesystem);
552 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()));
555 auto schemaOpt =
format->Inspect(fullpath);
556 auto physicalSchema = schemaOpt;
557 std::vector<std::shared_ptr<arrow::Field>>
fields;
558 for (
auto& original : (*schemaOpt)->fields()) {
559 if (original->name().ends_with(
"_size")) {
562 fields.push_back(original);
564 auto datasetSchema = std::make_shared<arrow::Schema>(
fields);
566 auto fragment =
format->MakeFragment(fullpath, {}, *physicalSchema);
578 f2b->fill(datasetSchema,
format);
586 if (inputFiles.size() == 1 && !inputFiles[0].empty() && inputFiles[0][0] ==
'@') {
589 for (
auto inputFile : inputFiles) {
599 for (
auto fn : mdefaultInputFiles) {
602 mdefaultInputFiles.clear();
603 mdefaultDataInputDescriptor =
nullptr;
605 for (
auto fn : mdataInputDescriptors) {
608 mdataInputDescriptors.clear();
613 mdataInputDescriptors.clear();
614 mdefaultInputFiles.clear();
615 mFilenameRegex = std::string(
"");
620 if (mdefaultDataInputDescriptor) {
621 delete mdefaultDataInputDescriptor;
628 mdefaultDataInputDescriptor->
tablename =
"any";
629 mdefaultDataInputDescriptor->
treename =
"any";
638 FILE*
f = fopen(fnjson.c_str(),
"r");
640 LOGP(error,
"Could not open JSON file \"{}\"!", fnjson);
645 char readBuffer[65536];
646 FileReadStream inputStream(
f, readBuffer,
sizeof(readBuffer));
650 jsonDoc.ParseStream(inputStream);
651 auto status = readJsonDocument(&jsonDoc);
659bool DataInputDirector::readJsonDocument(Document* jsonDoc)
662 std::string fileName(
"");
663 const char* itemName;
666 if (jsonDoc->HasParseError()) {
667 LOGP(error,
"Check the JSON document! There is a problem with the format!");
672 itemName =
"InputDirector";
673 const Value& didirItem = (*jsonDoc)[itemName];
674 if (!didirItem.IsObject()) {
675 LOGP(info,
"No \"{}\" object found in the JSON document!", itemName);
680 itemName =
"debugmode";
681 if (didirItem.HasMember(itemName)) {
682 if (didirItem[itemName].IsBool()) {
683 mDebugMode = (didirItem[itemName].GetBool());
685 LOGP(error,
"Check the JSON document! Item \"{}\" must be a boolean!", itemName);
695 PrettyWriter<StringBuffer> writer(
buffer);
696 didirItem.Accept(writer);
697 LOGP(info,
"InputDirector object: {}", std::string(
buffer.GetString()));
700 itemName =
"fileregex";
701 if (didirItem.HasMember(itemName)) {
702 if (didirItem[itemName].IsString()) {
705 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
710 itemName =
"resfiles";
711 if (didirItem.HasMember(itemName)) {
712 if (didirItem[itemName].IsString()) {
713 fileName = didirItem[itemName].GetString();
714 if (fileName.size() && fileName[0] ==
'@') {
715 fileName.erase(0, 1);
721 }
else if (didirItem[itemName].IsArray()) {
723 auto fns = didirItem[itemName].GetArray();
724 for (
auto& fn : fns) {
728 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string or an array!", itemName);
733 itemName =
"InputDescriptors";
734 if (didirItem.HasMember(itemName)) {
735 if (!didirItem[itemName].IsArray()) {
736 LOGP(error,
"Check the JSON document! Item \"{}\" must be an array!", itemName);
741 for (
auto& didescItem : didirItem[itemName].GetArray()) {
742 if (!didescItem.IsObject()) {
743 LOGP(error,
"Check the JSON document! \"{}\" must be objects!", itemName);
747 auto didesc =
new DataInputDescriptor(mAlienSupport, 0, mContext);
748 didesc->setDefaultInputfiles(&mdefaultInputFiles);
751 if (didescItem.HasMember(itemName)) {
752 if (didescItem[itemName].IsString()) {
753 didesc->tablename = didescItem[itemName].GetString();
756 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
760 LOGP(error,
"Check the JSON document! Item \"{}\" is missing!", itemName);
764 itemName =
"treename";
765 if (didescItem.HasMember(itemName)) {
766 if (didescItem[itemName].IsString()) {
767 didesc->treename = didescItem[itemName].GetString();
769 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
774 didesc->treename =
m[2];
777 itemName =
"fileregex";
778 if (didescItem.HasMember(itemName)) {
779 if (didescItem[itemName].IsString()) {
780 if (didesc->getNumberInputfiles() == 0) {
781 didesc->setFilenamesRegex(didescItem[itemName].GetString());
784 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string!", itemName);
788 if (didesc->getNumberInputfiles() == 0) {
789 didesc->setFilenamesRegex(mFilenameRegexPtr);
793 itemName =
"resfiles";
794 if (didescItem.HasMember(itemName)) {
795 if (didescItem[itemName].IsString()) {
796 fileName = didescItem[itemName].GetString();
797 if (fileName.size() && fileName[0] ==
'@') {
798 didesc->setInputfilesFile(fileName.erase(0, 1));
800 if (didesc->getFilenamesRegexString().empty() ||
801 std::regex_match(fileName, didesc->getFilenamesRegex())) {
805 }
else if (didescItem[itemName].IsArray()) {
806 auto fns = didescItem[itemName].GetArray();
807 for (
auto& fn : fns) {
808 if (didesc->getFilenamesRegexString().empty() ||
809 std::regex_match(fn.GetString(), didesc->getFilenamesRegex())) {
814 LOGP(error,
"Check the JSON document! Item \"{}\" must be a string or an array!", itemName);
818 didesc->setInputfilesFile(minputfilesFilePtr);
822 if (didesc->fillInputfiles() > 0) {
823 mdataInputDescriptors.emplace_back(didesc);
826 LOGP(info,
"This DataInputDescriptor is ignored because its file list is empty!");
828 mAlienSupport &= didesc->isAlienSupportOn();
856 for (
auto didesc : mdataInputDescriptors) {
857 if (didesc->matcher->match(dh, context)) {
871 didesc = mdefaultDataInputDescriptor;
876 return didesc->getFileFolder(
counter, numTF, wantedLevel,
origin);
884 didesc = mdefaultDataInputDescriptor;
895 didesc = mdefaultDataInputDescriptor;
900 return didesc->getTimeFrameNumber(
counter, numTF, wantedLevel,
origin);
905 std::string treename;
910 treename = didesc->treename;
915 didesc = mdefaultDataInputDescriptor;
920 auto result = didesc->readTree(outputs, dh,
counter, numTF, treename, totalSizeCompressed, totalSizeUncompressed);
927 for (
auto didesc : mdataInputDescriptors) {
928 didesc->closeInputFile();
932bool DataInputDirector::isValid()
936 for (
auto didesc : mdataInputDescriptors) {
937 status &= didesc->getNumberInputfiles() == numberFiles;
946 for (
auto didesc : mdataInputDescriptors) {
947 status &= (didesc->getNumberInputfiles() <=
counter);
955 LOGP(info,
"DataInputDirector");
956 LOGP(info,
" Default input files file : {}", minputfilesFile);
957 LOGP(info,
" Default file name regex : {}", mFilenameRegex);
958 LOGP(info,
" Default file names : {}", mdefaultInputFiles.size());
959 for (
auto const& fn : mdefaultInputFiles) {
960 LOGP(info,
" {} {}", fn->fileName, fn->numberOfTimeFrames);
962 LOGP(info,
" Default DataInputDescriptor:");
963 mdefaultDataInputDescriptor->
printOut();
965 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)
GLuint const GLchar * name
typedef void(APIENTRYP PFNGLCULLFACEPROC)(GLenum mode)
GLint GLint GLsizei GLint GLenum format
std::string getTreeName(header::DataHeader dh)
Defining ITS Vertex 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
std::vector< std::pair< std::string, TFile * > > openFiles
int levelForOrigin(std::string_view origin) const
static std::vector< LoadablePlugin > parsePluginSpecString(char const *str)
Parse a comma separated list of <library>:<plugin-name> plugin declarations.
std::vector< RootObjectReadingCapability > capabilities