The beauty production cross section for deep inelastic scattering events with at least one hard jet in the Breit frame together with a muon has been measured, for photon virtualities Q^2 > 2 GeV^2, with the ZEUS detector at HERA using integrated luminosity of 72 pb^-1. The total visible cross section is sigma_b-bbar (ep -> e jet mu X) = 40.9 +- 5.7 (stat.) +6.0 -4.4 (syst.) pb. The next-to-leading order QCD prediction lies about 2.5 standard deviations below the data. The differential cross sections are in general consistent with the NLO QCD predictions: however at low values of Q^2, Bjorken x, and muon transverse momentum, and high values of jet transverse energy and muon pseudorapidity, the prediction is about two standard deviations below the data.

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.

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.

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.

Diffractive photoproduction of D*+/-(2010) mesons was measured with the ZEUS detector at the ep collider HERA, using an integrated luminosity of 78.6 pb^{-1}. The D* mesons were reconstructed in the kinematic range: transverse momentum p_T(D*) > 1.9 GeV and pseudorapidity |eta(D*)| < 1.6, using the decay D*+ -> D0 pi+_s followed by D0 -> K- pi+ (+c.c.). Diffractive events were identified by a large gap in pseudorapidity between the produced hadronic state and the outgoing proton. Cross sections are reported for photon-proton centre-of-mass energies in the range 130 < W < 300 GeV and for photon virtualities Q^2 < 1 GeV^2, in two ranges of the Pomeron fractional momentum x_pom < 0.035 and x_pom < 0.01. The relative contribution of diffractive events to the inclusive D*+/-(2010) photoproduction cross section is about 6%. The data are in agreement with perturbative QCD calculations based on various parameterisations of diffractive parton distribution functions. The results are consistent with diffractive QCD factorisation.

Charm production in deep inelastic scattering has been measured with the ZEUS detector at HERA using an integrated luminosity of 120 pb^{-1}. The hadronic decay channels D^{+} -> K^{0}_{S} pi^{+}, Lambda_{c}^{+} -> p K^{0}_{S} and Lambda_{c}^{+} -> Lambda pi^{+}, and their charge conjugates, were reconstructed. The presence of a neutral strange hadron in the final state reduces the combinatorial background and extends the measured sensitivity into the low transverse momentum region. The kinematic range is 0 < p_{T}(D^{+}, Lambda_{c}^{+}) < 10 GeV, |eta(D^{+}, Lambda_{c}^{+})| < 1.6, 1.5 < Q^{2} < 1000 GeV^{2} and 0.02 < y < 0.7. Inclusive and differential cross sections for the production of D^{+} mesons are compared to next-to-leading-order QCD predictions. The fraction of c quarks hadronising into Lambda_{c}^{+} baryons is extracted.

The photoproduction of beauty quarks in events with two jets and a muon has been measured with the ZEUS detector at HERA using an integrated luminosity of 110 pb$^{- 1}$. The fraction of jets containing b quarks was extracted from the transverse momentum distribution of the muon relative to the closest jet. Differential cross sections for beauty production as a function of the transverse momentum and pseudorapidity of the muon, of the associated jet and of $x_{\gamma}^{jets}$, the fraction of the photon's momentum participating in the hard process, are compared with MC models and QCD predictions made at next-to-leading order. The latter give a good description of the data.

The photoproduction of $D^{*\pm} (2010)$ mesons associated with a leading neutron has been observed with the ZEUS detector in $ep$ collisions at HERA using an integrated luminosity of 80 pb$^{-1}$. The neutron carries a large fraction, {$x_L>0.2$}, of the incoming proton beam energy and is detected at very small production angles, {$\theta_n<0.8$ mrad}, an indication of peripheral scattering. The $D^*$ meson is centrally produced with pseudorapidity {$|\eta|<1.5$}, and has a transverse momentum {$p_{\it T} > 1.9$ GeV}, which is large compared to the average transverse momentum of the neutron of 0.22 GeV. The ratio of neutron-tagged to inclusive $D^*$ production is $8.85\pm 0.93({\rm stat.})^{+0.48}_{-0.61}({\rm syst.})\%$ in the photon-proton center-of-mass energy range {$130 <W<280$ GeV}. The data suggest that the presence of a hard scale enhances the fraction of events with a leading neutron in the final state.

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)$.