Photon–jet correlations in deep-inelastic scattering

The reaction e+p→γ+jet+X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e + p \rightarrow \gamma + \mathrm{jet} + X$$\end{document} is studied in QCD at the next-to-leading order. Previous studies on inclusive distributions showed good agreement with ZEUS data. To obtain a finer understanding of the dynamics of the reaction, several correlation functions are evaluated for ZEUS kinematics.


Abstract
The reaction e + p → γ + jet + X is studied in QCD at the next-to-leading order. Previous studies on inclusive distributions showed good agreement with ZEUS data. To obtain a finer understanding of the dynamics of the reaction, several correlation functions are evaluated for ZEUS kinematics.
This note is the continuation of the study of the deep inelastic scattering (DIS) reaction e + p → γ + jet + X with a photon and a jet in the final state. In a preceding paper [1] we calculated the next-to-leading order (NLO) QCD cross sections describing the photoproduction of a largep ⊥ photon ( p γ ⊥ , η γ ) accompanied by a jet. In particular, we calculated the transverse momentum and rapidity distributions of the photon with the jet constrained in a given large range as well as the jet transverse momentum spectrum in events with a detected photon. In the present note we concentrate on the distributions in the following correlation variables: the laboratory frame E e , E p are respectively the incident energies of the electron and the proton and y is the DIS inelasticity), η = η γ − η jet and φ = φ electron −φ γ . Note thatx p , as defined above, does not correspond to x p , the fraction of momentum carried by a parton in the proton, since we are not in a collinear frame. The correlations should provide a finer understanding of the underlying production mechanism than the inclusive cross sections.
Originally the NLO calculation [2,3] was performed in the γ * − p center of mass (γ * is the virtual photon) in which a large scale is provided by a large value of p * γ ⊥ , the final γ transverse momentum in that frame. Production cross sections of hadrons and jets were studied in this frame and successfully described by the NLO calculation [4][5][6]. On the other hand, in the ZEUS experiment [7] the reaction e + p → γ + jet + X is studied in the laboratory frame for a e-mail: patrick.aurenche@lapth.cnrs.fr b e-mail: michel.fontannaz@th.u-psud.fr which the original NLO calculation must be adapted. Indeed a large p γ ⊥ in the laboratory does not necessarily correspond to a large p * γ ⊥ in the γ * − p frame. Therefore a cut-off E * ⊥cut must be introduced in the calculation with p * γ ⊥ > E * ⊥cut in order to remain in a perturbative domain. In the preceding paper we considered the values E * ⊥cut = 2.5 GeV and E * ⊥cut = 0.5 GeV. Good agreement was found between theory and the ZEUS data with the cut-off E * ⊥cut = 2.5 GeV. We also found that there are kinematical domains in which the NLO cross sections are little dependent on E * ⊥cut , for instance for small ratios (Q/ p γ ⊥ ) 2 or (Q/ p jet ⊥ ) 2 . This result demonstrates that the theoretical calculations are almost model independent on the corresponding domains. A detailed discussion of this problem can be found in Ref. [1].
In this note we pursue the study of the reaction e + p → γ + jet + X by calculating the cross sections dσ/dx γ , dσ/dxx p , dσ/d η and dσ/d φ with the cut-offs E * ⊥cut = 2.5 GeV and E * ⊥cut = 0.5 GeV, and, following ZEUS [8], we consider two Q 2 -ranges, namely 10 < Q 2 < 350 GeV 2 and 10 < Q 2 < 30 GeV 2 . As in Ref. [1] we adopt the ZEUS kinematics [7] with √ s = 319 GeV; the photon momentum has to lie in the ranges 4 < p γ ⊥ < 15 GeV and −0.7 < η γ < 0.9. For the jet momentum we have 2.5 < p jet ⊥ < 35 GeV and −1.5 < η jet < 1.8. Constraints on the final electron are E e > 10 GeV, θ e > 140 • (the z axis is pointing toward the proton direction). As in the ZEUS experiment, the photon is isolated using the democratic k ⊥ -algorithm [9] where the photon is treated on the same footing as partons. We refer to Ref. [1] for details on the isolation criteria and other parameters and conventions used in the calculation. We do not consider here the emission of photons by the electron. Therefore, this cross section should be added to those calculated in this paper to reconstruct the experimental cross section measured by the ZEUS collaboration.
We recall that the cross section is the sum of four building blocks. The "direct" component where the initial virtual photon is coupled to the hard partonic process and the "resolved" component where it interacts via its structure function. In each of these the final real photon couples directly to the hard process or is a fragment of a parton produced at large transverse momentum (fragmentation component). Due to the photon isolation criteria the fragmentation components are small, typically 10-20% of the total depending on the kinematics. We use the virtual photon structure function presented in Ref. [10] and the CTEQ6M parton distributions in the proton [11]. The fragmentation function is that of the BFG collaboration (set II) [12]. All the components are calculated at NLO and details can be found in Ref. [13]. We work in the MS scheme for factorisation and renormalisation with The shapes of the curves corresponding to the small domain 10 < Q 2 < 30 GeV 2 (red curves) and to 10 < Q 2 < 350 GeV 2 are similar, but the curves of the small Q 2 -domain are a factor 2-3 smaller than those of the full Q 2 range. For each Q 2 domain two values of E * ⊥cut have been used, E * ⊥cut = .5 GeV and E * ⊥cut = 2.5 GeV. It clearly appears that the low Q 2 data are globally less sensitive to the transverse momentum cut-off E * ⊥cut : indeed, in the laboratory, the transverse momentum of the virtual photon is Q 2 (1 − y) and, therefore, at small Q 2 the values of p γ ⊥ and For the large Q 2 case the effect of the E * ⊥cut -variation is irregular and strongly depends on the values of the kinematical variables. Therefore the shapes of the distributions are sensitive to E * ⊥cut . But there are kinematical domains in which the distributions are little affected by E * ⊥cut . For instance, dσ/dx γ is rather independent on E * ⊥cut for x γ 0.7. This can be understood in the following way: a small value of x γ corresponds to a large value of y and consequently a small transverse momentum of the initial virtual photon. Also the parton with a small x γ 1 carries away a small transverse momentum. These two effects explain the weak sensitivity of dσ/dx γ to E * ⊥cut for x γ 0.7. Like- 1 The Born direct contribution corresponds to x γ = 1.
wise, the positive η = η γ − η jet domain is rather stable: in this domain the cross section is dominated by the direct processes, less sensitive to E * ⊥cut than the resolved contributions. As for dσ/dx p a similar reason leads to a stability region for largex p . A stable behavior is observable for dσ/d φ in the domain φ < 140 • , which can be explained by kinematical reasons. At φ = 0, for instance, p * γ ⊥ is always larger than p γ ⊥ > 4 GeV and an E * ⊥cut is not necessary. When there is a large sensitivity, we found in Ref. [1] that inclusive data were well described with the choice E * ⊥cut = 2.5 GeV. Concerning the sensitivity of the cross sections to the factorisation and renormalisation scales, we have studied it for one representative point of the distributions in the case of 10 < Q 2 < 30 GeV 2 . For the renormalisation and the proton factorisation we use μ = M = c Q 2 + p * γ ⊥ 2 while in the virtual photon structure function we take M F = Q 2 + (c p * γ ⊥ ) 2 . This last choice is natural since M F cannot be smaller than Q 2 in the photon structure function [10]. The sample points are 0.6 < x γ < 0.7, 0.6 < η < 1.3 and 80 • < φ < 110 • , domains in which the cross sections are almost insensitive to E * ⊥cut . We observe that in the range 0.5 < c < 2. the cross section varies by at most 15%. For example, for the distribution in x γ this variation corresponds to the thickness of the line in Fig. 1, while for the distribution in η γ − η jet it is displayed on Fig. 2, at the value η γ − η jet = 0.95. We notice that this relative stability is the result of a huge compensation between the direct and the resolved terms whose size varies by a factor 2 to 4 under the scales changes.
Finally, concerning the accuracy of the theoretical calculations, let us discuss another point besides the sensitivity to the cut E * ⊥cut and to the scale variations. When calculating the higher order corrections to the resolved component, we assume that the photon structure function, proportional to ln p * 2 ⊥ +Q 2 Q 2 , is large, so that we can neglect terms not proportional to this logarithm. We have no way to estimate the importance of such terms in the ZEUS experiment in which the condition Q 2 / p * 2 ⊥ 1 is not always verified. To test the validity of the approximation it would be interesting to compare theory and data in the two ranges 10 < Q 2 < 30 GeV 2 and 30 < Q 2 < 350 GeV 2 [15].
In conclusion, in the ZEUS experimental configuration for observables constructed with the laboratory kinematics we have a good theoretical stability under changes of cutoffs in a relatively large domain of the variables x γ < 0.7, η > 0 and φ < 130 • when 10 < Q 2 < 350 GeV 2 . Outside these ranges the physics becomes sensitive to nonperturbative effects. In the small Q 2 domain the QCD predictions are independent on the cut-offs in the above quoted ranges and rather insensitive elsewhere. As for the scale uncertainties they are less than ±8% when probing the standard scale ranges. On comparison of these NLO predictions with preliminary ZEUS data there is good agreement after the photon emission from the electron is taken into account [15].