# Measurement of 𝐷^{∗±} Cross Sections and the Charm Contribution to the Structure Function of the Proton in Deep Inelastic Scattering at HERA

## 1 INTRODUCTION

The first HERA measurements of the charm contribution to the proton structure function, $`F_2^{c\overline{c}}`$ were reported by the H1 and ZEUS Collaborations from an analysis of $`D^\pm `$ production in their 1994 data sets . The results were consistent with Photon Gluon Fusion (PGF) being the dominant mechanism for $`D^\pm `$ production in $`e^+p`$ Deep Inelastic Scattering (DIS). If this is the case, this type of measurements are sensitive to the gluons in the proton. In addition, they can provide a test of the universality of the parton distribution functions (pdf’s), namely, whether pdf’s extracted from the inclusive measurement of the proton structure function, $`F_2`$, can be used as input for calculations of more exclusive processes as charm production.

Here we present a study of $`D^\pm `$ production using the 1996 and 1997 data corresponding to an integrated luminosity of 37 pb<sup>-1</sup>. More than tenfold larger data sample, together with the modifications of the ZEUS detector made for the 1996 and 1997 operation allow an extension of the kinematic range to both smaller and larger $`Q^2`$. The $`D^\pm `$ is tagged via $`D^+(D^0K^{}\pi ^+)\pi ^+`$(+ c.c.) and $`D^+(D^0K^{}\pi ^+\pi ^+\pi ^{})\pi ^+`$(+ c.c.) decay channels, refered to as $`K2\pi `$ and $`K4\pi `$ respectively.

## 2 $`D^\pm `$ CROSS SECTIONS

The measured $`D^\pm `$ cross section using the $`K2\pi (K4\pi )`$ final state, in the region $`1.5(2.5)<p_\mathrm{T}(D^{})<15GeV`$, $`|\eta (D^{})|<1.5`$ is $`\sigma `$$`\left(e^+pe^+D^\pm X\right)`$ $`=8.31\pm 0.31\left(stat\right)_{0.50}^{+0.30}\left(sys\right)`$ nb($`3.65\pm 0.36(stat)_{0.41}^{+0.20}(sys)`$ nb). Figure 1 shows the differential $`D^{}`$ cross sections in the restricted $`Q^2`$, $`y`$, $`p_\mathrm{T}(D^{})`$ and $`\eta (D^{})`$ region as functions of $`log_{10}(Q^2)`$, $`log_{10}(x)`$, $`W`$, $`p_\mathrm{T}(D^{})`$, $`\eta (D^{})`$ and $`x(D^{})`$$`=2p^{}(D^{})/W`$, where $`p^{}(D^{})`$ is the momentum in the $`\gamma ^{}`$-proton CMS frame. The results using each decay channel can be directly compared in the $`p_\mathrm{T}(D^{})`$ differential cross section. The agreement is satisfactory.

## 3 COMPARISON WITH NLO QCD

We compute NLO QCD calculations with a semi-inclusive Monte Carlo generator HVQDIS for heavy quark production and subsequent fragmentation to $`D^\pm `$ via a Peterson fragmentation function  . This generator is based on NLO calculations  in the three flavor number scheme (TFNS), in which only light quarks ($`u,d,s`$) are included in the initial state proton. Heavy quarks are produced exclusively by the convolution of the light flavours and the gluon with the massive matrix elements and coefficient functions calculated previously .

We use as input pdf ZEUS94 . The QCD renormalization and factorization scales are set to $`\sqrt{Q^2+4m_{c}^{}{}_{}{}^{2}}`$. $`m_c`$ is varied between 1.3 and 1.5 GeV. $`f(cD^{})=0.222`$ is taken from $`e^+e^{}`$ measurements . The error on this quantity introduces a normalization uncertainty of $`9`$%. Finally, the Peterson fragmentation parameter is set to $`ϵ=0.035`$. The NLO QCD predictions (figure 1) are in reasonable agreement with the data, except in the $`\eta (D^{})`$ distribution, where the measurements show a shift into the positive $`\eta `$ region (proton direction) with respect to the prediction. Also a softer charm fragmentation is favored by the data.

### 3.1 Fragmentation Effects

In Monte Carlo fragmentation models like JETSET or HERWIG, a forward shift in the $`D^{}`$ direction w.r.t. that of the original charm quark is produced during the fragmentation due to the interaction of the charm quark with the proton remnant via either strings or soft gluon radiatio

To investigate how this affects the NLO QCD predictions we have reweighted a LO Monte Carlo for charm production, RAPGAP (it uses JETSET for the fragmentation) in a way such that at the stage of the hard interaction it will reproduce exactly the HVQDIS results for the $`p_t(c)`$ and $`\eta (c)`$ differential cross sections. The predictions from this NLO reweighted RAPGAP Monte Carlo are shown in (figure 1). They provide a better descripion of the data, especially in the $`\eta (D^{})`$ and $`x(D^{})`$ differential cross sections This result suggest that the small disagreement found with HVQDIS come from the fact that the Peterson function can not account for all the charm quark fragmentation effects present at HERA, in particular the interaction with the remnant.

## 4 $`F_2^{c\overline{c}}`$ EXTRACTION

The procedure for extracting $`F_2^{c\overline{c}}`$ starts by measuring the $`D^\pm `$ production cross section in the restricted $`p_\mathrm{T}(D^{})`$, $`\eta (D^{})`$ region in bins of $`Q^2`$ and $`y`$. The extrapolation to the full $`p_\mathrm{T}(D^{})`$, $`\eta (D^{})`$ region is done in the following way:

$$F_{2}^{c\overline{c}}{}_{}{}^{\text{m}}(x_i,Q_i^2)=\frac{\sigma _i^\text{m}(e^+pD^{}X)}{\sigma _i^\text{t}(e^+pD^{}X)}F_{2}^{c\overline{c}}{}_{}{}^{\text{t}}(x_i,Q_i^2)$$
(1)   
where $`x_i`$, $`Q_i^2`$ is the center of gravity of the bin $`i`$ and ‘m’ and ‘t’ denote ‘measured’ and ‘theoretical’ respectively. $`F_{2}^{c\overline{c}}{}_{}{}^{t}`$ is taken from ZEUS94  fit. For $`\sigma _i^t`$ we use the reweighted Monte Carlo. A number of assumptions are implicitly done in this procedure:

* The TFNS is valid,
* $`F_L^{c\overline{c}}`$ is negligible($`<1\%`$ of $`F_2`$ in our $`Q^2`$, $`y`$ region from calculations based on ),
* the value of $`f(cD^{})`$ measured in $`e^+e^{}`$ is valid also at HERA,
* the cross section outside the restricted region is well described by NLO QCD.

Figure 2 shows the measured $`F_2^{c\overline{c}}`$ after combining the results from both decay channels. Compared to our previous study we have extended the kinematic range to $`Q^2`$ as low as 1.8 GeV<sup>2</sup> and up to 130 GeV<sup>2</sup> and the errors are reduced substantialy. $`F_2^{c\overline{c}}`$ exhibits a steep rise with decreasing $`x`$ at constant $`Q^2`$. From a comparison with the ZEUS94 parametrization we determine that $`F_2^{c\overline{c}}`$ accounts for $`<10`$% of $`F_2`$ at low $`Q^2`$ and $`x510^4`$ and $`30`$% of $`F_2`$ for $`Q^2`$$`>11`$ GeV<sup>2</sup> at the lowest $`x`$ measured.

## 5 SUMMARY

We have presented a charm analysis in DIS using the combined 1996 and 1997 data sample. Charm was tagged with $`D^{}`$ mesons decaying into two decays ($`K2\pi `$ and $`K4\pi `$). In the experimentally accessible region, differential $`D^\pm `$ cross sections are in reasonable agreement with NLO QCD calculations of charm production in the TFNS using a pdf extracted from an inclusive measurement of $`F_2`$. This represents a successful test of the universality of the pdf’s.

Small disagreements in the $`\eta (D^{})`$ and $`x(D^{})`$ distributions show that the fragmentation a la Peterson can not account for all the charm quark fragmentation effects present at HERA.

Using these calculations to extrapolate outside the measured $`p_\mathrm{T}(D^{})`$, $`\eta (D^{})`$ region, $`F_2^{c\overline{c}}`$ was extracted. $`F_2^{c\overline{c}}`$ is rising steeply with decreasing $`x`$ at costant $`Q^2`$. It amounts to $`25\%`$ of $`F_2`$ at low $`x`$, $`Q^2>10`$.