The cross section for high-E_T dijet production in photoproduction has been measured with the ZEUS detector at HERA using an integrated luminosity of 81.8 pb-1. The events were required to have a virtuality of the incoming photon, Q^2, of less than 1 GeV^2 and a photon-proton centre-of-mass energy in the range 142 < W < 293 GeV. Events were selected if at least two jets satisfied the transverse-energy requirements of E_T(jet1) > 20 GeV and E_T(jet2) > 15 GeV and pseudorapidity requirements of -1 < eta(jet1,2) < 3, with at least one of the jets satisfying -1 < eta(jet) < 2.5. The measurements show sensitivity to the parton distributions in the photon and proton and effects beyond next-to-leading order in QCD. Hence these data can be used to constrain further the parton densities in the proton and photon.

The cross sections for inelastic photoproduction of J/psi and psi^prime mesons have been measured in ep collisions with the ZEUS detector at HERA, using an integrated luminosity of 468 pb-1 collected in the period 1996--2007. The psi^prime to J/psi cross section ratio was measured in the range 0.55 < z < 0.9 and 60 < W < 190 GeV as a function of W, z and p_T. Here W denotes the photon-proton centre-of-mass energy, z is the fraction of the incident photon energy carried by the meson and p_T is the transverse momentum of the meson with respect to the beam axis. The J/psi cross sections were measured for 0.1 < z < 0.9, 60 < W < 240 GeV and p_T > 1 GeV. Theoretical predictions within the non-relativistic QCD framework including NLO colour--singlet and colour--octet contributions were compared to the data, as were predictions based on the k_T--factorisation approach.

Charm production in deep inelastic scattering has been measured with the ZEUS detector at HERA using an integrated luminosity of 82 pb^{-1}. Charm has been tagged by reconstructing D^{*+}, D^0, D^{+} and D_s^+ (+ c.c.) charm mesons. The charm hadrons were measured in the kinematic range p_T(D^{*+},D^0,D^{+}) > 3 GeV, p_T(D_s^+)>2 GeV and |\eta(D)| < 1.6 for 1.5 < Q^2 < 1000 GeV^2 and 0.02 < y < 0.7. The production cross sections were used to extract charm fragmentation ratios and the fraction of c quarks hadronising into a particular charm meson in the kinematic range considered. The cross sections were compared to the predictions of next-to-leading-order QCD, and extrapolated to the full kinematic region in p_T(D) and \eta(D) in order to determine the open-charm contribution, F_2^{c\bar{c}}(x,Q^2), to the proton structure function F_2.

Inclusive jet cross sections in photoproduction for events containing a $D^*$ meson have been measured with the ZEUS detector at HERA using an integrated luminosity of $78.6 {\rm pb}^{-1}$. The events were required to have a virtuality of the incoming photon, $Q^2$, of less than 1 GeV$^2$, and a photon-proton centre-of-mass energy in the range $130&lt;W_{\gamma p}&lt;280 {\rm GeV}$. The measurements are compared with next-to-leading-order (NLO) QCD calculations. Good agreement is found with the NLO calculations over most of the measured kinematic region. Requiring a second jet in the event allowed a more detailed comparison with QCD calculations. The measured dijet cross sections are also compared to Monte Carlo (MC) models which incorporate leading-order matrix elements followed by parton showers and hadronisation. The NLO QCD predictions are in general agreement with the data although differences have been isolated to regions where contributions from higher orders are expected to be significant. The MC models give a better description than the NLO predictions of the shape of the measured cross sections.

Mean values and differential distributions of event-shape variables have been studied in neutral current deep inelastic scattering using an integrated {luminosity} of 82.2 pb$^{-1}$ collected with the ZEUS detector at HERA. The kinematic range was $80 &lt; Q^2 &lt; 20 480\gev^2$ and $0.0024 &lt; x &lt; 0.6$, where $Q^2$ is the virtuality of the exchanged boson and $x$ is the Bjorken variable. The data are compared with a model based on a combination of next-to-leading-order QCD calculations with next-to-leading-logarithm corrections and the Dokshitzer-Webber non-perturbative power corrections. The power-correction method provides a reasonable description of the data for all event-shape variables studied. Nevertheless, the lack of consistency of the determination of $\alpha_s$ and of the non-perturbative parameter of the model, $\albar$, suggests the importance of higher-order processes that are not yet included in the model.

The dissociation of virtual photons, $\gamma^{\star} p \to X p$, in events with a large rapidity gap between $X$ and the outgoing proton, as well as in events in which the leading proton was directly measured, has been studied with the ZEUS detector at HERA. The data cover photon virtualities $Q^2>2$ GeV$^2$ and $\gamma^{\star} p$ centre-of-mass energies $40&lt;W&lt;240$ GeV, with $M_X>2$ GeV, where $M_X$ is the mass of the hadronic final state, $X$. Leading protons were detected in the ZEUS leading proton spectrometer. The cross section is presented as a function of $t$, the squared four-momentum transfer at the proton vertex and $\Phi$, the azimuthal angle between the positron scattering plane and the proton scattering plane. It is also shown as a function of $Q^2$ and $\xpom$, the fraction of the proton's momentum carried by the diffractive exchange, as well as $\beta$, the Bjorken variable defined with respect to the diffractive exchange.

Exclusive rho^0 electroproduction at HERA has been studied with the ZEUS detector using 120 pb^{-1} of integrated luminosity collected during 1996-2000. The analysis was carried out in the kinematic range of photon virtuality 2 < Q^2 < 160 GeV$^2, and gamma^* p centre-of-mass energy 32 < W < 180 GeV. The results include the Q^2 and W dependence of the gamma^* p --> rho^0 p cross section and the distribution of the squared-four-momentum transfer to the proton. The helicity analysis of the decay-matrix elements of the rho^0 was used to study the ratio of the gamma^* p cross section for longitudinal and transverse photon as a function of Q^2 and W. Finally, an effective Pomeron trajectory was extracted. The results are compared to various theoretical predictions.

Charged particle production has been studied in neutral current deep inelastic ep scattering with the ZEUS detector at HERA using an integrated luminosity of 0.44 fb^-1. Distributions of scaled momenta in the Breit frame are presented for particles in the current fragmentation region. The evolution of these spectra with the photon virtuality, Q^2, is described in the kinematic region 10&lt;Q^2&lt;41000 GeV^2. Next-to-leading-order and modified leading-log-approximation QCD calculations as well as predictions from Monte Carlo models are compared to the data. The results are also compared to e+e- annihilation data. The dependences of the pseudorapidity distribution of the particles on Q^2 and on the energy in the \gamma p system, W, are presented and interpreted in the context of the hypothesis of limiting fragmentation.

The semi-inclusive reaction e+ p -> e+ X p was studied with the ZEUS detector at HERA using an integrated luminosity of 12.8 pb-1. The final-state proton, which was detected with the ZEUS leading proton spectrometer, carried a large fraction of the incoming proton energy, xL>0.32, and its transverse momentum squared satisfied pT^2<0.5 GeV^2/ the exchanged photon virtuality, Q^2, was greater than 3 GeV^2 and the range of the masses of the photon-proton system was 45<W<225 GeV. The leading proton production cross section and rates are presented as a function of xL, pT^2, Q^2 and the Bjorken scaling variable, x.

A combination is presented of the inclusive deep inelastic cross sections measured by the H1 and ZEUS Collaborations in neutral and charged current unpolarised ep scattering at HERA during the period 1994-2000. The data span six orders of magnitude in negative four-momentum-transfer squared, Q^2, and in Bjorken x. The combination method used takes the correlations of systematic uncertainties into account, resulting in an improved accuracy. The combined data are the sole input in a NLO QCD analysis which determines a new set of parton distributions HERAPDF1.0 with small experimental uncertainties. This set includes an estimate of the model and parametrisation uncertainties of the fit result.