software/JDB/JPMTThreshold.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <vector>
#include <map>

#include "TF1.h"

#include "JDB/JDB.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JDBToolkit.hh"
#include "JDB/JPMTThreshold.hh"
#include "JDB/JDetectorIntegration.hh"

#include "JDetector/JPMTParametersMap.hh"

#include "JSystem/JStat.hh"

#include "JSupport/JMeta.hh"

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


/**
 * \file
 * Auxiliary program to set PMT thresholds according DAQ configuration of a given data taking run.
 * \author mdejong
 */
int main(const int argc, const char * const argv[])
{
  using namespace std;
  using namespace JPP;

  JServer  server;
  string   usr;
  string   pwd;
  string   cookie;
  string   detid;
  int      run;
  string   pmtFile;
  string   formula;
  int      debug;

  try {

    JParser<> zap("Auxiliary program to set PMT thresholds according DAQ configuration of a given data taking run.");

    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);
    zap['r'] = make_field(run);
    zap['P'] = make_field(pmtFile,    "PMT parameters file")  =  "";
    zap['F'] = make_field(formula,    "transfer function")    =  "((x * 1.0)/255.0 * 1.6 + 0.8 - 0.989) / 0.44"; // "((x * 1.0)/255.0 * 1.6 + 0.8 - 1) / 0.44";
    zap['d'] = make_field(debug)      = 1;

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


  JPMTParametersMap parameters;

  if (pmtFile != "" && getFileStatus(pmtFile.c_str())) {
    parameters.load(pmtFile.c_str());
  }


  TF1 f1("user", formula.c_str());

  if (!f1.IsValid()) {
    FATAL("Formula <" << formula << "> invalid." << endl);
  }


  typedef vector<JDetectorIntegration> detector_type;

  JPMTThreshold getPMTThreshold;
  detector_type detector;


  try {

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

    const int id = getDetector<int>(detid);

    getPMTThreshold.configure(id, run);

    DEBUG(getPMTThreshold);

    ResultSet& rs = getResultSet(getTable<JDetectorIntegration>(), getSelector<JDetectorIntegration>(id));

    for (JDetectorIntegration parameters; rs >> parameters; ) {
      detector.push_back(parameters);
    }

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


  for (vector<JDetectorIntegration>::const_iterator i = detector.begin(); i != detector.end(); ++i) {

    if (i->PMTID != -1) {

      try {

        const JPMTThreshold::result_type threshold = getPMTThreshold(i->PMTUPI);

        const double value = f1.Eval((double) threshold.value);

        if (debug >= debug_t) {

          cout << FILL(4,'0') << i->DUID         << FILL() << '.' 
               << FILL(2,'0') << i->FLOORID      << FILL() << '[' 
               << FILL(2,'0') << i->CABLEPOS     << FILL() << ']'
               << ' '
               << setw(3)     << threshold.value << ' '
               << FIXED(5,2)  << value        << ' ';

          if (threshold.is_default) {
            cout << "* (" << FIXED(5,2) << JPMTParameters().threshold << ")";
          }

          cout << endl;
        }

        if (!threshold.is_default) {

          const JPMTIdentifier id(i->DOMID, i->CABLEPOS);

          parameters[id].threshold = value;
        }
      }
      catch(const exception& error) {
        WARNING("Error at " 
                << FILL(4,'0') << i->DUID         << FILL() << '.' 
                << FILL(2,'0') << i->FLOORID      << FILL() << '[' 
                << FILL(2,'0') << i->CABLEPOS     << FILL() << ']' 
                << ' '
                << error.what() << endl);
      }
    }
  }

  if (pmtFile != "") {

    parameters.comment.add(JMeta(argc, argv));

    parameters.store(pmtFile.c_str());

  } else {

    cout << parameters << endl;
  }
}