Jet Substructure at The Large Hadron Collider: 10.1103/revmodphys.91.045003
Jet Substructure at The Large Hadron Collider: 10.1103/revmodphys.91.045003
Mario Campanelli
University College London,
UK
Chris Delitzsch
University of Arizona,
USA
Philip Harris
Massachusetts Institute of Technology,
USA
Deepak Kar
University of Witwatersrand,
South Africa
Yuta Takahashi
Universität Zürich,
Switzerland
Marcel Vos
IFIC Valencia, Spain
Jet substructure has emerged to play a central role at the Large Hadron Collider (LHC),
where it has provided numerous innovative ways to search for new physics and to probe
the Standard Model, particularly in extreme regions of phase space. In this article
we focus on a review of the development and use of state-of-the-art jet substructure
techniques by the ATLAS and CMS experiments.
CONTENTS A. Inputs 6
Contents 1 B. Calibration 7
2
1 Sometimes the rapidity (y) is used and sometimes the pseudo- 2 For a theoretical introduction to jets, we recommend the reviews
rapidity (η) is used depending on the application. See Ref. (28) in Refs. (32; 33) as well as the theory companion this experimen-
for a detailed discussion. tal review, Ref. (16).
4
compensating nature of hadron calorimeters, suppression ployed in ATLAS and CMS analyses are discussed in
of electronic noise, tracking inefficiencies, dead material section IV. In section V we review jet grooming tech-
in front of calorimeters, the influence of pile-up and other niques in use in experimental analyses and discuss their
effects. While the calibration of the total jet energy impact on jet substructure observables. A special em-
scale is an important aspect in all analyses using jets, phasis is given on the jet mass calibration and jet mass
the precise knowledge of the jet mass scale and the de- measurements in different final states. Measurements of
tector response to jet substructure observables and jet other jet substructure distributions are described as well.
tagging algorithms is specific to jet substructure analy- One of the key developments within the field of jet sub-
ses. Calibrating the jet energy scale results in a change structure are tagging algorithms, which are described in
of the magnitude of the jet’s four-momentum, where the detail in section VI. Theoretical and experimental devel-
jet mass scale comprises an additional degree of freedom opments have resulted in large performance gains of sub-
that can not be constrained by the typical methods of structure taggers in the last years, relevant for a large
balancing a jet with a well-calibrated reference object. number of present and future physics analyses. We high-
The jet mass scale is usually calibrated using jets from light the main developments and improvements and give
fully-merged, highly boosted W → q q̄ 0 decays, facilitat- an overview of relevant experimental studies. The use of
ing a calibration of the peak position in the jet mass jet substructure taggers in existing cross section measure-
distribution. Measurements of the jet mass distribution ments is reviewed in section VII. So far, the major ben-
from light quark and gluon jets, as well as from fully- eficiaries of jet substructure methods have been analyses
hadronic highly-boosted W , Z and t decays allow for in search for new physical phenomena. We review the
precise tests of the modelling of perturbative and non- application of these methods to searches for new physics
perturbative effects in jet production. Similar measure- in section VIII and conclude in section IX.
ments can also be used to study the detector response to
jet substructure observables and their modelling in sim-
ulation. A mis-modelling of variables used for tagging, II. ATLAS AND CMS DETECTORS
either in the detector simulation or on the level of the
underlying physics, can result in a wrong estimation of The ATLAS (45) and CMS (46) detectors are designed
the tagging efficiency or the misidentification rate, with to observe leptons, photons, and hadrons resulting from
important consequences for measurements. In order to LHC pp collisions. The physics of the hard reaction
overcome this limitation, measurements of tagging effi- takes place at the point of collision (the primary ver-
ciencies and misidentification rates are performed in sam- tex) within the beam pipe. Beyond the beam pipe4 ,
ples enriched with the particle decays in question. While at 4.4 cm (3.3 cm) in CMS (ATLAS), the first cylindri-
these measurements do not help to understand the cause cal layer of detectors encountered are silicon pixels and
of the mis-modelling or to improve the description of jet strips for identification of charged particles. CMS pro-
substructure distributions, these can be used to correct vides a 3.8 T magnetic field via a solenoid positioned
the efficiencies in simulation. It is these measurements outside the silicon tracking detector, the Electromag-
that have enabled the use of jet substructure taggers in netic Calorimeter (ECAL) and most of the Hadronic
numerous physics analyses since the beginning of data Calorimeter (HCAL). ATLAS has an additional tracking
taking at the LHC. The increased statistics from a data layer composed of straw drift tubes (Transition Radia-
sample corresponding to about 150 fb−1 per experiment tion Tracking or TRT), with a 2 T magnetic field encom-
at a centre-of-mass energy of 13 TeV can now be used to passing the silicon and TRT detectors, while the ECAL
improve our understanding of the detector response to and HCAL are situated outside the solenoidal magnet.
jet substructure algorithms, the underlying physics and The calorimeters are surrounded by muon spectrometers
the performance differences of taggers. These studies and which build the outermost part of the ATLAS and CMS
measurements represent the continuation of an exciting detectors. Both detectors are nearly hermetic and can
physics program at the LHC in a field which reached therefore measure the missing transverse momentum.
its adolescence in the past few years. In the years to The energy and momentum ranges and resolutions for
come, the field of jet substructure will evolve and mature the barrel regions5 of ATLAS and CMS are shown in
through precision measurements and the exploration of table I along with the measurement granularity, which
unknown territory.
We begin this review with a brief overview of the AT-
LAS and CMS detectors in section II, followed by a de- 4 The LHC collaborations are continuously working to improve
scription of the input to jet reconstruction and jet cal- the detectors; the numbers given here are for the detectors that
operated in 2015-2017. Before and after this time, the exact
ibration in section III. An important aspect of jet re-
values are not the same as reported here.
construction at the LHC, and jet substructure in par- 5 For example, the ATLAS ECAL barrel covers the pseudorapidity
ticular, are algorithms to mitigate the effects of pile-up. range |η| < 1.475, the end-caps cover 1.375 < |η| < 3.2 and the
Recent experimental advancements and algorithms em- forward ECAL layer extends the coverage up to |η| < 4.9. The
6
limits the angular resolution. The better energy res- ATLAS CMS
olution of the CMS ECAL is due to the use of lead
tungstate (PbWO4 ) crystals, as opposed to the Liquid Tracking
Argon (LAr) used by ATLAS. The differences in the AT- 1/pT resolution 0.05% × 0.02%×pT / GeV⊕
LAS and CMS calorimeter designs are a result of the pT / GeV⊕1% (47) 0.8% (48)
different ranking of priorities decided by the two collab- d0 resolution 20 (49) 20 (48)
orations; ATLAS chose a radiation-hard technology with (µm)
sufficient resolution in a fine sampling LAr calorimeter,
ECAL
while CMS prioritized the excellent resolution of a total √ √
absorption crystal calorimeter (the focus was Higgs mass E resolution 10%/ E ⊕ 3%/ E ⊕
reconstruction), and accepted the accompanying limita- 0.2% (45) 12%/E ⊕
0.3% (46)
tions in radiation-hardness associated with this technol-
ogy. The CMS ECAL crystal response varies under ir- granularity 0.025 × 0.025 0.017 × 0.017
radiation, which is partially recovered in a few hours at HCAL
room temperature. √ √
E resolution 50%/ E ⊕ 100%/ E ⊕
The ATLAS ECAL is segmented into three (two) longi- 5% (45) 5% (50)
tudinal layers for |η| < 2.5 (|η| > 2.5). The granularity of
the ATLAS ECAL in table I refers to its second layer (as granularity 0.1 × 0.1 0.087 × 0.087
most of the electromagnetic energy is deposited there);
the first layer has a finer granularity in η. The multi- TABLE I ATLAS and CMS detectors in the barrel regions.
ple layers allow for a finer granularity than the cell size The granularity is in pseudorapidity and azimuth (η × φ) and
d0 is the transverse impact parameter resolution with respect
in any of the individual layers, being advantageous over to the beam-line. The tracker momentum resolution is from
a laterally segmented calorimeter, and additionally pro- muons while the d0 resolution is from generic charged particles
vide pointing information. The difference between AT- (mostly pions) in tt̄ events. The ECAL energy resolution
LAS and CMS for the HCAL resolution is particularly is presented for electrons. The granularity for the ATLAS
large at higher energies: a 1 TeV jet has σ(E) E ∼ 2% in calorimeters are for the middle layers only, which collect the
σ(E) largest amount of energy. For the ATLAS EM calorimeter,
ATLAS, in contrast to E ∼ 5% in CMS. This is one the innermost layer has ∆η = 0.0031 for γ/π 0 separation.
reason why CMS fully adapted a particle flow technique
since the beginning of the LHC (see section III.A below).
III. JET RECONSTRUCTION ation of particle-flow (54), most of this review will focus
on calorimeter-only jets as they are still the most widely
A. Inputs used setup. ATLAS benefits less than CMS from particle
flow because of its weaker magnetic field and longitudi-
Both experiments have dedicated algorithms to recon- nally segmented calorimeter.
struct particle kinematics from calorimeter and tracker ATLAS and CMS combine calorimeter cells using topo-
information designed to minimize the fake rate, max- logical clusters (53; 55). These clusters are three dimen-
imize the efficiency, and minimize the bias and reso- sional in ATLAS as a result of the longitudinal segmen-
lution of the particle candidate parameters. As there tation. Cluster seeds are started from highly significant
is no algorithm that can simultaneously optimize all of energy (high cell signal to average electronic ⊕ pileup
these objectives, the various approaches trade off op- noise) deposits which are combined (or split) based on
timality under one metric for improvements under an- the distribution of the significance of energy in nearby
other. ATLAS and CMS have also developed different cells. Calorimeter-cell clusters in CMS are obtained us-
algorithms that cater to the experiment’s hardware as ing a Gaussian-mixture model, which results in one or
well as the collaboration’s goals for the tradeoffs. By more calorimeter clusters within each topological cluster.
default, CMS combines tracker and calorimeter informa- HCAL clusters can be split according to the number and
tion into unified particle flow objects as inputs to jet re- energy distribution of associated ECAL clusters. Clus-
construction (51; 52; 53). ATLAS has traditionally used ter splitting is critical to achieve a better estimate of the
calorimeter-only information for jet reconstruction, with spatial energy distribution as input to jet substructure
tracking information used to augment/enhance the per- algorithms (56; 57).
formance. While ATLAS is current migrating to a vari- The topological clusters are calibrated using simula-
tions to account for the non-compensating calorimeter re-
sponse to hadrons, signal losses due to energy deposited
CMS ECAL barrel covers |η| < 1.48, the end-caps extend the in inactive detector material and signal losses on clus-
coverage up to |η| < 3. ter boundaries caused by the topological clustering al-
7
Energy resolution
classification of clusters as hadronic or electromagnetic
0.6 CMS Anti-kT, R = 0.4 Calo
in origin based on the energy and position of the clus-
ter, the longitudinal depth (λclus ) and normalized signal Simulation |ηRef| < 1.3 PF
energy density; hadronic showers tend to occur deeper
in the calorimeter and be less dense (55). Charged and
neutral pions are used to derive this classification and 0.4
calibration, called the Local Cell Weighting (LCW). In
CMS, dedicated ECAL (based on photons) and HCAL
(based on neutral kaons) calibrations are combined to
account for energy and |η|-dependent non-linearities in
the hadron calorimeter response (53). Both ATLAS and 0.2
CMS validate the performance of these calibrations with
single particle studies in data (53; 58).
Different strategies are used by ATLAS and CMS to
reconstruct tracks from their inner detectors. ATLAS fo- 0
cuses first on maintaining a high efficiency with a rather
20 100 200 1000
inclusive first pass through inner detector hits. A second pRef (GeV)
T
step known as ambiguity solving reduces the fake rate. In
contrast, CMS uses a sequential approach with multiple FIG. 3 Jet energy resolution for particle flow (red, lower line)
passes through the remaining inner detector hits. With and calorimeter-only (blue, upper line) jets in the barrel re-
each pass, the efficiency increases while maintaining a gion in CMS simulation, with no pile-up, as a function of the
pT of the reference jet. Taken from (53).
low fake rate. Both procedures are effective at identi-
fying about 90% of charged pions above 1 GeV with a
percent-level (or smaller) fake rate. Lower momentum
particles can be reconstructed, at the cost of a higher ter envelope in the less granular HCAL. Tracks with a
fake rate and lower efficiency. Due to its weaker mag- pT uncertainty in excess of the calorimetric energy reso-
netic field, ATLAS is able to reach low track momentum lution expected for charged hadrons are masked, which
of 100 MeV for physics analysis (59), although most jet allows the rate of misreconstructed tracks at large pT to
substructure measurements and searches use a threshold be reduced.
of 500 MeV. In contrast, the momentum resolution in
The ECAL and HCAL clusters not linked to any track
CMS is excellent up to higher momenta than in ATLAS.
give rise to photons and neutral hadrons. Charged
The TRT can be used to improve the momentum reso-
hadrons are created from the remaining ECAL and
lution of high pT tracks (60), but the weaker magnetic
HCAL clusters, linked to tracks. If the calibrated calori-
field despite a comparable inner detector radius is a fun-
metric energy is compatible with the corresponding track
damental limitation.
momenta under the charged-pion hypothesis, no neutral
Both experiments have implemented dedicated strate- particles are created. Otherwise, the excess energy is in-
gies for track reconstruction in high density environments terpreted to originate from photons and neutral hadrons
such as the core of high pT jets. In such environments, for deposits in the ECAL and HCAL, respectively. The
pixel and strip clusters can merge resulting in a loss in particle flow algorithm in ATLAS is similar to the one
tracking efficiency and degraded resolution. ATLAS has used by CMS and is described in more detail in Ref. (54).
implemented a stacked neural network (NN) approach to
examine pixel clusters to identify multi-particle clusters, The combination of tracking and calorimetric measure-
estimate the position of the particles passing through ments results in an optimal input for jet substructure
the clusters, and also predict the residual resolution of measurements, making use of the superior angular reso-
the position estimates (61; 62; 63; 64; 65). CMS has lution from the tracking detector and calibrated calorime-
introduced a dedicated tracking step in which a cluster ter clusters. Once the calibrated PF objects are clustered
splitting procedure attempts to split merged clusters ex- into jets, their relative momenta and angular distances
ploiting the information of the jet direction, predicting are kept constant, and only the total energy response
the expected cluster shape and charge (66). of jets is corrected with factorized JES calibrations (see
For particle flow in CMS, tracks and calibrated clusters section III.B).
are combined taking the tracking and calorimeter resolu- The particle flow algorithm improves the energy reso-
tions into account. First, a link is created between tracks lution as shown in figure 3. A similar performance gain is
in the central tracker and calorimeter clusters. Links are observed in ATLAS, but the weaker magnetic field means
also created between clusters in the ECAL and HCAL, that the point where calorimetery and tracking are com-
when the cluster position in the ECAL is within the clus- parable is lower (about 100 GeV).
8
0.1
Fractional JES uncertainty
19.7 fb-1 (8 TeV)
Data 2015, s = 13 TeV ATLAS 6
0 2
20 30 40 102 2×102 103 2×103
jet
p [GeV]
T
1
FIG. 4 ATLAS jet energy scale uncertainty. Adapted
from (68).
0
20 100 200 1000
p (GeV)
B. Calibration T
the whole jet, the individual particles within the jet can
I
19.7 fb−1 (8 TeV) be classified to whether they belong to the actual jet or to
ρ (GeV) 30 the underlying pile-up. Charged particles leave tracks in
CMS high granularity tracking detectors at the heart of multi-
25 purpose detectors like ATLAS and CMS and can be sep-
ρ(µ) = 1.00 + 0.54 µ + 0.0005 µ2 arated based on their longitudinal position ẑ (along the
ρ(µ) = 1.01 + 0.50 µ + 0.0007 µ2
beamline) within the luminous region (see figure 7). The
20 charged hadron subtraction (CHS) (74) method identi-
fies each pile-up track individually. Used in concert with
15 particle flow concepts which attempt to identify each par-
ticle in the event uniquely, CHS can effectively remove
all charged pile-up radiation from the event, including
10 calorimeter signals that are linked to tracks through the
particle flow algorithm. Identification of pile-up jets,
formed predominantly from the energy of one or many
5
Data pile-up vertices, is another technique for removing pile-
MC up using charged particles; by determining the fraction
0 of energy of the jet from the primary vertex, one can
0 10 20 30 40 distinguish such pile-up jets from the PV jets (75; 77).
µ The two methods discussed above can be combined.
First the more precise CHS method subtracts the pile-
FIG. 6 Average pile-up contribution to the jet pT , ρ, as a up contribution from charged particles; in a second step,
function of the average number of pile-up interactions per the remaining contributions from neutral particles are
bunch crossing, µ, for data (circles) and simulation (dia-
monds) at the CMS experiment. Taken from (44).
removed with the area subtraction method.
In a more advanced approach, local, topological infor-
mation is used, as QCD radiation from pile-up vertices is
even higher values in Run 3, and culminating at the high often uncorrelated and soft. It and can thus be removed
luminosity LHC (HL-LHC) reaching up to hµi = 140 − based on the local energy profile, i.e. if the radiation is
200. not consistent with hard scattering radiation from the
Pile-up typically leaves about 0.5 GeV of energy in the PV. This can be done in the transverse plane η, φ and
detector per unit area (η, φ), per pile-up vertex; the ef- also as a function of radiation depth. The jet grooming
fects of this are present in all aspects of LHC physics, technique is such an example to clean the jet of soft and
from detector design and software performance to the fi- wide-angle radiation which incidentally removes pile-up
nal sensitivity of measurements and searches. radiation. It is discussed in more detail in section V.A.
Topoclustering (55), used by the ATLAS Collaboration,
is deployed at the formation of clusters in the calorimeter
B. Mitigation Methods requiring radiation to have a certain topological profile.
In the forward region, where no tracking information is
Properties of pile-up interactions are exploited to dis- available, jet shapes and topological correlations can be
criminate pile-up particles from particles originating from used to identify pile-up (78).
the primary vertex, or to remove energy contributions While the above methods have been successfully de-
from pile-up to the individual jet. ployed in the LHC experiments, they each have some
Pile-up can be approximated as a spatially uniform de- deficiencies as well; ideally, one would hope to effectively
position of energy. The so-called area subtraction uses a combine all pile-up mitigation handles in order to maxi-
pile-up pT density per unit area estimator, ρ, and defines mally distinguish pile-up from PV radiation and to re-
a jet catchment area, A, to remove energy that is assumed move pile-up at the most granular level possible, i.e.
to originate from pile-up interaction. This approach cor- at the particle or constituent level, in order to be as
rects the jet in the following way: pcorrT = porig
T − ρA. generic as possible. For example, while area subtrac-
An example of ρ is shown in figure 6. There are many tion is very effective for correcting the jet pT , it is not
subtleties in defining both ρ and A, which are discussed used to mitigate the pile-up dependence of jet substruc-
in e.g. Refs. (73; 74; 75). An extension to this method is ture observables as it is only able to correctly remove
shape subtraction (76), where randomly distributed ghost pile-up contributions on average. In fact, jet substruc-
particles are used to calculate a jet shape’s sensitivity to ture variables are among the most difficult to correct for
pile-up, which can then be corrected for non-uniformities pile-up because they are so reliant on radiation profiles.
in the spatial distribution of pile-up particles. A number of hybrid methods have been proposed op-
Instead of a global, collective, treatment of pile-up for erating at the event constituent level. One example is
10
arbitrary units
Pythia QCD PF+PUPPI
35000 <∆m>=-0.6 GeV
Anti-kT (R=0.8) RMS=10.9 GeV
<nPU> = 40
30000 PF
200 GeV < p < 600 GeV <∆m>=13.7 GeV
T
|η| < 2.5 RMS=17.9 GeV
25000 PF+CHS
<∆m>=-6.2 GeV
RMS=14.6 GeV
20000
PF(Cleansing)
<∆m>=-0.8 GeV
RMS=12.7 GeV
15000
PF+CHS([Link].)
<∆m>=0.5 GeV
10000 RMS=13.7 GeV
5000
0
-100 -80 -60 -40 -20 0 20 40 60 80 100
mreco - mgen(GeV)
V. JET SUBSTRUCTURE METHODS AND ATLAS performed a broad study of the relative per-
OBSERVABLES formance of different grooming techniques for boson-
tagging (86; 94; 95), top-tagging (96; 97) and SM
A. Jet Grooming measurements (98; 99), using the removal of pile-up-
dependence, the jet mass resolution, and the tagging effi-
Jet grooming techniques have seen a particularly high ciency versus background rejection as performance met-
level of interest from the experimental and theoretical rics. The ‘standard’ grooming procedure adopted by
communities alike. Jet grooming is an additional ‘post- ATLAS is trimming with fcut = 0.05 for boson tagging
processing’ treatment of large radius jets, an extra step in both Run 1 (Rsub = 0.3) and Run 2 (Rsub = 0.2).
used to remove unwanted soft radiation and to allow The trimming algorithm with the same parameters was
the underlying hard substructure associated with a two- adopted for top tagging, along with several other tech-
prong (e.g. W boson) or three-prong (e.g. top quark) niques (see section VI.C). Another technique currently in
decay to be identified more efficiently. use by ATLAS is the reclustering of small-R jets (100),
In particular, grooming is the systematic removal of which uses fully-calibrated anti-kT , R = 0.4 jets as inputs
radiation from within a jet, often targeting soft and wide to the anti-kT algorithm with a larger distance parameter
angle radiation. There are a variety of techniques and (typically R = 1.0). This has proven a popular method
each one has tunable parameters which are chosen to in ATLAS analyses due to the flexibility of optimizing
suite the particular needs of the application. The three the jet distance parameter depending on the considered
main algorithms used by ATLAS and CMS are trim- phase-space of the analysis (101; 102; 103). A recent
ming (90), pruning (91), and soft drop (92). In each study of in-situ measurements (104) (including ‘closeby’
of these cases, the constituents of a jet are re-clustered effects) confirm that the differences between data and
and soft/wide angle radiation is rejected in this process. simulation observed with reclustered jets are indeed cov-
For trimming, the kT algorithm is used to re-cluster and ered by simply propagating the uncertainties associated
the radius parameter of the re-clustering is called Rsub . with the input anti-kT , R = 0.4 jets.
Those smaller-radius jets with a momentum fraction CMS studied a large number of grooming techniques
f < fcut are removed to produce the trimmed jet. The in the context of boosted boson-tagging (57; 105), top-
two other algorithms impose a condition on each 2 → 1 tagging (106; 107) and SM measurements (108; 109).
clustering step, by going backwards in the sequence in During Run 1 the grooming techniques were used to-
which the particles were combined in the re-clustering. gether with charged-hadron subtraction for pile-up mit-
The transverse momentum fraction of the softer particle igation (see section IV). All groomers studied showed
to the merged system, z = min(pT,1 , pT,2 )/(pT,1 + pT,2 ), reasonable or good agreement between data and simu-
is a natural choice for determining the scale of the soft lation and the pruning algorithm (R = 0.8, zp = 0.1
radiation, and the angular distance ∆R between the two and dp = 0.5) showed the best performance for boson
particles for identifying wide-angle radiation. The differ- tagging (105). For Run 2, soft drop (zcut = 0.1 and β
ence between pruning and soft drop lies in the way how = 0) is used for jets with R = 0.8 in jet substructure
particles and their combinations get rejected based on analyses in CMS together with the pile-up removal al-
the values of z and ∆R. For pruning, the softer particle gorithm PUPPI (80) (see section IV). Soft drop jets in
of the 2 → 1 clustering step is discarded if z < zp and combination with PUPPI show a similar performance as
∆R < dp . For soft drop, the softer particle is discarded if pruning when comparing signal efficiency versus back-
z < zcut (∆R/R)β , where zcut and the angular exponent ground rejection (87; 107), but allow for better theoret-
β are free parameters6 . ical control. While grooming techniques were found to
The role of grooming has traditionally satisfied two improve the performance (higher background rejection
purposes in ATLAS, being the mitigation of pile-up ef- at fixed signal efficiency) of the jet mass, N -subjettiness
fects on jets, and the removal of soft/wide-angle radia- ratios (110; 111) were found to perform better without
tion. The particle flow algorithm employed in CMS in grooming for boosted boson tagging (105). For top-
conjunction with CHS or PUPPI allows for a correction tagging applications, however, soft drop groomed N -
for pile-up effects. This reduces the usefulness of groom- subjettiness ratios improved the performance with re-
ing for pile-up mitigation, but retains its advantage for spect to ungroomed ones for jets with pT < 400 GeV.
the removal of soft/wide-angle radiation. For higher pT jets there was no significant gain observed
with grooming for N -subjettiness ratios (107).
Events / 5 GeV
surement of jet mass requires detection of the deposited ATLAS Preliminary
4000
energy with a granularity that is finer than the size of a s = 8 TeV, 20.3 fb-1
jet. The mass of a jet can only be estimated if the en- p > 200 GeV, |η| < 2.0
T
the jet constituents. For jet substructure techniques that Fitted Simulation
Particle-level Simulation
rely on the rejection of soft particles, it is also important 2000
to be able to reconstruct particles with low pT separately
from harder particles in a jet.
1000
The jet mass response distribution Rreco is constructed
from the calibrated, reconstructed jet mass Mreco divided
by the particle-level jet mass Mtrue . The mass response 0
distribution is calculated in bins of reconstructed jet 1.1
Data/MC
Detector-level Simulation
Fitted Simulation
pT,reco and ηreco . In ATLAS, the Jet Mass Scale (JMS) is
1
defined as the mean of this response distribution. The Jet
Mass Resolution (JMR) is then defined as half the 68% 0.9
0 50 100 150 200
interquartile range (IQnR) of the response distribution, Jet Mass [GeV]
as
FIG. 9 The trimmed jet mass before (detector-level) and after
(fitted) determining s and r. The particle-level distribution
r = 0.5 × 68% IQnR(Rreco ). (2) is shown for comparison. Jets are required to have pT > 200
GeV. Adapted from Ref. (113).
This is robust to large non-Gaussian tails but, if the dis-
tribution is Gaussian, is equal to its 1σ width. The frac-
tional JMR is expressed as the JMR divided by the me-
dian of the response distribution. uncertainties are due to the theoretical modeling of jet
ATLAS has recently developed a data-driven approach fragmentation and the cluster energy scale.
to extract the JMS and JMR from an enriched sample As the forward-folding method is currently restricted
of boosted tt events, however the method can also be to jets with pT < 350 and 500 GeV for boosted W
extended to other final states. This forward-folding ap- bosons and top quarks, respectively, the results are com-
proach folds the particle-level mass spectra by a modified bined with the so-called Rtrk method which constrains
response function such that the JMS in a given bin of the mass scale by comparing the calorimeter jet mass
particle-level jet mass and reconstructed jet pT is scaled to the mass calculated from track jets and extends up
by the scale parameter s and the JMR is scaled by the to pT = 3000 GeV (114). The Rtrk method can also be
resolution parameter r: generalized to other variables and is used in ATLAS to
m,pT
constrain the pT scale of large-R jets as well as to derive
Mfold = s × Mreco + (Mreco − hMreco i)(r − s)). (3) systematic uncertainties on jet substructure variables.
The values of r and s for which the Mfold distribution best The concept of a Track-Assisted Mass for trimmed,
matches the data are extracted from a 2 dimensional χ2 large-R jets has been studied in ATLAS (112) to main-
fit as shown in figure 9 and detailed in Ref. (112; 113). tain performance for highly boosted particles due to the
With the forward-folding approach, the JMS and JMR limited granularity of the calorimeter. The track-assisted
for hadronically decaying boosted W bosons with pT & mass is defined as:
200 GeV are determined with 2–3% and 20% system-
atic uncertainties, respectively (see figure 10). As the pcalo
jet mass and its detector-response depend on kinematics mTA = T
× mtrack , (4)
ptrack
T
and jet substructure, the measurement was repeated dif-
ferentially with an increased luminosity for boosted W
and top quarks in Ref. (114). It will be important to where pcalo
T is the transverse momentum of the calorime-
extend the technique to other final states in the future. ter jet, ptrack
T is the transverse momentum of the four-
This may require hybrid data/simulation methods. A vector sum of tracks associated to the calorimeter jet,
detailed study of the various contributions to the JMS and mtrack is the invariant mass of this four-vector sum,
and JMR has been performed in context of the soft drop where the track mass is set to the pion mass mπ . The
mass measurement (98), described in section VII.A.1, by track-assisted mass exploits the excellent angular resolu-
propagating experimental uncertainties on the inputs to tion of the tracking detector and the ratio pcalo
T to ptrack
T
the jet reconstruction to the jet mass. The dominating corrects for charged-to-neutral fluctuations. The Com-
13
0.3
1.2 0.2
0.15
1
0.1 Calorimeter mass
Track assisted mass
Combined mass
0.8
0.05
500 1000 1500 2000 2500
8 TeV, 20.3 fb -1 Data 1σ stat. uncertainty Truth jet p [GeV]
T
0.6 13 TeV, 3.2 fb -1 Data stat ⊕ syst. uncertainty
FIG. 11 The ATLAS combined jet mass resolution. Adapted
0.8 0.9 1 1.1 1.2 from (112).
Relative Jet Mass Scale
Events / (5 GeV)
300
and resolution (r) for trimmed anti-kT , R = 1.0 calorimeter CMS CMS data tt (unmerged)
jets from the 2012 and 2015 ATLAS datasets√ and the 1σ sta- Preliminary Data fit Single top
tistical and total uncertainty ellipses. The s = 8 and 13 250 MC fit W+jets
TeV selections are similar, although the trimming definition
slightly changed between Runs (Rsub = 0.3 to Rsub = 0.2). tt (merged) WW/WZ/ZZ
Adapted from Ref. (112). 200
CMS Simulation 13 TeV one for the combined mass in ATLAS (figure 11), even
JMR 0.3
though quark/gluon jets are compared with W /Z-jets
0.25 p bins and very different technologies are used to reconstruct
T
260-350 GeV the jet mass.
0.2 650-760 GeV
0.15 900-1000 GeV
1200-1300 GeV C. Other Jet Substructure Observables
0.1
Events
Events / 0.02 1.2 ATLAS Preliminary
Data 2015+2016 2500 Data
-1
Pythia8 dijet (× 0.71) s = 13 TeV, 36.5 fb Merged top t t
1 Herwig++ dijet (× 1.32) Trimmed anti-k t R =1.0 Merged W t t
Dijet Selection
2000
W+jets (× 50) Merged QB tt
0.8 Z+jets (× 50)
pT > 450 GeV Not merged tt
mcomb > 100 GeV 1500 W+Jets
all-had tt (× 50)
Single t
0.6 Stat. uncert.
DY+Jets
Stat. ⊕ syst. uncert. 1000 QCD
0.4
AK8 PUPPI
500 p > 400 GeV
T
0.2
DATA / MC
1.5 1.5
MC(sig. + bkg.)
1.25
1
Data
1
0.75 0.5
0.50
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
Probe jet τ3/ τ2
Leading Large-R Jet τ WTA
32
Normalized Entries
ATLAS Simulation Preliminary
Events/0.02
400
CMS Z+jets, Pythia8 s = 13 TeV
350 Preliminary
Data quark
0.20 Anti k EM+JES R=0.4
| | < 2.1
300
gluon undefined Quark Jet
0.15 Gluon Jet
250 80 GeV < p < 100 GeV 50 < pT < 100 GeV
T
| η| < 2.0 400 < pT < 500 GeV
200 0.10 1200 < pT < 1500 GeV
150
0.05
100
50
0.000 10 20 30 40 50 60
0
2
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
ntrack
1.8
1.6
1.4 FIG. 17 The distribution of the number of tracks inside jets
DATA
MC
1.2
1
0.8 for quark and gluon jets in multiple jet pT ranges. Repro-
0.6
0.4
0.2
duced from Ref. (135).
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
pD
T
quark-jet efficiency
e.g. (57; 142; 143; 144; 145; 146; 147)). Additionally, it CMS
0.9 Simulation Preliminary
has been shown that an improved W tagger can be con-
dijets, Pythia8
structed by utilizing q/g discrimination on subjets (57). 0.8
Both ATLAS and CMS have developed likelihood-
0.7
based discriminants for explicit q/g tagging. The dis-
criminants are constructed from variables sensitive to the 0.6
radiation pattern of quark and gluon jets, also taking into 0.5
account differences between light (uds) and heavy flavor
(cb) quark jets, where the latter are more similar to gluon 0.4 0 < |η| < 1.3, 80 < p < 100 GeV
T
jets. ATLAS uses the number of tracks ntrk as an approx- 0.3 0 < |η| < 1.3, 40 < p < 50 GeV
T
imation for the number of jet constituents and the jet 2.7 < |η| < 3, 40 < p < 50 GeV
0.2 T
width (138) while CMS utilizes the number of particle-
flow constituents nconst , the jet axes and fragmentation 0.1
functions (87). Since the distributions of these variables
0
depend on η, pT , and ρ, the likelihood discriminators are 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
constructed differentially with respect to these variables. gluon-jet rejection
In Run 2, ATLAS also introduced a simple and robust
FIG. 18 The CMS q/g tagging performance in simulation for
tagger using solely ntrack (135), which has the advantage
two bins in jet pT and two bins in jet |η|. Reproduced from
of a much-simplified uncertainty derivation. Ref. (87).
Figure 18 shows the CMS q/g tagging performance in
simulation. The q/g label is obtained through a matching
of jets on the detector level to outgoing partons from the figure 19 is much more pronounced for MC than for data
matrix-element calculation. For a 50% gluon or quark indicates that the simulation over-predicts the difference
efficiency, the misidentification rate (quark or gluon) is between quark and gluon jets. In contrast, Herwig (not
about 10%. This performance depends slightly on the shown) tends to underestimate the performance observed
jet pT , in part because the particle multiplicity increases in data.
with pT (and therefore the performance improves). Out- Multiple samples with a different (but known) q/g
side the tracking acceptance (|η| & 2.5), q/g tagging sig- composition can be used to extract the distribution of
nificantly degrades due to the coarse calorimeter granu- q/g tagging observables. ATLAS and CMS have both
larity and increased pile-up sensitivity. used Z/γ+jets and dijet samples, which are enriched in
ATLAS (138) and CMS (141) are also actively study- quark and gluon jets, respectively. The extracted average
ing sophisticated approaches based on modern machine ntrack from data is shown using this method in figure 20.
learning. While these methods hold great promise for As expected, gluon jets have more particles on average
their power and flexibility, simple combinations of a small than quark jets and the multiplicity distribution increases
number of features often achieves a similar performance. with jet pT .
Machine learning architecture design and input optimiza- The Run 2 ATLAS tagger is based entirely on dijets,
tion are still an active area of research and development. exploiting the rapidity dependence of the q/g fraction to
extract the track multiplicity separately for quarks and
The modeling of q/g discriminating observables is a key gluons. A Run 1 measurement is used to constrain the
concern for tagging applications. Typically, Pythia (148; particle-level modeling, and dedicated track reconstruc-
149) tends to describe quarks better than Herwig (150; tion uncertainties are used to complement the particle-
151), whereas the opposite is observed for gluons. Pythia level uncertainty with a Run 2 detector-level uncertainty.
tends to overestimate the q/g tagging performance with The uncertainties on q/g tagging are 2-5% over a wide
respect to data, as illustrated quite strikingly in figure 19. range of 200 GeV . pT . 1 TeV at a working point of
This figure shows that gluon jets tend to have more tracks 60% quark jet efficiency (135). The template-based cal-
and have a broader radiation pattern relative to quark ibration can also be used to directly construct the q/g
jets9 . The fact that the hot spot in the bottom left of tagger in data; however, when more than two observables
are used to construct the tagger, it becomes impractical
to extract the high-dimensional templates.
9
The likelihood-based discriminant used for q/g tagging
The jet flavor is obtained as the type of the highest energy
parton from the event record inside the jet cone. This gives
in CMS in Run 2 is calibrated with a template-based
nearly the same result as the CMS definition discussed above for fit using two discriminant distributions obtained from a
the two leading jets in a 2 → 2 calculation, but also works well Z+jets and a dijet sample. The different quark and gluon
for additional jets in the event. fractions in each bin of the discriminant distributions
18
20
<ntrk>
1 ATLAS Preliminary
L = q/(q+g)
ntrk 18 18
s = 8 TeV 20.3 fb-1
16
0.9 |η| < 0.8
16 ATLAS 14
Discriminant for Data-Driven Tagger
∫ anti-k R=0.4, |η| < 0.8 0.8
-1
12
14 L dt = 4.7 fb , s = 7 TeV
10
t
160 GeV<p <210 GeV
0.7 8
12 T
6
0.6 4 Extracted Validation
10 2
0.5 Gluons Light Quarks
8 1.2 100 200 300 400 500
Validation
Extracted
0.4 1 Jet PT (GeV)
6 0.8
0.3 100 200 300 400 500
4 Jet PT (GeV)
0.2
2 FIG. 20 The average track multiplicity in ATLAS for
0.1
Z/γ+jets (quark-enriched) and dijets (gluon-enriched). The
0 0 dashed lines indicate the measurement on the validation
0 0.05 0.1 0.15
samples : Z/γ+2-jets (quark-enriched) and trijets (gluon-
Track Width enriched). Reproduced from Ref. (137).
1
L = q/(q+g)
ntrk
&T
QCD
GeV, the decay products of a boosted W and Z boson ATLAS s = 8 TeV
1 / ∈G
anti-k t R=1.0 jets Trimmed (f =5%,Rsub =0.2)
are typically only captured by a jet of radius R ∼ 1, while |η| < 1.2, 350 GeV < p < 500 GeV
cut
a number of analyses, most notably searches for diboson Monte Carlo Predicted
ATLAS Simulation Preliminary ysis). In order to achieve a flat signal efficiency, AT-
1 / ∈G
Efficiency
1/(W Efficiency)
1.4 CHS
65 GeV < MPruned < 105 GeV
PUPPI
65 GeV < MSoftdrop < 105 GeV
ATLAS Simulation
1.2 65 GeV < MPruned < 105 GeV + τ21 ≤ 0.45
CHS
103
1 65 GeV < MSoftdrop < 105 GeV + τDDT
PUPPI
21 ≤ 0.52 No discrimination
+
M (Jet Mass)
0.8 Q (Jet Charge)
B (b-tagging)
102 M+B
0.6
M+Q
Q+B
0.4 M+Q+B
W-jet, AK R = 0.8
10
0.2 p > 200 GeV
T
|η| ≤ 2.4 GeV
0
500 1000 1500 2000 1
Jet p (GeV)
T
CMS Simulation Preliminary 13 TeV 0 0.2 0.4 0.6 0.8 1
0.3
Mistag rate
65 GeV <
CHS
MPruned < 105 GeV Z Efficiency
PUPPI
65 GeV < MSoftdrop < 105 GeV
MPruned < 105 GeV + τ21 ≤ 0.45
CHS
65 GeV <
65 GeV < MSoftdrop < 105 GeV + τ21 ≤ 0.4
PUPPI FIG. 24 Background rejection versus efficiency for discrimi-
0.2 65 GeV < MSoftdrop < 105 GeV + τDDT
PUPPI
21
≤ 0.52 nating Z boson jets from W boson jets for various jet observ-
QCD, Pythia8 ables and their combinations. Reproduced from Ref. (173).
p > 200 GeV
T
|η| ≤ 2.4 GeV
W and Z bosons. This has been exploited in a number ble shower histories that can lead to the observed
of algorithms, which usually aim at an optimal perfor- leading final state anti-kT , R = 1.0 jet are calcu-
mance in a particular kinematic regime. Flavor tagging lated. Each shower history is assigned a probability
also plays a key role for top tagging, which offers its own weight factor based on the aforementioned consid-
challenges because the b jet from the b quark may not be erations (to be signal-like or background-like), then
isolated from the radiation resulting from the associated a likelihood ratio χ(pN ) is constructed, and the
W boson decay. Due to the heavier mass of the top quark log χ(pN ) is used as the discriminating substruc-
compared with the electroweak bosons, top tagging must ture variable. For top quark tagging, efficiencies
also operate in a moderate boost regime where the decay of 80% with misidentification rates of 50% for jets
products may not all be contained inside a single jet with with 500 < pT < 1000 GeV were observed. The
R . 1.0. efficiency increases with increasing jet pT , where a
The techniques for tagging boosted top quarks have plateau is reached for pT > 2000 GeV, with efficien-
evolved as fairly complex methods in comparison to the cies of ∼ 80% at 10% misidentification rate. Re-
V taggers; these techniques include: cently, the Shower Deconstruction algorithm was
optimized for top quarks with pT > 800 GeV in
(a) The Johns Hopkins / CMS top tagger (CMSTT) context of the W 0 to tb hadronic search (182) by
(15) was designed for tagging top quarks with pT using exclusive kT subjets.
> 1 TeV. The algorithm is based on a decomposi-
tion of the primary jet into up to four subjets by In addition to the dedicated techniques described
reversing the CA clustering sequence. It has been above, simpler algorithms using grooming and substruc-
adapted by the CMS Collaboration (175; 176), and ture similar to V tagging methods have been investigated
was adopted as the standard top-tagging algorithm by ATLAS. A performance study at 7 TeV (179) inves-
in CMS in Run 1, where it was typically used in tigated a variety of performance metrics relating to the
the region of pT > 400 GeV, with an average iden- usage of groomed jets. Different grooming algorithms
tification efficiency of 38% at 3% misidentification were investigated for their resilience to pile-up and mass
rate (106). resolution. It was concluded that trimmed anti-kT jets
(b) The HEPTopTagger (HTT) (177; 178) was de- with a distance parameter of 1.0 and trimming parame-
signed to target ttH production in the H → bb de- ters of Rsub = 0.3 and fcut = 0.05 were a good candidate
cay channel. In ttH production the top quark pT for a one-fits-all large-R jet definition. This jet definition
distribution peaks around 150 GeV and is steeply became standard in ATLAS for W /Z/H and top quark
falling towards increasing pT , where it is already tagging in Run 1. A later ATLAS study (96) investi-
an order of magnitude smaller at pT ∼ 400 GeV. gated the various methods available for tagging hadronic,
This results in a requirement of non-zero signal ef- highly boosted top√quarks. The so-called √ Tagger V has
ficiency already at pT ≈ 200 GeV, where the top Mjet >p100 GeV, d12 > 40 GeV and d23 > 20 GeV,
quark decay is only moderately boosted. The HTT where dij is the kT -splitting scale (12). The efficiency
achieves this with a large jet distance parameter of versus rejection is shown for various taggers in figure 25.
1.5 and a sequence of declustering, filtering and re- The difference between√Taggers III and V is the addi-
clustering of the original CA jet. The performance tional requirement on d23 in Tagger V. At √ efficiencies
of the HTT was studied by the ATLAS and CMS smaller than 45%, the W 0 tagger, based on d12 and the
Collaborations on data with a center-of-mass en- N -subjettiness ratios τ21 and τ32 , has better background
√ rejection than Taggers III and V. ATLAS also tested the
ergy s = 7 and 8 TeV (96; 106; 179). Efficiencies
of 10% with misidentification rates of 0.5% for jets HTT and Shower Deconstruction (183), which have been
with 200 < pT < 250 GeV were observed. The ef- found to have good background rejection (larger than 50)
ficiency increases with increasing jet pT , where a for efficiency values smaller than about 35%. However,
plateau is reached for pT > 400 GeV, with efficien- similar as for the CMS experiment, the background ef-
cies of approximately 40% at 3% misidentification ficiencies of the two taggers show a significant rise with
rate, very similar to the performance achieved with increasing pT .
the CMSTT. CMS has focused on enhancing the performance of
CMSTT and HTT by identifying observables which carry
(c) Shower Deconstruction (180; 181) was designed to discriminatory power, but have only small or moderate
be analogous to running a parton shower Monte correlations with the observables used in the main algo-
Carlo generator in reverse, where emission and rithm. Typically, correlation coefficients of about 0.3 or
decay probabilities at each vertex, color connec- less are required for noticeable improvement when aug-
tions, and kinematic requirements are considered. menting an algorithm with additional variables. Exam-
Small-radius (generally R = 0.2) subjets are re- ples for discriminating variables which fulfill this are N -
constructed with the CA algorithm and all possi- subjettiness ratios, energy correlation functions and their
23
103 CMS
Background rejection ATLAS
SD (C/A 0.2 subjets)
εB
350 GeV < pT < 400 GeV
HEPTopTagger (C/A 1.5)
800 < p < 1000 GeV, |η| < 1.5
d12 T
trimmed mass ∆R(top,parton) < 0.6
2
10
flat p and η
10−1
Tagger V (scan d23)
d23
τ32
Tagger I
10 Tagger II
Tagger III
10−2 CMSTT min. m
CMSTT top m
Tagger IV Filtered (r=0.2, n=3) m
HTT V2 fRec
Tagger V
HTT V2 m
W' Top Tagger Pruned (z=0.1, rcut=0.5) m
1
0 0.5 1 10−3 Q-jet volatility
Softdrop (z=0.1, β=0) m
Top quark tagging efficiency Softdrop (z=0.2, β=1) m
Trimmed (r=0.2, f=0.03) m
Ungroomed τ3/τ2
FIG. 25 Top quark tagging efficiency versus background re- log(χ) (R=0.2)
jection for various substructure variables and combinations in 10−40 0.2 0.4 0.6 0.8 1
ATLAS. Taken from Ref. (96).
εS
Tagging efficiency
ATLAS SD anti-kt R=1.0
tional b-tagged jet) to ensure that events contain a fully- 0.8
∫ L dt = 20.3 fb ,-1 C/A R=0.2 subjets
s = 8 TeV
merged top quark decay in a single large-R jet. This 0.7
correlated data-MC systematic uncertainty
can never be fully achieved, as no requirements on the 0.6 0.7 < |η| < 2.0
substructure of the large-R jet can be imposed without 0.5
Data/Sim.
ciency measurement by using simulated events, as done in 1.5
350 400 450 500
a study by the ATLAS Collaboration (96), with the draw- 1
0.5 Large jet p [GeV]
back of relying on a specific simulation and the ambigu- 350 400 450 T 500
ous definition of a fully-merged top quark decay. By not Large-R jet p [GeV]
T
correcting for non-merged top quark decays, efficiency
0.05
Mistag rate
values are obtained smaller than the ones suggested by 0.045
ATLAS
ROC curve studies, see for example (107). Instead of 0.04 ∫ L dt = 20.3 fb ,-1
s = 8 TeV
by light-flavor jets. The disadvantage of this approach is 0.005 Pythia Syst. uncertainty
Multi-jet rejection
and V H production, searches for boosted mono-H, or ATLAS Simulation Preliminary
vector-like quark searches in the tH and bH final states. 105 calo
p > 250 GeV, No mjet selection
T
Because of the large predicted branching fraction for the Double b-tag
Asymm. b-tag (70% wp)
H → bb̄ decay of about 58%, its coupling to b quarks is 104
Single b-tag
one of the most interesting to study. For a large frac- Leading subjet b-tag
103
tion of Higgs bosons with pT > 300 GeV, the two b quark
jets merge into a single jet for a jet distance parameter of 102
R = 0.8 or 1.0, as used in CMS and ATLAS, respectively.
Several phenomenological studies have explored H → bb̄ 10
tagging algorithms using jet substructure, though ulti-
mately the optimal performance comes from using a com- 1
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
bination of substructure information and the track and Higgs-jet efficiency
vertex information related to the B hadron lifetime.
The approaches to identify boosted H → bb̄ candidates FIG. 30 The rejection of inclusive multijets versus Higgs-jet
that have been explored (and used) at CMS and ATLAS efficiency using all large-R jets with pT > 250 GeV for sin-
include: gle, double, asymmetric, and leading subjet b tagging require-
ments. Taken from Ref. (195).
(a) Subjet b tagging (174; 191; 192; 193; 194; 195),
where ‘standard’ b tagging is applied to each of
the subjets (the standard for CMS is the CSVv2 case of asymmetric b tagging, also due to the 70% b tag-
algorithm (196), and for ATLAS is MV2c20 (190)). ging working point requirement on one of the track jets.
Tagging b-jets in dense environments is of particu- The CMS double-b tagging algorithm (189; 198) at-
lar importance here, and was studied by ATLAS in tempts to fully exploit the strong correlations between
Ref. (197). In CMS subjets with R = 0.4 are clus- the b hadron flight directions and the energy flows of
tered with the kT algorithm using the constituents the two subjets, while adapting the variables used in the
of the large-R jet, while for ATLAS track jets with CSVv2 algorithm. The flexibility of the double-b tagger
a radius of 0.2 are matched to the large-R jet us- is ensured by avoiding a strong performance dependence
ing the ghost-association technique. At high pT on the jet pT and mass.
the subjets start to overlap causing the standard b With the double-b tagger, at the same signal efficiency,
tagging techniques to break down due to double- the misidentification rate is uniformly lower by about a
counting of tracks and secondary vertices when factor of two compared to the subjet b tagging approach.
computing the subjet b tag discriminants. Given the different kinematic properties expected for a
bb̄ pair originating from the decay of a massive resonance
(b) Double-b tagging (189; 195; 198), where in AT- compared to gluon splitting, the misidentification rate for
LAS, the term double-b tagging means that the the gluon splitting background reduces from 60% to 50%
two leading pT track jets must pass the same b at 80% signal efficiency and from 20% to 10% at 35%
tagging requirement. In CMS, the double-b tag- signal efficiency. At high pT , even larger performance
ger (189; 198) uses the N -subjettiness axes and the improvements are observed, which is an important gain
pruned anti-kT , R = 0.8 jet mass with a window for searches for heavy resonances, where very high pT
of 50 < M < 200 GeV to reduce the multijet back- jets are expected. In figure 31 the signal efficiencies and
ground. misidentification rates for the double-b tagger are shown
as a function of jet pT for three operating points: loose,
The Higgs-jet efficiency versus the inclusive multijet
medium and tight, which correspond to 80%, 70% and
rejection are shown in figure 30 for ATLAS subjet b tag-
35% signal efficiency, respectively, for a jet pT of about
ging, where the performance curves are shown for double-
1000 GeV. The misidentification rate is mostly flat across
b tagging, leading subjet b tagging, and asymmetric b tag-
the pT range considered while the signal efficiency de-
ging11 requirements. None of the curves reach a Higgs-jet
creases with increasing pT , as expected from the degra-
efficiency of 100% due to the imperfect efficiency to re-
dation of the tracking performance inside high pT jets.
construct the track jets needed for b tagging and, in the
Due to the small cross section of producing events with
boosted H → bb̄ or Z → bb̄ jets, the efficiency of the
11 Asymmetric b tagging means that among the two leading pT
ATLAS and CMS Higgs identification algorithms is mea-
track jets, the track jet with the largest b tagging weight must sured using QCD multijet events enriched in jets from
pass the fixed 70% b tagging working point threshold, while the gluon splitting, g → bb̄ with a topology similar to that of
b tagging requirement of the other jet is varied. boosted H → bb̄ jets.
27
CMS Simulation Preliminary (13 TeV) One of the major backgrounds for analyses selecting
1
boosted H or Z bosons decaying to bb̄ is tt̄ production.
Mistagging efficiency
double-b-tag
Subjet CSVv2
The misidentification rate for boosted top quark jets fak-
Fatjet CSVv2
ing H jets was measured in data by CMS (189; 198) in
AK8
70 < m < 200 GeV , 300 < p < 500 GeV enriched data samples of lepton+jets tt̄ events.
T
10−1 As previously discussed, for high pT of the Higgs bo-
son, the two subjets from b quarks start overlapping and
the performance of identifying the subjets as fixed-radius
track jets decreases significantly. To improve the perfor-
10−2
mance of the ATLAS standard H → bb̄ identification
algorithm for searches that require the presence of high
pT Higgs bosons, the ATLAS Collaboration studied alter-
native methods like the use of variable-radius track jets,
exclusive kT subjets, calorimeter subjets reconstructed in
10−3 the center-of-mass frame of the Higgs jet candidate (199)
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Tagging efficiency (H→bb) and the combination of three jet shape and jet substruc-
ture variables into a multivariate discriminator (200).
FIG. 31 The misidentification rate for inclusive multijets ver- For highly boosted Higgs bosons, these reconstruction
sus Higgs-jet efficiency using jets with 300 < pT < 500 GeV techniques significantly outperform the usage of fixed-
and pruned jet mass 70 < m < 200 GeV for three different b radius track jets.
tagging requirements. Taken from Ref. (198).
1.6
1.4 calibration factors derived from simulation and using in
1.2
1 situ techniques (from boosted W bosons) to correct the
0.8
0.6 jet energy and mass scale. Furthermore the jet energy
0.4 and mass are smeared in MC simulation to match the
0.2
0 20 40 60 80 100 120 140 160 180 200
Jet Mass [GeV]
resolution measurements in data. Various sources of sys-
tematic uncertainties, categorized as experimental and
FIG. 32 Normalized differential cross section as a function theoretical uncertainties, that impact the jet mass mea-
of the jet mass for CA jets with R = 1.2 after splitting and surement are taken into account. While CMS evaluated
filtering, taken from Ref. (201). the effect of the jet energy and mass scale uncertainties
on the measurement by varying the energy and mass by
their respective uncertainties, ATLAS evaluated the ex-
formed by the ATLAS Collaboration with a dataset cor- perimental uncertainties based on the accuracy of the
responding to 35 pb−1 of 7 TeV pp collisions (201). Both modelling of the topological cluster energies and posi-
the cross section for groomed and ungroomed CA R = 1.2 tions as well their reconstruction efficiency. Theoretical
jets was measured separately to gain sensitivity to both uncertainties on the physics model are taken into account
the hard and soft jet physics and to gain a deeper under- by comparing the response matrix for various MC gener-
standing of the various effects involved in QCD radiation. ators.
For the ungroomed jet mass, large discrepancies were ob- The comparison of the normalized cross section with
served in the tails of the mass distribution between the two analytical calculations as measured by CMS is shown
predictions from the MC event generators Pythia and in figure 33. ATLAS measured instead the log10 ρ2 distri-
Herwig++, and the data, whereas the core of the mass bution, shown in figure 34, where ρ is the ratio of the soft
distribution agreed within approximately 20% over the drop jet mass to the ungroomed jet pT . Both measure-
considered pT range. The largest discrepancies occur at ments are compared to calculations at next-to-leading
low jet masses which is sensitive to the underlying event order with next-to-leading-logarithm and leading order
description, hadronization model and pile-up effects. The with next-to-next-to-leading-logarithm accuracy. Good
normalized cross section after applying the split filtering agreement between the data and the predictions is ob-
algorithm (14) is shown in figure 32 with the mass drop served in resummation regime −3.7 < log10 ρ2 < −1.7.
parameters µfrac = 0.67 and yfilt = 0.09, and a filter- For higher jet masses, where fixed-order effects play an
ing parameter of Rfilt = min(0.3, ∆R/2). After removing important role, the NLO+NLL calculation provides a
soft radiation from the jet which is difficult to model, better description than the LO+NNLL calculation.
the MC prediction is in excellent agreement with the In addition to generic QCD jets, the jet mass has also
data within statistical precision. The CMS Collaboration been measured for boosted top quarks in lepton+jets
performed a similar measurement with anti-kT R = 0.7 tt̄ events collected by the CMS Collaboration at 8
jets using various grooming
√ techniques in selected dijet TeV (116). This measurement is the first jet mass distri-
events using 5 fb−1 of s = 7 TeV data and found as bution unfolded at the particle level probing three prong
well that the agreement between data and the MC predic- decays. Large-R jets are reconstructed with the CA al-
tion improves significantly after grooming techniques are gorithm using a distance parameter of 1.2. The larger
applied (108). Furthermore a measurement of the cross value of R in this measurement compared to the default
section was performed in V +jet final states which over- R = 0.8 applied for top tagging applications in CMS
all show a slightly better data/MC agreement than that is due to an optimization of of statistical precision ver-
observed in dijet events suggesting that the simulation of sus the width of the jet mass distribution at the parti-
quark jets is better than for gluon jets. cle level and the JMR. The number of fully-merged top
The CMS (ATLAS) Collaboration measured the quarks grows with increasing R, but so does the width of
double-differential jet cross section in balanced dijet the jet mass distribution and the suceptibility to pile-up
29
d2σ (1/GeV) CMS 2.3 fb-1 (13 TeV) 19.7 fb-1 (8 TeV)
0.016
dmjet GeV
650 < p < 760 GeV 2.5
fb
T
0.014 Data CMS Data
Stat. + syst. unc.
0.012
T
Stat. unc.
dσ/dp dmg dp
MADGRAPH+PYTHIA
dσ
0.01 PYTHIA8
HERWIG++
2
0.008
MC@NLO+HERWIG
POWHEG+PYTHIA8
T
Frye et al
1
0.006 POWHEG+PYTHIA
Marzani et al
1.5
0.004
0.002
0 1
2
Theory
Data
1
0
20 30 40 100 200 1000
0.5
Groomed jet mass mg (GeV)
0
FIG. 33 Normalized differential cross section as a function of 150 200 250 300 350
the mass for jets groomed with the soft drop algorithm in data Leading-jet m [GeV]
jet
and for two theoretical calculations. Taken from Ref. (109).
Data
0.6 s= 13 TeV, 32.9 fb-1
ATLAS Pythia 8.1
Pythia 8.1
lead Sherpa 2.1
2.1 Taken from Ref. (116).
T
0.50.6 s= drop,
Soft 13 TeV,
β = 1,32.9 -1
zcut =fb0.1
Herwig++ 2.7
Herwig++ 2.7
T
soft drop
/p
0.4 LO+NNLL
[(m
T NLO+NLL+NP
0.5 NLO+NLL+NP
19.7 fb-1 (8 TeV)
Ratio to Data (1 / σresum) d σ / d log
cut
GeV-1
0.4
10
The normalized mass distribution from boosted top uncertainties are the track pT resolution and the choice
quarks, shown in figure 36, can be used to extract the of MC generator used to construct the response matrix
top quark mass. The normalized distribution is used (Pythia 6 versus Herwig++) for the CMS Collaboration
since only the shape can be reliably calculated, and it has whereas the uncertainties on the unfolding procedure, the
the additional benefit that systematic uncertainties par- jet energy resolution at low pT and uncertainties on the
tially cancel. The top quark mass is measured to be mt = tracking at high pT dominate the measurement of the
170.8±6.0 (stat)±2.8 (sys) ±4.6 (model)±4.0 (theo) GeV ATLAS Collaboration. The unfolded jet charge distri-
in agreement with top quark mass measurement in re- bution (κ = 0.6) of the leading jet in data is compared
solved tt̄ events (see e.g. Refs. (202; 203; 204; 205)), al- to the prediction from Powheg+Pythia8 (PH+P8) and
beit with a much larger uncertainty. This constitutes Powheg+Herwig++ (PH+HPP) in figure 38. The dif-
a proof-of-principle, presenting the possibility to extract ferent hadronization and fragmentation model used by
a fundamental SM parameter from a jet mass distribu- Pythia8 and Herwig++ have the largest impact on the
tion. This is of particular interest, as ambiguities arise in jet charge distribution. Variations of the jet charge can
the interpretation of traditional mt measurements (206), also be observed for different PDF sets however the ef-
which can be circumvented by measurements and analyt- fect of the relative flavor fraction in the dijet samples
ical calculations in the highly-boosted
√ regime (207; 208). is significantly smaller than the choice of the showering
Future measurements at s = 13 TeV will allow for a and fragmentation model. It was further found that the
higher statistical precision and, in combination with jet predicted jet charge distribution has a significant depen-
grooming and pile-up mitigation techniques, lead to a dence on the chosen value of αs that describes final state
large improvement in the total precision of the measure- radiation whereas it is insensitive to NLO QCD effect in
ment. Measurements at higher jet pT will facilitate com- the matrix element calculation, color-reconnection and
parisons with analytical calculations. multiple parton interactions. These findings are consis-
tent between the ATLAS and CMS Collaboration.
In addition to studying the sensitivity to various non-
2. Jet Charge perturbative aspects of hadronization and parton distri-
bution functions, the jet charge measurement by ATLAS
The jet charge (209; 210) is defined as the energy includes the first direct comparison of a jet substructure
weighted sum of the electric charges of the jet con- quantity with a perturbative calculation at the LHC. As
stituents it is not collinear safe, the average jet charge is not cal-
X pT,i κ culable. However, the pT dependence for a particular jet
Qκ = qi , (9) type has been calculated (212; 213). A new technique
pT,J
i∈J
was introduced in Ref. (209) to separately extract the
where qi is the electric charge of particle i and the free average up and down quark jet charge. For a fixed pT ,
parameter κ that controls the sensitivity to soft particles the more forward of the two dijets has a higher energy
within the jet. The ATLAS (CMS) Collaboration mea- and is therefore more likely to be the scattering parton
sured the jet charge for different values of κ using anti-kT with a higher momentum fraction of the proton. In turn,
jets with a radius parameter of R = 0.4 (R = 0.5) in a the higher momentum fraction parton is most likely to be
sample of dijet events. The ATLAS Collaboration distin- a valence quark. Therefore, the fraction of up quark jets
guishes between the two leading jets using the pseudora- is higher for the more forward dijet than the more central
pidity instead of the pT to avoid cases where the lead- dijet. Assuming further that the jet charge is entirely de-
ing particle-level jet is reconstructed as the sub-leading termined by the jet pT and parton origin, one can then
detector-level jet due to the jet energy resolution and solve a system of equations to extract the average up and
to gain sensitivity to different jet flavors. The average down quark jet charge in each bin of jet pT :
jet charge at detector- and particle-level for the more
forward of the leading jets and for κ = 0.5 is shown hQfJ i = fuf hQuJ i + fdf hQdJ i
in figure 37. Due to the increasing fraction of scatter- hQcJ i = fuc hQuJ i + fdc hQdJ i, (10)
ing valence up quark jets (up quark charge > 0), the
average jet charge increases with pT . The difference of where f = forward, c = central, u = up and d = down.
the average jet charge distribution at detector-level and As expected (though not an input), the average up quark
particle-level in figure 37 shows that the unfolding cor- charge is positive and the average down quark charge
rections are large and growing at high pT , due to the loss is negative; furthermore, the latter is roughly half the
of charged-particle tracks inside jets as a result of track former in absolute value. The pT dependence of hQu,d J i
merging. The average jet charge as predicted by Pythia are fit with a logarithmic scale violating term c: hQii =
8 (148) using the Perugia tunes (211) is smaller than that hQi0 (1 + cκ ln(pT,i /pT,0 )), where i represents the pT bin.
observed in data due to a well-known over-estimation of Figure 39 shows the measured and predicted values of
the multiplicity inside jets. The dominating systematic cκ . The uncertainties are large, but there is an indication
31
0.15
0.05 0
Detector-level Data
Detector-level Pythia
Particle-level Pythia
-0.05
0
1.2
Detector/Particle Pythia/Data
1
-0.1 2012 Data
0.8
1.5 QCD (Leading Power)
1
-0.15
0.5
0 500 1000 1500 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Jet p [GeV] κ
T
FIG. 37 The detector- and particle-level average jet charge FIG. 39 The measured and predicted value of the average
as a function of jet pT . Reproduced from Ref. (209). jet charge scale violation parameter cκ . Reproduced from
Ref. (209).
-1
CMS 19.7 fb (8 TeV)
1.4
p > 400 GeV
angularity (201; 225).
1.2 T
|η| < 1.5
[1/e]
Data
1 PH + P8 (CT10)
0.6
PH + HPP (CT10)
0.8 B. Measurements with Jet Substructure
1/N dN/dQ
PH + P8 (HERAPDF)
0.6
0.4 While measurements of jet substructure observables
0.2 such as jet mass, jet charge and event shape variables
0
−1.5 −1 −0.5 0 0.5 1 1.5
have been discussed in section V, the following sections
Leading-jet Qκ =0.6 [e] present measurements of other quantities through the ex-
1.1 ploitation of jet substructure techniques such as top tag-
Data
MC
1 ging.
0.9
−1.5 −1 −0.5 0 0.5 1 1.5
1.5
0.5
−1.5 −1 −0.5 0 0.5 1 1.5
The selection cuts applied in traditional tt cross sec-
FIG. 38 Unfolded jet charge distribution for κ = 0.6 in data tion measurements (226; 227; 228; 229; 230; 231; 232)
and MC prediction. Taken from Ref. (210) are chosen to maximize the acceptance and minimize the
. associated uncertainties on the fiducial and total cross
section measurements. The fiducial region is such that
events with top pT below 100 GeV and above 600 GeV
that c < 0 and ∂c/∂κ < 0, as predicted. are under-represented, with the former caused by trig-
ger and reconstruction efficiencies and the latter by colli-
mated decays from large Lorentz boosts. This is evident
3. Other Jet Substructure Observables from figure 40, where a drop in selection efficiency below
100 GeV and above 600 GeV is apparent. This results in
The ATLAS and CMS Collaborations have performed a small number of events being selected with high top
further precision measurements of hadronic jet substruc- quark pT , as seen in the ATLAS Run 1 (7 TeV) mea-
ture in pp collisions, correcting for acceptance and resolu- surement shown in figure 41. This means that a very
tion such as jet and event shapes (214; 215; 216; 217; 218; interesting region in terms of new physics is the least well-
219), charged particle multiplicities (217; 220; 221), the measured. Despite often having similar signal efficiencies
jet fragmentation functions (222; 223), color flow (224) to resolved reconstructed techniques, boosted top tagging
and kT splitting scales, N -subjettiness ratios as well as techniques allow for more precise measurements at high
further substructure variables such as Planar Flow and pT due to their higher background rejection.
32
10
dσtt/dp t [fb/GeV]
GeV
Fiducial phase-space
s = 7 TeV
1
104 Parton-level hadronic top-quark Data
dEvents Particle-level th POWHEG+PYTHIA CT10+hdamp=∞
POWHEG+PYTHIA CT10+hdamp=mtop
T
T
103
dp
ATLAS Simulation
1 POWHEG+PYTHIA HERAPDF+hdamp=∞
POWHEG+PYTHIA POWHEG+PYTHIA HERAPDF+hdamp=mtop
2
10
10-1
10
ATLAS
s = 8 TeV, 20.3 fb-1
1
-2
10
10-1
2
Pred. / Data
300 400 500 600 700 800 900 1000 1100 1200
10-2 1.5
1.5 1
Normalised ratio
0.5
300 400 500 600 700 800 900 1000 1100 1200
1
Particle top-jet candidate p [GeV]
T
0.5
0 100 200 300 400 500 600 700 800 FIG. 42 Particle-level differential tt̄ cross section measure-
Top p [GeV]
T ment for two different PDF sets and choices of the hdamp
parameters. Taken from Ref. (233).
Data
stat. ⊕ syst. uncert.
in figure 42 to the predictions of several MC generators
POWHEG+PYTHIA normalized to the NNLO+NNLL inclusive cross section.
dp (th)
dσfid
POWHEG+HERWIG
10-2 Overall good agreement is observed, but a harder pT
T
POWHEG(HERAPDF)+PYTHIA
MC@NLO+HERWIG spectrum is predicted by the simulation than observed
in data with larger discrepancies at high pT . The dif-
ferential cross section measurement is also compared to
10-3 predictions from Powheg+Pythia using either the HER-
APDF (234) or CT10 (235) PDF set and two differ-
ent values for the resummation damping factor hdamp ,
hdamp = mtop and hdamp = ∞. The best data/MC
10-4 agreement is observed when using the HERAPDF set
1.80 200 400 600 800
and hdamp = mtop . For each of the settings, the trend
Expected/Data
T
t,2
pt,1 > 500 GeV, p > 350 GeV MG5_aMC@NLO+Py8
MC@NLO+Herwig6 T T
Sherpa 2.2.1
T
tt
Stat. ⊕ Syst. Unc.
10−3
10-1
10−4
Theory / Data
1.5
1
10−5
Prediction
0.5 1.5
400 500 600 700 800 900 1000 1100 1200
Data
1
Particle-level t jet p (GeV)
T 0.5
0 100 200 300 400 500 600 700 800
FIG. 43 Comparison of the particle-level differential tt̄ cross ptt [GeV]
T
section as a function of the jet pT to three different MC gen-
erators. Taken from Ref. (236). FIG. 44 The normalized differential cross section as a func-
tion of the tt̄ pT as measured by ATLAS in the all-hadronic
channel at 13 TeV. Taken from Ref. (238).
amount of background from QCD dijet production. The
cross section is determined from a maximum likelihood fit
to the jet mass distributions for signal-enriched and sig- cross section as a function of the top jet pT and rapidity
nal depleted regions. This allows for a simultaneous ex- is in good agreement with the different MC predictions.
traction of the tt̄ cross section and the QCD background. Larger discrepancies are observed for the pT of the tt̄
The measurement is in agreement with the results from system as shown in figure 44. The measurement is dom-
the lepton+jets final states, but has somewhat larger sta- inated by the systematic uncertainties on the jet energy,
tistical uncertainties of up to about 40% in the highest mass and substructure scale of the large-R jets, alterna-
pT bin with 0.8 < pT < 1.2 TeV. tive parton shower model and the uncertainties on the b
√
The increased s at Run 2 of the LHC offers the pos- jet identification.
sibility for more precise differential tt̄ cross section mea-
surements in the highly-boosted regime. The tt̄ produc-
tion cross section increased by more than a factor√ of ten 2. W/Z/H Cross Sections
for top quark pT > 400 GeV when going from s = 8 TeV
to 13 TeV. The cross section of boosted W and Z boson pro-
A first measurement based on 3.2 fb−1 of 13 TeV data duction was measured by ATLAS in 4.6 fb−1 of 7 TeV
in the lepton+jets channel has been performed by AT- pp collisions (99). The hadronically decaying W and Z
LAS (232). The measurement extends to pT of 1.5 TeV bosons are reconstructed as one single ungroomed anti-
and a similar trend as at 8 TeV is observed between the kT R = 0.6 jet with pT > 320 GeV, |η| < 1.9 and masses
data and the simulation at high pT . A newer mea- ranging between 50 and 140 GeV. The W and Z sig-
surement of the tt̄ differential cross section in the all- nal is enhanced over the dominating QCD background
hadronic channel is performed by the ATLAS Collabora- by constructing a likelihood discriminant from three sub-
tion with 36.1 fb−1 of 13 TeV data (238). The measure- structure variables; thrust minor (19; 20), sphericity (18)
ment uses trimmed anti-kT R = 1.0 jets with Rsub = 0.2 and aplanarity (18), resulting in a signal efficiency of 56%
and fcut = 0.05. To obtain a flat signal efficiency and a background rejection of 89%. The jet mass distri-
of 50% and a quark/gluon rejection of approximately bution after subtracting the expected background from
17 (10) for pT = 500 (1000) GeV, pT dependent cri- tt̄ events is shown in figure 45. A binned maximum likeli-
teria are applied on the jet mass and τ32 . Further- hood fit to the jet mass distribution is used to extract the
more the two top-tagged large-R jets are required to W/Z jet signal yield and to calculate the inclusive cross
have a b tagged small-R jet within ∆R < 1.0. The section. Only the combined W + Z cross section mea-
event selection results in a signal-to-background ratio of surement is performed in this analysis due to the limited
approximately 3:1. The measured fiducial phase-space jet mass resolution. The combined W +Z cross section is
cross section is σ = 292 ± 7 (stat) ± 76 (sys) fb com- measured to be σW +Z = 8.5 ± 0.8 (stat.) ±1.5 (syst.) pb
pared to the Powheg+Pythia8 prediction of 384 ± 36 fb and is in agreement with the Standard Model prediction
at NNLO+NNLL. The measured normalized differential of σW +Z = 5.1 ± 0.5 pb within 2 standard deviations.
34
Events / 7 GeV
Data 2011 8000
ATLAS 3000
Data − multijet − tt
10
5
σData
0
50 100 150 200 250 0
−5
40 60 80 100 120 140 160 180 200
Jet Mass [GeV] mSD (GeV)
FIG. 45 Binned maximum likelihood fit to the jet mass distri- FIG. 46 Soft drop jet mass mSD of anti-kT R = 0.8 jets
bution in data for selected W/Z events reconstructed as one in data and for the dominating background processes; mul-
single ungroomed anti-kT R = 0.6 jet. Taken from ref. (99). tijet production and W/Z+jets events. Jets are required to
pass criteria on N21 and to be identified as double-b jets by
the double-b tagger introduced in section VI.D. Taken from
The dominating systematic uncertainties are the jet mass Ref. (241).
resolution and the choice of the QCD background PDF.
The signal significance was furthermore studied when us-
ing groomed jets instead of ungroomed jets. Without an by utilizing the double-b tagger. The W/Z+jets back-
optimization of the analysis for groomed jets, similar sig- ground is estimated from MC simulation and the shape
nificances were observed for groomed and ungroomed jets of the multijet background is determined in a validation
as expected due to the low number of pile-up vertices in region in data with lower values of the double-b tagger
the 7 TeV dataset. discriminator. The soft drop mass distribution of the
As discussed in section VI.D the SM Higgs boson leading jet is shown in figure 46 with a clear resonant
decays with approximately 58% into bb̄. However the structure at the mass of the W and Z boson. The SM
H → bb̄ decay in the resolved channel can only be stud- background processes and the potential signal from SM
ied in associated production with either a vector boson H → bb̄ production are estimated simultaneously. The
(W/Z) (239; 240), top quarks, or via the vector-boson- observed (expected) significance for the H → bb̄ process
fusion production mechanism due to the overwhelming is 1.5(0.7)σ. The measured cross section for the Z+jets
multijet background. To search for H → bb̄ in the gluon- process is 0.85 ± 0.16 (stat.) +1.0
−0.4 (syst.) pb which is in
gluon fusion production mode with an additional high- agreement with the SM prediction of 1.09±0.11 pb. This
pT jet, jet substructure techniques can be employed to is the first observation of Z → bb̄ in the single jet topol-
suppress the enormous multijet background. The CMS ogy.
Collaboration performed a search for the SM Higgs bo- The ATLAS Collaboration also measured the high pT
son using a dijet topology with 35.9 fb−1 of 13 TeV pp Z → bb̄ cross section using two nearby b tagged anti-kT
collisions (241). The analysis uses anti-kT R = 0.8 jets R = 0.4 jets (instead of one large-radius jet) in 19.5 fb−1
corrected with the PUPPI algorithm to reduce the effects of 8 TeV pp collisions (242). The measured fiducial cross
from pile-up, and modified with the soft drop algorithm section was determined to be σZ→bb̄ = 2.02 ± 0.33 pb
(β = 0, zcut = 0.1) to mitigate the effects from the under- which is in excellent agreement with the next-to-leading-
lying event and soft/wide-angle radiation. At least one order theoretical predictions.
large-R jet with pT > 450 GeV and |η| < 2.5 is required.
To distinguish the two prong structure of a jet containing
the full H → bb̄ decay from quark- or gluon-initiated jets, VIII. SEARCHES FOR NEW PHYSICS
the N21 variable, calculated from the generalized energy
correlation functions, is exploited. To ensure a flat QCD Jet substructure methods have been successfully ap-
background rejection of 26% over the considered mass plied in a large variety of searches for physics beyond
and pT range, a decorrelation procedure (167) is applied the SM. The respective exclusion limits are substantially
to N21 . The multijet background is further suppressed improved through the application of these methods. In
35
Data-Fit
[Link]
288; 289; 290; 291; 292; 293; 294; 295; 296; 297; 298; 299; Dijet invariant mass (GeV)
0
300; 301; 302; 303; 304; 305; 306; 307).
−2
1500 2000 2500 3000
Dijet invariant mass (GeV)
A. Diboson Resonances
FIG. 47 Dijet invariant mass distribution in the high purity
Several new physics models predict resonances cou- W Z category of the fully hadronic W W /W Z/ZZ resonance
pling strongly to vector bosons to play a role in the can- search. The fit under the background-only hypothesis is over-
cellation of large corrections to the Higgs mass. These layed. Taken from Ref. (115).
models include extensions of the SM Higgs doublet,
where the simplest realizations are two-Higgs-doublet
models (308) with heavy, neutral Higgs bosons, which have one HP V jet, and are divided into HP and LP
can have large branching fractions to top quarks and events, depending on whether the other V jet is of high
W /Z/H bosons. Alternatives are composite Higgs mod- or low purity. To further suppress the large QCD mul-
els (309; 310; 311; 312; 313; 314; 315; 316) or Randall- tijet background a requirement on the dijet kinematics
Sundrum Kaluza-Klein models (317; 318; 319; 320). |η1 − η2 | < 1.3 is applied.
Searches for new resonances generally focus at high The background is estimated from a sig-
masses with m > 1 TeV such that the SM bosons re- nal+background fit with the function dm dN
jj
= (mjj /P√
0
s)P1
,
ceive high Lorentz boosts. In more than 60% of the where P0 is a normalization parameter and P1 is a pa-
cases, W /Z/H bosons decay into a quark anti-quark rameter describing the shape. This parametrization
pair, which makes the reconstruction of such decays with has been tested and validated on simulated events and
jet substructure techniques an essential ingredient for on data in a control region. As shown in figure 47 the
these searches. In the following, the analysis strategies data in the signal region is well described by the fit
and results from CMS and ATLAS using pp collision data function. Figure 47 also shows that no excess over the
√
with s = 13 TeV are discussed. background-only hypothesis is observed.
The searches for diboson resonances are performed A similar analysis has been performed by the ATLAS
in semi-leptonic (257; 262) and fully hadronic final Collaboration (165). In this analysis events are required
states (115; 163; 164; 165). As the methods of jet sub- to have at least two large-R jets with pT > 200 GeV in
structure analyses exhibit their full strength in hadronic the pseudo-rapidity range |η| < 2.0. These jets are re-
final states, the following discussion gives a summary and constructed with the anti-kT algorithm with a radius pa-
comparison of the ATLAS and CMS results in the search rameter R = 1.0. The trimming algorithm is applied
for W /Z resonances in hadronic final states only. using kT subjets with R = 0.2. The rapidity separation
In an analysis performed by the CMS Collabora- between the two leading jets has to satisfy |∆y12 | < 1.2.
tion (115) events with two anti-kT jets with R = 0.8, cor- The large-R jet mass is computed from the Combined
rected with the PUPPI algorithm, and 65 < msoft drop < Mass (see section V.B), and is required to be within a
105 GeV are selected. The jet is considered to be a W bo- window of the expected W or Z mass value. The window
son candidate if the mass is in the range 65–85 GeV, while width varies from 22 to 40 GeV depending on the jet pT .
it is a Z boson candidate if the mass is in the range 85– In addition, the D2β=1 variable is used to select jets with
105 GeV. This leads to the three signal categories W W , a two-prong structure.
ZZ and W Z. The jets are further categorized according Similar as in the CMS analysis, the background is es-
to τ21 into high purity (HP, τ21 < 0.35) and low purity timated by fitting the dijet invariant mass distribution
(LP, 0.35 < τ21 < 0.75). Events are always required to with the parametric form dn dx = p1 (1−x)
p2+ξp3 p3
x , where
36
Events
CMS Data CMS
Other 106 QCD (MC)
2 Simulation
10 tt tt
Z’ 2.0 TeV, 1% width 105 TT→tHtH (500 GeV/c2)
102
1
10
Data / Bkg
1.5
1
1
10-1
0.5
0 1000 2000 3000 0 1 ≥2
Mtt [GeV] Multiplicity of H tags
FIG. 49 Invariant mass of the reconstructed tt̄-pair in data FIG. 50 Multiplicity of CA15 jets which fulfill the Higgs tag-
and simulation for the lepton+jets channel in the category ging criteria. The solid histograms represent the simulated
with one top-tagged jet, taken from Ref. (248). background processes (tt and QCD multijet). The hatched
error bands show the statistical uncertainty of the simulated
events. Taken from Ref. (266).
No significant excess above the predicted background
is observed in the measured tt̄ invariant mass spectrum.
Figure 49 shows the mtt̄ spectrum in the analysis cate-
The first search for VLQs in the all-hadronic final state
gory with the highest S/B fraction.
(266) targeted the T → tH decay mode. The CA algo-
Depending on the model, narrow tt̄ resonances are ex-
rithm with a large size parameter of R = 1.5 was applied
cluded for masses less than approximately 4 TeV. The
to cluster top quarks and Higgs bosons in single large
exclusion limits are weaker for scenarios with large width
jets. To identify the origin of the large CA jets a top
of the resonance.
tagging algorithm (HEPTopTagger) and a Higgs tagging
algorithm based on subjet-b tagging (see section VI.D)
are used. This was the first time these two algorithms
C. Vector-like Quarks
have been applied in a data analysis by the CMS Collab-
oration. Two subjets must be b tagged and their invari-
Vector-like quarks (VLQs) are predicted by a va-
ant mass must be greater than 60 GeV to fulfill the Higgs
riety of theories introducing a mechanism that stabi-
tagging requirement. The multiplicity of these Higgs tags
lizes the mass of the Higgs particle. Such theories in-
is shown in figure 50 which demonstrates that both the
clude little Higgs models (325; 326), models with ex-
QCD multijet and the tt backgrounds can be suppressed
tra dimensions (327; 328), and composite Higgs mod-
by several orders of magnitudes.
els (327; 328; 329). As VLQs are expected to have large
masses and have top quarks and vector-bosons as decay Extensive use of substructure methods has also been
products, jet substructure analyses have been applied in made by the ATLAS Collaboration, in particular for the
many searches for VLQs. search for single production of VLQs. The single pro-
The first search for VLQs using jet substructure duction modes may have higher cross sections than pair
methods was an inclusive search for pair-produced T production depending on the VLQ mass and the coupling
quarks (265). As VLQs may have many decay modes parameters (330). ATLAS performed an analysis (277)
(T → bW , T → tZ, T → tH, B → tW , B → bZ, where the VLQ is searched for in the decay mode with a
B → bH), a large variety of final states needs to be W boson and a top quark (B → tW ). Final states with
explored. For this reason, an inclusive search has been at least one lepton are considered, where either the W
performed without the attempt to reconstruct a specific boson or the top quark appear in a boosted configura-
decay chain. The CA algorithm was used with a dis- tion. They are identified by the application of a jet mass
tance parameter R = 0.8 (CA8 jets). Boosted W jets are requirement (m > 50 GeV) on a trimmed large-R anti-kT
identified based on the mass of the CA8 jet while boosted jet with a distance parameter R = 1.0.
top jets are identified with the CMSTT, described in sec- A different strategy is followed in another ATLAS
tion VIII.B. search (278), where the decay into the bW final state
38
Events / 5 GeV
sumed to decay leptonically, no boosted hadronic W or 2500
Data W(qq)+jets (×3)
CMS
top quark decays are present. Therefore, the analysis Total SM pred. Z(qq)+jets (×3)
uses a veto on the presence of massive (m > 70 GeV), 2000 Multijet pred. t/tt(qq)+jets (×3)
Z'(qq), g =0.17, mZ'=135 GeV
trimmed large-R anti-kT jets with R = 1.0, to suppress q'
Data/Prediction
1.2 60 80 100 120 140 160 180 200 220 240
1
0.8 mSD (GeV)
D. Leptophobic Z 0 50 100 150 200 250
mSD (GeV)
Besides resonaces coupling to heavy SM particles, there FIG. 51 Soft drop jet mass of anti-kt R = 0.8 jets in data and
exist predictions for resonances that couple to quarks and for the dominating background processes; multijet production
gluons (331; 332; 333; 334), including simplified Dark and W/Z+jets events. Taken from Ref. (160).
Matter (DM) models in which resonances couple only to
quarks and DM particles (335; 336; 337). When the new ×10
3
Events / 2 GeV
ATLAS Data
it can be boosted when produced in association with 140
W + jets
s = 13 TeV, 36.1 fb-1
initial-state radiation and thus entirely captured by a sin- 120 Z + jets
W/Z validation
gle large-radius jet (160; 281). Searching in this mode can Jet channel
Background est.
100 Bkg. stat. uncert.
significantly extend the sensitivity of the existing search
Bkg. syst. uncert.
program, where resolved low-mass resonance searches 80
typically degrade due to high trigger thresholds and the
60
enormous QCD multijet background.
Both ATLAS and CMS have used this strategy to look 40
the jet mass spectrum due to large correlation between 50 60 70 80 90 100 110
the jet mass and substructure variables, a decorrelation Large-R jet mass [GeV]
with the DDT method is applied. Data-driven techniques
are used to determine the dominating background from
QCD multijet production. Subdominant processes such FIG. 52 Trimmed jet mass distribution anti-kT R = 1.0 jets
as W/Z+jets events are estimated from MC simulation. in data and for the dominating background processes. Taken
The jet mass distributions of the large-R jet is shown in from Ref. (281).
figure 51 and 52 for the CMS and ATLAS analyses, re-
spectively. No evidence for a resonant structure on top
of the SM background is observed. methods are used to exploit the details of hadronic ac-
tivity detectable by modern particle detectors such as
ATLAS and CMS, and precision Standard Model mea-
IX. CONCLUSIONS surements and searches for physics beyond the Standard
Model at both these experiments increasingly rely on one
Jet substructure is the term used to describe the cal- or more of the tools developed by the jet substructure
culations, algorithms, and analysis techniques developed community. With increasingly sophisticated hardware
over the last decade and reviewed in this article. These and software capabilities, jet substructure techniques of
39
the future will grow in complexity and utility, further [14] J. M. Butterworth, A. R. Davison, M. Rubin and G. P.
empowering the exploration of the subnuclear properties Salam, Jet substructure as a new Higgs search channel
of nature. at the LHC, Phys. Rev. Lett. 100 (2008) 242001,
[0802.2470].
[15] D. E. Kaplan, K. Rehermann, M. D. Schwartz and
B. Tweedie, Top Tagging: A method for identifying
ACKNOWLEDGMENTS boosted hadronically decaying top quarks, Phys. Rev.
Lett. 101 (2008) 142001, [0806.0848].
Much of the work in this field in recent years has been [16] A. J. Larkoski, I. Moult and B. Nachman, “Jet
galvanized by the Boost Workshop Series (338; 339; 340), substructure at the Large Hadron Collider: A review of
which continues to inspire fruitful collaborations between recent advances in theory and machine learning.” 2017.
[17] G. Hanson, G. S. Abrams, A. M. Boyarski,
experimentalists and theorists.
M. Breidenbach, F. Bulos, W. Chinowsky et al.,
The editors thank CERN and the ATLAS and CMS Evidence for jet structure in hadron production by
Collaborations, the participants and organizers of the e+ e− annihilation, Phys. Rev. Lett. 35 (1975) 1609.
Boost Workshops held in Zurich 2016 (341) and Buffalo [18] J. D. Bjorken and S. J. Brodsky, Statistical model for
2017 (342) for discussions and input, and Jon Butter- electron-positron annihilation into hadrons, Phys. Rev.
worth for suggesting this jet substructure review article. D1 (Mar, 1970) 1416.
We also thank Andrew Larkoski and Ian Moult for the [19] S. Brandt, C. Peyrou, R. Sosnowski and
A. Wroblewski, The principal axis of jets — an
collaboration on the theoretical review.
attempt to analyse high-energy collisions as two-body
processes, Phys. Lett. 12 (1964) 57.
[20] E. Farhi, Quantum chromodynamics test for jets, Phys.
REFERENCES Rev. Lett. 39 (1977) 1587.
[21] G. Sterman and S. Weinberg, Jets from Quantum
[1] JADE Collaboration, Observation of planar three jet Chromodynamics, Phys. Rev. Lett. 39 (Dec, 1977)
events in e+ e− annihilation and evidence for gluon 1436.
bremsstrahlung, Phys. Lett. B91 (1980) 142. [22] G. P. Salam and G. Soyez, A practical seedless
[2] PLUTO Collaboration, Evidence for gluon infrared-safe cone jet algorithm, JHEP 05 (2007) 086,
bremsstrahlung in e+ e− annihilations at high energies, [0704.0292].
Phys. Lett. B86 (1979) 418. [23] JADE Collaboration, Experimental studies on multijet
[3] MARK-J Collaboration, Discovery of Three Jet Events production in e+ e− annihilation at PETRA energies,
and a Test of Quantum Chromodynamics at PETRA Z. Phys. C33 (1986) 23.
Energies, Phys. Rev. Lett. 43 (1979) 830. [24] JADE Collaboration, Experimental investigation of the
[4] TASSO Collaboration, Evidence for planar events in energy dependence of the strong coupling strength,
e+ e− annihilation at high energies, Phys. Lett. B86 Phys. Lett. B213 (1988) 235.
(1979) 243. [25] S. Catani, Y. Dokshitzer, M. Olsson, G. Turnock and
[5] D0 Collaboration, Observation of the top quark, Phys. B. Webber, New clustering algorithm for multijet cross
Rev. Lett. 74 (1995) 2632, [hep-ex/9503003]. sections in e+ e− annihilation, Phys. Lett. B269 (1991)
[6] CDF Collaboration, Observation of top quark 432.
production in p̄p collisions with the Collider Detector [26] S. Catani, Y. Dokshitzer and B. Webber, The
at Fermilab, Phys. Rev. Lett. 74 (1995) 2626, k⊥ -clustering algorithm for jets in deep inelastic
[hep-ex/9503002]. scattering and hadron collisions, Phys. Lett. B285
[7] ATLAS Collaboration, Observation of Higgs boson (1992) 291.
production in association with a top quark pair at the [27] S. Catani, Y. Dokshitzer, M. Seymour and B. Webber,
LHC with the ATLAS detector, Phys. Lett. B784 Longitudinally-invariant k⊥ -clustering algorithms for
(2018) 173, [1806.00425]. hadron-hadron collisions, Nucl. Phys. B406 (1993)
[8] CMS Collaboration, Observation of ttH production, 187.
Phys. Rev. Lett. 120 (2018) 231801, [1804.02610]. [28] J. Gallicchio and Y.-T. Chien, “Quit Using
[9] ATLAS Collaboration, Observation of H → bb̄ decays Pseudorapidity, Transverse Energy, and Massless
and V H production with the ATLAS detector, Phys. Constituents.” 2018.
Lett. B786 (2018) 59–86, [1808.08238]. [29] Y. L. Dokshitzer, G. D. Leder, S. Moretti and B. R.
[10] CMS Collaboration, Observation of Higgs boson decay Webber, Better jet clustering algorithms, JHEP 08
to bottom quarks, Phys. Rev. Lett. 121 (2018) 121801, (1997) 001, [hep-ph/9707323].
[1808.08242]. [30] M. Wobisch and T. Wengler, Hadronization corrections
[11] M. H. Seymour, Searches for new particles using cone to jet cross-sections in deep inelastic scattering, in
and cluster jet algorithms: A comparative study, Z. Monte Carlo generators for HERA physics.
Phys. C62 (1994) 127. Proceedings, Workshop, Hamburg, Germany,
[12] J. M. Butterworth, B. E. Cox and J. R. Forshaw, W W 1998-1999, p. 270, 1998, hep-ph/9907280.
scattering at the CERN LHC, Phys. Rev. D65 (2002) [31] M. Cacciari, G. P. Salam and G. Soyez, The anti-kt jet
096014, [hep-ph/0201098]. clustering algorithm, JHEP 04 (2008) 063,
[13] K. Agashe, A. Belyaev, T. Krupovnickas, G. Perez and [0802.1189].
J. Virzi, LHC Signals from Warped Extra Dimensions, [32] G. P. Salam, Towards Jetography, Eur. Phys. J. C67
Phys. Rev. D77 (2008) 015003, [hep-ph/0612015]. (2010) 637, [0906.1833].
40
[33] R. K. Ellis, W. J. Stirling and B. R. Webber, QCD performance using particle flow with the ATLAS
and collider physics, Camb. Monogr. Part. Phys. Nucl. Detector, Eur. Phys. J. C77 (2017) 466, [1703.10485].
Phys. Cosmol. 8 (1996) 1–435. [55] ATLAS Collaboration, Topological cell clustering in
[34] S. D. Ellis, Z. Kunszt and D. E. Soper, Jets at hadron the ATLAS calorimeters and its performance in LHC
colliders at order αs3 : A look inside, Phys. Rev. Lett. Run 1, Eur. Phys. J. C77 (2017) 490, [1603.02934].
69 (1992) 3615, [hep-ph/9208249]. [56] ATLAS Collaboration, Improving jet substructure
[35] OPAL Collaboration, A direct observation of quark - performance in atlas using track-caloclusters, ATLAS
gluon jet differences at LEP, Phys. Lett. B265 (1991) PUB Note ATL-PHYS-PUB-2017-015, 2017.
462. [57] CMS Collaboration, V Tagging Observables and
[36] ALEPH Collaboration, Study of the subjet structure of Correlations, CMS Physics Analysis Summary
quark and gluon jets, Phys. Lett. B346 (1995) 389. CMS-PAS-JME-14-002, 2014.
[37] DELPHI Collaboration, Energy dependence of the [58] ATLAS Collaboration, A measurement of the
differences between the quark and gluon jet calorimeter response to single hadrons and
fragmentation, Z. Phys. C70 (1996) 179. determination of the jet energy scale uncertainty using
[38] L3 Collaboration, Ks0 and Λ production in quark and LHC Run-1 pp-collision data with the ATLAS detector,
gluon jets at LEP, Phys. Lett. B407 (1997) 389. Eur. Phys. J. C77 (2017) 26, [1607.08842].
[Erratum: Phys. Lett. B427, 409 (1998)]. [59] ATLAS Collaboration, Charged-particle
√ distributions
[39] CDF Collaboration,
√ A Measurement of jet shapes in pp̄ at low transverse momentum in s = 13 TeV pp
collisions at s = 1.8 TeV, Phys. Rev. Lett. 70 (1993) interactions measured with the ATLAS detector at the
713. LHC, Eur. Phys. J. C76 (2016) 502, [1606.01133].
[40] D0 Collaboration, Transverse √ energy distributions [60] ATLAS Collaboration, Performance of the ATLAS
within jets in pp̄ collisions at s = 1.8 TeV, Phys. Transition Radiation Tracker in Run 1 of the LHC:
Lett. B357 (1995) 500. tracker properties, JINST 12 (2017) P05002,
[41] ZEUS Collaboration, Measurement of jet shapes in [1702.06473].
photoproduction at HERA, Eur. Phys. J. C2 (1998) [61] ATLAS Collaboration, A neural network clustering
61, [hep-ex/9710002]. algorithm for the ATLAS silicon pixel detector, JINST
[42] ZEUS Collaboration, Measurement of jet shapes in 9 (2014) P09009, [1406.7690].
high Q2 deep inelastic scattering at HERA, Eur. Phys. [62] ATLAS Collaboration, Performance of the ATLAS
J. C8 (1999) 367, [hep-ex/9804001]. Track Reconstruction Algorithms in Dense
[43] H1 Collaboration, Measurement of internal jet Environments in LHC Run 2, Eur. Phys. J. C77
structure in dijet production in deep inelastic (2017) 673, [1704.07983].
scattering at HERA, Nucl. Phys. B545 (1999) 3, [63] ATLAS Collaboration, Measurement of track
[hep-ex/9901010]. reconstruction inefficiencies in the core of√jets via pixel
[44] CMS Collaboration, Jet energy scale and resolution in dE/dx with the ATLAS experiment using s = 13 TeV
the CMS experiment in pp collisions at 8 TeV, JINST pp collision data, ATLAS PUB Note
12 (2017) P02014, [1607.03663]. ATL-PHYS-PUB-2016-007, 2016.
[45] ATLAS Collaboration, The ATLAS Experiment at the [64] ATLAS Collaboration, Modelling of Track
CERN Large Hadron Collider, JINST 3 (2008) S08003. Reconstruction
√ Inside Jets with the 2016 ATLAS
[46] CMS Collaboration, The CMS experiment at the s = 13 TeV pp Dataset, ATLAS PUB Note
CERN LHC, JINST 3 (2008) S08004. ATL-PHYS-PUB-2017-016, 2017.
[47] ATLAS Collaboration, The ATLAS Inner Detector [65] ATLAS Collaboration, Measurement of performance of
commissioning and calibration, Eur. Phys. J. C70 the pixel neural network√ clustering algorithm of the
(2010) 787–821, [1004.5293]. ATLAS experiment at s = 13 TeV, ATLAS PUB
[48] CMS Collaboration, Description and performance of Note ATL-PHYS-PUB-2015-044, 2015.
track and primary-vertex reconstruction with the CMS [66] CMS Collaboration, High pT jet tracking with JetCore
tracker, JINST 9 (2014) P10009, [1405.6569]. tracking and pixel cluster splitting, CMS Detector
[49] ATLAS IBL Collaboration, Production and Integration Performance Summary , 2014.
of the ATLAS Insertable B-Layer, JINST 13 (2018) [67] ATLAS Collaboration, Jet energy measurement and its
T05008, [1803.00844]. systematic
√ uncertainty in proton-proton collisions at
[50] CMS Collaboration, The CMS hadron calorimeter s = 7 TeV with the ATLAS detector, Eur. Phys. J.
project: Technical Design Report, Tech. Rep. C75 (2015) 17, [1406.0076].
CERN-LHCC-97-031, Geneva, 1997. [68] ATLAS Collaboration, Jet energy scale measurements
[51] CMS Collaboration, Particle-Flow Event and their systematic
√ uncertainties in proton-proton
Reconstruction in CMS and Performance for Jets, collisions at s = 13 TeV with the ATLAS detector,
Taus, and MET, CMS Physics Analysis Summary Phys. Rev. D96 (2017) 072002, [1703.09665].
CMS-PAS-PFT-09-001, 2009. [69] CMS Collaboration, Determination of Jet Energy
[52] CMS Collaboration, Commissioning of the Calibration and Transverse Momentum Resolution in
Particle-Flow reconstruction in Minimum-Bias and Jet CMS, JINST 6 (2011) P11002, [1107.4277].
Events from pp Collisions at 7 TeV, CMS Physics [70] A. Cukierman and B. Nachman, Mathematical
Analysis Summary CMS-PAS-PFT-10-002, 2010. Properties of Numerical Inversion for Jet Calibrations,
[53] CMS Collaboration, Particle-flow reconstruction and Nucl. Instrum. Meth. A858 (2017) 1, [1609.05195].
global event description with the CMS detector, JINST [71] ATLAS Collaboration, Jet global sequential corrections
12 (2017) P10003, [1706.04965]. with
√ the ATLAS detector in proton-proton collisions at
[54] ATLAS Collaboration, Jet reconstruction and s = 8 TeV, ATLAS CONF Note
41
W/Z+jet events, JHEP 05 (2013) 090, [1303.4811]. [129] S. Bright-Thonney and B. Nachman, Investigating the
[109] CMS Collaboration, Measurements of the differential Topology Dependence of Quark and Gluon Jets, JHEP
jet cross section as a function of the jet√mass in dijet 03 (2019) 098, [1810.05653].
events from proton-proton collisions at s = 13 TeV, [130] J. Gallicchio and M. D. Schwartz, Quark and gluon
JHEP 11 (2018) 113, [1807.05974]. tagging at the LHC, Phys. Rev. Lett. 107 (2011)
[110] J. Thaler and K. Van Tilburg, Identifying boosted 172001, [1106.3076].
objects with N-subjettiness, JHEP 03 (2011) 015, [131] A. J. Larkoski, J. Thaler and W. J. Waalewijn,
[1011.2268]. Gaining (mutual) information about quark/gluon
[111] J. Thaler and K. Van Tilburg, Maximizing boosted top discrimination, JHEP 11 (2014) 129, [1408.3122].
identification by minimizing N-subjettiness, JHEP 02 [132] C. Frye, A. J. Larkoski, J. Thaler and K. Zhou,
(2012) 093, [1108.2701]. Casimir meets poisson: Improved quark/gluon
[112] ATLAS Collaboration, Jet mass reconstruction with discrimination with counting observables, JHEP 09
the ATLAS detector in early Run 2 data, ATLAS (2017) 083, [1704.06266].
CONF Note ATLAS-CONF-2016-035, 2016. [133] P. T. Komiske, E. M. Metodiev and M. D. Schwartz,
[113] ATLAS Collaboration, Measurement√of large radius jet Deep learning in color: Towards automated
mass reconstruction performance at s = 8 TeV using quark/gluon jet discrimination, JHEP 01 (2016) 110,
the ATLAS detector, ATLAS CONF Note [1612.01551].
ATLAS-CONF-2016-008, 2016. [134] P. Gras, S. Höche, D. Kar, A. Larkoski, L. Lönnblad,
[114] ATLAS Collaboration, In-situ measurements of S. Plätzer et al., Systematics of quark/gluon tagging,
large-radius jet reconstruction performance, ATLAS JHEP 07 (2017) 091, [1704.03878].
CONF Note ATLAS-CONF-2017-063, 2017. [135] ATLAS Collaboration, Quark versus gluon jet tagging
[115] CMS Collaboration, Search for massive resonances using charged particle multiplicity with the ATLAS
decaying into W √ W , W Z, ZZ, qW , and qZ with dijet detector, ATLAS PUB Note
final states at s = 13TeV, Phys. Rev. D97 (2018) ATL-PHYS-PUB-2017-009, 2017.
072006, [1708.05379]. [136] ATLAS Collaboration, Light-quark√and gluon jet
[116] CMS Collaboration, Measurement of the jet mass √ in discrimination in pp collisions at s = 7 TeV with the
highly boosted tt events from pp collisions at s = 8 ATLAS detector, Eur. Phys. J. C74 (2014) 3023,
TeV, Eur. Phys. J. C77 (2017) 467, [1703.06330]. [1405.6583].
[117] A. J. Larkoski, G. P. Salam and J. Thaler, Energy [137] ATLAS Collaboration, Discrimination √ of light quark
correlation functions for jet substructure, JHEP 06 and gluon jets in pp collisions at s = 8 TeV with the
(2013) 108, [1305.0007]. ATLAS detector, ATLAS CONF Note
[118] A. J. Larkoski, I. Moult and D. Neill, Power counting ATLAS-CONF-2016-034, 2016.
to better jet observables, JHEP 12 (2014) 009, [138] ATLAS Collaboration, Quark versus gluon jet tagging
[1409.6298]. using jet images with the ATLAS detector, ATLAS
[119] A. J. Larkoski, I. Moult and D. Neill, Analytic boosted PUB Note ATL-PHYS-PUB-2017-017, 2017.
boson discrimination, JHEP 05 (2016) 117, [139] CMS Collaboration, Performance of quark/gluon
[1507.03018]. discrimination in 8 TeV pp data, CMS Physics
[120] I. Moult, L. Necib and J. Thaler, New angles on Analysis Summary CMS-PAS-JME-13-002, 2013.
energy correlation functions, JHEP 12 (2016) 153, [140] CMS Collaboration, Performance of quark/gluon
[1609.07483]. discrimination in 13 TeV data, CMS Detector
[121] CMS Collaboration, W and top tagging scale factors, Performance Summary CMS-DP-2016-070, 2016.
CMS Detector Performance Summary [141] CMS Collaboration, New developments for jet
CMS-DP-2017-026, 2017. substructure reconstruction in CMS, CMS Detector
[122] ATLAS Collaboration, Top and boson tagger Performance Summary CMS-DP-2017-027, 2017.
background performance - Moriond 2017, Tech. Rep. [142] CMS Collaboration, Search for a Higgs boson in the
ATL-JETM-2017-005, 2017. decay
√ channel H → ZZ (∗) → q q̄`− l+ in pp collisions at
[123] G. Altarelli and G. Parisi, Asymptotic freedom in s = 7 TeV, JHEP 1204 (2012) 036, [1202.1416].
parton language, Nucl. Phys. B126 (1977) 298. [143] CMS Collaboration, Measurement of the hadronic
[124] A. Banfi, G. P. Salam and G. Zanderighi, Infrared safe activity in events with a Z and two jets and extraction
definition of jet flavor, Eur. Phys. J. C47 (2006) 113, of the cross section for the electroweak
√ production of a
[hep-ph/0601139]. Z with two jets in pp collisions at s = 7 TeV, JHEP
[125] A. Buckley and C. Pollard, QCD-aware partonic jet 10 (2013) 062, [1305.7389].
clustering for truth-jet flavour labelling, Eur. Phys. J. [144] CMS Collaboration, Search for the standard model
C76 (2016) 71, [1507.00508]. Higgs boson produced through vector boson fusion and
[126] E. M. Metodiev and J. Thaler, Jet topics: decaying to bb, Phys. Rev. D92 (2015) 032008,
Disentangling quarks and gluons at colliders, Phys. [1506.01010].
Rev. Lett. 120 (2018) 241602, [1802.00008]. [145] CMS Collaboration, Search for direct production of
[127] P. T. Komiske, E. M. Metodiev and J. Thaler, An supersymmetric partners of the top quark√in the all-jets
operational definition of quark and gluon jets, JHEP final state in proton-proton collisions at s = 13 TeV,
11 (2018) 059, [1809.01140]. JHEP 10 (2017) 005, [1707.03316].
[128] C. Frye, A. J. Larkoski, M. D. Schwartz and K. Yan, [146] ATLAS Collaboration, Search for the Standard Model
Factorization for groomed jet substructure beyond the Higgs boson produced by vector-boson fusion and
√
next-to-leading logarithm, JHEP 07 (2016) 064, decaying to bottom quarks in s = 8 TeV pp collisions
[1603.09338].
43
√
with the ATLAS detector, JHEP 11 (2016) 112, with boson-tagged jets in pp collisions at s = 13 TeV
[1606.02181]. with the ATLAS detector, Phys. Lett. B777 (2018)
[147] ATLAS Collaboration, Search for high-mass diboson 91–113, [1708.04445].
resonances with √ boson-tagged jets in proton-proton [166] ATLAS Collaboration, Search for W W/W Z resonance
collisions at s = 8 TeV with the ATLAS detector, production
√ in `νqq final states in pp collisions at
JHEP 12 (2015) 055, [1506.00962]. s = 13 TeV with the ATLAS detector, ATLAS
[148] T. Sjostrand, S. Mrenna and P. Z. Skands, A Brief CONF Note ATLAS-CONF-2017-051, 2017.
Introduction to PYTHIA 8.1, [Link]. [167] J. Dolen, P. Harris, S. Marzani, S. Rappoccio and
178 (2008) 852–867, [0710.3820]. N. Tran, Thinking outside the ROCs: Designing
[149] T. Sjostrand, S. Mrenna and P. Z. Skands, PYTHIA Decorrelated Taggers (DDT) for jet substructure,
6.4 Physics and Manual, JHEP 05 (2006) 026, JHEP 05 (2016) 156, [1603.00027].
[hep-ph/0603175]. [168] C. Shimmin, P. Sadowski, P. Baldi, E. Weik,
[150] J. Bellm et al., Herwig 7.0/Herwig++ 3.0 release note, D. Whiteson, E. Goul et al., Decorrelated jet
Eur. Phys. J. C76 (2016) 196, [1512.01178]. substructure tagging using adversarial neural networks,
[151] M. Bahr et al., Herwig++ Physics and Manual, Eur. Phys. Rev. D96 (2017) 074034, [1703.03507].
Phys. J. C58 (2008) 639–707, [0803.0883]. [169] J. A. Aguilar-Saavedra, J. H. Collins and R. K.
[152] L. M. Dery, B. Nachman, F. Rubbo and Mishra, A generic anti-QCD jet tagger, JHEP 11
A. Schwartzman, Weakly supervised classification in (2017) 163, [1709.01087].
high energy physics, JHEP 05 (2017) 145, [170] I. Moult, B. Nachman and D. Neill, Convolved
[1702.00414]. substructure: Analytically decorrelating jet substructure
[153] E. M. Metodiev, B. Nachman and J. Thaler, observables, JHEP 05 (2018) 002, [1710.06859].
Classification without labels: Learning from mixed [171] D. Krohn, J. Thaler and L.-T. Wang, Jets with
samples in high energy physics, JHEP 10 (2017) 174, variable R, JHEP 06 (2009) 059, [0903.0392].
[1708.02949]. [172] ATLAS Collaboration, Boosted object tagging with
[154] S. D. Ellis, A. Hornig, T. S. Roy, D. Krohn and M. D. variable-R jets in the ATLAS detector, ATLAS PUB
Schwartz, Qjets: A non-deterministic approach to Note ATL-PHYS-PUB-2016-013, 2016.
tree-based jet substructure, Phys. Rev. Lett. 108 (2012) [173] ATLAS Collaboration, A new method to distinguish
182003, [1201.1914]. hadronically decaying boosted Z bosons from W bosons
[155] L. G. Almeida, S. J. Lee, G. Perez, G. F. Sterman, using the ATLAS detector, Eur. Phys. J. C76 (2016)
I. Sung and J. Virzi, Substructure of high-pT jets at the 238, [1509.04939].
LHC, Phys. Rev. D79 (2009) 074017, [0807.0234]. [174] ATLAS Collaboration, Flavor tagging with track jets
[156] J. Thaler and L.-T. Wang, Strategies to Identify in boosted topologies with the ATLAS detector, ATLAS
Boosted Tops, JHEP 07 (2008) 092, [0806.0023]. PUB Note ATL-PHYS-PUB-2014-013, 2014.
[157] J. Gallicchio and M. D. Schwartz, Seeing in color: Jet [175] CMS Collaboration, A Cambridge-Aachen (C-A) based
superstructure, Phys. Rev. Lett. 105 (2010) 022001, jet algorithm for boosted top-jet tagging, CMS Physics
[1001.5027]. Analysis Summary CMS-PAS-JME-09-001, 2009.
[158] ATLAS Collaboration, Identification of [176] CMS Collaboration, Jet substructure algorithms, CMS
hadronically-decaying W bosons and top quarks using Physics Analysis Summary CMS-PAS-JME-10-013,
high-level features as input to boosted decision
√ trees 2011.
and deep neural networks in ATLAS at s = 13 TeV, [177] T. Plehn, G. P. Salam and M. Spannowsky, Fat jets
ATLAS PUB Note ATL-PHYS-PUB-2017-004, 2017. for a light Higgs, Phys. Rev. Lett. 104 (2010) 111801,
[159] A. J. Larkoski, D. Neill and J. Thaler, Jet shapes with [0910.5472].
the broadening axis, JHEP 04 (2014) 017, [1401.2158]. [178] T. Plehn, M. Spannowsky, M. Takeuchi and
[160] CMS Collaboration, Search for low mass vector D. Zerwas, Stop reconstruction with tagged tops, JHEP
resonances decaying into quark-antiquark
√ pairs in 10 (2010) 078, [1006.2833].
proton-proton collisions at s = 13 TeV, JHEP 01 [179] ATLAS Collaboration, Performance of jet substructure
(2018) 097, [1710.00159]. techniques
√ for large-R jets in proton-proton collisions
[161] CMS Collaboration, Search for massive resonances at s = 7 TeV using the ATLAS detector, JHEP 09
decaying into pairs √ of boosted bosons in semi-leptonic (2013) 076, [1306.4945].
final states at s = 8 TeV, JHEP 08 (2014) 174, [180] D. E. Soper and M. Spannowsky, Finding physics
[1405.3447]. signals with shower deconstruction, Phys. Rev. D84
[162] CMS Collaboration, Search for massive W H √ (2011) 074002, [1102.3480].
resonances decaying into the `νbb final state at s = 8 [181] D. E. Soper and M. Spannowsky, Finding top quarks
TeV, Eur. Phys. J. C76 (2016) 237, [1601.06431]. with shower deconstruction, Phys. Rev. D87 (2013)
[163] CMS Collaboration, Search for massive resonances 054012, [1211.3140].
decaying into√WW, WZ or ZZ bosons in proton-proton [182] ATLAS Collaboration, Search for W 0 → tb decays √ in
collisions at s = 13 TeV, JHEP 03 (2017) 162, the hadronic final state using pp collisions at s = 13
[1612.09159]. TeV with the ATLAS detector, Phys. Lett. B781
[164] ATLAS Collaboration, Search for heavy resonances (2018) 327, [1801.07893].
decaying to a W or Z boson and a Higgs √ boson in the [183] ATLAS Collaboration, Performance of top quark and
q q̄ (0) bb̄ final state in pp collisions at s = 13 TeV with W boson tagging in Run 2 with ATLAS, ATLAS
the ATLAS detector, Phys. Lett. B774 (2017) 494–515, CONF Note ATLAS-CONF-2017-064, 2017.
[1707.06958]. [184] G. Kasieczka, T. Plehn, T. Schell, T. Strebler and
[165] ATLAS Collaboration, Search for diboson resonances G. P. Salam, Resonance searches with an updated top
44