The production of beauty and charm quarks in ep interactions has been studied with the ZEUS detector at HERA for exchanged four-momentum squared 5 < Q^2 < 1000 GeV^2 using an integrated luminosity of 354 pb^{-1}. The beauty and charm content in events with at least one jet have been extracted using the invariant mass of charged tracks associated with secondary vertices and the decay-length significance of these vertices. Differential cross sections as a function of Q^2, Bjorken x, jet transverse energy and pseudorapidity were measured and compared with next-to-leading-order QCD calculations. The beauty and charm contributions to the proton structure functions were extracted from the double-differential cross section as a function of x and Q^2. The running beauty-quark mass, m_b at the scale m_b, was determined from a QCD fit at next-to-leading order to HERA data for the first time and found to be 4.07 \pm 0.14 (fit} ^{+0.01}_{-0.07} (mod.) ^{+0.05}_{-0.00} (param.) ^{+0.08}_{-0.05} (theo) GeV.

The reduced cross sections for $e^{+}p$ deep inelastic scattering have been measured with the ZEUS detector at HERA at three different centre-of-mass energies, $318$, $251$ and $225$ GeV. The cross sections, measured double differentially in Bjorken $x$ and the virtuality, $Q^2$, were obtained in the region $0.13\ \leq\ y\ \leq\ 0.75$, where $y$ denotes the inelasticity and $5\ \leq\ Q^2\ \leq\ 110$ GeV$^2$. The proton structure functions $F_2$ and $F_L$ were extracted from the measured cross sections.

The neutral current e+/-p cross section has been measured up to values of Bjorken x of approximately 1 with the ZEUS detector at HERA using an integrated luminosity of 187 inv. pb of e-p and 142 inv. pb of e+p collisions at sqrt(s) = 318GeV. Differential cross sections in x and Q2, the exchanged boson virtuality, are presented for Q2 geq 725GeV2. An improved reconstruction method and greatly increased amount of data allows a finer binning in the high-x region of the neutral current cross section and leads to a measurement with much improved precision compared to a similar earlier analysis. The measurements are compared to Standard Model expectations based on a variety of recent parton distribution functions.

The production of $D^{*\pm}$ mesons in deep inelastic $ep$ scattering has been measured for exchanged photon virtualities $ 5<Q^2<1000 \gev^2 $, using an integrated luminosity of 363 pb$^{-1}$ with the ZEUS detector at HERA. Differential cross sections have been measured and compared to next-to-leading-order QCD calculations. The cross-sections are used to extract the charm contribution to the proton structure functions, expressed in terms of the reduced charm cross section, $\sigma_{\rm red}^{c\bar{c}}$. Theoretical calculations based on fits to inclusive HERA data are compared to the results.

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.

Measurements of neutral current cross sections for deep inelastic scattering in e+p collisions at HERA with a longitudinally polarised positron beam are presented. The single-differential cross-sections d(sigma)/dQ2, d(sigma)/dx and d(sigma)/dy and the reduced cross-section were measured in the kinematic region Q2 > 185 GeV2 and y < 0.9, where Q2 is the four-momentum transfer squared, x the Bjorken scaling variable, and y the inelasticity of the interaction. The measurements were performed separately for positively and negatively polarised positron beams. The measurements are based on an integrated luminosity of 135.5 pb-1 collected with the ZEUS detector in 2006 and 2007 at a centre-of-mass energy of 318 GeV. The structure functions F3 and F3(gamma)Z were determined by combining the e+p results presented in this paper with previously published e-p neutral current results. The asymmetry parameter A+ is used to demonstrate the parity violation predicted in electroweak interactions. The measurements are well described by the predictions of the Standard Model.

For the first time, differential inclusive-jet cross sections have been measured in neutral current deep inelastic ep scattering using the anti-kT and SIScone algorithms. The measurements were made for boson virtualities Q^2 > 125 GeV^2 with the ZEUS detector at HERA using an integrated luminosity of 82 pb^-1 and the jets were identified in the Breit frame. The performance and suitability of the jet algorithms for their use in hadron-like reactions were investigated by comparing the measurements to those performed with the kT algorithm. Next-to-leading-order QCD calculations give a good description of the measurements. Measurements of the ratios of cross sections using different jet algorithms are also presented; the measured ratios are well described by calculations including up to O(alphas^3) terms. Values of alphas(Mz) were extracted from the data; the results are compatible with and have similar precision to the value extracted from the kT analysis.

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.

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.

The proton-dissociative diffractive photoproduction of J/psi mesons has been studied in ep collisions with the ZEUS detector at HERA using an integrated luminosity of 112 pb^-1. The cross section is presented as a function of the photon-proton centre-of-mass energy and of the squared four-momentum transfer at the proton vertex. The results are compared to perturbative QCD calculations.