Abstract

Prediction of production cross section (where ) in proton-proton collisions at  TeV is estimated up to next-to-next-to-leading order (NNLO) in perturbative QCD including next-to-leading order (NLO) electroweak (EW) corrections. The total inclusive boson production cross section times leptonic branching ratio, within the invariant mass window  GeV, is predicted using NNLO HERAPDF2.0 at NNLO QCD and NLO EW as (PDF) () (scale) (parameterization and model). Theoretical prediction of the fiducial cross section is further computed with the latest modern PDF models (CT14, MMHT2014, NNPDF3.0, HERAPDF2.0, and ABM12) at NNLO for QCD and NLO for EW. The central values of the predictions are based on DYNNLO 1.5 program and the uncertainties are extracted using FEWZ 3.1 program. In addition, the cross section is also calculated as functions of and scales. The choice of and for scale variation uncertainty is further discussed in detail.

1. Introduction

The Large Hadron Collider (LHC) has played a crucial role in recent particle physics discoveries and is currently the main tool for exploring TeV scale physics. The LHC is presently the biggest and most powerful proton-proton (pp) collider and allows to test Standard Model (SM) predictions in different center of mass () energies. The LHC started to run with a collision energy of 7 TeV in 2010 and the energy was raised to 8 TeV in 2012. The time period between 2010 and 2012 is called Run I. Then, there was a break between 2012 and 2015 at the LHC for the maintenance and upgrade of the collider. The LHC restarted to run at  TeV (Run II) early in 2015 in order to optimize the delivery of particle collisions for physics research. The second run of the LHC saw an unprecedented data set collected in 2015 and 2016 with more expected in 2017 which opens new possibilities for searches for new physics and precision measurements. The collider was designed to run at a maximum collision energy of 14 TeV; therefore, 14 TeV predictions are considered in this paper. More details about the LHC and its schedule could be found in [1].

Electroweak boson production at LHC provides a major testing ground for quantum chromodynamic (QCD) and electroweak (EW) processes. boson production is one of the most prominent examples of electroweak events. It has large cross section and clean experimental signature, which allow for important precision tests of SM. The process is also very valuable for the calibration of the detector [2, 3] and provides substantial input for the parton distribution function (PDF) determination [46].

Inclusive boson production cross section was previously measured by the ATLAS and CMS Collaborations at the LHC at  TeV [7, 8],  TeV [9], and  TeV [1012]. ATLAS and CMS publications showed that there is a good agreement between the observed and NNLO predicted results. Therefore, the theoretical predictions of higher-order corrections are essential work to validate the experimental results.

Theoretical prediction of boson is available up to next-to-next-to-leading order (NNLO) in perturbative QCD and to next-to-leading order (NLO) for EW processes. The cross section definition of a weak boson in pp collisions is given by (1), and the process is illustrated in Figure 1. Here , are the incoming partons with the momentum fraction , represent PDFs, and stands for partonic cross section. PDFs are known with limited precision, and partonic cross section is calculable to high orders (QCD, EW). Total and fiducial cross section predictions of events in pp collisions are presented providing PDF, strong coupling (), scale, model, and parameterization uncertainties based on DYNNLO 1.5 [13] and FEWZ 3.1 [14]. The central values of the predictions are calculated using DYNNLO 1.5, and the theoretical uncertainties and the impact of NLO EW correction on the cross section are extracted using FEWZ 3.1 MC generator program. The disagreement between FEWZ 3.1 and DYNNLO 1.5 is previously reported as smaller than 1 [15]. The selections used in this analysis are the same as those used in the latest ATLAS study at 13 TeV [10]. The total and fiducial cross section are estimated up to NNLO QCD including NLO EW corrections with the most recent PDF sets (HERAPDF2.0 [16], CT14 [17], NNPDF 3.0 [18], MMHT 2014 [19], and ABM12 [20]).

Another topic investigated in this paper is the parameter choices for scale uncertainty. The scale uncertainty depends on the choice of renormalization () and factorization () scales. Therefore, two different approaches for scale uncertainty calculation are also studied using HERAPDF2.0, as well as the cross section dependence of and scales used in the perturbative calculation. The standard convention of considering the range is followed, setting . Besides the standard convention, the possibility of varying independently the value of renormalization and factorization scale is investigated, and the main findings are discussed.

2. Total and Fiducial Cross Section

The total inclusive boson production cross section times leptonic branching ratio, within the invariant mass window  GeV, is estimated at  TeV in pp collisions. A fiducial region is further defined and boson production cross section is computed. For the fiducial cross section, the leptons are required to have and  GeV. The invariant mass window is kept as  GeV. All presented results are for born level (pre-QED FSR) leptons. HERAPDF2.0, CT14, NNPDF3.0, MMHT2014, and ABM12 PDF models are used with the corresponding value of and  GeV. For completeness, we first compute the predictions for 13 TeV and the predicted values are compared with ATLAS collaboration result [10]. Then, we rerun the code for 14 TeV center of mass energy.

Figure 2 presents the 13 TeV and 14 TeV predictions we performed. The upper two figures prove that there is a good agreement between experimental measurements and theoretical predictions. Particularly, HERAPDF2.0 prediction at NNLO best describes the measurement. The solid thick lines on upper figures represent the central value of ATLAS measurement. The dashed and dotted lines show the quadrature sum of the luminosity and systematic uncertainties and only systematic uncertainty, respectively. Black points stand for the central values of the predictions. The green bands present only scale uncertainty and the yellow bands illustrate the quadrature sum of PDF, scale, and for MMHT2014, NNPDF3.0, CT14, and ABM12 PDF predictions if they are available. The green band for HERAPDF2.0 is the quadrature sum of PDF, scale, , parameterization, and model uncertainties.

The PDF uncertainty is estimated following closely the prescription of PDF4LHC working group [21, 22]. These uncertainties are evaluated with five different NNLO PDFs for the fiducial cross section: HERAPDF2.0, CT14, NNPDF3.0, MMHT2014, and ABM12. The PDF uncertainty of the CT14 PDF set is rescaled from 90% CL to 68% CL to match the other sets. The PDF uncertainty of the total cross section is estimated using LO, NLO, and NNLO HERAPDF2.0 PDF sets at LO, NLO, and NNLO QCD, respectively. The uncertainty is estimated varying by 0.001. Then, the uncertainties are symmetrized by taking the bigger value from estimated up and down errors. The scale uncertainties are calculated varying the value of and scales. The default value of is set to , which is taken as 91.1876 GeV. The scales are changed by factors of two, , where . The maximum value of the variation is taken as the scale uncertainty, which makes the scale uncertainty symmetric. To calculate the parameterization and model uncertainties, HERAPDF20_NLO_VAR and HERAPDF20_NNLO_VAR PDF sets are used. Each PDF set has 13 eigenvectors and each variation is treated one by one to estimate the model and parameterization errors. First, model and parameterization errors are calculated separately, and, then, they are summed in quadrature. More details about the calculation of model and parameterization errors can be found in [23]. After the calculation of all considered uncertainties, the total uncertainty is determined by adding the individual uncertainties in quadrature. The numerical values are provided in Results.

3. Renormalization and Factorization Scales

The source of the renormalization scale is ultraviolet divergence, which is one of the most important aspects of field theories, and the factorization scale arises from collinear (infrared) divergence. Proper treatment of these divergences is crucial to understand the physics processes and perturbation theory. Any wrong choice of and scales may lead to an unstable perturbative behavior. The details and definition of the divergences and the scales can be found in a well written paper [24].

The hadronic cross section production for Drell-Yan production or other LHC processes is already defined in (1). The partonic cross section () can be written in terms of LO, NLO, and NNLO as follows:The cross section statement turns into (3), which is the cross section definition in terms of LO, NLO, NNLO, and so on.We may also express the cross section in terms of analytical expression of these scales.where is the experimental energy scale and it is set to in our case. The hard scattering function is presented as an analytical expression of and and it can be further expanded as follows:where , , and terms give the LO, NLO, and NNLO contributions, respectively. term represents the contributions beyond NNLO calculations. FEWZ and DYNNLO programs let us set the factorization and renormalization scales different to calculate the cross sections, so they successfully serve the purpose of our paper.

Figure 3 shows the total cross section dependence on and scales. The values are estimated at LO, NLO, and NNLO QCD using DYNNLO 1.5 MC generator program interfaced with LO, NLO, and NNLO HERAPDF2.0 PDF sets. is a renormalization scale dependent parameter and does not appear at leading order since LO boson diagram is purely electroweak. Based on (3), it starts to take a role at NLO. This feature can be proven by plotting total cross section as a function of . Figure 3(b) presents LO, NLO, and NNLO total cross section predictions as a function of . Here is fixed to while is varied. It is clearly seen that LO QCD prediction is constant for different values. NLO and NNLO QCD predictions decrease by increase of . NNLO prediction is always higher than NLO prediction. In Figure 3(a),    and scales are varied by setting the condition: . The cross section increases by increase of and . The increase is the least for NNLO QCD prediction and the most for LO. After some points, the LO and NLO QCD predictions become higher than NNLO QCD prediction. As can be seen from (3), PDFs are dependent on parameter, and scale starts to contribute the calculations at LO. Figure 3(c)   illustrates the cross section variation by . Here is set to while varies. The findings are similar to the one observed for the variation of scales with condition of .

There are two approaches to calculate the scale uncertainty. The first one is the standard convention researchers often used. It is considered that is multiplied by half or 2, where . Then, the difference with the default value is taken as the scale uncertainty (see (6)). The second way is varying independently the value of renormalization and factorization scales. The maximum variation is assigned as scale uncertainty (see (7)).where and could be 1/2, 1, or 2. We have total 9 variations including the default one (, ). Those are listed as (1, 1), (1, 2), (1, 1/2), (2, 1), (2, 2), (2, 1/2), (1/2, 1), (1/2, 2), and (1/2, 1/2). The scale uncertainty calculated at NLO and NNLO QCD is illustrated by Figure 4. The absolute value of each variation is also given in Table 1. The renormalization and factorization scales have different origins as stated above. Therefore, they should be varied independently by a factor of two. It is however found that large logarithms arise when the two scales differ by a factor of 4, and this gives spuriously large uncertainties, so these special configurations, that is, , should be removed from the calculation of the scale uncertainties. For example, at NNLO QCD, the scale uncertainty for two approaches is as follows: and . The uncertainty with is about two times bigger than the standard convention method. At NLO QCD, the scale uncertainty for two approaches is as follows: and . This time, is almost three times bigger than .

4. Results

The total and fiducial cross sections of inclusive boson production are given in Tables 2 and 3. The numbers are obtained with LO, NLO, and NNLO QCD and NLO EW corrections as stated before. Here the total uncertainty is the quadrature sum of the PDF, , scale, model, and parameterization errors. All are calculable for NNLO and NLO HERAPDF2.0 PDF sets but model and parameterization errors are not available for LO HERAPDF2.0 set. PDF, and scale uncertainties are calculable for NNLO MMHT2014, NNPDF3.0, and CT14 PDF sets. Only PDF and scale uncertainties are available for NNLO ABM12 PDF set.

The cross section values of full phase-space are presented in Table 2. Based on these numerical values, it can be safely stated that the scale uncertainty is getting smaller from LO to NNLO QCD as expected. The scale uncertainty at NNLO is 25.75 of the calculated total uncertainty and the scale uncertainty contributions to the total uncertainty at NLO and LO are 84.61 and 99.58, respectively. Table 3 provides the fiducial cross section numbers of inclusive production of boson in pp collisions at NNLO QCD and NLO EW. These numbers are produced using five different PDF models: HERAPDF2.0, MMHT2014, NNPDF3.0, CT14, and ABM12. The table proves that the predictions based on different PDF models are consistent with each other. The CT14 prediction has the biggest total uncertainty. The MMHT2014 prediction has the biggest scale uncertainty but overall all PDF models have the same level of scale uncertainties. ABM12 has the smallest total uncertainty since only PDF and scale uncertainties are calculable for ABM12 predictions. As can be seen from these numbers, the scale parameters are crucial for the QCD calculations.

The numerical values of cross sections based on variation of and are given in Table 4. The table also provides the cross section ratios at different orders, which is so-called -Factor. Table 4 is already illustrated in Figure 3. In Figure 3(a), LO prediction has the sharpest slope while NNLO prediction has the smallest slope. This illustrates why the scale uncertainty is the biggest at LO and the smallest at NNLO. If we take two points, for example, 0.5 and 2, on the prediction line and take the difference between these points and default value, which is 1, we will have the biggest difference at LO due to having this biggest slope. Another remark should be underlined here that the predicted cross section value as a function of is not changing as much as it changes by . Even LO QCD prediction does not show any reaction to the change of as discussed before.

The ratio of cross sections provided in Table 4 is highly dependent on factorization and renormalization scales. These numbers are extracted using HERAPDF2.0. Under normal conditions, an increase of cross section by increase of QCD orders is expected for most of the physics processes. That means the provided ratios in the table need to be above 1. However, / ratio drops below 1 in and , conditions. NLO prediction in the corresponding values will provide an overestimate of the NNLO cross section. It is also observed that LO predictions are getting greater than NNLO predictions at high points. There is a trend of decrease of / ratio in three conditions provided in Table 4. This is also true for / ratio in conditions of and , . The same behavior is not observed on / while the condition is set to , . All these outcomes draw a conclusion that wrong choices of renormalization and factorization scales may lead to bad perturbative behaviors.

5. Conclusion

The main findings of these studies are reiterated here. First, total and fiducial cross section predictions at 13 TeV are compared with ATLAS study at 13 TeV to validate our code. It is observed that there is a good agreement between the predicted and observed results at 13 TeV and the comparisons showed that HERAPDF2.0 best describes the experimental results. The 14 TeV predictions are further computed by running the code used at 13 TeV comparison. 14 TeV fiducial cross section predictions are presented at NNLO QCD including NLO EW corrections for five different PDF models: HERAPDF2.0, MMHT2014, NNPDF3.0, CT14, and ABM12. All PDF predictions are in well agreement with each other. Total cross section predictions at 14 TeV are also reported using HERAPDF2.0 at LO, NLO, and NNLO QCD including NLO EW corrections. In the calculations, PDF, , scale, model, and parameterization errors are considered. It is observed that the scale uncertainty is a dominant uncertainty at LO and NLO calculations and its contribution to the total uncertainty decreases from LO to NNLO.

Another topic studied here is the choice of and parameters. The obtained results show that the cross section increases by increase of and decreases by increase of . Varying the two scales simultaneously leads to a compensation of the two different behaviors. As a result, the scale dependence of cross section is mostly driven by the factorization scale. These scale dependence characteristics are reduced when high-order corrections are included. Another important issue observed here is that NLO calculations at certain values provide an overestimate of the NNLO cross sections. The results also present a similar overestimation for LO predictions for high values. For certain renormalization and factorization values, LO predictions are higher than NNLO predictions although an increase of cross section from LO to NNLO is expected. and choices on the determination of scale uncertainty are further investigated. and scales are altered independently, where and . Here could be 1, 1/2, and 2. We set 9 variations as listed (, ): (1, 1), (1, 2), (1, 1/2), (2, 1), (2, 2), (2, 1/2), (1/2, 1), (1/2, 2), and (1/2, 1/2). We found that there is very large uncertainty when . Based on the numbers presented in Table 1, it is found that the scale uncertainty is not altered significantly if the uncertainty at is excluded. All these findings draw a conclusion that there is no right choice but smart choice; however, any wrong choice may lead to an unexpected perturbative behavior.

Conflicts of Interest

The authors declare that they have no conflicts of interest.