JDestructDetector.cc File Reference

Auxiliary program to decompose detector to separate calibrations. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <fstream>
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorCalibration.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorAddressMapToolkit.hh"
#include "Jeep/JeepToolkit.hh"
#include "Jeep/JPrint.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 decompose detector to separate calibrations.

Author

mdejong

Definition in file JDestructDetector.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 23 of file JDestructDetector.cc.

{
  using namespace std;
  using namespace JPP;

  string   detectorFile;
  string   outputFile;
  double   precision;
  int      debug;

  try {

    JParser<> zap("Auxiliary program to decompose detector to separate calibrations.");

    zap['a'] = make_field(detectorFile, "detector file");
    zap['o'] = make_field(outputFile,   "detector calibrations file (json format)") = "";
    zap['p'] = make_field(precision)     = 1.0e-3;
    zap['d'] = make_field(debug)         = 1;

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


  JDetector detector;

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

  if (!hasDetectorAddressMap(detector.getID())) {
    FATAL("No detector address map for detector identier " << detector.getID() << endl);
  }

  const JDetectorAddressMap& demo = getDetectorAddressMap(detector.getID());


  // Test compatibility of modules with reference(s).

  for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {

    if (module->getFloor() != 0) {

      JModule buffer = getModule(demo.get(module->getID()), module->getID(), module->getLocation());

      const JRotation3D R = getRotation(buffer, *module);

      buffer.rotate(R);

      if (!JModule::compare(buffer, *module, precision)) {
        ERROR("Module " << setw(10) <<  module->getID() << ' ' << module->getLocation() << " incompatible with reference." << endl);
      }
    }
  }


  JPMTCalibration  calibration;
  JModulePosition  position;
  JModuleRotation  rotation;  

  for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {

    const JModule& buffer = getModule(demo.get(module->getID()), module->getID(), module->getLocation());

    {
      position.push_back(JModulePosition_t(module->getID(), getPosition(buffer, *module)));
    }
    if (module->getFloor() != 0) {
      rotation.push_back(JModuleRotation_t(module->getID(), getRotation(buffer, *module)));
    }

    for (JModule::const_iterator pmt = module->begin(); pmt != module->end(); ++pmt) {
      calibration.push_back(JPMTCalibration_t(pmt->getID(), getCalibration(JCalibration(), *pmt)));
    }
  }


  json js;

  js[Result_t] = OK_t;

  js[Data_t][0][PMTT0s_t]       = json(calibration);
  js[Data_t][0][DOMPositions_t] = json(position);
  js[Data_t][0][DOMRotations_t] = json(rotation);

  if (outputFile != "") {

    ostream* out = open<ostream>(outputFile.c_str());

    *out << setw(2) << setprecision(8);
    *out << js;

    close(out);
  }
}