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);
}