Hadron-jet correlations in high-energy hadronic collisions at the LHC

The inclusive production at the LHC of a charged light hadron and of a jet, featuring a wide separation in rapidity, is suggested as a new probe process for the investigation of the BFKL mechanism of resummation of energy logarithms in the QCD perturbative series. We present some predictions, tailored on the CMS and CASTOR acceptances, for the cross section averaged over the azimuthal angle between the identified jet and hadron and for azimuthal correlations.


Introduction
The LHC record energy, as well as the good resolution in azimuthal angles of the particle detectors, offer a unique opportunity to test a wide class of predictions of perturbative QCD. These include the so called Mueller-Navelet jet production [1], i.e. the inclusive production of two jets featuring a large rapidity separation between them, for which a wealth of theoretical analyses were produced in the last years [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16], and the somewhat related process where two identified, charged light hadrons well separated in rapidity are inclusively produced [17,18], instead of jets.
A common feature to these two processes is that their highenergy behavior is dominated by those final-state configurations where the produced particles are strongly ordered in rapidity, the tagged objects (jets or identified hadrons) being the two extrema in the rapidity tower, thus yielding a number a e-mail: ad.bolognino@unical.it b e-mail: francescogiovanni.celiberto@fis.unical.it c e-mail: d-ivanov@math.nsc.ru d e-mail: maherali8932@gmail.com e e-mail: alessandro.papa@fis.unical.it of energy logarithms growing with the number of produced particles. Such energy logarithms are so large to compensate the smallness of the coupling α s , so the perturbative series must be properly resummed.
The theoretical framework for the resummation of energy logs for these two processes, as well as for any semi-hard process in perturbative QCD, is provided by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach [19][20][21][22], whereby the resummation of all terms proportional to (α s ln(s)) n , the so called leading logarithmic approximation or LLA, and that of all terms proportional to α s (α s ln(s)) n , the next-toleading approximation or NLA, can be systematically carried out. The bottom line of the BFKL formalism is that azimuthal coefficients of the Fourier expansion of the cross section differential in the variables of the tagged objects over the relative azimuthal angle take the very simple form of a convolution between two impact factors, describing the transition from each colliding proton to the respective final state tagged object, and a process-independent Green's function. The BFKL Green's function obeys an integral equation, whose kernel is known at the next-to-leading order (NLO) both for forward scattering (i.e. for t = 0 and color singlet in the t-channel) [23,24] and for any fixed (not growing with energy) momentum transfer t and any possible two-gluon color state in the t-channel [25][26][27][28][29][30][31].
The impact factors for the proton to forward jet transition (the so called "jet vertices") are known up to the NLO for several jet selection algorithms [32][33][34][35][36]. The jet vertex, in its turn, can be expressed, within leading-twist collinear factorization, as the convolution of the parton distribution function (PDF) of the colliding proton, obeying the standard DGLAP evolution [37][38][39], with the hard process describ-ing the transition from the parton emitted by the proton to the forward jet in the final state. Two such jet vertices must be convoluted with the BFKL Green's function to theoretically describe the Mueller-Navelet jet production. The main aim is to calculate cross sections and azimuthal angle correlations [40,41] between the two measured jets, i.e. average values of cos (nφ), where n is an integer and φ is the angle in the azimuthal plane between the direction of one jet and the direction opposite to the other jet, and ratios of two such cosines [42,43].
Also the impact factors for the proton to identified hadron transition are known up to the NLO [44] and can be expressed, within leading-twist collinear factorization, as the convolution of the parton distribution function (PDF) of the colliding proton with the hard process describing the transition from the parton emitted by the proton to a finalstate parton and with the fragmentation function (FF) for that parton to the desired hadron. Two such hadron vertices must be convoluted with the BFKL Green's function to theoretically describe the above-mentioned inclusive hadronhadron production and finally get predictions for cross sections and azimuthal angle correlations, similarly to the case of jets.
Within the same formalism, other interesting processes have been proposed as a testfield for BFKL dynamics at the LHC, namely the inclusive production of three or four jets, well separated in rapidity from each other [45][46][47][48][49], the inclusive detection of two heavy quark-antiquark pairs, separated in rapidity, in the collision of two real (or quasi-real) photons [50], and the inclusive tag of a forward J/ -meson and a very backward jet at the LHC [51].
On the experimental side the situation is as follows: the CMS Collaboration [52] has presented the first measurements of the azimuthal correlation of the Mueller-Navelet jets at √ s = 7 TeV, but further experimental studies of the Mueller-Navelet jets are expected at higher LHC energies and larger rapidity intervals, including also the effects of using asymmetrical cuts for the jet transverse momenta. No experimental analyses have yet appeared on azimuthal correlation between two rapidity-separated identified light hadrons. The reason for that could be that events with identified hadrons in the final state, carrying transverse momenta of the order of, say, 5 GeV or larger, fall into the class of what experimentalists call "minimum bias events", which represent the main background in high-luminosity runs at a collider. They would be better studied in low-luminosity, dedicated, runs.
In this paper we want to introduce a new process which could serve as a probe of BFKL dynamics: the inclusive hadron-jet production in proton-proton collisions, when a charged light hadron: π ± , K ± , p (p) and a jet with high transverse momenta, separated by a large interval of rapidity, are produced together with an undetected hadronic system X (see Fig. 1 for a schematic view). The process (1) has many common features with the inclusive J/meson plus backward jet production, considered recently in Ref. [51]. From the experimental side, the detection of the J/ -meson looks rather appealing. But, from the theory side, there are more uncertainties in this case in comparison to our proposal. The J/ -meson production impact factor was considered in LO; moreover, several production mechanisms in the frame of NRQCD were discussed. Instead, the light hadron impact factor is well defined in collinear factorization and it is known in NLO. Previous experience in BFKL calculations for various processes at LHC shows that the account of NLO corrections to the impact factors leads both to a considerable change of predictions and to a big reduction of the theoretical uncertainties.
The theoretical task to build predictions for cross section and azimuthal correlations for our process is embarrassingly simple: one should simply replace one of the two jet impact factor entering the Mueller-Navelet formulas with the vertex for the proton-to-hadron transition. From the theoretical point of view, this process is definitely an easy target, since all the needed building blocks are available, with NLO accuracy.
Yet, we believe that there are some good reasons for building numerical predictions for this process and submitting them to the attention of both experimentalists and theorists: • the BFKL resummation implies certain factorization structure for the predicted observables: the latter are calculated as a convolution of the universal BFKL Green's function with the process dependent impact factors, which resembles the factorization in Regge theory. It is important to test this picture experimentally, considering all possible processes for which the full NLO BFKL description is available. • In Refs. [9,10,13,53] it was discussed, in the context of Mueller-Navelet jet production, that using asymmetric cuts for the transverse momenta of the tagged jets suppresses the Born term, present only for back-to-back jets, thus enhancing the effects of the additional undetected hard gluon radiation and making therefore more visible the impact of the BFKL resummation, with respect to the fixed-order (DGLAP) contribution. For the process we are considering here this kind asymmetry would be naturally imposed by the completely different nature of the two tagged objects: the identified jet should have transverse momentum not smaller than 20 GeV or so, whereas the minimum hadron transverse momentum can be as small as 5 GeV.

Fig. 1 Inclusive hadroproduction of a charged light hadron and of a jet
• For the process under consideration only one hadron in the final state should be identified, instead of two as in the hadron-hadron inclusive production, the other identified object being a jet with a typically much larger transverse momentum. This should facilitate the mining of these events out of the minimum-bias ones. • From the theoretical point of view one can use this process to compare models for FFs or for jet algorithms, handling expressions which are linear in the corresponding functions and not quadratic as it would be, respectively, in the hadron-hadron and in the Mueller-Navelet jet case.
The summary of the paper is as follows: in Sect. 1 we present the theoretical framework and sketch the derivation of our predictions; in Sect. 2 we show and discuss the results of our numerical analysis; finally, in Sect. 3, we draw our conclusions and give some outlook.

Theoretical framework
The final state configuration of the inclusive process under consideration is schematically represented in Fig. 1, where a charged light hadron (k H , y H ) and a jet (k J , y J ) are detected, featuring a large rapidity separation, together with an undetected system of hadrons. For the sake of definiteness, we will consider the case where the hadron rapidity y H is larger than the jet one y J , so that Y ≡ y H − y J is always positive. This implies that, for most of the considered values of Y , the hadron is forward and the jet is backward.
The hadron and the jet are also required to possess large transverse momenta, k 2 H ∼ k 2 J 2 QCD . The protons' momenta p 1 and p 2 are taken as Sudakov vectors satisfying p 2 1 = p 2 2 = 0 and 2( p 1 p 2 ) = s, so that the momenta of the final-state objects can be decomposed as In the center-of-mass system, the hadron/jet longitudinal momentum fractions x H,J are connected to the respective rapidities through the relations y H = 1 2 ln , here the space part of the four-vector p 1 being taken positive.
In QCD collinear factorization the cross section of the process (1) reads where the r, s indices specify the parton types (quarks q = u, d, s, c, b; antiquarksq =ū,d,s,c,b; or gluon g), f r,s (x, μ F ) denote the initial proton PDFs; x 1,2 are the longitudinal fractions of the partons involved in the hard subprocess, while μ F is the factorization scale; dσ r,s ŝ is the partonic cross section andŝ ≡ x 1 x 2 s is the squared centerof-mass energy of the parton-parton collision subprocess.
In the BFKL approach the cross section can be presented (see Ref. [4] for the details of the derivation) as the Fourier sum of the azimuthal coefficients C n , having so with φ H,J the hadron/jet azimuthal angles, while y H,J and k H,J are their rapidities and transverse momenta, respectively. The φ-averaged cross section C 0 and the other coefficients C n =0 are given by 1 1 In Ref. [18], on the last line of Eq. (5), which is closely related to this formula for C n , it was mistakenly written 2 ln k 2 1 k 2 2 instead of ln k 2 1 k 2 2 , although the numerical results presented there were obtained using the correct formula. Hereᾱ is the first coefficient of the QCD β-function, where n f is the number of active flavors, is the leading-order (LO) BFKL characteristic function, c H (n, ν) is the LO forward hadron impact factor in the νrepresentation, given as an integral in the parton fraction x, containing the PDFs of the gluon and of the different quark/antiquark flavors in the proton, and the FFs of the detected hadron, and the f (ν) function is defined by The remaining objects are the hadron/jet NLO impact factor corrections in the ν-representation, c H,J (n, ν, | k H,J |, x H,J ), their expressions being given in Eqs. (4.58)-(4.65) of Ref. [44] and in Eq. (36) of Ref. [4], respectively.

Integration over the final-state phase space
In order to match the actual LHC kinematic cuts, we integrate the coefficients over the phase space for two final-state objects and keep fixed the rapidity interval, Y , between the hadron and the jet: We consider two distinct ranges for the final-state objects:  [52]. For the sake of brevity, we will refer to this choice as the CMS-jet configuration; • a hadron always detected inside CMS in the range given above, together with a very backward jet tagged by CAS-TOR. In this peculiar, CASTOR-jet configuration, the jet lies in the typical range of the CASTOR experimental analyses, i.e. k min 17.68 GeV for √ s = 13 TeV (CASTOR-jet). The rapidity interval, Y , is taken to be positive: 0 < Y ≤ y max H − y min J . Two center-of-mass energies, √ s = 7 and 13 TeV, are taken into account in the CMS-jet configuration, while we give predictions for √ s = 13 TeV in the CASTORjet case.
In our calculations we use the MMHT 2014 NLO PDF set [55] with two different NLO parametrizations for hadron for a related discussion). In the results presented below, we sum over the production of forward charged light hadrons: π ± , K ± , p,p.

Scale optimization
To fix the renormalization scale μ R , which can be arbitrarily chosen within the NLA, we adopt the BLM [58][59][60][61] approach, which has become a quite common choice for semihard processes. We first perform a finite renormalization from the MS to the physical MOM scheme, whose definition is related to the 3-gluon vertex being a key ingredient of the BFKL approach and get with T = T β + T conf , 2.3439 and ξ is the gauge parameter of the MOM scheme, fixed at zero in the following. Then, the "optimal" BLM scale μ BLM R is the value of μ R that makes the β 0 -dependent part in the expression for the observable of interest vanish. In Ref. [12] some of us showed that terms proportional to the QCD β 0 -function are present not only in the NLA BFKL kernel, but also in the expressions for the NLA impact factor. This leads to a non-universality of the BLM scale and to its dependence on the energy of the process.
Finally, the condition for the BLM scale setting was found to be The term in the r.h.s. of Eq. (14) (proportional to α MOM s ) originates from the NLA part of the kernel, while the remaining ones come from the NLA corrections to the hadron/jet vertices.
In order to find the values of the BLM scales, we introduce the ratios of the BLM to the "natural" scale suggested by the kinematic of the process, μ N = | k H || k J |, so that m R = μ BLM R /μ N , and look for the values of m R which solve Eq. (14).
We finally plug these scales into our expression for the integrated coefficients in the BLM scheme (for the derivation see Ref. [12]): The coefficient C 0 gives the φ-averaged cross section, while the ratios R n0 ≡ C n /C 0 = cos(nφ) determine the values of the mean cosines, or azimuthal correlations, of the produced hadron and jet. In Eq. (15),χ(n, ν) is the eigenvalue of NLA BFKL kernel [62] and its expression is given, e.g. in Eq. (23) H,J are the NLA parts of the hadron/jet vertices (see Ref. [12]).
We set the factorization scale μ F equal to the renormalization scale μ R , as assumed by the MMHT 2014 PDF.
All calculations are done in the MOM scheme. For comparison, we present results for the φ-averaged cross section C 0 in the MS scheme, as implemented in Eq. (11).
In the latter case, we choose natural values for μ R , i.e. μ R = μ N ≡ | k H || k J |, and two different values of the factorization scale, (μ F ) 1,2 = | k H,J |, depending on which of the two vertices is considered. We checked that the effect of using natural values also for μ F , i.e. μ F = μ N , is negligible with respect to our two-value choice.

Used tools and uncertainty estimation
All numerical calculations were done using JETHAD, a Fortran code we recently developed, suited for the computation of cross sections and related observables for two-body final-state processes, and offering also support in the study of multi-body final-state reactions. In order to perform numerical integrations, we interfaced JETHAD with specific CERN program libraries [63] and with Cuba library integrators [64,65]. We made extensive use of the CERNLIB  The most relevant uncertainty comes from the numerical 4-dimensional integration over the two transverse momenta | k H,J |, the hadron rapidity y H , and over ν. Its effect was directly estimated by Dadmul integration routine [63]. The other three sources of uncertainty, which are respectively: the one-dimensional integration over the parton fraction x needed to perform the convolution between PDFs and FFs in the LO/NLO hadron impact factors, the one-dimensional integration over the longitudinal momentum fraction ζ in the NLO hadron/jet impact factor corrections, and the upper cutoff in the numerical integrations over | k H,J | and ν, are negligible with respect to the first one. For this reason the error bands of all predictions presented in this work are just those given by the Dadmul routine.

Discussion
In Fig. 2 we present our results at natural scales for the φaveraged cross section C 0 at √ s = 7 and 13 TeV in the CMSjet kinematic configuration. We can see that the NLO corrections become larger and larger at increasing Y , an expected phenomenon in the BFKL approach.
In Figs. 3 and 4, predictions with the BLM scale optimization for C 0 and several R nm ≡ C n /C m ratios with the jet tagged inside the CMS detector are shown for √ s = 7 and 13 TeV, respectively. Here the benefit of the use of BLM optimization appears, since the LLA and NLA predictions for C 0 are now comparable, a sign of stabilization of the perturbative series. The trend of ratios of the form R n0 is the standard one and indicates increasing azimuthal decorrelation between the jet and the hadron as Y goes up, with the NLA predictions systematically above the LLA ones, as it was also observed in Mueller-Navelet jets and in the hadron-hadron case. The ratios R 21 and R 32 seem to be almost insensitive to the NLO corrections.
Panels in Fig. 5 show results with BLM scale optimization for C 0 and several R nm ratios in the CASTOR-jet configuration at √ s = 13 TeV. They exhibit some new and, to some extent, unexpected features: (i) the two parametrizations for the FFs lead to clearly distinct predictions, (ii) cos φ exceeds one at the smaller values for Y , a clearly unphysical effect. The reason for these phenomena could reside in the fact that, the lower values for Y in the CASTORjet case are obtained for negative values of the hadron rapidity, i.e. in final-state configurations where both jet and hadron are backward.
Finally, in Fig. 6 we compare the φ-averaged cross section C 0 in different NLA BFKL processes: Mueller-Navelet jet, hadron-jet and hadron-hadron production, for μ F = μ BLM R , at √ s = 7 and 13 TeV, and Y ≤ 7.1 in the CMS-jet case. The hadron-hadron cross section, with the kinematical cuts adopted, dominates over the jet-jet one by an order of magnitude, with the hadron-jet cross section lying, not surprisingly, in-between.

Summary
In this paper we have proposed a new candidate probe of BFKL dynamics at the LHC in the process for the inclusive production of an identified charged light hadron and a jet, separated by a large rapidity gap.
We have given some arguments that this process, though being a naive hybridization of two already well studied ones, presents some own characteristics which can make it worthy of consideration in future analyses at the LHC.
In view of that, we have provided some theoretical predictions, with next-to-leading accuracy, for the cross section averaged over the azimuthal angle between the identified jet and hadron and for ratios of the azimuthal coefficients.
The trends observed in the distributions over the rapidity interval between the jet and the hadron are not different from the cases of Mueller-Navelet jets and hadron-hadron, when the jet is detected by CMS, whereas some new features have appeared when the jet is seen by CASTOR, which deserve further investigation.