Measurement of Di(cid:11)ractive Scattering of Photons with Large Momentum Transfer at HERA

The (cid:12)rst measurement of di(cid:11)ractive scattering of quasi-real photons with large momentum transfer (cid:13)p ! (cid:13)Y , where Y is the proton dissociative system, is made using the H1 detector at HERA. Single di(cid:11)erential cross sections are measured as a function of W , the incident photon-proton centre of mass energy, and t , the square of the four-momentum transferred at the proton vertex. The W dependence is well described by a perturbative QCD model using a leading logarithmic approximation of the BFKL evolution, whereas the measured j t j dependence is harder than predicted.

The study at the ep collider HERA of exclusive diffractive processes in the presence of a hard scale provides insight into the parton dynamics of the diffractive exchange. The four-momentum squared transferred at the proton vertex, t, provides a relevant scale to investigate the application of perturbative Quantum Chromodynamics (pQCD) for |t| Λ 2 QCD [2]. The diffractive photonproton scattering, γp → γY at large |t| gives an experimentally clean and almost fully perturbatively calculable process to study the parton dynamics. Its cross section measurement [3] is performed at HERA by studying the reaction e + p → e + γY in the photoproduction regime (initial photon virtualities Q 2 < 0.01 GeV 2 ) and at 4 < |t| < 36 GeV 2 , with a large rapidity gap between the final state photon and the proton dissociative system Y (as illustrated in Fig. 1) and the γp centre of mass energy in the range 175 < W < 247 GeV.
Diffractive photon scattering can be modelled in the proton rest frame by the fluctuation of the incoming photon into a qq pair, which is then involved in a hard interaction with the proton via the exchange of a colour singlet state. In the leading logarithmic approximation (LLA), the colour singlet exchange is modelled by the effective exchange of a gluon ladder (Fig. 2), described at sufficiently low values of Bjorken x by the BFKL [4] approach. In the LLA BFKL model the gluon ladder couples to a single parton within the proton. Due to the quasi-real nature of the incoming photon, the transverse momentum of the final state photon, P γ T , is entirely transferred by the gluon ladder to the parton in the proton. The separation in rapidity between the parton scattered by the gluon ladder and the final state photon is given by ∆η log(ŝ/(P γ T ) 2 ), whereŝ is the invariant mass of the system formed by the incoming photon and the struck parton. In addition to the usual DIS kinematic variables, the longitudinal momentum fraction of the diffractive exchange with respect to the proton, x I P and the elasticity of the γp interaction, y I P , are defined as (Following the notation in Fig. 1) 2 Experimental procedure The data, corresponding to an integrated luminosity of 46.2 pb −1 , were collected with the H1 detector during the 1999 − 2000 running period, when positrons of energy E e = 27.6 GeV collided with protons of energy E p = 920 GeV. The detailed description of the H1 detector can be found in [5]. The quantity W is calculated from the measured energy of the scattered positron, E e , using the relation W 2 ys s(1 − E e /E e ) with the relative resolution 4%, where s = (k + p) 2 is the ep centre of mass energy squared. The scattered positron is detected in a small calorimeter 33 m downstream the main detector. The photon candidate is detected as a cluster in the electromagnetic section of the calorimeter, with no associated track and is furthermore required to have an energy E γ > 8 GeV and a transverse momentum P γ T > 2 GeV. The kinematic variable |t| is reconstructed as |t| (P γ T ) 2 , with the resolution of 11%. The hadronic final state Y is reconstructed using a combination of tracking and calorimetric information. The event inelasticity of the γp interaction is reconstructed as y I P where the summation is performed on all detected hadronic final state particles in the event. Diffractive events are selected by requiring y I P < 0.05, ensuring a large pseudorapidity gap, ∆η, between the photon and the proton dissociative system Y , since y I P e −∆η .
After applying all selection criteria, described in a more detail in [3,6], 240 events remain in the data sample.

Models
The HERWIG 6.4 MC event generator is used to simulate the diffractive high |t| photon scattering using the LLA BFKL model [7,8,9]. At leading logarithmic accuracy there are two independent free parameters in the BFKL calculation: the value of the strong coupling α s and the scale, c, which defines the leading logarithms in the expansion of the BFKL amplitude, ln(W 2 /(c |t|)). In exclusive production of vector mesons, the scale parameter c is related to the vector meson mass. In the case of diffractive photon scattering, the unknown scale results in the absence of a normalisation prediction of the cross section. In

DIS 2009
the calculations considered here, the running of α s as a function of the scale is ignored and will henceforth be referred to as α BF KL s . The choice of α BF KL s = 0.17 is used for the simulation for this measurement.
In the asymptotic approximation of the calculations [9], the W dependence of the cross section follows a power-law, σ(W ) ∼ W 4ω0 , where the exponent is given by the choice of α BF KL s by ω 0 = (3α BF KL s /π) 4 ln 2. For the t dependence of the cross section, an approximate power-law behaviour is predicted of the form dσ/d|t| ∼ |t| −n , where n depends on the parton density functions of the proton and on the value of α BF KL s . In order to describe the data, the t dependence of the diffractive photon scattering simulation was weighted by a factor |t| 0.73 . This weighted HERWIG prediction is used only to correct the data for resolution and acceptance effects.
Possible sources of background were estimated using MC simulations and two of these sources were identified as non-negligible. The background from inclusive diffractive photoproduction events (ep → eXY , where the two hadronic final states are separated by a rapidity gap) is simulated using the PHOJET MC event generator. This background contributes when a single electromagnetic particle fakes the photon candidate and the process is estimated to amount to 3% of the measured cross section. The background from electron pair production (ep → ee + e − X) is modelled using the GRAPE event generator. This process contributes to the selection if two of the leptons fake signals from the photon candidate and the scattered electron, whereas the remaining lepton escapes detection. This background contributes 4% of the measured cross section.

Cross section measurement
The ep → eγY differential cross sections are defined using the formula: where N data is the number of observed events corrected for trigger efficiency, N bgr the expected contribution from background events, L the integrated luminosity, A the signal acceptance and ∆W and ∆t the bin widths in W and t, respectively. The γp → γY differential cross section is then extracted from the ep cross section using: where the photon flux, Γ(W ), is integrated over the range Q 2 < 0.01 GeV 2 according to the modified Weizsäcker-Williams approximation [10]. The γp cross section is obtained by modelling σ γp→γY as a power-law in W , whose parameters are iteratively adjusted to reproduce the measured W dependence of the ep cross section. The differential γp cross section in |t| is then extracted from the ep cross section by correcting for the effect of the photon flux over the visible W range. The γp cross section as a function of W is obtained by first integrating the above equation over the |t| range, and then correcting for the effect of the photon flux in each bin in W . The systematic error on the measurement stems from experimental uncertainties and from model dependences. They are calculated using the weighted HERWIG simulation of the signal process. The uncertainty on the PHOJET MC normalisation and the size of the  variation of the model dependence on x I P , |t| and M Y , are estimated from the measured distributions in the data: the variation corresponds to the a range where the model still describes the data. Each source of systematic error is varied in the weighted HERWIG Monte Carlo within its uncertainty. The largest systematic error on the cross sections comes from the uncertainty on the x I P and M Y dependences. The total systematic error on the W and |t| dependence of the cross section varies from 10% to 17% and from 8% to 14%, respectively. An additional global uncertainty of 4% arises from the γp cross section extraction procedure. The total systematic errors are comparable to or smaller than the statistical errors.

Results and conclusions
The cross sections are measured in the domain 175 < W < 247 GeV, 4 < |t| < 36 GeV 2 , y I P < 0.05 and Q 2 < 0.01 GeV 2 . The γp → γY cross section as a function of W is shown in Fig. 3a. A power-law dependence of the form σ ∼ W δ is fitted to the measured cross section, the fit yields δ = 2.73 ± 1.02 (stat.) +0.56 −0.78 (syst.) with χ 2 /n.d.f. = 2.7/2. The steep rise of the cross section with W is usually interpreted as an indication of the presence of a hard sub-process in the diffractive interaction and of the applicability of perturbative QCD. The present δ value, measured at an average |t| value of 6.1 GeV 2 , is compatible with that measured by H1 in diffractive J/ψ photoproduction of δ = 1.29 ± 0.23(stat.) ± 0.16(syst.) [11] at an average |t| of 6.9 GeV 2 . The Pomeron intercepts associated with these δ values correspond to the strongest energy dependences measured in diffractive processes.
The γp cross section differential in |t|, at W = 219 GeV, is shown in Fig. 3b, together with a fit of the form dσ/dt ∼ |t| −n , where the fitted n = 2.60 ± 0.19 (stat.) +0.03 −0.08 (syst.) with χ 2 /n.d.f. = 1.6/1. The |t| dependence is harder than that measured by H1 in the diffractive photoproduction of J/ψ mesons at large |t| [11]. The measured cross sections are compared to predictions of the LLA BFKL model, using the HERWIG Monte Carlo with no |t| weighting. The predictions are normalised to the integrated measured cross section, as the normalisation is not predicted by the LLA BFKL calculation. The W dependence of the cross section is well described by the LLA BFKL prediction. Using δ = 4 ω 0 = 4 (3α F it S /π) 4 ln 2, the measured W dependence leads to α F it S = 0.26 ± 0.10 (stat.) +0.05 −0.07 (syst.). Predictions are shown in Figure 3 also for the values α BF KL s = 0.14 and 0.37, corresponding to one standard deviation of α F it S . As shown in Figure 3b, the LLA BFKL calculation for α BF KL s = 0.14, 0.26 and 0.37, all of which give a reasonable description of the W dependence, predict steeper |t| distributions than is measured in the data. The same effect cannot be established for the exclusive ρ measurement [12], where the measured t range is limited to |t| < 8 GeV 2 , although an underestimate of the cross section was observed at the largest values of |t|. The present situation is in contrast with the analysis of J/ψ production [11,13], where the |t| dependence was found to be well described by the LLA BFKL prediction over a similar range in t.