JEditPMTParameters.cc File Reference

Auxiliary program to edit PMT parameters map. More...

#include <string>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <vector>
#include "JLang/JLangToolkit.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JPMTParametersMap.hh"
#include "JDetector/JPMTReadoutAddress.hh"
#include "JDetector/JPMTPhysicalAddress.hh"
#include "JDetector/JModuleAddressMap.hh"
#include "JDetector/JDetectorAddressMap.hh"
#include "JDetector/JDetectorAddressMapToolkit.hh"
#include "JDetector/JPMTAnalogueSignalProcessor.hh"
#include "JSupport/JMeta.hh"
#include "JTools/JRange.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Auxiliary program to edit PMT parameters map.

Syntax:

    -A "<PMT physical address> (set|add|sub|mul|div) <key> <value>"
    -M "<PMT identifier>       (set|add|sub|mul|div) <key> <value>"

The PMT physical address corresponds to the data structure JDETECTOR::JPMTPhysicalAddress and the PMT identifier to JDETECTOR::JPMTIdentifier.
The key corresponds to one of the data members of the JDETECTOR::JPMTParameters data structure.

Note that in the absence of option -a, the detector identifier should be specified using option -D so to obtain the correct PMT address mapping.

Multiple options -A will be processed in order of appearance.

Author

mdejong

Definition in file JEditPMTParameters.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 191 of file JEditPMTParameters.cc.

{
  using namespace std;
  using namespace JPP;

  typedef JRange<double>                 JRange_t;

  string               detectorFile;
  int                  detectorID;
  JPMTParametersMap    parameters;
  vector< JModifier<JPMTPhysicalAddress> >  mod;
  vector< JModifier<JPMTIdentifier> >       daq;
  double               mu;
  JRange_t             T_ns;
  string               outputFile;
  int                  debug;
  
  try {

    JParser<> zap("Auxiliary program to edit PMT parameters map.");

    zap['a'] = make_field(detectorFile,     "detector file.")                                                   = "";
    zap['D'] = make_field(detectorID,       "detector identifier (in absence of detector file).")               = 0;
    zap['P'] = make_field(parameters,       "PMT simulation data (or corresponding file name)")                 = JPARSER::initialised();
    zap['A'] = make_field(mod,              "PMT parameter modifier by physical address (e.g. B1).")            = JPARSER::initialised();
    zap['M'] = make_field(daq,              "PMT parameter modifier by DAQ address (e.g. <module> <channel>.")  = JPARSER::initialised();
    zap['E'] = make_field(mu,               "expectation value for npe given two-fold coincidence")             = 0.0;
    zap['T'] = make_field(T_ns,             "time-over-threshold range.")                                       = JRange_t();
    zap['o'] = make_field(outputFile,       "output file.");
    zap['d'] = make_field(debug)        = 3;

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


  if (detectorFile != "") {

    // Setting default PMT parameters for given detector.
    
    JDetector detector;

    try {
      load(detectorFile, detector);
    }
    catch(const JException& error) {
      FATAL(error);
    }

    if        (detectorID == 0) {

      detectorID = detector.getID();
      
    } else if (detectorID != detector.getID()) {

      FATAL("Inconsistent detector identifier " << detectorID << " != " << detector.getID() << endl);
    }
      
    
    for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {

      for (unsigned int pmt = 0; pmt != module->size(); ++pmt) {

        const JPMTIdentifier id(module->getID(), pmt);

        if (parameters.find(id) == parameters.end()) {

          DEBUG("Setting default parameters for PMT " << id << endl);
          
          parameters[id] = parameters.getDefaultPMTParameters();
        }
      }
    }
  }

  if (!mod.empty()) {

    if (!hasDetectorAddressMap(detectorID)) {
      FATAL("Invalid detector identifier " << detectorID << endl);
    }
    
    const JDetectorAddressMap& demo = getDetectorAddressMap(detectorID);
    
    for (JPMTParametersMap::iterator ps = parameters.begin(); ps != parameters.end(); ++ps) {

      const JPMTPhysicalAddress& address = demo.get(ps->first);

      for (vector< JModifier<JPMTPhysicalAddress> >::const_iterator i = mod.begin(); i != mod.end(); ++i) {
        
        if (compare(i->address, address)) {

          DEBUG("Modifying parameters for PMT " << ps->first << ' ' << i->action << ' ' << i->key << ' ' << i->value << endl);

          if (!i->apply(ps->second)) {
            ERROR("No valid action: " << *i << endl);
          }
        }
      }
    }
  }
  
  if (!daq.empty()) {

    for (JPMTParametersMap::iterator ps = parameters.begin(); ps != parameters.end(); ++ps) {

      for (vector< JModifier<JPMTIdentifier> >::const_iterator i = daq.begin(); i != daq.end(); ++i) {
        
        if (compare(ps->first, i->address)) {

          DEBUG("Modifying parameters for PMT " << ps->first << ' ' << i->action << ' ' << i->key << ' ' << i->value << endl);

          if (!i->apply(ps->second)) {
            ERROR("No valid action: " << *i << endl);
          }
        }
      }
    }
  }


  if        (mu > 0.0) {
    
    DEBUG("Correct measured QE for two-hit probability " << mu << endl);

    try {
      parameters.convertHitProbabilityToQE(mu);
    }
    catch(const JException& error) {
      FATAL(error.what());
    }
    
  } else if (mu < 0.0) {

    FATAL("Invalid expection value for two-hit probability " << mu << endl);
  }
  
  
  if (T_ns != JRange_t()) {

    DEBUG("Correct measured QE for time-over-threshold range " << T_ns << endl);

    const int NPE = 1;

    for (JPMTParametersMap::iterator i = parameters.begin(); i != parameters.end(); ++i) {

      const JPMTAnalogueSignalProcessor cpu(i->second);
      
      i->second.QE *= (cpu.getIntegralOfChargeProbability(i->second.threshold,
                                                          cpu.getNPE(T_ns.getUpperLimit(), NPE),
                                                          NPE)
                       /
                       cpu.getIntegralOfChargeProbability(cpu.getNPE(T_ns.getLowerLimit(), NPE),
                                                          cpu.getNPE(T_ns.getUpperLimit(), NPE),
                                                          NPE));
    }
  }
  
  
  if (outputFile != "") {

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

    ofstream out(outputFile.c_str());

    out << parameters << endl;

    out.close();
  }
}