main125
Back to index.
// main125.cc is a part of the PYTHIA event generator.
// Copyright (C) 2025 Torbjorn Sjostrand.
// PYTHIA is licenced under the GNU GPL v2 or later, see COPYING for details.
// Please respect the MCnet Guidelines, see GUIDELINES for details.
// Keywords:
//            Basic usage
//            LHE file
// This is a simple test program. It shows how PYTHIA 8 can write
// a Les Houches Event File v. 3.0 with hadron-level events, based on
// its process-level input events. Event weights are carried along.
// Note: HepMC is more commonly used for hadron-level output than LHEF.
#include "Pythia8/Pythia.h"
using namespace Pythia8;
//==========================================================================
int main() {
  // Set up for external input of LHEF 3.0 process-level events.
  Pythia pythia;
  pythia.readString("Beams:frameType = 4");
  pythia.readString("Beams:LHEF = wbj_lhef3.lhe");
  // If Pythia fails to initialize, exit with error.
  if (!pythia.init()) return 1;
  // Create and open file for LHEF 3.0 hadron-level output.
  LHEF3FromPythia8 myLHEF3(&pythia.event, &pythia.info);
  myLHEF3.openLHEF("weakbosons.lhe");
  // Write out initialization info on the file.
  myLHEF3.setInit();
  // Event generation loop.
  for (int iEvent = 0; iEvent < 10; ++iEvent) {
    // Generate next event.
    if (!pythia.next()) {
      if( pythia.info.atEndOfFile() ) break;
      else continue;
    }
    // Store and write event info.
    myLHEF3.setEvent();
  } // End loop over events to generate.
  // Statistics: full printout.
  pythia.stat();
  // Write endtag. Overwrite initialization info with new cross sections.
  myLHEF3.closeLHEF(true);
  // Done.
  return 0;
}