REVIEW
THERMUS -- A Thermal Model Package for ROOT
High Energy Physics - Phenomenology · Submitted 15 Jul 2004 · 61 references
Pith's verdict: A solid, honestly documented tool paper whose headline feature — per-yield decay chains and efficiencies — is asserted but never actually exercised in the worked example.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
A ROOT toolkit computes thermal-model particle yields in three statistical ensembles.
The reading
What carries the argument
The load-bearing mechanism is the ensemble-specific correction factor: in the canonical and strangeness-canonical treatments, the Boltzmann-approximation density, energy density and pressure of a species are the grand-canonical values with the canonically conserved chemical potential set to zero, multiplied by a factor such as $Z_{B-B_i,S-S_i,Q-Q_i}/Z_{B,S,Q}$ or $Z_{S-S_i}/Z_S$. THERMUS computes those factors, stores them in ensemble-specific thermal-model classes, and uses them to populate a hash table of densities keyed by PDG ID. Around that core it wraps particle and decay bookkeeping, a Broyden solver for constraints, and TMinuit for fits, with each yield able to point to its own particle set and decay chain.
Extended reading notes
Core claim
The paper's central claim is that a single ROOT-integrated package can carry out the full chain of a statistical-thermal analysis: define a hadron-resonance gas, choose an ensemble, compute primordial particle, energy, entropy densities and pressure, add resonance-decay contributions to stable particles, impose constraints such as $B/2Q$, strangeness or $E/N$, and fit thermal parameters to experimental yields with the ROOT implementation of MINUIT. The package is organised as a hierarchy of classes for particles, particle sets, parameter sets, thermal particles, thermal models and thermal fits, with ensemble-specific subclasses. Its practical contribution is not a new thermal-model formula but the assembly of the cited formulas into a reusable, interactive tool with per-yield decay chains and efficiency settings.
Load-bearing premise
The implementation reproduces the cited canonical and strangeness-canonical formalisms correctly, and the shipped default particle list, mass cutoff and rescaled branching ratios are adequate; the paper does not benchmark the package against measured yields, using fictitious ratios in its worked fit.
Editorial extensions
If this is right
- A user can go from a particle list and trial thermal parameters to primordial densities, decay feed-down and fitted freeze-out parameters without leaving ROOT.
- The same code base can treat high-energy collisions grand-canonically and low-energy collisions with strangeness treated canonically, so ensemble choices become a switch rather than a separate code.
- Constraints on $B/2Q$, strangeness, $E/N$, total baryon density or $s/T^3$ can be imposed self-consistently while fitting, which makes the thermal parameters satisfy the collision's conservation conditions.
- Quantum statistics, resonance width and excluded-volume corrections can be toggled, allowing systematic studies of how each effect changes the extracted parameters.
- Per-yield decay chains and efficiencies let the model match the feed-down corrections of a specific experimental data set rather than a single universal chain.
Reading between the lines
- If the package reproduces published fits, it could become a common baseline for comparing thermal-model results across experiments, provided the particle-list version and decay files are recorded with each fit.
- The per-yield decay-chain machinery could be extended to propagate branching-ratio and efficiency uncertainties into the fitted parameters; the paper does not quantify that error budget.
- A direct validation would be to recompute the canonical correction factors independently, for example by numerical integration of the partition-function integrals, and compare them with THERMUS across a grid of $T$, volume and chemical potentials.
- Because the data structures already allow charm and beauty quantum numbers, extending the default particle list to charmed hadrons would let the same workflow test $\gamma_C$ at collider energies, although that extension is not demonstrated here.
Original abstract
THERMUS is a package of C++ classes and functions allowing statistical-thermal model analyses of particle production in relativistic heavy-ion collisions to be performed within the ROOT framework of analysis. Calculations are possible within three statistical ensembles; a grand-canonical treatment of the conserved charges B, S and Q, a fully canonical treatment of the conserved charges, and a mixed-canonical ensemble combining a canonical treatment of strangeness with a grand-canonical treatment of baryon number and electric charge. THERMUS allows for the assignment of decay chains and detector efficiencies specific to each particle yield, which enables sensible fitting of model parameters to experimental data.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
THERMUS is a package of C++ classes and functions allowing statistical-thermal model analyses of particle production in relativistic heavy-ion collisions to be performed within the ROOT framework of analysis. Calculations are possible within three statistical ensembles; a grand-canonical treatment of the conserved charges B, S and Q, a fully canonical treatment of the conserved charges, and a mixed-canonical ensemble combining a canonical treatment of strangeness with a grand-canonical treatment of baryon number and electric charge. THERMUS allows for the assignment of decay chains and detector efficiencies specific to each particle yield, which enables sensible fitting of model parameters to experimental data.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | THERMUS is a package of C++ classes and functions allowing statistical-thermal model analyses of particle production in relativistic heavy-ion collisions to be performed within the ROOT framework of analysis. | contribution 0.97 | 0.53 | 0.19 | asserted 1.00 |
| s1 | Calculations are possible within three statistical ensembles; a grand-canonical treatment of the conserved charges B, S and Q, a fully canonical treatment of the conserved charges, and a mixed-canonical ensemble combining a canonical treatment of strangeness with a grand-canonical treatment of baryon number and electric charge. | method 0.85 | 0.45 | 0.07 | hedged possibility 0.59 |
| s2 | THERMUS allows for the assignment of decay chains and detector efficiencies specific to each particle yield, which enables sensible fitting of model parameters to experimental data. | method 0.64 | 0.51 | 0.15 | asserted 0.98 |
What was examined, and what was not
- arXiv listing abstract v2
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:hep-ph/0407174 by the authors.
References
How later papers used it
62 later papers cite it
9 later papers by other authors are listed below: 9 used its method or tool. 3 more papers by the same authors are listed: 3 used its method or tool.
Uses its method or tool 12 of 25 papers with this label · 3 by the same authors
-
Obtained its proton and antiproton results using the THERMUS code.
More uses for Thermal Models (2026) · 0.99 · paraphrased: the licence does not allow quoting
-
Did all its fits using the THERMUS thermal code with particle ratios calculated from mid-rapidity yields.
Energy and System Size Dependence of Strangeness Production from SPS to RHIC (2007) · 1.00 · paraphrased: the licence does not allow quoting
-
Used the THERMUS package to include the contribution from resonance decays, with the chemical freeze-out temperature and baryon chemical potential determined in another work.
Probing QCD critical fluctuations from the yield ratio of strange hadrons in relativistic heavy-ion collisions (2019) · 0.97 · paraphrased: the licence does not allow quoting
-
Used the thermal statistical model THERMUS to estimate the relative yield of primary particles.
Global spin polarization of multistrange hyperons and feed-down effect in heavy-ion collisions (2021) · 0.99 · paraphrased: the licence does not allow quoting
-
Used the Baryon Canonical ensemble approach as implemented in the THERMUS package to perform the fitting to the data.
Canonical treatment of strangeness and light nuclei production (2025) · 1.00 · paraphrased: the licence does not allow quoting
-
Utilized the open code THERMUS ver.2.1 to calculate the entropy density and baryon number density at various temperatures.
Phase diagram of hot and dense QCD constrained by the Statistical Model (2010) · 1.00 · paraphrased: the licence does not allow quoting
-
Utilised the implementation in the THERMUS code for its study.
Multiplicity dependence of light flavour hadron production at LHC energies in the strangeness canonical suppression picture (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
We have performed our analysis for grand canonical, canonical and strangeness canonical ensembles using a thermal model called THERMUS [ref] .
Multiplicity dependence freeze-out scenarios in pp collisions at $\sqrt{s}$ = 7 TeV (2021) · Multiplicity dependence freeze-out scenarios in p+p collision at LHC · 0.86
-
Determined the yield of particle A at a point on the trajectory with the statistical thermal model THERMUS.
Search for the QCD Critical Point by Transverse Velocity Dependence of Anti-deuteron to Deuteron Ratio (2018) · 1.00 · paraphrased: the licence does not allow quoting
By the same authors
-
Did all its calculations using the latest version of THERMUS.
A Comparison of p-p, p-Pb, Pb-Pb Collisions in the Thermal Model: Multiplicity Dependence of Thermal Parameters (2018) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Calculated the excluded volume correction factor and the entropy density using the THERMUS package.
Viscosity in the excluded volume hadron gas model (2007) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Used an extended version of the THERMUS framework for its calculations.
Multiplicity fluctuations in relativistic nuclear collisions (2007) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
How the literature uses it
50 arXiv papers cite this in 72 sentences
Papers whose main result relates to this one's
States the same result 1 paper
-
THERMUS is a package of C++ classes and functions allowing statistical-thermal model analyses of particle production in relativistic heavy-ion collisions to be performed within the ROOT framework of analysis.
THERMUS · 1.00
Papers citing it · 32 sentences without a confident label
Uses its method or tool 21 papers
-
All fits were done using the THERMUS thermal code [ref] , using particle ratios calculated from mid-rapidity yields (dN/dy).
Energy and System Size Dependence of Strangeness Production from SPS to RHIC · Energy Dependence · 1.00
-
However, corrections due to the weak-decay feed-down for protons still need to be performed. for STAR's $p$ and $p̅$ at were obtained using the THERMUS code [ref] .
More uses for Thermal Models · Data and analysis details · 0.99
-
The ratio $𝒪_K-$Ξ$--$ in this model is then given by [equation] To include the contribution from resonance decays, we use the THERMUS package [ref] with the chemical freeze-out temperature and baryon chemical potential determined in Ref. [ref] .
Probing QCD critical fluctuations from the yield ratio of strange hadrons in relativistic heavy-ion collisions · Introduction · 0.97
Uses its data or values 4 papers
-
The masses and branching ratios used are as given in [ref] .
Novel scheme for parametrizing the chemical freeze-out surface in Heavy Ion Collision Experiments · Introduction · 0.99
-
HIJING [ref] simulations indicate that $F$ = 0.67 and if the scales with the yield, THERMUS [ref] indicates that $F$ = 0.532.
Measurement of transverse energy at midrapidity in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV · Correction factors and · 0.87
-
The masses and branching ratios used are as given in [ref] .
Systematics of chemical freeze-out parameters in heavy-ion collision experiments · Data Analysis · 0.96
Compares with it 3 papers
-
In this paper, we implement this technique on published hadron yields and compare the extracted chemical potential ratios with the results of thermal-model fit obtained by STAR [ref] using THERMUS [ref] .
More uses for Thermal Models · Introduction · 0.98
-
Finally, the values and behaviour of $$ as a function of centrality are in very good agreement with previous findings, especially with the ones [ref] calculated with the THERMUS [ref] package.
Chemical freeze-out in ultra-relativistic heavy ion collisions at sqrt(s)_NN = 130 and 200 GeV · Comparison with previous analyses · 0.96
-
Figure [ref] shows the /(no contribution considered here), / and / ratios compared to calculations from grand-canonical thermal models [ref] , which were found to satisfactorily describe central Pb-Pb data at = 2.76 TeV.
Multiplicity dependence of (multi-)strange hadron production in proton-proton collisions at $\sqrt{s}$ = 13 TeV · Results and discussion · 0.80
Cites it as background 5 papers
-
Strange hadron production in heavy-ion collisions is currently usually described in the framework of statistical-hadronisation (or thermal) models [ref] .
Multiplicity dependence of (multi-)strange hadron production in proton-proton collisions at $\sqrt{s}$ = 13 TeV · Introduction · 0.90
-
Examples of commonly used hadron lists can be found in [ref] where, in some instances, additional states from the quark model have been included [ref] .
Updated Hadron List for Transport Simulations of Heavy-Ion Collisions · Introduction · 0.85
-
In canonical ensemble (CE) all the charges are exactly conserved (CE) [ref] whereas in strangeness canonical ensemble (SCE) [ref] only strangeness is exactly conserved.
Multiplicity dependence freeze-out scenarios in pp collisions at $\sqrt{s}$ = 7 TeV · Canonical and strangeness canonical ensemble · 0.96
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.