Multi-Regge kinematics and azimuthal angle observables for inclusive four-jet production

We evaluate differential cross sections for production of four jets in multi-Regge kinematics at a hadron collider. The main focus lies on azimuthal angle dependences. As in previous studies, the ratios of correlation functions of products of cosines of azimuthal angle differences among the tagged jets offer us the cleanest quantities to compare with experimental data. The calculations are based on the jet production from a single BFKL ladder with a convolution of three BFKL Green functions where we always have two forward/backward jets tagged in the final state. We also demand the tagging of two further jets in more central regions of the detectors with a relative separation in rapidity from each other, plus the inclusive production of an arbitrary number of mini-jets. We show that dependences on the transverse momenta and rapidity of the two central jets can be a distinct signal of the onset of BFKL dynamics.


Introduction
The study of the asymptotic behavior of scattering amplitudes in the limit of high center-of-mass energy is an active area of research for particle phenomenology. The Large Hadron Collider (LHC) is producing an abundance of data allowing for the study of very exclusive observables with stringent cuts in the final state. One of the key points for understanding multi-jet production at high energies is the multi-Regge kinematics (MRK). MRK is the kinematics that, by presupposing a strong ordering in rapidity for the final state jets, allows for large logarithms in the colliding energy to be present in all orders of the perturbative expansion. This fact alone calls for a resummation framework of the large logarithms in energy.
In the high energy (Regge) limit, the common basis for the perturbative description of a hard process in QCD is the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach, at leading (LL) [1,2,3,4,5,6] and next-to-leading (NLL) [7,8] accuracy. This approach offers a resummation of those enhanced terms in MRK in regions of phase space where a fixed order calculation might not be enough. This formalism has been successfully applied to lepton-hadron Deep Inelastic Scattering at HERA (see, e.g. [9,10]) to describe quite inclusive processes which are not that suitable though if one is interested in discriminating between BFKL dynamics and other resummation programs. At the LHC, however, it is possible to investigate processes with much more exclusive final states which could, in principle, be only described by the BFKL framework. This would allow us to precisely determine the applicability window of the framework.
With this idea in mind there has been a lot of recent activity in the study of the so-called Mueller-Navelet jets [11], i.e. the inclusive hadro-production of two forward jets with large and similar transverse momenta and a big relative separation in rapidity Y , proportional to √ s, and with associated inclusive mini-jet radiation. Interesting observables associated to this process are the azimuthal angle (θ) correlations cos (n θ) of the two tagged jets, and it has been shown [12,13,14,15] that the further gluon radiation manifests as a fast decrease of these functions with Y . However, these observables are strongly affected by collinear effects [16,17] due to their dependence on the n = 0 Fourier component in θ of the BFKL kernel, which is strongly dependent on collinear radiation. In order to remove this problem, new observables were proposed [16,17] which are independent from the n = 0 contribution: the ratios R M N = cos (M θ) / cos (N θ) . They have been calculated at NLL [18,19,20,21,22] and show a very good agreement with experimental data at the LHC.
Nevertheless, Mueller-Navelet configurations are still too inclusive to study MRK with precision. A step toward other observables capable to pin down the MRK dynamics in much more detail has been taken in [23], where a new study was proposed that demands the tagging of a third, central in rapidity, jet within the usual Mueller-Navelet configuration. It is important to remain within the general Mueller-Navelet setup because having two forward/backward jets allows the use of collinear factorization which is in a better theoretical control than k t -factorization. Since a unique footprint of BFKL physics is its azimuthal angle dependence, the main new observables studied in [23] are the ratios with φ 1 and φ 2 being the azimuthal angle difference respectively between the first and the second (central) jet and between this one and the third jet. These observables depend strongly on the p t and less strongly on the rapidity of the central jet, and this information can be used to probe characteristic properties of the BFKL ladder in a very precise way. The present work, original results of which are presented in the next Section, is a natural continuation of [23], by allowing the production of a second central jet, thus making it possible to define more differential distributions in the transverse momenta, azimuthal angles and rapidities of the two central jets, for fixed values of the four momenta of the forward jets.

Inclusive four-jet production
We now present the analysis of events with two forward/backward jets together with two more central jets, all of them well separated in rapidity from each other, making use of the BFKL formalism to describe the associated inclusive multi-jet emission. The two tagged forward/backward jets A and B have transverse momentum k A,B , azimuthal angle ϑ A,B and rapidity Y A,B , while the pair of tagged more central jets are characterized, respectively, by k 1,2 , ϑ 1,2 and y 1,2 . The differential cross section on these latter variables can be written in the form Here we have introduced the rapidity ordering characteristic of MRK: Y A > y 1 > y 2 > Y B ; and k 2 1 , k 2 2 lie above the experimental resolution scale. ϕ are BFKL gluon Green functions normalized to ϕ ( p, q, 0) = δ (2) Following the course taken in Ref. [23], our goal is to define and study the behavior of observables for which the BFKL approach will show distinct features with respect to other formalisms and, if possible, are also quite insensitive to higher-order corrections. We start with the study of a quantity similar to the usual Mueller-Navelet case such that we integrate over the azimuthal angles of the two central jets and over the difference in azimuthal angle between the two forward jets, wherẽ (ψ is the logarithmic derivative of Euler's gamma function). The associated experimental observable corresponds to the mean value of the cosine of ∆θ = ϑ A − ϑ B − π in the recorded events: In order to improve the perturbative stability of our predictions (see [24] for a related discussion) it is convenient to remove the contribution from the zero conformal spin (which corresponds to the index n = 0 in Eq. (4)) by defining the ratios where we consider M, N as positive integers. Our next step now is to propose new observables, different from those characteristic of the Mueller-Navelet case though still related to azimuthal angle projections. We thus define with Ω m,n,l = +∞ 0 dp In order to drastically reduce the dependence on collinear configurations we can remove the zero conformal spin contribution by defining the following ratios: with integer M, N, L, P, Q, R > 0.
It is now possible to numerically investigate many different momenta configurations. In order to cover two characteristic cases, namely k A ∼ k B and k A < k B (or equivalently k A > k B ) we choose the following two fixed configurations for the transverse momenta of the forward jets: (k A , k B ) = (40, 50) and (k A , k B ) = (30, 60) GeV. We also fix the rapidities of the four tagged jets to the values Y A = 9, y 1 = 6, Y 2 = 3, and Y B = 0 whereas the two inner jets can have transverse momenta in the range 20 < k 1,2 < 80 GeV.
In Fig. 1 we present our results for the normalized coefficients C 111 , C 112 , C 121 and C 122 after they are divided by their respective maximum. We find that the distributions are quite similar for the two configurations here chosen ((k A , k B ) = (40, 50), (30, 60) GeV) apart from the coefficient C 121 which is quite more negative for the latter configuration when the transverse momentum of the first central jet, k 1 , is low. Further coefficients, normalized as above, are calculated in Fig. 2 for the cases C 211 , C 212 , C 221 and C 222 . Again they are rather similar with the exception of C 221 at low p t of one of the centrals jets with largest rapidity.
Since these coefficients change sign on the parameter space here studied, it is clear that for the associated ratios R M N L P QR there will be some lines of singularities. We have investigated R 121 212 , R 212 211 and R 221 222 in Fig. 3. In this case the configurations (k A , k B ) = (40, 50), (30, 60) GeV behave quite differently. This is due to the variation of the position of the zeroes of those coefficients C M N P chosen as denominators in these quantities. It would be very interesting to test if these singularity lines are present in any form in the LHC experimental data. A further set of ratios, R 111 112 , R 111 122 , R 112 122 and R 222 211 , with their characteristic singular lines, is presented in Fig. 4. In general, we have found a very weak dependence on variations of the rapidity of the more central jets y 1,2 for all the observables here presented.
We have used both Fortran and Mathematica for the numerical computation of the ratios R M N L P QR . We made extensive use of the integration routine Vegas [25] as implemented in the Cuba library [26,27]. Furthermore, we used the Quadpack library [28] and a slightly modified version of the Psi [29] routine.
We conclude here our numerical analysis. We have offered several interesting observables where BFKL effects could have sizable effects. More detailed calculations are needed, including the introduction of parton distribution functions effects and higher order terms in jet vertices and Green functions. However, we argue that the bulk of the relevant contributions are already contained in the calculations here discussed, in particular, for the ratios R M N L P QR . It will    be very important to compare against the BFKL Monte Carlo code BFKLex [30,31,32,33,34,35,36] as well as to calculate the same quantities with other, more conventional, approaches [37,38,39] in order to gauge if they differ from our results. This includes those analysis where the four-jet predictions stem from two independent gluon ladders [40,41].

Summary & Outlook
We have presented new observables to study four-jet production at hadron colliders in terms of its azimuthal angle dependences. These correspond to the ratios of correlation functions of products of cosines of azimuthal angle differences among the tagged jets. We used a single BFKL ladder approach, with inclusive production of two forward/backward and two further, more central, tagged jets. The dependence on the transverse momenta and rapidities of the two central jets is a distinct signal of BFKL dynamics. For future works, more accurate analysis are needed: introduction of parton distribution functions, higher order effect and realistic experimental cuts. It is also pressing to calculate our proposed observables using other approaches not based on the BFKL approach and to test how they can differ from our predictions. Last but not least, we encourage our experimental colleagues to analyze these observables in recent and future LHC data.