JPrintTuneHV.cc

Go to the documentation of this file.

#include <iostream>
#include <iomanip>
#include <map>

#include "km3net-dataformat/online/JDAQ.hh"

#include "JSystem/JDate.hh"

#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

#include "JLang/JPredicate.hh"

#include "JSupport/JFilenameSupportkit.hh"

#include "JDetector/JDetectorCalibration.hh"

#include "JCalibrate/JFitToT.hh"

#include "JDB/JDB.hh"
#include "JDB/JDBToolkit.hh"
#include "JDB/JDBincludes.hh"
#include "JDB/JProductRouter.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JDetectorIntegration.hh"
#include "JDB/JDetectorIntegration_t.hh"

#include "JSon/JSon.hh"


/**
 * \file
 *
 * Auxiliary program to print HV-tuning output.
 * \author bjung
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  using namespace KM3NETDAQ;

  string usr;
  string pwd;
  string cookie;

  string inputFile;
  string outputFile;
  
  int    detid;
  int    debug;

  try {

    JParser<> zap("Auxiliary program to print high-voltage tuning results.");

    zap['u'] = make_field(usr)                                                         = "";
    zap['!'] = make_field(pwd)                                                         = "";
    zap['C'] = make_field(cookie)                                                      = "";
    zap['f'] = make_field(inputFile,    "input file (ROOT or JSON output of JTuneHV)");
    zap['o'] = make_field(outputFile,   "output file (for writing failed evaluations") = "";
    zap['D'] = make_field(detid,        "detector identifier");
    zap['d'] = make_field(debug,        "debug")                                       = 2;

    zap(argc, argv);
    
  } catch(const exception &error) {
    
    FATAL(error.what() << endl);
  }

  JDateAndTime           timer;
  JPersons               person;
  JDetectorIntegration_t detector;

  try {

    JDB::reset(usr, pwd, cookie);

    {
      JSelector  selector = getSelector<JPersons>(JDB::get()->User());
      ResultSet& rs       = getResultSet(getTable<JPersons>(), selector);
      rs >> person;
      rs.Close();
    }

    {
      ResultSet& rs = getResultSet(getTable<JDetectorIntegration_t>());
    
      if (! (rs >> detector)) {
        THROW(JDatabaseException, "Error reading " << getTable<JDetectorIntegration_t>());
      }
    }

  } catch (const exception& error) {
    
    FATAL(error.what() << endl);
  }

  detector.configure(getDetector(detid));
  
  const JProductRouter router(detector, PBS::PMT_SEQUENCES);

  string testType;
  size_t Ntotal = 0;
  
  map<int, JHVCalibration> failureMap;
  
  if (isJSONFile(inputFile.c_str())) {
    
    json js;

    ifstream ifs(inputFile.c_str());

    ifs >> js;
    ifs.close();

    JHVCalibration HVcals;
    
    json::const_iterator data = js.find(Tests_t);

    if (data != js.end()) {
    
      HVcals = data->get<JHVCalibration>();
    
    } else {

      ERROR("HV calibration data could not be found!");
    }

    json::const_iterator DBtype = js.find(Test_t + Type_t);

    if (DBtype != js.end()) {

      testType = DBtype->get<string>();
      
    } else {

      WARNING("No database test type specified.");
    }
    
    NOTICE(LEFT (30) << "UPI"  << CENTER(35) << "(identifier / location)" <<
           RIGHT(20) << "HV"   << RIGHT (10)  << "gain" << RIGHT(10) << "status" << endl);
    
    for (JHVCalibration::const_iterator it = HVcals.cbegin(); it != HVcals.cend(); ++it) {

      const JUPI_t&      pmtUPI   = it->getUPI();
      const JUPI_t&      domUPI   = detector.trace(pmtUPI, PBS::DOM).container.getUPI();
      const JLocation_t& location = router.getLocation(pmtUPI);

      int domID = -1;
      
      const JDetectorIntegration_t::range_type range = detector.get(domUPI);

      for (JDetectorIntegration_t::range_const_iterator i = range.first; i != range.second; ++i) {

        const JUPI_t& upi = detector[i->second].content.getUPI();
          
        if (upi.getPBS() == PBS::CLB) {

          domID = getCLBID(upi);

          if (it->result != OK_t) {
            failureMap[domID].push_back(*it);
          }
          
          break;
        }
      }

      NOTICE(LEFT(30) << pmtUPI << "(a.k.a. " << domID << '.' << FILL(2,'0') << location.position << " / " << location << "):" <<
             right << FIXED(20,1) << it->supplyVoltage << FIXED(10,1) << it->PMTgain << RIGHT(10) << it->result << endl);

      ++Ntotal;
    }

  } else {

    ERROR(inputFile << " is not a JSON file." << endl);
  }


  size_t Nfailed  = 0;
  size_t Nmissing = 0;

  JHVCalibration failures;
  
  NOTICE(endl << FILL(105, '-') << " List of failures" << setfill(' ') << endl);
  
  for (map<int, JHVCalibration>::const_iterator i = failureMap.cbegin(); i != failureMap.cend(); ++i) {

    Nfailed += i->second.size();
    
    if (i->second.size() == NUMBER_OF_PMTS) {
      
      Nmissing += NUMBER_OF_PMTS;
      WARNING("No successful calibrations found for module " << i->first << " (dead module?)" << endl);
      
    } else {

      for (JHVCalibration::const_iterator j = i->second.cbegin(); j != i->second.cend(); ++j) {

        const JUPI_t&      pmtUPI   = j->getUPI();
        const JLocation_t& location = router.getLocation(pmtUPI);
        
        NOTICE(LEFT(30) << pmtUPI << "(a.k.a. " << i->first << '.' << FILL(2,'0') << location.position << " / " << location << "):" <<
               right << FIXED(20,1) << j->supplyVoltage << FIXED(10,1) << j->PMTgain << RIGHT(10) << j->result << endl);

        failures.push_back(*j);
      }
    }
  }
  
  NOTICE(endl << FILL(105, '-') << " SUMMARY" << setfill(' ') << endl);
  NOTICE(LEFT(40)  << "Number of evaluated PMTs:" << RIGHT(20) << Ntotal << endl);
  NOTICE(LEFT(40)  << "Number of successful evaluations:" << RIGHT(20) << Ntotal - Nfailed << endl);
  NOTICE(LEFT(40)  << "Number of failed evaluations:" << RIGHT(20) << Nfailed << endl);
  NOTICE(RIGHT(20) << '-' << LEFT(20) << " in missing modules:" << RIGHT(20) << Nmissing << endl);
  NOTICE(RIGHT(20) << '-' << LEFT(20) << " other:" << RIGHT(20) << Nfailed - Nmissing << endl << endl);

  if (!outputFile.empty()) {
  
    json js;

    js[User_t]          = person.LOGIN;
    js[Location_t]      = person.LOCATIONID;
    js[StartTime_t]     = timer.toString();
    js[EndTime_t]       = timer().toString();
    js[Test_t + Type_t] = testType;
    js[Tests_t]         = json(failures);

    ofstream ofs(outputFile.c_str());

    ofs << setw(2) << setprecision(8);
    ofs << js;

    ofs.close();
  }

  return 0;
}