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<W_{\gamma p}<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.

Cross sections for e^+p neutral current deep inelastic scattering have been measured at a centre-of-mass energy of sqrt{s}=318 GeV with the ZEUS detector at HERA using an integrated luminosity of 63.2 pb^-1. The double-differential cross section, d^2sigma/dxdQ^2, is presented for 200 GeV^2 < Q^2 < 30000 GeV^2 and for 0.005 < x < 0.65. The single-differential cross-sections dsigma/dQ^2, dsigma/dx and dsigma/dy are presented for Q^2 > 200 GeV^2. The effect of Z-boson exchange is seen in dsigma/dx measured for Q^2 > 10000 GeV^2. The data presented here were combined with ZEUS e^+p neutral current data taken at sqrt{s}=300 GeV and the structure function F_2^{em} was extracted. All results agree well with the predictions of the Standard Model.

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.

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.

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.

Differential cross sections for dijet photoproduction in association with a leading neutron using the reaction e^+ + p --> e^+ + n + jet + jet + X_r have been measured with the ZEUS detector at HERA using an integrated luminosity of 6.4 pb^{-1}. The fraction of dijet events with a leading neutron in the final state was studied as a function of the jet kinematic variables. The cross sections were measured for jet transverse energies E^{jet}_T > 6 GeV, neutron energy E_n > 400 GeV, and neutron production angle theta_n < 0.8 mrad. The data are broadly consistent with factorization of the lepton and hadron vertices and with a simple one-pion-exchange model.

Jet cross sections were measured in charged current deep inelastic e+-p scattering at high boson virtualities Q^2 with the ZEUS detector at HERA II using an integrated luminosity of 0.36 fb^-1. Differential cross sections are presented for inclusive-jet production as functions of Q^2, Bjorken x and the jet transverse energy and pseudorapidity. The dijet invariant mass cross section is also presented. Observation of three- and four-jet events in charged-current e+-p processes is reported for the first time. The predictions of next-to-leading-order (NLO) QCD calculations are compared to the measurements. The measured inclusive-jet cross sections are well described in shape and normalization by the NLO predictions. The data have the potential to constrain the u and d valence quark distributions in the proton if included as input to global fits.

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.

The production and semi-leptonic decay of heavy quarks have been studied in the photoproduction process $e^+p -> e^+ + {dijet} + e^- + X with the ZEUS detector at HERA using an integrated luminosity of 38.5 ${\rm pb^{-1}}$. Events with photon-proton centre-of-mass energies, $W_{\gamma p}$, between 134 and 269 GeV and a photon virtuality, Q^2, less than 1 ${\rm GeV^2}$ were selected requiring at least two jets of transverse energy $E_T^{\rm jet1(2)} >7(6)$ GeV and an electron in the final state. The electrons were identified by employing the ionisation energy loss measurement. The contribution of beauty quarks was determined using the transverse momentum of the electron relative to the axis of the closest jet, $p_T^{\rm rel}$. The data, after background subtraction, were fit with a Monte Carlo simulation including beauty and charm decays. The measured beauty cross section was extrapolated to the parton level with the b quark restricted to the region of transverse momentum $p_T^{b} > p_T^{\rm min} =$ 5 GeV and pseudorapidity $|\eta^{b}| <$ 2. The extrapolated cross section is $1.6 \pm 0.4 (stat.)^{+0.3}_{-0.5} (syst.) ^{+0.2}_{-0.4} (ext.) {nb}$. The result is compared to a perturbative QCD calculation performed to next-to-leading order.