Azimuthal correlations of high transverse momentum jets at next-to-leading order in the parton branching method

The azimuthal correlation, Δϕ12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \phi _{12}$$\end{document}, of high transverse momentum jets in pp collisions at s=13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=13$$\end{document} TeV is studied by applying PB-TMD distributions to NLO calculations via MCatNLO together with the PB-TMD parton shower. A very good description of the cross section as a function of Δϕ12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \phi _{12}$$\end{document} is observed. In the back-to-back region of Δϕ12→π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta \phi _{12}}\rightarrow \pi $$\end{document}, a very good agreement is observed with the PB-TMD Set 2 distributions while significant deviations are obtained with the PB-TMD Set 1 distributions. Set 1 uses the evolution scale while Set 2 uses transverse momentum as an argument in αs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _\mathrm {s}$$\end{document}, and the above observation therefore confirms the importance of an appropriate soft-gluon coupling in angular ordered parton evolution. The total uncertainties of the predictions are dominated by the scale uncertainties of the matrix element, while the uncertainties coming from the PB-TMDs and the corresponding PB-TMD shower are very small. The Δϕ12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \phi _{12}$$\end{document} measurements are also compared with predictions using MCatNLO together Pythia8, illustrating the importance of details of the parton shower evolution.


Introduction
The description of the cross section of high p T jets in protonproton (pp) collisions is one of the most important tests of predictions obtained in Quantum Chromodynamics (QCD), and much progress has been achieved in the description of inclusive jets [1-12] by applying next-to-leading (NLO) [13-16] and next-to-next-to-leading-order (NNLO) calculations [17][18][19][20]. In multijet production, the azimuthal angle φ 12 between the two highest transverse momentum p T -jets is an inclusive measurement of additional jet radiation. At leading order (LO) in strong coupling α s , where only two jets are present, the jets are produced back-to-back, with φ 12 = π , while a deviation from this back-to-back configuration indi-cates the presence of additional jets, and only higher-order calculations can describe the observations. The azimuthal correlation between two jets has been measured in pp collisions at a center-of-mass energy of √ s = 1.96 TeV by the D0 collaboration [21,22] and in pp collisions by the ATLAS Collaboration at √ s = 7 TeV [23] and by the CMS Collaboration at √ s = 7, 8, and 13 TeV [24][25][26][27]. When measurements of azimuthal correlations of dijets are compared with LO or NLO computations supplemented by parton showers, deviations of 50% are observed in the medium φ 12 region even at NLO (see e.g. [25,26]), which requires a more detailed understanding. In the φ 12 → π region, deviations of up to 10 % are observed [27], significantly larger than the experimental uncertainties.
Since initial state parton radiation moves the jets away from the φ 12 = π region, it is appropriate to investigate the implications of transverse momentum dependent parton densities (TMDs [28]) in the description of the φ 12 measurements. Kinematic effects of the initial-state transverse momenta in the interpretation of jet measurements were pointed out in [29,30]. The region φ 12 → π is especially sensitive to soft multi-gluon emissions, for which QCD resummation is needed. Calculations at leadinglogarithm have been obtained in Ref. [31]. A calculation based on TMD distributions is found in Refs. [32,33] and further investigated in [34,35]. However, in the region φ 12 → π soft-gluon effects are expected which lead to socalled factorization-breaking [36][37][38]. An indirect strategy to explore the potential impact of these effects is to compare calculations which assume factorization with high-precision measurements.
The Parton Branching (PB)-method [39,40] allows one to determine TMD parton distributions. With these PB-TMD distributions a very good description of the Drell-Yan process [41] is achieved at the LHC [42] as well as at lower energies [43]. Drell-Yan lepton pair production in association with jet final states is also well described using the TMD jet merging [44]. In Ref. [45] it is shown that Z + b production is also well described. TMD parton distributions have been used together with off-shell matrix elements at the lowest order in Refs. [46][47][48] showing a reasonably good description of the measurements.
In this article we investigate in detail highp T dijet production by applying the PB formulation of TMD evolution together with NLO calculations of the hard scattering process in the MadGraph5_aMC@NLO [49] framework. We first give a very brief recap of the PB distributions in Sect. 2. In Sect. 3 we describe how TMDs and TMD parton showers are included in the Monte Carlo generator Cascade3 [50]. We discuss predictions obtained by applying fixed-order NLO perturbative calculations and study the region where soft gluon resummation becomes important. We show predictions using PB-TMDs together with TMD parton shower in Sect. 4.
We compare these predictions with the one using the Pythia8 parton shower. We finally give conclusions in Sect. 5.

PB TMDs
The PB method [39,40] provides a solution of evolution equations for collinear and TMD parton distributions. The equations are solved by applying the concept of resolvable and non-resolvable branchings with Sudakov form factors providing the probability to evolve from one scale to another without resolvable branching. The method is described in Refs. [40,51].
For the numerical calculations we use the NLO parton distribution sets, PB-NLO-2018-Set 1 and PB-NLO-2018-Set 2, as obtained in Ref. [51] from a fit to inclusive deep inelastic scattering precision measurements at HERA [52]. Both the collinear and TMD distributions are available in TMDlib [53,54], including uncertainty bands. PB-NLO-2018-Set 1 corresponds at collinear level to HERA-PDF 2.0 NLO [52], while PB-NLO-2018-Set 2 uses transverse momentum (instead of the evolution scale in Set1) for the scale in the running coupling α s which corresponds to the angular ordering of soft gluon emissions in the initial-state parton evolution [55][56][57][58].
In Fig. 1 the distributions of the collinear densities from Set 1 and Set 2 are shown for up-quarks at evolution scales of μ = 100 and 1000 GeV , typical for multi-jet production described below. The collinear densities are also available in a format compatible with LHAPDF [59], and can be used in calculations of physical processes at NLO. In Fig. 2 we show the TMD distributions for up-quarks at x = 0.01 and μ = 10 and 100 GeV . The differences between Set 1 and Set 2 are clearly visible in the small k T -region.
The uncertainties of the TMD distributions include both experimental and model uncertainties, as determined in Ref. [51]. In general, it is observed that those uncertainties are small; for k T > 1 GeV they are of the order of 2-3 %.

Multijet production
The predictions for multijet production at NLO are obtained using the MadGraph5_aMC@NLO [49] framework. We used MadGraph5_aMC@NLO in two different modes: one is the fixed NLO mode, in which only partonic events are produced, without parton shower and hadronization, and the other one is the real MC@NLO mode, in which infrared subtraction terms are included to avoid double counting of parton emissions between matrix-element and parton-shower calculations, so that events need to be supplemented with a parton shower (or with PB TMD evolution) in order to produce a physical cross section. Fixed NLO dijet production is calculated within the Mad-Graph5_aMC@NLO framework. Technically, in the fixed NLO mode, MadGraph5_aMC@NLO (version 2.9.3) produces event files with the partonic configuration in LHE format [60] which can be processed through Cascade3 [50] combining events and counter events (due to infrared subtraction) so that they are treated as one event for a proper calculation of statistical uncertainties. In the fixed NLO mode, the MadGraph5_aMC@NLO event record is kept without any modification. Processing through Cascade3 has a significant advantage that a fixed NLO calculation can be obtained making use of all the analyses coded in Rivet [61].
In the MC@NLO mode, subtraction terms are included which depend on the parton shower used. For the PB-TMDs and the PB-TMD parton shower we use Herwig6 [62,63] subtraction terms, as already applied in Z and Drell-Yan analyses [42,43], motivated by the angular ordering in the PB evolution. MadGraph5_aMC@NLO (version 2.6.4, hereafter labeled MCatNLO) [49] together with the NLO PB parton distributions with α s (M Z ) = 0.118 is used for NLO calculation of dijet production. The matching scale μ m , which limits the contribution from PB-TMDs and TMD showers (μ m =SCALUP included in the LHE file), guarantees that the overlap with real emissions from the matrix element is minimized according to the subtraction of counterterms in the MC@NLO method. The factorization and renormalization scale in MCatNLO is set to μ R,F = 1 2 i m 2 i + p 2 t,i , where the index i runs over all particles in the matrix element final state. This scale is also used in the PB-TMD parton distribution A(x, k T , μ).
In Cascade3, as described in detail in Ref. [50], the transverse momentum of the initial state partons is calculated according to the distribution of k T provided by the PB-TMD A(x, k T , μ) at given longitudinal-momentum fraction x and evolution scale μ. This transverse momentum is used for the initial state partons provided by MCatNLO, and their longitudinal momentum is adjusted such that the mass and the rapidity of the dijet system is conserved, similar to what has been done in the Drell-Yan case [43]. The initial state TMD parton shower is included in a backward evolution scheme, respecting all parameters and constraints from the PB-TMD. The kinematics of the hard process are not changed by the shower, after the k T from the TMD is included. The final state parton shower is obtained with the corresponding method implemented in Pythia6 [64], by vetoing emissions which do not satisfy angular ordering (MSTJ(42)=2).
In Fig. 3 we show results for the transverse momentum distribution of the dijet system p T,12 and the azimuthal correlation φ 12 between the two leading jets as obtained from the MCatNLO calculation at fixed NLO (blue curve), at the level including subtraction terms (LHE level, green curve) and after inclusion of PB-TMDs (red curve). One can clearly observe the rising cross section of the fixed NLO calculation towards small p T,12 (or at large φ 12 ). This is the region in p T,12 and φ 12 where the subtraction terms are relevant and a physical prediction is obtained when PB-TMDs and parton showers are included. The jets are defined with the  anti-k T jet-algorithm [65], as implemented in the FASTJET package [66], with a distance parameter of R=0.4 and a transverse momentum p T > 200 GeV. The use of jets (instead of partons) is the reason for the tail towards small φ 12 in the MCatNLO(LHE) and MCatNLO(fNLO) calculation.

Azimuthal correlations in multijet production
We next apply the framework described in the previous section, based on the matching of PB-TMDs with NLO matrix elements, to describe the measurement of azimuthal correlations φ 12 obtained by CMS at √ s = 13 TeV [26] and in the back-to-back region ( φ 12 → π ) [27]. Only leading jets with a transverse momentum of p > 1000 GeV, where the jets appear very collimated. We apply the collinear and TMD set PB-NLO-2018-Set 2, unless explicitly specified, with running coupling α s (m Z ) = 0.118. We may estimate the theoretical uncertainties on the predictions by considering two kinds of uncertainties: those that come from variation of the arbitrary scales that appear in the various factors that enter the jet cross section, and those that come from the determination of the TMD parton distributions and showers. The former include Shown are the uncertainties coming from the scale variation (as described in the text) as well as the uncertainties coming from the TMD the renormalization scale in the strong coupling, the factorization scale used in the parton distribution and the matching scale to combine the matrix element and PB TMD. The latter include both experimental and model uncertainties in the TMD extraction. As regards the scale variations, we present results corresponding to the 7-point scheme variation around the central values for the renormalization and factorization scale (avoiding the extreme cases of variation). We have studied the dependence on the matching scale μ m and found that is within the band of variation of factorization and renormalization scales. The experimental and model uncertainties on the determination of the TMD distributions as described in [51] are included. Shown are the uncertainties coming from the scale variation (as described in the text) as well as the uncertainties coming from the TMD In Fig. 4 we show a comparison of the measurement by CMS [26] for different values of p leading T with the calculation MCatNLO+CAS3 including PB-TMDs, parton shower, and hadronization. The uncertainties from scale variation and TMD determination are shown separately.
In Fig. 5 the measured φ 12 distribution [27] in the backto-back region is compared with the prediction MCatNLO +CAS3.
In general, the measurements are very well described, especially in the back-to-back region. The scale uncertainty is significantly larger than the TMD uncertainty, especially in the low p leading T region. A difference between the measure- In Figs. 6 and 7, the predictions using PB-NLO-2018-Set 1 are compared with those from PB-NLO-2018-Set 2 and with the measurements. The difference between Set 1 and Set 2 becomes significant in the back-to-back region, which is sensitive to the low k T -region of the TMD. As already observed in the case of Z-boson production in Ref. [42], Set 2 with the transverse momentum as a scale for α s , which is required from angular ordering conditions, allows a much better description of the measurement. It has been explicitly checked that the choice of the collinear parton density function (in contrast to the choice of the TMD densities) does not matter for the φ 12 distributions, since they are normalized. The region of low φ 12 in Figs. 4 and 6 is not well described with an NLO dijet matrix element calculation supplemented with TMD densities and TMD parton shower because in the low φ 12 region higher-order hard emissions play a significant role. It has been shown in [67] that the inclusion of higher order matrix elements with the new TMD merging method of Ref. [44] leads to a very good description of the low φ 12 region.
In Fig. 8 predictions obtained with MCatNLO+Pythia8 are compared with MCatNLO+CAS3. In the calculation of MCatNLO+Pythia8, the Pythia8 subtraction terms are used and the NNPDF3.0 [68] parton density and tune CUETP8M1 [69] are applied. The uncertainties of the PYTHIA prediction are derived by combining the fixed-order scale variation from MCatNLO with renormalization scale variations in the parton shower. We use the method of [70] together with the guidelines of [71] to obtain consistent scale variations where possible. In particular, this means that the renormalization scale variation at fixed order and in the parton shower are fully correlated. 1 The factorization scale variation is only applied at fixed order, as argued in [71]. We observe a significant dependence on the matching scale μ m , the details of matching in case of dijets needs further investigation.
Shown in Fig. 8 is also the contribution from multiparton interactions, which is very small for jets with p leading T > 200 GeV. The prediction obtained with MCatNLO+Pythia8 is in all φ 12 regions different from the measurement and MCatNLO+CAS3, illustrating the role of the treatment of parton showers.
In conclusion, the predictions of MCatNLO+CAS3 are in reasonable agreement with the measurements in the larger φ 12 regions, where the contribution from higher order matrix elements is small. In the back-to-back region ( φ 12 → π ) the predictions obtained with PB-TMDs and parton shower are in good agreement with the measurement. The uncertainties of the predictions are dominated by the scale uncertainties of the matrix element calculations, while the PB-TMD and TMD shower uncertainties are very small, as they are directly coming from the uncertainties of the PB-TMDs. No uncertainties, in addition to those from the PB-TMD, come from the PB-TMD parton shower.

Conclusion
We have investigated the azimuthal correlation of high transverse momentum jets in pp collisions at √ s = 13 TeV by applying PB-TMD distributions to NLO calculations via MCatNLO. We use the same PB-TMDs and MCatNLO calculations as we have used for Z -production at LHC energies in Ref. [42]. A very good description of the cross section as a function of φ 12 is observed. In the back-to-back region of φ 12 → π a very good agreement is observed with PB-TMD Set 2 distributions [51] while significant differences are obtained with PB-TMD Set 1 distributions, which use the evolution scale as an argument in α s . This observation confirms the importance of consistently handling the softgluon coupling in angular ordered parton evolution.
The uncertainties of the predictions are dominated by the scale uncertainties of the matrix element, while uncertainties coming from the PB-TMDs and the corresponding PB-TMD shower are very small. No other uncertainties, in addition to those of the PB-TMD, come from the PB-TMD shower, since it is directly correlated with the PB-TMD density.
We have also investigated predictions using MCatNLO with Pythia8 to illustrate the importance of details of the parton shower.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permit-ted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .