KVExpDB.cpp

//Created by KVClassFactory on Tue Jul 12 11:43:52 2016
//Author: bonnet,,,
 
#include "KVExpDB.h"
#include "KVDBSystem.h"
#include "KVNumberList.h"
#include "TSystem.h"
#include <KVFileReader.h>
#include <KVNucleus.h>
#include <iostream>
#include "KVUnownedList.h"
#include "KVEnv.h"
using namespace std;
 
KVExpDB* gExpDB = nullptr;
 
ClassImp(KVExpDB)
 
 
 
void KVExpDB::init()
{
   //default initialisations
 
   fRuns = AddTable("Runs", "List of available runs");
   fRuns->SetDefaultFormat("Run %d"); // default format for run names
   fSystems = AddTable("Systems", "List of available systems");
   index_multiplier = KVBase::GetDataSetEnv(fDataSet, "DataSet.RunFileIndexMultiplier.raw", 1);
   with_index_multiplier = index_multiplier != 1;
   raw_files_have_index = KVBase::GetDataSetEnv(fDataSet, "DataSet.RunFileName.raw", "").Contains("%I");
}
 
 
 
 
void KVExpDB::read_beams_targets_reactions()
{
   // Read infos on beams, targets, and reactions in files
   //~~~
   //Beams.dat
   //Targets.dat
   //Reactions.dat
   //~~~
   //(default filenames defined by config variables `EXPDB.Beams`, `EXPDB.Targets`, `EXPDB.Reactions`)
   //and set up the corresponding KVTarget and KVDBSystem objects, adding to the database
   //table `Systems`.
   //
   //###Beams
   //Format of file is
   //~~~
   //Beams: B1 B2 [...]
   //
   //B1.Ion: 129Xe      [nature of accelerated ion]
   //B1.Energy: 49.9976 [energy in MeV/A]
   //
   //B2.Ion: [...]
   //~~~
   //
   //###Targets
   //Format of file is
   //~~~
   //Targets: T1 T2 [...]
   //
   //T1.Material: 197Au
   //T1.AreaDensity: 0.66  [in mg/cm2]
   //[T1.Angle: 45]        [angle to beam in degrees, if != 0]
   //[T1.Name: Au-660]     [optional name]
   //
   //T2.Layers: 2
   //T2.Layer.1.Material: 12C
   //T2.Layer.1.AreaDensity: 0.02 [in mg/cm2]
   //T2.Layer.2.Material: 238U
   //T2.Layer.2.AreaDensity: 0.4  [in mg/cm2]
   //T2.TargetNucleus: 238U       [used for reaction kinematics]
   //~~~
   //
   //###Reactions
   //Format of file is
   //~~~
   //Reactions: R1 R2 [...]
   //
   //R1.Name: Alpha source
   //R1.Runlist: 1-12
   //
   //R2.Beam: B1     [beam defined in Beams.dat]
   //R2.Target: T1   [target defined in Targets.dat]
   //R2.Runlist: 13-25
   //~~~
 
   TString beam_fp, target_fp, reaction_fp;
   if (!(FindCalibFile("Beams", beam_fp) && FindCalibFile("Targets", target_fp) && FindCalibFile("Reactions", reaction_fp))) {
      Error("read_beams_targets_reactions",
            "One or more files not found:\nFull path to %s = %s\nFull path to %s = %s\nFull path to %s = %s",
            GetCalibFileName("Beams").Data(), beam_fp.Data(),
            GetCalibFileName("Targets").Data(), target_fp.Data(),
            GetCalibFileName("Reactions").Data(), reaction_fp.Data());
      return;
   }
 
   Info("read_beams_targets_reactions", "Reading files:\nFull path to %s = %s\nFull path to %s = %s\nFull path to %s = %s",
        GetCalibFileName("Beams").Data(), beam_fp.Data(),
        GetCalibFileName("Targets").Data(), target_fp.Data(),
        GetCalibFileName("Reactions").Data(), reaction_fp.Data());
 
   KVEnv beams_v(beam_fp), targets_v(target_fp), reactions_v(reaction_fp);
 
   auto make_target = [&](const KVString & react) -> KVTarget * {
      if (!reactions_v.HasValue(react, "Target"))
         return nullptr;
 
      auto targ = reactions_v.GetValueOf(react, "Target").Default<KVString>();
 
      auto nlayers = targets_v.GetValueOf(targ, "Layers").Default(1);
 
      auto T = new KVTarget;
      if (nlayers > 1)
      {
         for (int ilay = 1; ilay <= nlayers; ++ilay) {
            T->AddLayer(targets_v.GetValueOf(targ, "Layer", ilay, "Material").Default<KVString>(),
                        targets_v.GetValueOf(targ, "Layer", ilay, "AreaDensity").Default(0.));
         }
      }
      else if (targets_v.HasValue(targ, "Material"))
      {
         T->AddLayer(targets_v.GetValueOf(targ, "Material").Default<KVString>(),
                     targets_v.GetValueOf(targ, "AreaDensity").Default(0.));
      }
 
      T->SetAngleToBeam(targets_v.GetValueOf(targ, "Angle").Default(0.));
      if (targets_v.HasValue(targ, "Name"))
         T->SetName(targets_v.GetValueOf(targ, "Name").Default<KVString>());
 
      return T;
   };
   auto get_beam_properties = [&](const KVString & react) {
      auto beam = reactions_v.GetValueOf(react, "Beam").Default<KVString>();
      return KVNucleus(beams_v.GetValueOf(beam, "Ion").Default<KVString>(),
                       beams_v.GetValueOf(beam, "Energy").Default(0.));
   };
   auto get_target_properties = [&](const KVString & react) {
      auto targ = reactions_v.GetValueOf(react, "Target").Default<KVString>();
      if (targets_v.HasValue(targ, "TargetNucleus"))
         return KVNucleus(targets_v.GetValueOf(targ, "TargetNucleus").Default<KVString>());
      return KVNucleus(targets_v.GetValueOf(targ, "Material").Default<KVString>());
   };
   auto get_target_name = [&](const KVString & react) {
      auto targ = reactions_v.GetValueOf(react, "Target").Default<KVString>();
      return targets_v.GetValueOf(targ, "Name").Default<KVString>("");
   };
 
   KVString reactions = reactions_v.GetValue("Reactions", "");
   reactions.Begin(" ");
   while (!reactions.End()) {
      auto R = reactions.Next();
 
      KVDBSystem* sys = nullptr;
 
      if (!reactions_v.HasValue(R, "Beam") && !reactions_v.HasValue(R, "Target"))
         sys = new KVDBSystem(reactions_v.GetValueOf(R, "Name").Default<KVString>());
      else
         sys = new KVDBSystem(get_beam_properties(R), get_target_properties(R), get_target_name(R));
 
      // system with a defined target?
      auto targ = make_target(R);
      if (targ)
         sys->SetTarget(targ);
 
      AddSystem(sys);
 
      sys->SetRuns(reactions_v.GetValueOf(R, "Runlist").Default<KVNumberList>(), index_multiplier);
   }
}
 
 
 
 
KVExpDB::KVExpDB()
   : KVDataBase()
{
   // Default constructor
}
 
 
 
 
 
KVExpDB::KVExpDB(const Char_t* name)
   : KVDataBase(name), fDataSet(name)
{
   // Constructor inherited from KVDataBase
   init();
}
 
 
 
 
 
KVExpDB::KVExpDB(const Char_t* name, const Char_t* title)
   : KVDataBase(name, title), fDataSet(name)
{
   // Constructor inherited from KVDataBase
   init();
}
 
 
 
 
 
KVExpDB::~KVExpDB()
{
   // Destructor
   if (gExpDB == this) gExpDB = nullptr;
}
 
 
 
 
void KVExpDB::LinkRecordToRunRange(KVDBRecord* rec, UInt_t first_run,
                                   UInt_t last_run)
{
   //If the KVDBRecord 'rec' (i.e. set of calibration parameters, reaction system, etc.) is
   //associated to, or valid for, a range of runs, we use this method in order to link the record
   //and the runs. The list of associated runs will be kept with the record, and each of the runs
   //will have a link to the record.
 
   std::map<int, int> one_shot; //make sure each run is only linked once
   for (UInt_t ii = first_run; ii <= last_run; ii++) {
      Int_t rr = ii;
      if (rr > index_multiplier) rr /= index_multiplier;
      if (!one_shot[rr]) LinkRecordToRun(rec, ii);
      ++one_shot[rr];
   }
}
 
 
 
void KVExpDB::LinkRecordToRunRange(KVDBRecord* rec, const KVNumberList& nl)
{
   //If the KVDBRecord 'rec' (i.e. set of calibration parameters, reaction system, etc.) is
   //associated to, or valid for, a range of runs, we use this method in order to link the record
   //and the runs. The list of associated runs will be kept with the record, and each of the runs
   //will have a link to the record.
   nl.Begin();
   std::map<int, int> one_shot; //make sure each run is only linked once
   while (!nl.End()) {
      Int_t rr = nl.Next();
      if (rr > index_multiplier) rr /= index_multiplier;
      if (!one_shot[rr]) LinkRecordToRun(rec, rr);
      ++one_shot[rr];
   }
}
 
 
 
 
void KVExpDB::LinkRecordToRun(KVDBRecord* rec, Int_t rnumber)
{
 
   KVDBRecord* run = 0;
   if ((run = fRuns->GetRecord(rnumber)))
      rec->AddLink("Runs", run);
 
}
 
 
 
 
void KVExpDB::LinkRecordToRunRanges(KVDBRecord* rec, UInt_t rr_number,
                                    UInt_t run_ranges[][2])
{
   //Call LinkRecordToRunRange for a set of run ranges stored in the two-dimensional array
   //in the following way:
   //      run_ranges[0][0] = first run of first run range
   //      run_ranges[0][1] = last run of first run range
   //      run_ranges[1][0] = first run of second run range etc. etc.
   //rr_number is the number of run ranges in the array
 
   for (UInt_t i = 0; i < rr_number; i++) {
      LinkRecordToRunRange(rec, run_ranges[i][0], run_ranges[i][1]);
   }
}
 
 
 
 
void KVExpDB::LinkListToRunRanges(TList* list, UInt_t rr_number,
                                  UInt_t run_ranges[][2])
{
   //Link the records contained in the list to the set of runs (see LinkRecordToRunRanges).
 
   if (!list) {
      Error("LinkListToRunRanges",
            "NULL pointer passed for parameter TList");
      return;
   }
   if (list->GetSize() == 0) {
      Error("LinkListToRunRanges(TList*,UInt_t,UInt_t*)",
            "The list is empty");
      return;
   }
   TIter next(list);
   KVDBRecord* rec;
   while ((rec = (KVDBRecord*) next())) {
      LinkRecordToRunRanges(rec, rr_number, run_ranges);
   }
}
 
 
 
void KVExpDB::LinkListToRunRange(TList* list, const KVNumberList& nl)
{
   //Link the records contained in the list to the set of runs (see LinkRecordToRunRanges).
 
   if (!list) {
      Error("LinkListToRunRange",
            "NULL pointer passed for parameter TList");
      return;
   }
   if (list->GetSize() == 0) {
      Error("LinkListToRunRange(TList*,KVNumberList)",
            "The list is empty");
      return;
   }
   TIter next(list);
   KVDBRecord* rec;
   while ((rec = (KVDBRecord*) next())) {
      LinkRecordToRunRange(rec, nl);
   }
}
 
 
 
 
void KVExpDB::ReadSystemList()
{
   // Set up all informations on the reactions studied during the experiment,
   // store as KVDBSystem objects in the database table `Systems`.
   //
   // For information on necessary input files and formats, see read_beams_targets_reactions()
 
   TString fp;
   if (FindCalibFile("Reactions", fp)) {
      Info("ReadSystemList", "Reading reaction parameters...");
      read_beams_targets_reactions();
   }
   else { // old-style 'Systems.dat' file
      //There are 2 formats for the description of systems:
      //
      //+129Xe + natSn 25 MeV/A           '+' indicates beginning of system description
      //129 54 119 50 0.1 25.0            A,Z of projectile and target, target thickness (mg/cm2), beam energy (MeV/A)
      //Run Range : 614 636               runs associated with system
      //Run Range : 638 647               runs associated with system
      //
      //This is sufficient in the simple case where the experiment has a single
      //layer target oriented perpendicular to the beam. However, for more
      //complicated targets we can specify as follows :
      //
      //+155Gd + 238U 36 MeV/A
      //155 64 238 92 0.1 36.0
      //Target : 3 0.0                    target with 3 layers, angle 0 degrees
      //C 0.02                            1st layer : carbon, 20 g/cm2
      //238U 0.1                          2nd layer : uranium-238, 100 g/cm2
      //C 0.023                           3rd layer : carbon, 23 g/cm2
      //Run Range : 770 804
      //
      //Lines beginning '#' are comments.
      std::ifstream fin;
      if (OpenCalibFile("Systems", fin)) {
         Info("ReadSystemList()", "Reading Systems parameters ...");
 
         TString line;
 
         char next_char = fin.peek();
         while (next_char != '+' && fin.good()) {
            line.ReadLine(fin, kFALSE);
            next_char = fin.peek();
         }
 
         while (fin.good() && !fin.eof() && next_char == '+') {
            KVDBSystem* sys = new KVDBSystem("NEW SYSTEM");
            AddSystem(sys);
            sys->Load(fin, index_multiplier);
            next_char = fin.peek();
         }
         fin.close();
      }
      else {
         Error("ReadSystemList()", "Could not open file %s",
               GetCalibFileName("Systems").Data());
      }
   }
   // if any runs are not associated with any system
   // we create an 'unknown' system and associate it to all runs
   KVDBSystem* sys = 0;
   TIter nextRun(GetRuns());
   KVDBRun* run;
   while ((run = (KVDBRun*)nextRun())) {
      if (!run->GetSystem()) {
         if (!sys) {
            sys = new KVDBSystem("[unknown]");
            AddSystem(sys);
         }
         sys->AddRun(run);
      }
   }
 
   // rehash the record table now that all names are set
   fSystems->Rehash();
}
 
 
 
 
void KVExpDB::WriteSystemsFile() const
{
   //Write the 'Systems.dat' file for this database.
   //The actual name of the file is given by the value of the environment variable
   //[dataset_name].INDRADB.Systems (if it exists), otherwise the value of
   //INDRADB.Systems is used. The file is written in the
   //$KVROOT/[dataset_directory] directory.
 
   ofstream sysfile;
   KVBase::SearchAndOpenKVFile(GetDBEnv("Systems"), sysfile, fDataSet.Data());
   TIter next(GetSystems());
   KVDBSystem* sys;
   TDatime now;
   sysfile << "# " << GetDBEnv("Systems") << " file written by "
           << ClassName() << "::WriteSystemsFile on " << now.AsString() << endl;
   cout << GetDBEnv("Systems") << " file written by "
        << ClassName() << "::WriteSystemsFile on " << now.AsString() << endl;
   while ((sys = (KVDBSystem*)next())) {
      if (strcmp(sys->GetName(), "[unknown]")) { //do not write dummy 'unknown' system
         sys->Save(sysfile);
         sysfile << endl;
      }
   }
   sysfile.close();
}
 
 
 
 
void KVExpDB::Save(const Char_t* what)
{
   //Save (in the appropriate text file) the informations on:
   // what = "Systems" : write Systems.dat file
   // what = "Runlist" : write Runlist.csv
   if (!strcmp(what, "Systems")) WriteSystemsFile();
   else if (!strcmp(what, "Runlist")) WriteRunListFile();
}
 
 
 
 
 
void KVExpDB::WriteRunListFile() const
{
   //Write a file containing a line describing each run in the database.
   //The delimiter symbol used in each line is '|' by default.
   //The first line of the file will be a header description, given by calling
   //KVDBRun::WriteRunListHeader() for the first run in the database.
   //Then we call KVDBRun::WriteRunListLine() for each run.
   //These are virtual methods redefined by child classes of KVDBRun.
 
   ofstream rlistf;
   KVBase::SearchAndOpenKVFile(GetDBEnv("Runlist"), rlistf, fDataSet.Data());
   TDatime now;
   rlistf << "# " << GetDBEnv("Runlist") << " file written by "
          << ClassName() << "::WriteRunListFile on " << now.AsString() << endl;
   cout << GetDBEnv("Runlist") << " file written by "
        << ClassName() << "::WriteRunListFile on " << now.AsString() << endl;
   if (GetRuns() && GetRuns()->GetEntries() > 0) {
      TIter next_run(GetRuns());
      //write header in file
      ((KVDBRun*) GetRuns()->At(0))->WriteRunListHeader(rlistf, GetDBEnv("Runlist.Separator")[0]);
      KVDBRun* run;
      while ((run = (KVDBRun*) next_run())) {
 
         run->WriteRunListLine(rlistf, GetDBEnv("Runlist.Separator")[0]);
 
      }
   }
   else {
      Warning("WriteRunListFile()", "run list is empty !!!");
   }
   rlistf.close();
}
 
 
 
 
Bool_t KVExpDB::OpenCalibFile(const Char_t* type, ifstream& fs) const
{
   //Find and open calibration parameter file of given type. Return kTRUE if all OK.
   //types are defined in $KVROOT/KVFiles/.kvrootrc by lines such as (use INDRA as example)
   //
   //# Default name for file describing systems for each dataset.
   //INDRADB.Systems:     Systems.dat
   //
   //A file with the given name will be looked for in the dataset calibration file
   //directory given by GetDataSetDir()
   //
   //Filenames specific to a given dataset may also be defined:
   //
   //INDRA_camp5.INDRADB.Pedestals:      Pedestals5.dat
   //
   //where 'INDRA_camp5' is the name of the dataset in question.
 
   return KVBase::SearchAndOpenKVFile(GetCalibFileName(type), fs, GetDataSetDir());
}
 
 
 
 
Bool_t KVExpDB::FindCalibFile(const Char_t* type, TString& fullpath, const TString& array_name) const
{
   //Find calibration parameter file of given type. Return kTRUE if all OK.
   //In this case fullpath contains the full path to the file.
   //
   //Types are defined in $KVROOT/KVFiles/.kvrootrc by lines such as (use INDRA as example)
   //
   //~~~~
   //# Default name for file describing systems for each dataset.
   //INDRADB.Systems:     Systems.dat
   //~~~~
   //
   //A file with the given name will be looked for in the dataset calibration file
   //directory given by GetDataSetDir()
   //
   //Filenames specific to a given dataset may also be defined:
   //
   //~~~~
   //INDRA_camp5.INDRADB.Pedestals:      Pedestals5.dat
   //~~~~
   //
   //where 'INDRA_camp5' is the name of the dataset in question.
   //
   //If 'array_name' is given and no file is found with this filename,
   //we try looking for a file called 'array_name.filename', i.e. if a filename
   //
   //~~~~
   //INDRADB.Pressures:  ChIoPressures.dat
   //~~~~
   //
   //is defined, but no ChIoPressures.dat file exists for the dataset, if array_name="INDRA"
   //then we look for a file called INDRA.ChIoPressures.dat
 
   auto cal_file_name = GetCalibFileName(type);
   if (cal_file_name.IsNull()) return kFALSE;
   if (KVBase::SearchKVFile(cal_file_name, fullpath, GetDataSetDir()))
      return kTRUE;
   if (!array_name.IsNull()) {
      cal_file_name.Prepend(Form("%s.", array_name.Data()));
      return KVBase::SearchKVFile(cal_file_name, fullpath, GetDataSetDir());
   }
   return kFALSE;
}
 
 
 
 
 
void KVExpDB::PrintRuns(KVNumberList& nl) const
{
   // Print compact listing of runs in the number list like this:
   //
   // root [9] gIndraDB->PrintRuns("8100-8120")
   // RUN     SYSTEM                          TRIGGER         EVENTS          COMMENTS
   // ------------------------------------------------------------------------------------------------------------------
   // 8100    129Xe + 58Ni 8 MeV/A            M>=2            968673
   // 8101    129Xe + 58Ni 8 MeV/A            M>=2            969166
   // 8102    129Xe + 58Ni 8 MeV/A            M>=2            960772
   // 8103    129Xe + 58Ni 8 MeV/A            M>=2            970029
   // 8104    129Xe + 58Ni 8 MeV/A            M>=2            502992          disjonction ht chassis 1
   // 8105    129Xe + 58Ni 8 MeV/A            M>=2            957015          intensite augmentee a 200 pA
 
   printf("RUN\tSYSTEM\t\t\t\tTRIGGER\t\tEVENTS\t\tCOMMENTS\n");
   printf("------------------------------------------------------------------------------------------------------------------\n");
   nl.Begin();
   while (!nl.End()) {
      KVDBRun* run = GetDBRun(nl.Next());
      if (!run) continue;
      printf("%4d\t%-30s\t%s\t\t%llu\t\t%s\n",
             run->GetNumber(), (run->GetSystem() ? run->GetSystem()->GetName() : "            "), run->GetTriggerString(),
             run->GetEvents(), run->GetComments());
   }
}
 
 
 
 
void KVExpDB::cd()
{
   gExpDB = this;
}
 
 
 
 
 
KVExpDB* KVExpDB::MakeDataBase(const Char_t* name, const Char_t* datasetdir, Bool_t minimal)
{
   // Static function which will create and 'Build' the database object corresponding to 'name'
   // defined as 'Plugin' classes for KVExpDB:
   //
   //~~~
   //      +Plugin.KVExpDB:    [name]    [class]   [library]    "class(const char*)"
   //~~~
   //
   // The 'name' corresponds to the name of the dataset, it is the argument given to the
   // constructor and stored in member variable fDataSet.
   //
   // \param[in] minimal if =true, only a minimal database with runs/systems infos will be built
 
   TPluginHandler* ph;
   if (!(ph = KVBase::LoadPlugin("KVExpDB", name))) {
      return 0;
   }
   //execute constructor/macro for database
   auto mda = (KVExpDB*) ph->ExecPlugin(1, name);
   mda->SetDataSetDir(datasetdir);
 
   mda->Build(minimal);
 
   return mda;
}
 
 
 
 
ULong64_t KVExpDB::GetTotalEvents(int first_run, int last_run) const
{
   // Return total number of events in range [first_run,last_run]
   // (if last_run=-1, go to last known run)
 
   ULong64_t total = 0;
   TIter next(GetRuns());
   KVDBRun* dbr;
   while ((dbr = (KVDBRun*)next())) {
      if (dbr->GetNumber() >= first_run) {
         if ((last_run > 0 && dbr->GetNumber() <= last_run)
               || last_run == -1) total += dbr->GetEvents();
      }
   }
   return total;
}
 
 
 
 
ULong64_t KVExpDB::GetTotalEvents(const KVString& system) const
{
   // Return total number of events for given system
 
   if (!GetSystem(system)) {
      Error("GetTotalEvents", "No system with name : %s", system.Data());
      return 0;
   }
   ULong64_t total = 0;
   for (auto dbr : GetSystem(system)->GetRuns()) {
      total += dynamic_cast<KVDBRun*>(dbr)->GetEvents();
   }
   return total;
}
 
 
 
 
run_index_t KVExpDB::GetRunIndexFromGlobalRunFileNumber(int glob_runfile) const
{
   // \param[in] glob_runfile global runfile number of a runfile (given by KVDBRunFile::GetGlobalRunFileNumber())
   // \returns run/index information of runfile with given global runfile number
   //
   // \note only _good_ runfiles (with KVDBRunFile::IsBad() == kFALSE) have a global runfile number,
   // the _bad_ runfiles are not included in the list
 
   if (fMapGlobRunFileNumberToRunIndex.empty()) {
      // have to regenerate map on first call after reading database in from file.
      // this is because no Streamer can be generated for run_index_t and so the
      // structure cannot be written to disk
      TIter it_run(GetRuns());
      KVDBRun* r;
      while ((r = (KVDBRun*)it_run())) {
         for (auto& ri : r->GetRunIndexList()) {
            auto& rf = GetDBRunFile(ri);
            if (!rf.IsBad())
               fMapGlobRunFileNumberToRunIndex[rf.GetGlobalRunFileNumber()] = ri;
         }
      }
   }
   return fMapGlobRunFileNumberToRunIndex[glob_runfile];
}
 
 
 
 
run_index_list KVExpDB::SetRunIndexListFromString(const TString& rlist) const
{
   // \param[in] rlist string containing run numbers, runfile identifiers (run+index), and ranges of both
   // \returns run_index_list containing the sequence of all _good_ runfiles corresponding to the input string
   //
   // \note any _bad_ runfiles (with KVDBRunFile::IsBad() == kTRUE) are excluded.
   //
   // The string can contain individual run numbers, runfile (runindex) identifiers,
   // either in sequences (separated by spaces and/or commas), or as part of a range,
   // i.e. two values separated by '-'.
   //
   // A number without a '.index' index suffix is interpreted as a run number,
   // and corresponds to all runfiles of the given run. A number _with_ '.index'
   // suffix is interpreted as meaning the given individual runfile of the corresponding
   // run.
   //
   // A number with a '.0' (zero) index is interpreted as the first runfile of the
   // given run (note that normally such a runfile is represented by the run number alone,
   // without a suffix).
   //
   // Examples:
   //
   //~~~
   //"40-43"    => all good runfiles from runs 40, 41, 42 & 43
   //
   //"40-43.0"  => all good runfiles from runs 40, 41 & 42, plus the first runfile of run 43
   //
   //"40.3-42.7" => all good runfiles from 40.3 to 42.7
   //~~~
 
   run_index_list rilist;
 
   KVString _the_list(rlist);
   _the_list.Begin(" ,");
   while (!_the_list.End()) {
      KVString item = _the_list.Next(kTRUE);
      if (item.GetNValues("-") == 2) {
         // this is a range
         item.Begin("-");
         std::vector<int> limits(2);
         for (int i = 0; i < 2; ++i) {
            auto titbit = item.Next(kTRUE);
            if (titbit.Contains(".")) {
               // this is a runindex identifier
               limits[i] = GetDBRunFile(titbit).GetGlobalRunFileNumber();
            }
            else {
               // this is a run number
               // if at beginning of range, we start from first good file
               // if at end of range, we include everything up to last good file
               run_index_t ri(titbit);
               if (i == 0)
                  limits[i] = GetDBRun(ri.run())->GetFirstGoodFile().value().get().GetGlobalRunFileNumber();
               else limits[i] = GetDBRun(ri.run())->GetLastGoodFile().value().get().GetGlobalRunFileNumber();
            }
         }
         for (int ri = limits[0]; ri <= limits[1]; ++ri)
            rilist.Add(GetRunIndexFromGlobalRunFileNumber(ri));
      }
      else {
         // individual run number (no index) or runfile identifier
         if (item.Contains(".")) {
            // this is a runindex identifier
            rilist.Add(run_index_t{item});
         }
         else {
            // run number: add all files from first to last (good) file of run
            auto run_num = run_index_t{item}.run();
            for (int i = GetDBRun(run_num)->GetFirstGoodFile().value().get().GetGlobalRunFileNumber();
                  i <= GetDBRun(run_num)->GetLastGoodFile().value().get().GetGlobalRunFileNumber(); ++i)
               rilist.Add(GetRunIndexFromGlobalRunFileNumber(i));
         }
      }
   }
   return rilist;
}
 
 
 
 
 
TString KVExpDB::GetDBEnv(const Char_t* type) const
{
   //Will look for
   //~~~
   //gEnv->GetValue("name_of_dataset.fDBType.type")
   //~~~
   //then
   //~~~
   //gEnv->GetValue("fDBType.type")
   //~~~
   //if no dataset-specific value is found,
   //then
   //~~~
   //gEnv->GetValue("EXPDB.type")
   //~~~
   //if no database-specific value is found
 
   if (fDataSet == "") return "";
   auto p = KVBase::GetDataSetEnv(fDataSet, Form("%s.%s", fDBType.Data(), type), "");
   if (!p.Length()) return KVBase::GetDataSetEnv(fDataSet, Form("EXPDB.%s", type), "");
   return p;
}
 
 
 
 
void KVExpDB::ReadComments()
{
   // Looks for file with name given by one of the following variables:
   //
   //    [DBtype].Comments
   //    [dataset].[DBtype].Comments
   //
   // and opens it to read and add comments on runs.
   // Format of file is:
   //
   //    run=3830-3836 | really amazing data in these runs
   //
 
   TString comments_file = GetCalibFileName("Comments");
   if (comments_file == "") return;
   TString fullpath;
   if (!FindCalibFile("Comments", fullpath)) return;
   Info("ReadComments", "Reading run comments in file %s...", fullpath.Data());
 
   KVFileReader fr;
   if (!fr.OpenFileToRead(fullpath)) {
      Error("ReadComments", "Problem opening file %s", fullpath.Data());
      return;
   }
 
   while (fr.IsOK()) {
      fr.ReadLine("|");
      if (fr.GetCurrentLine().BeginsWith("#")) {
 
      }
      else if (fr.GetCurrentLine() == "") {
 
      }
      else {
         if (fr.GetNparRead() == 2) {
            KVString srun(fr.GetReadPar(0));
            srun.Begin("=");
            srun.Next();
            KVNumberList lruns(srun.Next());
            KVString comments(fr.GetReadPar(1));
            lruns.Begin();
            while (!lruns.End()) {
               Int_t run = lruns.Next() / index_multiplier;
               KVDBRun* dbrun = GetDBRun(run);
               if (dbrun)
                  dbrun->SetComments(comments.Data());
            }
         }
      }
   }
}