MEASUREMENT AND PHENOMENOLOGY OF THE PROTON STRUCTURE FUNCTION F FROM ZEUS AT HERA A. Quadt Department of Physics, Particle Physics, Keble Road, Oford OX 3RH, England E-mail: quadt@mail.desy.de arxiv:hep-ph/9836v 3 Nov 998 Measurents of the proton structure function F in the Q range.6 7 GeV from ZEUS 995 shifted verte data and Q.5 GeV from 996 and 997 ZEUS data are presented. From the former and other ZEUS F data the slopes df /dlnq at fied and df /dln(/) at fied Q are derived. F data at Q below.9 GeV are described successfully by a combination of generalised vector meson dominance and Regge theory. Using a NLO QCD fit the gluon density in the proton is etracted in the range 3 5 < <.7 from ZEUS 994 and 995 data. For Q GeV it is found that the q q sea distribution is still rising at small whereas the gluon distribution id strongly suppressed. It is shown that these observations may be understood from the behaviour of the F and df /dlnq data thselves. Introduction Measurents of the low and medium Q a neutral current (NC) deep inelastic scattering (DIS) cross sections atherahaverevealedthe rapidriseofthe protonstructure function F as Bjorken- decreases below. At low Q down to. GeV ZEUS data allows study of the transition region as Q in which perturbative QCD (pqcd) must break down. At high Q, NC DIS measurents are sensitive to details of the QCD evolution of parton densities, electroweak couplings and the propagator mass of the Z gauge boson. Furthermore, such measurents allow the searches for physics beyond the Standard Model, such as resonance searches or contact interactions. Phenomenology of F at low and low Q. Phenomenology of the low Q region The primarypurposeis to usenlo DGLAP QCDon the one hand and the simplest non-perturbative models on the other to eplorethe Q transitionregionand through probing their limitations to shed light on how the pqcd description of F breaks down. One way to understand the rise in F at low is advocated by Glück, Reya and Vogt (GRV94) who argue that the starting scale for the evolution of the parton densities should be very low (.3 GeV ) and at the starting scale the parton density functions should be non-singular. The observed rise in F, with a parameterisation valid above Q GeV, is then generated dynamically. On the other hand, at low one might epect that the standard NLO Q evolution given by the DGLAP equations breaks down because of the large ln(/) terms that are not included. Such terms are taken into account by the BFKL formalism, which in a the negative of the square of the four-momentum transfer between the positron and the proton leading order predicts a rising F at low. The rise comes from a singular gluon density, g λ, with λ in the range.3 to -.5. Clearly accurate eperimental results on F at low and the implied value of λ are of great interest. At some low value of Q pqcd will break down and non-perturbative models must be used to describe the data. At low and largeγ p centre-of-massenergy, W Q /, the total γ p cross-section is given by σ γ p tot (W,Q ) = σ T +σ L = 4π α Q F (,Q ) () where σ T and σ L are the cross-sections for transversely and longitudinally polarised virtual photons respectively. Two non-perturbative approaches are considered, the generalised vector meson dominance model (GVMD) and a Regge-type two component Pomeron+Reggeon approacha la Donnachieand Landshoff(DL) to giveagood description of hadron-hadron and photoproduction total cross-section data.. Measurent of F with Shifted Verte Data The shifted verte data correspond to an integrated luminosity of 36 nb taken in a special running period, in which the nominal interaction point was offset in the proton beam direction by +7 cm, away from the detecting calorimeter. Compared to the earlier shifted verte analysis, for the 995 data taking period the calorimeter modules above and below the beam were moved closer to the beam, thus etending the shifted verte Q range down to.6 GeV. The double differential cross-section for single virtual-boson echange in DIS is given by d σ ddq = πα [ Y+ Q 4 F y F L Y F 3 ] (+δr )() ] [( y)+ y F (+δ r ), (3) +R πα Q 4
F F F F.5.5.5.5 ZEUS 995 Q =. GeV Q =.5 GeV Q =. GeV Q =.5 GeV Q =.3 GeV Q =.4 GeV Q =.5 GeV Q =.65 GeV Q =.9 GeV Q =.3 GeV Q =.9 GeV Q =.5 GeV Q = 3.5 GeV -5 - Q = 4.5 GeV -5 - Q = 6. GeV -5 - ZEUS BPC95 ZEUS SVX95 ZEUS94 E665 H SVX95 ZEUSREGGE ZEUSQCD Figure : Low Q F data for different Q bins together with the ZEUSDL style Regge model fit to the ZEUS BPC95 data. Atlarger Q values the ZEUS NLO QCD fit is also shown. where R is related to the longitudinal structure function F L by R = F L /(F F L ) and δ r gives the radiative corrections to the Born cross-section, which in this kinatic region is at most %. The parity violating term F 3 arising from the Z echange is negligible in the Q range of this analysis. Further details about the data analysis can be found in ref.. Fig. shows the results for F as a function of in bins of Q between.65 and 6 GeV (ZEUS SVX95) together with ZEUS F measurents at very low Q =..65 GeV (ZEUS BPC95) and at larger Q those from the ZEUS94. There is good agreent between the different ZEUS data sets in the region of overlap. Also shown are data from the shifted verte measurents by H (H SVX95) and fied target data from E665. The steep increase of F at low observed in the higher Q bins softens at the lower Q values of this analysis. The curves shown will be discussed later in the tet..3 The low Q region We first give an overview of the low Q region, Q < 5 GeV, taking ZEUS SVX95, BPC95 and ZEUS94 F data. Using Eq. we calculate σ γ p tot values from the F data. The DL model predicts that the cross-section rises slowly with energy W λ, λ = α P.8 and this behaviour ses to be followed by the data at very low Q. Above Q =.65 GeV, the DL model predicts a shallower rise of the cross-section than the data ehibit. For Q values of around GeV and above, the GRV94 curves describe the qualitative behaviour of the data, namely the increasing rise of σ γ p tot with W, as Q increases. This suggests that the perturbative QCD calculations can account for a significant fraction of the cross-section at the larger Q values. For the rainder of this section we concentrate on non-perturbative descriptions of the ZEUS BPC95 data. Since BPC95 data are binned in Q and y we first rewrite the double differential cross-section of Eq. 3 as d σ dydq = Γ (σ T +ǫσ L ) where σ L = Q 4π α F L and σ T has been defined by Eq.. The virtual photon has flu factor Γ and polarisation ǫ. Keeping only the continuum states in the GVMD at a fied W the longitudinal and transverse γ p cross-section are related to the corresponding photoproduction cross-section σ γp by [ M σ L (W,Q ) = ξ Q ln M +Q M ] M M σ γp +Q (W ) σ T (W,Q M ) = M +Qσγp (W ) (4) where the parameter ξ is the ratio σ V p L /σv p T for vector meson (V) proton scattering and M is the effective vector meson mass. Neither ξ nor M are given by the model and they are either determined from a fit to data or by other approaches. As we do not have much sensitivity to ξ and it is small (. -.4) we set it here to zero. We thus have 9 parameters to be determined by fitting the BPC data to the simplified GVMD epression F = Q M σ γp M +Q 4π α in 8 bins of W between 4 and 5 GeV. The fit is reasonable and its quality might also be judged from the upper plot in Fig.. The value obtained for M is.53±.4(stat)±.9(sys). The resulting etrapolated values of σ γp areshown as a function of W in the lower plot of Fig., along with measurents from HERA and lower energy eperiments. The etrapolated BPC data lie somewhat above the direct measurents from HERA. They are also above the cross section prediction of the DL model. It should be clearly understood that the σ γp data derived from the BPC95data are not a measurent of the total photoproduction cross-section but the result of a physically motivated ansatz. The simple GVMD approach just described gives a concise account of the Q dependence of the BPC data but it says nothing about the energy dependence of σ γp. To eplore this aspect of the data we turn to a two component Regge model σ γp tot (W ) = A R (W ) αr +A P (W ) αp where P and R denote the Pomeron and Reggeon contributions. The Reggeon intercept α R is fied to the value.5whichiscompatiblewiththe originaldlvalueandby the re-evaluation of Cudell et al. With such an intercept
σ T +εσ L (µb) (scaled) σ tot γp (µb) 6 4 8 6 4 4 8 6 4 ZEUS 995 W = 5 GeV ( 8) W = 33 GeV ( 7) W = GeV ( 6) W = 9 GeV ( 5) W = 73 GeV ( 4) W = 53 GeV ( 3) W = 34 GeV ( ) W = 4 GeV ( )...3.4.5.6.7.8.9 Reggeon+Pomeron Pomeron (BPC only) DL ZEUS BPC 995 (etrapolated) ZEUS, H γp low W γp Q (GeV ) 3 4 5 W (GeV ) Figure : Upper plot: ZEUS BPC95 measurents of the total cross-section σ T +ǫσ L in bins of W and the GVMD fit to the data. Lower plot: σ γp tot as a function of W. The ZEUS BPC95 points are those from the GVMD etrapolation of σ γp. the Reggeon contribution is negligible at HERA energies. Fitting the etrapolated BPC95 data alone yields a value.4±.(stat)for α P. Fitting both terms to the real photoproduction data (with W > 3 GeV ) and BPC95 data yields α P =.±.(stat). Including in addition the two original measurents from HERA as well gives α P =.±.(stat). All these values of α P are larger than the value of.8 used originally by DL, but we note that the best estimate of Cudell et al. is.964 +.5.94, which within the errors is consistent with our result. The final step in the analysisof the BPC data is to combine the GVMD fitted Q dependence with the Regge model energy dependence σ γ p tot = ( M M +Q ) (A R (W ) αr +A P (W ) αp ). The parameters M and α R are fied to their previous values of.53 and.5, respectively. The 3 raining parameters are determined by fitting to real photoproduction data and the original BPC data. The description of the lowq F data givenby this DL stylodel is shown in Fig.. Data in the BPC region Q <.65 GeV is well described. At larger Q values the curves fall below the data. Also shown in Fig. for Q > 6 GeV are the results of a NLO QCD fit (full line) as described in Sec..5..4 F slopes: dlnf /dln(/);df /dlnq To quantify the behaviour of F as a function of Q and at low we calculate the two slopes dlnf /dln(/);df /dlnq from the ZEUS SVX95, BPC95 and ZEUS94 data sets. At a fied value of Q and at small the behaviour of F can be characterised by F λ, with λ taking rather different values in the Regge and BFKL approaches. λ eff is calculated from horizontal slices of ZEUS F data between the y = HERA kinatic limit and a fied cut of <., here including E665 data. In a given Q bin is calculated from the mean value of ln(/) weighted by the statistical errors of the corresponding F values. The same procedure is applied to the theoretical curves shown for comparison. Figure3showstheasuredvaluesofλ eff asafunction of Q. From the Regge approach one would epect λ eff. and independent of Q. Data for Q < GeV is consistent with this epectation. The linked points labelled DL are calculated from the Donnachie- Landshoff fit and as epected from the discussion of the previous section are somewhat below the data. For Q > GeV, λ eff increases slowly to around.3 at Q values of 4 GeV. Qualitatively the tendency of λ eff to increase with Q is described by a number of pqcd approaches. The linked points labelled GRV94 are calculated from the NLO QCD GRV94 fit. Although the GRV94 prediction follows the trend of the data it tends to lie above the data, particularly in the Q range 3 GeV. λ eff.5.45.4.35.3.5..5..5.3.5.3.4.68.9.86 3.838 ZEUS 995.398.88.35.6.4.79.669.56.97.99.76 3.377.669 3.68 3.693 4.77 4.88 5.74 6.79 6.8-6.994 7.9 Q (GeV ) df /dlnq.6.5.4.3.. ZEUS 995 -. -6-5 -4-3 - - Figure 3: Left plot: dlnf /dln(/) as a function of Q calculated by fitting ZEUS and E665 F data in bins of Q. Right plot: df /dlnq as a function of calculated by fitting ZEUS F data in bins of. Within the framework of pqcd, at small the behaviour of F is largely determined by the behaviour of the sea quarks F S, whereas the df /dlnq is determined by the convolution of the splitting function P qg and the gluon density, df /dlnq α s P qg g. In order to study the scaling violations of F in more detail the logarithmic slope df /dlnq is derived from the data by fitting F = a+blnq in bins of fied. The statistical and systatic errors are determined as described above. 3
The results for df /dlnq as a function of are shown in Fig. 3. For values of down to 3 4, the slopes are increasing as decreases. At lower values of and Q, the slope decreases. Comparing the rapid increase in F at small with the behaviour of the df /dlnq, one is tpted to the naive conclusion that the underlying behaviour of the sea quark and gluon momentum distributions must be different at small, with the sea dominant and the gluon tending to zero. The failure of DL is in line with the earlier discussion. GRV94 does not follow the trend of the data when it turns over..5 NLO QCD fit to F data In perturbative QCD the scaling violations of the F structure function are caused by gluon brsstrahlung from quarks and quark pair creation from gluons. In the low domain accessible at HERA the latter process dominates the scaling violations. A QCD analysis of F structure functions measured at HERA therefore allows one to etract the gluon momentum density in the proton down to low values of. In this section we present NLO QCD fits to the ZEUS 994 nominal verte data and the SVX95 data of this paper. We are not attpting to include all available information on parton densities, but concentrating on what ZEUS data and their errors allow us to conclude about the gluon density at low. To constrain the fits at high proton and deuteron F structure function data from NMC and BCDMS are included. The kinatic range covered in this analysis is 3 5 < <.7 and < Q < 5 GeV. The QCD predictions for the F structure functions are obtained by solving the DGLAP evolution equations at NLO. These equations yield the quark and gluon momentum distributions at all values of Q provided they are given at some input scale Q. In this analysis we adopt the so-called fied flavour number sche where only three light flavours (u,d,s) contribute to the quark density in the proton. The corresponding structure functions F c and Fb are calculated from the photon-gluon fusion process including massive NLO corrections. The input valence distributions are taken from the parton distribution set MRS(R). As for MRS(R) we assume that the strange quark distribution is a given fraction K s =. of the sea at the scale Q = GeV. The gluon normalisation is fied by the momentum sum rule. The input value for the strong coupling constant is set to α s (MZ ) =.8 and the charm mass is taken to be m c =.5 GeV. In the QCD evolutions and the evaluation of the structure functions the renormalisation scale and mass factorisation scale are both set equal to Q. In the definition of the χ only statistical error are included and the relative normalisation of the data sets is fied at unity. The fit yields a good description of the data as shown in Fig.. We have also checked that the gluon obtained from this fit to scaling violations is in agreent with the recent ZEUS measurents of charm production and F c in deep inelastic scattering at HERA. Two types of systatic uncertainties have been considered in this analysis. HERA standard errors contain statistical error on the data, eperimental systatic uncertainties, relative normalisation of the different data sets and uncertainties on α s, the strange quark content of the proton and the charm mass. Parametrisation errors contain uncertainties from a χ definition including statistical and eperimental systatic errors, variations of the starting scale Q and an alternative, more fleible parametrisation of the gluon density using Chebycheff polynomials. The first type of errors amounts to 6% g/g at = 5 5, Q = 7 GeV, the second type yields 9.5% in g/g. 6 4 3 5 Q = GeV min Q = GeV g Σ ZEUS 995 NLO(MS) Q = 4 GeV min Q = GeV Q = 7 GeV Q = 7 GeV Q = GeV -4-3 - - g Σ Q = GeV NLO(MS) -4-3 - - Figure 4: The quark singlet momentum distribution, Σ (shaded), and the gluon momentum distribution, g() (hatched), as a function of at fied values of Q =, 7 and GeV. The error bands correspond to the quadratic sum of all error sources considered for each parton density. The three plots of Fig. 4 show the distribution for Σ and g as a function of for Q at, 7 and GeV. It can be seen that even at the smallest Q Σ is rising at small whereas the gluon distribution has become almost flat. These results give support to the naive conclusion of Sec..4, that the sea distribution dominates at low and Q. At Q = GeV the gluon distribution is poorly determined and can, within errors, be negative at low. 4
3 Measurent of the Proton Structure Function F from 996 and 997 data 3. Kinatics in Deep Inelastic Scattering Recalling the double differential NC cross-section(3), but nowincluding the corrections(δ L and δ 3 ) forf L and F 3 yields d σ ddq = πα Y + Q 4 F ( δ L δ 3 )(+δ r ) (5) Here the F structure function contains contributions from virtual photon and Z echange F overlap with the fied target eperiments; in the overlap region reasonable agreent has been found. The F scaling violation from this analysis and the fied target data are also shown in Figure 5. For Q > GeV the increase in statistics allows a measurent of F in smaller bins with respect to our previous measurent. Above Q = 8 GeV, the statistical error grows typically to 5-5% and dominates the total error. Overall our data are in agreent with our published data. ZEUS Preliminary 996-97 Q =.5 ZEUS 996-97 ZEUS 994 BCDMS, E665, NMC, SLAC Q =.7 Q = 3.5 F ZEUS Preliminary 996-97 Q = 35 Q = 45 Q = 6 F = F + Q + (Q +M Z )Fint Q 4 (Q +MZ )Fwk (6) where M Z is the mass of the Z and F, F wk and F int are the contributions to F due to photon echange, Z echange and γz interference respectively. In this analysiswedeterminedthestructurefunctionf using 996 and 997 data with an integrated luminosity of 6.8 pb and 7.4 pb, respectively. The selection and kinatic reconstruction of NC DIS events is based on an observed positron and the hadronic final state. For further details see ref.. F.5. -5-4 -3 - - Q = 4.5 Q = 6.5 Q = 8.5 Q = Q = Q = 5 Q = 8 Q =. -5-4 -3 - - Q = 7-5 -4-3 - - ZEUS Preliminary 996-97 Q = 8 Q = Q = 5 Q = 7 Q = 9 Q = Q = 5 Q = Q = 5 Q = 35-5 -4-3 - - Q = 45-5 -4-3 - - Q = 65-5 -4-3 - - 3. Results Monte Carlo samples are used to estimate the acceptance, migration, radiative corrections, electroweak corrections and background contributions. F is then determined based on a bin-by-bin unfolding. The resulting statistical error, including the Monte Carlo statistics, ranges from % below Q = GeV to 5-6% at Q 8 GeV. The systatic uncertainties have been estimated by varying the selection cuts, efficiencies and reconstruction techniques and redetermining the cross section including background estimates. Potential error source such as possible detector misalignment, event verte reconstruction, calorimeter energy scale, positron identification efficiency, background contributions and hadronic energy flow have been considered. The total systatic uncertainty amounts to 3-4% ecept at low and high y, where it grows to %. At the present preliminary state of the analysis we estimate an overall normalisation uncertainty of 3%. The resulting F is shown as a function of for fied Q in Figure 5. Results from our previous analysis, and from fied target eperiments are also shown for comparison. At low Q the rise in F for is measured with improved precision. The coverage in has also been etended to higher, yielding etended.5..5..5 Q = Q = 3 Q = 5 Q = 8 Q = -3 - - -3 - - Q = -3 - - Figure 5: Top and bottom left plots: F versus for fied Q. Bottom right plot: F as a function of Q for fied. References [] ZEUS Results on the Measurent and Phenomenology of F at Low and Low Q, DESY 98- (August 998) submitted to The European Physical Journal. [] Measurent of the Proton Structure Function F in e + p Collisions at HERA, Submitted paper to the XXIX International Conference on High Energy Physics, Vancouver, July 3-9, 998. 5