Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] adxl_v2_0.tar.gz(2787 Kbytes)
Manuscript Title: THERMINATOR 2: THERMal heavy IoN generATOR 2
Authors: Mikolaj Chojnacki, Adam Kisiel, Wojciech Florkowski, Wojciech Broniowski
Program title: THERMINATOR 2
Catalogue identifier: ADXL_v2_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 183(2012)746
Programming language: C++ with the CERN ROOT libraries, BASH shell.
Computer: Any with a C++ compiler and the CERN ROOT environment, ver. 5.26 or later, tested with Intel Core2 Duo CPU E8400 @ 3 GHz, 4 GB RAM.
Operating system:
Linux Ubuntu 10.10 x64 (gcc 4.4.5) ROOT 5.26
Linux Ubuntu 11.04 x64 (gcc Ubuntu/Linaro 4.5.2-8ubuntu4) ROOT 5.30/00 (compiled from source)
Linux CentOS 5.2 (gcc Red Hat 4.1.2-42) ROOT 5.30/00 (compiled from source)
Mac OS X 10.6.8 (i686-apple-darwin10-g++-4.2.1) ROOT 5.30/00 (for Mac OS X 10.6 x86-64 with gcc 4.2.1)
cygwin-1.7.9-1 (gcc gcc4-g++-4.3.4-4) ROOT 5.30/00 (for cygwin gcc 4.3).
  • 30 MB therm2 events
  • 150 MB therm2 femto
Keywords: Relativistic heavy-ion collisions, Statistical hadronization, Monte Carlo event generator, LHC, RHIC, SPS, FAIR, NICA.
PACS: 25.75.-q, 25.75.Dw, 25.75.Ld.
Classification: 11.2.

External routines: CERN ROOT (http://root.cern.ch/drupal/)

Does the new version supersede the previous version?: Yes

Nature of problem:
Particle production via statistical hadronization in relativistic heavy-ion collisions

Solution method:
Monte-Carlo simulation, analyzed with ROOT

Reasons for new version:
The increased functionality of the code contains the following important features.
The input of any shape of the freeze-out hypersurface and the expansion velocity field, including the 3+1 dimensional profiles, in particular those generated externally with the various popular hydrodynamic codes. The hypersurfaces may have variable thermal parameters, which allows for studies departing significantly from the mid-rapidity region. We include a library of standard sets of hypersurfaces and velocity profiles describing the RHIC Au+Au and the LHC Pb+Pb data. A separate code, FEMTO-THERMINATOR, is provided to carry out the analysis of femtoscopic correlations.

Summary of revisions:
THERMINATOR 2 incorporates major revisions to encompass the enhanced functionality.
  • Classes: The Integrator class has been expanded and a new subgroup of classes defined
    • Model and abstract class: These classes are responsible for the physical models of the freeze-out process. The functionality and readability of the code has been substantially increased by implementing each freeze-out model in a different class.
    • The Hypersurface class was added to handle the input form hydrodynamic codes. The hydro input is passed to the program as a lattice of the freeze-out hypersurface. That information is stored in the .xml files.
  • Input: THERMINATOR 2 programs are now controlled by *.ini type files. The programs parameters and the freeze-out model parameters are now in separate ini files.
  • Output:
    • The event files generated by the therm2_events program are not backward compatible with the previous version. The event*.root file structure was expanded with two new TTree structures. From the particle entry it is possible to back-trace the whole cascade.
    • Event text output is now optional.
    • The ROOT macros produce the *.eps figures with physics results, e.g. the pT-spectra, the elliptic-flow coefficient, rapidity distributions, etc.
    • The THERMINATOR HBT package creates the ROOT files femto*.root (therm2_femto) amd hbtfit*.root (therm2_hbtfit).
  • Directory structure: The directory structure has been reorganised. Source code resides in the build directory. The freeze-out model input files, event files, ROOT macros are stored separately.
    The THERMINATOR 2 system, after installation, is able to run on a cluster.
  • Scripts: The package contains a few BASH scripts helpful when running e.g. on a cluster the whole system can be executed via a single script.

Additional comments:
Typical data file size: default configuration
  • 45 MB/500 events
  • 35 MB/correlation file (one kT bin)
  • 45 kB/fit file (projections and fits)

Running time:
Default configuration at 3GHz.
  • primordial multiplicities 70 min (calculated only once per case)
  • 8 min/500 events
  • 10 min - draw all figures
  • 25 min/one kT bin in the HBT analysis with 5000 events