software/JDB/JCheckHV.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
#include <limits>
#include <vector>
#include <map>
 
#include "JDB/JDB.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JDBToolkit.hh"
#include "JDB/JDetectorIntegration.hh"
#include "JDB/JDatalogNumbers.hh"
#include "JDB/JPMTHV.hh"
#include "JDB/JLocation_t.hh"
#include "JDB/JSupport.hh"
 
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * \file
 *
 * Auxiliary program to check HV.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
 
  JServer         server;
  string          usr;
  string          pwd;
  string          cookie;
  string          detid;
  int             run;
  long long int   numberOfRows;
  double          precision;
  int             qaqc;
  int             debug;
 
  try {
 
    JParser<> zap("Auxiliary program to check HV."\
                  "\nFor detailed ROOT output, use corresponding example.");
    
    zap['s'] = make_field(server)        = getServernames();
    zap['u'] = make_field(usr)           = "";
    zap['!'] = make_field(pwd)           = "";
    zap['C'] = make_field(cookie)        = "";
    zap['D'] = make_field(detid,        "Detector identifier");
    zap['R'] = make_field(run,          "Run");
    zap['n'] = make_field(numberOfRows, "Maximal number of rows")         = numeric_limits<long long int>::max();
    zap['e'] = make_field(precision)     = 0.5;
    zap['Q'] = make_field(qaqc)          = 0;
    zap['d'] = make_field(debug)         = 2;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  map<JUPI_t, size_t> counts;
  map<JUPI_t, size_t> errors;
  map<JUPI_t, size_t> aliens;
 
  try {
 
    JDB::reset(usr, pwd, cookie);
 
    const int ID = getDetector<int>   (detid);
    detid        = getDetector<string>(detid);
  
    map<JUPI_t, JLocation_t> detector;
 
    {
      ResultSet& rs = getResultSet(getTable<JDetectorIntegration>(), getSelector<JDetectorIntegration>(ID));
 
      for (JDetectorIntegration parameters; rs >> parameters; ) {
 
        if (parameters.PMTUPI.getPBS() == PBS::PMT) {
          detector[parameters.PMTUPI] = JLocation_t(parameters.DUID, parameters.FLOORID, parameters.CABLEPOS);
        }
      }
 
      rs.Close();
    }
    {    
      const JPMTHV HV(ID, run);
 
      JSelector selection(&JDatalogNumbers::PARAMETER_NAME, "pmt_highvolt");
 
      selection += getSelector<JDatalogNumbers>(detid, run, run);
 
      ResultSet& rs  = getResultSet(getTable<JDatalogNumbers>(), selection);
 
      long long int counter = 0;
 
      for (JDatalogNumbers parameters; rs >> parameters && counter != numberOfRows; ++counter) {
 
        STATUS(setw(10) << counter << '\r'); DEBUG(endl);
 
        const JUPI_t      upi      = parameters.SOURCE_NAME;
        const JLocation_t location = detector[upi];
 
        DEBUG(location << ' ' << FIXED(7,1) << HV(upi) << ' ' << FIXED(7,1) << parameters.DATA_VALUE << endl);
 
        if (location.is_valid()) {
 
          counts[upi] += 1;
 
          if (fabs(HV(upi).value - parameters.DATA_VALUE) > precision) {
 
            ERROR("HV " << setw(10) << parameters.UNIXTIME << ' ' << location << " set/read " << FIXED(7,1) << HV(upi) << "/" << FIXED(7,1) << parameters.DATA_VALUE << endl);
 
            errors[upi] += 1;
          }
          
        } else {
 
          aliens[upi] += 1;
        }
        
        STATUS(endl);
      }
 
      rs.Close();
    }
 
    if (debug >= warning_t) {
      for (map<JUPI_t, size_t>::const_iterator i = errors.begin(); i != errors.end(); ++i) {
        cout << "Error HV  " << setw(32) << left << i->first << right << ' ' << detector[i->first] << ' ' << setw(3) << i->second << '/' << setw(3) << counts[i->first] << endl;
      }
      for (map<JUPI_t, size_t>::const_iterator i = aliens.begin(); i != aliens.end(); ++i) {
        cout << "Alien UPI " << setw(32) << left << i->first << right << ' ' << setw(3) << i->second << endl;
      }
    }
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
  QAQC(setw(4) << errors.size() << endl);
 
  return 0;
}