Difference: PythiattH (1 vs. 11)

Revision 112012-09-27 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Changed:
<
<
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events. The relevant code can be found in the Glasgow ATLAS subversion repository under Pythiatth. Note that in order to run it you will need to have the correct Pythia version in a directory at the same level as the tth directory (the linktth_64.sh script specifies the path to the Pythia version to use). There is a compiled and ready to use Pythia version (6.4) in /afs/phas.gla.ac.uk/user/g/gkirby/public/pythia6426. Note that the older linktth.sh script is for use with pythia6222, which is no longer used. Information about studies that have been carried out using this code may be found at PythiattHWork.
>
>
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events. The relevant code can be found in the Glasgow ATLAS subversion repository under Pythiatth. Note that in order to run it you will need to have the correct Pythia version in a directory at the same level as the tth directory (the linktth_64.sh script specifies the path to the Pythia version to use). There is a compiled and ready to use Pythia version (6.4) in /afs/phas.gla.ac.uk/user/g/gkirby/public/pythia6426. Note that the older linktth.sh script is for use with pythia6222, which is no longer used. Information about studies that have been carried out using this code may be found at PythiattHWork. The "standard" version of this code generates semileptonic ttH events; a version to produce diletponic events is in progress, and will be included as a branch in the SVN repository.
 

Revision 102012-07-26 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Changed:
<
<
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events. The relevant code can be found in the Glasgow ATLAS subversion repository under Pythiatth. Note that in order to run it you will need to have the correct Pythia version in a directory called pythia6222, at the same level as the tth directory (this requirement is due to the linktth.sh script, which specifies the path to the Pythia version to use). This can be found in ~gkirby/Pythia_versions/pythia6222.
>
>
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events. The relevant code can be found in the Glasgow ATLAS subversion repository under Pythiatth. Note that in order to run it you will need to have the correct Pythia version in a directory at the same level as the tth directory (the linktth_64.sh script specifies the path to the Pythia version to use). There is a compiled and ready to use Pythia version (6.4) in /afs/phas.gla.ac.uk/user/g/gkirby/public/pythia6426. Note that the older linktth.sh script is for use with pythia6222, which is no longer used. Information about studies that have been carried out using this code may be found at PythiattHWork.
 
Changed:
<
<
The code consists primarily of three files, in addition to the .params file which specifies the run parameters. These files are:
>
>

Overview

The code can be found in the Pythiatth directory of the Glasgow "atlas" SVN repository.

The code consists primarily of four files, in addition to the .params file which specifies the run parameters. These files are:

 
  • src/tth.f -- the main Fortran source code,
Changed:
<
<
  • inc/tth.inc -- contains variable declarations
  • inc/tth_inl.inc -- containing inline functions
>
>
  • inc/tth.inc -- contains variable declarations for the Fortran code
  • inc/tth_inl.inc -- containing inline functions for the Fortran code
 
  • inc/pycomm.inc -- inline functions for Pythia variables
Changed:
<
<
The shell script test/linktth.sh is used to compile the code. It creates two executable files: tth and tthsf (i.e. using structure functions).
>
>
The shell script test/linktth_64.sh is used to compile the code. It creates an executable file called tth (6.2 also has tthsf, i.e. using structure functions).
  A run may be initiated in the test area by using the command
Changed:
<
<
./tthsf < [filename].params > [output logfile]
>
>
./tth < [filename].params > [output logfile]
 
Changed:
<
<
for a given parameters file. The code is set up to dump unsatisfactory events in the output (logfile) for debugging purposes.
>
>
for a given parameters file. The code is set up to dump unsatisfactory events and other select information in the output (logfile) for debugging purposes.
 

Params Files

Line: 27 to 31
  14000D0 5000 10 500 p p
Changed:
<
<
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 (e.g. p pbar would be used for Tevatron runs).
>
>
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 (e.g. p pbar would be used for Tevatron runs).
 
Changed:
<
<

Pythia Code

>
>
The .params files also contain switches to control particle decay modes, Q^2 scale parameters, particle mass settings (in particular, PMAS(25,1) controls the Higgs mass). Consult the Pythia 6.4 manual for a full description of Pythia parameters.

Code Structure

 

tth.f

Changed:
<
<
This contains the main function (program higgs) and the subroutines it calls. The main function generates the required number of events (as specified by the .params file) and tests them to ensure that they comply with specific criteria (for debugging/sanity checking purposes). It also creates an output text file containing specific information about the generated events (the Boolean variable writeout is used to control which events are considered acceptable for this output file). The format of this output file is such that it can be parsed to create an ntuple using the relevant tools.
>
>
This contains the main function ( program higgs) and the subroutines it calls. The main function generates the required number of events (as specified by the .params file) and tests them to ensure that they comply with specific criteria (for debugging/sanity checking purposes). It also creates an output text file called PythiaOutput.dat containing specific information about the generated events, such as neural net variables or event shape variables (the Boolean variable writeout is used to control which events are considered acceptable for this output file). The format of this output file is such that it can be parsed to create an ntuple using the relevant tools.
 

tth.inc

Changed:
<
<
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.
>
>
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. New variables must be declared in the main body of the file and then listed under the appropriate block below. All new variables should be zeroed before being used in the main Fortran code.
 

tth_inl.inc

Changed:
<
<
>
>
Contains definitions of inline functions for the tth code.
 

pycomm_inl.inc

Changed:
<
<
This file contains inline functions relating to Pythia internal variables (e.g. the energy and momenta of particles). It assigns these values to a set of variables (declared in the file) which can then be used in the main function (e.g. for writeout of a text file containing event information suitable for building an ntuple).
>
>
This file contains inline functions relating to Pythia internal variables (e.g. the energy and momenta of particles). It assigns these values to a set of variables (declared in the file) which can then be used in the main function (e.g. for writeout of a text file containing event information suitable for building an ntuple). Thus, pt_r(particle index) will return the pt of the relevant particle.
 

Changes from VBF

Line: 52 to 56
  The code was initially designed to generate vector boson (WW) fusion events, and it has been adapted in stages to generate ttH events 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, with a view to amending the code to produce a text output file containing the necessary information (in the appropriate format) to prepare root ntuples for histogramming and analysis.
Deleted:
<
<
 

Particle Number Checking

Changed:
<
<
The code uses a number of independent checks to ensure that the events that Pythia is generating have the anticipated structure. In particular, it is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar. The code has been developed to focus specifically on H-->b/bbar, but can be modified for studies involving other Higgs decay modes.
>
>
The code uses a number of independent checks to ensure that the events that Pythia is generating have the anticipated structure. In particular, it is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar. The code has been developed to focus specifically on semileptonic ttH with H-->b/bbar, but can be modified for studies involving other Higgs production and decay modes.
 

Particle Parentage Checking

Changed:
<
<
It is also important to check that the particles in each event correspond to the expected pattern of decays (i.e. parent/daughter particles). For instance, one expects the lepton in a semileptonic ttH process to appear as a result of the weak decay of a W (or from H->ll) and this condition can be imposed on events, so together with the independent counting of particles it is possible to impose the condition (for example) that any acceptable event must contain exactly one lepton with a W parent (of appropriate charge).
>
>
It is also important to check that the particles in each event correspond to the expected pattern of decays (i.e. parent/daughter particles). For instance, one expects the lepton in a semileptonic ttH process to appear as a result of the weak decay of a W (or from H->ll) and this condition can be imposed on events, so together with the independent counting of particles it is possible to impose the condition (for example) that any acceptable event must contain exactly one lepton with a W parent (of appropriate charge).

Revision 92011-02-17 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Changed:
<
<
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events.
>
>
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events. The relevant code can be found in the Glasgow ATLAS subversion repository under Pythiatth. Note that in order to run it you will need to have the correct Pythia version in a directory called pythia6222, at the same level as the tth directory (this requirement is due to the linktth.sh script, which specifies the path to the Pythia version to use). This can be found in ~gkirby/Pythia_versions/pythia6222.
 
Line: 11 to 11
 
  • src/tth.f -- the main Fortran source code,
  • inc/tth.inc -- contains variable declarations
  • inc/tth_inl.inc -- containing inline functions
Changed:
<
<
  • pycomm.inc --
>
>
  • inc/pycomm.inc -- inline functions for Pythia variables
 
Changed:
<
<
The shell script test/linktth is used to compile the code.
>
>
The shell script test/linktth.sh is used to compile the code. It creates two executable files: tth and tthsf (i.e. using structure functions).
  A run may be initiated in the test area by using the command

./tthsf < [filename].params > [output logfile]

Changed:
<
<
for a given parameters file. The code is set up to dump unsatisfactory events in the output logfile for debugging purposes.
>
>
for a given parameters file. The code is set up to dump unsatisfactory events in the output (logfile) for debugging purposes.
 

Params Files

Line: 33 to 33
 

tth.f

Changed:
<
<
>
>
This contains the main function (program higgs) and the subroutines it calls. The main function generates the required number of events (as specified by the .params file) and tests them to ensure that they comply with specific criteria (for debugging/sanity checking purposes). It also creates an output text file containing specific information about the generated events (the Boolean variable writeout is used to control which events are considered acceptable for this output file). The format of this output file is such that it can be parsed to create an ntuple using the relevant tools.
 

tth.inc

Changed:
<
<
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.
>
>
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

Added:
>
>

pycomm_inl.inc

This file contains inline functions relating to Pythia internal variables (e.g. the energy and momenta of particles). It assigns these values to a set of variables (declared in the file) which can then be used in the main function (e.g. for writeout of a text file containing event information suitable for building an ntuple).

 

Changes from VBF

Line: 54 to 55
 

Particle Number Checking

Changed:
<
<
The code uses a number of independent checks to ensure that the events that Pythia is generating have the appropriate structure. In particular, it is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar.
>
>
The code uses a number of independent checks to ensure that the events that Pythia is generating have the anticipated structure. In particular, it is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar. The code has been developed to focus specifically on H-->b/bbar, but can be modified for studies involving other Higgs decay modes.
 

Particle Parentage Checking

Revision 82010-11-01 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Line: 27 to 27
  14000D0 5000 10 500 p p
Changed:
<
<
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.
>
>
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 (e.g. p pbar would be used for Tevatron runs).
 

Pythia Code

Line: 49 to 49
 

Changes from VBF

Changed:
<
<
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.
>
>
The code was initially designed to generate vector boson (WW) fusion events, and it has been adapted in stages to generate ttH events 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, with a view to amending the code to produce a text output file containing the necessary information (in the appropriate format) to prepare root ntuples for histogramming and analysis.
 

Particle Number Checking

Changed:
<
<
It is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar.
>
>
The code uses a number of independent checks to ensure that the events that Pythia is generating have the appropriate structure. In particular, it is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar.
 

Particle Parentage Checking

\ No newline at end of file
Added:
>
>
It is also important to check that the particles in each event correspond to the expected pattern of decays (i.e. parent/daughter particles). For instance, one expects the lepton in a semileptonic ttH process to appear as a result of the weak decay of a W (or from H->ll) and this condition can be imposed on events, so together with the independent counting of particles it is possible to impose the condition (for example) that any acceptable event must contain exactly one lepton with a W parent (of appropriate charge).
 \ No newline at end of file

Revision 72010-10-22 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Line: 54 to 54
 

Particle Number Checking

Added:
>
>
It is important to check that the numbers of particles of each kind in an event are correct. For semileptonic ttH we expect to find one top, one tbar, one each of Wplus and Wminus (or two, for H->bbar), one lepton and its associated neutrino (or three leptons, for H->leptons), one (or two) each of b and bbar.
 

Particle Parentage Checking

Revision 62010-10-20 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Line: 19 to 19
  ./tthsf < [filename].params > [output logfile]
Changed:
<
<
for a given parameters file.
>
>
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

Line: 37 to 45
 

tth_inl.inc

Deleted:
<
<

Params Files

 
Deleted:
<
<
Run information is specified in a .params file. The top line of a .params file should look like:
 
Changed:
<
<
14000D0 5000 10 500 p p
>
>

Changes from VBF

 
Changed:
<
<
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.
>
>
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

Revision 52010-10-19 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Line: 11 to 11
 
  • src/tth.f -- the main Fortran source code,
  • inc/tth.inc -- contains variable declarations
  • inc/tth_inl.inc -- containing inline functions
Added:
>
>
  • pycomm.inc --
 
Changed:
<
<
The shell script test/linktth is used to compile the code
>
>
The shell script test/linktth is used to compile the code.
  A run may be initiated in the test area by using the command

Revision 42010-10-19 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Line: 12 to 12
 
  • inc/tth.inc -- contains variable declarations
  • inc/tth_inl.inc -- containing inline functions
Changed:
<
<
The script test/linktth is used to compile the code
>
>
The shell script test/linktth is used to compile the code
  A run may be initiated in the test area by using the command
Changed:
<
<
./tthsf < v[number].params > [output logfile]
>
>
./tthsf < [filename].params > [output logfile]
  for a given parameters file.
Line: 28 to 28
 

tth.inc

Added:
>
>
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

\ No newline at end of file
Added:
>
>

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.

Revision 32010-10-13 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Line: 20 to 20
  for a given parameters file .
Added:
>
>

Pythia Code

 

tth.f

Revision 22010-10-13 - GavinKirby

Line: 1 to 1
 
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29
Deleted:
<
<
This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events
 \ No newline at end of file
Added:
>
>
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

The script test/linktth is used to compile the code

A run may be initiated in the test area by using the command

./tthsf < v[number].params > [output logfile]

for a given parameters file .

tth.f

tth.inc

tth_inl.inc

Revision 12010-09-29 - GavinKirby

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="PythiattH"
-- GavinKirby - 2010-09-29

This page documents the use of Pythia to generate ttH (Higgs-strahlung from top/tbar) events

 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback