The photon PDF of the proton

We show how the photon input parton distribution function (PDF) may be calculated with good accuracy and used in an extended DGLAP global parton analysis in which the photon is treated as an additional point-like parton. The uncertainty of the input photon PDF is relatively small, since the major part of the distribution (which is produced by the coherent emission of the photon from a proton that remains intact) is well known. We present the expected photon PDFs and compare the predictions with ZEUS data for isolated photon electroproduction at negative rapidities.


Introduction
Precise parton distribution functions (PDFs) are an essential ingredient in analysing data obtained at high energy hadron colliders, such as the Tevatron and the LHC. In perturbative QCD the PDFs are now known at next-to-next-to leading order. Indeed, with the current precision, it is important to investigate the effect of electroweak corrections in hadron collider physics. In particular, the QED contributions have large logarithmic terms, up to αlog(Q 2 /m 2 q ), arising from photons emitted from the incoming quark lines. At high Q 2 scales these corrections should be resummed. Fortunately the QCD factorisation theorem also applies to QED, and so the photon-induced logarithms can be absorbed into the PDFs, just as the α s logQ 2 terms are summed by DGLAP evolution. As a consequence the normal DGLAP equations are slightly modified and a photon PDF of the proton, γ p (x, Q 2 ), is generated. Thus, for example, (at leading order in both α s and α) we have an extra equation for the evolution of the photon PDF ∂γ (x, Q 2 ) Similarly, the photon PDF, γ , contributes to the evolution ∂q i /∂logQ 2 via the P qγ splitting.

Existing determinations of γ p compared to this work
Indeed, with the present level of precision, it has become topical to treat the photon as one of the point-like partons inside the nucleon and to account for this QED effect explicitly in the global parton analysis. This approach was first followed 10 years ago by MRST (2004) [1]. Recently it has been used by the NNPDF [2] and CTEQ [3] groups. The central issue is the choice of input distributions for the photon PDFs of the proton and neutron. In the original MRST study it was assumed that the starting distributions are given by one-photon emission off valence (constituent) quarks in the leading logarithm approximation. For example, for the photon PDF of the proton the starting distribution was taken to be 1 where u 0 and d 0 are the valence-like distributions of the proton, and where the current quark masses were used. The most direct measurement of the photon PDF at the time of the MRST analysis [4] It appeared to be wide-angle scattering of the photon by an electron beam via the process ep → eγ X , where the final state electron and photon are produced with equal and opposite large transverse momentum.
The subprocess is QED Compton scattering, 2 eγ → eγ , for which the contribution to the cross section is where μ is the factorisation scale. MRST [1] predicted a cross section in agreement with the only measurement of this process available at that time [5].
The NNPDF [2] and CTEQ [3] groups use a different approach to MRST. They parametrise the input photon PDFs, γ (x, Q 2 0 ), and attempt to determine the parameters from the global data, along with the quark and gluon PDFs. 3 Unfortunately the present data are not of sufficient accuracy to provide a reasonable determination of the photon input.
The NNPDF collaboration [2] used freely parametrised (without bias) starting distributions, including the photon PDFs, and constrain the photon PDFs mainly from the Drell-Yan (low-mass, on-shell W and Z production and high-mass) LHC data. There is expected to be the most sensitivity to the low-mass Drell-Yan data [9]. However, the uncertainties observed in the resulting photon PDFs are huge, especially at low x.
The preliminary CTEQ analysis [3] proceeds differently. CTEQ keep a similar theoretical form of the distributions γ (x, Q 2 0 ) to that proposed by MRST, but with an arbitrary normalisation parameter, which is expressed as the momentum fraction, p 0 (γ ), carried by the input photon. They find that the constraint coming from the energy-momentum sum rule is weak (allowing p 0 (γ ) to range up to 5 %), while to fit the updated ZEUS data for ep → eγ X [10] requires p 0 (γ ) ∼ 0.1-0.2 %, using the valence quark induced input (3) and allowing for the extra normalisation parameter.
Unlike the above analyses, here we emphasise that the major part of the input, γ p (x, Q 2 0 ), especially at low x, comes from the coherent emission of the photon from the 'elastic' proton, which can be calculated theoretically with good accuracy. The process is sketched in the left diagram of Fig. 1.
Indeed, by definition, inclusive PDFs include the contribution from all possible final states. With some probability, w c , after photon emission, the coherence between the partons in initial nucleon wave function may be not destroyed and the final state will be just the same (intact) nucleon, as shown in the left diagram of Fig. 1. This probability is given by the square of the form factor F 1 , that is, w c = F 2 1 (t). At large photon virtuality, |t|, it falls down rapidly and may be considered as a negligible power correction to the DGLAP equation. However, at the beginning of DGLAP evolution (at a low |t|) this coherent contribution is important, and The coherent and incoherent contributions to the photon PDF, γ p (x, Q 2 ), corresponding, respectively, to photon emission directly from the proton and from a quark it should be considered separately; see (6) below. In order to avoid double counting, the remaining (incoherent) contribution must be multiplied by the 'remaining' probability, ; see (9) below. Note that the previous analyses are based only on incoherent emission from individual quarks within the proton. The uncertainty on our determination γ p (x, comes essentially only from the relatively small contribution γ p incoh which, from a hadron viewpoint, actually corresponds to the QED excitations of the proton into higher-mass states. However, here, in Sect. 3 we adopt the quark viewpoint and calculate this contribution as the incoherent emission of photons from quarks within the proton. This contribution turns out to be relatively small. Therefore, since γ p coh may be calculated with good accuracy, it means that the uncertainty in the theoretically determined photon input distributions is expected to be small; essentially coming from the uncertainty in the 'extrapolation' of valence quarks needed to estimate the contribution to γ p incoh coming from the region below the starting scale Q 0 . We will quantify this uncertainty below.
We summarise the discussion of this section in Table 1.

Improved input distributions for the photon PDFs
Here we will follow the MRST approach, but we will use much improved starting distributions for the photon PDFs of the proton and neutron. Indeed, we have where N = p, n. As discussed above, and sketched in Fig. 1, the contribution γ p coh is caused by coherent photon emission from the proton that remains intact, whereas γ incoh is due to non-coherent emission from individual quarks. The coherent emission from the proton is given by [11] where q t is the transverse momentum of the emitted photon and F 1 is the electromagnetic proton form factor corresponding to γ μ at the vertex. For the neutron we have For the non-coherent emission from individual quark lines we use an improved form of (3) where accounts for the fact that the lowest possible proton excitation is the -isobar. The factor (1 − F 2 1 ) in (9) is the probability to have no intact proton in the final state. We have to exclude an intact proton as its contribution is calculated separately in (6).
In (9), m q = m d when convoluted with d 0 , and m q = m u when convoluted with u 0 . 4 In this contribution we use the current quark masses. Here the quark distribution u 0 = u valence + 2u sea is frozen for Q < Q 0 at its value at Q 0 . The same is true for the other quarks-d, s. A similar expression holds for γ n incoh , with 4/9 ↔ 1/9 and F p 1 → F n 1 . In this way we get an upper limit for the non-coherent contribution to the photon input distributions. The other extreme is to take for u 0 and d 0 just the non-relativistic quark model expectation with 4 To be precise, we replace the integral dt/(t − m 2 q ) by with t min given by (10).
for the proton and to use constituent quark masses m q = 300-350 MeV.
The optimum estimate of the non-coherent contribution to the photon PDF input is probably to take a, physicsmotivated, linear interpolation between the two limits. That is, to use in (9) with m q = m current + m eff (t), where the 'effective' constituent quark mass is parametrised by a simplified formula of the form with m(0) = 345 MeV and slope b = 1.4 GeV −1 (see, for example, Fig. 4 in [12], where the light quark in the instanton vacuum was studied).
In general, one may also account for the -isobar excitation. In the latter case, we have to add to (6) γ coh , which is also of the form of (6), but with 5 For the contribution Also when including the contribution we have to replace in (9) where here |t| is given by (15). In addition, it is possible to include a coherent contribution caused by the anomalous magnetic moment of the proton, described by the proton form factor F 2 . These non-logarithmic corrections will reduce the remaining incoherent contribution and therefore decrease the final uncertainty in the input γ p (x, Q 2 0 ). However, since they do not change the result noticeably, we do not consider these possibilities here.
Moreover, recall that in the case of photon emission arising from the γ μ vertex of the form factor F 1 , we get [ignoring the t min factor x 2 m 2 p in (6)] the same logarithmic ( dq 2 t /q 2 t ) distribution as that in conventional DGLAP evolution, which justifies the strong q t ordering: q 2 t Q 2 0 . 6 On the contrary, the vertex of the -isobar excitation (or the proton spinflip contribution driven by the form factor F 2 ) contains an additional factor q t in the numerator. Therefore here the typical value of photon virtuality will be |t| ∼ Q 2 0 . Thus one cannot use, at the next step of the evolution, the conventional splitting function P qγ calculated under the assumption that the incoming parton (photon) has a virtuality much smaller than that at the next step of the evolution (and thus can be considered as on-mass-shell). To be safe here, we do not consider the contributions of the baryon resonances excitations and of the proton spin-flip, as separate parts of our photon input, but instead we include them in the whole incoherent input described in terms of the radiation from quarks.
Our approach is close to that proposed in [13] where the coherent contribution from an intact proton was calculated using the Weizsäcker-Williams approximation accounting for both form factors, F 1 and F 2 , and assuming zero incoherent input at the starting scale Q 2 = Q 2 0 . In comparison with [13], we exclude the spin-flip contribution caused by F 2 (since it does not have the conventional DGLAP logarithmic structure), but we do not ignore the initial non-coherent part at Q 2 = Q 2 0 formed by radiation from individual quarks. The difference is especially clear for the photon distribution inside a neutron.

Results for the photon PDF
In Figs. 2 and 3 we show, respectively, the input distributions of the photon PDF for the proton and neutron at Q 2 0 = 1 GeV 2 , as calculated in the previous section, together with the photon PDF determined from the evolution up to Q 2 = 200 GeV 2 using (1), but with NLO partons. We see that the major part of the input photon PDF of the proton is generated by coherent emission of the photon from an 'elastic' proton, while the non-coherent contribution only enlarges this value by about 20 %. Actually, the figure shows three curves for the inclusion of γ p incoh , corresponding to the two limits of the input distribution together with their 'interpolated average', shown by the continuous curve. We also note that for a low x ∼ 10 −2 -10 −3 the evolution to Q 2 = 200 GeV 2 increases the photon density by about a factor of two; whereas for x = 10 −4 the increase is about a factor of three. Considering the two terms involving P γ γ and P γ q of (1), which contribute to the evolution of the photon PDF, we note that the first term, 6 Since the q t dependence of the photon flux produced by the coherent emission (via the γ μ vertex of the F 1 form factor) has the same logarithmic form as that for the LO DGLAP partons, we have the possibility to treat these photons as a new sort of point-like massless partons, and we hence include them in the conventional DGLAP evolution. The photon PDF of the neutron at input (Q 2 0 = 1 GeV 2 ) and after evolution up to Q 2 = 200 GeV 2 . At each Q 2 value, the three curves correspond to the upper and lower estimates of γ n incoh , together with continuous ('interpolating') curve obtained from (12) P γ γ , decreases γ p only slightly (<1 % in the evolution up to Q 2 = 200 GeV 2 ), whereas the growth comes from the P γ q term-the photons emitted by quarks. We have the same Here we evolve using MSTW [14] NLO quarks. 7 At first sight we might expect the contribution generated by gluons, via the gluon-photon splitting, to be important, due to the large gluon PDF, especially at low x. However, the value of the NLO splitting function P (1) γ g (z) is rather small. 8 Moreover, P (1) γ g (z) is negative at large z. As a result, the inclusion of gluons enlarges γ p (x, Q 2 ) by <2 %, and we ignore this effect.
Recall that for the photon PDF of the neutron we have γ n coh (x, Q 2 0 ) = 0, and so the input is given entirely by γ n incoh (x, Q 2 0 ); see Fig. 3. However, the increase in γ n in the evolution up to Q 2 = 200 GeV 2 [which is driven by the final term in (1)] is comparable to that for γ p . In Fig. 4 we compare γ p and γ n with the other PDFs at Q 2 = 20 GeV 2 .
In Fig. 5 we show how γ p evolves as Q 2 increases from input Q 2 0 = 1 to Q 2 = 20, 200 and 10 4 GeV 2 , together with the behaviour predicted 9 by MRST (2004) [1] input. Figure  5 indicates that at very large Q 2 the behaviour of γ p (x, Q 2 ) becomes stable and relatively insensitive to the form taken 7 We find that the use of the updated CPdeut parton set of MMSTWW [15] makes a negligible difference. 8 We extract the splitting function as the term proportional to C F N f from the known P (1) gg (z) splitting [16]. 9 Actually we use the MRST formulation, but with NLO MSTW partons [14], which make little change to the behaviour of γ p . for the input. The reason is that the evolution of (1) is driven by P γ q ⊗ q i and that the quark PDFs are well determined.
The comparison of the predictions by MRST [1] and of this work, may, at first sight, appear surprising. MRST is purely based on the incoherent contribution, γ p incoh , whereas here the prediction comes dominantly from the coherent contribution (dashed-dotted curve) with a small addition from γ p incoh . The explanation is as follows. The incoherent contribution determined by MRST should be suppressed by t min , (10), and by [1− F 2 1 (t)] of (9). In the present work, the coherent emission from the proton is added. The above two effects (that is, the suppression of the incoherent contribution and the inclusion of the coherent emission) partly compensate each other. However, indeed at large x, where |t min | is large, MRST goes above the present input, while at low x, where |t min | (xm N ) 2 is even less than the current quark mass, our input exceeds the MRST curve.

Comparison with ep → eγ X data
To probe the photon PDF experimentally we should consider a hard subprocess where the photon distribution dominates. A good example is the inclusive electroproduction of isolated photons observed in the direction of the incoming electron. This reaction is mediated by the eγ → eγ hard subprocess and its cross section is described by (4), and sketched in Fig. 6a. It is well known that this Compton scattering process is sharply peaked in the backward direction. Therefore the outgoing photon should be observed at high negative rapidity, η γ , at angles close to the electron beam. Indeed it was already shown [17,18] that this component (called LL) already dominates for 10 η γ < −1.
Unfortunately, the present experimental data are quite limited and the lowest rapidity bin is −0.7 < η γ < −0.3. We compare our theoretical prediction for the LL component with the data in Tables 2 and 3. As can be seen from the comparison in the first bin, the predicted cross section is close to the measured value. The result depends on the choice of factorisation scale. We present results for μ = E γ T and E γ T /2 to indicate the scale dependence.
At large η γ , the contribution of the Compton-induced process decreases rapidly. In this domain, inclusive isolated photons are mainly produced by quarks; see Fig. 6b.
Note that in our theoretical calculation of ep → eγ X we have accounted for the angular, the E T and the other experimental cuts imposed by the ZEUS collaboration [10], but we have no possibility to include the photon isolation criteria. Therefore the observed cross section corresponding to the LL process should be lower than our prediction. The photon PDF at input Q 2 0 = 1 GeV 2 and after evolution up to Q 2 = 20, 200, and 10 4 GeV 2 . The notation of the curves are as in Fig. 2. We have included the prediction of MRST (2004) [1] for comparison  Table 2 The second column gives the values of the ep → eγ X cross section measured by the ZEUS collaboration [10] in different rapidity, η γ , intervals. The final two columns show the contribution to the cross section arising from the LL process of Fig. 6a for two different choices of the factorisation scale μ  Table 3 The second column gives the values of the ep → eγ X cross section measured by the ZEUS collaboration [10] in different E γ T intervals. The final two columns show the contribution to the cross section arising from the LL process of Fig. 6a  In Table 3 we compare our prediction of the LL contribution with E T dependence of the measured cross section. However, now the data were collected over a large rapidity interval: −0.7 < η γ < 0.9. Here the quark contribution is important, and the LL subprocess describes only about half of the cross section.

Conclusions
We have demonstrated that the major part of the photon input PDF of the proton (caused by the coherent emission of the photon that does not destroy the proton) can be calculated with good accuracy. This strongly reduces the possible uncertainties in the QED part of an extended global parton analysis which includes the photon as a point-like parton. In this way, we evaluate the expected photon PDFs by DGLAP evolution with LO QED splittings and NLO MSTW quarks. Note that the further step of including the photon-to-quark splitting will introduce a small (suppressed by the α QED /2π factor) violation of isospin symmetry in the 'singlet' PDF, in particular u d = d n . The resulting photon distributions agree with data for the electroproduction of isolated photons, ep → eγ X , at negative rapidities where the cross section is dominated by the eγ → eγ hard subprocess.