getAAnetFiles.cc

Go to the documentation of this file.

#include <string>
#include <iostream>
#include <iomanip>
#include <map>

#include "km3net-dataformat/offline/Head.hh"
#include "JLang/JVectorize.hh"
#include "JAAnet/JHead.hh"
#include "JAAnet/JHeadToolkit.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JSupport.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

namespace {

  using JAANET::JHead;
  
  struct JWeight_t :
    public std::map<std::string, JHead>
  {
    /**
     * Default constructor.
     */
    JWeight_t()
    {
      using namespace JPP;

      (*this)["MUPAGE"]      = getMUPAGEHeader();
      (*this)["GENHEN"]      = getGenhenHeader();
      (*this)["GENIE"]       = getGenieHeader();
      (*this)["SIRENE"]      = getSireneHeader();
      (*this)["KM3"]         = getKM3Header();
      (*this)["KM3SIM"]      = getKM3SimHeader();
      (*this)["CORSIKA_P"]   = getCorsikaHeader( 0,14);   // proton is single nucleon 
      (*this)["CORSIKA_HE"]  = getCorsikaHeader( 4, 2);
      (*this)["CORSIKA_C"]   = getCorsikaHeader(12, 6);
      (*this)["CORSIKA_0"]   = getCorsikaHeader(16, 8);
      (*this)["CORSIKA_FE"]  = getCorsikaHeader(56,26);
      (*this)["DAQ"]         = getDAQHeader();
    }
  };
}


/**
 * \file
 * Auxiliary program to select Monte Carlo files according generator.
 *
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
  
  JWeight_t weight;

  JMultipleFileScanner_t inputFile;
  string         header;
  bool           option;
  int            debug;

  try { 

    JParser<> zap("Auxiliary program to select Monte Carlo files according generator.");
    
    zap['f'] = make_field(inputFile);
    zap['H'] = make_field(header)          = get_keys(weight);
    zap['O'] = make_field(option);
    zap['d'] = make_field(debug)           = 1;

    zap['H'] = JPARSER::not_initialised();

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

  DEBUG(weight[header]);
  
  cout << getAAnetFiles(inputFile, weight[header], option);
}