HELAC-PHEGAS is a multi-purpose, tree-level  event generator which is the only existing implementation of the algorithm based on Dyson-Schwinger equations. It can be used to efficiently obtain helicity amplitudes and total cross sections  for arbitrary multiparticle processes in the Standard Model. The program can generate all processes with 10 or more final state particles with full off-shell and finite width effects taking into account naturally both, spin and color correlations. All electroweak vertices in both Feynman and unitary gauges have been included. The program incorporates the possibility to use extended numerical precision. It is also possible to get numerical predictions for amplitudes for arbitrary processes for phase-space points provided by the user.

For the phase space generation the PHEGAS program is used which has been the first implementation of a complete automated algorithm of multi-channel phase-space mappings for arbitrary number of external particles. It also uses a self-adaptation procedure that reshapes the generated phase-space density, including also the PARNI  algorithm for the efficient integration over parton distribution functions.

The program generates parton-level events with all necessary information, in the most recent Les Houches Accord format, for the study of any process at hadron and lepton colliders. Moreover, the emission of secondary partons  as well as the translation of the emerging partons into primordial hadrons is performed by interfacing HELAC-PHEGAS to PYTHIA.  In addition, the so called MLM matching procedure is used to combine fixed order tree level matrix element calculations with parton showers in order to avoid double counting of jet configurations.

The program has been built with future extensions in mind. The design includes all mechanisms for expanding the program with new capabilities without having to make major changes to the main structure. Because of this it is straightforward to add new physics effects by including new models and new couplings. Moreover, a tested upgrade with high color charge processes treated via Monte Carlo summation over color is available.

The code is written in Fortran 95 and has been developed using the
following compilers under Linux: Lahey Fujitsu, Intel Fortran, GNU gfortran and g95. The stand-alone codes do not need any external library. However, LHAPDF is supported as an option.

The use of the code is very easy. The user has just to run a shell script that produces the file with the un-weighted events in the Les Houches format. Parameters, like physical constants, number of un-weighted events or phase-space cuts can be entered, in a very transparent way, in the user.inp file. Default values for all of them exist in the default.inp file. The script also supports running under the LSF batch system, so that each sub-process runs on a different computer on the LSF cluster and it automatically collects all information, without any intervention from the user. It has been successfully tested at lxplus at CERN.

:: top ::


Costas G. Papadopoulos  
:: top ::


  lahey 95
  Intel Fortran

:: top ::

Software Download

   HELAC-PHEGAS  Curent Version 2.0.1 (29 January 2015): helac-phegas.2.0.1.tgz

Version 2.0.0 (8 October 2014): 

Version 1.2.2 (13 December 2010):  helac-phegas.1.2.2.tgz

Version 1.2.1 (7 July 2010):  helac-phegas.1.2.1.tgz

Version 1.2.0 (6 May 2009):  helac-phegas.1.2.0.tar.gz

Version 1.1.0 (18 January 2008):  tar file    

Version 1.0.0 (17 October 2007):
  tar file

Version 0.1.0 (24 July 2006):  tar file

:: top ::


The following papers should be cited while using HELAC-PHEGAS package

Helac-Phegas: A Generator for all parton level processes
A. Cafarella, C. G. Papadopoulos, M. Worek
Comput. Phys. Commun. 180 (2009) 1941

: A Phase space generator for automatic cross-section computation
C. G. Papadopoulos
Comput. Phys. Commun. 137 (2001) 247

HELAC: A Package to compute electroweak helicity amplitudes
A. Kanaki, C. G. Papadopoulos
Comput. Phys. Commun. 132 (2000) 306

:: top ::

News & Developments & Bug Report & Validation

Changes in HELAC-PHEGAS Version 2.0.1

     One point phase space option corrected. Example for gg-> ttbb process included (see the user.inp file and the mom file, where the phase
        space point is given) 

Changes in HELAC-PHEGAS Version 2.0.0

     Simple way to select the order of the strong and electroweak coupling constants for a given process has been implemented. As an example
        the top anti-top production in di-lepton channel with all contributions of the order of alpha_s = 2 alpha_ew = 4  is given in the user.inp file

Changes in HELAC-PHEGAS Version 1.2.2

     Summation mode for high number of subprocesses with LHA files corrected
        Tested on pp -> gamma + 5jets process 

Changes in HELAC-PHEGAS Version 1.2.1

     Optimization for the phase space corrected for the massive cases
        - on-shell production of top, W and Z 

     Bugs fixed in the summation mode for high number of subprocesses with unweighting
        - thanks to Tim Adye !
     pT_miss cut added for TeVatron and LHC case as well as for e+e- collisions
     Some minor changes in the output format

Changes in HELAC-PHEGAS Version 1.2.0

     lhaglue.f file has been removed from the package to allow HELAC to run with the newest
        LHAPDF lhaglue interface within the lhapdf-5.7.0
     Several bugs fixed
     New version of PARNI included
Changes in HELAC-PHEGAS Version 1.1.0

     Scripts have been added and some files have been modified to allow running HELAC with
        the LSF system

     A bug in main_mc.f has been fixed
     Shell changed from sh to bash
     "special" directory added to allow removal of some vertices in particular situations
     Many unnecessary files have been removed
     README file has been updated
     Some minor changes in the output format

Validation for current HELAC-PHEGAS version (helac-phegas.1.2.1.tgz)

     tT production @ TeVatron (CTEQ6L1 from LHAPDF version 5.7.0, scale = mtop)
     tT production @ LHC (CTEQ6L1 from LHAPDF version 5.7.0, scale = mtop)
     tT production @ e+e- (scale = mtop)
     e+ve + 1j @ LHC (CTEQ6L1 stand-alone, scale = mZ)
     e+e- + 1j  @ LHC (CTEQ6L1 stand-alone, scale = mZ)

:: top ::

Contact us

If you have a question, comment, suggestion or bug report, please e-mail us at:



:: top ::

Last modified by Malgorzata Worek
Thu, 29 January, 2015