JAttachPMTs.cc

Go to the documentation of this file.

#include <string>
 
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JSupport/JMeta.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * \file
 * 
 * Auxiliary program to attach PMTs to detector.
 *
 * \author mdejong
 */
int main(int argc, char **argv)
{
 
  using namespace std;
  using namespace JPP;
  
  string    detectorFile_a;
  string    detectorFile_b;
  string    outputFile;
  bool      squash;
  int       debug;
  
  try { 
 
    JParser<> zap("Auxiliary program to attach PMTs to detector.");
    
    zap['a'] = make_field(detectorFile_a, "detector file w/o  PMTs");
    zap['f'] = make_field(detectorFile_b, "detector file with PMTs");
    zap['o'] = make_field(outputFile,     "detector file with PMTs");
    zap['q'] = make_field(squash,         "squash meta data");
    zap['d'] = make_field(debug)          = 0;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
   
  JDetector detector_a;
  JDetector detector_b;
 
  try {
    load(detectorFile_a, detector_a);
    load(detectorFile_b, detector_b);
  }
  catch(const JException& error) {
    FATAL(error);
  }
 
  if (squash) {
    detector_a.comment.clear();
  }
 
  detector_a.comment.add(JMeta(argc,argv));
 
  if (detector_a.setToLatestVersion()) {
    NOTICE("Set detector version to " << detector_a.getVersion() << endl);
  }  
 
  for (JDetector::iterator module_b = detector_b.begin(); module_b != detector_b.end(); ++module_b) {
 
    for (JDetector::iterator module_a = detector_a.begin(); module_a != detector_a.end(); ++module_a) {
 
      if (module_a->getID() == module_b->getID()) {
 
        DEBUG("Attach PMTs in module " << setw(10) << module_a->getID() << endl);
 
        module_b->set ( module_a->getPosition());
        module_b->swap(*module_a);
 
        break;
      }
    }
  }
 
  try {
    store(outputFile, detector_a);
  }
  catch(const JException& error) {
    FATAL(error);
  }
}