Durham Research Online Deposited in DRO: 31 May 2017 Version of attached �le: Published Version Peer-review status of attached �le: Peer-reviewed Citation for published item: Harland-Lang, L. A. and Martin, A. D. and Motylinski, P. and Thorne, R. S. (2016) 'Charm and beauty quark masses in the MMHT2014 global PDF analysis.', European physical journal C., 76 (1). p. 10. Further information on publisher's website: https://doi.org/10.1140/epjc/s10052-015-3843-5 Publisher's copyright statement: c© The Author(s) 2016 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3 Additional information: Use policy The full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-pro�t purposes provided that: • a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders. Please consult the full DRO policy for further details. Durham University Library, Stockton Road, Durham DH1 3LY, United Kingdom Tel : +44 (0)191 334 3042 | Fax : +44 (0)191 334 2971 https://dro.dur.ac.uk https://www.dur.ac.uk https://doi.org/10.1140/epjc/s10052-015-3843-5 http://dro.dur.ac.uk/21920/ https://dro.dur.ac.uk/policies/usepolicy.pdf https://dro.dur.ac.uk Eur. Phys. J. C (2016) 76:10 DOI 10.1140/epjc/s10052-015-3843-5 Regular Article - Theoretical Physics Charm and beauty quark masses in the MMHT2014 global PDF analysis L. A. Harland-Lang1, A. D. Martin2, P. Motylinski1, R. S. Thorne1,a 1 Department of Physics and Astronomy, University College London, London WC1E 6BT, UK 2 Institute for Particle Physics Phenomenology, Durham University, Durham DH1 3LE, UK Received: 13 October 2015 / Accepted: 11 December 2015 / Published online: 6 January 2016 © The Author(s) 2016. This article is published with open access at Springerlink.com Abstract We investigate the variation in the MMHT2014 PDFs when we allow the heavy-quark masses mc and mb to vary away from their default values. We make PDF sets avail- able in steps of �mc = 0.05 GeV and �mb = 0.25 GeV, and present the variation in the PDFs and in the predictions. We examine the comparison to the HERA data on charm and beauty structure functions and note that in each case the heavy-quark data, and the inclusive data, have a slight pref- erence for lower masses than our default values. We provide PDF sets with three and four active quark flavours, as well as the standard value of five flavours. We use the pole mass definition of the quark masses, as in the default MMHT2014 analysis, but briefly comment on the MS definition. 1 Introduction Over the past few years there has been a significant improve- ment both in the precision and in the variety of the data for deep-inelastic and related hard-scattering processes. Since the MSTW2008 analysis [1] we have seen the appearance of the HERA combined H1 and ZEUS data on the total [2] and also on the charm structure functions [3], together with a variety of new hadron-collider data sets from the LHC, and in the form of updated Tevatron data (for full references see [4]). Additionally, the procedures used in the global PDF analyses of data have been improved, allowing the parton distributions of the proton to be determined with more pre- cision and with more confidence. This allows us to improve predictions for Standard Model signals and to model Stan- dard Model backgrounds to possible experimental signals of New Physics more accurately. One area that now needs care- ful attention, at the present level of accuracy, is the treatment of the masses of the charm and beauty quarks, mc and mb, in the global analyses. Here we extend the recent MMHT2014 a e-mail: thorne@hep.ucl.ac.uk global PDF analysis [4] to study the dependence of the PDFs, and the quality of the comparison to data, under variations of these masses away from their default values of mc = 1.4 GeV and mb = 4.75 GeV, as well as the resulting predictions for processes at the LHC. We make available central PDF sets for a variety of masses, namely mc = 1.15–1.55 GeV in steps of 0.05 GeV and mb = 4.25–5.25 GeV in steps of 0.25 GeV. We also make available the standard MMHT2014 PDFs, and the sets with varied masses in the three and four flavour number schemes. 2 Dependence on the heavy-quark masses 2.1 Choice of range of heavy-quark masses In the study of heavy-quark masses that accompanied the MSTW2008 PDFs [5] we varied the charm and beauty quark masses, defined in the pole mass scheme, from mc = 1.05 to 1.75 GeV and mb = 4 to 5.5 GeV. This was a very generous range of masses, and it was not clear that there was a demand for PDFs at the extreme limits. Hence, this time we are a little more restrictive, and study the effects of varying mc from 1.15 to 1.55 GeV, in steps of 0.05 GeV, and of varying mb from 4.25 to 5.25 GeV in steps of 0.25 GeV. Part of the reason for this is that the values are constrained by the comparison to data, though for both charm and beauty the preferred values are at the lower end of the range, as we will show. However, there is also the constraint from other determinations of the quark masses. These are generally quoted in the MS scheme, and in [6] are given as mc(mc) = (1.275 ± 0.025) GeV and mb(mb) = (4.18 ± 0.03) GeV. The transformation to the pole mass definition is not well defined due to the diverging series, i.e. there is a renormalon ambiguity of ∼0.1–0.2 GeV. The series is less convergent for the charm quark, due to the lower scale in the coupling, but the renormalon ambiguity cancels in difference between the charm and beauty masses. 123 http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-015-3843-5&domain=pdf mailto:thorne@hep.ucl.ac.uk 10 Page 2 of 14 Eur. Phys. J. C (2016) 76 :10 Indeed, we obtain m pole b −m pole c = 3.4 GeV with a very small uncertainty [7, 8]. Using the perturbative expression for the conversion of the beauty mass, and the relationship between the beauty and charm mass, as shown in [5], we obtain m pole c = 1.5 ± 0.2 GeV and mpoleb = 4.9 ± 0.2 GeV. (1) This disfavours mc ≤ 1.2–1.3 GeV and mb ≤ 4.6–4.7 GeV. As the fit quality prefers values in this region, or lower, we allow some values a little lower than this. In the upper direc- tion the fit quality clearly deteriorates, so our upper values are not far beyond the central values quoted above. There is some indication from PDF fits for a slightly lower mpole than that suggested by the use of the perturbative series out to the order at which it starts to show lack of convergence. We now consider the variation with mc and mb in more detail. 2.2 Dependence on mc We repeat the global analysis in [4] for values of mc = 1.15– 1.55 GeV in steps of 0.05 GeV. As in [4] we use the “optimal” version [9] of the TR’ general-mass variable-flavour-number scheme GM-VFNS [10]. This is smoother near the transi- tion point, which we define to be at Q2 = μ2 = m2c , than the original version, so has a slight tendency to prefer lower masses—the older version growing a little more quickly at low scales, which could be countered by increasing the mass. We also assume all heavy flavour is generated by evolution from the gluon and light quarks, i.e. there is no intrinsic heavy flavour. We perform the analysis with αS(M 2 Z ) left as a free parameter in the fit at both NLO and NNLO, but also use our fixed default values of the coupling of αS(M 2 Z ) = 0.118 and 0.120 at NLO and αS(M 2 Z ) = 0.118 at NNLO. Unlike the MSTW2008 study [5] we will concentrate on the results and PDFs with fixed coupling, as the standard MMHT PDFs were made available at these values. We present results in terms of the χ 2 for the total set of data in the global fit and for just the data on the reduced cross section, σ̃ cc̄, for open charm production at HERA [3]. This is shown at NLO with αS(M 2 Z ) = 0.120 in Fig. 1. The variation in the quality of the fit to the HERA combined charm cross section data is relatively slight, less than the variation in the fit to the separate H1 and ZEUS data used in [5]. This is presum- ably due to the use of the full information now available on correlated systematics, which allows movement of the data relative to the theory with only a moderate penalty in χ 2. The lower variation is also likely due in part to the improved flavour scheme. Despite the fairly small variation in χ 2 the charm data clearly prefer a value close to mc = 1.35 GeV, near our default value of mc = 1.4 GeV. However, there is more variation in the fit quality to the global data set, with a clear preference for values close to mc = 1.2 GeV. The dete- rioration is clearly such as to make values of mc > 1.5 GeV strongly disfavoured. The main constraint comes from the inclusive HERA cross section data, but there is also some preference for a low value of the mass from NMC structure function data, where the data for x ∼ 0.01 and Q2 ∼ 4 GeV2 is sensitive to the turn-on of the charm contribution to the structure function. Overall, there is some element of tension between the preferred value for the global fit and the fit to charm data. We do not attempt to make a rigorous determi- nation of the best value of the mass or its uncertainty, as provided in [11] for example, as we believe there are more precise and better controlled methods for this. However, a rough indication of the uncertainty could be obtained from the χ 2 profiles by treating mc in the same manner as the stan- dard PDF eigenvectors and applying the dynamic tolerance procedure. In this case the appropriate tolerance, obtained mc [GeV] χ2 σ̃cc (52 pts), NLO, αS(M 2Z) = 0.120 . 1.551.51.451.41.351.31.251.21.15 110 100 90 80 70 60 mc [GeV] χ2global (2996 pts), NLO, αS(M 2 Z) = 0.120 . 1.551.51.451.41.351.31.251.21.15 3360 3340 3320 3300 3280 3260 3240 Fig. 1 The quality of the fit versus the quark mass mc at NLO with αS(M2Z ) = 0.120 for (left) the reduced cross section for charm production σ̃ cc̄ for the combined H1 and ZEUS data and (right) the full global fit 123 Eur. Phys. J. C (2016) 76 :10 Page 3 of 14 10 mc [GeV] χ2 σ̃cc (52 pts), NLO, αS(M 2Z) = 0.118 . 1.551.51.451.41.351.31.251.21.15 110 100 90 80 70 60 mc [GeV] χ2global (2996 pts), NLO, αS(M 2 Z) = 0.118 . 1.551.51.451.41.351.31.251.21.15 3360 3340 3320 3300 3280 3260 3240 Fig. 2 The quality of the fit versus the quark mass mc at NLO with αS(M2Z ) = 0.118 for (left) the reduced cross section for charm production σ̃ cc̄ for the combined H1 and ZEUS data and (right) the full global fit Table 1 The quality of the fit versus the quark mass mc at NLO with αS(M 2 Z ) left as a free parameter mc (GeV) χ 2 global χ 2 σ̃ cc̄ αs(M 2 Z ) 2996 pts 52 pts 1.15 3239 75 0.1190 1.2 3237 73 0.1192 1.25 3239 71 0.1194 1.3 3245 70 0.1195 1.35 3254 70 0.1196 1.4 3268 71 0.1198 1.45 3283 73 0.1200 1.5 3303 76 0.1201 1.55 3327 81 0.1202 by assuming the charm cross section data is the dominant constraint, would be of the order T = √ �χ 2 ≈ 2.5.1 The analogous results for αS(M 2 Z ) = 0.118 and αS(M2Z ) left free are shown in Fig. 2 and Table 1, respectively, where in the latter case the corresponding αs(M 2 Z ) values are shown as well. For αS(M 2 Z ) = 0.118 the picture is much the same as for αS(M 2 Z ) = 0.120 except that the fit to charm data is marginally better, while the global fit is a little worse, though more-so for higher masses. The results with free αS(M 2 Z ) are consistent with this, with the preferred value of αS(M 2 Z ) falling slightly with lower values of mc. However, the values of mc preferred by charm data and the full data sets are much 1 As discussed in [1], for a 68 % confidence level uncertainty we insist the fit quality to a given data set deteriorates by no more than the width of the χ 2 distribution for N points, roughly √ N/2 multiplied by the χ 2 per point for the best fit. For the charm cross section data this is ≈√52/2 × 1.3. the same as for fixed coupling—the values of the χ 2 just being a little lower in general. The results of the same analysis at NNLO are shown for αS(M 2 Z ) = 0.118 and αS(M2Z ) left free in Fig. 3 and Table 2, respectively, where again in the latter case the corresponding αs(M 2 Z ) values are shown. Broadly speaking, the results are similar to those at NLO, but with lower values of mc preferred and where the χ 2 variation is greater for the inclusive data than for the charm cross section data. However, in this case there is essentially no tension at all between the inclusive and charm data, with both χ 2 values minimising very near to mc = 1.25 GeV—this lower preferred value for the charm data meaning that the fit quality at mc = 1.55 GeV has deteriorated more than at NLO. The picture is exactly the same for fixed and free strong coupling, with the values of χ 2 simply being a little lower when αS(M 2 Z ) is left free, since the best fit value of the coupling is a little below 0.118, particularly for low mc. 2.3 Dependence on mb We repeat essentially the same procedure for varying val- ues of mb in the range 4.25–5.25 GeV in steps of 0.25 GeV. However, this time there were no data on the beauty contri- bution to the cross section included in the standard global fit [4]. In the previous heavy-quark analysis [5] we compared to beauty cross section data from H1 [12]. This placed a weak constraint on the value of mb but had negligible constraint on the PDFs for fixed mb. Hence, we did not include these data in the updated global fit [4]. There are now also data of comparable precision from ZEUS [13], and we will include both these data sets in future global fits. In this article we study the quality of the comparison to these data to predic- 123 10 Page 4 of 14 Eur. Phys. J. C (2016) 76 :10 mc [GeV] χ2 σ̃cc (52 pts), NNLO, αS(M 2Z) = 0.118 . 1.551.51.451.41.351.31.251.21.15 110 100 90 80 70 60 mc [GeV] χ2global (2663 pts), NNLO, αS(M 2 Z) = 0.118 . 1.551.51.451.41.351.31.251.21.15 2770 2760 2750 2740 2730 2720 2710 2700 Fig. 3 The quality of the fit versus the quark mass mc at NNLO with αS(M2Z ) = 0.118 for (left) the reduced cross section for charm production σ̃ cc̄ for the combined H1 and ZEUS data and (right) the full global fit Table 2 The quality of the fit versus the quark mass mc at NNLO with αS(M 2 Z ) left free mc (GeV) χ 2 global χ 2 σ̃ cc̄ αs(M 2 Z ) 2663 pts 52 pts 1.15 2703 78 0.1164 1.2 2699 76 0.1166 1.25 2698 75 0.1167 1.3 2701 76 0.1169 1.35 2707 78 0.1171 1.4 2717 82 0.1172 1.45 2729 88 0.1173 1.5 2749 96 0.1173 1.55 2769 105 0.1175 tions obtained using the MMHT PDFs with different values of mb. The data themselves are not included in the fit, i.e. we use predictions from the PDFs, as they still provide negligible direct constraint. The results for the NLO PDFs with αS(M 2 Z ) = 0.120 and 0.118 are shown in Figs. 4 and 5, respectively. The picture for the data in the global fit (not including the σ̃ bb̄ data) is slightly different in the two cases: for αS(M 2 Z ) = 0.120 there is a fairly weak tendency to prefer lower values of mb, similar to the results in [5], but for αS(M 2 Z ) = 0.118 the global fit prefers a value of between 4.5 and 5.0 GeV. For the predictions for the beauty cross section data, however, the picture is similar in the two cases, and low values of mb ∼ 4.4–4.5 GeV are preferred. The results for the NNLO fit with αS(M 2 Z ) = 0.118 are shown in Fig. 6. As can be seen the global fit is fairly weakly dependent on mb, though more than for αS(M 2 Z ) = 0.120 at NLO, and prefers a value lower than mb = 4.25 GeV. As in the NLO case the χ 2 for the prediction for σ̃ bb̄ is better for lower values of mb. The slightly larger variations in the quality of the global fit with varying mb compared to [5] is perhaps due to the greater precision of the inclusive HERA cross section data used in this analysis, and to the fact that the CMS double-differential Drell–Yan data [14] has some sensitivity to the value of mb due to the induced variation in sea quark flavour composition for low scales. The previous analysis preferred a value of mb ∼ 4.75 GeV for the com- parison to the H1 beauty data. However, the definition of the general-mass variable-number scheme has improved since this previous analysis, being smoother near to the transition point Q2 = m2b, and including an improvement to the approx- imation for the O(α3S) contribution at low Q2 at NNLO, so some changes are not surprising. Another important differ- ence is in the treatment of the correlated experimental errors, which we now take as being multiplicative. The result within exactly the same framework, but with the experimental errors on the HERA beauty data instead treated as additive is also shown in Fig. 6 and a higher value of mb ∼ 4.75 GeV is clearly preferred. Similar results are seen in the NLO fits. In Fig. 7 the comparison to the (unshifted) HERA beauty data for different values of mb at NNLO is shown. At low Q2 and for ZEUS data in particular, the curves for lower mb are clearly a better fit to unshifted data. How- ever, the low-m2b predictions do significantly overshoot some of the unshifted data points. These predictions will work better with the multiplicative definition of uncertainties as the size of the correlated uncertainties then scales with the prediction, not the data point (as would be the case in the additive definition), or equivalently, if data are nor- malised up to match theory, then so is the uncorrelated uncertainty. 123 Eur. Phys. J. C (2016) 76 :10 Page 5 of 14 10 mb [GeV] χ2 σ̃bb (29 pts), NLO, αS(M 2Z) = 0.120 . 5.254.84.64.4 50 45 40 35 30 25 20 15 mb [GeV] χ2global (2996 pts), NLO, αS(M 2 Z) = 0.120 . 5.254.84.64.4 3310 3300 3290 3280 3270 3260 Fig. 4 The quality of the fit versus the quark mass mb at NLO with αS(M2Z ) = 0.120 for (left) the reduced cross section for beauty production σ̃ bb̄ for the H1 and ZEUS data and (right) the global fit, not including the beauty data mb [GeV] χ2 σ̃bb (29 pts), NLO, αS(M 2Z) = 0.118 . 5.254.84.64.4 50 45 40 35 30 25 20 15 mb [GeV] χ2global (2996 pts), NLO, αS(M 2 Z) = 0.118 . 5.254.84.64.4 3310 3300 3290 3280 3270 3260 Fig. 5 The quality of the fit versus the quark mass mb at NLO with αS(M2Z ) = 0.118 for (left) the reduced cross section for beauty production σ̃ bb̄ for the H1 and ZEUS data and (right) the global fit, not including the beauty data 2.4 Changes in the PDFs We show how the NLO PDFs for mc = 1.25 and 1.55 GeV compare to the central PDFs in Figs. 8 and 9. Results are very similar at NNLO, though more complicated to inter- pret for the charm distribution at low Q2 due to the non-zero transition matrix element at Q2 = m2c in this case. We see at Q2 = 4 GeV2 (that is, close to the transition point Q2 = m2c ) that the change in the gluon is well within its uncertainty band, though there is a slight increase at smaller x with higher mc (and vice versa) such that extra gluon quickens the evo- lution of the structure function which is suppressed by larger mass. Similarly the light quark singlet distribution increases slightly near the transition point for larger mc to make up for the smaller charm contribution to structure functions, and this is maintained, helped by the increased gluon, at larger scales. In both cases, however, the changes are within uncertainties for these moderate variations in mc. The charm distribution increases at low Q2 for decreasing mc, and vice versa, sim- ply due to increased evolution length ln(Q2/m2c). As men- tioned before we have identified the transition point at which heavy flavour evolution begins with the quark mass. This has the advantage that the boundary condition for evolution is zero up to NLO (with our further assumption that there is no intrinsic charm), though there is a finite O(α2S) boundary condition at NNLO in the GM-VFNS, available in [15]. In principle the results on the charm distribution at relatively low scales, such as that in Fig. 8 are sensitive to these defini- 123 10 Page 6 of 14 Eur. Phys. J. C (2016) 76 :10 mb [GeV] χ2 σ̃bb (29 pts), NNLO, αS(M 2Z) = 0.118 . 5.254.84.64.4 50 45 40 35 30 25 20 15 mb [GeV] χ2global (2663 pts), NNLO, αS(M 2 Z) = 0.118 . 5.254.84.64.4 2728 2726 2724 2722 2720 2718 2716 2714 2712 mb [GeV] χ2 σ̃bb (29 pts), NNLO, αS(M 2Z) = 0.118, additive errorrs . 5.254.84.64.4 50 45 40 35 30 25 20 15 Fig. 6 The quality of the fit versus the quark mass mb at NNLO with αS(M 2 Z ) = 0.118 for (left) the reduced cross section for beauty pro- duction σ̃ bb̄ for the H1 and ZEUS data and (right) the global fit, not including the beauty data. Recall that in the MMHT analysis the exper- imental errors are treated multiplicatively. The lower plot shows the χ 2 profile if the errors in the HERA beauty data were to be treated additively tions at finite order, though as the order in QCD increases the correction for changes due to different choices of the tran- sition point arising from the corresponding changes in the boundary conditions become smaller and smaller, ambigu- ities always being of higher order than the calculation. At scales typical of most of LHC physics, however, the rela- tive change in evolution length for the charm distribution is much reduced, as are the residual effects of choices relating to the choice of the transition point and intrinsic charm. At these scales the change in the charm distribution is of the same general size as the PDF uncertainty for fixed mc, as seen in Fig. 9. We also note that the charm structure func- tion at these high scales is reasonably well represented by the charm distribution, while at low scales, certainty including Q2 = 4 GeV2, this is not true. Indeed at NNLO the boundary condition for the charm distribution is negative at very low x if the transition point is m2c , but this is more than compensated for by the gluon and light quark initiated cross section. As noted in [9], use of a zero mass scheme becomes unfeasible at NNLO. The dependence on the heavy-quark cross section at low scales relative to the mass is much better gauged from Fig. 7. The relative changes in the gluon and light quarks for variations in mb are significantly reduced due to the much smaller impact of the beauty contribution to the structure functions from the charge-squared weighting, as can be seen in Figs. 10 and 11, where we show NLO PDFs for mb = 4.25 and 5.25 GeV. At Q2 = 40 GeV2 ∼ 2m2b the relative change 123 Eur. Phys. J. C (2016) 76 :10 Page 7 of 14 10 ZEUS H1 mb = 5.25 GeV mb = 5.0 GeV mb = 4.75 GeV mb = 4.5 GeV mb = 4.25 GeV x σ̃bb, Q2 = 5 GeV2 . 0.010.0010.0001 0.004 0.003 0.002 0.001 0 x σ̃bb, Q2 = 6.5 GeV2 . 0.010.0010.0001 0.006 0.005 0.004 0.003 0.002 0.001 0 x σ̃bb, Q2 = 12 GeV2 . 0.010.001 0.01 0.008 0.006 0.004 0.002 0 x σ̃bb, Q2 = 25 GeV2 . 0.010.001 0.02 0.015 0.01 0.005 0 x σ̃bb, Q2 = 30 GeV2 . 0.10.010.001 0.016 0.012 0.008 0.004 0 x σ̃bb, Q2 = 60 GeV2 . 0.10.010.001 0.025 0.02 0.015 0.01 0.005 0 x σ̃bb, Q2 = 80 GeV2 . 0.10.010.001 0.05 0.04 0.03 0.02 0.01 0 x σ̃bb, Q2 = 160 GeV2 . 0.10.01 0.04 0.03 0.02 0.01 0 x σ̃bb, Q2 = 200 GeV2 . 0.10.01 0.035 0.03 0.025 0.02 0.015 0.01 0.005 0 x σ̃bb, Q2 = 600 GeV2 . 0.10.01 0.03 0.025 0.02 0.015 0.01 0.005 0 x σ̃bb, Q2 = 650 GeV2 . 0.10.01 0.03 0.025 0.02 0.015 0.01 0.005 0 x σ̃bb, Q2 = 2000 GeV2 . 0.1 0.02 0.015 0.01 0.005 0 Fig. 7 The (unshifted) HERA σ̃ (bb̄) data versus x at 12 different values of Q2, namely Q2 = 5, 6.5, . . . , 2000 GeV2; the H1 [12] and ZEUS [13] data are shown as solid circular and clear trian- gular points, respectively. The curves are the NNLO predictions for five different values of mb, namely, in descending order, mb = 4.25, 4.5, 4.75, 5, 5.25 GeV in the beauty distribution for a ∼10 % change in the mass is similar to that for the same type of variation for mc. How- ever, the extent to which this remains at Q2 = 104 GeV2 is much greater than the charm case due to the smaller evolution length. 3 Effect on benchmark cross sections In this section we show the variation with mc and mb for cross sections at the Tevatron, and for 7 and 14 TeV at the LHC. Variations for 8 and 13 TeV will be very similar to 123 10 Page 8 of 14 Eur. Phys. J. C (2016) 76 :10 −20 −10 0 10 20 0.0001 0.001 0.01 0.1 Gluon (NLO), percentage difference at Q2 = 4 GeV2 x MMHT2014 mc = 1.25 GeV mc = 1.55 GeV −20 −10 0 10 20 0.0001 0.001 0.01 0.1 Light quarks (NLO), percentage difference at Q2 = 4 GeV2 x MMHT2014 mc = 1.25 GeV mc = 1.55 GeV −60 −40 −20 0 20 40 60 0.0001 0.001 0.01 0.1 Charm quark (NLO), percentage difference at Q2 = 4 GeV2 x MMHT2014 mc = 1.25 GeV mc = 1.55 GeV Fig. 8 The mc dependence of the gluon, light-quark singlet and charm distributions at NLO for Q2 = 4 GeV2, compared to the standard MMHT2014 distributions with mc = 1.4 GeV and mb = 4.75 GeV those at 7 and 14 TeV, respectively. We calculate the cross sections for W and Z boson, Higgs boson via gluon–gluon fusion and top-quark pair production. To calculate the cross section we use the same procedure as was used in [4, 16]. That is, for W, Z and Higgs production we use the code provided by W.J. Stirling, based on the calculation in [17– 19], and for top pair production we use the procedure and code of [20]. Here our primary aim is not to present definitive predictions or to compare in detail to other PDF sets, as both these results are frequently provided in the literature with very specific choices of codes, scales and parameters which may differ from those used here. Rather, our main objective is to illustrate the relative influence of varying mc and mb for these benchmark processes. We show the predictions for the default MMHT2014 PDFs, with PDF uncertainties, and the relative changes due to changing mc from 1.25 to 1.55 GeV, and mb from 4.25 to 5.25 GeV, i.e. changing the default values by approximately 10 % in each case. The dependence of the benchmark pre- dictions on the value of mc in Tables 3, 4 and 5 reflects the behaviour of the gluon with √ s. The changes in cross sec- tion generally scale linearly in variation of masses away from the default values to a good approximation, although for mb, where the cross section sensitivity to the mass choice is often small, this is less true, and in some cases the cross section is even found to decrease or increase in both directions away from the best fit mass. We begin with the predictions for the W and Z production cross sections. The results at NNLO are shown in Table 3. The PDF uncertainties on the cross sections are 2 % at the Teva- tron and slightly smaller at the LHC—the lower beam energy at the Tevatron meaning the cross sections have more contri- bution from higher x where the PDF uncertainties increase. The mc variation is at most about 0.4 % at the Tevatron and is 0.5–1 % at the LHC, being larger at 14 TeV. The results at NLO are very similar. In Table 4 we show the analogous results for the top-quark pair production cross section. At the Tevatron the PDFs are probed in the region x ≈ 0.4/1.96 ≈ 0.2, and the main pro- duction is from the qq̄ channel. At the LHC the dominant pro- duction at higher energies (and with a proton–proton rather than proton–antiproton collider) is gluon–gluon fusion, with the central x value probed being x ≈ 0.4/7 ≈ 0.06 at 7 TeV, and x ≈ 0.4/14 ≈ 0.03 at 14 TeV. The PDF uncertainties on the cross sections are nearly 3 % at the Tevatron, similar for 7 TeV at the LHC, but a little smaller at 14 TeV as there is less sensitivity to the high-x gluon. The mc variation are slightly less than 1 % at the Tevatron and for 7 TeV at the 123 Eur. Phys. J. C (2016) 76 :10 Page 9 of 14 10 −10 0 10 0.0001 0.001 0.01 0.1 Gluon (NLO), percentage difference at Q2 = 104 GeV2 x MMHT2014 mc = 1.25 GeV mc = 1.55 GeV −10 0 10 0.0001 0.001 0.01 0.1 Light quarks (NLO), percentage difference at Q2 = 104 GeV2 x MMHT2014 mc = 1.25 GeV mc = 1.55 GeV −20 −10 0 10 20 0.0001 0.001 0.01 0.1 Charm quark (NLO), percentage difference at Q2 = 104 GeV2 x MMHT2014 mc = 1.25 GeV mc = 1.55 GeV Fig. 9 The mc dependence of the gluon, light-quark singlet and charm distributions at NLO for Q2 = 104 GeV2, compared to the standard MMHT2014 distributions with mc = 1.4 GeV and mb = 4.75 GeV LHC, but rather lower at 14 TeV since the x probed is near the fixed point for the gluon (see Fig. 9). In Table 5 we show the uncertainties in the rate of Higgs boson production from gluon–gluon fusion. At the Tevatron the dominant x range probed, i.e. x ≈ 0.125/1.96 ≈ 0.06, corresponds to a region where the gluon distribution falls as mc increases and at the LHC where x ≈ 0.01–0.02 at central rapidity the gluon increases as mc increases, though at 7 TeV we are only just below the fixed point. At the Tevatron the resultant uncertainty is ∼0.7 %. At the LHC at 7 TeV it is in the opposite direction but only ∼0.1 %, whereas at 14 TeV it has increased to near 0.5 %. As in [5] we recommend that in order to estimate the total uncertainty due to PDFs and the quark masses it is best to add the variation due to the variation in quark mass in quadrature with the PDF uncertainty, or the PDF+αS uncertainty, if the αS uncertainty is also used. 4 PDFs in three- and four-flavour-number schemes In our default studies we work in a general-mass variable- flavour-number scheme (GM-VFNS) with a maximum of five active flavours. This means that we start at our input scale of Q20 = 1 GeV2 with three active light flavours. At the transition point m2c the charm quark starts evolution and then at m2b the beauty quark also starts evolution. The evo- lution is in terms of massless splitting functions, and at high Q2 the contribution from charm and bottom quarks lose all mass dependence other than that in the boundary conditions at the chosen transition point. The explicit mass dependence is included at lower scales, but falls away like inverse powers as Q2/m2c,b → ∞. We do not currently ever consider the top quark as a parton. We could alternatively keep the information about the heavy quarks only in the coefficient functions, i.e. the heavy quarks would only be generated in the final state. This is called a fixed-flavour-number scheme (FFNS). One example would be where neither charm and beauty exist as partons. This would be a three-flavour FFNS. An alternative would be to turn on charm evolution but never allow beauty to be treated as a parton. This is often called a four-flavour FFNS. We will use this notation, but strictly speaking it is a GM- VFNS with a maximum of four active flavours. One might produce the partons for the three- and four- flavour FFNS by performing global fits in these schemes. 123 10 Page 10 of 14 Eur. Phys. J. C (2016) 76 :10 −10 0 10 0.0001 0.001 0.01 0.1 Gluon (NLO), percentage difference at Q2 = 40 GeV2 x MMHT2014 mb = 4.25 GeV mb = 5.25 GeV −10 0 10 0.0001 0.001 0.01 0.1 Light quarks (NLO), percentage difference at Q2 = 40 GeV2 x MMHT2014 mb = 4.25 GeV mb = 5.25 GeV −75 −50 −25 0 25 50 75 0.0001 0.001 0.01 0.1 Bottom quark (NLO), percentage difference at Q2 = 40 GeV2 x MMHT2014 mb = 4.25 GeV mb = 5.25 GeV Fig. 10 The mb dependence of the gluon, light-quark singlet and charm distributions at NLO for Q2 = 40 GeV2, compared to the standard MMHT2014 distributions with mc = 1.4 GeV and mb = 4.75 GeV However, it was argued in [21] that the fit to structure func- tion data is not optimal in these schemes. Indeed, evidence for this has been provided in [9, 22, 23]. Moreover, much of the data (for example, on inclusive jets and W, Z production at hadron colliders) is not known to NNLO in these schemes, and is very largely at scales where mc,b are relatively very small. So it is clear that the GM-VFNS are more appropri- ate. Hence, in [24] it was decided to make available PDFs in the three- and four-flavour schemes simply by using the input PDFs obtained in the GM-VFNS, but with evolution of the beauty quark, or both the beauty and the charm quark, turned off. This procedure was continued in [5] and is the common choice for PDF groups who fit using a GM-VFNS but make PDFs available with a maximum of three or four active flavours. Hence, here, we continue to make this choice for the MMHT2014 PDFs. We make PDFs available with a maximum of three or four active flavours for the NLO central PDFs and their uncertainty eigenvectors for both the standard choices of α n f,max=5 S (M 2 Z ) of 0.118 and 0.120, and for the NNLO central PDF and the uncertainty eigenvectors for the standard choice of α n f,max=5 S (M 2 Z ) of 0.118. We also provide PDF sets with αS(M 2 Z ) displaced by 0.001 from these default values, so as to assist with the calculation of αs uncertainties in the different flavour schemes. Finally, we make available PDF sets with different values of mc and mb in the different fixed-flavour schemes. By default, when the charm or beauty quark evolution is turned off, we also turn off the contribution of the same quark to the running coupling. This is because most calcu- lations use this convention when these quarks are entirely final state particles. This results in the coupling running more quickly. So if the coupling at Q20 is chosen so that α n f,max=5 S (M 2 Z ) ≈ 0.118, then we find that α n f,max=3 S (M 2 Z ) ≈ 0.105 and α n f,max=4 S (M 2 Z ) ≈ 0.113. There are sometimes cases where a set of PDFs with no beauty quark but with five-flavour running coupling is desired, e.g. [25]. After the publication of [5], PDF sets with this definition were made available. Here we make available PDFs for the central sets together with their eigenvectors with a maximum of four active flavours, but the beauty quark included in the running of the coupling. This type of PDF has also been considered very recently in [26]. The variation of the PDFs defined with a maximum num- ber of three and four flavours, compared to our default of five flavours, is shown at Q2 = 104 GeV2 in Fig. 12 for NNLO 123 Eur. Phys. J. C (2016) 76 :10 Page 11 of 14 10 −10 0 10 0.0001 0.001 0.01 0.1 Gluon (NLO), percentage difference at Q2 = 104 GeV2 x MMHT2014 mb = 4.25 GeV mb = 5.25 GeV −5 0 5 0.0001 0.001 0.01 0.1 Light quarks (NLO), percentage difference at Q2 = 104 GeV2 x MMHT2014 mb = 4.25 GeV mb = 5.25 GeV −10 0 10 0.0001 0.001 0.01 0.1 Bottom quark (NLO), percentage difference at Q2 = 104 GeV2 x MMHT2014 mb = 4.25 GeV mb = 5.25 GeV Fig. 11 The mb dependence of the gluon, light-quark singlet and charm distributions at NLO for Q2 = 104 GeV2, compared to the standard MMHT2014 distributions with mc = 1.4 GeV and mb = 4.75 GeV Table 3 Predictions for W ± and Z cross sections (in nb), including leptonic branching, obtained with the NNLO MMHT2014 parton sets. The PDF uncertainties and mc and mb variations are also shown, where the mc variation corresponds to ±0.15 GeV and the mb variation cor- responds to ±0.5 GeV, i.e. about 10 % in each case σ PDF unc. mc var. mb var. W Tevatron (1.96 TeV) 2.78 +0.0017−0.056 ( +2.0 % −2.0 % ) +0.0017 −0.0086 ( +0.061 % −0.31 % ) −0.00092 −0.0015 ( −0.033 % −0.052 % ) Z Tevatron (1.96 TeV) 0.256 +0.0052−0.0046 ( +2.0 % −1.8 % ) +0.00042 −0.0011 ( +0.16 % −0.43 % ) −0.00029 −0.000016 ( −0.11 % −0.0059 % ) W + LHC (7 TeV) 6.20 +0.103−0.092 ( +1.7 % −1.5 % ) +0.029 −0.040 ( +0.48 % −0.64 % ) +0.0043 −0.014 ( +0.070 % −0.22 % ) W − LHC (7 TeV) 4.31 +0.067−0.076 ( +1.6 % −1.8 % ) +0.019 −0.022 ( +0.44 % −0.51 % ) +0.0059 −0.0091 ( +0.14 % −0.21 % ) Z LHC (7 TeV) 0.964 +0.014−0.013 ( +1.5 % −1.3 % ) +0.0074 −0.0088 ( +0.77 % −0.92 % ) −0.00096 −0.00038 ( −0.10 % −0.039 % ) W + LHC (14 TeV) 12.5 +0.22−0.18 ( +1.8 % −1.4 % ) +0.091 −0.12 ( +0.73 % −0.93 % ) +0.0087 −0.037 ( +0.069 % −0.30 % ) W − LHC (14 TeV) 9.3 +0.15−0.14 ( +1.6 % −1.5 % ) +0.064 −0.075 ( +0.69 % −0.81 % ) +0.012 −0.029 ( +0.13 % −0.31 % ) Z LHC (14 TeV) 2.06 +0.035−0.030 ( +1.7 % −1.5 % ) +0.021 −0.025 ( +1.03 % −1.2 % ) −0.0035 −0.0013 ( −0.17 % −0.062 % ) PDFs. The general form of the differences are discussed in detail in section 4 of [5] and are primarily due to two effects. For fewer active quarks there is less gluon branching, so the gluon is larger if the flavour number is smaller. Also, as Q2 increases the coupling gets smaller for fewer active quarks, so evolution is generally slower, which means par- tons decrease less quickly for large x and grow less quickly at small x. The latter effect dominates for quark evolution, while for the gluon the two effects compete at small x. For the case where the maximum number of flavours is 4, but 123 10 Page 12 of 14 Eur. Phys. J. C (2016) 76 :10 Table 4 Predictions for tt cross sections (in nb), obtained with the NNLO MMHT2014 parton sets. The PDF uncertainties and mc and mb variations are also shown, where the mc variation corresponds to ±0.15 GeV and the mb variation corresponds to ±0.5 GeV σ PDF unc. mc var. mb var. tt Tevatron (1.96 TeV) 7.5 +0.21−0.20 ( +2.8 % −2.7 % ) −0.059 +0.077 ( −0.78 % +1.0 % ) +0.0088 +0.0015 ( +0.12 % +0.20 % ) tt LHC (7 TeV) 176 +3.9−5.5 ( +2.2 % −3.1 % ) −1.1 +1.4 ( −0.60 % +0.77 % ) +0.77 −0.009 ( +0.44 % −0.0051 % ) tt LHC (14 TeV) 970 +16−20 ( +1.6 % −2.1 % ) −3.0 +3.1 ( −0.31 % +0.32 % ) +3.1 −1.7 ( −0.32 % +0.17 % ) Table 5 Predictions for the Higgs boson cross sections (in nb), obtained with the NNLO MMHT 2014 parton sets. The PDF uncertainties and mc and mb variations are also shown, where the mc variation corresponds to ±0.15 GeV and the mb variation corresponds to ±0.5 GeV σ PDF unc. mc var. mb var. Higgs Tevatron (1.96 TeV) 0.87 +0.024−0.030 ( +2.7 % −3.4 % ) −0.0060 +0.0070 ( −0.68 % +0.79 % ) +0.0042 −0.0011 ( +0.48 % −0.13 % ) Higgs LHC (7 TeV) 14.6 +0.21−0.29 ( +1.4 % −2.0 % ) +0.025 −0.019 ( +0.17 % −0.13 % ) +0.049 −0.044 ( +0.34 % −0.30 % ) Higgs LHC (14 TeV) 47.7 +0.63−0.88 ( +1.3 % −1.8 % ) +0.27 −0.22 ( +0.57 % −0.48 % ) +0.16 −0.16 ( +0.34 % −0.33 % ) 0.9 1 1.1 1.2 1.3 0.0001 0.001 0.01 0.1 f3(x, Q2)/f5(x, Q2), Q2 = 104 GeV 2 x Gluon Up Down Strange Charm Anti–up Anti–down 0.9 1 1.1 1.2 1.3 0.0001 0.001 0.01 0.1 f4(x, Q2)/f5(x, Q2), Q2 = 104 GeV 2 x Gluon Up Down Strange Charm Anti–up Anti–down 0.9 1 1.1 1.2 1.3 0.0001 0.001 0.01 0.1 f4,αs(5)(x, Q 2)/f5(x, Q2), Q2 = 104 GeV 2 x Gluon Up Down Strange Charm Anti–up Anti–down Fig. 12 The ratio of the different fixed-flavour PDFs to the standard five flavour PDFs at NNLO and at Q2 = 104 GeV2. The three and four flavour schemes are show in the top left and right plots, while the four flavour scheme with five flavours in the running of αS is shown in the bottom plot 123 Eur. Phys. J. C (2016) 76 :10 Page 13 of 14 10 the coupling has five-flavour evolution, the overwhelming effect is that the gluon is larger—effectively replacing the missing beauty quarks in the momentum sum rule. However, the increase in the gluon is maximal at small x, where the increased coupling compared to the case where we use the four-flavour coupling leading to increased loss of gluons at high x from evolution. 5 Renormalisation schemes At present most PDF fitting groups, including the most recent updates [4, 27–29], use the pole mass definition for the heavy quarks. Hence, we have remained with this definition in our investigation of quark mass dependence in this article. The analyses in [30, 31] use the MS definition, following the developments in [32]. The latter analyses perform their fits in the fixed-flavour-number scheme (FFNS), while all the others groups use a general-mass variable-flavour-number scheme. There is no fundamental obstacle to switching between the two renormalisation schemes using either approach. The mass dependence in a GM-VFNS appears in entirety from the FFNS coefficient functions and in the transition matrix elements which set the boundary conditions for the (mass- less) evolution of the charm and beauty quarks. These are used along with the FFNS coefficient functions to define the GM-VFNS coefficient functions which tend to the massless versions as m2c,b/Q 2 → 0. Under a change in the definition of the quark mass2 mpole = m(μR)(1 + αS(μ2R)d1(μ2R) + · · · ) (2) the coefficient functions and transition matrix elements can be transformed from one mass scheme to the other straight- forwardly, as illustrated in Eq. (8) of [32], and the mass in GM-VFNS defined in the MS renormalisation scheme. However, there is more sensitivity to the definition of the mass in a FFNS at given order than in a GM-VFNS. At LO there is no mass scheme dependence in the same way that there is no renormalisation scheme dependence of any sort. At NLO in the FFNS the variation of the LO O(αS) coeffi- cient function under the change in Eq. (2) leads to a change in the NLO O(αS) coefficient function. Some NLO GM-VFNS definitions (e.g. the SACOT(χ ) [33] and the FONLL-A [34]) only use the FFNS coefficient functions at O(αS). The tran- sition matrix element for heavy-quark evolution in an NLO GM-VFNS is also defined at O(αS) (and indeed is zero with the standard choice μF = mc,b), so neither depend on the mass definition, and the NLO GM-VFNS is independent of the mass scheme [35]. Some NLO GM-VFNS definitions do use the O(α2S) FFNS coefficient functions. Hence, these will contain some 2 Note that d1(μ2R) = 4/3π if μR = m. dependence on the mass scheme. However, in the original TR [36] and then the TR’ [10] schemes this contribution is frozen at Q2 = m2c,b, so becomes relatively very small at high Q2. In the “optimal” TR’ scheme [9], and in the FONLL-B, the dependence falls away like m2c,b/Q 2 (in the former case the whole O(α2S) coefficient function is weighted by m2c,b/Q2, while in the FONLL-B scheme the subtraction means that only the massless limit of the O(α2S) coefficient function remains as m2c,b/Q 2 → 0). Hence, the dependence on the mass scheme is more limited than in the FFNS at NLO, and is particularly small. Indeed, in all but the original TR and TR’ schemes, there is no dependence at high Q2. At NNLO the mass scheme dependence in the FFNS enters in the O(α2S) and O(α3S) coefficient functions. In a GM- VFNS it now enters in the O(α2S) coefficient functions at low scales, and in boundary conditions for evolution, which gives effects which persist to all scales. If the GM-VFNS uses the O(α3S) coefficient functions these will also give mass scheme dependent effects at low Q2. However, the expressions for the O(α3S) coefficient functions are themselves still approx- imations [37]. Hence, at present it does not seem too important whether the pole mass or MS renormalisation scheme is used in a GM-VFNS (indeed in [27] the pole mass scheme is used, but the MS values for the masses are taken). Nevertheless, in the future it is probably ideal to settle on the MS mass, since the value of this is quite precisely determined in many experiments, which is not true of the pole mass. At the same time it will also be desirable for different PDF groups to agree on a common value of mc and mb (there is no agreement at present). 6 Conclusions The main purpose of this article is to present and make available PDF sets in the framework used to produce the MMHT2014 PDFs, but with differing values of the charm and beauty quark mass. We do not make a determination of the optimum values of these masses, but we do investigate and note the effect the mass variation has on the quality of the fits to the data, concentrating on the HERA cross section data with charm or beauty in the final state. We note that for both the charm and beauty quarks a lower mass than our default values of mc = 1.4 GeV and mb = 4.75 GeV is pre- ferred, although these are roughly the values of pole masses one would expect from conversion from the values measured in the MS scheme. This suggests that in the future it may be better to use the MS definition, though this is currently not the practice in global fits using a GM-VFNS—perhaps because, as we discuss, the mass scheme dependence has less effect in these schemes than for the FFNS. We also make PDFs avail- 123 10 Page 14 of 14 Eur. Phys. J. C (2016) 76 :10 able with a maximum of three or four active quark flavours. The PDF sets obtained for different quark masses and for different active quark flavours can be found at [38] and will be available from the LHAPDF library [39]. We investigate the variation of the PDFs and the predicted cross sections for standard processes at the LHC (and Teva- tron) corresponding to these variations in heavy-quark mass. For reasonable variations of mc the effects are small, but not insignificant, compared to PDF uncertainties. For vari- ations in mb the effect is smaller, and largely insignificant, except for the beauty distribution itself, which can vary more than its uncertainty at a fixed value of mb; see, in particu- lar Fig. 10. Hence, currently the uncertainties on PDFs due to quark masses are not hugely important, but need to be improved in the future for very high precision predictions at hadron colliders. Acknowledgments We particularly thank W. J. Stirling and G. Watt for numerous discussions on PDFs and for previous work without which this study would not be possible. We would like to thank Mandy Cooper- Sarkar, Albert de Roeck, Stefano Forte, Joey Huston, Pavel Nadolsky and Juan Rojo for various discussions on the relation between PDFs and quark masses. We would like to thank A. Geiser, A. Gizhko and K. Wichmann for help regarding the treatment of uncertainties for ZEUS beauty cross sections. This work is supported partly by the London Cen- tre for Terauniverse Studies (LCTS), using funding from the European Research Council via the Advanced Investigator Grant 267352. RST would also like to thank the IPPP, Durham, for the award of a Research Associateship held while most of this work was performed. We thank the Science and Technology Facilities Council (STFC) for support via Grant Awards ST/J000515/1 and ST/L000377/1. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. References 1. A.D. Martin, W.J. Stirling, R.S. Thorne, G. Watt, Eur. Phys. J. C 63, 189 (2009). arXiv:0901.0002 2. H1 and ZEUS Collaboration, F. Aaron et al., JHEP 1001, 109 (2010). arXiv:0911.0884 3. H1 Collaboration, ZEUS Collaboration, H. Abramowicz et al., Eur. Phys. J. C 73, 2311 (2013). arXiv:1211.1182 4. L. Harland-Lang, A. Martin, P. Motylinski, R. Thorne, Eur. Phys. J. C 75, 204 (2015). arXiv:1412.3989 5. A. Martin, W. Stirling, R. Thorne, G. Watt, Eur. Phys. J. C 70, 51 (2010). arXiv:1007.2624 6. Particle Data Group, K. Olive et al., Chin. Phys. C 38, 090001 (2014) 7. C.W. Bauer, Z. Ligeti, M. Luke, A.V. Manohar, M. Trott, Phys. Rev. D 70, 094017 (2004). arXiv:hep-ph/0408002 8. A.H. Hoang, A.V. Manohar, Phys. Lett. B 633, 526 (2006). arXiv:hep-ph/0509195 9. R.S. Thorne, Phys. Rev. D 86, 074017 (2012). arXiv:1201.6180 10. R.S. Thorne, Phys. Rev. D 73, 054019 (2006). arXiv:hep-ph/0601245 11. S. Alekhin, J. Blümlein, K. Daum, K. Lipka, S. Moch, Phys. Lett. B 720, 172 (2013). arXiv:1212.2355 12. H1, F.D. Aaron et al., Eur. Phys. J. C 65, 89 (2010). arXiv:0907.2643 13. ZEUS, H. Abramowicz et al., JHEP 09, 127 (2014). arXiv:1405.6915 14. CMS Collaboration, S. Chatrchyan et al., JHEP 1312, 030 (2013). arXiv:1310.7291 15. M. Buza, Y. Matiounine, J. Smith, W.L. van Neerven, Eur. Phys. J. C 1, 301 (1998). arXiv:hep-ph/9612398 16. L.A. Harland-Lang, A.D. Martin, P. Motylinski, R. S. Thorne, Eur. Phys. J. C 75, 435 (2015). doi:10.1140/epjc/s10052-015-3630-3 17. R. Hamberg, W. van Neerven, T. Matsuura, Nucl. Phys. B 359, 343 (1991) 18. R.V. Harlander, W.B. Kilgore, Phys. Rev. Lett. 88, 201801 (2002). arXiv:hep-ph/0201206 19. A. Djouadi, M. Spira, P. Zerwas, Phys. Lett. B 264, 440 (1991) 20. M. Czakon, P. Fiedler, A. Mitov, Phys. Rev. Lett. 110, 252004 (2013). arXiv:1303.6254 21. R.S. Thorne, W.K. Tung, in PQCD Formulations with Heavy Quark Masses and Global Analysis, Proceedings, Workshop: HERA and the LHC Workshop Series on the Implications of HERA for LHC Physics (2008). arXiv:0809.0714 22. R. Thorne, Eur. Phys. J. C 74, 2014 (2958). arXiv:1402.3536 23. The NNPDF Collaboration, R.D. Ball et al., Phys. Lett. B 723, 330 (2013). arXiv:1303.1189 24. A.D. Martin, W.J. Stirling, R.S. Thorne, Phys. Lett. B 636, 259 (2006). arXiv:hep-ph/0603143 25. F. Cascioli, P. Maierhöfer, N. Moretti, S. Pozzorini, F. Siegert, Phys. Lett. B 734, 210 (2014). arXiv:1309.5912 26. V. Bertone, S. Carrazza, J. Rojo, in Doped Parton Distributions, 27th Rencontres de Blois on Particle Physics and Cosmology Blois, France (2015). arXiv:1509.04022 27. The NNPDF Collaboration, R.D. Ball et al. (2014). arXiv:1410.8849 28. S. Dulat et al. (2015). arXiv:1506.07443 29. ZEUS, H1, H. Abramowicz et al. (2015). arXiv:1506.06042 30. S. Alekhin, J. Bluemlein, S. Moch, Phys. Rev. D 89, 054028 (2014). arXiv:1310.3059 31. P. Jimenez-Delgado, E. Reya, Phys. Rev. D 89, 074049 (2014). arXiv:1403.1852 32. S. Alekhin, S. Moch, Phys. Lett. B 699, 345 (2011). arXiv:1011.5790 33. W.-K. Tung, S. Kretzer, C. Schmidt, J. Phys. G 28, 983 (2002). arXiv:hep-ph/0110247 34. S. Forte, E. Laenen, P. Nason, J. Rojo, Nucl. Phys. B 834, 116 (2010). arXiv:1001.2312 35. V. Bertone, J. Rojo (2012). arXiv:1212.0741. [AIP Conf. Proc. 1523, 51 (2012)] 36. R.S. Thorne, R.G. Roberts, Phys. Rev. D 57, 6871 (1998). arXiv:hep-ph/9709442 37. H. Kawamura, N.A. Lo Presti, S. Moch, A. Vogt, Nucl. Phys. B 864, 399 (2012). arXiv:1205.5727 38. http://www.hep.ucl.ac.uk/mmht/. Accessed 4 Jan 2016 39. http://lhapdf.hepforge.org. Accessed 4 Jan 2016 123 http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ http://arxiv.org/abs/0901.0002 http://arxiv.org/abs/0911.0884 http://arxiv.org/abs/1211.1182 http://arxiv.org/abs/1412.3989 http://arxiv.org/abs/1007.2624 http://arxiv.org/abs/hep-ph/0408002 http://arxiv.org/abs/hep-ph/0509195 http://arxiv.org/abs/1201.6180 http://arxiv.org/abs/hep-ph/0601245 http://arxiv.org/abs/1212.2355 http://arxiv.org/abs/0907.2643 http://arxiv.org/abs/1405.6915 http://arxiv.org/abs/1310.7291 http://arxiv.org/abs/hep-ph/9612398 http://dx.doi.org/10.1140/epjc/s10052-015-3630-3 http://arxiv.org/abs/hep-ph/0201206 http://arxiv.org/abs/1303.6254 http://arxiv.org/abs/0809.0714 http://arxiv.org/abs/1402.3536 http://arxiv.org/abs/1303.1189 http://arxiv.org/abs/hep-ph/0603143 http://arxiv.org/abs/1309.5912 http://arxiv.org/abs/1509.04022 http://arxiv.org/abs/1410.8849 http://arxiv.org/abs/1506.07443 http://arxiv.org/abs/1506.06042 http://arxiv.org/abs/1310.3059 http://arxiv.org/abs/1403.1852 http://arxiv.org/abs/1011.5790 http://arxiv.org/abs/hep-ph/0110247 http://arxiv.org/abs/1001.2312 http://arxiv.org/abs/1212.0741 http://arxiv.org/abs/hep-ph/9709442 http://arxiv.org/abs/1205.5727 http://www.hep.ucl.ac.uk/mmht/ http://lhapdf.hepforge.org