software/JAcoustics/JToA.cc

Go to the documentation of this file.

#include <iomanip>
#include <iostream>
 
#include "JAcoustics/JToA.hh"
#include "JAcoustics/JSupport.hh"
#include "JDB/JDB.hh"
#include "JDB/JDBToolkit.hh"
#include "JDB/JDatabaseObjectIterator.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSupport.hh"
#include "JDB/JToAshort.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JMeta.hh"
#include "Jeep/JMessage.hh"
#include "Jeep/JParser.hh"
#include "dbclient/KM3NeTDBClient.h"
 
namespace JROOT {
 
  // make ROOT output file flat
 
  template<> class JRootCreateFlatTree<JACOUSTICS::JToA> : public std::true_type {};
}
 
/**
 * \file
 *
 * Auxiliary program to convert acoustics data from data base into ROOT format.
 *
 * Example usage :
 * <pre>
 * JToA -@ "minrun = 16514"      \
 *      -@ "maxrun = 16514"      \
 *      -@ "detid  = D1ORCA015"  \
 *      -o KM3NeT_00000146_00016514_toa.root
 * </pre>
 * \author mdejong
 */
int main(int argc, char** argv)
{
  using namespace std;
  using namespace JPP;
 
  JFileRecorder<JTYPELIST<JToA, JMeta>::typelist> outputFile;
 
  JServer   server;
  string    usr;
  string    pwd;
  string    cookie;
  JSelector selection;
  int       debug;
 
  try {
 
    JParser<> zap("Auxiliary program to convert acoustics data from data base into ROOT format.");
 
    // clang-format off
    zap['o'] = make_field(outputFile);
    zap['s'] = make_field(server)     = getServernames();
    zap['u'] = make_field(usr)        = "";
    zap['!'] = make_field(pwd)        = "";
    zap['C'] = make_field(cookie)     = "";
    zap['@'] = make_field(selection)  = JPARSER::initialised();
    zap['d'] = make_field(debug)      = 1;
    // clang-format on
 
    zap(argc, argv);
  }
  catch (const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  try {
 
    JDB::reset(usr, pwd, cookie);
 
    JDBToolkit::initialise(getDetector);
 
    outputFile.open();
 
    outputFile.put(JMeta(argc, argv));
 
    JToA toa;
 
    for (JDatabaseObjectIterator<JToAshort> in(JToAshort::getName(), selection); in.hasNext();) {
 
      const JToAshort* p = in.next();
 
      toa.DETID                = getDetector(p->DETID);
      toa.RUN                  = p->RUN;
      toa.DOMID                = p->DOMID;
      toa.WAVEFORMID           = p->EMITTERID;
      toa.TOA_NS               = p->getAbsoluteToA_ns();
      toa.SECONDS              = p->getDAQFrameSeconds();
      toa.TICKS                = p->getDAQFrameTicks();
      toa.QUALITYFACTOR        = p->QUALITYFACTOR;
      toa.QUALITYNORMALISATION = p->QUALITYNORMALISATION;
 
      outputFile.put(toa);
    }
 
    outputFile.close();
  }
  catch (const exception& error) {
    FATAL(error.what() << endl);
  }
}