JPMTAngularAcceptance.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <limits>
 
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
 
/**
 * \author mdejong
 */
 
/**
 * Auxiliary program to convert PMT QE(angle) data.
 */
int main(int argc, char **argv)
{
  using namespace std;
 
  string  inputFile;
  int     debug;
 
  try {
 
    JParser<> zap("Auxiliary program to convert PMT QE(angle) data.");
 
    zap['f'] = make_field(inputFile);
    zap['d'] = make_field(debug)          = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  
  const string DEG("deg");
 
  
  if (inputFile != "") {
 
    ifstream in(inputFile.c_str());
 
    while (in.peek() == '#') {
      in.ignore(numeric_limits<streamsize>::max(), '\n');
    }
 
    // header
    
    vector<double> energy;
    vector<double> wavelength;
    
    string buffer, key;
 
    if (getline(in,buffer)) {
          
      istringstream is(buffer);
 
      is >> key;
      
      for (double x; is >> x; ) {
        energy.push_back(x);
      }
    }
 
    if (getline(in,buffer)) {
          
      istringstream is(buffer);
 
      is >> key;
      
      for (double x; is >> x; ) {
        wavelength.push_back(x);
      }
    }
 
    //in.ignore(numeric_limits<streamsize>::max(), '\n');
 
 
    // data
    
    while (getline(in,buffer)) {
 
      istringstream is(buffer);
 
      // first column
      /*
        is >> key;
 
        const size_t pos = key.find(DEG);
 
        if (pos != string::npos) {
        key.replace(pos, DEG.size(), "");
        }
 
        double angle;
      
        istringstream(key) >> angle;
      */
      double ct;
 
      is >> ct;
 
      // following columns
      
      size_t i = 0;
      
      for (double y ; is >> y; ++i) {
 
        if (i != energy.size()) {
          cout << "      (*this)"
               << noshowpos
            //<< "[" << FIXED(5,1) << angle     << "]"
            //<< "[" << FIXED(3,1) << energy[i] << "]"
               << "[" << FIXED(5,2) << -ct           << "]"
               << "[" << FIXED(3,1) << wavelength[i] << "]"
               << " = "
               << showpos
               << FIXED(9,5) << y
               << ";" << endl;
        } else {
          //FATAL("Inconsistent data " << i << " != " << energy.size() << endl);
          FATAL("Inconsistent data " << i << " != " << wavelength.size() << endl);
        }
      }
 
      if (i != energy.size()) {
        //FATAL("Inconsistent data " << i << " != " << energy.size() << endl);
        FATAL("Inconsistent data " << i << " != " << wavelength.size() << endl);
      }
    }
 
    in.close();
  }
}