getPythia.cc

Go to the documentation of this file.

#include <iostream>
#include <iomanip>
 
#include "JSirene/pythia.hh"
#include "JPhysics/JPhysicsToolkit.hh"
#include "JAAnet/JPDB.hh"
#include "JLang/JManip.hh"
 
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
 
/**
 * \file
 * Main program to get equivalent electro-magnetic energy.
 *
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;
 
  int    type;
  double E_GeV;
  int    debug;
 
  try {
 
    JParser<> zap("Main program to get equivalent electro-magnetic energy.");
 
    zap['P'] = make_field(type,     "PDG particle type");
    zap['E'] = make_field(E_GeV,    "energy [GeV]");
    zap['d'] = make_field(debug)             = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
  double E = E_GeV;
  
  try {
    E = getKineticEnergy(E, JPDB::getInstance().getPDG(type).mass);
  }
  catch(const exception& error) {
    ERROR(error.what() << endl);
  }
  
  E = pythia(type, E);
 
  cout << SCIENTIFIC(12,3) << E << endl;
}