A Brief Introduction to PYTHIA 8.1
Pith reviewed 2026-05-13 11:38 UTC · model grok-4.3
The pith
PYTHIA 8.1 is a complete C++ rewrite of the event generator with new physics models for LHC studies.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
PYTHIA 8 represents a complete rewrite in C++. The current release is the first main one after this transition, and does not yet in every respect replace the old code. It does contain some new physics aspects, on the other hand, that should make it an attractive option especially for LHC physics studies.
What carries the argument
The PYTHIA 8 event generator, a library of hard processes and models for initial- and final-state parton showers, multiple parton-parton interactions, beam remnants, string fragmentation and particle decays.
Load-bearing premise
The new C++ code and added physics features will prove stable enough and match data well enough for the community to adopt them despite the incomplete replacement of the Fortran version.
What would settle it
If detailed comparisons with LHC data reveal systematic discrepancies that were absent in the Fortran version, or if users continue to rely exclusively on the older code for production studies, the claim that the new release is attractive would be falsified.
read the original abstract
The PYTHIA program is a standard tool for the generation of high-energy collisions, comprising a coherent set of physics models for the evolution from a few-body hard process to a complex multihadronic final state. It contains a library of hard processes and models for initial- and final-state parton showers, multiple parton-parton interactions, beam remnants, string fragmentation and particle decays. It also has a set of utilities and interfaces to external programs. While previous versions were written in Fortran, PYTHIA 8 represents a complete rewrite in C++. The current release is the first main one after this transition, and does not yet in every respect replace the old code. It does contain some new physics aspects, on the other hand, that should make it an attractive option especially for LHC physics studies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript is a brief introduction to PYTHIA 8.1, a complete rewrite in C++ of the PYTHIA event generator. It describes the core components including libraries of hard processes, initial- and final-state parton showers, multiple parton interactions, beam remnants, string fragmentation, particle decays, and utilities with external interfaces. The text states that this first major post-transition release does not yet fully replace the prior Fortran version in every respect but incorporates new physics aspects that make it attractive for LHC studies.
Significance. If the descriptive claims hold, the work is significant because PYTHIA is a standard tool for high-energy collision simulations. The C++ implementation and added physics modules (updated showers and MPI) provide a modern platform that facilitates integration with contemporary analysis frameworks and supports LHC physics studies, with the manuscript correctly noting the incomplete replacement of legacy functionality.
minor comments (1)
- Abstract: the sentence beginning 'It does contain some new physics aspects, on the other hand' would read more clearly if 'on the other hand' were replaced by 'however' to improve flow without altering meaning.
Simulated Author's Rebuttal
We thank the referee for the positive assessment and recommendation to accept the manuscript. The referee's summary accurately reflects the scope of our brief introduction to PYTHIA 8.1, including its C++ implementation, core physics models, and current status relative to the Fortran version.
Circularity Check
No significant circularity; descriptive software release note
full rationale
This paper is a brief introduction to the PYTHIA 8.1 event generator and its C++ rewrite. It describes existing physics models (hard processes, showers, MPI, fragmentation) and new features without presenting any derivations, predictions, or first-principles results. No equations or claims reduce to fitted parameters or self-citations within the text; all statements concern implemented code verifiable by inspecting the released source and external data comparisons. The manuscript is self-contained as documentation.
Axiom & Free-Parameter Ledger
free parameters (1)
- parton-shower and fragmentation parameters
axioms (1)
- domain assumption Validity of perturbative QCD for hard processes and string fragmentation for hadronization
Forward citations
Cited by 60 Pith papers
-
Evidence for the decay $B^0_s\to\phi\eta'$
First evidence for B_s^0 to phi eta-prime decay with relative branching ratio (3.56 ± 0.79 ± 0.18 ± 0.06) x 10^{-2} and absolute branching fraction (0.66 ± 0.15 ± 0.03 ± 0.02) x 10^{-6}.
-
Observation of the charmless purely baryonic decay $\mathinner{\mathit{\Lambda}^0_b\!\to \mathit{\Lambda} p \overline{p}}$
First observation of Λ_b^0 → Λ p p-bar with 5.1σ significance and relative branching fraction (5.1 ± 1.3(stat) ± 0.3(syst)) × 10^{-2} to the reference mode Λ_b^0 → Λ K^+ K^-.
-
First evidence of the decay $B^+\to\pi^+ e^+ e^-$
First evidence of B⁺→π⁺e⁺e⁻ decay at 3.2σ significance with branching fraction (2.4 +0.9-0.8 +0.4-0.2)×10^{-8} consistent with SM prediction.
-
First evidence of the decay $B^+\to\pi^+ e^+ e^-$
LHCb observes 3.2 sigma evidence for B+ -> pi+ e+ e- with measured branching fraction (2.4 +0.9 -0.8 +0.4 -0.2) x 10^-8, consistent with Standard Model.
-
Angular analysis of the $B^+\to\pi^+\mu^+\mu^-$ decay
First measurement of A_FB and F_H in B+→π+μ+μ− decay is consistent with Standard Model predictions in both high- and low-mass dimuon regions.
-
Observation of a new excited charm-strange meson $D_{s1}(2933)^+$ in $B^0\to D^+ D^- K^+ \pi^-$ decays
A new charm-strange resonance D_s1(2933)^+ with J^P=1^+ is observed at >10 sigma in B^0 to D+ D- K+ pi- decays, with measured mass 2933 MeV and width 72 MeV.
-
Observation of a cross-section enhancement near the $t\bar{t}$ production threshold in $\sqrt{s}=13$ TeV $pp$ collisions with the ATLAS detector
ATLAS reports an 8+ sigma excess in ttbar production near threshold, consistent with NRQCD quasi-bound states and measuring 9.3 pb.
-
First evidence of the $B_s^0\rightarrow K^-\pi^+\gamma$ decay
LHCb reports the first experimental evidence for the Bs0 to K- pi+ gamma decay with 3.5 sigma significance and measures its branching fraction ratio to the B0 decay in two mass ranges.
-
Searches for $B^0\to K^+\pi^-\tau^+\tau^-$ and $B_s^0\to K^+K^-\tau^+\tau^-$ decays
LHCb reports the first upper limits on B0 → K+π−τ+τ− and Bs0 → K+K−τ+τ− branching fractions, with recast limits of 2.8×10−4 on B0 → K*(892)0 τ+τ− at 95% CL that improve prior bounds by an order of magnitude.
-
First observation of $CP$ violation and measurement of polarization in $B^+\to\rho(770)^0 K^*(892)^+$ decays
First observation of CP violation in B+ to rho(770)^0 K*(892)+ decays with A_CP = 0.507 and longitudinal polarization fraction f_L = 0.720 from 9 fb^-1 of LHCb data.
-
Nonflow Subtraction Beyond Two-Particle Correlations
A nonflow subtraction framework for m-particle cumulants is developed and tested in HIJING simulations for O+O and d+Au collisions.
-
Reweighting Adversarial Networks for Unbinned Unfolding
RANs generalize moment unfolding to full phase-space unbinned unfolding via detector-level Wasserstein critics without requiring support overlap or multiple iterations.
-
Search for electroweak scale dijet resonances in pile-up collisions at $\sqrt{s}=13$ TeV with the ATLAS detector
No significant excess found in low-mass dijet spectrum from pile-up collisions; exclusion limits set on Gaussian and simplified dark matter models with 1.3 pb^{-1} effective luminosity.
-
Study of the $B^0 \to \Lambda_c^+ \bar{\Lambda}_c^- K_S^0$ decay
Relative branching fraction B(B0 → Λc+ Λc- KS0)/B(B+ → Λc+ Λc- K+) measured as 0.53 ± 0.05 ± 0.05 with 3.9σ evidence for Ξc(2923)+ and Ξc(2939)+ resonances consistent with isospin partners.
-
Precise measurement of the CKM angle $\gamma$ with a novel approach
A joint fit to LHCb B to D h decays and BESIII quantum-correlated D Dbar data yields gamma = 71.3 plus or minus 5.0 degrees, the most precise measurement to date.
-
Measurement of the CKM angle $\gamma$ in $B^{\pm} \rightarrow D(\rightarrow K^{0}_{\rm S} h^{\prime+}h^{\prime-})h^{\pm}$ decays with a novel approach
A novel model-independent approach with per-event phase-space weights on combined BESIII and LHCb data measures the CKM angle γ as (71.3 ± 5.0)° in B± → D(→ K0S h'+h'-) h± decays.
-
Many Wrongs Make a Right: Leveraging Biased Simulations Towards Unbiased Parameter Inference
Template-Adapted Mixture Model uses many biased simulations for data-driven estimates of signal and background distributions, yielding unbiased signal fraction estimates with well-calibrated uncertainties.
-
Observation of the decay $\chi_{c1}(3872)\rightarrow J\mskip -3mu/\mskip -2mu\psi \mu^+\mu^-$
First observation of χ_c1(3872) → J/ψ μ⁺μ⁻ reported at 6.5σ significance with branching fraction ratio (1.68 ± 0.37) × 10^{-3} relative to the π⁺π⁻ mode.
-
Observation of a cross-section enhancement near the $t\bar{t}$ production threshold in $\sqrt{s}=13$ TeV $pp$ collisions with the ATLAS detector
ATLAS observes an excess over 8 sigma in ttbar production near threshold, consistent with NRQCD color-singlet S-wave quasi-bound states, with measured cross-section 9.3^{+1.4}_{-1.3} pb.
-
(3+1)D event-by-event pre-equilibrium dynamics in heavy-ion collisions
Extension of KoMPoST to (3+1)D non-boost-invariant pre-equilibrium dynamics via kinetic theory response functions, chained into a full McDIPPER+KoMPoST-3D+CLVisc+SMASH simulation chain to study sensitivity of longitud...
-
Measurement of the $W \to \mu \nu_\mu$ cross-sections as a function of the muon transverse momentum in $pp$ collisions at 5.02 TeV
Differential W to muon neutrino cross-section measurement at 5.02 TeV with proof-of-principle W mass extraction from corrected distributions.
-
The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations
MadGraph5_aMC@NLO automates tree-level, NLO, shower-matched, and merged cross-section computations for collider processes in a unified flexible framework.
-
Machine learning fully hadronic events with spectral functions
Spectral functions from two-point correlations serve as multiplicity-independent ML inputs and improve expected gluino mass reach by 150-250 GeV in a fully hadronic ttbar vs gluino benchmark.
-
Reweighting Underlying Event and Colour Reconnection parameter variations in Sherpa
A reweighting method is introduced to efficiently study and tune MPI and CR parameter variations in Sherpa using relative weights, demonstrated with LHC and Tevatron data.
-
"Hadron-in-fat-jet'' AI Tagging to Detect Rare Decays Such as $W^{\pm}\to\pi^{\pm}\gamma$
Proof-of-principle for hadron-in-fat-jet AI tagging that yields an expected 95% CL limit of B(W±→π±γ) < 2.78×10^{-5} at 450 fb^{-1}.
-
Search for soft unclustered energy patterns containing muons in the final state in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector
No significant excess observed in search for SUEP signatures with muons; exclusion limits set on mediator cross section times branching fraction down to 0.05 fb at 750 GeV.
-
NLO QCD and parton-shower effects for Higgs-boson production in association with a hard photon via vector-boson fusion
NLO QCD plus parton-shower matched implementation for VBF Higgs plus photon production, with studies showing small shower effects on Higgs observables and larger effects on sub-leading jets.
-
Uncovering Hidden Systematics in Neural Network Models for High Energy Physics
Neural networks for HEP tasks can be fooled at significant rates by subtle perturbations inside uncertainty envelopes, revealing hidden systematics not captured by conventional methods.
-
SemiCharmTag: a tool for Semileptonic Charm tagging
SemiCharmTag achieves a factor of ~4 signal-over-background improvement at 81% efficiency for Drell-Yan muons using secondary-vertex hadron tagging in LHCb proton-proton collision simulations.
-
Test of lepton flavour universality with $B^0\to K^{*0}\ell^+\ell^-$ decays at large dilepton invariant mass
R_K*0 is measured as 1.08^{+0.14}_{-0.12}(stat) ± 0.07(syst) for q² > 14 GeV²/c⁴ in B⁰ → K*⁰ ℓ⁺ℓ⁻ decays, consistent with the Standard Model.
-
Search for the lepton-flavour violating decays $B^+ \to \pi^+ \mu^\pm e^\mp$
No signal observed for B+ → π+ μ± e∓; branching fraction upper limit set at 1.8 × 10^{-9} at 90% CL.
-
Monte Carlo Event Generation with Continuous Normalizing Flows
Continuous normalizing flows improve unweighting efficiency in Monte Carlo event generation for high-jet-multiplicity collider processes by factors up to 184, with wall-time gains of about ten when combined with coupl...
-
Non-Thermal Production of Sexaquark Dark Matter
Non-thermal production via late-decaying reheatons can achieve the observed dark matter density for sexaquarks by controlling branching fractions and coalescence probabilities, unlike thermal freeze-out which underpro...
-
Measurement of the top-quark mass using decays with a $J/\psi$ meson at $\sqrt{s}=$13 TeV with the ATLAS detector
Top quark mass measured at 172.17 ± 1.56 GeV via unbinned maximum-likelihood fit to m(ℓ μ⁺μ⁻) in ATLAS Run 2 data using J/ψ decays.
-
Measurement of charged-hadron distributions in heavy-flavor jets in proton-proton collisions at $\sqrt{s}$=13 TeV
Charged-hadron distributions in heavy-flavor jets differ from light-quark jets in ways consistent with dead-cone suppression and hard fragmentation of the heavy hadron.
-
Search for $K_{\mathrm{S(L)}}^{0} \rightarrow \pi^{+}\pi^{-}\mu^{+}\mu^{-}$ decays at LHCb
No evidence for KS0 or KL0 to pi+ pi- mu+ mu- decays; first upper limits set at 1.4e-9 and 6.6e-7 (90% CL).
-
Observation of $W^{+}W^{-}\gamma$ production in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector and constraints on anomalous quartic gauge-boson couplings
ATLAS reports the observation of W⁺W⁻γ triboson production with 5.9σ significance and derives constraints on anomalous quartic gauge-boson couplings via effective field theory.
-
Search for Higgs bosons produced in association with a high-energy photon via vector-boson fusion and decaying to a pair of $b$-quarks in the ATLAS detector
ATLAS reports a new measurement of the Higgs signal strength in the VBF plus photon channel with H to bb decay, obtaining 0.2 ± 0.7 using an improved neural-network-based analysis on 133 fb^{-1} of data.
-
Electroweak diboson production in association with a high-mass dijet system in semileptonic final states from $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector
Electroweak diboson plus high-mass dijet production observed at 7.4 sigma with signal strength 1.28, plus first semileptonic-channel limits on S02, T0 and M0 Wilson coefficients.
-
Measurement of $J/\psi$-jet correlations in $pp$ and Pb+Pb collisions at $\sqrt{s_{\text{NN}}}=5.02$ TeV with the ATLAS detector
New measurements of isolated, non-isolated, and inclusive J/ψ yields, nuclear modification factors, and fractions up to 60 GeV in pp and Pb+Pb collisions at 5.02 TeV.
-
Toponium effects on quantum steering and Bell nonlocality of top quarks
Toponium strengthens the spin-singlet component in top-quark pairs, substantially enhancing entanglement and enabling observable quantum steering at 10 sigma and Bell nonlocality at 9 sigma near threshold with current...
-
Measurement of isolated photon plus two-jet correlations in Pb+Pb and $pp$ collisions at 5.02 TeV with ATLAS
ATLAS reports significant suppression (I_AA < 1) of photon-tagged two-jet yields in Pb+Pb versus pp collisions at 5.02 TeV across three observables, compared to JEWEL, JETSCAPE, and LBT models.
-
Bayesian inference constraints on jet quenching across centrality, beam energy, and observable classes in LHC heavy-ion collisions
Bayesian posteriors from JETSCAPE jet-quenching model are largely compatible across centrality but exhibit shifts across beam energy and observable class, with varying ability to predict complementary datasets.
-
Measurement of inclusive production of charmonium states in $b$-hadron decays via their decay into $\phi \phi$
LHCb reports branching fractions B(b→χ_c0,1,2 X) and B(b→η_c(2S)X)×B(η_c(2S)→φφ) plus the most precise η_c(1S) mass from φφ decays in 5.9 fb⁻¹ of data.
-
Hadronic Origin of Sub-PeV Gamma-Ray Emission from LHAASO J0621+3755
Hadronic proton-proton interactions in a superdiffusive one-zone model with proton luminosity 0.14 of spin-down power and diffusion index 1.05 reproduce the 7-200 TeV gamma-ray spectrum from PSR J0622+3749 halo at amb...
-
Precision measurement of CP violation and branching fractions in $B^{\pm} \to K^0_{\mathrm{S}} h^{\pm}$ $(h = \pi, K)$ decays and search for the rare decay $B_c^{\pm} \to K^0_{\mathrm{S}} K^{\pm}$
Most precise CP asymmetry and branching ratio measurements for B± → K0S h± decays with no significant signal observed for the rare Bc decay.
-
Observation and investigation of the $T_{c\bar{c}1}(4430)^{+}$ structure in $B^{+} \to \psi(2S) K_{\text{S}}^{0} \pi^{+}$ decays
First 4D amplitude analysis of B+ → ψ(2S) K_S0 π+ confirms the exotic T_{c c-bar 1}(4430)+ structure with properties consistent with prior observations in the isospin-related channel.
-
Measurement of the $W$-boson angular coefficients and transverse momentum in $pp$ collisions at $\sqrt{s}=$ 13 TeV with the ATLAS detector
ATLAS reports the first measurement of the complete angular coefficients and pT-differential cross sections for W+ and W- bosons in full lepton phase space at 13 TeV, finding agreement with QCD predictions up to order α_S².
-
DELPHES 3, A modular framework for fast simulation of a generic collider experiment
DELPHES 3 delivers a modular fast-simulation framework with particle-flow and pile-up features for reconstructing physics objects in collider detector studies.
-
Refined anti-proton and anti-deuteron fluxes from weak-scale Dark Matter
Updated tabulated interstellar fluxes of antiprotons and antideuterons from weak-scale DM are provided using new propagation schemes and inelastic cross-sections, showing greater robustness than prior calculations.
-
Probing anomalous quartic gauge couplings in same-sign $W$ boson scattering with polarization and spin correlation
Spin-correlation asymmetries in same-sign WW production yield sensitivity to anomalous WWWW couplings comparable to transverse-mass distributions, and their combination improves Wilson-coefficient limits while respect...
-
Energy-energy correlators inside single inclusive jets in heavy-ion collisions with CoLBT-hydro model
Updated CoLBT-hydro simulations with Q_M=2.0 GeV reproduce CMS in-jet EEC data, validate background subtraction, and show path-length and diffusion-wake effects.
-
Forward hadron production in pp collisions at LHC energies from an event generator based on the color glass condensate framework
MC-CGC simulations with rcBK evolution favor HERA-tuned MV^γ and MV^e initial conditions for LHCb forward hadron data, show kT factorization outperforming DHJ at mid-rapidity, and provide predictions for ALICE FoCal n...
-
Novel Machine Learning Methods to Improve Z Pole Integrated Luminosity at Future Colliders
Gradient boosted decision trees suppress diphoton backgrounds while adaptive symbolic memetic regression corrects beam deflection biases, reaching luminosity uncertainties below 10^{-4} and 5x10^{-6}.
-
$C\!P$ violation analysis of local and nonlocal amplitudes in the $\overline{B}^0 \to \overline{K}^{*0}\mu^+\mu^-$ decay
LHCb measures CP-violating Wilson coefficients in B0 -> K*0 mu+ mu- with higher precision than before and finds results consistent with the Standard Model.
-
Measurement of $\gamma$ using $B^{\pm}\rightarrow DK^{\pm}$ and $B^{\pm}\rightarrow D\pi^{\pm}$ decays with $D\rightarrow K_{\rm S}^{0}\pi^{+}\pi^{-}$ and $D\rightarrow K_{\rm S}^{0}K^{+}K^{-}$
The CKM angle γ is measured to be (68.1 ± 6.7)° from CP violation observed in the Dalitz plots of B± → DK± and B± → Dπ± decays with D → KS0π+π− and KS0K+K−.
-
Measurement of the $W$-boson production cross-sections in $pp$ collisions at $\sqrt{s}$ = 13 TeV in the forward region
LHCb measures forward W+ and W- production cross-sections of 1754.2 pb and 1178.1 pb at 13 TeV, agreeing with NNLO QCD predictions at higher precision than prior results.
-
Precision measurement of the muon charge asymmetry from $W$-boson decays in $pp$ collisions at $\sqrt{s}$ = 13 TeV in the forward region
Muon charge asymmetry from W decays is measured with highest precision in the forward region at 13 TeV and agrees with NNLO pQCD.
-
Muonphilic asymmetric dark matter at a future muon collider
Muonphilic portals to fermionic asymmetric dark matter are constrained by existing data and can be probed further by 3 and 10 TeV muon colliders.
-
Branching fraction measurement of the $\mathit{\Lambda} \to p \mu^- \overline{\nu}_{\mu}$ decay
Branching fraction B(Λ → p μ⁻ ν̄_μ) measured as (1.462 ± 0.016 ± 0.100 ± 0.011) × 10^{-4}, improving prior precision by a factor of two and yielding R^{μe} = 0.175 ± 0.012 consistent with the Standard Model.
Reference graph
Works this paper leans on
-
[1]
Sj¨ ostrand, Computer Physics Commun
T. Sj¨ ostrand, Computer Physics Commun. 27 (1982) 243, 28 (1983) 229, 39 (1986) 347; T. Sj¨ ostrand and M. Bengtsson, Computer Physics Commun. 43 (1987) 367
work page 1982
-
[2]
Bengtsson, Computer Physics Commun
H.-U. Bengtsson, Computer Physics Commun. 31 (1984) 323; H.-U. Bengtsson and G. Ingelman, Computer Physics Commun. 34 (1985) 251; H.-U. Bengtsson and T. Sj¨ ostrand, Computer Physics Commun . 46 (1987) 43; T. Sj¨ ostrand, Computer Physics Commun. 82 (1994) 74
work page 1984
-
[3]
T. Sj¨ ostrand, P. Ed´ en, C. Friberg, L. L¨ onnblad, G. Miu , S. Mrenna and E. Norrbin, Computer Physics Commun. 135 (2001) 238
work page 2001
-
[4]
T. Sj¨ ostrand, S. Mrenna and P. Skands, JHEP 05 (2006) 026 [hep-ph/0603175]
work page Pith review arXiv 2006
-
[5]
L¨ onnblad, Computer Physics Commun
L. L¨ onnblad, Computer Physics Commun. 118 (1999) 213; M. Bertini, L. L¨ onnblad and T. Sj¨ ostrand, Computer Physic s Commun. 134 (2001) 365
work page 1999
-
[6]
see webpage http://www.thep.lu.se/ThePEG/
-
[7]
S. Gieseke, A. Ribon, M.H. Seymour, P. Stephens and B.R. W ebber, JHEP 02 (2004) 005; see webpage http://hepforge.cedar.ac.uk/herwig/
work page 2004
-
[8]
Generic User Process Interface for Event Generators
E. Boos et al., in the Proceedings of the Workshop on Physi cs at TeV Colliders, Les Houches, France, 21 May - 1 Jun 2001 [hep-ph/0109068]
work page internal anchor Pith review arXiv 2001
-
[9]
Alwall et al., Computer Physics Comm
J. Alwall et al., Computer Physics Comm. 176 (2007) 300 26
work page 2007
-
[10]
M.R. Whalley, D. Bourilkov and R.C. Group, in ‘HERA and t he LHC’, eds. A. De Roeck and H. Jung, CERN-2005-014, p. 575 [hep-ph/0 508110]
work page 2005
- [11]
- [12]
- [13]
- [14]
-
[15]
B. Andersson, G. Gustafson, G. Ingelman and T. Sj¨ ostra nd, Phys. Rep. 97 (1983) 31
work page 1983
- [16]
- [17]
- [18]
- [19]
- [20]
- [21]
-
[22]
see webpage http://proj-clhep.web.cern.ch/proj-clhep/
- [23]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.