JPrintDetector.cc File Reference

#include <string>
#include <iostream>
#include <iomanip>
#include <set>
#include <iterator>
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JTools/JConstants.hh"
#include "JTools/JRange.hh"
#include "JGeometry3D/JGeometry3DToolkit.hh"
#include "JGeometry3D/JCylinder3D.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 print detector file in human friendly format.

Author

mdejong

Definition in file JPrintDetector.cc.

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 50 of file JPrintDetector.cc.

{
  using namespace std;
 
  string         detectorFile;
  string         option;
  int            debug;
 
  try {
 
    JParser<> zap("Auxiliary program to print detector file in human friendly format.");
 
    zap['a'] = make_field(detectorFile);
    zap['O'] = make_field(option)            = 
      default_t,
      header_t,
      pmts_t,
      modules_t,
      version_t,
      identifier_t,
      can_t,
      center_t,
      summary_t,
      VERSION_t,
      IDENTIFIER_t,
      CAN_t,
      CENTER_t,
      SUMMARY_t;
    zap['d'] = make_field(debug)             = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  using namespace JPP;
  
  JDetector detector;
 
  try {
    load(detectorFile, detector);
  }
  catch(const JException& error) {
    FATAL(error);
  }
 
 
  if        (option == default_t) {
 
    cout << detector << endl;
 
  } else if (option == header_t) {
 
    cout << detector.getProperties() << endl;
 
  } else if (option == modules_t) {
      
    for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
 
      double t0 = 0.0;
 
      for (JModule::const_iterator pmt = module->begin(); pmt != module->end(); ++pmt) {
        t0 += pmt->getT0();
      }
      
      t0 /= module->size();
      
      cout << "Module";
      cout << ' ' << noshowpos << setw(8) << right << module->getID();
      cout << ' ' << noshowpos << setw(3) << right << module->getString();
      cout << ' ' << noshowpos << setw(2) << right << module->getFloor();
      cout << ' ' << FIXED(7,2) << module->getX();
      cout << ' ' << FIXED(7,2) << module->getY();
      cout << ' ' << FIXED(7,2) << module->getZ();
      cout << ' ' << FIXED(8,2) << t0;
      cout << endl;
    }
    
  } else if (option == pmts_t) {
          
    for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
      for (JModule::const_iterator pmt = module->begin(); pmt != module->end(); ++pmt) {
        
        cout << "PMT";
        cout << ' ' << noshowpos << setw(8) << right << pmt->getID();
        cout << ' ' << FIXED(7,2) << pmt->getX();
        cout << ' ' << FIXED(7,2) << pmt->getY();
        cout << ' ' << FIXED(7,2) << pmt->getZ();
        cout << ' ' << FIXED(6,3) << pmt->getDX();
        cout << ' ' << FIXED(6,3) << pmt->getDY();
        cout << ' ' << FIXED(6,3) << pmt->getDZ();
        cout << ' ' << FIXED(8,2) << pmt->getT0();
        cout << ' ' << noshowpos << setw(8) << right << pmt->getStatus();
        cout << endl;
      }
    }
  
  } else if (option == version_t || 
             option == VERSION_t) {
 
    if (option == version_t) {
 
      cout << detector.getVersion() << endl;
 
    } else {
 
      cout << "set_variable DETECTOR_VERSION " << detector.getVersion() << ";" << endl;
    }
 
  } else if (option == identifier_t ||
             option == IDENTIFIER_t) {
 
    if (option == identifier_t) {
 
      cout << "Detector " << detector.getID() << endl;
 
    } else {
 
      cout << "set_variable DETECTOR_ID      " << detector.getID() << ";" << endl;
    }
 
  } else if (option == can_t ||
             option == CAN_t) {
 
    const JCylinder3D cylinder(detector.begin(), detector.end());
 
    const double      V = (cylinder.getZmax() - cylinder.getZmin()) * PI * cylinder.getRadius() * cylinder.getRadius(); 
 
    if (option == can_t) {
      
      cout << "X      = " << FIXED(7,1)       << cylinder.getX()      << endl;
      cout << "Y      = " << FIXED(7,1)       << cylinder.getY()      << endl;
      cout << "Zmin   = " << FIXED(7,1)       << cylinder.getZmin()   << endl;
      cout << "Zmax   = " << FIXED(7,1)       << cylinder.getZmax()   << endl;
      cout << "Radius = " << FIXED(7,1)       << cylinder.getRadius() << endl;
      cout << "Depth  = " << FIXED(7,1)       << detector.getUTMZ()   << endl;
      cout << "Volume = " << SCIENTIFIC(12,3) << V                    << endl;
 
    } else {
 
      cout << "set_variable CAN_X_M      " << FIXED(7,1)       << cylinder.getX()           << ";" << endl;
      cout << "set_variable CAN_Y_M      " << FIXED(7,1)       << cylinder.getY()           << ";" << endl;
      cout << "set_variable CAN_ZMIN_M   " << FIXED(7,1)       << cylinder.getZmin()        << ";" << endl;
      cout << "set_variable CAN_ZMAX_M   " << FIXED(7,1)       << cylinder.getZmax()        << ";" << endl;
      cout << "set_variable CAN_RADIUS_M " << FIXED(7,1)       << cylinder.getRadius()      << ";" << endl;
      cout << "set_variable CAN_DEPTH_M  " << FIXED(7,1)       << detector.getUTMZ()        << ";" << endl;
      cout << "set_variable CAN_VOLUME_M3" << SCIENTIFIC(12,3) << V                         << ";" << endl;
    }
  
  } else if (option == center_t || option == CENTER_t) {
      
    const JCenter3D center(detector.begin(), detector.end());
      
    if (option == center_t) {
 
      cout << "center = ";
      cout << showpos << FIXED(8,3) << center.getX() << ' ';
      cout << showpos << FIXED(8,3) << center.getY() << ' ';
      cout << showpos << FIXED(8,3) << center.getZ() << endl;
 
    } else {
 
      cout << "set_variable CENTER_X_M " << FIXED(7,1) << center.getX() << ";" << endl;
      cout << "set_variable CENTER_Y_M " << FIXED(7,1) << center.getY() << ";" << endl;
      cout << "set_variable CENTER_Z_M " << FIXED(7,1) << center.getZ() << ";" << endl;
    }
    
  } else if (option == summary_t ||
             option == SUMMARY_t) {
      
    const int numberOfStrings    = getNumberOfStrings(detector);
    const int numberOfFloors     = getNumberOfFloors (detector);
    const int numberOfModules    = getNumberOfModules(detector);
 
    typedef JTOOLS::JRange<int>   JRange_t;
 
    const JRange_t string(detector.begin(), detector.end(), &JModule::getString);
    const JRange_t floor (detector.begin(), detector.end(), &JModule::getFloor);
 
    set<int> strings;
 
    for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
      strings.insert(module->getString());
    }
 
    if (option == summary_t) {
        
      cout << "Number of strings = " << setw(4) << numberOfStrings << endl;
      cout << "Number of floors  = " << setw(4) << numberOfFloors  << endl;
      cout << "Number of modules = " << setw(4) << numberOfModules << endl;
      cout << "First string      = " << setw(4) << string.first    << endl;
      cout << "Last  string      = " << setw(4) << string.second   << endl;
      cout << "First floor       = " << setw(4) << floor .first    << endl;
      cout << "Last  floor       = " << setw(4) << floor .second   << endl;
 
    } else {
      
      cout << "set_variable NUMBER_OF_STRINGS " << setw(4) << numberOfStrings << ";" << endl;
      cout << "set_variable NUMBER_OF_FLOORS  " << setw(4) << numberOfFloors  << ";" << endl;
      cout << "set_variable NUMBER_OF_MODULES " << setw(4) << numberOfModules << ";" << endl;
      cout << "set_variable FIRST_STRING      " << setw(4) << string.first    << ";" << endl;
      cout << "set_variable LAST_STRING       " << setw(4) << string.second   << ";" << endl;
      cout << "set_variable FIRST_FLOOR       " << setw(4) << floor .first    << ";" << endl;
      cout << "set_variable LAST_FLOOR        " << setw(4) << floor .second   << ";" << endl;
      cout << "set_array    STRINGS           ";
      copy(strings.begin(), strings.end(), ostream_iterator<int>(cout, " "));
      cout << endl;
    }
  }
}