JEditDetector.cc File Reference

Auxiliary program to modify detector calibration. More...

#include <iostream>
#include <sstream>
#include <string>
#include <vector>
#include "TRandom3.h"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorHeader.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JPMTIdentifier.hh"
#include "JGeometry3D/JVector3D.hh"
#include "JTools/JConstants.hh"
#include "JLang/JException.hh"
#include "JSupport/JMeta.hh"
#include "Jeep/JeepToolkit.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 modify detector calibration.

Syntax:

    -M "<module identifier> (set|add|sub|randset|randadd|randsub) x0 [x1 [x2]]"
    -S "<string identifier> (set|add|sub|randset|randadd|randsub) x0 [x1 [x2]]"
    -M "<module identifier> (rot|randrot) phi"
    -S "<string identifier> (rot|randrot) phi"
    -M "<module identifier> (mul|randmul) factor"
    -S "<string identifier> (mul|randmul) factor"
    -P "<PMT identifier>    (set|reset)   PMT_STATUS"
    -@ "<key>=<value>[;<key>=<value>"

Options -M and -S refer to a module and a string, respectively.
If module identifier or string identifier is -1, the action takes place on all modules or strings in the detector, respectively.
The number of values apply to position or time calibration in the following way:

1. time calibration (t = x0)
2. position calibration (x = x0, y = x1, z = 0)
3. position calibration (x = x0, y = x1, z = x2)

The rotation is anti-clockwise around the z-axis.
The multiplication applies to the z-coordinates.

Note that the angles are defined in radians.
Option -@ refers to the header information.
The list of possible keys can be obtained using JPrintDetector.cc with option -O header.

Multiple options -M, -S or -@ will be processed in order of appearance.

Note that if the output file name is the same as the input file name, the original file will be overwritten.

Author

mdejong

Definition in file JEditDetector.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 360 of file JEditDetector.cc.

{
  using namespace std;
  using namespace JPP;

  string               inputFile;
  string               outputFile;
  vector<JModifier>    mod;
  vector<JModifier>    str;
  vector<string>       hdr;
  vector<JPMTModifier> pmt;
  int                  debug;
  
  try { 

    JParser<> zap("Auxiliary program to modify detector.");

    zap['a'] = make_field(inputFile);
    zap['M'] = make_field(mod)           = JPARSER::initialised();
    zap['S'] = make_field(str)           = JPARSER::initialised();
    zap['@'] = make_field(hdr)           = JPARSER::initialised();
    zap['P'] = make_field(pmt)           = JPARSER::initialised();
    zap['o'] = make_field(outputFile);
    zap['d'] = make_field(debug)         = 2;

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

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


  gRandom->SetSeed(0);


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


  if (!hdr.empty()) {

    JProperties helper = detector.getProperties();
    
    for (vector<string>::const_iterator i = hdr.begin(); i != hdr.end(); ++i) {
      
      istringstream is(*i);
      
      is >> helper;
    }
  }


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

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

      if (module->getID() == i->id || i->id == JModifier::WILD_CARD ){  

        DEBUG("Modifier"                                                                              << ' '
              << "(" << setw(3) << module->getString() << "," << setw(2) << module->getFloor() << ")" << ' '
              << setw(8) << module->getID()                                                           << ' '
              << "action" << ' ' << i->action << endl);
        
        if (!i->apply(*module)) {
          ERROR("No valid action: " << *i << endl);
        }
      }
    }
  }


  for (vector<JModifier>::const_iterator i = str.begin(); i != str.end(); ++i) {
  
    for (JDetector::iterator module = detector.begin(); module != detector.end(); ++module) {

      if (module->getString() == i->id || i->id == JModifier::WILD_CARD) {

        DEBUG("Modifier"                                                                              << ' '
              << "(" << setw(3) << module->getString() << "," << setw(2) << module->getFloor() << ")" << ' '
              << setw(8) << module->getID()                                                           << ' '
              << "action" << ' ' << i->action << endl);
        
        if (!i->apply(*module)) {
          ERROR("No valid action: " << *i << endl);
        }
      }
    }
  }


  for (vector<JPMTModifier>::const_iterator i = pmt.begin(); i != pmt.end(); ++i) {
  
    for (JDetector::iterator module = detector.begin(); module != detector.end(); ++module) {

      if (module->getID() == i->getModuleID() || i->getModuleID() == JModifier::WILD_CARD) {

        DEBUG("PMT modifier"                                                                      << ' '
              << "(" << setw(8) << module->getID() << "," << setw(2) << i->getPMTAddress() << ")" << ' '
              << "action" << ' ' << i->action << ' ' 
              << "value"  << ' ' << i->value  << endl);
        
        if (i->getPMTAddress() <  0                         ||
            i->getPMTAddress() >= getNumberOfPMTs(*module)  ||
            !i->apply(module->getPMT(i->getPMTAddress()))) {
          ERROR("No valid action: " << *i << endl);
        }
      }
    }
  }

  try {
    store(outputFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }
}