Yields for J/psi production in Cu+Cu collisions at sqrt (s_NN)= 200 GeV have been measured by the PHENIX experiment over the rapidity range |y| < 2.2 at transverse momenta from 0 to beyond 5 GeV/c. The invariant yield is obtained as a function of rapidity, transverse momentum and collision centrality, and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data provide greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, providing a key constraint that is needed for disentangling cold and hot nuclear matter effects.

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 at the Relativistic Heavy Ion Collider has measured low mass vector meson, $\omega$, $\rho$, and $\phi$, production through the dimuon decay channel at forward rapidity ($1.2<|y|<2.2$) in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The differential cross sections for these mesons are measured as a function of both $p_T$ and rapidity. We also report the integrated differential cross sections over $1<p_T<7$ GeV/$c$ and $1.2<|y|<2.2$: $d\sigma/dy(\omega+\rho\rightarrow\mu\mu) = 80 \pm 6 \mbox{(stat)} \pm 12 \mbox{(syst)}$ nb and $d\sigma/dy(\phi\rightarrow\mu\mu) = 27 \pm 3 \mbox{(stat)} \pm 4 \mbox{(syst)}$ nb. These results are compared with midrapidity measurements and calculations.

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 PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured J/psi production for rapidities 2.2 < y < 2.2 in Au+Au collisions at sqrt(s_NN) = 200 GeV. The J/psi invariant yield and nuclear modification factor R_AA as a function of centrality, transverse momentum and rapidity are reported. A suppression of J/psi relative to binary collision scaling of proton-proton reaction yields is observed. Models which describe the lower energy J/Psi data at the Super Proton Synchrotron (SPS) invoking only J/psi destruction based on the local medium density would predict a significantly larger suppression at RHIC and more suppression at mid rapidity than at forward rapidity. Both trends are contradicted by our data.

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.

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.

We report on the inclusive production cross sections of J/$\psi$, $\psi$(2S), $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S), measured at forward rapidity with the ALICE detector in pp collisions at a center-of-mass energy $\sqrt{s}=8$ TeV. The analysis is based on data collected at the LHC and corresponds to an integrated luminosity of 1.28 pb$^{-1}$. Quarkonia are reconstructed in the dimuon-decay channel. The differential production cross sections are measured as a function of the transverse momentum $p_{\rm T}$ and rapidity $y$, over the $p_{\rm T}$ ranges $0<p_{\rm T}<20$ GeV/$c$ for J/$\psi$, $0<p_{\rm T}<12$ GeV/$c$ for all other resonances, and for $2.5<y<4$. The cross sections, integrated over $p_{\rm T}$ and $y$, and assuming unpolarized quarkonia, are $\sigma_{{\rm J}/\psi} = 8.98\pm0.04\pm0.82$ $\mu$b, $\sigma_{\psi{\rm (2S)}} = 1.23\pm0.08\pm0.22$ $\mu$b, $\sigma_{\Upsilon{\rm(1S)}} = 71\pm6\pm7$ nb, $\sigma_{\Upsilon{\rm(2S)}} = 26\pm5\pm4$ nb and $\sigma_{\Upsilon{\rm(3S)}} = 9\pm4\pm1$ nb, where the first uncertainty is statistical and the second one is systematic. These values agree, within at most $1.4\sigma$, with measurements performed by the LHCb collaboration in the same rapidity range.

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.

The inclusive production cross sections at forward rapidity of J/$\psi$, $\psi$(2S), $\Upsilon$(1S) and $\Upsilon$(2S) are measured in pp collisions at $\sqrt{s} = 7$ TeV with the ALICE detector at the LHC. The analysis is based in a data sample corresponding to an integrated luminosity of 1.35 pb$^{-1}$. Quarkonia are reconstructed in the dimuon-decay channel and the signal yields are evaluated by fitting the $\mu^+\mu^-$ invariant mass distributions. The differential production cross sections are measured as a function of the transverse momentum $p_{\rm T}$ and rapidity $y$, over the ranges $0 < p_{\rm T} < 20$ GeV/$c$ for J/$\psi$, $0 < p_{\rm T} < 12$ GeV/$c$ for all other resonances and for $2.5 < y < 4$. The measured cross sections integrated over $p_{\rm T}$ and $y$, and assuming unpolarized quarkonia, are: $\sigma_{J/\psi} = 6.69 \pm 0.04 \pm 0.63$ $\mu$b, $\sigma_{\psi^{\prime}} = 1.13 \pm 0.07 \pm 0.14$ $\mu$b, $\sigma_{\Upsilon{\rm(1S)}} = 54.2 \pm 5.0 \pm 6.7$ nb and $\sigma_{\Upsilon{\rm (2S)}} = 18.4 \pm 3.7 \pm 2.2$ nb, where the first uncertainty is statistical and the second one is systematic. The results are compared to measurements performed by other LHC experiments and to theoretical models.