--
GavinKirby - 2010-09-29
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events.
The code consists primarily of three files, in addition to the .params file which specifies the run parameters. These files are:
- src/tth.f -- the main Fortran source code,
- inc/tth.inc -- contains variable declarations
- inc/tth_inl.inc -- containing inline functions
- pycomm.inc --
The shell script test/linktth is used to compile the code.
A run may be initiated in the test area by using the command
./tthsf < [filename].params > [output logfile]
for a given parameters file. The code is set up to dump unsatisfactory events in the output logfile for debugging purposes.
Params Files
Run information is specified in a .params file. The top line of a .params file should look like:
14000D0 5000 10 500 p p
In this case, 14000 is the centre of mass energy (in
GeV), 5000 is the number of events to generate, 10 is the number of events to printout at the start of the run, 500 is the subsequent printout interval, and "p p" is the initial state appropriate to the LHC.
Pythia Code
tth.f
tth.inc
This file contains variable declarations for the main program and each of the subroutines in tth.f: error counters, event particle counters and particle indices, as well as logical (Boolean) variables used in the main Fortran code.
tth_inl.inc
Changes from VBF
The code was initially designed to generate vector boson (WW) fusion events, and it has been adapted in stages to generate ttH instead, with some additional modifications. In particular, the histogramming aspect of the initial VBF code (which used the Pythia data to prepare an hbook file) has been disabled.
Particle Number Checking
Particle Parentage Checking