KVExpDB.cpp

Go to the documentation of this file.

//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 <iostream>
#include "KVUnownedList.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");
}
 
 
 
 
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.
 
   for (UInt_t ii = first_run; ii <= last_run; ii++) {
      LinkRecordToRun(rec, ii);
   }
}
 
 
 
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();
   while (!nl.End()) {
      Int_t rr = nl.Next();
      //Info("LinkRecordToRunRange","run number %d",rr);
      LinkRecordToRun(rec, 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;
 
   for (UInt_t ru_ra = 0; ru_ra < rr_number; ru_ra++) {
      UInt_t first_run = run_ranges[ru_ra][0];
      UInt_t last_run = run_ranges[ru_ra][1];
      for (UInt_t i = first_run; i <= last_run; i++) {
         KVDBRecord* run = GetDBRun(i);
         while ((rec = (KVDBRecord*) next())) {
            if (run)
               rec->AddLink("Runs", run);
         }
         next.Reset();
      }
   }
}
 
 
 
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()
{
   //Reads list of systems with associated run ranges, creates KVDBSystem
   //records for these systems, and links them to the appropriate KVDBRun
   //records using LinkListToRunRanges.
   //
   //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);
         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
{
   //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.
   return KVBase::SearchKVFile(GetCalibFileName(type), fullpath, GetDataSetDir());
}
 
 
 
 
 
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)
{
   //Static function which will create and 'Build' the database object corresponding to 'name'
   //These are defined as 'Plugin' objects in the file $KVROOT/KVFiles/.kvrootrc :
   //
   //      Plugin.KVExpDB:    INDRA_camp1    KVDataBase1     KVIndra    "KVDataBase1()"
   //      +Plugin.KVExpDB:    INDRA_camp2    KVDataBase2     KVIndra    "KVDataBase2()"
   //      +Plugin.KVExpDB:    INDRA_camp4    KVDataBase4     KVIndra    "KVDataBase4()"
   //      +Plugin.KVExpDB:    INDRA_camp5    KVDataBase5     KVIndra5    "KVDataBase5()"
   //
   //The 'name' ("INDRA_camp1" etc.) corresponds to the name of a dataset in $KVROOT/KVFiles/manip.list
   //This name is stored in member variable fDataSet.
   //The constructors/macros used have arguments (const Char_t* name)
 
   //does plugin exist for given name ?
   TPluginHandler* ph;
   if (!(ph = KVBase::LoadPlugin("KVExpDB", name))) {
      return 0;
   }
   //execute constructor/macro for database
   KVExpDB* mda = (KVExpDB*) ph->ExecPlugin(1, name);
   mda->SetDataSetDir(datasetdir);
 
   mda->Build();
 
   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;
   }
   TIter it(GetSystem(system)->GetRuns());
   ULong64_t total = 0;
   KVDBRun* dbr;
   while ((dbr = (KVDBRun*)it())) {
      total += dbr->GetEvents();
   }
   return total;
}
 
 
 
 
 
TString KVExpDB::GetDBEnv(const Char_t* type) const
{
   //Will look for gEnv->GetValue name "name_of_dataset.fDBType.type",
   //then "fDBType.type" if no dataset-specific value is found,
   //then "EXPDB.type" if no database-specific value is found
 
   if (fDataSet == "") return "";
   TString 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();
               KVDBRun* dbrun = GetDBRun(run);
               if (dbrun)
                  dbrun->SetComments(comments.Data());
            }
         }
      }
   }
}