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.

J/psi production has been measured in proton-proton collisions at sqrt(s)= 200 GeV over a wide rapidity and transverse momentum range by the PHENIX experiment at RHIC. Distributions of the rapidity and transverse momentum, along with measurements of the mean transverse momentum and total production cross section are presented and compared to available theoretical calculations. The total J/psi cross section is 3.99 +/- 0.61(stat) +/- 0.58(sys) +/- 0.40(abs) micro barns. The mean transverse momentum is 1.80 +/- 0.23(stat) +/- 0.16(sys) GeV/c.

Muon production at forward rapidity (1.5 < |\eta| < 1.8) has been measured by the PHENIX experiment over the transverse momentum range 1 < p_T \le 3 GeV/c in sqrt(s) = 200 GeV p+p collisions at the Relativistic Heavy Ion Collider. After statistically subtracting contributions from light hadron decays an excess remains which is attributed to the semileptonic decays of hadrons carrying heavy flavor, i.e. charm quarks or, at high p_T, bottom quarks. The resulting muon spectrum from heavy flavor decays is compared to PYTHIA and a next-to-leading order perturbative QCD calculation. PYTHIA is used to determine the charm quark spectrum that would produce the observed muon excess. The corresponding differential cross section for charm quark production at forward rapidity is determined to be d\sigmac c^bar)/dy|_(y=1.6)=0.243 +/- 0.013 (stat.) +/- 0.105 (data syst.) ^(+0.049(-0.087) (PYTHIA syst.) mb.

J/psi production in d+Au and p+p collisions at sqrt(s_NN) = 200 GeV has been measured by the PHENIX experiment at rapidities -2.2 < y < +2.4. The cross sections and nuclear dependence of J/\psi production versus rapidity, transverse momentum, and centrality are obtained and compared to lower energy p+A results and to theoretical models. The observed nuclear dependence in d+Au collisions is found to be modest, suggesting that the absorption in the final state is weak and the shadowing of the gluon distributions is small and consistent with Dokshitzer-Gribov-Lipatov-Altarelli-Parisi-based parameterizations that fit deep-inelastic scattering and Drell-Yan data at lower energies.

The Upsilon production cross section in proton-proton collisions at sqrt(s) = 7 TeV is measured using a data sample collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 3.1 +/- 0.3 inverse picobarns. Integrated over the rapidity range |y|<2, we find the product of the Upsilon(1S) production cross section and branching fraction to dimuons to be sigma(pp to Upsilon(1S) X) B(Upsilon(1S) to mu+ mu-) = 7.37 +/- 0.13^{+0.61}_{-0.42}\pm 0.81 nb, where the first uncertainty is statistical, the second is systematic, and the third is associated with the estimation of the integrated luminosity of the data sample. This cross section is obtained assuming unpolarized Upsilon(1S) production. If the Upsilon(1S) production polarization is fully transverse or fully longitudinal the cross section changes by about 20%. We also report the measurement of the Upsilon(1S), Upsilon(2S), and Upsilon(3S) differential cross sections as a function of transverse momentum and rapidity.

The first observation of $Z$ boson production in proton-lead collisions at a centre-of-mass energy per proton-nucleon pair of $\sqrt{s_{NN}}=5~\text{TeV}$ is presented. The data sample corresponds to an integrated luminosity of $1.6~\text{nb}^{-1}$ collected with the LHCb detector. The $Z$ candidates are reconstructed from pairs of oppositely charged muons with pseudorapidities between 2.0 and 4.5 and transverse momenta above $20~\text{GeV}/c$. The invariant dimuon mass is restricted to the range $60-120~\text{GeV}/c^2$. The $Z$ production cross-section is measured to be \begin{eqnarray*} \sigma_{Z\to\mu^+\mu^-}(\text{fwd})&=&13.5^{+5.4}_{-4.0}\text{(stat.)}\pm1.2\text{(syst.)}~\text{nb} \end{eqnarray*} in the direction of the proton beam and \begin{eqnarray*} \sigma_{Z\to\mu^+\mu^-}(\text{bwd}) & =&10.7^{+8.4}_{-5.1}\text{(stat.)}\pm1.0\text{(syst.)}~\text{nb} \end{eqnarray*} in the direction of the lead beam, where the first uncertainty is statistical and the second systematic.

We have studied J/psi production in pp collisions at $\sqrt{s}=7$ TeV at the LHC through its muon pair decay. The polar and azimuthal angle distributions of the decay muons were measured, and results on the J/$\psi$ polarization parameters $\lambda_{\theta}$ and $\lambda_\phi$ were obtained. The study was performed in the kinematic region 2.5<y<4, 2<$p_{\rm T}$<8 GeV/$c$, in the helicity and Collins-Soper reference frames. In both frames, the polarization parameters are compatible with zero, within uncertainties.

A measurement of the J/psi and psi(2S) production cross sections in pp collisions at sqrt(s)=7 TeV with the CMS experiment at the LHC is presented. The data sample corresponds to an integrated luminosity of 37 inverse picobarns. Using a fit to the invariant mass and decay length distributions, production cross sections have been measured separately for prompt and non-prompt charmonium states, as a function of the meson transverse momentum in several rapidity ranges. In addition, cross sections restricted to the acceptance of the CMS detector are given, which are not affected by the polarization of the charmonium states. The ratio of the differential production cross sections of the two states, where systematic uncertainties largely cancel, is also determined. The branching fraction of the inclusive B to psi(2S) X decay is extracted from the ratio of the non-prompt cross sections to be: BR(B to psi(2S) X) = (3.08 +/- 0.12(stat.+syst.) +/- 0.13(theor.) +/- 0.42(BR[PDG])) 10^-3

Using 1.8 fb-1 of pp collisions at a center-of-mass energy of 7 TeV recorded by the ATLAS detector at the Large Hadron Collider, we present measurements of the production cross sections of Upsilon(1S,2S,3S) mesons. Upsilon mesons are reconstructed using the di-muon decay mode. Total production cross sections for p_T<70 GeV and in the rapidity interval |Upsilon|<2.25 are measured to be 8.01+-0.02+-0.36+-0.31 nb, 2.05+-0.01+-0.12+-0.08 nb, 0.92+-0.01+-0.07+-0.04 nb respectively, with uncertainties separated into statistical, systematic, and luminosity measurement effects. In addition, differential cross section times di-muon branching fractions for Upsilon(1S), Upsilon(2S), and Upsilon(3S) as a function of Upsilon transverse momentum p_T and rapidity are presented. These cross sections are obtained assuming unpolarized production. If the production polarization is fully transverse or longitudinal with no azimuthal dependence in the helicity frame the cross section may vary by approximately +-20%. If a non-trivial azimuthal dependence is considered, integrated cross sections may be significantly enhanced by a factor of two or more. We compare our results to several theoretical models of Upsilon meson production, finding that none provide an accurate description of our data over the full range of Upsilon transverse momenta accessible with this dataset.

The production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons decaying into the dimuon final state is studied with the LHCb detector using a data sample corresponding to an integrated luminosity of 3.3 pb^{-1} collected in proton-proton collisions at a centre-of-mass energy of sqrt{s}=2.76 TeV. The differential production cross-sections times dimuon branching fractions are measured as functions of the Upsilon transverse momentum and rapidity, over the ranges p_T<15 GeV/c and 2.0<y<4.5. The total cross-sections in this kinematic region, assuming unpolarised production, are measured to be sigma(pp -> Upsilon(1S) X) x B(Upsilon(1S) -> mu+mu-) = 1.111 +/- 0.043 +/- 0.044 nb, sigma(pp -> Upsilon(2S) X) x B(Upsilon(2S) -> mu+mu-) = 0.264 +/- 0.023 +/- 0.011 nb, sigma(pp -> Upsilon(3S) X) x B(Upsilon(3S) -> mu+mu-) = 0.159 +/- 0.020 +/- 0.007 nb, where the first uncertainty is statistical and the second systematic.