JCheckHV.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
#include <limits>
#include <map>
 
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
 
#include "JDB/JDB.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JDBToolkit.hh"
#include "JDB/JDetectorIntegration.hh"
#include "JDB/JProductRouter.hh"
#include "JDB/JAllParams.hh"
#include "JDB/JDatalogNumbers.hh"
#include "JDB/JPMTHVRunSettings.hh"
#include "JDB/JSupport.hh"
 
#include "JGizmo/JManager.hh"
#include "JLang/JLangToolkit.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;
 
  string          usr;
  string          pwd;
  string          cookie;
  string          outputFile;
  string          detid;
  int             run;
  long long int   numberOfRows;
  double          precision;
  int             debug;
 
  try {
 
    JParser<> zap("Auxiliary program to check HV.");
    
    zap['u'] = make_field(usr)           = "";
    zap['!'] = make_field(pwd)           = "";
    zap['C'] = make_field(cookie)        = "";
    zap['o'] = make_field(outputFile)    = "hv.root";
    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['p'] = make_field(precision)     = 0.5;
    zap['d'] = make_field(debug)         = 2;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  JManager<JLocation_t, TH1D> manager(new TH1D("HV[%]", NULL, 200, -1500.0, +1500.0));
 
  try {
 
    JDB::reset(usr, pwd, cookie);
 
    if (is_integer(detid)) {
      detid = getDetector(to_value<int>(detid));
    }
 
    JDetectorIntegration detector;
 
    {
      DEBUG("Reading database table " << getTable<JDetectorIntegration>() << endl); 
 
      ResultSet& rs  = getResultSet(getTable<JDetectorIntegration>());
 
      if (! (rs >> detector)) {
        THROW(JDatabaseException, "Error reading " << getTable<JDetectorIntegration>());
      }
    }
 
    DEBUG("Number of integrated products (before) " << right << detector.size() << endl);
    
    detector.configure(detid);
    
    DEBUG("Number of integrated products (after)  " << right << detector.size() << endl);
 
    const JProductRouter router(detector);
    
    map<JLocation_t, double> HV;
    
    {
      ResultSet& rs  = getResultSet(getTable<JPMTHVRunSettings>(), getSelector<JPMTHVRunSettings>(detid, run));
 
      for (JPMTHVRunSettings parameters; rs >> parameters; ) {
 
        const JLocation_t location(parameters.DUID, parameters.FLOORID, parameters.PMTINTID);  
 
        HV[location] = parameters.HV_VALUE;
      }
      
      rs.Close();
    }
 
    if (debug >= debug_t) {
 
      cout << "HV settings:" << endl;
 
      for (map<JLocation_t, double>::const_iterator i = HV.begin(); i != HV.end(); ++i) {
        cout << i->first << ' ' << FIXED(6,1) << i->second << endl;
      }
    }
 
    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         upi      = parameters.SOURCE_NAME;
      const JLocation_t& location = router.getLocation(upi);
 
      DEBUG(location << ' ' << FIXED(6,1) << HV[location] << ' ' << FIXED(6,1) << parameters.DATA_VALUE << endl);
 
      if (location != JLocation_t()) {
 
        const double hv = HV[location] - parameters.DATA_VALUE ;
 
        manager[location]->Fill(hv);
 
        if (fabs(hv) > precision) {
          WARNING("HV " << location << " set/read " << FIXED(6,1) << HV[location] << "/" << FIXED(6,1) << parameters.DATA_VALUE << endl);
        }
      }
      STATUS(endl);
    }
  }
  catch(const exception& error) {
    FATAL(error.what() << endl);
  }
 
 
  TFile out(outputFile.c_str(), "recreate");
 
  out << manager;
 
  out.Write();
  out.Close();
}