Top-quark mass measurement in the all-hadronic $t\bar{t}$ decay channel at $\sqrt{s}$ = 8 TeV with the ATLAS detector

The top-quark mass is measured in the all-hadronic top-antitop quark decay channel using proton-proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 8 TeV with the ATLAS detector at the CERN Large Hadron Collider. The data set used in the analysis corresponds to an integrated luminosity of 20.2 fb$^{-1}$. The large multi-jet background is modelled using a data-driven method. The top-quark mass is obtained from template fits to the ratio of the three-jet to the dijet mass. The three-jet mass is obtained from the three jets assigned to the top quark decay. From these three jets the dijet mass is obtained using the two jets assigned to the W boson decay. The top-quark mass is measured to be 173.72 $\pm$ 0.55 (stat.) $\pm$ 1.01 (syst.) GeV.


Introduction
Of all known fundamental particles, the top quark has the largest mass. Its existence was predicted in 1973 by Kobayashi and Maskawa [1], and it was not observed directly until 1995, by the CDF and D0 experiments at the Tevatron [2,3]. Since 2010, top quarks have also been observed at the Large Hadron Collider (LHC) [4] at CERN. Due to the higher centre-of-mass energy, top quark production at the LHC is an order of magnitude larger than at the Tevatron. The large data sets of top-antitop quark (tt) pairs allow many precision studies and measurements of top quark properties. The Yukawa coupling of the top quark is predicted to be close to unity [5,6], suggesting that it may play a special role in electroweak symmetry breaking. In the Standard Model (SM), the top quark dominantly contributes to the quantum corrections to the Higgs self coupling [7,8]. Precise measurements of the top-quark mass (m top ) are therefore very important in probing the stability of the vacuum [9, 10], and contribute to searches for signs of physics beyond the SM.
Today the most precise individual measurement of m top is in the single-lepton decay channel of top-antitop quark pairs, where one top quark decays into a b-quark, a charged lepton and a neutrino and the other top quark decays into a b-quark and two u/d/c/s-quarks, performed by the CMS Collaboration, yielding a value of m top = 172.35 ± 0.16 (stat.) ± 0.48 (syst.) GeV [11]. The most precise measurement of m top in the dileptonic tt decay channel, where each of the top quarks decays into a b-quark, a charged lepton and its neutrino, is from the ATLAS Collaboration, yielding a value of m top = 172.99 ± 0.41 (stat.) ± 0.74 (syst.) GeV [12]. Further m top results are available in Refs. [13][14][15].
The top-quark mass measurement in the all-hadronic tt channel takes advantage of the largest branching ratio (46%) among the possible top quark decay channels [16]. The all-hadronic channel involves six jets at leading order, two originating from b-quarks and four originating from the two W boson hadronic decays. It is a challenging measurement because of the large multi-jet background arising from various quantum chromodynamics (QCD) processes, which can exceed the tt production by several orders of magnitude. However, all-hadronic tt events profit from having no neutrinos among the decay products, so that all four-momenta can be measured directly. The multi-jet background for the all-hadronic tt channel, while large, leads to different systematic uncertainties than in the case of the single-and dileptonic tt channels. Thus, all-hadronic analyses offer an opportunity to cross-check top-quark mass measurements performed in the other channels. The most recent measurements of m top in the all-hadronic channel were performed by the CMS Collaboration with m top = 172.32 ± 0.25 (stat.) ± 0.59 (syst.) GeV [11], and the ATLAS Collaboration with m top = 175.1 ± 1.4 (stat.) ± 1.2 (syst.) GeV [17].
This paper presents a top-quark mass measurement in the tt all-hadronic channel using data collected by the ATLAS experiment in 2012. The m top measurement is obtained from template fits to the distribution of the ratio of three-jet to dijet masses (R 3/2 = m j j j /m j j ), similarly to a previous measurement at √ s = 7 TeV [17]. The three-jet mass is obtained from the three jets assigned to the top quark decay. From the selected three jets the dijet mass is obtained using the two jets assigned to the W boson decay. The jet assignment is accomplished by using a χ 2 fit to the tt system, so there are two values of R 3/2 measured in each event. The observable R 3/2 employed in this analysis achieves a partial cancellation of systematic effects common to the masses of the reconstructed top quark and associated W boson, notably the significant uncertainty on the jet energy scale. Data-driven techniques are used to estimate the contribution from multi-jet background events. Data events are divided into several disjoint regions using two uncorrelated observables. The region containing the largest relative fraction of tt events is labeled the signal region. The background is estimated from the other regions, which determine the shape of the background distribution in the signal region.
The paper is organised as follows. After a brief description of the ATLAS detector in Section 2, the data and Monte Carlo (MC) samples used in the analysis are described in Section 3. The analysis event selection is further detailed in Section 4. Section 5 describes the method used to select the candidate four-momenta that comprise the reconstructed tt system. The estimation of the multi-jet background is detailed in Section 6. The method used to measure the top-quark mass and its uncertainties are reported in Sections 7, 8, and 9. The results of the measurement are presented in Section 10, and the analysis is summarised in Section 11.

ATLAS detector
The ATLAS detector [18] is a multi-purpose particle physics experiment with a forward-backward symmetric cylindrical geometry and near 4π coverage in solid angle 1. The inner tracking detector (ID) covers the pseudorapidity range |η| < 2.5, and consists of a silicon pixel detector, a silicon microstrip detector, and, for |η| < 2.0, a transition radiation tracker. The ID is surrounded by a thin superconducting solenoid providing a 2 T magnetic field. A high-granularity lead/liquid-argon (LAr) sampling electromagnetic calorimeter covers the region |η| < 3.2. A steel/scintillator-tile calorimeter provides hadronic coverage in the range |η| < 1.7. LAr technology is also used for the hadronic calorimeters in the endcap region 1.5 < |η| < 3.2 and for electromagnetic and hadronic measurements in the forward region up to |η| = 4.9. The muon spectrometer surrounds the calorimeters. It consists of three large air-core superconducting toroid systems, precision tracking chambers providing accurate muon tracking for |η| < 2.7, and additional detectors for triggering in the region |η| < 2.4.

Data and Monte Carlo simulation
This analysis is performed using the proton-proton (pp) collision data set at a centre-of-mass energy of √ s = 8 TeV collected with the ATLAS detector at the LHC. The data correspond to an integrated luminosity of 20.2 fb −1 . Samples of simulated MC events are used to optimise the analysis, to study the detector response and the efficiency to reconstruct tt events, to build signal template distributions used for fitting the top-quark mass, and to estimate systematic uncertainties. Most of the MC samples used in the analysis are based on a full simulation of the ATLAS detector [19] obtained using GEANT4 [20]. Some of the systematic uncertainties are studied using alternative tt samples processed through a faster ATLAS simulation (AFII) using parameterised showers in the calorimeters [21]. Additional simulated pp collisions generated with P [22] are overlaid to model the effects of additional collisions in the same and nearby bunch crossings (pile-up). All simulated events are processed using the same reconstruction algorithms and analysis chain as used for the data.
The nominal tt simulation sample is generated using the next-to-leading-order (NLO) MC program POWHEG-BOX [23-25] with the NLO parton distribution function (PDF) set CT10 [26, 27], interfaced 1 The coordinate system used to describe the ATLAS detector is briefly summarised here. The nominal interaction point is defined as the origin of the coordinate system, while the beam direction defines the z-axis and the x-y plane is transverse to the beam direction. The positive x-axis is defined as pointing from the interaction point to the centre of the LHC ring and the positive y-axis is defined as pointing upwards. The azimuthal angle φ is measured around the beam axis, and the polar angle θ is the angle from the beam axis. The pseudorapidity is defined as η = − ln tan(θ/2). The transverse momentum p T , the transverse energy E T , and the missing transverse momentum (E miss T ) are defined in the x-y plane unless stated otherwise. The distance ∆R in the η-φ angle space is defined as ∆R = (∆η) 2 + (∆φ) 2 . to P 6.427 [28] with a set of tuned parameters called the P 2012 tune [29] for parton shower, fragmentation and underlying-event modelling. For the construction of the signal templates, MC events are generated at five different assumed values of m top , between 167.5 and 177.5 GeV, in steps of 2.5 GeV. The full simulation of the ATLAS detector sample at 172.5 GeV has the largest number of generated events, and is used as the nominal signal sample. The h damp parameter [30], which regulates the high-p T radiation in POWHEG-BOX, is set to the same m top value as used in each of the generated POWHEG-BOX samples. All the simulated samples used to estimate systematic uncertainties are further described in Section 9.
All MC samples are normalised using the predicted top-antitop quark pair cross-section (σ tt ) at √ s = 8 TeV. For m top = 172.5 GeV, the next-to-next-to-leading-order cross-section of σ tt = 253 +13 −15 pb is calculated using the program T ++2.0 [31], which includes re-summation of next-to-next-to-leading logarithmic soft gluon terms.

Event selection
Events in this analysis are selected by a trigger that requires at least five jets with p T > 55 GeV. Only events with a well-reconstructed primary vertex formed by at least five tracks with p T > 400 MeV are considered for the analysis. Events with isolated electrons (muons) with E T > 25 GeV (p T > 20 GeV) and reconstructed in the central region of the detector within |η| < 2.5 are rejected. Both lepton types are identified using the tight working points as specified in Refs. [32,33]. Jets ( j) are reconstructed using the anti-k t algorithm with radius parameter R = 0.4 [34] employing topological clusters [35] in the calorimeter. These jets are calibrated to the hadronic energy scale as described in Refs. [36][37][38]. The four-vector of the highest-energy muon (µ) from among those matched within ∆R( j, µ) < 0.3 to a reconstructed jet, is added to the reconstructed jet four-vector. This is done to compensate for the energy losses in the calorimeter arising from semimuonic quark decays. In simulation this correction slightly improves both the jet energy response and resolution across the full range of jet energies.
To ensure that the selected events are in the plateau region of the trigger efficiency curve where the trigger efficiency in data is greater than 90%, at least five of the reconstructed central jets (within |η| < 2.5) are required to have p T > 60 GeV. Any additional jet is required to have p T > 25 GeV and |η| < 2.5. All selected jets in an event must be isolated; any pairing of two jets ( j i and j k ) reconstructed with the above criteria are required to not overlap within ∆ R( j i , j k ) < 0.6. Events with jets failing this isolation requirement are rejected.
Events containing neutrinos are removed by requiring E miss T < 60 GeV. The E miss T in an event is computed as the sum of a number of different terms [39,40]. Muons, electrons and jets are accounted for using the appropriate calibrations for each object. For each term considered, the missing transverse momentum is calculated as the negative sum of the calibrated reconstructed objects, projected onto the x and y directions.
For the final selection, events are kept if at least two of the six leading transverse momentum jets are identified as originating from a b-quark. Such jets are said to be b-tagged. A neural network trained on decay vertex properties [41] is used to identify these b-tagged jets. Because of the large number of c-quarks originating from the W boson decays in this analysis (on average one c-quark per tt event) a b-tagger trained to reject u/d/s-jets but also a large fraction of c-jets is used. Events with fewer than two b-tagged jets are used for the background estimate described in Section 6. The chosen working point for
In each event the two jets with leading b-tag weights (b i and b j ) are required to satisfy ∆φ(b i , b j ) > 1.5. The quantities b i and b j represent here the vectorial transverse momentum of a b-jet: b T,i and b T, j . This ∆φ cut is very powerful in rejecting combinatorial background events; most of these are true tt events where the incorrect jets are associated with the top quark. Finally, a cut is applied based on the azimuthal angle between b-jets and their associated W boson candidate: the average of the two angular separations for each event is required to satisfy ∆φ(b, W) < 2. Here the b, and the W are the vectorial transverse momentum of a b-jet and a W boson: b T and W T , identified by means of the three-jet combination that best fits the tt event hypothesis described in Section 5. This ∆φ cut rejects a large fraction of events from the multi-jet and combinatorial backgrounds, as well as events from non-all-hadronic tt decays. Events failing this final selection cut are, however retained for the purpose of modelling the multi-jet background, as detailed in Section 6. Table 1 summarises the yields obtained after each of the individual selection cuts. The χ 2 cut listed in Table 1 is described in Section 5. The number of b-tagged jets (N b tag ) and ∆φ(b, W) are the two observables used for the data-driven multi-jet background estimation, further detailed in Section 6.

tt reconstruction
In each event the tt final state is reconstructed using all the jets from the all-hadronic tt decay chain: To determine the top-quark mass in each tt event, a minimum-χ 2 approach is adopted, with the χ 2 defined as: (1) Here, two of the reconstructed jets are associated with the bottom-type quarks produced directly from the top quark and antitop quark decays (b 1

Multi-jet background estimation
The available MC generators for multi-jet production include only leading-order theory calculations for final states with up to six partons. Therefore, the dominant multi-jet background in this analysis is determined directly from the data. Two largely uncorrelated variables are used to divide the data events into four different regions, such that the background is determined in the control regions and extrapolated to the signal region. The two chosen observables are the N b tag in an event, and the ∆φ(b, W) variable, both described in Section 4. These have a correlation measured in data of ρ = −0.038. The value of N b tag in each event is determined from the leading six jets ordered by p T .
The four regions, labelled ABCD, are identified by defining two bins in the number of b-tagged jets, Table 2. The R 3/2 distributions are studied for each of the defined regions. Region D represents the signal region (SR), and contains the largest fraction of tt events (34.05%). Regions A, B, and C are the control regions (CR), and are dominated by multi-jet background events. Table 2 summarises the expected fractions of signal events in each of the four regions. Each signal fraction is estimated by comparing the total predicted number of signal events from tt simulation to the number of observed data events in each region.
To obtain an unbiased estimate of the number of background events in each considered CR, the signal contamination is removed using simulated tt events with m top = 172.5 GeV. The method validation and the template closure described in Section 8 show that the m top dependence of this signal subtraction is ABCD region and definition Estimated signal fraction Table 2: Definitions and signal fractions for each of the four regions used to estimate the multi-jet background. Region D is the signal region. The signal fraction with statistical uncertainty is estimated by comparing the total predicted number of signal events from tt simulation to the number of observed data events in each region.
significantly smaller than other uncertainties on the method, and is ignored. The estimated background in a given bin i of R 3/2 for SR D (N SR D background,i ) is given by: The background in a given bin i of the R 3/2 spectrum of CR B (N CR B background,i ) is estimated after subtraction of the signal contamination and after scaling by the ratio of the number events in control regions C (N CR C background ) and A (N CR A background ), also after signal removal. The signal contamination present in CR C comes from improperly reconstructed tt events which form a smoothly varying distribution in R 3/2 . This signal contribution in CR C is not relevant in the analysis, as this region only affects the normalisation of the distribution obtained for the multi-jet background, which is not used in the fit for m top described in Section 7. Figure 1 shows the distributions of the masses of the W boson (m j j ) and top quark (m j j j ) after applying the event selection, the χ 2 approach defined in Eq. (1), and using the data-driven multi-jet background method. In the figure, the reconstruction using MC events is said to be correct for one (or both) top quark(s) if each of the three jets ( j) selected by the reconstruction algorithm matches to each of the three quarks (q) within a ∆R( j, q) < 0.3, modulo the interchange of the two jets assigned to the hadronically decaying W boson. If at least one of the jets selected by the algorithm is not one of the three jets matched to the quarks, the top quark reconstruction is classified as incorrect. Finally, events where at least one quark is not matched uniquely to a reconstructed jet are classified as non-matched. The R 3/2 distribution obtained after using the data-driven multi-jet background estimation methods to determine the shape and normalisation is shown in Figure 2. In general, good agreement between data and prediction is observed in all the distributions.

Top-quark mass determination
To extract a measurement of the top-quark mass, a template method with a binned minimum-χ 2 approach is employed. For each tt event, two R 3/2 values are obtained, one for each top-quark mass measurement. To properly correct for the correlation between the two R 3/2 values in each event, the statistical uncertainty of m top returned from the final χ 2 fit described later in this section is scaled up by a factor √ 1 + ρ = 1.26, where ρ = 0.59 is the correlation factor as obtained from data. Signal and background templates binned  in R 3/2 are created using the simulated tt events described in Section 3, and the data-driven background distribution.
The top quark contribution is parameterised by a probability distribution function (pdf) which is the sum of a Novosibirsk function [43] and a Landau function [44]. These describe, respectively, the signal and the combinatorial background. As a first step, the R 3/2 distributions from the five tt simulation samples with differing m top are fitted separately to determine the six parameters for each template mass. The MC simulation shows that each of these parameters depends linearly on the input m top . In the next step, the parameters are fitted to obtain the offsets and slopes of the linear m top dependencies. These values are then used as inputs to a combined, simultaneous fit to all five R 3/2 distributions. In total 12 parameters are derived by the combined fit to determine the pdf. Figure 3 shows the R 3/2 distributions obtained using the tt MC samples based on the full simulation of the ATLAS detector and generated at three top-quark mass points: 167.5, 172.5, and 175 GeV. Results from the combined, simultaneous fit to all five R 3/2 distributions are superimposed. Shown are the functions describing the signal and combinatorial background, respectively, and their sum. The Novosibirsk mean and width parameters offer the strongest sensitivity to m top . Template distributions obtained simultaneously for three separate input values of m top (167.5, 172.5, and 177.5 GeV), highlighting the R 3/2 shape to m top , are shown in Figure 4.
The multi-jet background template distribution obtained from the output of the data-driven method described in Section 6 can be parameterised in a similar fashion. In this case the sum of a Gaussian function and a Landau function was found to be a suitable choice for the functional form. The background pdf requires five parameters.
As a final step in the parameterisation, in order to take properly into account the uncertainties and the correlations between the various signal and background shape parameters, a more generalised version of the χ 2 function is used. The final χ 2 fit, which uses matrix algebra to include non-diagonal covariance matrices, has the form: Here m top and F bkgd are the two parameters which are left to float. The shape of the fitted multi-jet background parameterisation is assumed to be independent of m top while the normalisation, controlled by a background fraction parameter, F bkgd , is obtained by fitting the data distribution. The F bkgd is defined within the fit range of the R 3/2 distribution: 1.5 ≤ R 3/2 < 3.5. The term n i in Eq.(3) corresponds to the number of entries in bin i in the R 3/2 data distribution, whereas µ i corresponds to the estimated total number of signal and background entries. The term V data is the N bin × N bin diagonal data covariance matrix with V ik = δ ik n i , which accounts for the statistical uncertainty in each bin i. Similarly, V signal and V bkgd are N bin × N bin non-diagonal covariance matrices which account for the signal and background shape parameterisation uncertainties and their correlations. In the R 3/2 distribution which has a total number of data entries N d , and a given bin width w bin , the number of estimated entries in bin i, µ i , is given by: where P S and P B are the probability density functions for the signal and background, respectively.
where δm top is the statistical uncertainty of the m top parameter obtained from the fit of the pseudo experiment. The correction that takes into account the correlation between two R 3/2 values in each event, described in Section 7, is not applied here, as the values of R 3/2 drawn for pseudo experiments are uncorrelated. The pull distribution for an unbiased measurement has a mean of zero and a standard deviation of unity. A fitted pull mean value of 0.19 ± 0.13 and a fitted pull width of 0.98 ± 0.01 are obtained, which shows that the uncertainty determination is unbiased.

Systematic uncertainties
This section outlines the various sources of systematic uncertainty in m top which are summarised in Table 3. All sources are treated as uncorrelated. Individual contributions are symmetrised and the total uncertainty is taken as the sum in quadrature of all contributions.
The majority of the systematic uncertainties are assessed by varying the tt MC sample to reflect the uncertainty from each of these sources. Pseudo experiments are constructed from the varied sample, which are then passed through the analysis chain; the change in the result relative to that obtained from the nominal MC sample is evaluated. Exceptions to this are described in the following subsections. To facilitate a combination with other results, each systematic uncertainty is assigned a statistical uncertainty, taking into account the statistical correlation of the considered samples. Following Ref. [46], the systematic uncertainties listed in Table 3  In what follows, each source of systematic uncertainty is briefly described. These are broken down into three categories. The first category, theory and modelling uncertainties, is associated with the simulation of the signal events. The second set of uncertainties is related to the analysis method. These involve uncertaintie due to the way that the analysis was performed, including the choice of a template method, the background modelling, and the final m top extraction procedure. Finally a third category, calibrationand detector-related uncertainties involves uncertainties coming from the standard calibrations of physics objects.  Table 3: Summary of all sources of statistical and systematic uncertainties in the measured values of the top-quark mass. Totals are evaluated by means of a sum in quadrature and assuming that all contributions are uncorrelated. The uncertainties are subdivided into three categories: theory and modelling uncertainties, method-related uncertainties, and calibration-or detector-related uncertainties, as described in the text. Adjacent to each of the quoted systematic variations in m top is its associated statistical uncertainty. The ABCDEF method is further described in Section 9.2 and in Ref.

Source of uncertainty
[17]. The quoted statistical uncertainty is corrected for the correlation between the two R 3/2 measurements of each event.  [50,51]. To determine the final systematic uncertainty, the quantity m top is calculated for each of the three sets, where m top is the measured value from the central reference sample of the corresponding PDF set. Half of the difference between the largest and the smallest of these m top values is quoted as the symmetrised final uncertainty, and is 0.04 GeV.

Initial-state and final-state radiation:
Varying the amount of initial-and final-state radiation (ISR and FSR) can have an impact on the number of reconstructed jets, which in turn can affect the overall measurement of the top-quark mass. In order to quantify the sensitivity of the measurement to ISR/FSR, two alternative POWHEG-BOX plus . Simulated tt AFII events were produced with an increased number of semi-hard MPI (P 2012 H ) in order to assess the potential impact on the final measurement. The absolute difference between the results of these pseudo experiments and the one using the nominal simulated AFII sample is quoted as the systematic uncertainty and is 0.13 GeV.
Colour reconnection: When simulating AFII signal events using P 6.427 [28] for the parton shower and hadronisation modelling, there is a tunable parameter associated with the colour reconnection strength due to the colour flow along parton lines in the strong-interaction hard-scattering process. An alternative AFII tt sample uses the P P 2012 CR tune [29], which corresponds to reduced colour reconnection strength. The absolute difference of 0.12 GeV between the results of these pseudo experiments and the average m top value obtained using the nominal P 6.427 tt events is quoted as the systematic uncertainty.

Signal and background parameterisation:
To extract m top as described in Section 7, the uncertainties in the shape parameters of the R 3/2 observable for the signal contributions are included in the N bin × N bin covariance matrices which enter into the χ 2 minimisation used to extract m top (see Eq. (3)). Omitting these contributions would yield a simplified definition of the χ 2 variable: which can be recognised as the standard definition of the χ 2 variable for a least-squares fit assuming only a diagonal covariance matrix. The fit to the data distribution is repeated using this simplified definition of the χ 2 variable. This results in a slightly modified returned value of the m top parameter and a smaller statistical uncertainty. The difference in quadrature of 0.09 GeV between the final statistical uncertainty returned from the original minimisation and this modified value is quoted as the uncertainty in the signal and background parameterisation.
Inclusion of non-all-hadronic tt background: A number of event selection requirements, such as the lepton veto and the requirement that E miss T < 60 GeV, result in a large suppression of background contributions arising from non-all-hadronic tt events. The estimated fractional contribution from such events in the final signal region is below 3%, and is not considered in the nominal case. Pseudo experiments are performed by drawing events from the nominal signal distribution but from a modified background, now consisting of QCD events and tt events with at least one leptonic W boson decay. The absolute difference of 0.06 GeV between the average m top value obtained in this way and that from the nominal case is quoted as a systematic uncertainty.

Variation in the number of control regions:
A variation of the background estimation procedure is considered in which six distinct regions, rather than four, are defined to estimate the multi-jet background. This is done by allowing three different values of N b tag : 0, 1, or ≥ 2. Events can then be separated into the six differing regions as in the nominal analysis. As in the nominal case the number of b-tagged jets in an event considers only the leading six jets, ordered by p T . The values of the second ABCD variable, ∆φ(b, W) , are unchanged from the nominal case. One reason for considering this alternative is that the inclusion of a larger number of control regions could potentially provide sensitivity to different physics processes. Additionally, the systematic uncertainty contribution arising from uncertainties in the b-tagging scale factors could differ between these methods.
With a total of six regions shown in Table 4 the background estimation technique remains similar to that using four regions.
The final SR is labelled F. The new region D, together with region B, is now used to predict the shape of the multi-jet background in SR F, whereas CR A, C, and E set the multi-jet background normalisation [17]. Pseudo experiments are performed by drawing background events from the modified multi-jet distribution in the final signal region. The absolute difference of 0.16 GeV between this and the nominal case is quoted as the systematic uncertainty.

Calibration-and detector-related uncertainties
Trigger efficiency: The trigger efficiency obtained using simulated signal events [23-25, 28, 29] is compared to an equivalent distribution obtained using data, which results in a small observed discrepancy. The data here are expected to consist primarily of multi-jet events. It is expected that some true kinematic differences give rise to the difference observed between the data and MC trigger efficiencies. In order to obtain a conservative uncertainty, it is assumed that the difference represents a mis-modelling of the data by the trigger simulation. The simulated events are assigned a p T -dependent trigger efficiency correction such that the corrected MC and data trigger efficiencies agree. Pseudo experiments are performed by drawing signal events from the modified R 3/2 distribution with the trigger SFs applied, and the 0.08 GeV absolute difference from the nominal case is quoted as a conservative uncertainty on m top due to the trigger efficiency.

Pile-up reweighting scale:
The distribution of the average number of interactions per bunch crossing, denoted by µ , is known to differ between data and simulation. Simulated events are reweighted so that µ matches the value observed in data. In order to assess the impact on the final result, pseudo experiments are performed in which the reweighting scale is shifted up and down according to its uncertainty, and the fit procedure is repeated. A negligible maximum change of 0.01 GeV in m top is found as the symmetrised up/down uncertainty.

Lepton and E miss T soft-term calibrations:
Uncertainties in the calibration scales and in the resolutions of the lepton (e/µ) four-vector objects [33,53,54] can potentially lead to small differences in the event selection or the jet-quark assignment in the top reconstruction algorithm. Similarly, small uncertainties in m top can be expected due to the uncertainties in the scale and resolution of the E miss T soft term [39,40]. The E miss T soft term is varied according to these uncertainties and pseudo experiments are performed with the modified MC events. In the case of the muon-related uncertainties, Gaussian smearing is performed to assess the impact on the final result. The maximum absolute deviation from the reference m top value is taken as the uncertainty in each case, and these are added in quadrature to obtain a single value of 0.02 GeV for all lepton-and E miss T -related scale and resolution uncertainties.

Flavour-tagging efficiencies:
In the validation of the flavour-tagging algorithms, the differences between tagging efficiencies and mis-tag rates evaluated in data and simulation are removed by applying scale factor (SF) weights to the simulated events. The uncertainties in the flavour-tagging SFs are calculated separately for the b-tagging SFs, the c/τ-tagging SFs, and the overall mis-tag SFs [42]. The uncertainties in the flavour-tagging SFs are split into various components. The full covariance matrix between the various bins of jet transverse momentum is built and decomposed into eigenvectors. Each eigenvector corresponds to an independent source of uncertainty, each with an upward and a downward fluctuation, and the resulting total systematic uncertainty is 0.10 GeV.

Jet energy scale:
The different contributions to the total JES uncertainty are estimated individually as described in Ref. [36]. For each component the resulting differences from the up and down variations, corresponding to one-standard-deviation relative to the nominal JES, are quoted separately. The total uncertainty for each contribution is taken as half of the absolute difference between the up and down variation. In case both the up and down variations result in a change in the parameter in the same direction, the largest absolute difference (either from the up or down variation) is taken as the symmetrised uncertainty. The total JES uncertainty is the sum in quadrature of all subcontributions, and is 0.60 GeV. This includes all but the b-jet energy scale contribution, which is quoted separately and discussed below. Jet energy resolution: An eigenvector decomposition strategy similar to that followed for the JES and the flavour-tagging systematic uncertainties is used for the determination of jet energy resolution (JER) systematic uncertainties [56]. The final quoted JER systematic uncertainty is 0.10 GeV.

Jet reconstruction efficiency:
A small difference between the jet reconstruction efficiencies measured in data and simulation was observed [37], and as this difference can affect the final measured m top value, a set of pseudo experiments are performed in which jets from simulated events are removed at random. The frequency of this is chosen such that the modified jet reconstruction efficiency in simulation matches the value measured in data. The analysis is repeated with this change and no significant difference is observed.

Measurement of m top
After applying the method described in Section 7 the top-quark mass is measured to be: The statistical error quoted in Eq. (7) is corrected for the correlation between the two R 3/2 measurements of each event, as discussed in Section 7. The systematic uncertainty quoted above is the sum in quadrature of all the systematic uncertainties described in Section 9 and summarised in Table 3. Figure 6 shows the

Conclusion
From the analysis of 20.2 fb −1 of data recorded with the ATLAS detector at the LHC at a pp centre-of-mass energy of 8 TeV, the top-quark mass has been measured in the all-hadronic decay channel of top-antitop quark pairs to be m top = 173.72 ± 0.55 (stat.) ± 1.01 (syst.) GeV.
This measurement is obtained from template fits to the R 3/2 observable, which is chosen due to its reduced dependence on the jet energy scale uncertainty.