Inclusive four-jet production at 7 and 13 TeV: Azimuthal profile in multi-Regge kinematics

Recently, new observables in LHC inclusive events with three tagged jets were proposed. Here, we extend that proposal to events with four tagged jets. The events are characterized by one jet in the forward direction, one in the backward direction with a large rapidity distance Y from the first one and two more jets tagged in more central regions of the detector. In our setup, non-tagged associated mini-jet multiplicity is present and needs to be accounted for by the inclusion of BFKL gluon Green functions. The projection of the cross section on azimuthal-angle components opens up the opportunity for defining new ratios of correlation functions of the azimuthal-angle differences among the tagged jets that can be used as probes of the BFKL dynamics.


Introduction
The large hadron collider (LHC) gives a unique opportunity to study high energy scattering in quantum chromodynamics (QCD). Jet production studies play a crucial role in high energy QCD phenomenology since the plethora of data makes possible the analysis of even more exclusive observables than usual. Here, we concentrate on four-jet production in the so-called multi-Regge kinematics. In an experimental setup containing final state jets with a large rapidity separation, the Balitsky-Fadin-Kuraev-Lipatov (BFKL) framework in the leading logarithmic (LL) [1][2][3][4][5][6][7][8][9][10][11] and next-toleading logarithmic (NLL) approximation [12,13] presents itself as a powerful tool to probe the dominant dynamics of the QCD high energy limit.
Arguably, jet production studies in the last decade searching for the onset of BFKL effects were mainly focused on Mueller-Navelet jets (dijets) [14] in hadronic colliders. They D. Gordo Gómez: 'la Caixa'-Severo Ochoa Scholar. a e-mail: chachamis@gmail.com correspond to the inclusive production of two jets with similar transverse sizes, k A,B , and a large rapidity difference Y = ln(x 1 x 2 s/(k A k B )). x 1,2 are the usual Bjorken parameters of the two partons that are linked to the jets and s is the center-of-mass-energy squared. A number of analyses [15][16][17][18][19][20][21] of the average values, cos (m φ) , for the azimuthalangle difference between the two tagged jets, φ, suggest mini-jet activity in the rapidity interval between the most forward and most backward jet which cannot be dismissed and which affects the azimuthal-angle difference of these jets. A downside is that collinear effects [22,23], having their origin at the zeroth component of the conformal spin, affect significantly the azimuthal-angle observables. Nevertheless, this collinear contamination can be mostly eliminated if the ratios of projections on azimuthal-angle observables R m n = cos (m φ) / cos (n φ) [22,23] (where m, n are integers and φ is the azimuthal angle between the two tagged jets) are considered instead. Moreover, the ratios offer a clearer signal of BFKL effects than the standard predictions for the growth of hadron structure functions F 2,L (well fitted within NLL approaches [24,25]). The confrontation of different NLL theoretical predictions for these ratios R m n [26][27][28][29][30][31][32][33][34][35][36][37][38][39] against LHC experimental data has been quite successful.
In Refs. [40,41], we proposed new observables associated to the inclusive production of three jets. We argued there that the new observables feature appealing attributes. First, in order to obtain data for the proposed three-jet analysis, one would only need to search for an additional central jet within any existing Mueller-Navelet analysis data set. Second, on more theoretical grounds, these observables probe fundamental characteristics of the BFKL ladder. They correspond to the ratios being the azimuthal angle of the forward, central,  and backward jet, respectively. Here, we extend our discussion to the case of four-jet events (different experimental analyses can be found in Refs. [42][43][44]). For the present study, we need to have one jet in the forward direction with rapidity Y A , one in the backward direction with rapidity Y B and both well separated in rapidity from the each other so that Y A − Y B is large, along with two extra jets tagged in more central regions of the detector. Additionally, the relative rapidity separation between any two neighboring jets cannot be very different from one third of Y A − Y B so that the kinematical configurations of the events actually follow the multi-Regge kinematics. Extending Eq. (1) to the partonic four-jet production, we studied in [45] different ratios of three cosines in numerator and denominator: where φ 1 , φ 2 , and φ 3 are the azimuthal-angle differences between jets neighboring in rapidity. These ratios allow for the study of even more differential distributions in the transverse momenta, azimuthal angles and rapidities of the two central jets as well as for detailed work in connection to multiple parton scattering [46][47][48][49][50][51][52].
In this paper, we define and study ratios of three cosines in numerator and denominator beyond the partonic level. We make use of the collinear factorization scheme to produce the two uttermost jets and we convolute the partonic differential cross section, which is described by the BFKL dynamics, with collinear parton distribution functions. We also include in our computation the forward "jet vertex" [53][54][55][56][57]. Three BFKL gluon Green functions link the outermost (Mueller-Navelet-like) jets with the more centrally produced ones. We integrate over the momenta of the four produced jets, using LHC kinematical cuts so that a comparison of our predictions with forthcoming experimental analyses of LHC data is possible. In the following section we will overview the main formulas and present our numerical results. We conclude with Sect. 3.

Hadronic level inclusive four-jet production in multi-Regge kinematics
We study (see Figs. 1 and 2) the production of one forward and one backward jet, both characterized by high transverse momenta k A,B and well separated in rapidity, together with two more jets produced in the central rapidity region and with possible associated mini-jet production: The cross section for the inclusive four-jet production process (3) reads in collinear factorization where α, β characterize the partons (gluon g; the parton distribution functions of the protons; x 1,2 represent the longitudinal fractions of the partons involved in the hard subprocess; dσ α,β ŝ, μ F is the partonic cross section for the production of jets andŝ ≡ x 1 x 2 s is the partonic squared center-of-mass energy (see Fig. 1). The cross section for the partonic hard subprocess dσ α,β features a dependence on BFKL dynamics keeping in mind that the emissions of mini-jet in the rapidity span between any two subsequentin-rapidity jets can be described by a forward gluon Green function ϕ.
Making use of the leading order approximation of the jet vertex [53], the cross section for the process (3) reads In order to follow a multi-Regge kinematics setup, we demand that the rapidities of the produced particles obey Y A > y 1 > y 2 > Y B , while k 2 1 and k 2 2 are well above the resolution scale of the detectors. x J A,B are the longitudinal momentum fractions of the two external jets, connected to the respective rapidities Y A,B by the relation The strong coupling isᾱ s = α s (μ R ) N c /π and ϕ are BFKL gluon Green functions following the normalization ϕ ( p, q, 0) = δ (2) ( p − q).
Based upon the work presented in Refs. [41,45], we seek observables for which the BFKL dynamics would surface in a distinct form. Moreover, we are interested in observables that should be rather insensitive to possible higher order corrections. Let us first define the following azimuthal- Then the related experimental observable we propose corresponds to the mean value (with M, N , L being positive integers) Fig. 2 A primitive lego plot depicting a four-jet event. k A is a forward jet with large positive rapidity Y A and azimuthal angle ϑ A , k B is a backward jet with large negative rapidity Y B and azimuthal angle ϑ B and k 1 and k 2 are two jets with azimuthal angles ϑ 1 and ϑ 2 , respectively, and rapidities y 1 and y 2 such

Rapidity
The numerator in Eq. (6) The quantity˜ m,n,l is simply a convolution of the BFKL gluon Green functions, originally defined in [45]: where (ψ is the logarithmic derivative of Euler's gamma function).
As we mentioned previously, we would like to consider quantities that are easily measured experimentally; moreover, we want to eliminate as much as possible any dependence on higher order corrections. Thus, we need to consider ratios 1 similar to Eq. (2) which are defined on a partonic level though. Therefore, in order to provide testable theoretical predictions against any current and forthcoming experimental data, we proceed in two steps. First of all, we impose LHC kinematical cuts by integrating C M N L over the momenta of the tagged jets. More precisely, where the rapidity Y A of the most forward jet k A is restricted to 0 < Y A < 4.7 and the rapidity Y B of the most backward jet k B is restricted to −4.7 < Y B < 0, while their difference Y ≡ Y A − Y B is kept fixed at definite values within the range 6.5 < Y < 9. Obviously, the last condition on the allowed values of Y makes both the integration ranges over Y A and Y B smaller than 4.7 units of rapidity. Second, we remove the zeroth conformal spin contribution responsible for any collinear contamination (contributions that originate at ϕ 0 ) and we minimize possible higher order effects by introducing the ratios To keep things simple, in both cuts, we set k 2 to be larger than all the other three-jet momenta and we only vary the range of k 1 . In the cut defined in Eq. (13), k 1 is smaller that all the other three-jet momenta whereas in the cut defined in Eq. (14), the allowed k 1 values overlap with the ranges of k A and k B . In the plots to follow, we plot the ratios for the cut defined in Eq. (13) with a red dot-dashed line and the ratios for the cut defined in Eq. (14) with a blue dashed line.
The numerical computation of the observables to be shown was done in Fortran. Mathematica was used for various cross-checks. We used the NLO MSTW 2008 PDF sets [58]  Vegas [59] as implemented in the Cuba library [60,61] was our main integration routine. We also made use of the The functional dependence of the ratios R M N L P Q R on the rapidity difference between k A and k B is rather smooth. We can further notice that there are ratios with an almost linear behavior with Y and with a rather small slope. To be specific, the ratios represented by the blue curve in Fig. 3 and the red curve in Figs. 4, 5, and 6 demonstrate this linear behavior in a striking fashion. Furthermore, whenever a ratio exhibits a linear dependence on Y (for a certain kinematical cut of k 1 ) at colliding energy 7 TeV, we observe that the ratio maintains On the other hand, there are configurations for which the functional dependence on Y is much stronger and far from linear. In Fig. 4, the blue curve on the top rises from ∼1.2 at Y = 6.5 to ∼6.8 at Y = 9, whereas in Fig. 6 on the top it drops from ∼(−1.5) to ∼(−4.8) for the same variation in Y . Generally, if for some ratio there is a strong functional dependence on Y for a k 1 of intermediate size (blue curve), this dependence is 'softened' at higher colliding energy (see plots in Figs. 4,5,6,8). However, for a k 1 of smaller size (red curve), we see that the functional dependence on Y gets stronger at 13 TeV (Figs. 3, 7, 8), unless of course it exhibits a linear behavior as was discussed in the previous paragraph.
In all plots presented in Figs. 3, 4, 5, 6, 7 and 8, there is no red or blue curve that changes sign in the interval 6.5 < Y < 9. Moreover, if a ratio R M N L P Q R is positive (negative) at 7 TeV, it will continue being positive (negative) at 13 TeV, disregarding the specific functional behavior on Y .
In contrast to our main observation in Ref. [41] where in general, for most of the observables R M N P Q there were no significant changes after increasing the colliding energy from 7 to 13 TeV, here we notice that, depending on the kinematical cut, an increase in the colliding energy may lead to a noticeable change to the shape of the functional Y dependence, e.g. the red curve in Fig. 3, the blue and red curve in Fig. 8. This is a very interesting point for the following reason. If a BFKLbased analysis for an observable dictates that the latter does not change much when the energy increases, this fact actually indicates that a kind of asymptotics has been reached, e.g. the slope of the gluon Green function plotted as a function of the rapidity for very large rapidities. In asymptotics, the dynamics is driven by pure BFKL effects whereas pre-asymptotic effects are negligible. In the present study, we have a mixed picture. We have ratios that do not really change when the energy increases and other ratios for which a higher colliding energy changes their functional dependence on Y. A crucial point that allows us to speak about pre-asymptotic effects, from which in itself one infers that BFKL is still the relevant dynamics, was outlined previously in this section: despite the fact that for some cases we see a different functional dependence on Y after raising the colliding energy, it is important to note that we observe no change of sign for any ratio R M N L P Q R . Therefore, the four-jet ratio observables we are studying here are more sensitive to pre-asymptotic effects than the related three-jet ratio observables studied in Ref. [41]. Nevertheless, by imposing different kinematical cuts one can change the degree of importance of these effects.
To conclude with, a carefully combined choice of cuts for the R M N L P Q R observables and a detailed confrontation between theoretical predictions and data may turn out to be an excellent way to probe deeper into the BFKL dynamics.

Summary and outlook
We presented a first phenomenological study for some aspects of the jets' azimuthal profile in LHC inclusive fourjet production within the BFKL resummation framework. Following up the work in Ref. [45], where a new set of BFKL probes was proposed for the LHC based on a partonic level study, we have calculated some of these observables here, after convoluting the previous results with parton distribution functions and imposing LHC kinematical cuts, at two different center-of-mass energies, √ s = 7, 13 TeV.
We have chosen an asymmetric kinematical cut with respect to the transverse momenta of the most forward (k A ) and most backward (k B ) jet which is arguably a more interesting kinematical configuration that a symmetric cut, since it allows for an easier distinction between BFKL and fixed order predictions [30,33]. The asymmetry was realized by imposing different lower limits to k A and k B (k min A = 35 GeV and k min B = 45 GeV). Additionally, we demanded k 2 to be larger than both k A and k B whereas the value of the transverse momentum k 1 was allowed to be either smaller than both k A and k B or overlapping the k A and k B range of values.
We have plotted six generalized-azimuthal-ratio observables, R 111 221 , R 112 111 , R 112 211 , R 212 111 , R 122 221 , R 221 112 , as a function of the rapidity distance Y between k A and k B for 6.5 < Y < 9. A smooth functional dependence of the ratios on Y appears to be the rule. It is noteworthy that the plots for ratios we presented exhibit in some cases considerable change when the colliding energy increases from 7 to 13 TeV. This tells us that pre-asymptotic effects do play a role for the azimuthal ratios in inclusive four-jet production. A comparison with predictions for these observables from fixed order analyses as well as from the BFKL inspired Monte Carlo BFKLex [64][65][66][67][68][69][70] seems to be the logical next step. Predictions from multi-purpose Monte Carlo tools should also be pursued.
We will conclude our discussion by stressing that it would be very interesting to have an experimental analysis for these observables using existing and future LHC data. We have the strong belief that such an analysis will be a big step forward to the direction of gauging the applicability, at present energies, of the BFKL dynamics in phenomenological studies.