Power Counting Energy Flow Polynomials

The Big Picture

Imagine trying to identify a person in a crowd using every possible facial measurement: the distance between each pair of freckles, every angle of every feature. You’d drown in data. Now imagine a physicist’s version of that problem, applied to the chaos of particle collisions at the Large Hadron Collider.

Every time protons smash together at CERN, the resulting debris sprays into tight, narrow streams of particles called jets. The internal patterns of these jets, their substructure, carry fingerprints of the fundamental interactions that created them.

For decades, physicists have built specialized measurements to read those fingerprints. More recently, the field has pushed toward complete mathematical toolkits: collections of measurements that, in principle, capture everything about a jet. One such toolkit, energy flow polynomials (EFPs), is mathematically guaranteed to form an overcomplete linear basis for jet observables. Think of EFPs as a systematically organized family of measurements, each probing a different combination of angular relationships between particles in the jet, arranged by their complexity.

Completeness comes at a steep price, though. By degree six (the sixth level of complexity) there are 314 distinct EFPs. Using all of them is computationally painful and conceptually murky.

This paper by Pedro Cal, Jesse Thaler, and Wouter J. Waalewijn cuts through that excess. Using a classical physics technique called power counting, they show that the EFP toolkit contains enormous redundancy for realistic jets, and that a dramatically smaller subset performs just as well.

Key Insight: By applying power counting arguments appropriate for quark and gluon jets, the authors reduce the EFP basis by more than a factor of six at degree six (from 212 down to 22 elements), without sacrificing machine learning performance.

How It Works

EFPs rest on an elegant idea: represent each observable as a graph. Each node corresponds to a particle in the jet; each edge encodes an angular correlation between particles. The degree of an EFP (the number of edges) controls how many angular correlations it probes. Computing all EFPs up to degree six means evaluating 314 such graphs over every particle pair, triplet, or quartet in a jet. Computational cost scales steeply with graph complexity.

Power counting is a tool borrowed from effective field theory, the physics framework that describes complex systems by focusing only on effects that matter at a given scale and systematically ignoring the rest. Not all contributions are equally important. By identifying what’s “big” and what’s “small” in typical jet configurations, you can organize observables by relative importance and discard subleading terms. The authors apply three power counting schemes:

• Strongly-ordered expansion: Assumes jet emissions are hierarchically ordered in both energy and angle, a realistic approximation for leading-logarithmic calculations.
• 1-collinear expansion: Keeps one energetic collinear emission (a particle flying nearly parallel to the jet axis) with full angular information, expanding around that configuration.
• 2-collinear expansion: Extends to two energetic collinear emissions, capturing next-to-leading logarithmic physics at a finer level of precision.

Within each scheme, power counting reveals linear relationships between EFPs. Two graphically distinct EFPs may be numerically equivalent, to the precision of the approximation, once you account for how jets actually form. The authors derive these relations analytically, then test them on millions of simulated jets from Pythia, a widely used program that simulates particle production and scattering after high-energy collisions.

The numbers speak for themselves. For the strongly-ordered basis, the degree-six EFP count drops from 212 to just 22. The 2-collinear basis needs 37 elements, richer because it retains more angular information, but still tiny compared to the full set.

The predicted linear relations hold up in Pythia with excellent numerical agreement.

Why It Matters

The immediate payoff is practical. The authors run quark/gluon tagging, distinguishing jets initiated by quarks from those initiated by gluons, using logistic regression on both the full EFP set and their reduced bases. The reduced bases match full-set performance using a fraction of the inputs. For regression tasks predicting continuous jet properties, the story is the same.

There’s a deeper gain here, too. EFPs were always intended as a complete basis, but completeness without structure is just a list. Power counting gives that list a hierarchy: it tells you which observables matter most and for which physics.

Machine learning alone can’t provide that kind of understanding. A neural network might implicitly learn to ignore redundant EFPs; power counting explains why they’re redundant and when that redundancy breaks down.

On the computational side, there’s a bonus in the 1-collinear case. Naively, computing an N-point EFP on M particles costs O(M^N), so the computation balloons as both particle count and observable complexity grow. Power counting lets the authors “cut open” high-complexity graphs and express them as products of simpler ones, reducing the tree-width (a graph-theory measure of computational complexity) of the underlying computation. Fewer floating-point operations per jet matters at the LHC, where jets are reconstructed millions of times per second.

The approach raises natural follow-up questions. Can power counting extend to multi-prong jets from top quarks or W bosons? Can the basis reductions inform better neural network architectures for jet physics? The authors focus on single-prong quark and gluon jets here and flag multi-prong extension as explicit future work.

Bottom Line: Power counting transforms the EFP basis from an unwieldy 314-element list into a compact, physically motivated toolkit of 22 to 37 elements, enabling faster computation, cleaner interpretability, and undiminished machine learning performance on jet classification tasks.