Phenomenology of single-inclusive jet production with jet radius and threshold resummation

Phenomenology of single-inclusive jet production with jet radius and threshold resummation

Xiaohui Liu xiliu@bnu.edu.cn Center of Advanced Quantum Studies, Department of Physics, Beijing Normal University, Beijing 100875, China    Sven-Olaf Moch sven-olaf.moch@desy.de II. Institut für Theoretische Physik, Universität Hamburg, Luruper Chaussee 149, D-22761 Hamburg, Germany    Felix Ringer fmringer@lbl.gov Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
Abstract

We perform a detailed study of inclusive jet production cross sections at the LHC and compare the QCD theory predictions based on the recently developed formalism for threshold and jet radius joint resummation at next-to-leading logarithmic accuracy to inclusive jet data collected by the CMS collaboration at and  TeV. We compute the cross sections at next-to-leading order in QCD with and without the joint resummation for different choices of jet radii and observe that the joint resummation leads to crucial improvements in the description of the data. Comprehensive studies with different parton distribution functions demonstrate the necessity of considering the joint resummation in fits of those functions based on the LHC jet data.

DESY 18-007

I Introduction

Long term persistence in achieving higher order calculations in perturbative Quantum Chromodynamics (pQCD) paves the way to the precision frontier at the Large Hadron Collider (LHC). With many Standard Model processes now being measured with an impressive accuracy at the LHC, theoretical predictions beyond next-to-leading order (NLO) in pQCD, nowadays considered the standard in phenomenological analyses, are often mandatory. During the past three years, there have been a burst of publications on complete next-to-next-to-leading order (NNLO) calculations for various hadro-production processes involving jets Boughezal et al. (2015a, b, c); Chen et al. (2015); Gehrmann-De Ridder et al. (2016); Boughezal et al. (2016a); Currie et al. (2017a, b); Campbell et al. (2017). The list of those processes includes the hadro-production of gauge bosons +jet () as well as single-inclusive jets and dijets, but it is limited to reactions at Born level due to the enormous computational complexity at NNLO. In particular the calculations for +jet production have already been shown to greatly improve the description of the available LHC data Boughezal et al. (2016b); Lindert et al. (2017); Boughezal et al. (2017).

For the hadro-production of jets at the LHC the experimental collaborations have provided very precise data for the single-inclusive jet production cross sections at all collider energies and differential in the jet transverse momentum and the rapidity . Specifically, ALICE Abelev et al. (2013), ATLAS Aad et al. (2013) and CMS Khachatryan et al. (2017a) have collected data at  TeV and ATLAS and CMS at  TeV Aad et al. (2015); Chatrchyan et al. (2014),  TeV Aaboud et al. (2017a); Khachatryan et al. (2017b) and  TeV Aaboud et al. (2017b); ATL (2017); Khachatryan et al. (2016). These data allow for important consistency tests of pQCD as well as a precise extraction of the value of the strong coupling constant  Britzger et al. (2017) and they provide very valuable constraints on parton distribution functions (PDFs) which govern the parton luminosity of the colliding initial protons Nocera and Ubiali (2017); Harland-Lang et al. (2017).

In order to fully utilize the available data, a precise understanding of the corresponding theoretical calculations within pQCD is very important. The current accuracy for fixed order pQCD predictions is NNLO where the coefficient is known in the leading-color approximation Currie et al. (2017a), i.e. for large values of for a general SU gauge group. Any additional corrections are parametrically suppressed as , so that the results of Currie et al. (2017a) are supposed to approximate the full NNLO calculation very well. Preliminary comparisons of those NNLO results with some of the LHC data, however, have not been entirely satisfactory. Refs. Currie et al. (2017c); ATL (2017) have shown that the agreement between theory and data heavily depends on the choices for the renormalization and factorization scales and . Moreover, for some natural scale choices, such as identifying and with the transverse momentum of the leading jet in the event, i.e. , the theory description of the data at NNLO deteriorates compared to NLO. This situation implies the existence of potentially large higher order corrections beyond fixed NNLO.

Improvements beyond fixed order in pQCD are possible by supplementing the fixed order calculations with resummation results where dominant classes of logarithmic corrections are summed up to all orders in the strong coupling constant. Recently, a joint resummation framework was developed Liu et al. (2017) that allows to resum both threshold and jet radius logarithms simultaneously. Threshold logarithms appear in the partonic cross section at -th order as where and . Here, is the invariant mass of the partonic system recoiling against the observed jet and is the partonic center-of-mass energy Kidonakis et al. (1998). Since these logarithms are integrated over the specified parton kinematics together with the steeply falling parton luminosity, threshold logarithms can dominate the entire cross section in a wide kinematic range. Instead, the jet radius is an external quantity and the dependence of the cross section is single-logarithmic with instead of double-logarithmic Dasgupta et al. (2015); Kang et al. (2016); Dai et al. (2016). The framework developed in Liu et al. (2017) addresses both these logarithmic corrections on the same footing and it was shown that numerically the threshold and the jet radius logarithmic terms do account for the dominant bulk of the NLO corrections. The explicit resummation of these logarithms to next-to-leading logarithmic (NLL) accuracy was also realized within the joint resummation framework derived in Liu et al. (2017) and the subsequent matching to fixed order NLO results leads to theory predictions at the combined NLO NLL accuracy. The approach of Liu et al. (2017), however, is not limited to this logarithmic accuracy and the framework is ready for a systematic extension to the next-to-next-to leading logarithmic (NNLL) accuracy which may then be matched to the available fixed order NNLO results to achieve a combined accuracy of NNLO NNLL. We leave the extension to NNLL for future work and instead focus here on the phenomenological results at NLO NLL accuracy.

In general, one expects competing effects from threshold and small- resummation. As it was observed in Kidonakis and Owens (2001); de Florian and Vogelsang (2007); Kumar and Moch (2014); de Florian et al. (2014) threshold resummation leads to an enhancement whereas small- resummation alone leads to a decrease of the cross section Dasgupta et al. (2016); Kang et al. (2016), see also Kang et al. (2018) for studies on jet angularities. Depending on the non-trivial interplay within the joint resummation framework, one or the other effect will dominate. For certain kinematics and values of , the two effects may even largely cancel out. In order to obtain a good understanding of the convergence of the perturbative series expansion it is important to disentangle these two effects. A closely related issue is the dependence of the fixed order and the resummed calculations on the renormalization and factorization scales and , collectively denoted by in the following. As it was pointed out in Currie et al. (2017c); ATL (2017), the fixed order results change significantly depending on whether the hard scale is chosen as of the individual jet or as the transverse momentum of the leading jet in the event. At the same time, the residual scale dependence is very small and even vanishes for some kinematic configurations. In Dasgupta et al. (2015, 2016); Kang et al. (2016), it was argued that this is generally an artifact of results at fixed order in perturbation theory. Here, we address this issue within the joint resummation formalism.

In this work, we provide a detailed comparison with LHC data and find that the inclusion of the resummation generally yields a much better description of those data. In addition, our studies highlight possible improvements that can be obtained by using a resummed calculation in fits of PDFs. The constraints from inclusive jet data on PDF fits are most significant for the gluon PDF in the large- region. In this endpoint region the cross sections from which PDFs are extracted can be subject to large logarithmic corrections that need to be taken into account to all orders. Improvements in the precision of the extracted PDFs eventually have direct impact on all PDF sensitive analyses at the LHC and recent progress on PDFs in the large- region has been made in Accardi et al. (2015); Bonvini et al. (2015).

The remainder of this work is organized as follows. In section II, we briefly review the theoretical framework of Liu et al. (2017). In section III, we present detailed phenomenological studies of the resummation effect and the scale dependence of the resummed cross section. We study cross section ratios for different jet radii to discriminate the predictive power of the NLL NLO and the NLO results. Finally, we present a comprehensive comparison to the inclusive jet data from the LHC together with the impact of different PDF sets. We conclude in section IV with a summary and an outlook.

Ii Theoretical framework

First we review the theoretical formalism which allows us to achieve the threshold and small- joint resummation used in this work. The resummation is based on the factorization theorem Liu et al. (2017) developed within the framework of the Soft Collinear Effective Theory Bauer et al. (2000, 2001); Bauer and Stewart (2001); Bauer et al. (2002), in which the single-inclusive jet cross section with jet transverse momentum , jet rapidity and a small anti- Cacciari et al. (2008) jet radius near the partonic threshold can be written as

(1)

where the partonic cross sections are further factorized as

(2)

In Eq. (1) the PDFs are denoted by which are evaluated at the momentum fractions and , where , and is the hadronic center-of-mass energy. The sum runs over all partonic channels initiating the subprocesses and runs over the collinear splitting history. The associated angular integrals are denoted by ‘’ Becher et al. (2016) to resum non-global logarithms Dasgupta and Salam (2001); Caron-Huot (2015); Becher et al. (2016); Larkoski et al. (2015, 2016); Neill (2017). Besides the jet , the partonic cross sections depend on the partonic kinematic variables , and with and . Here, are the momenta of the incoming partons and is the momentum of the parton which initiates the signal-jet.

The hard scattering functions in Eq. (II) are denoted by which are available to two loops Broggio et al. (2014). The inclusive jet function is also known to order  Becher and Neubert (2006a); Becher and Bell (2011) and the NLO jet function can be extracted from Ellis et al. (2010); Liu and Petriello (2013). The global soft function and the soft collinear Becher et al. (2016); Chien et al. (2015) function have been derived to NLO in Dai et al. (2017); Liu et al. (2017). The global soft function and the soft collinear function can be readily calculated to two loops following Becher et al. (2012) and Kelley et al. (2011); Boughezal et al. (2015d). All the functions are evolved from their natural scales to the common scale according to their renormalization group equations in order to obtain the NLL resummation used in this work.

The factorization formalism in Eq. (II) holds in the threshold regime in which and . To extend the region of validity, we combine the NLL resummed results with the NLO predictions using an additive matching procedure and define

(3)

Here, the logarithmically enhanced contributions at NLO are obtained within the resummation framework and denoted by . They are subtracted from the full NLO calculation and replaced by the NLL resummed results . For the phenomenological studies presented in the next section, we use as a default scale choice the leading jet transverse momentum for the fixed NLO calculations Gao et al. (2013). We vary the scales around the central scale up and down by a factor of two and take the maximal deviation as our NLO scale uncertainties. For the resummed results, we make the central scale choices , for the hard and the signal-jet functions, respectively, and we set , see also Becher and Neubert (2006b); Becher and Schwartz (2010). The other scales are determined in the seesaw way: and for the global soft and the soft collinear functions, respectively. Our uncertainty estimates are obtained by varying , , and independently by a factor of two around their central values while keeping the seesaw relations for and in terms of . The final scale uncertainty is obtained by taking the envelope of the scale variations.

[GeV]
[GeV]
Table 1: Cross section ratios of Eq. (4) for different jet radii at  TeV (top) and 13 TeV (bottom) for selected values of the the signal jet using the MMHT PDF set Harland-Lang et al. (2015) at NLO.
[GeV]
[GeV]
Table 2: Ratios of Eq. (5) for the cross sections at NLO and NLO NLL accuracy denoted by and , respectively, at  TeV (top) and 13 TeV (bottom) for selected values of the the signal jet using the MMHT PDF set Harland-Lang et al. (2015) at NLO. The results at  TeV include NP correction factors which are taken from Chatrchyan et al. (2014).

Iii Phenomenology

We start by studying the overall numerical impact of the joint threshold and small- resummation. We then continue by analyzing the scale dependence of the resummed cross section and provide a detailed comparison to LHC data. Finally, we study in detail the impact of different PDF sets. The two single-inclusive jet data sets

from CMS that we are comparing to throughout this section were taken at  TeV Chatrchyan et al. (2014) and at 13 TeV Khachatryan et al. (2016). For the  TeV data set, the jets were reconstructed using two different values of the jet radius, and covering a rapidity range of . Instead, for the  TeV data set, the jet radius parameters were chosen as and covering . For both data sets, the jets were reconstructed using the anti- algorithm Cacciari et al. (2008) and the transverse momentum of the identified jets ranges up to  TeV.

Figure 1: Cross section ratios of Eq. (4) with different jet radii at  TeV (left) and 13 TeV (right) using the MMHT PDF set Harland-Lang et al. (2015) at NLO.
Figure 2: Ratios of Eq. (5) with at  TeV for the cross sections at NLO NLL (left) and NLO (right) accuracy using the MMHT PDF set Harland-Lang et al. (2015) at NLO with NP correction factors which are taken from Chatrchyan et al. (2014). The red dots indicate the single-inclusive jet data for from CMS collected at  TeV with  Chatrchyan et al. (2014).
Figure 3: Ratios of Eq. (5) with at  TeV for the cross sections at NLO NLL (left) and NLO (right) accuracy using the MMHT PDF set Harland-Lang et al. (2015) at NLO.
Figure 4: The ratio for the CMS data  TeV of Chatrchyan et al. (2014) with (left) and (right) to the theoretical results at NLO (black) and at NLO NLL (red) accuracy using the MMHT PDF set Harland-Lang et al. (2015) at NLO.
Figure 5: Same as Fig. 4 for the CMS data of Khachatryan et al. (2016) at  TeV with (left) and (right).

iii.1 Numerical impact of the joint resummation

In Fig. 1, we plot the ratio

(4)

of the NLO NLL and the NLO cross sections for different jet radii as a function of the signal-jet and at both  TeV (left) and TeV (right). Results for selected values of are listed in Tab. 1. For all numerical calculations in this section we use the MMHT14 PDFs Harland-Lang et al. (2015) at NLO as an example. We find that for a large range of the jet , the joint resummation leads to a decrease of the NLO cross section. The effect is more pronounced for smaller values of , where the impact of the resummation becomes more noticeable and leads to a significant decrease of the cross section. For larger values of the jet , the threshold enhancement will compensate the contributions and can eventually result in an enhancement. This effect is most clearly illustrated in Fig. 1 and Tab. 1 for for around . A similar trend has been observed in recent work on threshold summation with a parton shower event generator including quantum interference Nagy and Soper (2017).

To illustrate the dependence of the cross section on the jet radius parameter we plot in Figs. 2 and 3 for both, fixed order NLO and NLO NLL resummed predictions, the respective ratios

(5)

as a function of the signal-jet and for the energies and  TeV. Nonperturbative (NP) correction factors are taken from Chatrchyan et al. (2014) and results for selected values of are listed in Tab. 2. The predicted dependence of the cross section on the jet radius parameter tends to be more pronounced when the effects of the joint resummation are considered. For the smaller value of , the ratio with resummation is smaller than the corresponding fixed order result , while the opposite trend is observed for larger radii , cf. Tab. 2. We emphasize that the ratios are quite insensitive to the chosen input PDFs. This leads to precise predictions of the NLL NLO and the NLO calculations which are experimentally well testable. In Fig. 2, we have performed a comparison of predictions with the inclusive jet data for shown in red dots measured by CMS Chatrchyan et al. (2014) at the LHC with at  TeV. The error bars represent the experimental uncertainties which are small due to the cancellation of systematic errors in the ratio. After considering the NP effects, the NLO NLL resummed result in Fig. 2 (left) agrees very well with the measurement while the fixed order NLO prediction in Fig. 2 (right) fails to describe the data for the entire range of jet considered. We note that those CMS data for are even larger than the fixed order NLO prediction , i.e. a situation which intuitively should be reversed.

As a further illustration of the resummation effects we compare the predicted central values at fixed NLO and at NLO NLL accuracy in Eq. (3) with the LHC data Chatrchyan et al. (2014); Khachatryan et al. (2016), collected at and  TeV, respectively. In Figs. 4 and 5, the LHC data in the rapidity bin is normalized to the theoretical predictions and displayed as a function of the signal-jet . For  TeV the NP effects have been included in the predictions in Fig. 4. It is clearly visible how the resummed predictions lead to an increase of the ratio of cross sections compared to the NLO result for all choices of jet radii, so that the NLO NLL results of Eq. (3) are in perfect agreement with the  TeV data Chatrchyan et al. (2014). For the  TeV data Khachatryan et al. (2016) with the choice of , the resummation improves the theory description as well, whereas the data for slightly overshoots the theory predictions. However, the experimental uncertainties of those data sets are still relatively large.

The observations presented here do neither depend significantly on the scales chosen as the leading jet transverse momentum nor on the PDFs. This will be quantified in detail in the following sections.

iii.2 Scale dependence

In this section we analyze in detail the scale dependence of the jointly resummed single-inclusive jet cross section. It is instructive to compare the obtained scale dependence to the case where only the logarithms in the jet size parameter are resummed Kang et al. (2016). In Fig. 6, we show the residual scale uncertainty of the jointly resummed cross section normalized to NLO. The scale band is obtained as discussed in the section II above. In addition, we show the cross section where only terms are resummed. In this case the scale band is obtained by varying only the hard scale and the jet scale by factors of two around their canonical choices. One observes a significant reduction of the residual scale dependence once also threshold resummation is taken into account. This observation holds true even though for the small- resummed calculation there are only two scales that are varied in order to estimate the QCD uncertainty whereas there are three separate scales that are all varied independently for the jointly resummed result. One also notices that the reduction of the scale uncertainty gets more pronounced at higher jet transverse momenta where threshold resummation is more relevant. In fact this behavior is generally expected for threshold resummed calculations and has been analyzed in more detail before in many instances, see for example the studies for Higgs boson hadro-production Catani et al. (2003); Ahrens et al. (2009). When approaching large , the joint resummation surpasses the small- resummed cross section due to the threshold enhancement.

An important caveat here is that the resummed calculation of Kang et al. (2016) can currently be performed only with the scale choice whereas in the threshold limit we always have . This difference is most relevant at small values of the jet transverse momentum and likely explains the difference of the central values of the two curves at small . On the other hand, it is interesting to note that for both scale choices, the resummation consistently leads to a suppression relative to the respective NLO calculation. In addition, the jointly resummed calculation is matched and normalized to the full NLO. Instead, the resummed calculation is using the narrow jet approximation. However, the differences are of order which are negligible for  Jager et al. (2004); Mukherjee and Vogelsang (2012).

Figure 6: Comparison of the scale dependence of the jointly resummed cross section (red) and the case where only small- logarithms are resummed to all orders (blue) Kang et al. (2016). As an example, we use the CT14 PDF set Dulat et al. (2015) at NLO with and both results are normalized to NLO. See text for more details.

iii.3 Comparison to LHC data

Now we move on to the comparison of the theoretical predictions with the CMS inclusive jet analyses for both and  TeV Chatrchyan et al. (2014); Khachatryan et al. (2016). Other data sets, such as those by CMS collected at  TeV Khachatryan et al. (2017a) and the one ATLAS at  TeV ATL (2017) have already been considered in Liu et al. (2017).

We start with  TeV following the CMS analysis Chatrchyan et al. (2014) and focus on the cross section data with the anti- jet radius , which we bin into 4 different rapidity regions: , , and . For each rapidity bin, we present the pQCD predictions at NLO NLL and at NLO accuracy based on the CT10 PDFs Lai et al. (2010) at NLO as in the original CMS analysis Chatrchyan et al. (2014).

Figs. 7 and 8 show the ratio of the CMS data to the theoretical predictions, that is for both NLO NLL and NLO accuracy. For both cases, also the NP effects as provided by CMS Chatrchyan et al. (2014) have been included in the perturbative calculations to convert the predictions from the parton level to the particle level. The yellow bands in Figs. 7 and 8 indicate the theoretical uncertainties from scale variations obtained as discussed in the previous section with the hard scale chosen as . The solid brown lines on the other hand indicate the experimental systematic errors, whereas the error bars on the data represent the experimental statistical errors Chatrchyan et al. (2014).

Figure 7: The ratio for the CMS data collected at  TeV Chatrchyan et al. (2014) with to the theoretical results using the CT10 PDF set Lai et al. (2010) at NLO. The error bars represent the experimental statistical errors and the solid (brown) lines the systematic ones. The band (yellow) indicates theoretical scale uncertainties. The NP corrections from Chatrchyan et al. (2014) have been included.
Figure 8: Same as Fig. 7 for the ratio .

In Fig. 7 we observe very good agreement with the data in all rapidity regions for the NLO NLL predictions where the NLL joint resummation is taken into account. In the high- region, the NLO NLL calculations still somewhat overestimate the CMS data. However, this can be further improved by switching from the CT10 PDFs Lai et al. (2010) to more recent PDF sets, as we will detail in the next section. In contrast, the NLO results in Fig. 8 are consistently larger than the inclusive jet data by an amount of in all rapidity bins. Thus, all predictions are lying along the lower boundary of the systematic errors (brown lines) in Fig. 8, except for the high tail region of the rapidity bin .

The theoretical uncertainty of the NLL NLO predictions in Fig. 7 is still large and comparable with the experimental errors. However, this can be reduced further in the future with the help improved accuracy for the resummation, i.e., upon resumming the relevant logarithms to NNLL accuracy together with matching to the available NNLO calculations.

Next we study the inclusive jet production cross section with at  TeV. The results are shown in Figs. 9 and 10, in which the ratio of the CMS data Khachatryan et al. (2016) to the cross sections and are displayed, respectively. Here, we have applied the CT14 PDF set Dulat et al. (2015) at NLO for both predictions and we note that the NP and the electroweak effects have not been included in this analysis. Again, the yellow band in Figs. 9 and 10 represents the theoretical scale uncertainties whereas the experimental systematic and statistical errors are shown as solid brown lines and the error bars, respectively. At present, the published CMS data at  TeV have larger statistical errors, since they are based on data samples corresponding to a relativley small integrated luminosities of 71 and 44 inverse picobarns, whereas the TeV data Chatrchyan et al. (2014) discussed above correspond to an integrated luminosity of 5.0 inverse femtobarns.

As shown in Fig. 9, the NLO NLL calculation leads to a good agreement of the ratio with unity in the region of central rapidities, but slightly overshoots it in the rapidity bin , although still being compatible within the errors. On the contrary, in Fig. 10 the ratio based on the NLO predictions systematically undershoots unity in the rapidity region , but it is still compatible within the quoted uncertainties. Better consistency of the NLO results with the CMS data is only observed in the rapidity region .

Figure 9: The ratio for CMS data collected at  TeV Khachatryan et al. (2016) with to the theoretical results using the PDFs of CT14 Dulat et al. (2015) at NLO. The error bars represent the experimental statistical errors and the solid (brown) lines the systematic ones. The band (yellow) indicates theoretical scale uncertainties.
Figure 10: Same as Fig. 9 for the ratio .

iii.4 Impact of different PDFs

We now proceed to study the sensitivity of different choices of PDFs in predicting the inclusive jet distributions. We benchmark our study using the CMS data at  TeV since the experimental errors of those data are relatively small. Besides the CT10 PDFs used above, we also consider the following alternative PDF extractions which are available in the literature to NLO and NNLO accuracy in pQCD: ABMP16 Alekhin et al. (2018, 2017), CT14 Dulat et al. (2015), HERAPDF2.0 Abramowicz et al. (2015), MMHT2014 Harland-Lang et al. (2015) and NNPDF3.1 Ball et al. (2017). In addition, we use the PDF set of Bonvini et al. (2015) obtained within the framework of NNPDF by fitting only data for the Drell-Yan (DY) process, deep-inelastic scattering (DIS) and top-quark hadro-production but including threshold resummation in all theory predictions for the hard scattering.

Figure 11: Cross sections at NLL NLO accuracy with using the central PDF sets HERAPDF2.0 Abramowicz et al. (2015), MMHT2014 Harland-Lang et al. (2015) and NNPDF3.1 Ball et al. (2017) normalized to the one with CT10 PDFs Lai et al. (2010) at NLO. The dashed lines (blue) indicate the PDF uncertainties for the NNPDF3.1 set. The CMS data collected at TeV Chatrchyan et al. (2014) with with their experimental statistical errors are displayed as dots (red).
Figure 12: Same as Fig. 11 for the PDF sets ABMP16 Alekhin et al. (2018), CT14 Dulat et al. (2015) and NNPDF3.1 Ball et al. (2017).
Figure 13: Same as Fig. 11 for cross sections at NLO in pQCD.
Figure 14: Same as Fig. 12 for cross sections at NLO in pQCD.
Figure 15: Same as Fig. 11 with the NNLO variant taken for all PDF sets ABMP16 Alekhin et al. (2017), CT14 Dulat et al. (2015), HERAPDF2.0 Abramowicz et al. (2015), MMHT2014 Harland-Lang et al. (2015) and NNPDF3.1 Ball et al. (2017).

The NLO NLL studies are presented in Figs. 11 and 12, and the NLO ones in Figs. 13 and 14, respectively. To maintain consistency, the NLO variants for all PDFs are used here, the value of the strong coupling is taken as provided by the respective PDF sets and the predictions are normalized to the one with CT10 PDFs at NLO to allow for comparisons with the CMS analysis Chatrchyan et al. (2014). In Figs. 11 and 12 we see that the predictions based on the NLL joint resummation achieve excellent agreement with the experimental data for the PDF sets CT14, MMHT2014 and NNPDF3.1, while the PDF set ABMP16 slightly undershoots the data for large values of . Only in the high- region of the bin with rapidities , the data are in excess of the theoretical predictions, but the discrepancy is still within the experimental error, see Fig. 7. The predictions with HERAPDF2.0 show a somewhat different trend. They are lower than the data in most of the region and only tend to agree with the data in the high- regime around  TeV. To illustrate the robustness of the observations further, we also plot the uncertainties for the NNPDF3.1 PDFs, which turn out to be quite small compared with the theoretical error from the scale uncertainty, see again Fig. 7. The uncertainties for other PDFs are of similar size.

The situation deteriorates substantially when only the NLO corrections are taken into account. For at NLO, as shown in Figs. 13 and 14, none of the PDFs do a good job in describing the CMS inclusive jet data. The predictions with ABMP16, CT14, MMHT2014 and NNPDF3.1 are all higher than the data in all rapidity bins in most of the jet regions. For  GeV those predictions tend to agree with or slightly undershoot the data for rapidities . For HERAPDF2.0 the NLO predictions are within the experimental errors of the inclusive jet data, but the shape of the dependence of the cross section differs slightly from the one of the data. Like in the NLO NLL case above, we also display the uncertainties for the NNPDF3.1 PDF sets in Figs. 13 and 14. As the PDF errors are found to be small, the theory predictions are rather stable against the uncertainties of current PDFs from global fits and the strong tension between the NLO theory and the data persists. We also note, that recent studies Nocera and Ubiali (2017); Harland-Lang et al. (2017) found it to be impossible to re-constrain the PDFs within a global analysis including inclusive jet data from the LHC when all current cross-correlations among different rapidity bins are consistently taken into account.

Fig. 15 shows again the comparison of the NLO NLL calculations with the CMS data at TeV but now using the NNLO variants of the PDF sets under study. This choice is reasonable to a certain extent, since the NLL resummation includes a dominant part of the full NNLO contributions. On the other hand, potentially large NNLO corrections, for instance possible large corrections from the complete two-loop virtual corrections are still missing in . We can see from Fig. 15 that the NNLO variants of CT14, MMHT2014, NNPDF3.1 get slightly shifted, but are well consistent with the CMS data, again except for the highest values of in the rapidity bin . In contrast, the predictions with the NNLO variant of HERAPDF2.0 are significantly shifted compared to the NLO one, cf. Fig. 11, and display now also good consistency with the CMS data. The NNLO variant of the AMBP16 PDFs predicts the correct shape, but it is lower than the data for all rapidity bins as a consequence of the lower value of compared to used by CT14, HERAPDF2.0, MMHT2014 or NNPDF3.1. This sensitivity to confirms again the great potential of inclusive jet cross section data for the determination of the strong coupling constant Britzger et al. (2017).

Finally, in Fig. 16, we display the results with the PDFs of Bonvini et al. (2015) which have been extracted from data for the DY process, DIS and top-quark hadro-production within the NNPDF framework. These PDFs are subject to improvements at large-, since the theory predictions for DIS and DY as well as for top-quark hadro-production include threshold resummation. For the PDF variant without threshold resummation (labeled as NNPDF30NLO in Fig. 16) both cross sections at NLO NLL and NLO accuracy, and , respectively, are shown. In addition to that, the NLO NLL results for the PDF variant with threshold resummation (labeled as NNPDF30NLL in Fig. 16) are presented as well. Overall, the NLO NLL predictions exhibit better agreement with the data compared to the NLO results, although the PDFs uncertainties of Bonvini et al. (2015) are substantially larger than the ones of global fits. Those large PDF uncertainties at large- and relevant scales of  TeV originate from the gluon PDF at and the light flavor PDFs at lower through the standard parton evolution. The findings in Fig. 16 underpin the necessity to carefully examine and analyze data which constrain those PDFs, including the need to delineate resummation effects from power corrections in the kinematic regions. It will be interesting to observe to what extend improvements can be made in the future in extractions of PDFs with threshold resummation when the inclusive jet data are included.

Figure 16: Same as in Figs. 11 and 13 for the cross sections and using the PDF sets NNPDF3.0 at NLO without and with resummation Bonvini et al. (2015), labeled as NNPDF30NLO and NNPDF30NLL in the plots, respectively. The dashed lines (blue) indicate the PDF uncertainties for with the set NNPDF30NLO.

Iv Summary and Conclusions

In this manuscript, we have provided a detailed study of pQCD calculations from first principles for cross sections of single-inclusive jet production at the LHC. We have performed comprehensive comparisons between the fixed NLO results and the NLL threshold and small- joint resummation improved calculations obtained recently, and we have achieved remarkable advances in perturbative predictions upon using the latter. In our studies, significant differences between the NLO and the NLO NLL joint resummation predictions have been observed in the kinematic regions of interest for the LHC analyses and we have found that these differences account for the discrepancy between the NLO predictions and the LHC data for the jet spectrum in various rapidity bins collected by the CMS experiment at various center-of-mass energies. Once the joint resummation has been included, a remarkable agreement was found between the QCD theory predictions and the LHC data in a large range of jet rapidities.

We have illustrated the impact of the joint resummation in a study of the jet radius ratios at both, and  TeV, which have the advantage of being largely independent of the PDFs and other residual theory uncertainties. At  TeV these jet radius ratios between and , i.e., , have been compared with a CMS analysis in Fig. 2 and overall we have found a significant improvement in the theoretical description of those data. While the predicted double-differential cross sections in and at NLO in pQCD are, for a given value of , systematically higher than the central values of those LHC data in all rapidity bins, arguably they still agree within the theoretical and experimental uncertainties. For the jet radius ratios , however, such consistency is definitely not the case due to the much reduced experimental uncertainties. The NLO predictions for completely miss those LHC data and also cannot be changed by considering different PDF sets since those effects largely cancel out in the jet radius ratios . Therefore, we conclude that the NLO NLL joint resummation is a crucial ingredient in order to achieve a good description of the  TeV jet data within pQCD. We have also presented predictions for the jet radius ratios at  TeV in Fig. 3. using different jet radii with a jet up to . These results will be useful for future experimental analyses of inclusive jet data.

Due to the great importance of the inclusive jet data for constraints on PDFs and determinations of the strong coupling we have also investigated in this study the impact of different PDF sets on the theoretical predictions. We have found that the NLO NLL predictions at  TeV based on the NLO variants of the PDF sets ABMP16, CT14, MMHT2014 and NNPDF3.1 or on the NNLO variant of HERAPDF2.0, respectively, describe the distributions remarkably well for the various rapidity bins. On the other hand, the tension with the CMS inclusive jet data and the pure NLO predictions persists for all those PDF sets and cannot be removed or relieved by selecting a specific PDF set. Therefore, global PDFs which also fit inclusive jet data from the LHC need to be based on theory predictions using the joint resummation for the single-inclusive jet production in order to avoid a possible bias in the PDF extraction due to missing large logarithms in the hard cross sections beyond NLO. We have also noticed in our studies that PDFs extracted with account of threshold resummation but without inclusive jet data have significantly increased PDF uncertainties compared to the ones from the global fits. We suggest to use the joint resummed calculations of the present article in the on-going efforts to refine those PDFs.

Finally, we note that although the NLO NLL calculations greatly improve the theoretical predictions, the associated scale uncertainties are still large and comparable with the current experimental errors. Therefore, in future studies it will be necessary to go beyond the currently achieved accuracy by matching the resummed results with the now available inclusive jet calculations at NNLO and by improving the logarithmic accuracy of the joint threshold and small- joint resummation to NNLL. Both these tasks are feasible within the present framework for resummation and will be subject of future work.

Acknowledgements.
We would like to thank Jun Gao and Katerina Lipka for helpful discussions and Marco Bonvini for providing the PDF set of Bonvini et al. (2015). X.L. would like to thank Zhong-bo Kang and UCLA department of Physics and Astronomy for the hospitality. X.L. is supported by the National Natural Science Foundation of China under Grants No. 11775023 and the Fundamental Research Funds for the Central Universities. S.M. acknowledges contract 05H15GUCC1 by BMBF. F.R. is supported by the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 and by the LDRD Program of Lawrence Berkeley National Laboratory. This research used resources of the Argonne Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC02-06CH11357.

Appendix A Cross sections at  TeV

We present the cross sections for the LHC in Tabs. A.1A.4, respectively, including the theory uncertainty arising from scale variations. At NLO the scale uncertainties have been obtained from the envelope of the variation around up and down by a factor of two, while the scale uncertainties at NLO NLL have been computed as described in section II. The values for the cross sections contain an additional error of not shown explicitly from the numerical integration of the NLO corrections, which dominates both and is correlated between and . The PDF sets used and other parameters are given in the table captions. We also note that for the small bins ( GeV) at  TeV, the threshold resummation may not be applicable anymore.

for LHC at
[GeV] [pb]  (
56 74 1.20 ( 15.06%) 1.17 ( 13.72%) 1.07 ( 14.10%) 8.99 ( 19.41%) 7.25 ( 25.95 %)
74 97 3.49 ( 10.61%) 3.28 ( 11.83%) 2.95 ( 12.10%) 2.47 ( 15.22%) 1.93 ( 21.42 %)
97 133 8.76 ( 9.53%) 8.23 ( 10.41%) 7.31 ( 10.91%) 6.04 ( 12.44%) 4.71 ( 17.60 %)
133 174 2.05 ( 8.50%) 1.93 ( 8.34%) 1.69 ( 8.24%) 1.37 ( 9.75%) 1.02 ( 15.06 %)
174 220 5.66 ( 7.96%) 5.27 ( 8.06%) 4.55 ( 8.41%) 3.63 ( 8.38%) 2.60 ( 13.37 %)
220 272 1.75 ( 7.24%) 1.63 ( 7.04%) 1.39 ( 7.10%) 1.08 ( 7.46%) 7.26 ( 12.22 %)
272 330 5.91 ( 6.86%) 5.42 ( 6.60%) 4.57 ( 6.22%) 3.43 ( 7.08%) 2.14 ( 11.40 %)
330 395 2.11 ( 7.07%) 1.93 ( 6.25%) 1.59 ( 5.80%) 1.15 ( 6.65%) 6.42 ( 10.81 %)
395 468 7.91 ( 6.93%) 7.16 ( 6.10%) 5.79 ( 5.18%) 3.95 ( 6.48%) 1.90 ( 10.41 %)
468 548 3.07 ( 7.22%) 2.75 ( 6.00%) 2.17 ( 4.45%) 1.37 ( 6.35%) 5.44 ( 10.00 %)
548 638 1.23 ( 7.03%) 1.09 ( 5.96%) 8.28 ( 3.99%) 4.78 ( 6.29%) 1.48 ( 9.61 %)
638 737 4.98 ( 7.02%) 4.36 ( 5.73%) 3.17 ( 3.51%) 1.62 ( 6.31%) 3.67 ( 9.15 %)
737 846 2.07 ( 6.86%) 1.77 ( 5.49%) 1.22 ( 3.35%) 5.32 ( 6.33%) 8.05 ( 8.77 %)
846 967 8.63 ( 6.93%) 7.23 ( 5.12%) 4.61 ( 3.45%) 1.66 ( 6.34%) 1.50 ( 8.41 %)
967 1101 3.58 ( 6.73%) 2.92 ( 4.81%) 1.70 ( 3.59%) 4.79 ( 6.34%) 2.19 ( 8.04 %)
1101 1248 1.49 ( 6.69%) 1.17 ( 4.58%) 6.05 ( 3.84%) 1.25 ( 6.39%) 2.21 ( 8.57 %)
1248 1410 6.09 ( 6.53%) 4.57 ( 4.23%) 2.05 ( 4.03%) 2.86 ( 6.42%) 1.18 ( 9.93 %)
1410 1588 2.45 ( 6.35%) 1.73 ( 3.84%) 6.49 ( 4.22%) 5.50 ( 6.48%) 1.75 ( 12.79 %)
1588 1784 9.60 ( 6.23%) 6.30 ( 3.50%) 1.88 ( 4.41%) 8.30 ( 6.55%) 5.45 ( 23.09 %)
1784 2000 3.60 ( 6.03%) 2.16 ( 3.16%) 4.85 ( 4.61%) 8.57 ( 6.64%) —–
2000 2238 1.28 ( 5.89%) 6.84 ( 2.86%) 1.07 ( 4.81%) 4.69 ( 7.07%) —–
2238 2500 4.25 ( 5.70%) 1.96 ( 2.59%) 1.91 ( 5.02%) 7.58 ( 8.68%) —–
2500 2787 1.29 ( 5.49%) 4.91 ( 2.79%) 2.57 ( 5.23%) 4.71 ( 12.81%) —–
2787 3103 3.51 ( 5.22%) 1.04 ( 3.03%) 2.25 ( 5.47%) —– —–
Table A.1: The double-differential cross sections in bins of and at  TeV with using the CT14 Dulat et al. (2015) PDFs at NLO. Theoretical uncertainties from the scale variation are given in parenthesis.
for LHC at
[GeV] [pb]  (
56 74 1.52 ( 2.60%) 1.48 ( 3.96%) 1.36 ( 5.17%) 1.17 ( 4.59%) 9.57 ( 7.17 %)
74 97 4.28 ( 6.03%) 4.05 ( 3.83%) 3.66 ( 5.55%) 3.11 ( 5.34%) 2.47 ( 5.42 %)
97 133 1.05 ( 4.66%) 9.94 ( 4.44%) 8.87 ( 5.48%) 7.38 ( 5.14%) 5.79 ( 5.85 %)
133 174 2.42 ( 5.02%) 2.28 ( 5.83%) 2.00 ( 5.45%) 1.63 ( 5.80%) 1.22 ( 6.05 %)
174 220 6.59 ( 5.65%) 6.15 ( 5.93%) 5.32 ( 4.82%) 4.24 ( 5.31%) 3.05 ( 7.15 %)
220 272 2.02 ( 5.26%) 1.88 ( 5.25%) 1.60 ( 4.75%) 1.24 ( 4.24%) 8.37 ( 4.44 %)
272 330 6.76 ( 4.80%) 6.20 ( 5.21%) 5.21 ( 5.07%) 3.91 ( 4.63%) 2.43 ( 5.70 %)
330 395 2.40 ( 4.98%) 2.18 ( 4.35%) 1.80 ( 4.90%) 1.30 ( 4.83%) 7.23 ( 5.90 %)
395 468 8.92 ( 4.87%) 8.07 ( 4.68%) 6.50 ( 3.62%) 4.42 ( 4.66%) 2.12 ( 6.01 %)
468 548 3.43 ( 4.11%) 3.08 ( 4.18%) 2.42 ( 4.23%) 1.53 ( 4.44%) 6.05 ( 6.08 %)
548 638 1.37 ( 4.23%) 1.21 ( 4.19%) 9.19 ( 4.46%) 5.29 ( 4.77%) 1.63 ( 6.91 %)
638 737 5.53 ( 4.43%) 4.83 ( 4.68%) 3.50 ( 4.83%) 1.78 ( 4.87%) 4.04 ( 8.56 %)
737 846 2.28 ( 5.00%) 1.96 ( 5.02%) 1.34 ( 5.30%) 5.84 ( 5.01%) 8.86 ( 9.05 %)
846 967 9.48 ( 5.16%) 7.94 ( 5.36%) 5.06 ( 5.38%) 1.82 ( 5.06%) 1.65 ( 11.51 %)
967 1101 3.92 ( 5.36%) 3.20 ( 5.65%) 1.86 ( 5.82%) 5.24 ( 5.36%) 2.40 ( 15.04 %)
1101 1248 1.62 ( 5.92%) 1.27 ( 5.94%) 6.60 ( 5.92%) 1.37 ( 5.87%) 2.41 ( 19.70 %)
1248 1410 6.62 ( 6.30%) 4.97 ( 6.35%) 2.24 ( 6.22%) 3.14 (