JCheckTriggerParameters.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <vector>
#include <map>
 
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/Evt.hh"
 
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQEventIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
 
#include "JTrigger/JTriggerParameters.hh"
 
#include "JSupport/JSingleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JTriggerParametersSupportkit.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
 
#include "JLang/JObjectMultiplexer.hh"
 
#include "JDB/JDB.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JDetectors.hh"
#include "JDB/JRuns.hh"
#include "JDB/JAllParams.hh"
#include "JDB/JRunsetupParams.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * \file
 * Auxiliary program to check trigger parameters in a given ROOT files against the values from the database.
 * \author lnauta
 */
int main(const int argc, const char * const argv[])
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;
 
  JSingleFileScanner<JDAQTypes_t> inputFile;
  JFileRecorder     <JAllTypes_t> outputFile;
  string   usr;
  string   pwd;
  string   cookie;
  int      debug;
 
  try {
 
    JParser<> zap("Auxiliary program to check trigger parameters in a given ROOT files against the values from the database.");
 
    zap['f'] = make_field(inputFile);
    zap['u'] = make_field(usr)        = "";
    zap['!'] = make_field(pwd)        = "";
    zap['C'] = make_field(cookie)     = "";
    zap['o'] = make_field(outputFile) = "";
    zap['d'] = make_field(debug)      = 1;
 
    zap(argc, argv);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  const double Dmax_m = 9999.0;
 
  JTriggerParameters parametersA;
  JTriggerParameters parametersB;
 
 
  DEBUG("Get trigger parameters from input file." << endl);
 
  try {
    parametersA = getTriggerParameters(inputFile);
  }
  catch(const JException& error) {
    FATAL("No trigger parameters from input." << endl);
  }
 
  parametersA.set(Dmax_m);
 
  DEBUG(parametersA);
 
 
  int detector_id = 0;
  int run_number  = 0;
 
  JObjectMultiplexer<JDAQTypes_t, JDAQHeader> in(inputFile);
 
  if (in.hasNext()) {
 
    detector_id = in.next()->getDetectorID();
    run_number  = in.next()->getRunNumber();
 
  } else {
 
    FATAL("File does not contain a valid run number. Aborting" << endl);
  }  
 
 
  try {
    JDB::reset(usr, pwd, cookie);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  string detid = "";
 
  try {
 
    JDetectors buffer;
 
    ResultSet& rs = getResultSet(getTable<JDetectors>(), getSelector<JDetectors>(detector_id));
 
    if (rs >> buffer) {
      detid = buffer.OID;
    }
 
    rs.Close();
 
    DEBUG("Detector " << detector_id << " -> " << detid << endl);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  string rs_oid;
 
  try {
 
    JRuns buffer;
 
    ResultSet& rs = getResultSet(getTable<JRuns>(), getSelector<JRuns>(detector_id, run_number));
 
    if (rs >> buffer) {
      rs_oid = buffer.RUNSETUPID;
    }
 
    rs.Close();
 
    DEBUG("Run setup " << rs_oid << endl);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  map<string, string> umap;
  map<string, string> zmap;
 
  try {
 
    JAllParams buffer;
 
    ResultSet& rs = getResultSet(getTable<JAllParams>(), getSelector<JAllParams>("DATAACQUISITION"));
 
    while (rs >> buffer) {
 
      umap[buffer.OID] = buffer.NAME;
 
      string::size_type pos = buffer.NAME.find(PREFIX_DATAFILTER);
 
      if (pos != string::npos) {
 
        string::size_type len = PREFIX_DATAFILTER.length();
 
        while (pos + len != buffer.NAME.length() && buffer.NAME[pos + len] != '=') {
          ++len;
        }
 
        zmap[buffer.OID] = buffer.NAME.erase(pos, len + 1);
      }
    }
 
    rs.Close();
 
    DEBUG("Parameters map " << zmap.size() << endl);
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  stringstream data;    
 
  try {
 
    JRunsetupParams buffer;
 
    ResultSet& rs = getResultSet(getTable<JRunsetupParams>(), getSelector<JRunsetupParams>(detid, rs_oid));
 
    while (rs >> buffer) {
 
      DEBUG(buffer.PAR_OID << "->" << umap[buffer.PAR_OID] << " = " << buffer.VALUE << endl);
 
      map<string, string>::const_iterator i = zmap.find(buffer.PAR_OID);
 
      if (i != zmap.end()) {
        data << i->second << "=" << buffer.VALUE << ";\n";
      }
    }
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
  DEBUG("Raw data" << endl << "<<" << endl << data.str() << ">>" << endl);
 
 
  data >> parametersB;
 
  parametersB.set(Dmax_m);  
  parametersB.actionAtFileRead();
 
  DEBUG(parametersB);
 
 
  if (outputFile.getFilename() != "") {
 
    DEBUG("Write to " << outputFile << "... " << flush);
 
    outputFile.open();
 
    outputFile.put(parametersB);
    
    JSingleFileScanner<JRemove<JAllTypes_t, JTriggerParameters>::typelist> io(inputFile);
 
    io >> outputFile;
 
    outputFile.close();
 
    DEBUG("OK" << endl);
  }
 
 
  const bool is_equal = (parametersA == parametersB);
 
  if (is_equal)
    NOTICE(GREEN << "The parameters for run " << run_number << " are equivalent."     << RESET << endl);
  else
    ERROR (RED   << "The parameters for run " << run_number << " are NOT equivalent." << RESET << endl);
 
  return (is_equal ? 0 : 1);
}