Determination of the strong coupling constant α s m Z from measurements of the total cross section for top-antitop-quark production.

We present a determination of the strong coupling constant α s m Z using inclusive top-quark pair production cross section measurements performed at the LHC and at the Tevatron. Following a procedure first applied by the CMS Collaboration, we extract individual values of α s m Z from measurements by different experiments at several centre-of-mass energies, using QCD predictions complete in NNLO perturbation theory, supplemented with NNLL approximations to all orders, and suitable sets of parton distribution functions. The determinations are then combined using a likelihood-based approach, where special emphasis is put on a consistent treatment of theoretical uncertainties and of correlations between various sources of systematic uncertainties. Our final combined result is α s m Z = 0 . 1177 - 0.0036 + 0.0034 .

Introduction

The strong coupling constant of Quantum Chromodynamics (QCD), , is, together with the quark masses, the main free parameter of the QCD Lagrangian. It enters into every process that involves the strong interaction and is the fundamental parameter of the perturbative expansion used in calculating cross sections for processes with large momentum transfers.

The strong coupling is a function of a renormalisation scale . Its dependence on is governed by renormalisation group equations [1, 2], however, its value at a given reference scale must be determined from experimental data. The current world average value for the coupling evaluated at the Z-boson mass scale, , as determined by the Particle Data Group (PDG), is  [3]. The world average incorporates information from a wide variety of experimental data and of methods to deduce from that data. It requires at least next-to-next-to-leading order (NNLO) accuracy in the perturbative expansions that are used.

Even with the precision that is quoted by the PDG, the uncertainty on contributes significantly to uncertainties on physical predictions for colliders. For example, it leads to about uncertainty on the gluon-fusion Higgs cross section, comparable with the largest of any of the other individual uncertainties [4]. Furthermore, while the bulk of the evidence points to values of the strong coupling that are compatible with , including precise lattice-QCD-based determinations, e.g. [3, 5–7], there are a handful determinations with small quoted uncertainties that suggest values that are several standard deviations below the world average. Notable cases are those from the thrust and C-parameter distributions in collisions, which yield and , respectively [8, 9],1 or the ABMP PDF fit [11], .

Of the various NNLO determinations of the strong coupling, so far only one is based on hadron-collider data, using a measurement of the top-quark pair production cross section () performed by the CMS Collaboration at a centre-of-mass energy TeV [12]. It yields . This extraction is intriguingly placed between the world average and the outlying low extractions, albeit compatible with both. However, it is based on a single, early and now outdated measurement of . It is of interest, therefore, to examine how it is affected by more recent precise measurements by the ATLAS and CMS Collaborations at CERN’s Large Hadron Collider (LHC) [13-17] as well as by a combination of measurements from the D0 and CDF Collaborations at the Tevatron [18].

In the course of our discussion, we will encounter issues related to the treatment of theoretical uncertainties and the choice of the parton distribution function (PDF) set that are of relevance more generally in the determination of the strong coupling and other fundamental constants (e.g. the top-quark mass) from collider data. Such studies may become increasingly widespread in the coming years, given the recent rapid progress in NNLO calculations, e.g. for vector-boson (e.g. Refs. [19, 20]) and inclusive jet distributions [21] at hadron colliders and jet distributions in Deep Inelastic Scattering (DIS) [22].

Determination of from cross section measurements

Theory prediction for the top-pair-production cross section

Theory predictions for the dependence of on are calculated using the program top++2.0 [23]. It provides the computation of the total cross section up to NNLO [24], with possible inclusion of soft-gluon resummation at next-to-next-to-leading logarithmic order (NNLL), as described in Refs. [25, 26].

The predicted cross section is evaluated setting both the renormalisation scale and the factorisation scale equal to the top-quark pole mass. The theoretical uncertainty associated with missing higher-order contributions is evaluated by independently varying and up and down by a factor of 2, under the constraint that . The scale uncertainties are modelled as corresponding to a confidence interval with a Gaussian-shaped uncertainty profile. This choice is more conservative than the (flat) confidence interval that is sometimes taken for scale variations and used, notably, in Ref. [12]. The latter choice leads to a scale uncertainty contribution that is smaller by a factor (the ratio of the standard deviations of the two uncertainty profiles). Note that a confidence interval for scale uncertainties is known to be inconsistent with the observation that a significant fraction of NNLO calculations is outside the scale uncertainty interval of the corresponding NLO calculation.2

A further choice that needs to be made is whether to include the NNLL threshold resummation for the cross section. This is a procedure that resums terms whose leading-logarithmic (LL) structure is , where and s is the squared centre-of-mass energy. When approaches one, i.e. when one approaches the threshold for production, N is proportional to ) and the threshold resummation is a necessity. However, at the LHC and even at the Tevatron, top-pair production is far from threshold and N is not especially large: for the dominant gluon–gluon production channel at LHC, for and ; while for the dominant production channel at the Tevatron, . Accordingly, there is debate within the community as to whether threshold resummation is called for. On one hand, one may argue that it brings about terms that have a certain physical meaning. On the other, one may argue that there is no reason why the terms brought by threshold resummation should dominate over other, neglected terms, and therefore it is more consistent to include just the fixed-order contributions, which are known exactly. We will take an agnostic approach to this question, carry out fits with and without NNLL resummation, and then average both the central values and the uncertainties in the two cases in order to obtain our final result.

The theory prediction for also depends on a choice of PDF set. Since that choice needs to be related to the data that we fit, we postpone our discussion of the PDF choice to Sect. 2.3.

Measurements of the top-pair production cross section

Our determination is performed using seven inputs, listed in Table 1. The six measurements at the LHC include three updated measurements by the CMS Collaboration at centre-of-mass energies of 7, 8 TeV [13] and 13 TeV [14]. These measurements were performed in the decay channel,3 where the W-bosons from the top-quark decays each themselves decay into a charged lepton and a neutrino, one of the W decays producing an electron, the other producing a muon. The measurements are based on data collected in the years of 2011, 2012 and 2015, respectively, with integrated luminosities of 5.0, 19.7, and 2.2 fb. From the ATLAS Collaboration, three similar measurements performed in the decay channel are included, based on datasets with integrated luminosities of 4.6, 20.3 and 3.2 fb for the 7, 8 TeV [16] and 13 TeV [17] centre-of-mass energies, respectively. A seventh input from the Tevatron collider [18] at a centre-of-mass energy of 1.96 TeV is included, which comprises a combination of measurements performed in multiple decay channels from both the CDF Collaboration and the D0 Collaboration.

Table 1

Cross sections and experimental uncertainties for the inputs that we use [13, 14, 16–18]. The LHC beam-energy uncertainties quoted in these references have been scaled down by a factor 6.6 in the light of the recent beam-energy calibration [30], which has a 0.1% uncertainty and coincides with the nominal energy within uncertainties. The original beam-energy-induced uncertainties corresponded to  [31]. The Tevatron beam-energy uncertainty is sufficiently small (cf. Ref. [32]) that no beam-energy uncertainty is quoted by CDF and D0 in the cross section measurements. The cross section and uncertainties listed here are adjusted to the top mass corresponding to the latest world average value computed by the Particle Data Group [3], . The “Exp. ” column corresponds to the uncertainty discussed in Sect. 2.4, signed such that the upper (lower) uncertainty corresponds to an increase (decrease) in

	σ_tt¯ (pb)	Statistical unc. (%)	Systematic unc. (%)	Luminosity unc. (%)	E_beam unc. (%)	Exp. m_t unc. (%)
ATLAS (7 TeV) [16]	182.5	1.7	2.3	2.0	0.3	+0.2-0.2
ATLAS (8 TeV) [16]	242.4	0.7	2.3	2.1	0.3	+0.2-0.2
ATLAS (13 TeV) [17]	816.3	1.0	3.3	2.3	0.2	+0.3-0.3
CMS (7 TeV) [13]	173.4	1.2	2.5	2.2	0.3	+0.2-0.2
CMS (8 TeV) [13]	244.1	0.6	2.4	2.6	0.3	+0.4-0.4
CMS (13 TeV) [14]	809.8	1.1	4.7	2.3	0.2	+0.8-0.8
Tevatron (1.96 TeV) [18]	7.52	2.7	3.9	2.8	0.0	+1.4-1.1

Choice of PDF

Several considerations arise in our choice of PDF. Firstly, we restrict our attention to recent global fits that are available through the LHAPDF interface [33]. Secondly, we require that the PDFs should be available for at least three values, so that we can correctly determine the dependence of the cross section in the context of that PDF. These two conditions limit us to the CT14 [34], MMHT2014 [35] and the NNPDF3.0 [36] series. Thirdly, we impose a requirement that the PDF should not have included data in its fitting procedure. As should be obvious qualitatively, and as we will discuss quantitatively elsewhere [37], using a PDF with top-data included would bias our fits.

Table 2 summarises what data has been included in each of these PDF sets, including both the default NNPDF30 set and NNPDF30_nolhc, obtained without LHC data. One sees that the two options that are available to us are CT14 and NNPDF30_nolhc.4

Table 2

Top-pair cross section data included in a selection of recent PDF fits. A “” (“”) indicates that the corresponding cross section measurement is (is not) included in the PDF fit. The specific sets of 7 and ATLAS and CMS data used in the fits do not always coincide with those that we list in Table 1. All the PDFs shown here predate the measurements

	Tevatron	ATLAS (7 TeV)	ATLAS (8 TeV)	CMS (7 TeV)	CMS (8 TeV)
CT14 [34]	-	-	-	-	-
MMHT2014 [35]	✓	✓	-	✓	✓
NNPDF30 [36]	-	✓	✓	✓	✓
NNPDF30_noLHC [36]	-	-	-	-	-

We use PDF uncertainties calculated at the 68% confidence level, following the error propagation prescription from the individual PDF groups. The uncertainties from the CT14 PDF set, which are provided at a 90% confidence level by default, are scaled by a factor of .

The predicted cross sections for both PDF sets, with NNLO and NNLO + NNLL calculations, are listed in Table 3. The cross sections are higher when including NNLL contributions. The scale uncertainties are in the range for the NNLO results and get reduced by between one third and one half when including NNLL terms. At LHC energies, the cross sections with NNPDF30_nolhc are about larger than those with CT14, however, the opposite pattern is seen at Tevatron. Finally, the PDF uncertainties are somewhat larger with NNPDF_nolhc than with CT14.

Table 3

Predicted cross sections and uncertainties for the PDF sets that we use [34, 36], as determined with the Top++ program [23] at a reference value of . The results are for and the “” column corresponds to the uncertainty discussed in Sect. 2.4, signed such that the upper (lower) uncertainty corresponds to an increase (decrease) in

	σtt¯pred(αsref) [pb]	PDF unc. [%]	Scale unc. [%]	m_t unc. [%]	dlnσ_tt¯(αsref)dlnα_s
CT14 (NNLO)
LHC (7 TeV)	172.7	-3.8+4.5	-6.5+4.1	+2.7-2.6	2.486
LHC (8 TeV)	246.7	-3.5+4.0	-6.3+3.9	+2.6-2.5	2.404
LHC (13 TeV)	807.3	-2.7+2.6	-5.6+3.5	+2.4-2.3	2.133
Tevatron (1.96 TeV)	7.3	-2.2+3.4	-5.5+3.8	+2.8-2.7	1.757
NNPDF30_nolhc (NNLO)
LHC (7 TeV)	174.8	-5.0+5.0	-6.5+4.1	+2.7-2.6	2.247
LHC (8 TeV)	249.7	-4.4+4.4	-6.3+3.9	+2.6-2.5	2.099
LHC (13 TeV)	816.2	-2.9+2.9	-5.6+3.5	+2.4-2.3	1.681
Tevatron (1.96 TeV)	7.2	-3.1+3.5	-5.5+3.8	+2.8-2.7	2.396
CT14 (NNLO + NNLL)
LHC (7 TeV)	177.9	-3.7+4.4	-3.5+2.6	+2.7-2.6	2.545
LHC (8 TeV)	253.6	-3.4+3.9	-3.5+2.6	+2.6-2.5	2.459
LHC (13 TeV)	825.9	-2.7+2.6	-3.6+2.4	+2.4-2.3	2.178
Tevatron (1.96 TeV)	7.4	-2.2+3.5	-2.9+1.6	+2.8-2.7	1.842
NNPDF30_nolhc (NNLO + NNLL)
LHC (7 TeV)	180.1	-5.0+4.9	-3.5+2.6	+2.7-2.6	2.296
LHC (8 TeV)	256.7	-4.4+4.3	-3.5+2.6	+2.6-2.5	2.147
LHC (13 TeV)	835.0	-2.8+2.8	-3.6+2.4	+2.4-2.3	1.722
Tevatron (1.96 TeV)	7.3	-3.2+3.6	-2.9+1.5	+2.8-2.7	2.476

To understand the final errors on the determination it is important also to examine how the predicted cross sections depend on , a result of the dependence both of the hard cross section and of the PDFs. This is shown in Fig. 1: points correspond to the values of for which the given PDF is available, and lines correspond to a fit for using a polynomial of . We use polynomials of degree 3 and 1, respectively, for the CT14 and NNPDF30_nolhc PDFs, chosen based on the available number of points and requirements of stability of the extrapolation beyond the available points. A steeper slope of the dependence (also quoted at in the last column of Table 3) leads to a smaller final error on for any given source of uncertainty on . For LHC energies, CT14 is generally steeper, while at the Tevatron it is NNPDF_nolhc that is steeper. Note also that CT14 curves have substantial curvature, and this will induce asymmetric uncertainties for , even in the case of uncertainties on the cross section that are symmetric.

Fig. 1

Predicted cross section as a function of . The points are the cross sections calculated using the Top++ program [23], and the line is our polynomial fit. The plot also includes horizontal lines corresponding to the central values of the measured cross sections, adjusted to correspond to the same top mass as the theory cross sections (), cf. Sect. 2.4

Top-mass dependence

The top-quark pole mass is taken to be GeV, which is consistent with the world average value computed by the Particle Data Group [3]. The experimentally measured cross section, , depends on through the acceptance corrections, whose parametrisation is given together with the individual measurements. The uncertainty on the experimentally measured cross section due to the top-quark pole mass is given in Table 1, where the uncertainty was calculated by shifting the top mass up and down by its uncertainty. An increase in the top mass leads to a decrease in the measured total cross section. This is because the experiments effectively measure a fiducial cross section (which is independent of ) and then extrapolate it to a total cross section by dividing by the acceptance for the fiducial cross section. For larger values of the acceptance is larger, since decay products are more likely to pass transverse momentum cuts, and so the resulting total cross section is lower. The theoretically predicted cross section, , also depends on , because of the structure of the underlying hard cross section and the x-dependence of the PDFs, cf. Table 3. It too decreases for an increase in the cross section, and this effect is larger than for the measured cross section.

To define a single error contribution associated with the top-mass uncertainty, it is convenient to absorb these different sources of dependence into an effective predicted cross section,where is the central value of the world average top mass. For , this effective predicted cross section coincides with the actual predicted one.

The final uncertainty on the effective predicted cross section associated with the error of on the world average top mass is then given byThis can be used in our determination in a manner similar to any of the theoretical and PDF uncertainties on the predicted cross section. To a good approximation, the final top-mass uncertainty on the effective cross section is equal to the difference between the percentage uncertainties in Tables 1 and 3.

Strong coupling determination procedure

In the determination of from , the theory prediction is treated as a Bayesian prior (one prior for any given value of ) and the experimental result as a likelihood function. The multiplication of these is the joint posterior probability function from which and its uncertainties are determined after marginalisation of . The procedure is mostly analogous to that used by the CMS Collaboration in Ref. [12].

The construction of the Bayesian prior from the theory dependence necessitates a single probability distribution function given all individual theory uncertainties. The three theory uncertainties are each interpreted as corresponding to an asymmetric Gaussian function: where is the predicted central value at a given value of , and is the positive (negative) uncertainty from a given theory uncertainty source. This function has the advantage that the integral normalises naturally to one, and that the integral from to corresponds to a 68% confidence interval. On average there is a 20% difference between and , up to a difference of about 85% for the most asymmetric uncertainty. The central value for corresponds to the median of the distribution.

The combined probability distribution function of the predicted cross section, , is computed by taking the numerical convolution of the individual asymmetric Gaussian functions:where the convolution is performed such that the probability distribution functions are centred around . While the individual uncertainty distributions contain a discontinuity at , the convolution is a smooth function. The dependence on of the width of the uncertainty band is neglected.5 The probability distribution function of the predicted cross section is multiplied by the probability distribution function of the measured cross section , yielding the joint Bayesian posterior in terms of and . The Bayesian confidence interval of can be computed through marginalisation of the posterior by integrating over :Here, is taken to be independent of . Technically a small dependence on is introduced in through the acceptance corrections; however, in the region of relevance around , the effect of this on the uncertainty of the cross section is below the percent level [12], and can thus be safely neglected. The marginalised joint posterior can be treated as a probability distribution function. The central value for the determination is taken to be the location of the peak of , and the uncertainty is extracted by computing the 68% confidence interval whose left and right bounds are at equal height.6 The procedure is illustrated in Fig. 2, showing the experimental and theory probability distribution functions and the unmarginalised posterior (Fig. 2a) as well as the marginalised posterior with extracted central value and uncertainties (Fig. 2b).

Fig. 2

a The central values and 1 deviations for the predicted cross section (, blue oblique lines) and the experimental cross section (, red horizontal lines) and the product of the probability distribution functions (green shading). The markers on the predicted cross section indicate the fit points from top++2.0. b Marginalisation of the joint posterior with Bayesian confidence interval

The combination of determinations from different experiments necessitates a breakdown of the total uncertainty into components that can be assigned to the individual uncertainty sources. To this end, the determination is repeated each time omitting a different uncertainty source, and the squared difference of the resulting uncertainty with respect to the total uncertainty is computed. A relative contribution to the total uncertainty is then computed per uncertainty source.

Individual results for per measurement

The results of our determination are listed for the CT14nnlo PDF set in Tables 4 and 6 and for the NNPDF30_nolhc PDF set in Tables 5 and 7.

Table 4

Results for the strong coupling evaluated at the Z-boson mass scale and individual uncertainty contributions. These are based on cross sections calculated at NNLO using the CT14nnlo series of PDFs

	Center	Stat.	Syst.	Lumi.	E_beam	PDF	Scale	m_t	Total
ATLAS (7 TeV)	0.1205	-0.0009+0.0007	-0.0012+0.0009	-0.0010+0.0008	-0.0001+0.0001	-0.0021+0.0015	-0.0021+0.0021	-0.0012+0.0009	-0.0036+0.0030
ATLAS (8 TeV)	0.1171	-0.0004+0.0003	-0.0014+0.0011	-0.0013+0.0010	-0.0002+0.0001	-0.0025+0.0017	-0.0026+0.0027	-0.0015+0.0011	-0.0044+0.0037
ATLAS (13 TeV)	0.1187	-0.0006+0.0006	-0.0021+0.0017	-0.0014+0.0012	-0.0001+0.0001	-0.0016+0.0014	-0.0024+0.0026	-0.0013+0.0011	-0.0041+0.0038
CMS (7 TeV)	0.1182	-0.0007+0.0005	-0.0014+0.0010	-0.0013+0.0009	-0.0002+0.0001	-0.0025+0.0017	-0.0025+0.0025	-0.0014+0.0010	-0.0043+0.0035
CMS (8 TeV)	0.1175	-0.0004+0.0003	-0.0015+0.0011	-0.0016+0.0012	-0.0001+0.0001	-0.0024+0.0017	-0.0026+0.0026	-0.0014+0.0010	-0.0044+0.0037
CMS (13 TeV)	0.1183	-0.0007+0.0006	-0.0030+0.0025	-0.0015+0.0013	-0.0001+0.0002	-0.0017+0.0014	-0.0025+0.0026	-0.0010+0.0009	-0.0047+0.0042
Tevatron (1.96 TeV)	0.1202	-0.0018+0.0013	-0.0026+0.0019	-0.0019+0.0014	-0.0000+0.0000	-0.0020+0.0014	-0.0027+0.0024	-0.0009+0.0006	-0.0050+0.0039

Table 6

As in Table. 4, but now using NNLO + NNLL cross sections with the CT14nnlo series of PDFs

	Center	Stat.	Syst.	Lumi.	E_beam	PDF	Scale	m_t	Total
ATLAS (7 TeV)	0.1192	-0.0009+0.0007	-0.0012+0.0010	-0.0010+0.0008	-0.0001+0.0001	-0.0021+0.0016	-0.0014+0.0012	-0.0012+0.0010	-0.0033+0.0027
ATLAS (8 TeV)	0.1158	-0.0004+0.0004	-0.0014+0.0011	-0.0013+0.0011	-0.0002+0.0001	-0.0025+0.0019	-0.0018+0.0016	-0.0015+0.0012	-0.0040+0.0032
ATLAS (13 TeV)	0.1175	-0.0006+0.0005	-0.0020+0.0018	-0.0014+0.0012	-0.0001+0.0001	-0.0016+0.0014	-0.0017+0.0017	-0.0013+0.0012	-0.0037+0.0033
CMS (7 TeV)	0.1168	-0.0007+0.0006	-0.0015+0.0011	-0.0013+0.0010	-0.0002+0.0001	-0.0026+0.0019	-0.0017+0.0014	-0.0015+0.0012	-0.0040+0.0031
CMS (8 TeV)	0.1162	-0.0004+0.0003	-0.0015+0.0012	-0.0016+0.0013	-0.0002+0.0001	-0.0024+0.0018	-0.0018+0.0016	-0.0014+0.0011	-0.0040+0.0032
CMS (13 TeV)	0.1171	-0.0007+0.0006	-0.0029+0.0025	-0.0015+0.0013	-0.0002+0.0001	-0.0017+0.0015	-0.0018+0.0017	-0.0011+0.0009	-0.0043+0.0038
Tevatron (1.96 TeV)	0.1188	-0.0017+0.0014	-0.0025+0.0021	-0.0018+0.0015	-0.0000+0.0000	-0.0020+0.0015	-0.0013+0.0011	-0.0009+0.0007	-0.0043+0.0035

Table 5

As in Table 4, but now using NNLO cross sections with the NNPDF30_nolhc series of PDFs

	Center	Stat.	Syst.	Lumi.	E_beam	PDF	Scale	m_t	Total
ATLAS (7 TeV)	0.1206	-0.0009+0.0009	-0.0013+0.0012	-0.0011+0.0010	-0.0001+0.0002	-0.0027+0.0025	-0.0025+0.0029	-0.0013+0.0012	-0.0043+0.0044
ATLAS (8 TeV)	0.1166	-0.0004+0.0004	-0.0013+0.0012	-0.0012+0.0011	-0.0002+0.0001	-0.0026+0.0024	-0.0026+0.0032	-0.0014+0.0013	-0.0043+0.0045
ATLAS (13 TeV)	0.1183	-0.0007+0.0007	-0.0024+0.0022	-0.0017+0.0016	-0.0001+0.0002	-0.0021+0.0020	-0.0029+0.0035	-0.0015+0.0015	-0.0049+0.0051
CMS (7 TeV)	0.1179	-0.0007+0.0006	-0.0013+0.0013	-0.0012+0.0011	-0.0001+0.0001	-0.0028+0.0025	-0.0025+0.0030	-0.0013+0.0012	-0.0044+0.0045
CMS (8 TeV)	0.1170	-0.0003+0.0003	-0.0014+0.0013	-0.0015+0.0014	-0.0002+0.0001	-0.0026+0.0024	-0.0026+0.0032	-0.0013+0.0012	-0.0044+0.0046
CMS (13 TeV)	0.1178	-0.0008+0.0008	-0.0034+0.0032	-0.0017+0.0016	-0.0002+0.0003	-0.0021+0.0020	-0.0029+0.0034	-0.0011+0.0011	-0.0054+0.0055
Tevatron (1.96 TeV)	0.1205	-0.0014+0.0013	-0.0020+0.0019	-0.0014+0.0014	-0.0000+0.0000	-0.0017+0.0015	-0.0021+0.0023	-0.0007+0.0007	-0.0040+0.0039

Table 7

As in Table 4, but now using NNLO + NNLL cross sections with the NNPDF30_nolhc series of PDFs

	Center	Stat.	Syst.	Lumi.	E_beam	PDF	Scale	m_t	Total
ATLAS (7 TeV)	0.1190	-0.0009+0.0009	-0.0012+0.0012	-0.0011+0.0010	-0.0001+0.0001	-0.0026+0.0025	-0.0015+0.0016	-0.0013+0.0012	-0.0037+0.0036
ATLAS (8 TeV)	0.1152	-0.0004+0.0004	-0.0013+0.0012	-0.0012+0.0011	-0.0001+0.0001	-0.0025+0.0024	-0.0017+0.0018	-0.0014+0.0013	-0.0037+0.0037
ATLAS (13 TeV)	0.1168	-0.0007+0.0007	-0.0023+0.0022	-0.0016+0.0015	-0.0002+0.0002	-0.0020+0.0019	-0.0020+0.0022	-0.0015+0.0015	-0.0043+0.0043
CMS (7 TeV)	0.1163	-0.0006+0.0006	-0.0013+0.0012	-0.0011+0.0011	-0.0001+0.0001	-0.0027+0.0026	-0.0016+0.0016	-0.0013+0.0012	-0.0038+0.0037
CMS (8 TeV)	0.1155	-0.0003+0.0003	-0.0013+0.0013	-0.0014+0.0014	-0.0001+0.0001	-0.0025+0.0024	-0.0017+0.0017	-0.0013+0.0012	-0.0038+0.0037
CMS (13 TeV)	0.1163	-0.0007+0.0008	-0.0032+0.0031	-0.0016+0.0015	-0.0002+0.0002	-0.0020+0.0019	-0.0020+0.0022	-0.0011+0.0011	-0.0048+0.0047
Tevatron (1.96 TeV)	0.1194	-0.0013+0.0013	-0.0019+0.0018	-0.0014+0.0013	-0.0000+0.0000	-0.0017+0.0016	-0.0010+0.0010	-0.0007+0.0007	-0.0034+0.0033

The individual determinations are all compatible with the world average to within uncertainties. The central values are rather similar to the CT14 and NNPDF sets. The largest individual sources of uncertainty on are the PDF uncertainties and the scale uncertainties. For the LHC determinations, the PDF uncertainties tend to be larger with NNPDF, in part a consequence of the larger uncertainties in the cross section in Table 3. However, the other uncertainties are also larger with NNPDF, because of its weaker dependence on .

The NNLO + NNLL determinations all have smaller results, consistent with the larger cross sections in Table 3. The scale uncertainties are also noticeably smaller. Other uncertainties are largely unchanged.

A final comment concerns the somewhat larger scale, and PDF uncertainties with the CT14 PDF for the CMS 7 TeV case as compared to the ATLAS 7 TeV case, or also ATLAS 8 TeV as compared to ATLAS 7 TeV. In general with the CT14 PDF, a smaller value of corresponds to larger uncertainties, because the dependence of the cross section is weaker for small values; cf. Fig. 1. Note, however, that the scale and other uncertainties on the cross section predictions have been evaluated only for the reference value of , and in general the question of how one should correlate uncertainties with the central value is a delicate one.7 Accordingly one should be wary of reading too much into the variation of uncertainties with the central value.

Combination of determinations

Correlation coefficients

A combination of measurements can strongly depend on the assumed or calculated correlations [39]. It is therefore necessary to carefully evaluate the correlation coefficients used for the combination. In the case of determinations many correlations can be reasonably motivated or computed. The correlation coefficients between individual measurements are motivated per uncertainty source.The uncertainties on the predicted cross sections (due to the PDF, the top-quark mass and the renormalisation and factorisation scale) are generally strongly correlated. The combination result strongly depends on the assumed correlation structure of these theoretical uncertainties if included in the combination, which is usually not known precisely in particular for the scale uncertainty. We therefore adopt a different procedure: The individual results are simultaneously shifted up and down by their respective total theory uncertainties, and the combination is re-evaluated. The difference between the upper and lower bounds and the original combination is taken to be the (asymmetric) theoretical uncertainty.

Statistical uncertainties are considered uncorrelated for all experimental inputs.

Systematic uncertainties are considered fully correlated only for measurements obtained with the same detector. This concerns the measurements performed by CMS and ATLAS at different centre-of-mass energies.

Uncertainties due to beam energy are fully correlated between ATLAS and CMS and are taken to be correlated across energies. The beam-energy uncertainty at the Tevatron was tiny and is neglected, as outlined in the caption of Table 1.

Uncertainties due to luminosity are partially correlated between ATLAS and CMS. The correlated component of the luminosity uncertainty stems from the uncertainty on the bunch current density and similarities in the Van der Meer scan fit model.

The correlated and uncorrelated uncertainties are estimated using the same principles as used for the top-quark-pair production cross section combinations between ATLAS and CMS at 7 and 8 TeV [40, 41], updated with the latest luminosity determinations [42-46]. The luminosity uncertainty (as a percentage of the top-quark-pair production cross section) is displayed in Table 8. The luminosity uncertainties on are taken to have the same correlation coefficient.

Table 8

Correlated, uncorrelated and total luminosity uncertainties with respect to the top-quark-pair production cross section (in percentages) [40–46]

s	Experiment	Corr.	Uncorr.	Total
7 TeV	ATLAS	0.46	1.72	1.78
	CMS	0.46	2.13	2.17
8 TeV	ATLAS	0.60	1.84	1.94
	CMS	0.68	2.50	2.59
13 TeV	ATLAS	0.36	2.29	2.32
	CMS	0.36	2.31	2.34

The impact of the alternative procedure of including also the theory uncertainties within a single combination is discussed in Appendix A.

Combining correlated measurements: likelihood-based approach

In order to combine the individual results, we opted for a likelihood-based approach [47].8 In this approach a global likelihood function is constructed from the probability distribution functions of individual determinations. Let us suppose we have measurements of the top cross section and associated determinations of . For each determination i, , we have uncorrelated error components, each specific to that determination. The magnitude of the kth uncorrelated error for determination i is labelled . We additionally have error components that are correlated across all determinations. For each of the correlated components, j, we introduce a nuisance parameter that is common across all measurements. Its impact on measurement i is governed by a coefficient . The full set of will be denoted .

The likelihood will be composed of a product of probability distribution functions (pdfs).9 For each nuisance parameter we will have one pdf, a Gaussian distribution with a standard deviation of one:

There will also be a pdf for each combination of measurement i and associated uncorrelated error . It is given byTo address the issue of errors that are not symmetric, we adopt the following prescription for the and : An overview of the values used for and is given in Appendix B. The probability distribution function of determination i including all uncorrelated uncertainties is then constructed by convolution: where the convolution is performed such that the probability distribution functions are centred around . The global likelihood function is constructed by multiplication of the probability distribution functions of the determinations and the nuisance parameters:In order to complete the formalism of a statistical test the test statistic q is introduced:Here L is maximised for variables that carry a hat and in general will take on different values from . The quantity is therefore the global maximum likelihood, and the ratio cannot be larger than one. The normalisation is such that q can be treated as -distributed with one degree of freedom.

Fig. 3

Combination results using the CT14 PDF set (NNLO in a and NNLO + NNLL in b) and the NNPDF3.0 noLHC PDF set (NNLO in c and NNLO + NNLL in d). The individual determinations and their uncertainties are shown in grey, where the darker shade represents the experimental uncertainties which enter into the combination. The test statistic q as a function of is plotted as a black line. The green line and band represent the central value of the combination and the confidence interval, respectively. The red band depicts the total combination uncertainty with scale, PDF and top-mass uncertainties included

The test statistic q is scanned over a range of values. The minimum of the scan, by construction at , is the maximum likelihood value for , and the confidence interval is extracted from the interval between the intersection points of the scan with . Any skewness of the parabola of the scan is due to the inclusion of asymmetric uncertainties. Figure 3 shows the scan and the corresponding combination results for each of the PDF sets.

Results and discussion

The combination procedure is performed for each of the two PDF sets taken into consideration at NNLO and at NNLO + NNLL separately. The combination results and their unweighted average are displayed numerically in Table 9, and graphically in Fig. 4.

Table 9

Combination results for all PDF sets taken into consideration, at NNLO and NNLO + NNLL

	Center	Stat.	Syst.	E_beam	Lumi.	m_t	PDF	Scale	Total
CT14 (NNLO)	0.1184	-0.0003+0.0003	-0.0007+0.0006	-0.0001+0.0001	-0.0006+0.0006	-0.0014+0.0010	-0.0023+0.0016	-0.0025+0.0025	-0.0038+0.0033
NNPDF30_nolhc (NNLO)	0.1182	-0.0003+0.0003	-0.0007+0.0007	-0.0000+0.0000	-0.0007+0.0007	-0.0013+0.0012	-0.0025+0.0023	-0.0026+0.0031	-0.0040+0.0042
CT14 (NNLO + NNLL)	0.1172	-0.0003+0.0003	-0.0007+0.0007	-0.0001+0.0001	-0.0007+0.0006	-0.0014+0.0011	-0.0023+0.0017	-0.0017+0.0015	-0.0033+0.0027
NNPDF30_nolhc (NNLO + NNLL)	0.1168	-0.0003+0.0003	-0.0007+0.0006	-0.0001+0.0001	-0.0007+0.0007	-0.0013+0.0012	-0.0024+0.0023	-0.0017+0.0018	-0.0034+0.0033
Average	0.1177	+0.0003+0.0003	+0.0007+0.0007	+0.0001+0.0001	+0.0007+0.0006	+0.0013+0.0012	+0.0024+0.0020	+0.0021+0.0022	-0.0036+0.0034

Fig. 4

Combination results for all PDF sets taken into consideration, at NNLO and NNLO + NNLL. The solid blue line is the unweighted average of the individual combination results, and the dashed blue lines represent the 68% confidence interval. The red and green bands are as in Fig. 3

There is no unique way to quote a final best estimate of based on the results obtained from the different PDF sets and QCD calculation choices (NNLO v. NNLO + NNLL). An unbiased approach for combining results from different PDFs, in line with the PDF4LHC recommendations [49], is to average without applying any further weighting. In accordance with that approach we take the straight average of the mean values and the uncertainties of the individual combinations. This coincides with the procedure for combining results from a single class of observables in Ref. [3]. The final result iswhich can be compared to the result of Ref. [12], . Our central value is larger mainly because recent measurements of the cross sections are higher than that used in Ref. [12], but also in part because of our choice to take the average of results from NNLO and NNLO + NNLL cross sections (a increase relative to just NNLO + NNLL). Our symmetrised uncertainty of is somewhat increased with respect to that of Ref. [12], (symmetrised). The difference in uncertainty is due to several choices. On one hand we have taken a smaller uncertainty on the top-quark mass, in line with the PDG determination. One the other hand, we have been somewhat more conservative in our treatment of theoretical and PDF uncertainties. Firstly, the choice of treating the scale uncertainties on as a confidence interval instead of a (flat) confidence interval increases the scale uncertainty component by roughly a factor of . Secondly, we have used an average of the uncertainties from NNLO and NNLO + NNLL cross section determinations, which also yields a larger uncertainty than using NNLO + NNLL cross section determinations only. Finally, the PDF sets used for the determination were chosen with minimisation of potential biases in mind, rather than the ones with smallest uncertainty.

Conclusions

We have used seven measurements of the top–antitop-quark production cross section at the LHC and the Tevatron in order to determine the strong coupling constant , using the CT14 PDF set and the NNPDF30_nolhc PDF set at NNLO and NNLO + N NLL. Overall, our determination of yields a value that is compatible with the world average value and uncertainties that are somewhat larger than the best individual determinations, though comparable with that from the electroweak precision data [50]. The largest uncertainties are associated with unknown higher-order contributions and PDF uncertainties.

Table 10

Combination results including also uncertainties from the PDF, the scale and the top mass in the combination. The first row corresponds to our approach of excluding correlated uncertainties from the combination. The results listed here are obtained using NNLO cross sections with the CT14nnlo series of PDFs

Uncertainties included in combination	Center	Combination uncertainty	Total uncertainty
–	0.1184	-0.0010+0.0009	-0.0038+0.0033
PDF	0.1191	-0.0018+0.0018	-0.0034+0.0033
PDF and m_t	0.1194	-0.0020+0.0020	-0.0032+0.0032
PDF, m_t and scale	0.1207	-0.0029+0.0030	-0.0029+0.0030

Table 11

As in Table 10, but now using NNLO + NNLL cross sections with the CT14nnlo series of PDFs

Uncertainties included in combination	Center	Combination uncertainty	Total uncertainty
–	0.1172	-0.0010+0.0010	-0.0033+0.0027
PDF	0.1180	-0.0020+0.0019	-0.0029+0.0027
PDF and m_t	0.1183	-0.0022+0.0022	-0.0027+0.0027
PDF, m_t and scale	0.1188	-0.0025+0.0025	-0.0025+0.0025

Table 12

As in Table 10, but now using NNLO cross sections with the NNPDF30_nolhc series of PDFs

Uncertainties included in combination	Center	Combination uncertainty	Total uncertainty
–	0.1182	-0.0010+0.0010	-0.0040+0.0042
PDF	0.1188	-0.0022+0.0023	-0.0037+0.0040
PDF and m_t	0.1190	-0.0024+0.0025	-0.0036+0.0040
PDF, m_t and scale	0.1200	-0.0036+0.0035	-0.0036+0.0035

Table 13

As in Table 10, but now using NNLO + NNLL cross sections with the NNPDF30_nolhc series of PDFs

Uncertainties included in combination	Center	Combination uncertainty	Total uncertainty
–	0.1168	-0.0010+0.0010	-0.0034+0.0033
PDF	0.1175	-0.0023+0.0023	-0.0031+0.0031
PDF and m_t	0.1178	-0.0025+0.0025	-0.0030+0.0030
PDF, m_t and scale	0.1182	-0.0028+0.0028	-0.0028+0.0028

Table 14

Overview of the uncertainty coefficients for the CT14 (NNLO) PDF set. The coefficients with a correspond to the coefficients of the correlated uncertainty sources, and those with a to the uncorrelated uncertainty sources

Exp.	δ^Syst.	δ^Lumi.	δ^E_beam	Δ^Stat.	Δ^Lumi.
ATLAS (13 TeV)	-0.0021+0.0017	-0.0002+0.0002	-0.0001+0.0001	-0.0006+0.0006	-0.0014+0.0012
ATLAS (8 TeV)	-0.0014+0.0011	-0.0004+0.0003	-0.0002+0.0001	-0.0004+0.0003	-0.0013+0.0009
ATLAS (7 TeV)	-0.0012+0.0009	-0.0003+0.0002	-0.0001+0.0001	-0.0009+0.0007	-0.0010+0.0007
CMS (13 TeV)	-0.0030+0.0025	-0.0002+0.0002	-0.0001+0.0002	-0.0007+0.0006	-0.0015+0.0013
CMS (8 TeV)	-0.0015+0.0011	-0.0004+0.0003	-0.0001+0.0001	-0.0004+0.0003	-0.0016+0.0012
CMS (7 TeV)	-0.0014+0.0010	-0.0003+0.0002	-0.0002+0.0001	-0.0007+0.0005	-0.0013+0.0009
Tevatron (1.96 TeV)	-0.0026+0.0019	–	–	-0.0018+0.0013	-0.0019+0.0014

Table 15

As in Table 14, but now using NNLO + NNLL cross sections with the CT14 PDF set

Exp.	δ^Syst.	δ^Lumi.	δ^E_beam	Δ^Stat.	Δ^Lumi.
ATLAS (13 TeV)	-0.0020+0.0018	-0.0002+0.0002	-0.0001+0.0001	-0.0006+0.0005	-0.0014+0.0012
ATLAS (8 TeV)	-0.0014+0.0011	-0.0004+0.0003	-0.0002+0.0001	-0.0004+0.0004	-0.0013+0.0010
ATLAS (7 TeV)	-0.0012+0.0010	-0.0003+0.0002	-0.0001+0.0001	-0.0009+0.0007	-0.0010+0.0008
CMS (13 TeV)	-0.0029+0.0025	-0.0002+0.0002	-0.0002+0.0001	-0.0007+0.0006	-0.0015+0.0013
CMS (8 TeV)	-0.0015+0.0012	-0.0004+0.0003	-0.0002+0.0001	-0.0004+0.0003	-0.0016+0.0012
CMS (7 TeV)	-0.0015+0.0011	-0.0003+0.0002	-0.0002+0.0001	-0.0007+0.0006	-0.0013+0.0010
Tevatron (1.96 TeV)	-0.0025+0.0021	–	–	-0.0017+0.0014	-0.0018+0.0015