A fresh look at factorization breaking in diffractive photoproduction of dijets at HERA at next-to-leading order QCD

We calculate the cross section of diffractive dijet photoproduction in ep scattering at next-to-leading order (NLO) of perturbative QCD (pQCD), which we supplement by a model of factorization breaking for the resolved-photon contribution. In this model, the suppression depends on the flavor and momentum fraction of the partons in the photon. We show that within experimental and theoretical uncertainties, the resulting approach provides a good description of the available HERA data in most of the bins. Hence, taken together with the observation that NLO pQCD explains well the data on diffractive photoproduction of open charm in ep scattering, our model of factorization breaking presents a viable alternative to the scheme based on the global suppression factor.


Introduction
One of the highlights of the physics results obtained at HERA was the measurement of inclusive diffraction γ * p → X p in lepton-proton deep-inelastic scattering (DIS). Combined analyses of the H1 and ZEUS experiments have been published in Ref. [1]. Contrary to the expectations that the probability of large rapidity gap events in DIS should be very small [2], it was found that diffractive events in DIS constitute approximately 10-15 % of the total cross section over a wide range of Q 2 . Furthermore, the QCD collinear factorization theorem for hard inclusive diffraction [3] allowed one to treat diffractive DIS on the same footing as inclusive DIS. One can thus first introduce universal diffractive parton distribution functions (dPDFs) and then determine them by fitting to the measured diffractive structure function [4][5][6]. The universality of the resulting dPDFs is confirmed by a e-mail: vguzey@pnpi.spb.ru b e-mail: michael.klasen@uni-muenster.de good agreement between the perturbative QCD (pQCD) calculations and the data on diffractive production of dijets [7] and open charm [8] in DIS.
At the same time, in diffractive photoproduction of dijets in ep scattering, based on the well-known factorization breaking in diffractive dijet production inp p collisions at the Tevatron [9][10][11], collinear factorization is generally not expected to hold [12][13][14]. However, the pattern of this factorization breaking remains an open question [15]; for recent reviews see, e.g., Refs. [16][17][18]. To summarize, while the recent H1 [19,20] and ZEUS [21] data on diffractive dijet photoproduction are largely consistent with each other after the renormalization of the ZEUS data, the next-to-leading order (NLO) perturbative QCD calculations overestimate the data by approximately 40-50 %. The theory and the data can be made consistent by introducing either a global suppression factor of 0.5 or a suppression factor of approximately 0.4 only for the resolved-photon contribution. The most recent H1 measurement of diffractive dijet photoproduction with a leading proton [22] is also consistent with the observation that NLO pQCD globally overestimates the data by the factor of 0.5.
Factorization breaking in diffractive dijet photoproduction is a result of soft inelastic photon interactions with the proton, which populate and thus partially destroy the finalstate rapidity gap. This effect is usually described in the literature by a rapidity gap survival factor S 2 ≤ 1. Since the magnitude of S 2 decreases with an increase of the interaction strength between the probe and the target, the pattern of the factorization breaking can be related to various components of the photon [17]. In the laboratory reference frame, the high-energy photon interacts with hadronic targets by fluctuating into various configurations (components) interacting with the target with different cross sections. These fluctuations contain both weakly interacting (the so-called point-like) components and the components interacting with large cross sections, which are of the order of the vector meson-proton cross sections. This general spacetime picture of photon-hadron interactions at high energies is usually realized in the framework of such approaches as the vector meson dominance (VMD) model and its generalizations [23,24] or the color dipole model [25,26]. It is also used in the language of collinear factorization, where the photon structure function and parton distribution functions (PDFs) are given by a sum of the resolved-photon contribution corresponding to the VMD part of the photon wave function and the point-like (inhomogeneous) term originating from the γ → qq splitting; see, e.g., Ref. [27]. Note that the direct-photon contribution to diffractive dijet photoproduction corresponds to the configurations interacting with very small cross sections of the order of 1/E 2 T (E T is the transverse jet energy), which preserves factorization.
Let us recall that the hadron (VMD) part of the photon PDFs contributes only for small x γ , whereas the point-like term gives the dominant contribution for large x γ . Here, x γ is the light-cone momentum fraction of a parton in the photon. Based on the arguments presented above, it is then natural to expect S 2 = 1 for the direct-photon contribution localized near x γ = 1, S 2 ≈ 0.34 for the hadron-like component of the photon at small x γ , and S 2 ≈ 0.53-0.75 for the gluon and quark contributions at large x γ corresponding to small, but non-negligible factorization breaking due to the point-like component of the resolved photon [17,18]. Note that in pQCD the separation of the direct and the resolvedphoton contributions is unambiguous only at leading order. At NLO, it becomes a matter of convention depending on the factorization scheme and the factorization scale [28][29][30].
Another important observation relevant to the possible pattern of factorization breaking is that the HERA data on diffractive photoproduction of open charm [31] agree with the NLO pQCD calculations [32,33], and hence no factorization breaking is required. This calls into question the assumption of a global suppression factor modeling factorization breaking and indicates that S 2 for the resolved-photon contribution may depend on the parton flavor. In particular, S 2 ≈ 1 for the charm quark distribution in the photon, which agrees with the observation that in the VMD model the J/ψ-proton cross section is of the order of a few mbarn, i.e. much smaller than the ρ-nucleon cross section [23,24].
In the present work we revisit the issue of factorization breaking in diffractive dijet photoproduction in ep scattering and perform NLO pQCD calculations of the corresponding cross sections, which we combine with a new flavor-dependent and momentum fraction-dependent scheme of factorization breaking for the resolved-photon contribution. We demonstrate that the results of our calculations provide a good description of the H1 [19,20] and ZEUS [21] data on diffractive dijet photoproduction in ep scattering at HERA, while simultaneously, by construction, not conflicting with the good pQCD description of diffractive photoproduction of open charm in ep scattering. Thus, next-to-leading order perturbative QCD coupled with the physically motivated assumption as regards the rapidity gap survival probability for the resolved-photon contribution and the effect of hadronization corrections provide a good description of all available HERA data on diffractive dijet photoproduction. This result reinforces the conclusion of Ref. [17].

New scenario for factorization breaking in diffractive photoproduction of dijets in ep
We performed next-to-leading order (NLO) pQCD calculations [18] of the cross sections of diffractive photoproduction of dijets in ep scattering ep → e + 2jets + X + Y using the kinematic conditions and cuts of the H1 [19,20] and ZEUS [21] measurements of this process. This process is illustrated in Fig. 1. As is well known, jet production with real photons involves direct interactions of the photon with quarks or gluons from the proton (or in our case from the pomeron) as well as resolved-photon contributions, leading to parton-parton interactions and an additional remnant jet coming from the photon. For the direct interactions, factorization is expected to be valid as in the case of DIS, whereas we expect it to fail for the resolved process as in hadron-hadron scattering. For this part of photoproduction one would therefore naively expect a similar suppression factor due to rescattering effects of the hadronic fluctuations of the photon. The expression for the cross section reads where X denotes the pomeron (and possibly photon) remnant jet(s); Y denotes either a proton or a low-mass proton excitation; S 2 i (x γ ) is the factor modeling factorization breaking; f γ /e (y) is the photon flux of the electron depending on the photon light-cone momentum fraction y; f i/γ (x γ , μ 2 ) is the PDF of parton i in the photon depending on the momentum fraction x γ and the factorization scale μ; f D (4) j/ p (x IP , z IP , t, μ 2 ) is the diffractive PDF of the proton, which depends on x IP (the momentum fraction carried by the diffractive exchange or "Pomeron") and z IP (the momen- The same as in Fig. 2, but for the H1 data with the high-E jet T cut [20] scheme [27] and the 2006 H1 proton diffractive PDFs (fit B) [5]. We explained in the Introduction that the space-time picture of high-energy photon-proton interactions suggests that in diffractive dijet photoproduction on the proton, QCD diffractive factorization holds for the direct-photon contribution and is broken for the resolved-photon contribution. Moreover, for the latter contribution, the factorization breaking is strongest at small x γ , small, but non-negligible, for large x γ , and depends on the parton flavor. In the framework of collinear QCD factorization, the decrease of factorizationbreaking effects in the resolved-photon contribution with an increase of x γ can be explained by the observation that based on the factorization of the collinear singularity, there should be a smooth transition from the resolvedphoton contribution at large x γ to the direct-photon contribution. Therefore, we model the effect of factorization breaking by introducing the following suppression factor of S 2 i (x γ ) for the resolved-photon contribution (i.e. for the photon PDFs) in Eq. (1): where i is the parton flavor; A q = 0.37-0.41 and A g = 0.19-0.24. The given ranges of values take into account the possible effective dependence of S 2 i (x γ ) on the hard resolution scale, where the first and the second values correspond to E The same as in Fig. 2, but for the ZEUS data [21] cuts of the experimental analyses are summarized in Table 1, where Q 2 refers to the photon virtuality, y its momentum fraction in the electron, E jet1(2) T are the leading and subleading transverse jet energies, and η jet1(2) their rapidities. In Figs. 2, 3, and 4, the thick red solid lines correspond to the calculation, when the renormalization and factorization scale μ is identified with the transverse energy of jet 1, μ = E jet1 T . The thin red dotted lines quantify the scale uncertainty of our NLO calculations and correspond to μ = 2E jet1 T (lower) and μ = E jet1 T /2 (upper). For comparison, we also show the unsuppressed predictions (assuming no factorization breaking) by the blue dot-dashed lines labeled "NLO, R = 1". Note that our theoretical calculations have been multiplied by the hadronization corrections in each bin [19][20][21]. In the different panels, the values of the cross section are shown as functions of the following variables; see, e.g., Ref. [19]: is the hadron-level estimator of the parton momentum fraction in the photon, where the sum runs over the hadronic final states i included in the jets and X refers to the full diffractive final state; z jets IP = jets (E i + P i,z )/(E X + P X,z ) is the estimator of the "Pomeron" momentum fraction carried by a parton; x IP = (E X + P X,z )/(2E beam p ) is the measured "Pomeron" momentum fraction, where E beam p is the proton beam energy; E jet1 T is the transverse energy of jet 1; M X is the invariant mass of the diffractive final state (two jets plus diffractively produced remnants of the photon and the "Pomeron"); M 12 is the invariant mass of the dijet system; η jets = (η 1 + η 2 )/2 and η jets | = |η 1 − η 2 | are the average and the relative jet rapidities, and W is the invariant photonproton energy.
One can see from Figs. 2, 3, and 4 that within theoretical uncertainties, NLO perturbative QCD combined with the model of factorization breaking of Eq. (2) provides a good description of the data for most of the bins. The quality of the description of the data is similar to that of Ref. [18], where factorization breaking is realized either by the global or the resolved-only, flavor-and x-independent suppression factors.
An inspection of Figs. 2, 3, and 4 shows that NLO pQCD correctly reproduces the shape of almost all considered distributions and only fails to explain the normaliza-  tion in some bins receiving a significant contribution from the direct-photon contribution, which is unsuppressed in our factorization-breaking scheme. The most notable example is the distribution at high x jets γ at low E jet T . It was hypothesized in Ref. [17] that additional sizable hadronization corrections along with bin migration effects, which are not included in our analysis, might help to improve the agreement between theory and data at large x jets γ . This hypothesis is supported by very similar observations in inclusive photoproduction at low E jet T [30]. Integrating our results for dσ/dE jet1 T over E jet1 T , we obtain the theoretical prediction for the integrated cross section σ tot NLO . Table 2 gives our results for σ tot NLO , the corresponding experimental cross sections σ tot data , and the ratios of the theoretical predictions to the measured values. In the presented values for R, we have added in quadrature the experimental and theoretical uncertainties. One can see from the results shown in the table that within combined experimental and theoretical uncertainties, NLO pQCD with our factorization-breaking scheme gives a good description of the integrated cross section of diffractive dijet photoproduction in ep scattering measured at HERA. Note that the central value of R for the high-E T cut data is somewhat larger than that for the low-E T cut, which can be explained by the fact that the unsuppressed charm quark con-tribution becomes more prominent due to the QCD evolution of the photon PDFs.

Conclusions
We calculated the cross sections of diffractive dijet photoproduction in ep scattering in HERA kinematics using NLO perturbative QCD and a scenario of factorization breaking, which assumes that only the resolved-photon contribution is suppressed. The suppression depended on the parton flavor and the light-cone momentum fraction of partons in the photon. It was absent for charm quarks, larger for gluons than for light quarks, and decreased with an increase of the parton momentum. This model for factorization breaking in diffractive QCD is based on the space-time picture of photon-hadron interactions and complies with the good pQCD description of diffractive photoproduction of open charm in ep scattering. We compared our results with the available H1 and ZEUS data and found that various measured distributions and the integrated cross sections can be reproduced by our calculations with good accuracy. This agreement allows us to advocate our model as a viable alternative to the purely phenomenological scheme based on a global suppression factor. for hospitality. The work of VG is partially supported by a grant of Deutscher Akademischer Austauschdienst (DAAD). The work of MK is partially supported by the BMBF Verbundprojekt 05H2015 through Grant 05H15PMCCA.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .