Determination of the strong coupling constant α s from transverse energy_energy correlations in multijet events at √s̅ = 8 TeV using the ATLAS detector

Measurements of transverse energy–energy correlations and their associated asymmetries in multi-jet events using the ATLAS detector at the LHC are presented. The data used correspond to √ s = 8 TeV proton–proton collisions with an integrated luminosity of 20.2 fb − 1 . The results are presented in bins of the scalar sum of the transverse momenta of the two leading jets, unfolded to the particle level and compared to the predictions from Monte Carlo simulations. A comparison with next-to-leading-order per-turbative QCD is also performed, showing excellent agreement within the uncertainties. From this comparison, the value of the strong coupling constant is extracted for different energy regimes, thus testing the running of α s (μ) predicted in QCD up to scales over 1 TeV. A global ﬁt to the transverse energy–energy correlation distributions yields α s ( m Z ) = 0 . 1162 ± 0 . 0011 (exp.) + 0 . 0084 − 0 . 0070 (theo.), while a global ﬁt to the asymmetry distributions yields a value of α s ( m Z ) = 0 . 1196 ± 0 . 0013 (exp.) + 0 . 0075 − 0 . 0045 (theo.).


Introduction
Experimental studies of the energy dependence of event shape variables have proved very useful in precision tests of quantum chromodynamics (QCD).Event shape variables have been measured in e + e − experiments from PETRA-PEP [1][2][3] to LEP-SLC [4][5][6][7] energies, at the ep collider HERA [8][9][10][11][12] as well as in hadron colliders from Tevatron [13] to LHC energies [14,15].Most event shape variables are based on the determination of the thrust's principal axis [16] or the sphericity tensor [17].A notable exception is given by the energy-energy correlations (EEC), originally proposed by Basham et al. [18], and measurements [19][20][21][22][23][24][25][26][27][28][29][30][31] of these have significantly improved the precision tests of perturbative QCD (pQCD).The EEC is defined as the energy-weighted angular distribution of hadron pairs produced in e + e − annihilation and, by construction, the EEC as well as its associated asymmetry (AEEC) are infrared safe.The second-order corrections to these functions were found to be significantly smaller [32][33][34][35] than for other event shape variables such as thrust.The transverse energy-energy correlation (TEEC) and its associated asymmetry (ATEEC) were proposed as the appropriate generalisation to hadron colliders in Ref. [36], where leading-order (LO) predictions were also presented.As a jet-based quantity, it makes use of the jet transverse energy E T = E sin θ since the energy alone is not Lorentz-invariant under longitudinal boosts along the beam direction.Here θ refers to the polar angle of the jet axis, while E is the jet energy. 1 The next-to-leading-order (NLO) corrections were obtained recently [37] by using the NLOJET++ program [38,39].They are found to be of moderate size so that the TEEC and ATEEC functions are well suited for precision tests of QCD, including a precise determination of the strong coupling constant α s .The TEEC is defined as [40] 1 σ dΣ d cos φ ≡ 1 σ i j dσ dx Ti dx T j d cos φ x Ti x T j dx Ti dx T j = 1 where the last expression is valid for a sample of N hard-scattering multi-jet events, labelled by the index A, and the indices i and j run over all jets in a given event.Here, x Ti is the fraction of transverse energy of jet i with respect to the total transverse energy, i.e. x Ti = E Ti / k E Tk , φ i j is the angle in the transverse plane between jet i and jet j and δ(x) is the Dirac delta function, which ensures φ = φ i j .
The associated asymmetry ATEEC is then defined as the difference between the forward (cos φ > 0) and the backward (cos φ < 0) parts of the TEEC, i.e.
1 σ Recently, the ATLAS Collaboration presented a measurement of the TEEC and ATEEC [41], where these observables were used for a determination of the strong coupling constant α s (m Z ) at an energy regime of Q = 305 GeV.This paper extends the previous measurement to higher energy scales up to values close to 1 TeV.The analysis consists in the measurement of the TEEC and ATEEC distributions in different energy regimes, determining α s (m Z ) in each of them, and using these determinations to test the running of α s predicted by the QCD β-function.Precise knowledge of the running of α s is not only important as a precision test of QCD at large scales but also as a test for new physics, as the existence of new coloured fermions would imply modifications to the β-function [42,43]. 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 upward.Cylindrical coordinates (r, ϕ) are used in the transverse plane, ϕ being the azimuthal angle around the beam pipe.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).

ATLAS detector
The ATLAS detector [44] is a multipurpose particle physics detector with a forward-backward symmetric cylindrical geometry and a solid angle coverage of almost 4π.The inner tracking system covers the pseudorapidity range |η| < 2.5.It consists of a silicon pixel detector, a silicon microstrip detector and, for |η| < 2.0, a transition radiation tracker.It is surrounded by a thin superconducting solenoid providing a 2 T magnetic field along the beam direction.A high-granularity liquid-argon sampling electromagnetic calorimeter covers the region |η| < 3.2.A steel/scintillator tile hadronic calorimeter provides coverage in the range |η| < 1.7.The endcap and forward regions, spanning 1.5 < |η| < 4.9, are instrumented with liquid-argon calorimeters for electromagnetic and hadronic measurements.The muon spectrometer surrounds the calorimeters.It consists of three large air-core superconducting toroid systems and separate trigger and high-precision tracking chambers providing accurate muon tracking for |η| < 2.7.
The trigger system [45] has three consecutive levels: level 1, level 2 and the event filter.The level 1 triggers are hardware-based and use coarse detector information to identify regions of interest, whereas the level 2 triggers are software-based and perform a fast online data reconstruction.Finally, the event filter uses reconstruction algorithms similar to the offline versions with the full detector granularity.

Monte Carlo simulation
Multi-jet production in pp collisions is described by the convolution of the production cross-sections for parton-parton scattering with the parton distribution functions (PDFs).Monte Carlo (MC) event generators differ in the approximations used to calculate the underlying short-distance QCD processes, in the way parton showers are built to take into account higher-order effects and in the fragmentation scheme responsible for long-distance effects.Pythia and Herwig++ event generators were used for the description of multi-jet production in pp collisions.These event generators differ in the modelling of the parton shower, hadronisation and underlying event.Pythia uses p T -ordered parton showers, in which the p T of the emitted parton is decreased in each step, while for the angle-ordered parton showers in Herwig++, the relevant scale is related to the angle between the emitted and the incoming parton.The generated events were processed with the ATLAS full detector simulation [46] based on Geant4 [47].
The baseline MC samples were generated using Pythia 8.160 [48] with the matrix elements for the underlying 2 → 2 processes calculated at LO using the CT10 LO PDFs [49] and matched to p T -ordered parton showers.A set of tuned parameters called the AU2CT10 tune [50] was used to model the underlying event (UE).The hadronisation follows the Lund string model [51].
A different set of samples were generated with Herwig++ 2.5.2 [52], using the LO CTEQ6L1 PDFs [53] and the CTEQ6L1-UE-EE-3 tune for the underlying event [54].Herwig++ uses angle-ordered parton showers, a cluster hadronisation scheme and the underlying-event parameterisation is given by Jimmy [55].Additional samples are generated using Sherpa 1.4.5 [56], which calculates matrix elements for 2 → N processes at LO, which are then convolved with the CT10 LO PDFs, and uses the CKKW [57] method for the parton shower matching.These samples were generated with up to three hard-scattering partons in the final state.
In order to compensate for the steeply falling p T spectrum, MC samples are generated in seven intervals of the leading-jet transverse momentum.Each of these samples contain of the order of 6 × 10 6 events for Pythia8 and 1.4 × 10 6 events for Herwig++ and Sherpa.All MC simulated samples described above are subject to a reweighting algorithm in order to match the average number of pp interactions per bunch-crossing observed in the data.The average number of interactions per bunch-crossing amounts to µ = 20.4 in data, and to µ = 22.0 in the MC simulation.

Data sample and jet calibration
The data used were recorded in 2012 at √ s = 8 TeV and collected using a single-jet trigger.It requires at least one jet, reconstructed with the anti-k t algorithm [58] with radius parameter R = 0.4 as implemented in FastJet [59].The jet transverse energy measured by the trigger system is required to be greater than 360 GeV at the trigger level.This trigger is fully efficient for values of the scalar sum of the calibrated transverse momenta of the two leading jets, p T1 + p T2 , denoted hereafter by H T2 , above 730 GeV.This is the lowest unprescaled trigger for the 2012 data-taking period, and the integrated luminosity of the full data sample is 20.2 fb −1 .Events are required to have at least one vertex, with two or more associated tracks with transverse momentum p T > 400 MeV.The vertex maximising p 2 T , where the sum is performed over tracks, is chosen as the primary vertex.In the analysis, jets are reconstructed with the same algorithm as used in the trigger, the anti-k t algorithm with radius parameter R = 0.4.The input objects to the jet algorithm are topological clusters of energy deposits in the calorimeters [60].The baseline calibration for these clusters corrects their energy using local hadronic calibration [61,62].The four-momentum of an uncalibrated jet is defined as the sum of the four-momenta of its constituent clusters, which are considered massless.Thus, the resulting jets are massive.However, the effect of this mass is marginal for jets in the kinematic range considered in this paper, as the difference between transverse energy and transverse momentum is at the per-mille level for these jets.The jet calibration procedure includes energy corrections for multiple pp interactions in the same or neighbouring bunch crossings, known as "pile-up", as well as angular corrections to ensure that the jet originates from the primary vertex.Effects due to energy losses in inactive material, shower leakage, the magnetic field, as well as inefficiencies in energy clustering and jet reconstruction, are taken into account.This is done using an MC-based correction, in bins of η and p T , derived from the relation of the reconstructed jet energy to the energy of the corresponding particle-level jet, not including muons or non-interacting particles.In a final step, an in situ calibration corrects for residual differences in the jet response between the MC simulation and the data using p T -balance techniques for dijet, γ+jet, Z+jet and multi-jet final states.The total jet energy scale (JES) uncertainty is given by a set of independent sources, correlated in p T .The uncertainty in the p T value of individual jets due to the JES increases from (1-4)% for |η| < 1.8 to 5% for 1.8 < |η| < 4.5 [63].The selected jets must fulfill p T > 100 GeV and |η| < 2.5.The two leading jets are further required to fulfil H T2 > 800 GeV.In addition, jets are required to satisfy quality criteria that reject beam-induced backgrounds (jet cleaning) [64].The number of selected events in data is 6.2 × 10 6 , with an average jet multiplicity N jet = 2.3.In order to study the dependence of the TEEC and ATEEC on the energy scale, and thus the running of the strong coupling, the data are further binned in H T2 .The binning is chosen as a compromise between reaching the highest available energy scales while keeping a sufficient statistical precision in the TEEC distributions, and thus in the determination of α s .Table 1 summarises this choice, as well as the number of events in each energy bin and the average value of the chosen scale Q = H T2 /2, obtained from detector-level data.

Results at the detector level
The data sample described in Section 4 is used to measure the TEEC and ATEEC functions.In order to study the kinematical dependence of such observables, and thus the running of the strong coupling with the energy scale involved in the hard process, the binning introduced in Table 1 is used.Figure 1 compares the TEEC and ATEEC distributions, measured in two of these bins, with the MC predictions from Pythia8, Herwig++ and Sherpa.
The TEEC distributions show two peaks in the regions close to the kinematical endpoints cos φ = ±1.The first one, at cos φ = −1 is due to the back-to-back configuration in two-jet events, which dominate the sample, while the second peak at cos φ = +1 is due to the self-correlations of one jet with itself.These self-correlations are included in Equation (1) and are necessary for the correct normalisation of the TEEC functions.The central regions of the TEEC distributions shown in Figure 1 are dominated by gluon radiation, which is decorrelated from the main event axis as predicted by QCD and measured in Refs.[65,66].Among the MC predictions considered here, Pythia8 and Sherpa are the ones which fit the data best, while Herwig++ shows significant discrepancies with the data.

Correction to particle level
In order to allow comparison with particle-level MC predictions, as well as NLO theoretical predictions, the detector-level distributions presented in Section 5 need to be corrected for detector effects.Particlelevel jets are reconstructed in the MC samples using the anti-k t algorithm with R = 0.4, applied to final-state particles with an average lifetime τ > 10 −11 s, including muons and neutrinos.The kinematical requirements for particle-level jets are the same as for the definition of TEEC / ATEEC at the detector level.
In the data, an unfolding procedure is used which relies on an iterative Bayesian unfolding method [67] as implemented in the RooUnfold program [68].The method makes use of a transfer matrix for each distribution, which takes into account any inefficiencies in the detector, as well as its finite resolution.The Pythia8 MC sample is used to determine the transfer matrices from the particle-level to detector-level TEEC distributions.Pairs of jets not entering the transfer matrices are accounted for using inefficiency correction factors.Figure 1: Detector-level distributions for the TEEC (top) and ATEEC functions (bottom) for the first and the last H T2 intervals chosen in this analysis, together with MC predictions from Pythia8, Herwig++ and Sherpa.The total uncertainty, including statistical and detector experimental sources, i.e. those not related to unfolding corrections, is also indicated using an error bar for the distributions and a green-shaded band for the ratios.The systematic uncertainties are discussed in Section 7.
The excellent azimuthal resolution of the ATLAS detector, together with the reduction of the energy scale and resolution effects by the weighting procedure involved in the definition of the TEEC function, are reflected in the fact that the transfer matrices have very small off-diagonal terms (smaller than 10%), leading to very small migrations between bins.The statistical uncertainty is propagated through the unfolding procedure by using pseudo-experiments.A set of 10 3 replicas is considered for each measured distribution by applying a Poisson-distributed fluctuation around the nominal measured distribution.Each of these replicas is unfolded using a fluctuated version of the transfer matrix, which produces the corresponding set of 10 3 replicas of the unfolded spectra.The statistical uncertainty is defined as the standard deviation of all replicas.

Systematic uncertainties
The dominant sources are those associated with the MC model used in the unfolding procedure and the JES uncertainty in the jet calibration procedure.
• Jet Energy Scale: The uncertainty in the jet calibration procedure [63] is propagated to the TEEC by varying each jet energy and transverse momentum by one standard deviation of each of the 67 nuisance parameters of the JES uncertainty, which depend on both the jet transverse momentum and pseudorapidity.The total JES uncertainty is evaluated as the sum in quadrature of all nuisance parameters, and amounts to 2%.• Jet Energy Resolution: The effect on the TEEC function of the jet energy resolution uncertainty [69] is estimated by smearing the energy and transverse momentum by a smearing factor depending on both p T and η.This amounts to approximately 1% in the TEEC distributions.
• Monte Carlo modelling: The modelling uncertainty is estimated by performing the unfolding procedure described in Section 6 with different MC approaches.The difference between the unfolded distributions using Pythia and Herwig++ defines the envelope of the uncertainty.This was crosschecked using the difference between Pythia and Sherpa, leading to similar results.This is the dominant experimental uncertainty for this measurement, being always below 5% for the TEEC distributions, and being larger for low H T2 .
• Unfolding: The mismodelling of the data made by the MC simulation is accounted for as an additional source of uncertainty.This is assessed by reweighting the transfer matrices so that the level of agreement between the detector-level projection and the data is enhanced.The modified detector-level distributions are then unfolded using the method described in Section 6.The difference between the modified particle-level distribution and the nominal one is then taken as the uncertainty.This uncertainty is smaller than 0.5% for the full cos φ range for all bins in H T2 .The impact of this uncertainty on the TEEC function is below 1%.
• Jet Angular Resolution: The uncertainty in the jet angular resolution is propagated to the TEEC measurements by smearing the azimuthal coordinate ϕ of each jet by 10% of the resolution in the MC simulation.This is motivated by the track-to-cluster matching studies done in Ref. [65].This impacts the TEEC measurement at the level of 0.5%.
• Jet cleaning: The modelling of the efficiency of the jet-cleaning cuts is considered as an additional source of experimental uncertainty.This is studied by tightening the jet cleaning-requirements in both data and MC simulation, and considering the double ratio between them.The differences are below 0.5%.
In order to mitigate statistical fluctuations, the resulting systematic uncertainties are smoothed using a Gaussian kernel algorithm.The impact of these systematic uncertainties is summarised in Figure 2, where the relative errors are shown for the TEEC and ATEEC distributions for each H T2 bin considered.

Experimental results
The results of the unfolding are compared with particle-level MC predictions, including the estimated systematic uncertainties.Figure 3 shows this comparison for the TEEC, while the ATEEC results are shown in Figure 4.The level of agreement seen here between data and MC simulation is similar to that at detector level.Pythia and Sherpa broadly describe the data, while the Herwig++ description is disfavoured.Figure 4: Particle-level distributions for the ATEEC functions in each of the H T2 intervals chosen in this analysis, together with MC predictions from Pythia8, Herwig++ and Sherpa.The total uncertainty, including statistical and other experimental sources is also indicated using an error bar for the distributions and a green-shaded band for the ratios.

Theoretical predictions
The theoretical predictions for the TEEC and ATEEC functions are calculated using perturbative QCD at NLO as implemented in NLOJET++ [38,39].Typically O(10 10 ) events are generated for the calculation.The partonic cross-sections, σ, are convolved with the NNLO PDF sets from MMHT 2014 [70], CT14 [71], NNPDF 3.0 [72] and HERAPDF 2.0 [73] using the LHAPDF6 package [74].The value of α s (m Z ) used in the partonic matrix-element calculation is chosen to be the same as that of the PDF.At leading order in α s , the TEEC function defined in Equation ( 1) can be expressed as where Σa 1 a 2 →b 1 b 2 b 3 is the partonic cross-section weighted by the fractions of transverse energy of the outgoing partons, x Ti x T j as in Equation (1); x i (i = 1, 2) are the fractional longitudinal momenta carried by the initial-state partons, f a 1 /p (x 1 ) and f a 2 /p (x 2 ) are the PDFs and ⊗ denotes a convolution over x 1 ,x 2 .
At O(α 4 s ), the numerator in Equation ( 2) entails calculations of the 2 → 3 partonic subprocesses at NLO accuracy, and the 2 → 4 partonic subprocesses at LO.In order to avoid the double collinear singularities appearing in the latter, the angular range is restricted to | cos φ| ≤ 0.92.This avoids calculating the twoloop virtual corrections to the 2 → 2 subprocesses.Thus, with the azimuthal angle cut, the denominator in Equation ( 2) includes the 2 → 2 and 2 → 3 subprocesses up to and including the O(α 3 s ) corrections.The nominal renormalisation and factorisation scales are defined as a function of the transverse momenta of the two leading jets as follows [75] This choice eases the comparison with the previous measurement at √ s = 7 TeV [41], where the renormalisation scale was the same.The relevant scale for the perturbative calculation is the renormalisation scale, as variations of the factorisation scale lead to small variations of the physical observable.The scale choice for the NLO pQCD templates used to extract α s as well as for the presentation of the measurement is not uniquely defined.The nominal scale choice, H T2 /2, used in this paper is based on previous publications [41,76].However, it should be noted that other scale choices, which explicitly take into account the kinematics of the third jet, are also viable options and can be considered in future measurements.The following comments are in order.The NLOJet++ calculations are performed in the limit of massless quarks.PDFs are based on the n f = 5 scheme.There is therefore a residual uncertainty due to the mass of the top quark.This is expected to be small since at LHC energies σ tt σ QCD .The correct treatment of top quark mass effects in the initial as well as in final state is not yet available.

Non-perturbative corrections
The pQCD predictions obtained using NLOJET++ are generated at the parton level only.In order to compare these predictions with the data, one needs to correct for non-perturbative (NP) effects, namely hadronisation and the underlying event.Here, doing this relies on bin-by-bin correction factors calculated as the ratio of the MC predictions for TEEC distributions with hadronisation and UE turned on to those with hadronisation and UE turned off.These factors, which are calculated using several MC models, are used to correct the pQCD prediction to the particle level by multiplying each bin of the theoretical distributions.Figure 5 shows the distributions of the factors for the TEEC as a function of cos φ and for two bins in the energy scale H T2 .They were calculated using several models, namely Pythia8 with the AU2 [77] and 4C tunes [78] and Herwig++ with the LHC-UE-EE-3-CTEQ6L1 and LHC-UE-EE-3-LOMOD tunes [54].From these four posibilities, Pythia8 with the AU2 tune is used for the nominal corrections.

Theoretical uncertainties
The theoretical uncertainties are divided into three classes: those corresponding to the renormalisation and factorisation scale variations, the ones corresponding to the PDF eigenvectors, and the ones for the non-perturbative corrections.• The theoretical uncertainty due to the choice of renormalisation and factorisation scales is defined as the envelope of all the variations of the TEEC and ATEEC distributions obtained by varying up and down the scales µ R , µ F by a factor of two, excluding those configurations in which both scales are varied in opposite directions.This is the dominant theoretical uncertainty in this measurement, which can reach 20% in the central region of the TEEC distributions.
• The parton distribution functions are varied following the set of eigenvectors / replicas provided by each PDF group [70][71][72][73].The propagation of the corresponding uncertainty to the TEEC and ATEEC is done following the recommendations for each particular set of distribution functions.The size of this uncertainty is around 1% for each TEEC bin.
• The uncertainty in the non-perturbative corrections is estimated as the envelope of all models used for the calculation of the correction factors in Figure 5.This uncertainty is around 1% for each of the TEEC bins considered in the NLO predictions, i.e. those with | cos φ| ≤ 0.92.
• The uncertainty due to α s is also considered for the comparison of the data with the theoretical predictions.This is estimated by varying α s by the uncertainty in its value for each PDF set, as indicated in Refs.[70][71][72][73].
The total theoretical uncertainty is obtained by adding these four theoretical uncertainties in quadrature.The total uncertainty can reach 20% for the central part of the TEEC, due to the large value of the scale uncertainty in this region.

Comparison of theoretical predictions and experimental results
The unfolded data obtained in Section 8 are compared to the pQCD predictions, once corrected for nonperturbative effects.Figures 6 and 7 show the ratios of the data to the theoretical predictions for the TEEC and ATEEC functions, respectively.The theoretical predictions were calculated, as a function of cos φ and for each of the H T2 bins considered, using the NNPDF 3.0 PDFs with α s (m Z ) = 0.1180.
-  -0.9 -0.From the comparisons in Figures 6 and 7, one can conclude that perturbative QCD correctly describes the data within the experimental and theoretical uncertainties.

Determination of α s and test of asymptotic freedom
From the comparisons made in the previous section, one can determine the strong coupling constant at the scale given by the pole mass of the Z boson, α s (m Z ), by considering the following χ 2 function where the theoretical predictions are varied according to In Equations ( 3) and (4), α s stands for α s (m Z ); x i is the value of the i-th point of the distribution as measured in data, while ∆x i is its statistical uncertainty.The statistical uncertainty in the theoretical predictions is also included as ∆ξ i , while σ (i) k is the relative value of the k-th source of systematic uncertainty in bin i.This technique takes into account the correlations between the different sources of systematic uncertainty discussed in Section 7 by introducing the nuisance parameters {λ k }, one for each source of experimental uncertainty.Thus, the minimum of the χ 2 function defined in Equation ( 3) is found in a 74-dimensional space, in which 73 correspond to nuisance parameters {λ i } and one to α s (m Z ).The method also requires an analytical expression for the dependence of the fitted observable on the strong coupling constant, which is given by ψ i (α s ) for bin i.For each PDF set, the corresponding α s (m Z ) variation range is considered and the theoretical prediction is obtained for each value of α s (m Z ).The functions ψ i (α s ) are then obtained by fitting the values of the TEEC (ATEEC) in each (H T2 , cos φ) bin to a second-order polynomial.For both the TEEC and ATEEC functions, the fits to extract α s (m Z ) are repeated separately for each H T2 interval, thus determining a value of α s (m Z ) for each energy bin.The theoretical uncertainties are determined by shifting the theory distributions by each of the uncertainties separately, recalculating the functions ψ i (α s ) and determining a new value of α s (m Z ).The uncertainty is determined by taking the difference between this value and the nominal one.Each of the obtained values of α s (m Z ) is then evolved to the corresponding measured scale using the NLO solution to the renormalisation group equation (RGE), given by where the coefficients β 0 and β 1 are given by and Λ is the QCD scale, determined in each case from the fitted value of α s (m Z ).Here, n f is the number of active flavours at the scale Q, i.e. the number of quarks with mass m < Q.Therefore, n f = 6 in the six bins considered in Table 1.When evolving α s (m Z ) to α s (Q), the proper transition rules for n f = 5 to n f = 6 are applied so that α s (Q) is a continuous function across quark thresholds.Finally, the results are combined by performing a global fit, where all bins are merged together.

Fits to individual TEEC functions
The values of α s (m Z ) obtained from fits to the TEEC function in each H T2 bin are summarised in Table 2.
The theoretical predictions used for this extraction use NNPDF 3.0 as the nominal PDF set.The values summarised in Table 2 are in good agreement with the 2016 world average value [79], as well as with previous measurements, in particular with previous extractions using LHC data [41,76,[80][81][82][83][84].
The values of the χ 2 indicate that agreement between the data and the theoretical predictions is good.The nuisance parameters for the TEEC fits are generally compatible with zero.One remarkable exception is the nuisance parameter associated to the modelling uncertainty, which deviates by half standard deviation with a very small error bar.This is an indication that these data can be used to further tune MC event generators which model multi-jet production.Figure 8 compares the data with the theoretical predictions after the fit, i.e.where the fitted values of α s (m Z ) and the nuisance parameters are already constrained.Table 3 shows the values of α s evolved from m Z to the corresponding scale Q using Equation (5).The appendix includes tables in which the values of α s (m Z ) obtained from the TEEC fits are extrapolated to different values of Q, given by the averages of kinematical quantities other than H T2 /2.

Global TEEC fit
The combination of the previous results is done by considering all the H T2 bins into a single, global fit.The result obtained using the NNPDF 3.0 PDF set has the largest PDF uncertainty and thus, in order to be conservative, it is the one quoted as the final value of α s (m Z ).
The impact of the correlations of the JES uncertainties on the result is studied by considering two additional correlation scenarios, one with stronger and one with weaker correlation assumptions [63].From the envelope of these results, an additional uncertainty of 0.0007 is assigned in order to cover this difference.
The results for α s (m Z ) are summarised in The results for α s from fits to the TEEC using different PDFs.The uncertainty referred to as NP is the one related to the non-perturbative corrections.The uncertainty labelled as 'mod' corresponds to the HERAPDF modelling and parameterisation uncertainty.
As a result of considering all the data, the experimental uncertainties are reduced with respect to the partial fits.Also, it should be noted that the values of α s extracted with different PDF sets show good agreement with each other within the PDF uncertainties, and are compatible with the latest world average value α s (m Z ) = 0.1181 ± 0.0011 [79].
A comparison of the results for α s from the global and partial fits is shown in Figure 9.In this figure, the results from previous experiments [41, 76, 80-83, 85, 86] are also shown, together with the world average band [79].Agreement between this result and the ones from other experiments is very good, even though the experimental uncertainties in this analysis are smaller than in previous measurements in hadron colliders.

Fits to individual ATEEC functions
The values of α s extracted from the fits to the measured ATEEC functions are summarised in Table 5, together with the values of the χ 2 functions at the minima.
The values extracted from the ATEEC show smaller scale uncertainties than their counterpart values from TEEC.This is understood to be due to the fact that the scale dependence is mitigated for the ATEEC distributions because, for the TEEC, this dependence shows some azimuthal symmetry.Also, it is important to note that the values of the χ 2 indicate excellent compatibility between the data and the theoretical predictions.Good agreement, within the scale uncertainty, is also observed between these values and the ones extracted from fits to the TEEC, as well as among themselves and with the current world average.The nuisance parameters are compatible with zero within one standard deviation.
As before, the values of α s (Q 2 ) at the scales of the measurement are obtained by evolving the values in Table 5 using Equation (5).The results are given in Table 6.As in the TEEC case, Figure 10 compares the data with the theoretical predictions after the fit.The appendix includes tables in which the values of α s (m Z ) obtained from the ATEEC fits are extrapolated to different values of Q, given by the averages of kinematic quantities other than H T2 /2.  ) obtained from fits to the ATEEC function for each H T2 interval using the NNPDF 3.0 parton distribution functions.The uncertainty referred to as NP is the one related to the non-perturbative corrections.

Global ATEEC fit
As before, the global value of α s (m Z ) is obtained from the combined fit of the ATEEC data in the six bins of H T2 .Again, the NNPDF 3.0 PDF set is used for the final result as it provides the most conservative choice.Also, as in the TEEC case, two additional correlation scenarios have been considered for the JES uncertainty.An additional uncertainty of 0.0003 is assigned in order to cover the differences.
The results are summarised in The values shown in Table 7 are in good agreement with the values in Table 4, obtained from fits to the TEEC functions.Also, it is important to note that the scale uncertainty is smaller in ATEEC fits than in TEEC fits.The values of the χ 2 function at the minima show excellent agreement between the data and the pQCD predictions.

Conclusion
The TEEC and ATEEC functions are measured in 20.2 fb −1 of pp collisions at a centre-of-mass energy √ s = 8 TeV using the ATLAS detector at the LHC.The data, binned in six intervals of the sum of transverse momenta of the two leading jets, H T2 = p T1 + p T2 , are corrected for detector effects and compared to the predictions of perturbative QCD, corrected for hadronisation and multi-parton interaction effects.The results show that the data are compatible with the theoretical predictions, within the uncertainties.
The data are used to determine the strong coupling constant α s and its evolution with the interaction scale Q = (p T1 + p T2 )/2 by means of a χ 2 fit to the theoretical predictions for both TEEC and ATEEC in each energy bin.Additionally, global fits to the TEEC and ATEEC data are performed, leading to α s (m Z ) = 0.1162 ± 0.0011 (exp.)+0.0076 −0.0061 (scale) ± 0.0018 (PDF) ± 0.0003 (NP), α s (m Z ) = 0.1196 ± 0.0013 (exp.)+0.0061 −0.0013 (scale) ± 0.0017 (PDF) ± 0.0004 (NP), respectively.Conservatively, the values obtained using the NNPDF 3.0 PDF set are chosen, as they provide the largest PDF uncertainty among the four PDF sets investigated.These two values are in good agreement with the determinations in previous experiments and with the current world average α s (m Z ) = 0.1181 ± 0.0011.The correlation coefficient between the two determinations is ρ = 0.60.
The present results are limited by the theoretical scale uncertainties, which amount to 6% of the value of α s (m Z ) in the case of the TEEC determination and to 4% in the case of the ATEEC.This uncertainty is expected to decrease as higher orders are calculated for the perturbative expansion.

Figure 2 :
Figure 2: Systematic uncertainties in the measured TEEC (top) and ATEEC distributions (bottom) for the first and the last bins in H T2 .The total uncertainty is below 5% in all bins of the TEEC distributions.

Figure 3 :
Figure3: Particle-level distributions for the TEEC functions in each of the H T2 intervals chosen in this analysis, together with MC predictions from Pythia8, Herwig++ and Sherpa.The total uncertainty, including statistical and other experimental sources is also indicated using an error bar for the distributions and a green-shaded band for the ratios.

Figure 5 :
Figure 5: Non-perturbative correction factors for TEEC in the first and last bins of H T2 as a function of cos φ.

Figure 6 :
Figure6: Ratios of the TEEC data in each H T2 bin to the NLO pQCD predictions obtained using the NNPDF 3.0 parton distribution functions, and corrected for non-perturbative effects.

Figure 7 :
Figure7: Ratios of the ATEEC data in each H T2 bin to the NLO pQCD predictions obtained using the NNPDF 3.0 parton distribution functions, and corrected for non-perturbative effects.

Figure 8 :
Figure 8: Comparison of the TEEC data and the theoretical predictions after the fit.The value of α s (m Z ) used in this comparison is fitted independently for each energy bin.

Figure 9 :
Figure9: Comparison of the values of α s (Q) obtained from fits to the TEEC functions at the energy scales given by H T2 /2 (red star points) with the uncertainty band from the global fit (orange full band) and the 2016 world average (green hatched band).Determinations from other experiments are also shown as data points.The error bars, as well as the orange full band, include all experimental and theoretical sources of uncertainty.The strong coupling constant is assumed to run according to the two-loop solution of the RGE.

Figure 10 :
Figure 10: Comparison of the ATEEC data and the theoretical predictions after the fit.The value of α s (m Z ) used in this comparison is fitted independently for each energy bin.

Table 1 :
Summary of the H T2 bins used in the analysis.The table shows the number of events falling into each energy bin together with the value of the scale Q at which the coupling constant α s is measured.

Table 2 :
Values of the strong coupling constant at the Z boson mass scale, α s (m Z ) obtained from fits to the TEEC function for each H T2 interval using the NNPDF 3.0 parton distribution functions.The values of the average scale Q for each energy bin are shown in the first column, while the values of the χ 2 function at the minimum are shown in the third column.The uncertainty referred to as NP is the one related to the non-perturbative corrections.

Table 4 :
Table 4 for each of the four PDF sets investigated in this analysis

Table 6 :
Values of the strong coupling constant at the measurement scales, α s (Q 2

Table 7 :
Table 7for the four sets of PDFs considered in the theoretical predictions.The results for α s from fits to the ATEEC using different PDFs.The uncertainty referred to as NP is the one related to the non-perturbative corrections.The uncertainty labelled as 'mod' corresponds to the HERAPDF modelling and parameterisation uncertainty.

Table 6
Comparison of the values of α s (Q) obtained from fits to the ATEEC functions at the energy scales given by H T2 /2 (red star points) with the uncertainty band from the global fit (orange full band) and the 2016 world average (green hatched band).Determinations from other experiments are also shown as data points.The error bars, as well as the orange full band, include all experimental and theoretical sources of uncertainty.The strong coupling constant is assumed to run according to the two-loop solution of the RGE.

Table 10 :
Values of α s , obtained from TEEC fits, evolved to the average value of transverse momentum for every pair of jets in each event, (p Ti + p T j /2 for each bin in H T2 .