The production of D*+-(2010) mesons in deep inelastic scattering has been measured in the ZEUS detector at HERA using an integrated luminosity of 37 pb^-1. The decay channels D*+ -> D0 pi+(+c.c.), with D0 -> K- pi+ or D0 ->K- pi- pi+ pi+, have been used to identify the D mesons. The e+p cross section for inclusive D*+- production with 1<Q^2<600 GeV^2 and 0.02<y<0.7 is 8.31 +- 0.31(stat.) +0.30-0.50(syst.) nb in the kinematic region 1.5< pT(D*+-)<15 GeV and |eta(D*+-)|<1.5. Differential cross sections are consistent with a next-to-leading-order perturbative-QCD calculation when using charm-fragmentation models which take into account the interaction of the charm quark with the proton remnant. The observed cross section is extrapolated to the full kinematic region in pT(D*+-) and eta(D*+-) in order to determine the charm contribution, F^ccbar_2(x,Q^2), to the proton structure function. The ratio F^ccbar_2/F_2 rises from ~10% at Q^2 ~1.8 GeV^2 to ~30% at Q^2 ~130 GeV^2 for x values in the range 10^-4 to 10-3.

Three- and four-jet final states have been measured in photoproduction at HERA using the ZEUS detector with an integrated luminosity of 121 pb^-1. The results are presented for jets with transverse energy E_T^jet>6 GeV and pseudorapidity |eta^jet|<2.4, in the kinematic region given by the virtuality of the photon Q^2<1 GeV^2 and the inelasticity 0.2<y<0.85 and in two mass regions defined as 25<M_nj<50 GeV and M_nj>50 GeV, where M_nj is the invariant mass of the n-jet system. The four-jet photoproduction cross section has been measured for the first time and represents the highest-order process studied at HERA. Both the three- and four-jet cross sections have been compared with leading-logarithmic parton-shower Monte Carlo models, with and without multi-parton interactions. The three-jet cross sections have been compared to an order(alpha alpha_s^2) perturbative QCD calculation.

The forward-jet cross section in deep inelastic ep scattering has been measured using the ZEUS detector at HERA with an integrated luminosity of 6.36 pb^-1. The jet cross section is presented as a function of jet transverse energy squared, E(T,jet)^2, and Q^2 in the kinematic ranges 10^-2<E(T,jet)^2/Q^2<10^2 and 2.5 10^-4<x<8.0 10^-2. Since the perturbative QCD predictions for this cross section are sensitive to the treatment of the log(E_T/Q)^2 terms, this measurement provides an important test. The measured cross section is compared to the predictions of a next-to-leading order pQCD calculation as well as to various leading-order Monte Carlo models. Whereas the predictions of all models agree with the measured cross section in the region of small E(T,Jet)^2/Q^2, only one model, which includes a resolved photon component, describes the data over the whole kinematic range.

Inclusive photoproduction of D*+- mesons has been measured for photon-proton centre-of-mass energies in the range 130 < W < 280 GeV and a photon virtuality Q^2 < 1 GeV^2. The data sample used corresponds to an integrated luminosity of 37 pb^-1. Total and differential cross sections as functions of the D* transverse momentum and pseudorapidity are presented in restricted kinematical regions and the data are compared with next-to-leading order (NLO) perturbative QCD calculations using the "massive charm" and "massless charm" schemes. The measured cross sections are generally above the NLO calculations, in particular in the forward (proton) direction. The large data sample also allows the study of dijet production associated with charm. A significant resolved as well as a direct photon component contribute to the cross section. Leading order QCD Monte Carlo calculations indicate that the resolved contribution arises from a significant charm component in the photon. A massive charm NLO parton level calculation yields lower cross sections compared to the measured results in a kinematic region where the resolved photon contribution is significant.

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 dijet and trijet production in deep inelastic $ep$ scattering has been measured for $10&lt;Q^2&lt;100$ GeV$^2$ and low Bjorken $x$, $10^{-4}&lt;x_{\rm Bj}&lt;10^{-2}$. The data were taken at the HERA $ep$ collider with centre-of-mass energy $\sqrt{s} = 318 \gev$ using the ZEUS detector and correspond to an integrated luminosity of $82 {\rm pb}^{-1}$. Jets were identified in the hadronic centre-of-mass (HCM) frame using the $k_{T}$ cluster algorithm in the longitudinally invariant inclusive mode. Measurements of dijet and trijet differential cross sections are presented as functions of $Q^2$, $x_{\rm Bj}$, jet transverse energy, and jet pseudorapidity. As a further examination of low-$x_{\rm Bj}$ dynamics, multi-differential cross sections as functions of the jet correlations in transverse momenta, azimuthal angles, and pseudorapidity are also presented. Calculations at $\mathcal{O}(\alpha_{s}^3)$ generally describe the trijet data well and improve the description of the dijet data compared to the calculation at $\mathcal{O}(\alpha_{s}^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.

The production of the neutral strange hadrons $K^{0}_{S}$, $\Lambda$ and $\bar{\Lambda}$ has been measured in $ep$ collisions at HERA using the ZEUS detector. Cross sections, baryon-to-meson ratios, relative yields of strange and charged light hadrons, $\Lambda$ ($\bar{\Lambda}$) asymmetry and polarization have been measured in three kinematic regions: $Q^2 > 25 \gev^2$: $5 < Q^2 < 25 \gev^2$: and in photoproduction ($Q^2 \simeq 0$). In photoproduction the presence of two hadronic jets, each with at least $5 \gev$ transverse energy, was required. The measurements agree in general with Monte Carlo models and are consistent with measurements made at $e^+ e^-$ colliders, except for an enhancement of baryon relative to meson production in photoproduction.

Jet substructure and differential cross sections for jets produced in the photoproduction and deep inelastic ep scattering regimes have been measured with the ZEUS detector at HERA using an integrated luminosity of 82.2 pb-1. The substructure of jets has been studied in terms of the jet shape and subjet multiplicity for jets with transverse energies Et(jet) > 17 GeV. The data are well described by the QCD calculations. The jet shape and subjet multiplicity are used to tag gluon- and quark-initiated jets. Jet cross sections as functions of Et(jet), jet pseudorapidity, the jet-jet scattering angle, dijet invariant mass and the fraction of the photon energy carried by the dijet system are presented for gluon- and quark-tagged jets. The data exhibit the behaviour expected from the underlying parton dynamics. A value of alphas(Mz) of alphas(Mz) = 0.1176 +-0.0009(stat.) -0.0026 +0.0009 (exp.) -0.0072 +0.0091 (th.) was extracted from the measurements of jet shapes in deep inelastic scattering.

Cross sections for e^-p neutral current deep inelastic scattering have been measured at a centre-of-mass energy of 318 GeV using an integrated luminosity of 15.9 pb^-1 collected with the ZEUS detector at HERA. Results on the double-differential cross-section d^2s/dxdQ^2 in the range 185 < Q^2 < 50000 GeV^2 and 0.0037 < x < 0.75, as well as the single-differential cross-sections ds/dQ^2, ds/dx and ds/dy for Q^2 > 200 GeV^2, are presented. To study the effect of Z-boson exchange, ds/dx has also been measured for Q^2 > 10000 GeV^2. The structure function xF_3 has been extracted by combining the e^-p results presented here with the recent ZEUS measurements of e^+p neutral current deep inelastic scattering. All results agree well with the predictions of the Standard Model.