A TLAS-CONF-2019-035 25 July 2019
ATLAS CONF Note
ATLAS-CONF-2019-035
July 24, 2019
Measurement of the Lund jet plane using charged particles with the ATLAS detector from 13 TeV
proton–proton collisions
The ATLAS Collaboration
The prevalence of hadronic jets at the Large Hadron Collider (LHC) requires that a deep understanding of jet formation and structure must be achieved in order to reach the highest levels of experimental and theoretical precision. There have been many measurements of jet substructure at the LHC and previous colliders, but the targeted observables mix physical effects from various origins. Based on a new proposal to factorize physical effects, this note presents a double-differential cross section measurement of the Lund jet plane using 139 fb
−1of
√ s =13 TeV pp data collected with the ATLAS detector. The measurement uses charged particles to achieve a fine angular resolution and is corrected for acceptance and detector effects. Several parton shower Monte Carlo models are compared with the data.
Multiple parton shower Monte Carlo simulations are compared with data to study the modeling of various physical effects across the plane.
© 2019 CERN for the benefit of the ATLAS Collaboration.
Reproduction of this article or parts of it is allowed as specified in the CC-BY-4.0 license.
High energy quarks and gluons are ubiquitous in collisions at the Large Hadron Collider (LHC), but cannot be directly observed due to the confining nature of the strong force: instead, collimated jets of particles emerge from quark and gluon scattering. The internal substructure of a jet encodes information about the originating parton. Jet dynamics are governed by Quantum Chromodynamics (QCD) with one coupling constant called α
s. In the soft gluon (‘eikonal’) picture of jet formation, a quark or gluon radiates a haze of relatively low energy and statistically-independent gluons [1, 2]. QCD is nearly scale-invariant, and so this emission pattern is approximately uniform in the Lund plane [3]: a two-dimensional space spanned by ln ( 1 /z) and ln ( 1 /θ) , where z is the momentum fraction of the (relatively soft) emitted gluon relative to the primary quark or gluon and θ is the opening angle of the emission. This description can be extended to higher orders in QCD and has served as the basis for many calculations of jet substructure observables (see Refs. [4–7]).
The Lund plane is a powerful representation to provide insight and predictions for the entire radiation pattern within jets. However, because it is constructed from quarks and gluons within a particular theoretical framework, it is not observable due to the effects of confinement. A recent proposal [8] described an approach to approximate the radiation from a primary high-energy quark or gluon by declustering jets. Jets are defined in-part by a clustering algorithm; the most popular algorithms sequentially combine pairs of proto-jets to produce a set of jets from an initial set of jet constituents [9]. The proposal from Ref. [8]
prescribes the reclustering of a jet’s constituents using the Cambridge/Aachen (C/A) [10, 11] recombination algorithm, which imposes an angularly-ordered hierarchy on the jet’s clustering history. Then, by traversing the C/A clustering history in reverse while following the harder proto-jet (i.e., by following the primary clutering sequence), one can approximate the Lund plane by treating the proto-jet pair at each step in the sequence as a proxy for an emission (represented by the softer proto-jet) and the jet core (represented by the harder proto-jet) in the original theoretical depiction. For each proto-jet pair, at each step in the C/A declustering sequence, an entry may be made in the approximate Lund plane (henceforth, the ‘Lund jet plane’) by defining a coordinate pair of ‘Lund jet plane observables’ ( log ( 1 /z), log (R/∆R)) with
z = p
emissionT
p
emissionT
+ p
coreT
and ∆R
2= ( y
emission− y
core)
2+ (φ
emission− φ
core)
2,
where p
Tis the transverse momentum
1, y is the rapidity, and R is the radius parameter used during jet reconstruction. Using this approach, individual jets are represented as a set of points within the Lund jet plane, and ensembles of jets may be studied by measuring the double-differential cross section in this space.
The structure of emissions, which may themselves be composite objects, is not considered in this analysis.
To leading-logarithmic (LL) order, the average density of emissions within the primary Lund jet plane is uniform [8]:
1 N
jetsd
2N
emissionsdzd∆R ∝ 1
z 1
∆R ,
where n
jetsis the jet multiplicity. This picture of the plane in terms of z and ∆R is selected in order to cleanly factorise the momentum and angular measurements, although other depictions drawn in terms of e.g. k
t= zθ are also valid representations.
This note presents a two-dimensional differential cross section measurement of the Lund jet plane which is corrected for distortions due to detector effects, using 139 fb
−1of
√ s = 13 TeV pp collision data
1
ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector
and the z-axis along the beam pipe. The x-axis points from the IP to the center of the LHC ring, and the y -axis points
upward. Cylindrical coordinates ( r , φ ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The
pseudorapidity is defined in terms of the polar angle as η = − ln tan ( polar angle / 2 ) .
collected by the ATLAS detector during Run 2 of the LHC. In contrast to previous measurements of jet substructure, in the context of QCD, the radiation from different physical processes such as initial state radiation, underlying event (UE), hadronization, and perturbative emissions are factorized as illustrated in Figure 1. As different regions are dominated by factorised processes, the differential cross section of the Lund jet plane can be simultaneously useful for tuning non-perturbative models and for constraining advanced parton shower (PS) Monte Carlo programs [12–15] and state-of-the-art QCD calculations of jet substructure [16–21] which have so far only been constrained by the jet mass [22, 23] (which is itself a diagonal ray in the Lund jet plane: ln 1 /z ∼ ln m
2/p
2T
− 2 ln R/∆R ). Furthermore, the number of emissions within regions of the Lund jet plane is calculable within QCD and provides provably optimal discrimination between quark and gluon jets [5].
ln(R/∆R)
ln(1/ z )
har d &
wid e- angle
soft- collinear
hard- collinear UE/MPI
Regions of the Lund jet plane
P( ) ~ dz z d∆R ∆R
non -per
t. (s ma ll k t : z θ <
Λ QCD )
Figure 1: Schematic representation of the Lund jet plane, indicating the source of emissions which typically populate each region. Areas dominated by radiation in QCD originating from the underlying event (UE) and multi-parton interactions (MPI), soft-collinear splittings, hard-collinear splittings and hard, wide-angle radiation are indicated.
The ‘non-perturbative’ region lies above the diagonal band, and is populated by emissions with small k
t∼ zθ .
The ATLAS detector [24–26] is a general-purpose particle detector which provides nearly 4 π coverage in
solid angle. The inner tracking detector (ID) is inside a 2 T magnetic field and is designed to measure
charged-particle trajectories up to |η| = 2 . 5. The innermost component of the ID is a pixelated silicon
detector with fine granularity that is able to resolve ambiguities inside the dense hit environment of jet
cores [27], surrounded by silicon strip and transition radiation tracking detectors. Beyond the ID are
electromagnetic and hadronic calorimeters, from which topologically-connected clusters of cells [28] are clustered to form jets with the anti- k
talgorithm with radius parameter R = 0 . 4 [29, 30]. The energy scale of these jets is calibrated so that on average, the detector-level jet energy is the same as that of the corresponding particle-level jets [31]. The rate of jet production is too large to record every multijet event;
therefore, only jets with a large enough p
Tare saved as part of the ATLAS two-level trigger system [32, 33].
Events which are used in this measurement are selected online using a set of single-jet triggers which are fully efficient throughout the Run 2 dataset for jets used in this analysis. A well-balanced dijet system is required to be present, where the leading jet satisfies p
leadT
> 500 GeV and the system satisfies p
leadT
< 1 . 5 × p
subleadT
. Both jets must be within |η | < 2 . 1 in order to ensure that they are completely within the ID acceptance. Approximately 115 million reconstructed jets satisfy these selection criteria.
Particle-level charged hadrons and their reconstructed tracks are used for this measurement because individual particle trajectories can be precisely identified with the ID. As the Lund jet plane observables are dimensionless and isospin is an approximate symmetry of the strong force, the difference between the considered observables constructed using either all particles or charged particles following the analysis selections is small. ID tracks within ∆R = 0 . 4 of the axis of selected jets are used to construct the Lund jet plane observables: tracks associated with each jet are clustered using the C/A algorithm, and the plane is populated by iteratively declustering the clustering history. These tracks are required to have a p
T> 0 . 5 GeV, to be associated with the primary vertex of the event, and to satisfy “loose”
reconstruction quality criteria [34]. The measured Lund jet plane cross section is presented with 21 bins in ln ( 1 /z) between ln ( 1 / 0 . 5 ) and 9 . 4 × ln ( 1 / 0 . 5 ) and 15 bins in ln (R/∆R) between 0 . 0 and 5 . 0. The log (R/∆R) range is set to span the distance scale between the jet radius (log ( 0 . 4 / 0 . 4 ) = 0) and a minimum scale which is comparable to the pixel pitch (log ( 0 . 4 / 5 ) ∼ 50 microns / 33 mm ∼ 0 . 0015 ). The maximum z is 0.5 and the minimum is 500 MeV /p
jetT
, which is 0 . 001 for jets at the threshold of p
jetT
= 500 GeV.
A set of dijet events are simulated in order to perform the unfolding as well as to compare with the corrected data. The nominal set is simulated using Pythia 8.186 [35, 36] with the NNPDF2.3 LO [37] set of parton distribution functions (PDFs), a p
T-ordered PS, Lund string hadronization [38, 39], and the A14 set of tuned parameters (tune) [40]. Additional sets include: two sets of Pythia 8.230 [41] with the NNPDF2.3 LO PDF set and the A14 tune that are either using the Pythia LO matrix elements or NLO matrix elements from Powheg [42–45], Sherpa 2.1.1 [46] with the CT10LO PDF set, a p
T-ordered PS [47], and a cluster-based hadronization model [48], two sets of Sherpa 2.2.5 with the CT14NNLO PDF set [49]
and either a cluster hadronization model or the Lund string model, and two sets of Herwig 7.1.3 [15, 50, 51] with the MMHT2014NLO PDF set [52] and either the default angularly-ordered PS or a dipole PS and a cluster hadronization [48]. Additional details of these samples may be found in Ref. [53]. The Pythia 8.186 and Sherpa 2.1.1 events are passed through the ATLAS detector simulation [54] based on Geant 4 [55]. The effect of multiple pp interactions in the same and neighbouring bunch crossings (‘pileup’) is modelled in these samples by overlaying minimum-bias pp collisions, which are generated using Pythia 8 with the A3 set of tuned parameters [56] and NNPDF2.3 LO PDF set. The distribution of pileup is reweighted to match that of the data.
Selected data are unfolded in order to correct for detector bias, resolution, and acceptance effects using an Iterative Bayesian (IB) technique [57] and four iterations, as implemented in the RooUnfold framework [58].
The number of iterations was chosen to minimize the total uncertainty. The unfolding procedure corrects the
Lund jet plane constructed from detector-level objects to particle-level, where jets and charged particles are
defined similarly to detector-level: jets are clustered using the same anti- k
talgorithm with detector-stable
( cτ > 10 mm) hard-scatter hadrons as inputs. The same kinematic requirements are imposed on the
particle-level jet as for detector-level jets; charged particles with p
T> 0 . 5 GeV within ∆R < 0 . 4 of the particle-level jets are used to populate the particle-level Lund jet plane.
Emissions at detector-level and particle-level are matched in η − φ to construct the response matrix. This matching is unique, with the closest match within ∆R = 0 . 1 taking precedence. Unmatched emissions from tracks not due to a single charged particle (detector-level) and from non-reconstructed charged particles (particle-level) are accounted for with acceptance and efficiency corrections, before and after the response matrix portion of the IB method. These corrections are about 20% for the lower left quadrant of the plane and increase to 90% in the most collinear bins and to 50% in the lowest z bins. Despite the large corrections, the correction factors are robust to all systematic variations within 10%, as described below. The matching procedure follows the order of the C/A declustering, starting from the widest-angle emissions and iterating towards the jet core. An uncertainty on this matching procedure is established by repeating the unfolding procedure and iterating through the C/A declustering sequence in reverse, and taking the difference in the unfolded result as an uncertainty. This uncertainty is less than 1% except in the highest log ( R/∆R) bin.
For matched emissions, the log ( 1 /z) and log (R/∆R) migrations are largely independent of each other.
Furthermore, since the differential distribution is slowly varying across the Lund jet plane, the purity and efficiencies are approximately the same. The log (R/∆R) migrations in a given log ( 1 /z) bin are about 50%
for the smallest opening angles and decrease to under 10% for the widest angles. The log ( 1 /z) migrations decrease from about 50% for the softest emissions to about 20% for the hardest ones, with some degradation for the softest emissions at small opening angles. Migrations for both observables are nearly symmetric except for log (R/∆R) & 3, where harder to resolve small opening angles are asymmetrically measured as larger ones and . 10% of the time emissions are mis-matched between particle-level and detector-level and therefore measured with the wrong log ( 1 /z) . While the log (R/∆R) migrations are nearly the same when log ( 1 /z) migrates by one bin, the log ( 1 /z) migrations increase by about 30% when log (R/∆R) migrates by one bin.
The unfolded distribution is normalized by the number of jets that pass the event selection, and so this measurement is made insensitive to the total jet cross section. Following normalization, the integral of the Lund jet plane is the average number of emissions within the fiducial region studied.
Systematic uncertainties on the measured cross section distribution may be broadly classified as one of two types: experimental uncertainties related to the measurement and reconstruction of the objects used to perform this measurement, and theoretical uncertainties related to the modeling of jet fragmentation.
These uncertainties are evaluated by modifying some aspect of the unfolding procedure and taking the difference with respect to the nominal result.
Experimental uncertainties from various sources are considered and found to be comparable to those arising from theoretical sources in some regions of the Lund jet plane. Uncertainties on the jet energy scale and resolution are measured using several simulation-based and in situ techniques [31]. These uncertainties cause the migration of jets into or out of the fiducial acceptance which causes an uncertainty on the measured observables to typically be above 3%, reaching up to 7% in certain regions of the plane.
Uncertainties related to the track reconstruction both for isolated tracks and those within dense environments are considered, which modify the measured p
Tof individual tracks or remove them completely [27, 59].
These uncertainties are typically small, contributing less than a half-percent to the total uncertainty.
A data-driven non-closure uncertainty is determined by unfolding the nominal particle-level spectrum following a reweighting based on the corresponding simulated detector-level spectrum and the data [60].
This uncertainty is less than 1% except for the smallest opening angles, where it can reach 5%. Other
experimental uncertainties related to the stability of the measurement to the effects of pileup are considered, and are found to be negligible.
Theoretical uncertainties predominantly arise due to the limited accuracy of jet fragmentation modelling.
Variations in jet fragmentation can impact the response matrix as well as the unfolding correction factors.
The impact of this uncertainty on the measurement is evaluated by repeating the unfolding procedure with variations of the Monte Carlo generators used. As Pythia 8.186 is used to perform the nominal unfolding, the impact of constructing the unfolding matrix using simulated Sherpa 2.2.1 dijet events is evaluated. This alternative matrix is used to unfold the nominal reconstructed Pythia 8.186 sample, and differences between these results and the nominal particle-level Pythia 8.186 Lund jet plane are taken as an uncertainty, which is typically between 1–10% depending on the considered region (larger in the soft-collinear) and is the largest single source of uncertainty associated with this measurement.
The measured Lund jet plane following the unfolding procedure is shown in Figure 2. A triangular region is populated uniformly by perturbative emissions within the ensemble of jets. A large number of emissions can be found at the transition to the non-perturbative regime, as α
s( k
t) is enhanced for small values of k
t. Non-perturbative emissions beyond this region, for which k
t∼ Λ
QCD, are suppressed. The average number of emissions in the measured region (integral of the Lund jet plane) is measured to be 7.75 emissions while the corresponding value for Pythia 8.186 is 8.13 emissions and for Sherpa 2.2.1 is 7.51 emissions. A total uncertainty of ± 0 . 11 is associated with this measured value, found by propagating all systematic uncertainties from the unfolded measurement in an uncorrelated and symmetrized manner. The associated statistical uncertainty on this value is negligible. While this bracketing of the data by Pythia 8.186 and Sherpa 2.2.1 in the particle multiplicity inside jets has been observed before [61], it has never been decomposed into a perturbative and non-perturbative component.
Figure 3 shows data from four selected horizontal and vertical slices through the Lund jet plane. These data are compared to predictions from multiple Monte Carlo generators. While no one prediction describes the data in all regions, the Herwig 7.1.3 angularly-ordered prediction provides the best description across most of the plane. The Pythia 8.186 and Sherpa 2.2.1 predictions are not shown, but are consistent with the Pythia 8.230 and Sherpa 2.2.5 (Cluster Had.) predictions, respectively. The differences between the shower algorithms implemented in Herwig 7.1.3 are present only at low values of z∆R while the differences between the hadronization algorithms implemented in Sherpa 2.2.5 are present mostly at large z∆R . Similarly the Pythia and Sherpa +string predictions are nearly the same in the high z∆R region. Since the Lund jet plane is reconstructed for small-radius jets, there is no difference between the Powheg+Pythia and Pythia predictions. These observations demonstrate that the Lund jet plane is able to isolate various effects and can therefore provide useful input to both perturbative and non-perturbative model development and tuning.
In summary, a measurement of the substructure of jets initiated by quarks and gluons as represented by
the Lund jet plane is reported, which was performed in 139 fb
−1of LHC pp collisions recorded by the
ATLAS detector during Run 2 of the LHC. This measurement is performed on an inclusive selection of dijet
events, with a leading jet p
T> 500 GeV. Selected jets are reconstructed from topological clusters using
the anti- k
talgorithm with R = 0 . 4, and their associated charged-particle tracks are used to construct the
observables of interest. The measured data are presented as an unfolded double-differential cross section,
and compared to several Monte Carlo generators. The average number of emissions entering the region
of the Lund jet plane is measured to be 7 . 75 ± 0 . 11. The measured data disagree with predictions from
simulation throughout the measured parameter space, and are a first demonstration of the Lund jet plane’s
potential to isolate various effects and to provide useful input to both perturbative and non-perturbative
model development and tuning.
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 2
) R ∆ / R ) d ln( z d ln(1/ / N d 1/N
emissionsjets0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆ R / R ln(
1 2 3 4 5 6
) z ln(1/
ATLAS Preliminary s = 13 TeV, 139 fb
-1−2
10
−1
10 )
core Tp +
emission Tp / (
emission Tp = z
−2 1
10
10
−(emission, core)
∆ R =
∆ R
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
Relative Uncertainty
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆ R / R ln(
1 2 3 4 5 6
) z ln(1/
ATLAS Preliminary s = 13 TeV, 139 fb
-1−2
10
−1
10 )
core Tp +
emission Tp / (
emission Tp = z
−2 1
10
10
−(emission, core)
∆ R =
∆ R
Figure 2: Top: The Lund jet plane as measured using jets in 13 TeV pp collision data, corrected to particle-level.
The inner set of axes indicate the coordinates of the Lund jet plane itself, while the outer set indicate corresponding
values of z and ∆R . Bottom: The total relative uncertainty (experimental, statistical and related to Monte Carlo
modeling effects) as a function of the Lund jet plane observables.
1 2 3 4 5 6 ) z ln(1/
0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 )z / d ln(1/emissions) dNjets(1/N
ATLAS Preliminary = 13 TeV, 139 fb-1
s
) < 1.00
∆R / R 0.67 < ln(
Data Pythia 8.230 Powheg + Pythia 8.230 Sherpa 2.2.5 (Cluster Had.) Sherpa 2.2.5 (String Had.) Herwig 7.1.3 (Dipole Shower) Herwig 7.1.3 (Angular Shower)
1 2 3 4 5 6
) z ln(1/
0.6 0.81 1.2 1.4
Ratio to Data
(a)
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆R / R ln(
0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 )R∆/R / d ln(emissions) dNjets(1/N
ATLAS Preliminary = 13 TeV, 139 fb-1
s
) < 1.25 z 0.97 < ln(1/
Data Pythia 8.230 Powheg + Pythia 8.230 Sherpa 2.2.5 (Cluster Had.) Sherpa 2.2.5 (String Had.) Herwig 7.1.3 (Dipole Shower) Herwig 7.1.3 (Angular Shower)
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆R / R ln(
0.6 0.81 1.2 1.4
Ratio to Data
(b)
1 2 3 4 5 6
) z ln(1/
0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 )z / d ln(1/emissions) dNjets(1/N
ATLAS Preliminary = 13 TeV, 139 fb-1
s
) < 3.67
∆R / R 3.33 < ln(
Data Pythia 8.230 Powheg + Pythia 8.230 Sherpa 2.2.5 (Cluster Had.) Sherpa 2.2.5 (String Had.) Herwig 7.1.3 (Dipole Shower) Herwig 7.1.3 (Angular Shower)
1 2 3 4 5 6
) z ln(1/
0.6 0.81 1.2 1.4
Ratio to Data
(c)
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆R / R ln(
0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 )R∆/R / d ln(emissions) dNjets(1/N
ATLAS Preliminary = 13 TeV, 139 fb-1
s
) < 3.47 z 3.19 < ln(1/
Data Pythia 8.230 Powheg + Pythia 8.230 Sherpa 2.2.5 (Cluster Had.) Sherpa 2.2.5 (String Had.) Herwig 7.1.3 (Dipole Shower) Herwig 7.1.3 (Angular Shower)
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆R / R ln(
0.6 0.81 1.2 1.4
Ratio to Data
(d)
Figure 3: Representative horizontal and vertical slices through the Lund jet plane. Unfolded data are compared to
particle-level simulation from several Monte Carlo generators. The uncertainty band includes all sources of systematic
and statistical uncertainty. The inset triangle illustrates which slice of the plane is depicted: 0 . 67 < ln (R/∆R) < 1 . 00
(a), 0 . 97 < ln ( 1 /z) < 1 . 25 (b), 3 . 33 < ln (R/∆R) < 3 . 67 (c), and 3 . 19 < ln ( 1 /z) < 3 . 47 (d).
Appendix
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 2 ) R ∆ / R ) d ln( z d ln(1/ / N d 1/N emissions jets
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆ R / R ln(
1 2 3 4 5 6
) z ln(1/
ATLAS Simulation Preliminary s = 13 TeV, 139 fb
-1−2
10
−1
10 ) core T p + emission T p / ( emission T p = z
−2 1
10
10
−(emission, core)
∆ R =
∆ R
Figure 4: The Lund jet plane as simulated by the Pythia v8.186 Monte Carlo generator. corrected to particle-level.
The inner set of axes indicate the coordinates of the Lund jet plane itself, while the outer set indicate corresponding
values of z and ∆R .
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08
) R ∆ / R ) d ln( z d ln(1/ / emissions N 2 d jets 1/N
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
)
∆ R / R ln(
1 2 3 4 5 6
) z ln(1/
ATLAS Simulation Preliminary s = 13 TeV, 139 fb
-1−2
10
−1
10 ) core T p + emission T p / ( emission T p = z
−2 1