---
title: Jets, Hadronization, and Testing QCD
module: Quantum Chromodynamics
moduleNumber: 8
lessonNumber: 4
order: 804
summary: >
  Quarks and gluons produced in a collision fragment into collimated sprays of
  hadrons — jets — whose directions track the underlying partons. This lesson reads
  two-jet events as the quark and antiquark of electron-positron annihilation,
  three-jet events as direct evidence of the radiated gluon, and the hadronization
  step as the flux tube breaking into color singlets. Jet algorithms and event-shape
  variables turn the pattern into precision measurements of alpha_s.
topics: [Quantum Chromodynamics]
draft: false
sources:
  - book: Griffiths
    ref: "Ch. 9 — Quantum Chromodynamics, §9.5"
  - book: Thomson
    ref: "Ch. 10 §10.6 (jets and the measurement of alpha_s)"
  - book: Perkins
    ref: "Ch. 6 (quarks and gluons in hard collisions; jets)"
  - book: Tong
    ref: "The Standard Model (Cambridge Part III), §3.1 (asymptotic freedom and hadronic final states)"
---

[Confinement](/particle-physics/qcd/asymptotic-freedom-confinement) forbids a free
quark, yet high-energy collisions manifestly produce quarks and gluons. The
resolution is that a struck parton does not travel far before its
[color flux tube](/particle-physics/qcd/asymptotic-freedom-confinement) breaks into
hadrons, and because the fragmentation is soft it barely deflects the parton's
momentum. The result is a **jet**: a narrow cone of hadrons whose total momentum
reconstructs the parton that produced it. Jets make the quarks and gluons of the
Lagrangian visible as collimated sprays, and their rates and shapes are among the most
precise tests of QCD. This lesson reads the jet structure of electron-positron
annihilation, identifies the gluon in three-jet events, and extracts $\alpha_s$ from
the pattern.

Throughout, $\hbar = c = 1$.

## Two-jet events

The reference reaction is $e^+ e^- \to \gamma^\ast \to q\bar q$, the quark version of
the muon-pair process of the [QED tree
lesson](/particle-physics/qed/qed-tree-processes). The quark and antiquark emerge
back to back in the center-of-momentum frame, each fragmenting into a jet. Two
observations confirm the partonic origin. First, the two jets are collinear and
opposite, as required by momentum conservation for a two-body partonic final state.
Second, the jet axis follows the angular distribution

$$
\frac{\d\sigma}{\d\Omega} \propto 1 + \cos^2\theta ,
$$

the distribution for producing a spin-$\tfrac12$ pair from a virtual photon. A
spin-$0$ pair would give $\sin^2\theta$ instead. The measured $1 + \cos^2\theta$
confirms that the jets originate from spin-$\tfrac12$ quarks.[^grif-jets]

The total annihilation rate into hadrons, summed over the accessible flavors and
colors, is the $R$-ratio measured in the
[quark model](/particle-physics/quark-model/color-confinement-exotics); the two-jet
angular shape confirms the spin of the constituents that $R$ counts.

$$
% caption: A two-jet event in electron-positron annihilation. The virtual photon
% produces a quark and antiquark back to back; each fragments into a collimated spray
% of hadrons along the original parton direction. The jet axis follows the
% one-plus-cosine-squared distribution characteristic of a spin-half pair.
\begin{tikzpicture}[>=stealth, font=\footnotesize, scale=1.0]
  \definecolor{acc}{HTML}{4A6FA5}
  % beams
  \draw[black, thick, ->] (-3.6,0) -- (-2.2,0);
  \node[black, left, font=\scriptsize] at (-3.6,0) {$e^-$};
  \draw[black, thick, ->] (3.6,0) -- (2.2,0);
  \node[black, right, font=\scriptsize] at (3.6,0) {$e^+$};
  \fill (0,0) circle (2.4pt);
  % jet 1 (up-right cone)
  \foreach \a in {28,34,40,46,52}
    \draw[thick] (0,0) -- ({2.6*cos(\a)},{2.6*sin(\a)});
  \node at (2.4,2.1) {jet};
  % jet 2 (down-left cone)
  \foreach \a in {208,214,220,226,232}
    \draw[thick] (0,0) -- ({2.6*cos(\a)},{2.6*sin(\a)});
  \node at (-2.4,-2.1) {jet};
\end{tikzpicture}
$$

## Three-jet events and the gluon

A quark can radiate a hard gluon before hadronizing, $e^+ e^- \to q\bar q g$, and when
the gluon is energetic and well separated it fragments into a jet of its own. The
event then has three jets in a plane. Three-jet events were observed at the PETRA
collider at DESY in 1979 and provided the first direct evidence for the gluon: their
rate is fixed by $\alpha_s$, and their angular correlations match a spin-$1$ radiated
quantum rather than a spin-$0$ one.[^thom-gluon] The fraction of three-jet to two-jet
events measures the coupling directly,

$$
\frac{\sigma_{3\text{-jet}}}{\sigma_{2\text{-jet}}} \sim \frac{\alpha_s}{\pi},
$$

the QCD analogue of hard-photon bremsstrahlung in QED. The gluon jet is
distinguishable statistically: gluons carry the adjoint color charge $C_A = 3$ against
the quark's $C_F = \tfrac43$, so gluon jets radiate more and are broader with higher
hadron multiplicity.

$$
% caption: A three-jet event. One of the quarks radiates a hard gluon before
% hadronizing, and the gluon fragments into its own jet. The three jets lie in a
% plane; the rate relative to two-jet events is proportional to alpha_s, and the
% observation at PETRA in 1979 was the first direct evidence for the gluon.
\begin{tikzpicture}[>=stealth, font=\footnotesize, scale=1.0]
  \definecolor{acc}{HTML}{4A6FA5}
  \fill (0,0) circle (2.4pt);
  % quark jet (right)
  \foreach \a in {-14,-8,-2,4,10}
    \draw[thick] (0,0) -- ({2.7*cos(\a)},{2.7*sin(\a)});
  \node at (2.9,0.2) {quark};
  % antiquark jet (upper-left)
  \foreach \a in {128,134,140,146,152}
    \draw[thick] (0,0) -- ({2.7*cos(\a)},{2.7*sin(\a)});
  \node at (-2.4,1.9) {antiquark};
  % gluon jet (lower-left), the highlight
  \foreach \a in {236,242,248,254,260}
    \draw[acc, thick] (0,0) -- ({2.5*cos(\a)},{2.5*sin(\a)});
  \node[acc] at (-1.3,-2.35) {gluon};
\end{tikzpicture}
$$

## Fragmentation and hadronization

Between the hard parton and the observed hadrons lies **hadronization**, the
nonperturbative step in which color-charged partons convert into color-singlet
hadrons. It cannot be computed in perturbation theory, but its gross features follow
from the flux-tube picture. As the quark and antiquark separate, the tube between them
stretches and breaks by pair creation, again and again, until all the energy is stored
in low-momentum hadrons. Each break produces a new $q\bar q$ pair, and the fragments
combine into mesons and baryons streaming along the original parton direction.

The distribution of hadron momenta is summarized by a **fragmentation function**
$D_i^h(z)$, the probability that parton $i$ yields a hadron $h$ carrying a fraction $z$
of the parton momentum. Fragmentation functions are measured, not calculated, and
like the [parton distributions](/particle-physics/qcd/deep-inelastic-scattering-partons)
they evolve with scale by DGLAP equations. The guiding empirical fact is **local
parton-hadron duality**: the flow of hadrons closely tracks the flow of the underlying
partons, so the jet momentum reproduces the parton momentum up to soft corrections.

$$
% caption: From parton to hadrons. A hard quark radiates a cascade of gluons, which
% split into quark-antiquark pairs (the parton shower); at the confinement scale the
% partons combine into color-singlet hadrons flowing along the original direction.
% The perturbative shower is calculable; the final hadronization is not, but it
% preserves the parton's momentum and direction.
\begin{tikzpicture}[>=stealth, font=\footnotesize, scale=1.0]
  \definecolor{acc}{HTML}{4A6FA5}
  % initiating quark
  \draw[thick, ->] (-4.6,0) -- (-3.5,0);
  \draw[thick] (-3.5,0) -- (-2.8,0);
  \node[left, font=\scriptsize] at (-4.6,0) {quark};
  % branchings
  \draw[thick] (-2.8,0) -- (-1.6,0.6);
  \draw[thick] (-2.8,0) -- (-1.6,-0.5);
  \draw[thick] (-1.6,0.6) -- (-0.4,0.9);
  \draw[thick] (-1.6,0.6) -- (-0.4,0.3);
  \draw[thick] (-1.6,-0.5) -- (-0.4,-0.2);
  \draw[thick] (-1.6,-0.5) -- (-0.4,-0.8);
  \fill (-2.8,0) circle (1.4pt);
  \fill (-1.6,0.6) circle (1.4pt);
  \fill (-1.6,-0.5) circle (1.4pt);
  \node[black, font=\scriptsize] at (-1.6,-1.5) {parton shower};
  % hadronization boundary
  \draw[black, densely dashed] (0.3,-1.3) -- (0.3,1.3);
  \node[black, font=\scriptsize] at (0.3,1.55) {confinement scale};
  % hadrons
  \foreach \yy in {0.9,0.4,-0.05,-0.5,-0.95}
    {\draw[acc, thick, ->] (0.7,\yy) -- (2.1,\yy);
     \draw[acc, thick, fill=acc!12] (2.25,\yy) circle (2.6pt);}
  \node[black, font=\scriptsize] at (2.2,-1.5) {hadrons};
\end{tikzpicture}
$$

## Jet algorithms and event shapes

Turning a spray of hadrons into a definite jet requires a rule. A **jet algorithm**
clusters particles into jets by a distance measure; the modern sequential-recombination
algorithms merge the closest pair repeatedly until well-separated jets remain. The
requirement that a jet definition be **infrared and collinear safe** — unchanged when a
particle emits a soft gluon or splits collinearly — is what makes the jet cross section
calculable in perturbation theory.

An alternative to counting jets is to characterize the whole event by a continuous
**event-shape variable**. The most-used is the **thrust**,

$$
T = \max_{\hat n} \frac{\sum_i \lvert \vec p_i\cdot\hat n\rvert}{\sum_i \lvert \vec p_i\rvert},
$$

maximized over the axis $\hat n$. A perfect two-jet (pencil-like) event has $T \to 1$;
a spherical multi-jet event has $T \to \tfrac12$. The distribution of $T$ is calculable
in QCD, and its shape depends on $\alpha_s$: more gluon radiation broadens events and
pushes $T$ away from $1$.

$$
% caption: The thrust distribution. Most events are close to two-jet, piling up near
% thrust equal to one; gluon radiation produces the tail toward lower thrust. The
% shape of the tail is proportional to alpha_s, so fitting it measures the strong
% coupling. A larger coupling shifts weight toward broader, lower-thrust events.
\begin{tikzpicture}[font=\footnotesize, scale=1.0]
  \definecolor{acc}{HTML}{4A6FA5}
  \draw[->, black] (0,0) -- (8.4,0) node[right, font=\scriptsize] {thrust T};
  \draw[->, black] (0,0) -- (0,4.0) node[above, font=\scriptsize] {events};
  \node[below, black, font=\scriptsize] at (1.0,-0.05) {0.5};
  \node[below, black, font=\scriptsize] at (7.6,-0.05) {1.0};
  % distribution: rises toward T=1, sharp peak near right edge
  \draw[acc, very thick]
    (1.0,0.35) .. controls (3.0,0.6) and (5.0,1.0) .. (6.6,2.0)
    .. controls (7.1,2.7) and (7.35,3.4) .. (7.5,2.2)
    .. controls (7.55,1.4) and (7.6,0.6) .. (7.62,0.1);
  \node[black, font=\scriptsize, align=left] at (2.8,1.7) {gluon-radiation\\ tail};
  \draw[black, ->] (3.4,1.4) -- (4.6,0.9);
  \node[black, font=\scriptsize] at (6.9,3.1) {two-jet peak};
\end{tikzpicture}
$$

## Precision alpha_s

Jet rates and event shapes give some of the most precise determinations of the strong
coupling. The three-jet rate, the thrust distribution, and related observables are
each proportional to $\alpha_s$ at leading order, computed to higher orders in
perturbation theory, and fit to the data at the $Z$ mass. Combined with the
[deep-inelastic](/particle-physics/qcd/deep-inelastic-scattering-partons) and $\tau$-decay
determinations at other scales, they populate the running-coupling curve and yield the
world average

$$
\alpha_s\!\left(M_Z^2\right) = 0.1180 \pm 0.0009 .
$$

That a single coupling, extracted from jet shapes at $M_Z$, from DIS at a few GeV, and
from $\tau$ decays below that, all lie on the one asymptotic-freedom curve is the
quantitative triumph of QCD as the theory of the strong force.[^pdg-jets]

[^grif-jets]: Two-jet production in $e^+e^-$ annihilation, the $1 + \cos^2\theta$ angular distribution as evidence for spin-$\tfrac12$ quarks, and the jet picture of hadronization are treated in Griffiths, _Introduction to Elementary Particles_, §9.5, and Thomson, _Modern Particle Physics_, Ch. 10.

[^thom-gluon]: The discovery of three-jet events and the gluon at PETRA (TASSO, MARK-J, PLUTO, JADE) in 1979, and the gluon spin determination from three-jet angular correlations, are described in Thomson, Ch. 10, and Perkins, _Introduction to High Energy Physics_, Ch. 6. The relative rate $\sigma_{3\text{-jet}}/\sigma_{2\text{-jet}} \sim \alpha_s/\pi$ follows from the $q\bar q g$ matrix element.

[^pdg-jets]: The extraction of $\alpha_s$ from jet rates and event shapes, the thrust variable, and infrared-collinear safety are discussed in Thomson, Ch. 10, and Perkins, Ch. 6. The world average $\alpha_s(M_Z^2) = 0.1180 \pm 0.0009$ is from the Particle Data Group review of QCD, [pdg.lbl.gov](https://pdg.lbl.gov).