createPMTMap.cc File Reference

Auxiliary program to create map for reference PMTs. More...

#include <string>
#include <iostream>
#include <iomanip>
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleRouter.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
#include "JTrigger/JTriggerToolkit.hh"

Go to the source code of this file.

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

Detailed Description

Auxiliary program to create map for reference PMTs.

the option

    -T "s" 

or

    -T{"s1","s2"} 

has to be used with a vector of strings. This vector contains either the list of reference PMT identities or the keyword odd. If the keyword odd is used, all the PMTs with an odd identity are listed in the output file. The parity is calculated as following:

        (floor id - 1)*31 + PMT id

Author

acreusot

Definition in file createPMTMap.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 37 of file createPMTMap.cc.

{
  using namespace std;
  using namespace JPP;
  //using namespace KM3NETDAQ;

  string         detectorFile;
  string         outputFile;
  vector<string> mapType;
  int            debug;

  try {
    JParser<> zap("Auxiliary program to create PMT map files.");
    zap['a'] = make_field(detectorFile);
    zap['o'] = make_field(outputFile);
    zap['T'] = make_field(mapType);
    zap['d'] = make_field(debug)          = 1;
    zap(argc, argv);
  }

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

  fstream outstr(outputFile.c_str(), ofstream::out);
  JDetector detector;
  try {
    load(detectorFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }
  for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
    const int floorId = module->getFloor();
    const int domId = module->getID();
    for (unsigned int t = 0; t < mapType.size(); ++t) {
      if (mapType[t] == "odd") {
        for (unsigned int pmt = 0; pmt < NUMBER_OF_PMTS; ++pmt) {
          const int id = (floorId - 1)*NUMBER_OF_PMTS + pmt;
          if (id % 2 != 0) {
            outstr << domId << ' ' << pmt << endl;
          }
        }
      } else {
        outstr << domId << ' ' << mapType[t] << endl;
      } 
    } 
  }
  outstr.close();

}