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.
-
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.
-
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.
-
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.
-
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.
-
Octet scalars shaping LHC distributions in 4-jet final states
A color-octet scalar Θ pair-produced via gluons and decaying to quark pairs can account for the CMS excess in equal-mass dijet pairs at 0.95 TeV with a 65 fb cross section for the real case.
-
Measurement of high-mass $t\bar{t}\ell^{+}\ell^{-}$ production and lepton flavour universality-inspired effective field theory interpretations at $\sqrt{s}=13$ TeV with the ATLAS detector
ATLAS reports a measurement of high-mass ttbar ll production with no observed deviation from the Standard Model and derives EFT constraints on four-fermion operators from 140 fb^{-1} of Run 2 data.
-
Search for Higgs boson exotic decays into Lorentz-boosted light bosons in the four-$\tau$ final state at $\sqrt{s}=13$ TeV with the ATLAS detector
Search for H->aa->4tau (4<m_a<15 GeV) with boosted di-tau reconstruction finds no excess and sets 95% CL upper limits of 0.03-0.10 on (sigma_H/sigma_SM)*BR(H->aa->4tau).
-
Luminosity determination in $pp$ collisions at $\sqrt{s}=13$ TeV using the ATLAS detector at the LHC
The integrated luminosity of the ATLAS Run 2 pp dataset at 13 TeV is determined to be 140.1 ± 1.2 fb^{-1}.
-
TMD parton distributions from parton showers
Effective TMD parton distributions are derived from PYTHIA8 and HERWIG6 parton showers and compared with those from the Parton Branching method.
-
Impacts of radiative corrections on measurements of lepton flavour universality in $B \to D \ell \nu_{\ell}$ decays
Radiative corrections introduce up to 7% bias on R(D) measurements in LHCb-like pseudo-experiments when using PHOTOS for soft-photon simulation.
-
An Introduction to PYTHIA 8.2
PYTHIA 8.2 is a mature C++ event generator that combines hard processes, parton showers, multiparton interactions, and string fragmentation into a complete simulation framework for high-energy collisions.
-
USMEFT as a tool for discovery of universal new physics at high luminosity LHC
Universal SMEFT fits to pseudo-data from neutral and charged Drell-Yan processes at HL-LHC can detect universal new physics and extract its properties stably across EFT truncation orders.
-
Search for the decays $B_{(s)}^0\to J/\psi\gamma$ at LHCb
Upper limits of 2.9×10^{-6} for B_s^0 and 2.5×10^{-6} for B^0 on the branching fractions to J/ψγ at 90% CL, with the B_s limit improved by a factor of 2.5.
-
Search for signatures of electroweakinos with photons, jets, and large missing transverse momentum in $\sqrt{s}=13$ TeV pp collisions with the ATLAS detector
No significant excess observed; gauginos excluded up to 1.2 TeV in specific supersymmetric decay scenarios.
-
FastJet user manual
FastJet is a C++ package providing implementations of sequential recombination jet algorithms, cone algorithms via plugins, jet substructure tools, and pileup estimation for pp and e+e- collisions.
-
VBSCan Thessaloniki 2018 Workshop Summary
The document reports the first year of activity of the VBSCan COST Action network on vector-boson scattering phenomenology and experiments from a 2018 workshop.
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.