Measurement of isolated photon production in deep inelastic ep scattering

Isolated photon production in deep inelastic ep scattering has been measured with the ZEUS detector at HERA using an integrated luminosity of 320 pb. Measurements were made in the isolated-photon transverse-energy and pseudorapidity ranges 4 < E T < 15 GeV and −0.7 < η < 0.9 for exchanged photon virtualities, Q, in the range 10 < Q < 350 GeV and for invariant masses of the hadronic system WX > 5 GeV. Differential cross sections are presented for inclusive isolated photon production as functions of Q, x, E T and η . Leadinglogarithm parton-shower Monte Carlo simulations and perturbative QCD predictions give a reasonable description of the data over most of the kinematic range. The ZEUS Collaboration S. Chekanov, M. Derrick, S. Magill, B. Musgrave, D. Nicholass, J. Repond, R. Yoshida Argonne National Laboratory, Argonne, Illinois 60439-4815, USA n M.C.K. Mattingly Andrews University, Berrien Springs, Michigan 49104-0380, USA P. Antonioli, G. Bari, L. Bellagamba, D. Boscherini, A. Bruni, G. Bruni, F. Cindolo, M. Corradi, G. Iacobucci, A. Margotti, R. Nania, A. Polini INFN Bologna, Bologna, Italy e S. Antonelli, M. Basile, M. Bindi, L. Cifarelli, A. Contin, S. De Pasquale, G. Sartorelli, A. Zichichi University and INFN Bologna, Bologna, Italy e D. Bartsch, I. Brock, H. Hartmann, E. Hilger, H.-P. Jakob, M. Jüngst, A.E. Nuncio-Quiroz, E. Paul, U. Samson, V. Schönberg, R. Shehzadi, M. Wlasenko Physikalisches Institut der Universität Bonn, Bonn, Germany b J.D. Morris H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom m M. Kaur, P. Kaur, I. Singh Panjab University, Department of Physics, Chandigarh, India M. Capua, S. Fazio, A. Mastroberardino, M. Schioppa, G. Susinno, E. Tassi Calabria University, Physics Department and INFN, Cosenza, Italy e J.Y. Kim Chonnam National University, Kwangju, South Korea Z.A. Ibrahim, F. Mohamad Idris, B. Kamaluddin, W.A.T. Wan Abdullah Jabatan Fizik, Universiti Malaya, 50603 Kuala Lumpur, Malaysia r Y. Ning, Z. Ren, F. Sciulli Nevis Laboratories, Columbia University, Irvington on Hudson, New York 10027, USA o J. Chwastowski, A. Eskreys, J. Figiel, A. Galas, K. Olkiewicz, B. Pawlik, P. Stopa, L. Zawiejski The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Poland i L. Adamczyk, T. Bo ld, I. Grabowska-Bo ld, D. Kisielewska, J. Lukasik, M. Przybycień, L. Suszycki Faculty of Physics and Applied Computer Science, AGH-University of Science and Technology, Cracow, Poland p


DIS Photon emission Quark radiation Isolated photon
Isolated photon production in deep inelastic ep scattering has been measured with the ZEUS detector at HERA using an integrated luminosity of 320 pb −1 . Measurements were made in the isolated-photon transverse-energy and pseudorapidity ranges 4 < E γ T < 15 GeV and −0.7 < η γ < 0.9 for exchanged photon virtualities, Q 2 , in the range 10 < Q 2 < 350 GeV 2 and for invariant masses of the hadronic system W X > 5 GeV. Differential cross sections are presented for inclusive isolated photon production as functions of Q 2 , x, E γ T and η γ . Leading-logarithm parton-shower Monte Carlo simulations and perturbative QCD predictions give a reasonable description of the data over most of the kinematic range.

Introduction
In the study of high-energy collisions involving hadrons, events in which an isolated high-energy photon is observed provide a direct probe of the underlying partonic process, since the emission of these photons is unaffected by parton hadronisation. * Corresponding author.
E-mail address: tobias.haas@desy.de (T. Haas). 1 Also affiliated with University College London, United Kingdom. 2 Supported by the US Department of Energy. 3 Supported by the Italian National Institute for Nuclear Physics (INFN). 4 Now at University of Salerno, Italy. 5 Supported by the German Federal Ministry for Education and Research (BMBF), under contract Nos. 05 HZ6PDA, 05 HZ6GUA, 05 HZ6VFA and 05 HZ4KHA. 6 Now at Queen Mary University of London, United Kingdom. 7 Supported by the Science and Technology Facilities Council, UK. 8 Also working at Max Planck Institute, Munich, Germany. 9 Also Senior Alexander von Humboldt Research Fellow at Hamburg University, Institute of Experimental Physics, Hamburg, Germany. 10 Supported by Chonnam National University, South Korea, in 2009. 11 Supported by an FRGS grant from the Malaysian government. 12 Supported by the US National Science Foundation. Any opinion, findings and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation. 13   Isolated high-energy photon production has been studied in a number of fixed-target and hadron-collider experiments [1]. Previous ZEUS and H1 publications have also reported the production of isolated photons in photoproduction [2][3][4][5][6], in which the exchanged photon is quasi-real ( Q 2 ≈ 0), and deep inelastic scattering (DIS) [7,8], in which Q 2 ≈ GeV 2 .
Isolated photons are produced in DIS at lowest order in QCD as shown in Fig. 1. Photons produced by radiation from an incoming or outgoing quark are called "prompt"; an additional class of high-energy photons comprises those radiated from the incoming or outgoing lepton. In this Letter, results are presented from a new inclusive measurement of isolated photon production in neutral current DIS. The data provide a test of perturbative QCD in a kinematic region with two hard scales: Q 2 , the exchanged photon virtuality, and E γ T , the transverse energy of the emitted photon. Compared to the previous ZEUS publication [7], the kinematic reach extends to lower values of Q 2 and to higher values of E γ T . The statistical precision is also improved. Leading-logarithm parton-shower Monte Carlo (MC) and perturbative QCD predictions are compared to the measurements. The cross sections for isolated photon production in DIS have been

Experimental set-up
The measurements are based on a data sample corresponding to an integrated luminosity of 320 ± 131 ± 3 pb −1 of e + p data and 189 ± 5 pb −1 of e − p data 57 with centre-of-mass energy √ s = 318 GeV.
A detailed description of the ZEUS detector can be found elsewhere [13]. Charged particles were tracked in the central tracking detector (CTD) [14] and a silicon micro vertex detector (MVD) [15] which operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The high-resolution uranium-scintillator calorimeter (CAL) [16] consisted of three parts: the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. The BCAL covers the pseudorapidity range −0.74 to 1.01 as seen from the nominal interaction point. The FCAL and RCAL extend the range to −3.5 to 4.0. The smallest subdivision of the CAL was called a cell.
The barrel electromagnetic calorimeter (BEMC) cells had a pointing geometry aimed at the nominal interaction point, with a cross section approximately 5 × 20 cm 2 , with the finer granularity in the Z -direction. 58 This fine granularity allows the use of shower-shape distributions to distinguish isolated photons from the products of neutral meson decays such as π 0 → γ γ .
A three-level trigger system was used to select events online [17] by requiring well isolated electromagnetic deposits in the CAL.
The luminosity was measured using the Bethe-Heitler reaction ep → eγ p by a luminosity detector which consisted of two independent systems: a lead-scintillator calorimeter [18] and a magnetic spectrometer [19].

Event selection and reconstruction
Events were selected offline by requiring a scattered-electron candidate, identified using a neural network [20]. The candidates were required to have a polar angle in the range 139.8 • < θ e < 171.9 • in order to ensure that they were well measured in the RCAL. The impact point ( X, Y ) of the candidate on the surface of the RCAL was required to lie outside the region (±15 cm, ±15 cm) centred on (0, 0) to ensure well understood acceptance. The energy of the candidate, E e , was required to be larger than 10 GeV. The kinematic quantities Q 2 and x were reconstructed from the scattered electron by means of the relationships Q 2 = −(k − k ) 2 and is the four-momentum of the incoming (outgoing) lepton and P is the four-momentum of the incoming proton. The kinematic region 10 < Q 2 < 350 GeV 2 was selected.
To reduce backgrounds from non-ep collisions, events were required to have a reconstructed vertex position, Z vtx , within the range |Z vtx | < 40 cm and to have 35 is the energy of the ith CAL cell, θ i is its polar angle and the sum runs over all cells [21]. At least one reconstructed track, well separated from the electron, was required, ensuring some hadronic activity which suppressed deeply virtual Compton scattering (DVCS) [22] to a negligible level.
Photon candidates were identified as CAL energy-flow objects (EFOs) [23] for which at least 90% of the reconstructed energy was measured in the BEMC. EFOs with wider electromagnetic showers than are typical of a single photon were accepted to allow evaluation of backgrounds. The reconstructed transverse energy of the EFO, E γ T , was required to lie within the range 4 < E γ T < 15 GeV and the pseudorapidity, η γ , had to satisfy −0.7 < η γ < 0.9. The upper limit on the reconstructed transverse energy was selected to en-57 Hereafter 'electron' refers to both electrons and positrons unless otherwise specified. 58 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the proton beam direction, referred to as the "forward direction", and the X axis pointing towards the centre of HERA. The coordinate origin is at the nominal interaction point. sure that the shower shapes from background and signal remained distinguishable.
To reduce the background from photons and neutral mesons within jets, the EFO was required to be isolated from reconstructed tracks and hadronic activity. Isolation from tracks was initially is the distance to the nearest reconstructed track with momentum greater than 250 MeV in the η-φ plane, where φ is the azimuthal angle. Jet reconstruction was performed on all EFOs in the event, including the electron and photon candidates, using the k T cluster algorithm [24] in the longitudinally invariant inclusive mode [25] with R parameter set to 1.0. Further isolation was imposed by requiring that the photon-candidate EFO possessed at least 90% of the total energy of the jet of which it formed a part.
Each event was required to contain both an electron and a photon candidate. The invariant mass of the hadronic system, W X , is then defined by W 2 where p γ is the four-vector of the outgoing photon. A total of 15 699 events were selected; at this stage the sample was dominated by background events. The largest source of background was neutral current (NC) DIS events where a genuine electron candidate was found in the RCAL and neutral mesons, such as π 0 and η, decaying to photons, produced a photon-candidate EFO in the BEMC.

Theory
Two theoretical predictions are compared to the measurements presented in this Letter. In the approach of GGP [10], the contributions to the scattering cross section for ep → eγ X are calculated at order α 3 in the electromagnetic coupling. One of these contributions comes from the radiation of a photon from the quark line (called QQ photons; Fig. 1a, b) and a second from the radiation from the lepton line (called LL photons; Fig. 1c, d). In addition to QQ and LL photons, the interference term between photon emission from the lepton and quark lines, called LQ photons by GGP, is evaluated. For the kinematic region considered here, where the outgoing photon is well separated from both outgoing electron and quark, the interference term gives only a 3% effect on the cross section. This effect is further reduced to ≈ 1% when e + p and e − p data are combined as the LQ term changes sign when e − is replaced by e + . The QQ contribution includes both wide-angle photon emission and the leading q → qγ fragmentation term. GGP have chosen to use CTEQ6L leading-order parton distribution functions [26]. The factorisation scales used are Q for QQ events and max(Q , μ F ,min ) for LL events where μ F ,min = 1 GeV. Parton-tohadron corrections were not made, in view of technical issues in relating 2 → 2 and 2 → 3 topologies, following the advice of the GGP authors. We note that others have taken a different view [8]. A naïve study indicated the likely effect to be a reduction after hadronisation in predicted inclusive cross-sections of order 15%.
In the approach of MRST [12,27], a partonic photon component of the proton, γ p , is introduced as a consequence of including QED corrections in the parton distribution functions. This leads to ep interactions taking place via QED Compton scattering, γ p e → γ e.
A measurement of the isolated high-energy photon production cross section therefore provides a constraint on the photon density in the proton. The model includes the collinearly divergent LL contribution, which is enhanced relative to that of GGP by the DGLAP resummation due to the inclusion of QED Compton scattering. The QQ component is not included in the MRST model, in which the transverse momentum of the scattered electron is expected to balance approximately that of the isolated photon. In the analysis presented here, such a constraint was not imposed. The theoretical uncertainties in the models have been estimated by varying the factorisation scales by a factor two.
Since the MRST cross sections include the LL contribution of GGP to a good approximation, but exclude the QQ, an improved prediction can be constructed by summing the MRST cross section and the QQ cross section from GGP [27,28]. The theory uncertainties are of the same order as those of the individual QQ and LL components.

Monte Carlo event simulation
The MC program Pythia 6.416 [29] was used to simulate prompt-photon emission for the study of the event-reconstruction efficiency. In Pythia, this process is simulated as a DIS process with additional photon radiation from the quark line to account for QQ photons. Radiation from the lepton is not simulated in this Pythia sample.
The LL photons radiated at large angles from the incoming or outgoing electron were simulated using the generator Djangoh 6 [30], an interface to the MC program Heracles 4.6.6 [31]; higherorder QCD effects were included using the colour dipole model of Ariadne 4.12 [32]. Hadronisation of the partonic final state was performed by Jetset 7.4 [33]. The small LQ contribution was neglected.
The NC DIS background was simulated using Djangoh 6, within the same framework as the LL events. This provided a realistic spectrum of mesons and overlapping clusters with well modelled kinematic distributions and hence was preferred to single-particle MC samples for backgrounds, such as were used in the previous ZEUS publication [7].
The MC samples described above contained only events in which W X was larger than 5 GeV. Isolated photons can also be produced at values of W X less than 5 GeV in 'elastic' and 'quasielastic' processes (ep → epγ ) such as DVCS and Bethe-Heitler photon production. Such events were simulated using the GenDVCS [34] and Grape-Compton [35] generators. The contribution of these elastic processes was negligible after the selections described in Section 3.
The generated MC events were passed through the ZEUS detector and trigger simulation programs based on Geant 3.21 [36]. They were reconstructed and analysed by the same programs as the data. In addition to the full-event simulations, MC samples of single particles (photons and neutral mesons) were generated and used to study the MC description of electromagnetic showering in the BEMC.

Extraction of the photon signal
The event sample selected according to the criteria in Section 3 was dominated by background; thus the photon signal was extracted statistically following the approach used in previous ZEUS analyses [2][3][4]7].
The photon signal was extracted from the background using BEMC energy-cluster shapes. Two shape variables were considered: of the centre of the ith cell, Z cluster is the centroid of the EFO cluster, w cell is the width of the cell in the Z direction, E i is the energy recorded in the cell and the sum runs over all BEMC cells in the EFO; • the ratio f max of the highest energy deposited in any one BEMC cell in the EFO to the total EFO BEMC energy.
The distributions of δ Z and f max (after the requirement δ Z < 0.8) in the data and the MC are shown in Fig. 2. The MC LL and QQ distributions have been corrected in each two-dimensional (η, E T ) bin using factors derived from the difference between simulated and real DIS electron data. The δ Z distribution exhibits a double-peaked structure with the first peak at ≈ 0.1, associated with the signal, and a second peak at ≈ 0.5, dominated by the π 0 → γ γ background. The f max distribution shows a single peak at ≈ 0.9 corresponding to the photon signal, and has a shoulder extending down to ≈ 0.5, which is dominated by the hadronic background.
The number of isolated-photon events contributing to Fig. 2 and in each cross-section bin was determined by a χ 2 fit to the δ Z distribution in the range 0 < δ Z < 0.8 using the LL and QQ signal and background MC distributions as described in Section 5. By treating the LL and QQ photons separately, one automatically takes account of their differing hadronic activity (resulting in significantly different acceptances) and their differing (η, E T ) distributions (resulting in different bin migrations due to finite measuring precision).
In performing the fit, the LL contribution was kept constant at its MC-predicted value and the other components were varied. Of the 15 699 events selected, 4164 ± 168 correspond to the extracted signal (LL and QQ). The scale factor resulting from the global fit for the QQ photons in Fig. 2 was 1.6; this factor was used for all the plots comparing MC to data. The fitted global scale factor for the hadronic background was 1.0. The signal fraction in the crosssection bins varied from 21% to 62%. In all cross-section bins, the χ 2 /n.d.f. of the fits was 2.1 or smaller.
For a given observable Y , the production cross section was determined using: where N(γ QQ ) is the number of QQ photons extracted from the fit, Y is the bin width, L is the total integrated luminosity, σ MC LL is the predicted cross section for LL photons from Djangoh, and A QQ is the acceptance correction for QQ photons. The value of A QQ was calculated using Monte Carlo from the ratio of the number of events generated to those reconstructed in a given bin. It varied between 1.2 and 1.7 from bin to bin.
The fits employed in this analysis were performed using δ Z because of the larger difference in shape between signal and background for this quantity. Fits in terms of the f max distributions were performed as a cross-check and gave similar results. As a further cross-check, an algorithm from the previous ZEUS publication [7], which selects wider electromagnetic clusters as photon candidates, was used. This proved to be more sensitive to the modelling of calorimeter backgrounds. In every case where a satisfactory fit was obtained, good agreement with the principal method was found. The corrections to the MC photon-signal energy-cluster shapes gave changes to the results within the statistical uncertainties and were not further considered [37].

Systematic uncertainties
The following sources of systematic uncertainty were investigated [37]: • the energy scale of the electromagnetic calorimeter (EMC) was varied by its known scale uncertainty of ±2% causing variations in the measured cross sections of typically less than ±2%; • the dependence on the modelling of the hadronic background by Ariadne was investigated by varying the upper limit for the δ Z fit in the range 0.6-1.0, giving variations that were typically ±5% but up to +12% and −14% in the most forward η γ and highest-x bins respectively.
The following sources of systematic uncertainty were also investigated and found to be negligible compared to the statistical uncertainty [37]: • variation of the EMC energy-fraction cut for the photon candidate EFO by ±5%; • variation of the Z vtx cut by ±5 cm; • variation of the upper and lower cuts on δ by ±3 GeV; • variation of the R cut used for track isolation by ±0.1; • variation of the track-momentum cut used in calculating track isolation by ±100 MeV; • variation of the LL-signal component by ±5%.
All the uncertainties listed above were added in quadrature to give separate positive and negative systematic uncertainties in each bin. The uncertainty of 2.6% on the luminosity measurement was not included in the differential cross sections but included in the integrated cross sections.

Results
The cross section for inclusive isolated photon production, ep → eγ X , was measured in the kinematic region defined by: 10 < Q 2 < 350 GeV 2 , W X > 5 GeV, E e > 10 GeV, 139.8 • < θ e < 171.8 • , −0.7 < η γ < 0.9 and 4 < E γ T < 15 GeV, with isolation such that at least 90% of the energy of the jet containing the photon belongs to the photon, where jets were formed according to the k T algorithm with R parameter set 1.0. The measured integrated cross section is      Fig. 3 and given in Tables 1-4. It can be seen that the cross section decreases with increasing E γ T , η γ , Q 2 and x. The predictions for the sum of the expected LL contribution from Djangoh and a factor of approximately 1.6 times the expected QQ contribution from Pythia agree well with the measurements, except for some differences at the lowest Q 2 (and correspondingly lowest x). The theoretical predictions described in Section 4 are compared to the measurements in Fig. 4. The predictions from GGP describe the shape of the E γ T and η γ distributions well, but their central value typically lies 20% below the measured cross sections. The calculations fail to reproduce the shape in Q 2 ; a similar observation was made by H1 [8]. As with the MC comparison, the measured cross section is larger than the theoretical prediction; this is also reflected in an excess of data over theory at low x.
The MRST predictions mostly fall below the measured differential cross sections. However, they lie close to the measurements at large values of Q 2 and x, for backward η γ and for high values of E γ T , where the LL cross section is expected to be a substantial fraction of the total. Also included in Fig. 4 is the sum of MRST and QQ of GGP; it gives an improved description of the data over much of the range of the kinematic variables.   5 shows the measured dσ /dη γ compared to previous measurements from ZEUS [7] and H1 [8] for the restricted range Q 2 > 35 GeV 2 and 5 < E γ T < 10 GeV. The results are consistent but the uncertainty in the present measurement is smaller. The symbols are mutually displaced for clarity.

Conclusions
Inclusive isolated photon production has been measured in deep inelastic scattering using the ZEUS detector at HERA using an integrated luminosity of 320 pb −1 . Differential cross sections as functions of several kinematic variables are presented for 10 < Q 2 < 350 GeV 2 and W X > 5 GeV in the pseudorapidity range −0.7 < η γ < 0.9 for photon transverse energies in the range 4 < E γ T < 15 GeV. The order α 3 predictions of Gehrmann-de Ridder et al. reproduce the shapes of the experimental results as functions of transverse energy and pseudorapidity, but are lower than the measurements at low Q 2 and low x. The predictions of Martin et al. mostly fall below the measured cross sections but are close in the kinematic regions where lepton emission is expected to be dominant. An improved description of the data is obtained by appropriately combining the two predictions, suggesting a need for further calculations to exploit the full potential of the measurements.