The production of $\Upsilon (nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{-1}$. The $\Upsilon (nS)$ mesons are reconstructed through their decays into two opposite-sign muons. The measurements comprise the differential production cross-sections of the $\Upsilon (1S)$ and $\Upsilon (2S)$ states, their forward-to-backward ratios and nuclear modification factors, performed as a function of the transverse momentum $p_{\mathrm{T}}$ and rapidity in the nucleon-nucleon centre-of-mass frame $y^*$ of the $\Upsilon (nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($-5.0<y^*<-2.5$) for $p$Pb (Pb$p$) collisions. In addition, production cross-sections for $\Upsilon (3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon (nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.

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.

The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < p_T < 4.0 GeV/c) of single electrons as a function of centrality in Au+Au collisions at sqrt(s_NN) = 200 GeV. Contributions to the raw spectra from photon conversions and Dalitz decays of light neutral mesons are measured by introducing a thin (1.7% X_0) converter into the PHENIX acceptance and are statistically removed. The subtracted ``non-photonic'' electron spectra are primarily due to the semi-leptonic decays of hadrons containing heavy quarks (charm and bottom). For all centralities, charm production is found to scale with the nuclear overlap function, T_AA. For minimum-bias collisions the charm cross section per binary collision is N_cc^bar/T_AA = 622 +/- 57 (stat.) +/- 160 (sys.) microbarns.

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 PHENIX experiement has measured the electron-positron pair mass spectrum from 0 to 8 GeV/c^2 in p+p collisions at sqrt(s)=200 GeV. The contributions from light meson decays to e^+e^- pairs have been determined based on measurements of hadron production cross sections by PHENIX. They account for nearly all e^+e^- pairs in the mass region below 1 GeV/c^2. The e^+e^- pair yield remaining after subtracting these contributions is dominated by semileptonic decays of charmed hadrons correlated through flavor conservation. Using the spectral shape predicted by PYTHIA, we estimate the charm production cross section to be 544 +/- 39(stat) +/- 142(syst) +/- 200(model) \mu b, which is consistent with QCD calculations and measurements of single leptons by PHENIX.

J/Psi production in p+p collisions at sqrt(s) = 200 GeV has been Measured in the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over a rapidity range of -2.2 < y < 2.2 and a transverse momentum range of 0 < pT < 9 GeV/c. The statistics available allow a detailed measurement of both the pT and rapidity distributions and are sufficient to constrain production models. The total cross section times branching ratio determined for J/Psi production is B_{ll} sigma_pp^J/psi = 178 +/- 3(stat) +/- 53(syst) +/- 18(norm) nb.

The momentum distribution of electrons from decays of heavy flavor (charm and beauty) for midrapidity |y| < 0.35 in p+p collisions at sqrt(s) = 200 GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over the transverse momentum range 0.3 < p_T < 9 GeV/c. Two independent methods have been used to determine the heavy flavor yields, and the results are in good agreement with each other. A fixed-order-plus-next-to-leading-log pQCD calculation agrees with the data within the theoretical and experimental uncertainties, with the data/theory ratio of 1.72 +/- 0.02^stat +/- 0.19^sys for 0.3 < p_T < 9 GeV/c. The total charm production cross section at this energy has also been deduced to be sigma_(c c^bar) = 567 +/- 57^stat +/- 224^sys micro barns.

The momentum distribution of electrons from semi-leptonic decays of charm and bottom for mid-rapidity |y|<0.35 in p+p collisions at sqrt(s)=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over the transverse momentum range 2 < p_T < 7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D^bar --> e^{+/-} K^{-/+} X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in p_T. A fixed-order-plus-next-to-leading-log (FONLL) perturbative quantum chromodynamics (pQCD) calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is \sigma_{b\b^bar}= 3.2 ^{+1.2}_{-1.1}(stat) ^{+1.4}_{-1.3}(syst) micro b.

At the electron-proton collider HERA the inclusive $D~{*\pm}$ meson photoproduction cross section has been measured with the H1 detector in two different, but partly overlapping, kinematical regions. For the first, where $\langle W_{\gamma p}\rangle \approx 200$\ GeV and $Q~2 < 0.01\,\gev~2$, the result is $\sigma(\gamma p \rightarrow c \bar{c} X) = (13.2 \pm 2.2 ~{+2.1}_{-1.7}\, ~{+9.9}_{-4.8})\,\mu b$. The second measurement for $Q~2 < 4\,\gev~2$ yields $\sigma(\gamma p \rightarrow c \bar{c} X) = ( 9.3 \pm 2.1 ~{+1.9}_{-1.8}\, ~{+6.9}_ {-3.2} )\,\mu b$ at $\langle W_{\gamma p}\rangle \approx 142$\,GeV and $\sigma(\gamma p \rightarrow c \bar{c} X) = ( 20.6 \pm 5.5 ~ {+4.3}_{-3.9}\, ~{+15.4}_{-7.2})\,\mu b$ at $\langle W_{\gamma p} \rangle \approx 230$\,GeV, respectively. The third error accounts for an additional uncertainty due to the proton and photon parton density parametrizations. Differential cross sections are presented as a function of the $D~{*\pm}$ transverse momentum and rapidity. The results compare reasonably well with next-to-leading order QCD calculations. Evidence for diffractive photoproduction of charm quark is presented.