Measurement of jet-substructure observables in top quark, W boson and light jet production in proton-proton collisions at p s = 13 TeV with the ATLAS detector

: A measurement of jet substructure observables is presented using data collected in 2016 by the ATLAS experiment at the LHC with proton-proton collisions at p s = 13 TeV. Large-radius jets groomed with the trimming and soft-drop algorithms are studied. Dedicated event selections are used to study jets produced by light quarks or gluons, and hadronically decaying top quarks and W bosons. The observables measured are sensitive to substructure, and therefore are typically used for tagging large-radius jets from boosted massive particles. These include the energy correlation functions and the N -subjettiness variables. The number of subjets and the Les Houches angularity are also considered. The distributions of the substructure variables, corrected for detector e(cid:11)ects, are compared to the predictions of various Monte Carlo event generators. They are also compared between the large-radius jets originating from light quarks or gluons, and hadronically decaying top quarks and W bosons.


Introduction
Increasing the centre-of-mass energy of proton-proton (pp) collisions from 7 and 8 TeV in Run 1 to 13 TeV in Run 2 of the Large Hadron Collider (LHC) leads to a larger fraction of heavy particles such as top quarks, vector bosons and Higgs bosons being produced with large transverse momenta.This large transverse momentum leads to collimated decay products.They are usually reconstructed in a large-radius jet, whose internal (sub)structure shows interesting features that can be used to identify the origin of the jet [1,2].
Grooming algorithms are employed to reduce the effects of soft, uncorrelated radiation contained in the large-radius jet.Jet substructure observables are then used to classify the jets originating from different particles.This process is called jet tagging and the algorithms are referred to as taggers.Most of the grooming algorithms and jet substructure observables were developed on the basis of theoretical calculations or Monte Carlo (MC) simulation programs and then they are applied to data in searches for new physics signatures.Thus it is important to measure these observables and compare the results with predictions from MC simulation.This will also help in optimising and developing current and future substructure taggers.
The ATLAS Collaboration has performed measurements of jet mass and substructure variables at the pp centre-of-mass energies of √ s = 7, 8 and 13 TeV [3][4][5][6][7][8][9] in inclusive jet events, and the CMS Collaboration has performed measurements of jet mass in dijet, W/Z boson, and t t events [10][11][12] at √ s = 7 and 13 TeV.This paper presents measurements of substructure variables in large-radius jets produced in inclusive dijet events and in t t events at √ s = 13 TeV using 33 fb −1 of data collected in 2016 by the ATLAS experiment.In this analysis, the lepton+jets decay mode of t t events is selected, where one W boson decays into a muon and a neutrino, and the other W boson decays into a pair of quarks.Then the large-radius jets are separated into those that contain all the decay products of a hadronically top quark and those containing only hadronic W boson decay products.
The contents of this paper are organised as follows.First, a description of the ATLAS detector is presented in Section 2 and then the MC samples used in the analysis are discussed in Section 3. In Section 4, event and object selections are summarised.The measured jet substructure observables are defined in Section 5.The background estimation is described in Section 6 and the systematic uncertainties are assessed in Section 7. In Section 8, detector-level mass and p T distributions corresponding to selected large-radii jets are shown, and the unfolding is described in Section 9. Finally, the unfolded results are presented in Section 10, and the conclusions in Section 11.

ATLAS detector
The ATLAS experiment uses a multipurpose particle detector [13,14] with a forward-backward symmetric cylindrical geometry and a near 4π coverage in solid angle.1It consists of an inner tracking detector (ID) surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, electromagnetic (EM) and hadron calorimeters, and a muon spectrometer.The ID consists of silicon pixel, silicon microstrip, and straw-tube transition-radiation tracking detectors, covering the pseudorapidity range |η| < 2.5.The 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe.The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards.Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the z-axis.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).An angular separation between two objects is defined as ∆R ≡ (∆η) 2 + (∆φ) 2 , where ∆η and ∆φ are the separations in η and φ.Momentum in the transverse plane is denoted by p T .
calorimeter system covers the pseudorapidity range |η| < 4.9.Electromagnetic calorimetry is performed with barrel and endcap high-granularity lead/liquid-argon (LAr) sampling calorimeters, within the region |η| < 3.2.There is an additional thin LAr presampler covering |η| < 1.8, to correct for energy loss in material upstream of the calorimeters.For |η| < 2.5, the LAr calorimeters are divided into three layers in depth.Hadronic calorimetry is performed with a steel/scintillator-tile calorimeter, segmented into three barrel structures within |η| < 1.7, and two copper/LAr hadronic endcap calorimeters, which cover the region 1.5 < |η| < 3.2.The forward solid angle up to |η| = 4.9 is covered by copper/LAr and tungsten/LAr calorimeter modules, which are optimised for energy measurements of electrons/photons and hadrons, respectively.The muon spectrometer consists of separate trigger and high-precision tracking chambers that measure the deflection of muons in a magnetic field generated by superconducting air-core toroids.
The ATLAS detector selects events using a tiered trigger system [15].The first level is implemented in custom electronics.The second level is implemented in software running on a general-purpose processor farm which processes the events and reduces the rate of recorded events to 1 kHz.

Monte Carlo samples
Simulated events are used to optimise the event selection, correct the data for detector effects and estimate systematic uncertainties.The predictions of different phenomenological models implemented in the Monte Carlo (MC) generators are compared with the data corrected to the particle level (i.e.observables constructed from final-state particles within the detector acceptance).
The generators used to produce the samples are listed in Table 1.The dijet, t t and single-top-quark samples are considered to be signal processes in this analysis, corresponding to the dedicated selections.The background is estimated using Z/W+jets and diboson samples.The tt samples are scaled to next-to-next-to-leading order (NNLO) in perturbative QCD, including soft-gluon resummation to nextto-next-to-leading-log order (NNLL) [16] in cross-section, assuming a top quark mass m t = 172.5 GeV.The P model [17] resummation damping parameter, h damp , which controls the matching of matrix elements to parton showers and regulates the high-p T radiation, was set to 1.5m t .The single-top-quark [18][19][20][21][22][23] and W/Z samples [24] are scaled to the NNLO theoretical cross-sections.
The predicted shape of jet substructure distributions depends on the modelling of final-state radiation (FSR), and fragmentation and hadronisation, as well as on the merging/matching between matrix element (ME) and parton shower (PS) generators.The P 8 and the S generators use a dipole shower ordered in transverse momentum, with the Lund string [25] and cluster hadronisation model [26] respectively.The H 7 generator uses an angle-ordered shower, with the cluster hadronisation model.For comparison purposes in dijet events, a sample was generated with S using the string hadronisation model.
The MC samples were processed through the full ATLAS detector simulation [41] based on Geant4 [42], and then reconstructed and analysed using the same procedure and software that are used for the data.Additional pp collisions generated by P 8, with parameter values set to the A2 tune [43] and using the MSTW2008 [44] PDF set, were overlaid to simulate the effects of additional collisions from the same and nearby bunch crossings (pile-up), with a distribution of the number of extra collisions matching that of data.

Object and event selection
This analysis uses pp collision data at √ s = 13 TeV collected by the ATLAS detector in 2016, that satisfy a number of criteria to ensure that the ATLAS detector was in good operating condition.All selected events must have at least one vertex with at least two associated tracks with p T > 400 MeV.The vertex with the highest p 2 T,track , where p T,track is the transverse momentum of a track associated with the vertex, is chosen as the primary vertex.
Jets are reconstructed from the EM-scale or locally-calibrated topological energy clusters [45] in both the EM and hadronic calorimeters using the anti-k t algorithm [46] with a radius parameter of R = 0.4 or R = 1.0, referred to as small-radius and large-radius jets respectively.These clusters are assumed to be massless when computing the jet four-vectors and substructure variables.A trimming algorithm [47] is employed for the large-radius jets to mitigate the impact of initial-state radiation, underlying-event activity, and pile-up.Trimming removes subjets of radius R sub = 0.2 with p i T /p jet T < f cut , where p i T is the transverse momentum of the i th subjet, p jet T is the transverse momentum of the jet under consideration, and f cut = 0.05.All large-radius jets used in this paper are trimmed before applying the selection criteria.The energies of jets are calibrated by applying p T -and rapidity-dependent corrections derived from Monte Carlo simulation with additional in situ correction factors determined from data [48].
In order to reduce the contamination by small-radius jets originating from pile-up, a requirement is imposed on the output of the Jet Vertex Tagger (JVT) [49].The JVT algorithm is a multivariate algorithm that uses tracking information to reject jets which do not originate from the primary vertex, and is applied to jets with p T < 60 GeV and |η| < 2.4.Small-radius jets containing b-hadrons are tagged using a neural-network-based algorithm [50][51][52] that combines information from the track impact parameters, secondary vertex location, and decay topology inside the jets.The operating point corresponds to an overall 70% b-tagging efficiency in simulated t t events, and to a probability of mis-tagging light-flavour jets of approximately 1%.
Muons are reconstructed from high-quality muon spectrometer track segments matched to ID tracks.Muons with a transverse momentum greater than 30 GeV and within |η| < 2.5 are selected if the associated track has a longitudinal impact parameter |z 0 sin(θ)| < 0.5 mm and a transverse impact parameter significance |d 0 |/σ(d 0 )| < 3. The impact parameter d 0 is measured relative to the beam line.The muon candidates are also required to be isolated from nearby hadronic activity [53].The muon isolation criteria remove muons that lie a distance ∆R(µ, jet) < 0.04 + 10 GeV/p T,µ from a small-radius jet axis, where p T,µ is the p T of the muon.Since muons deposit energy in the calorimeters, an overlap removal procedure is applied in order to avoid double counting of leptons and small-radius jets.
Electrons are reconstructed from energy deposits measured in the EM calorimeter which are matched to ID tracks.They are required to be isolated from nearby hadronic activity by using a set of p T -and η-dependent criteria based on calorimeter and track information as described in Ref. [54].Their selection also requires p T > 30 GeV and |η| < 2.5, excluding the region 1.37 < |η| < 1.52 which corresponds to the transition region between the barrel and end-cap calorimeters.Photon candidates are reconstructed from clusters of energy deposited in the EM calorimeter, and must have p T > 30 GeV and |η| < 2.5.Photon identification is based primarily on shower shapes in the calorimeter [55].
The missing transverse momentum, with magnitude E miss T , is calculated as the negative vectorial sum of the transverse momenta of calibrated photons, electrons, muons and jets associated with the primary vertex [56].The transverse mass of the leptonically decaying W boson, m W T , is defined using the absolute value of In order to examine large-radius jets originating from light quarks and gluons, from top quarks and from W bosons, three event selections are defined.These are referred to as dijet, top and W selections, and are indicative of the origin of the large-radius jet.
In the dijet selection, the events are accepted by a single-large-radius-jet trigger that becomes fully efficient for jets with p T > 400 GeV.The offline dijet selection requires a leading trimmed large-radius jet with p T > 450 GeV and |η| < 1.5, and at least one other trimmed large-radius jet with p T > 200 GeV and |η| < 2.5, and rejects the event if an electron or muon is present.
For both the top and W selections, events are collected with a set of single-muon triggers that become fully efficient for muon p T > 28 GeV.The top quarks and the W bosons are identified from their decay products.A geometrical separation between the decay products of the two top quark candidates is required.Additional requirements are applied to separate large-radius jets containing all decay products of the top quark from those where the large-radius jet only contains the hadronic W boson decays, with the b-tagged small-radius jet reconstructed independently.These form the top selection and the W selection respectively.The selections are described in Table 2.After these requirements the data sample contains about 3.2 × 10 7 events in the dijet selection, and roughly 6800 and 4500 events in the top and W selection respectively.
Particle-level observables in Monte Carlo simulation are constructed from stable particles, defined as those with proper lifetimes cτ 10 mm.Muons at particle level are dressed by including contributions from photons with an angular distance ∆R < 0.1 from the muon.Particle-level jets do not include muons or neutrinos.Particle-level b-tagging is performed by requiring a prompt b-hadron to be ghost-associated [57] with the jet.Table 2: Summary of object event selections for detector-level and particle-level dijet and t t events."Leptonic top" refers to the top quark that decays into a leptonically decaying W boson, while "b-tagged jet" refers to small-radius jets that pass a b-tagging requirement.The top and W selections are common up to the requirement on the leptonic top, then they differ on the requirements on the leading-p T trimmed large-radius jet.All selections are inclusive, unless otherwise mentioned.

Definition of the jet observables
All large-radius jets are trimmed before being used in the selections, and subsequently only the leading trimmed large-radius is considered in the analysis.Then the large-radius jet constructed from the original constituents of the selected jet before the trimming step is groomed using the soft-drop algorithm, and the jet substructure observables studies are constructed from that soft-dropped large-radius jet.
Soft-drop [58,59] is an extension of the original split-filtering technique [60] and relies on reclustering the jet constituents using the angle-ordered Cambridge-Aachen jet algorithm and then sequentially considering each splitting in order to remove soft and wide-angle radiation.At each step the jet is split into two proto-jets.The removal of proto-jets in a splitting is controlled by two parameters: a measure of the energy balance of the pair, z cut , and the significance of the angular separation of the proto-jets, β SD .These are used to define the soft-drop condition: where R 12 is the angular distance between the two proto-jets and R is the radius of the large jet.In this analysis, values of z cut = 0.1 and β SD = 0.0 are used, based on previous ATLAS studies [8].An important feature of soft-drop is that groomed observables are analytically calculable to high-order resummation accuracy [61][62][63].
The following substructure variables are measured in this analysis: • Number of subjets with p T > 10 GeV, reconstructed from the selected large-radius jet constituents using the k t algorithm [64] with R = 0.2.
• Generalised angularities defined as: , where z i is the transverse momentum fraction p i T /p jet T of the constituent, and θ i is the angle of the i th constituent of the jet relative to the jet axis, normalised by the jet radius.The exponents κ and β LHA probe different aspects of the jet fragmentation.The (κ = 1, β LHA = 0.5) variant is termed the Les Houches angularity (LHA) [65] and used in this analysis.It is an infrared-safe version of the jet-shape angularity, and provides a measure of the broadness of a jet.
• Energy correlation functions ECF2 and ECF3 [66], and related ratios C 2 , D 2 [67].The 1-point, 2-point and 3-point energy correlation functions for a jet J are given by: , where the parameter β ECF weights the angular separation of the jet constituents.In the above functions, the sum is over the i constituents in the jet J, such that the 1-point correlation function ECF1 is approximately the jet p T .Likewise, if one takes β ECF = 2, the 2-point correlation functions scale as the mass of a particle undergoing a two-body decay in collider coordinates.In this analysis, β ECF = 1 is used, and for brevity, β ECF is not explicitly mentioned hereafter.
The ratios of some of these quantities (written in an abbreviated form) are defined as : The observables e 2 and e 3 are measured, and are later referred to as ECF2 norm and ECF3 norm .These ratios are then used to generate the variable C 2 [66], and its modified version D 2 [65,67], which have been shown to be particularly useful in identifying two-body structures within jets [68].The C 2 and D 2 variables as defined below are measured in this analysis: • Ratios of N-subjettiness [69], τ 21 and τ 32 .The N-subjettiness describes to what degree the substructure of a given jet is compatible with being composed of N or fewer subjets.
In order to calculate τ N , first N subjet axes are defined within the jet by using the exclusive k t algorithm, where the jet reconstruction continues until a desired number of jets are found.The 0-, 1-, 2-,and 3-subjettiness are defined as: where ∆R is the angular distance between constituent i and the jet axis, a i , and ∆R a,n is the angular distance between constituent i and the axis of the n th subjet.The term R is the radius parameter of the jet.The parameter β NS gives a weight to the angular separation of the jet constituents.In the studies presented here, the value of β NS = 1 is used.In the above functions, the sum is performed over the constituents i in the jet J, and a normalisation factor τ 0 (Eq.( 1a)) is used.The ratios of the N-subjettiness functions, τ 21 = τ 2 /τ 1 and τ 32 = τ 3 /τ 2 have been shown to be particularly useful in identifying two-body and three-body structures within jets.
Studies presented in Ref.
[70] have shown that an alternative axis definition can increase the discrimination power of these variables.The winner-takes-all (WTA) axis uses the direction of the hardest constituent in the subjet obtained from the exclusive k t algorithm instead of the subjet axis, such that the distance measure ∆R a 1 ,i changes in the calculation.In this analysis, the same observables calculated with the WTA axis definition, τ WTA 21 and τ WTA 32 , are used.

Data-driven background estimation
The largest non-t t contributions to the W and top selections come from the W+jets and single-top processes.
Additionally non-prompt and mis-reconstructed muons are a separate source of background for the top and W selections. Contributions from other processes were considered and found to be negligible.A data-driven method, following Ref.
[71], is used to estimate the contribution from the W+jets process while the single-top process is considered part of the signal.
At the LHC the production rate of W + +jets events is larger than that of W − +jets due to the higher density of u-quarks than d-quarks in the proton.This results in more events with positively charged leptons.Other processes do not contribute significantly to this charge asymmetry.The data are used to derive scale factors that correct the normalisation and flavour fraction given by the MC simulation [72].
Normalisation scale factors are determined by comparing the charge asymmetry in data with the asymmetry estimated by simulation.Contributions to the asymmetry from other processes are estimated by simulation and subtracted.A selection that contains the full top and W selection criteria without any b-tagging requirements is initially used.The total number of W+jets events in data, N W + + N W − , is given by where r MC is the ratio of the number of events with positive muons to the number of events with negative muons obtained from the MC simulation while D + and D − are the number of events with positive and negative muons in data, respectively, after using simulation to subtract the estimated background contribution of all processes other than W+jets.From the above equation the scale factor C A is extracted which is defined as the ratio of W+jets events evaluated from data to the number predicted by the simulation.
Scale factors correcting the relative fractions of W bosons produced in association with jets of different flavour are also estimated using data.The fractions of W + b b, W + c c, W + c and W+light-quark events are initially estimated from simulation in a selection without the b-tagging requirements, which corresponds to the selection mentioned in Table 2 without the ∆R requirement imposed during the top and W selections.
A system of three equations is used to fit the fractions estimated from simulation to the selection with full b-tagging requirements: where f bb , f cc , f c and f light are flavour factors estimated from simulation while K bb , K cc , K c and K light are the respective correction factors.The corresponding number of events estimated by simulation with positive (negative) leptons are given by and N +(−) light .The terms D W ± are the expected numbers of W+jets events with positively or negatively charged leptons in the data.An iterative process is used to find the K flavour correction factors which are used to correct the associated f flavour fractions used in the calculation of C A .The correction factors are determined by inverting Eq. ( 2) and then the process is repeated with a new C A calculated using the corrected flavour fractions.This process is repeated 10 times and further iterations produce negligible changes in C A .This process is repeated individually for all variables in the top and W selections since, depending on the substructure of the selected large-radius jet, events can fall out of the acceptance for a subset of the variables.The final calculated scale factors are, however, consistent across both selections.These scale factors are 0.84 ± 0.02, where the uncertainty is statistical, and the overall contribution to the final selections is shown in Table 3.In order to determine the uncertainty in the shape of the subtracted W+jets distribution, the contribution from an alternative MC generator (MG5_aMC@NLO+P 8 as opposed to default S ) was used.Both MC samples were scaled to the estimated number of events and the envelope of the shape difference was taken as an uncertainty.
There is also a contribution from events where a jet is misreconstructed as a muon or when a non-prompt muon is misidentified as a prompt muon which satisfies the selection criteria.This contribution is estimated using the matrix method, comparing the yields of muons and non-prompt muons that pass a loose selection with the yields of those that pass a tight selection.The efficiency for real muon selection (ε real ) is measured using a tag-and-probe method with muons from Z → µµ events.The efficiency for misreconstructed muon selection (ε fake ) is measured in control regions dominated by background from multijet processes, after using simulation to subtract the contribution of other processes.Event weights are computed using the above efficiencies, which are parameterised in the kinematics of the event.The weight for event i, where the muons satisfy the loose criteria, is given by where δ i equals unity if the muon in event i satisfies the tight criteria and zero otherwise.The background estimate in a given bin is therefore the total sum of weights in that bin.The estimated contributions to the yield from misreconstructed or non-prompt muons for the top and W selections are shown in Table 3.These corrections have very little effect on the shape of the distributions considered.

Large-radius jet uncertainties
As jets are built from topological clusters reconstructed in the calorimeter, systematic uncertainties in the jet substructure observables are calculated using a bottom-up approach applied to the clusters forming each jet [8].The following components of the uncertainty are considered: • Cluster reconstruction efficiency (CE): Accounts for low energy particles that fail to seed a cluster based on the fraction of inner-detector tracks matched to no clusters in low µ data.The uncertainty is the observed difference between simulation and data.Since the efficiency reaches 100% for cluster energy above 2.5 GeV, no uncertainty is assumed above this value.
• Cluster energy scale variation (CESu/CESd): The cluster energy scale is determined by studying clusters matched to isolated tracks in data events with low pile-up.A fit of the E/p distribution is used to extract an overall energy scale.The uncertainty in the scale is given by taking the difference of the ratio of the scales calculated in data and simulation from unity.Clusters are independently scaled up and down and the resulting variations in observables are added in quadrature.
• Cluster energy smearing (CES): The difference in quadrature of the width of the E/p distribution measured in data and given by simulation is defined as the uncertainty in the energy resolution.The cluster energies are smeared by this value and the effect on the observables is taken as an uncertainty.
• Cluster angular resolution (CAR): The radial distance between clusters and their matched tracks (extrapolated to the corresponding calorimeter layer) is measured in bins of η and as a function of E, to account for the resolution in various regions of the calorimeter.A conservative uncertainty of 5 mrad is used to smear cluster positions.
Uncertainties in the jet p T and mass are derived by the Rtrk method [73], comparing the variables calculated using the energy deposited in the calorimeter with those using the momenta of charged-particle tracks.The largest effect on the majority of measured distributions comes from cluster energy smearing for the top and W selections, typically around 8% but can be as high as 16% in some regions.The other cluster uncertainty components contribute between 1% and 6% in the statistically significant part of the distributions for the top and W selections.For the dijet selection, the typical values are between 2% and 4% for all observables, but reach 10% in some bins.The dominant large-radius jet uncertainties for a subset of variables are shown in Figure 1.
In addition to the above uncertainties the sensitivity of the measured distributions to other detector effects was considered.This are summarised as follows: • Energy scaling correlation scheme: applying the variations to clusters with different kinematics and with different properties, assuming them to be uncorrelated.
• Since the cluster energy calibration is based on pion energy deposition, additional tests are carried out to account for the different energy deposited by non-pion hadrons, such as K L , and the impact on the distributions under study.
• Cluster merging and splitting: topo-clusters can be split or merged during the clustering procedure and this process can be sensitive to noise fluctuations.
In all cases, very conservative variations were applied in order to ensure that the distributions considered were not sensitive to the above effects.For the majority of the distributions the observed variations due to other detector effects were smaller than the cluster uncertainties.However, it was found that N-subjettiness variables in the dijet selection had shifts of about 50% when some of the cluster merging and splitting variations were applied.While these variations were conservative, in order to ensure that no systematic uncertainties are being underestimated the N-subjettiness variables and their ratios were not used in the dijet selection.

Other sources of uncertainties
Systematic uncertainties are also derived for other reconstructed objects which are considered in the top and the W selections [74].Uncertainties associated with small-radius jets, b-tagged jets, reconstructed muons and E miss T are all considered and are found to be subdominant.The theory normalisation uncertainties are also found to be negligible.
Finally, uncertainties in the shape of the subtracted W+jets component are derived by comparing, for each variable, the shapes obtained using the nominal MC sample and an alternative sample, as listed in Table 1.The envelope is taken as an uncertainty in the subtracted shape, and results in uncertainties which are smaller than 1%.The uncertainties due to signal modelling in MC generators are accounted for in unfolding, as described in Section.9.

Detector-level results
The distributions of the trimmed large-radius jet mass and p T at detector level are shown in Figure 2 for dijet, top and W selections.The peaks in the distributions due to the top and W masses are clearly visible.In general, good agreement is observed between data and simulation for the distribution of transverse momenta, while a shift is observed for the distributions of mass.This is a known effect [2], due to the lack of in situ calibrations of jet mass, and to jet mass scale uncertainties in the detector-level plots.

Unfolding
The measured distributions are unfolded to correct for detector effects.The Iterative Bayesian (IB) unfolding method [75] with three iterations (as implemented in RooUnfold [76]) is used to correct detector-level data to particle level, as defined in Section 4. Response matrices (a ji ) for each distribution are derived from MC simulation and used in order to estimate the probability for a given event at particle level (T), contributing to bin i, to be reconstructed in a given detector-level (D) bin j, also defined as P(D j |T i ).Rather than using a simple matrix inversion, IB unfolding uses a probabilistic approach.In order to do this, the unfolding matrix (θ i j ) is defined such that the number of events in a particle-level bin, T i , is given by where d j is the number of data events measured in bin j.Using Bayes' theorem, one can define the unfolding matrix as: .
where P(T i ) is the input prior.The unfolding matrix can therefore be constructed using the response matrix obtained from simulation.After corrections are applied for detector acceptance and reconstruction efficiency, Eq. ( 3) can be used to perform the unfolding.To ensure that the final distributions are not biased by the shape predicted by simulation the process is iterated, each subsequent iteration using the previous estimate for the final corrected distribution as P(T i ).The number of iterations is chosen such that differences between multiple subsequent iterations are smaller than data-driven cross-closure uncertainties, described below.
The consistency of the unfolding procedure was tested using several closure and cross-closure tests.
• MC closure: a test where the distributions from the nominal MC generator are unfolded using the nominal method.Uncertainties are found to be negligible.
• Cross-closure: accounts for modelling differences between two different MC generators.The distributions from an alternative generator are unfolded using the nominal method and the differences account for differences in the predicted shape.These result in the largest uncertainties and are typically around 5% in the dijet selection and around 14% in the top and W selections, depending on the observable and the bin.
• Data-driven cross-closure: accounts for the sensitivity of the unfolding method to differences between the shape of the observable seen in data and in simulation.The detector-level and particle-level distributions are reweighted such that the corresponding detector-level distributions match the data.These reweighted distributions are unfolded using the nominal method and uncertainties are estimated as the differences between the reweighted particle-level and unfolded distributions.

Particle-level results
The results are presented in two sets of distributions: substructure observables in data are compared with MC predictions.and distributions measured in data corresponding to different selections are compared with each other.For the latter, it must be noted that the comparisons are performed in different large-radius jet p T ranges; however, in each instance the most inclusive selection is used.They are indicative of different substructures of the large-radius jets according to their origin even with somewhat different kinematic ranges.All plots with soft-drop grooming are shown; the trimmed versions have very similar characteristics.
In Figure 3, the subjet multiplicity inside the large-radius jets from the three different selections is compared with different MC predictions, and the data are compared between the three selections.While for the dijet selection most events have one subjet, for the top selection and W selection the distributions peak at three and two subjets respectively, as expected.In both cases a non-negligible fraction of events have more subjets, indicating the presence of semi-hard gluon radiation.In the W selection, the instances with one subject are few, while for the top selection, some fraction of events have two subjets, indicating non-containment of the top quark decay products.For the dijet selection, P 8 and S describe the data the best, while for the top selection and W selection, there is more spread among MC predictions.Predictions from H 7 are very different from data for the dijet selection, a trend which is consistent across all observables.The difference between the different hadronisation models used in S is negligible.Although these observables depend on hadronisation modelling, it can be inferred that both models can be tuned to give a good description of data.
In Figure 4, the Les Houches angularity (LHA) is compared between large-radius jets for the three selections and with MC model predictions.For the dijet selection, all models except H 7 describe the data, while for the top and W selections, the level of agreement between all models and data is worse, and the peaks of the distributions in the models are shifted relative to those in data.While in the case of the top and W selections the shapes are similar, the distribution for the dijet selection peaks at the lowest value.This indicates that the additional radiation in quark/gluon jets is soft, with little activity away from the large-radius jet axis, while for the large-radius jets from top quarks and W bosons, there are hard emissions separated by appreciable angles.
In Figure 5, a comparison of C 2 among the three different selections with MC is presented, as well as a comparisons of data and MC predictions for each selection.For the dijet selection, all models except H 7 describe the data well, while for the top and W selections, the models predict shapes that differ from data, with P +H 7 performing somewhat worse than the rest.The three distributions have distinct peaks, corresponding to their substructure.The value of C 2 increases as the number of subjets inside the large-radius jets increases.
In Figure 6, comparisons of the data with MC predictions for D 2 reveal some interesting features.For the dijet selection, most of the models describe the data well, and for the top selection the some differences can be seen.For the W selection, all MC predictions have a peak shifted relative to data, suggesting that the models are overestimating gluon radiation.The distributions in data for the three selections are also compared in Figure 6 (bottom right), where peaks at different values are observed.
The distributions of ECF2 norm , as shown in Figure 7 for the different selections, can discriminate between events with two and three prong decays as opposed to one prong decay.Similarly to C 2 , for the dijet selection, all models except H 7 describe the data well, while for the top and W selections, the models predict shapes that differ somewhat from data, with agreement being worse for the W selection case.
The modelling of ECF3 norm in the dijet selection is better for P 8 than for the other generators, as shown in Figure 8.For the top and W selections, none of the models describe the shape of the data distribution well, with noticeable differences at low values.The three different selections again show distinct shapes.(bottom right) in data are compared between the soft-dropped large-radius jets for the two selections mentioned above.The subscript WTA indicates that WTA axis was used in calculating these observables.The shaded bands represent the total uncertainty, while the error bars show the statistical uncertainty, except in the bottom plots, where the shaded areas represent the total uncertainty.

Conclusions
A measurement of jet substructure observables using groomed large-radius jets from light quarks or gluons, hadronically decaying top quarks and W bosons is presented using 33 fb −1 of √ s = 13 TeV proton-proton collision data taken with the ATLAS detector at the LHC.The data discriminate between the various MC models probed.In general, P 8 for light-quark/gluon large-radius jet observables, and P +P 8, S as well as MG5_aMC@NLO+P 8 for top quark and W boson large-radius jet observables, describe the data better than other models.The different hadronisation models in S in the djiet selection result in similar predictions.For most observables, H 7 in the dijet selection, and P +H 7 in the top and W selections do not describe the data well.These measurements will be useful in improving the modelling of these substructure variables in MC generators.Since searches that utilise boosted topologies use these observables, or combinations of them, in tagging large-radius jets, a better modelling of them will help to increase the sensitivity of such searches.[74] ATLAS Collaboration, Search for heavy particles decaying into top-quark pairs using lepton-plus-jets events in proton-proton collisions at

Figure 1 :
Figure 1: Bin-by-bin systematic uncertainties due to large-radius jet reconstruction uncertainties associated with cluster, Rtrk and jet mass calibrations in the dijet (top) and W (bottom) selections for the soft-drop groomed Les Houches angularity variable (left) and the normalised ECF2 variable (right).

Figure 2 :
Figure 2: Comparison of detector-level distributions in data and MC simulation for trimmed large-radius jets for dijet (top row), top (middle row), and W (bottom row) selections.For the top and W selections, jet mass requirements have not been applied.The mass is shown in the left column, while the transverse momentum is in the right column.The shaded bands represent the combined statistical and systematic uncertainty.Contributions from dominant backgrounds are shown for the top and W selections, while the smaller contributions from other processes are grouped under other backgrounds.

Figure 3 :
Figure 3: Subjet multiplicity distributions compared with different MC predictions for soft-dropped large-radius jets from dijet (top left), top (top right), and W (bottom left) selections.For the dijet selection, S is tested with two different hadronisation models.Data are compared between the soft-dropped large-radius jets for the three selections mentioned above (bottom right).The shaded bands represent the total uncertainty, while the error bars show the statistical uncertainty, except in the bottom right plot, where the shaded areas represent the total uncertainty.

Figure 4 :Figure 5 :Figure 6 :
Figure 4: Les Houches angularity is compared with different MC predictions for soft-dropped large-radius jets from dijet (top left), top (top right), and W (bottom left) selections.For the dijet selection, S is tested with two different hadronisation models.Data are compared between the soft-dropped large-radius jets for the three selections mentioned above (bottom right).The shaded bands represent the total uncertainty, while the error bars show the statistical uncertainty, except in the bottom right plot, where the shaded areas represent the total uncertainty.

Figure 7 :Figure 8 :
Figure 7: The distributions of ECF2 norm compared with different MC predictions for soft-dropped large-radius jets from dijet (top left), top (top right), and W (bottom left) selections.For the dijet selection, S is tested with two different hadronisation models.Data are compared between the soft-dropped large-radius jets for the three selections mentioned above (bottom right).The subscript norm indicates that normalised versions of ECF2 norm are used.The shaded bands represent the total uncertainty, while the error bars show the statistical uncertainty, except in the bottom right plot, where the shaded areas represent the total uncertainty.

Figure 9 :
Figure 9: The distributions of τ WTA 21 (left) and τ WTA 32 (right) are compared with different MC predictions for large-radius jets from top (top row) and W (bottom row) selections.The distributions of τ WTA 21 (bottom left) and τ WTA 32

Table 1 :
Main features of the Monte Carlo models used to simulate signal and background samples, and to produce predictions to be compared with data.The nominal samples listed are used for comparisons with corrected data at particle level as well.For convenience, MG5_aMC@NLO is referred to as MG5_aMC in Figures3-9.

Table 3 :
Contributions from background processes which are subtracted in the top and W selections.The uncertainties are statistical only.
13 TeV with the ATLAS detector, Eur.Phys.J. C 78 (2018) 565, arXiv: 1804.10823[hep-ex].Also at California State University, East Bay; United States of America.c Also at Centre for High Performance Computing, CSIR Campus, Rosebank, Cape Town; South Africa.d Also at CERN, Geneva; Switzerland.e Also at CPPM, Aix-Marseille Université, CNRS/IN2P3, Marseille; France.f Also at Département de Physique Nucléaire et Corpusculaire, Université de Genève, Genève; Switzerland.g Also at Departament de Fisica de la Universitat Autonoma de Barcelona, Barcelona; Spain.h Also at Departamento de Física, Instituto Superior Técnico, Universidade de Lisboa, Lisboa; Portugal.i Also at Department of Applied Physics and Astronomy, University of Sharjah, Sharjah; United Arab Emirates.j Also at Department of Financial and Management Engineering, University of the Aegean, Chios; Greece.k Also at Department of Physics and Astronomy, University of Louisville, Louisville, KY; United States of America.l Also at Department of Physics and Astronomy, University of Sheffield, Sheffield; United Kingdom.m Also at Department of Physics, California State University, Fresno CA; United States of America.n Also at Department of Physics, California State University, Sacramento CA; United States of America.o Also at Department of Physics, King's College London, London; United Kingdom.p Also at Department of Physics, St. Petersburg State Polytechnical University, St. Petersburg; Russia.q Also at Department of Physics, Stanford University, Stanford CA; United States of America.r Also at Department of Physics, University of Fribourg, Fribourg; Switzerland.s Also at Department of Physics, University of Michigan, Ann Arbor MI; United States of America.t Also at Giresun University, Faculty of Engineering, Giresun; Turkey.u Also at Graduate School of Science, Osaka University, Osaka; Japan.v Also at Hellenic Open University, Patras; Greece.w Also at Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest; Romania.x Also at Institucio Catalana de Recerca i Estudis Avancats, ICREA, Barcelona; Spain.y Also at Institut für Experimentalphysik, Universität Hamburg, Hamburg; Germany.z Also at Institute for Mathematics, Astrophysics and Particle Physics, Radboud University Nijmegen/Nikhef, Nijmegen; Netherlands.aa Also at Institute for Nuclear Research and Nuclear Energy (INRNE) of the Bulgarian Academy of Sciences, Sofia; Bulgaria.ab Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest; Hungary.ac Also at Institute of Particle Physics (IPP); Canada.ad Also at Institute of Physics, Academia Sinica, Taipei; Taiwan.ae Also at Institute of Physics, Azerbaijan Academy of Sciences, Baku; Azerbaijan. a f Also at Institute of Theoretical Physics, Ilia State University, Tbilisi; Georgia.ag Also at Instituto de Física Teórica de la Universidad Autónoma de Madrid; Spain.ah Also at Istanbul University, Dept. of Physics, Istanbul; Turkey.Also at Joint Institute for Nuclear Research, Dubna; Russia.a j Also at LAL, Université Paris-Sud, CNRS/IN2P3, Université Paris-Saclay, Orsay; France.ak Also at Louisiana Tech University, Ruston LA; United States of America.al Also at LPNHE, Sorbonne Université, Paris Diderot Sorbonne Paris Cité, CNRS/IN2P3, Paris; France.am Also at Manhattan College, New York NY; United States of America.an Also at Moscow Institute of Physics and Technology State University, Dolgoprudny; Russia.ao Also at National Research Nuclear University MEPhI, Moscow; Russia.ap Also at Physics Dept, University of South Africa, Pretoria; South Africa.aq Also at Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Freiburg; Germany.ar Also at School of Physics, Sun Yat-sen University, Guangzhou; China.as Also at The City College of New York, New York NY; United States of America.at Also at The Collaborative Innovation Center of Quantum Matter (CICQM), Beijing; China.au Also at Tomsk State University, Tomsk, and Moscow Institute of Physics and Technology State University, Dolgoprudny; Russia.av Also at TRIUMF, Vancouver BC; Canada.aw Also at Universita di Napoli Parthenope, Napoli; Italy.
b ai * Deceased