To perform simulation of parton showers and/or hadronization one have to process the PEGASUS output in Les Houches Event format (*.lhe) with Monte-Carlo event generator CASCADE.

NOTE: to reconstruct the TMD initial state parton showers, the TMD gluon density from CCFM or PB families should be chosen.

To produce an Les Houches Event file with PEGASUS please follow the manual. Shortly, one should choose the process by double-clicking, limit the phase space in order to record only necessary events, and check the “Generate Les Houches Event file” box in the lower left corner. One then starts the calculation with the Start button (or F5). The produced file will be created in the PEGASUS folder and is ready to use with the CASCADE. Here we describe the procedure for 2.4 version.

As a basic file one can take the steer_pp-LHEin file, provided with the distributive. Just check, that the following lines are set as needed:

NEVENT -1
'IPRO'     1       0        -1      ! read in LHE file 

Other parameters define, in particular, whether you turn the fragmentation and showers on or off. You can set them as you need. Importantly, one has to set the IGLU parameter corresponding to the TMD gluon density, used in the LHE file calculation. The TMD number can be found in the LHE file, it corresponds to the TMDLib numbering scheme. PEGASUS uses the following of them (see Overview). The ICCF parameter should be defined according to the used TMD gluon distribution.

The path to the LHE file to be processes is defined by CLHE.

NOTE: ILHE parameter should be set to be 0 for the off-shell partons.

Typically one produces a HepMC output file with CASCADE. One defines the name of the file with

export HEPMCOUT=name_of_the_output_file

Then one starts the CASCADE with

export CASEED=some_random_number
./cascade < steer_pp-LHEin

The produced HepMC file can be then processed to form necessary cross sections. Note that the weighted events from PEGASUS are normalized by CASCADE, so one should multiply the weights by the total cross section during processing. The total cross section can be found in the CASCADE output or in the LHE file.