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.

We report measurements of single- and double- spin asymmetries for $W^{\pm}$ and $Z/\gamma^*$ boson production in longitudinally polarized $p+p$ collisions at $\sqrt{s} = 510$ GeV by the STAR experiment at RHIC. The asymmetries for $W^{\pm}$ were measured as a function of the decay lepton pseudorapidity, which provides a theoretically clean probe of the proton's polarized quark distributions at the scale of the $W$ mass. The results are compared to theoretical predictions, constrained by recent polarized deep inelastic scattering measurements, and show a preference for a sizable, positive up antiquark polarization in the range $0.05<x<0.2$.

The invariant differential cross section for inclusive neutral pion production in p+p collisions at sqrt(s_NN) = 200 GeV has been measured at mid-rapidity |eta| < 0.35 over the range 1 < p_T <~ 14 GeV/c by the PHENIX experiment at RHIC. Predictions of next-to-leading order perturbative QCD calculations are consistent with these measurements. The precision of our result is sufficient to differentiate between prevailing gluon-to-pion fragmentation functions.

We report a measurement of the longitudinal double-spin asymmetry A_LL and the differential cross section for inclusive Pi0 production at midrapidity in polarized proton collisions at sqrt(s) = 200 GeV. The cross section was measured over a transverse momentum range of 1 < p_T < 17 GeV/c and found to be in good agreement with a next-to-leading order perturbative QCD calculation. The longitudinal double-spin asymmetry was measured in the range of 3.7 < p_T < 11 GeV/c and excludes a maximal positive gluon polarization in the proton. The mean transverse momentum fraction of Pi0's in their parent jets was found to be around 0.7 for electromagnetically triggered events.

We report a new high-precision measurement of the mid-rapidity inclusive jet longitudinal double-spin asymmetry, $A_{LL}$, in polarized $pp$ collisions at center-of-mass energy $\sqrt{s}=200$ GeV. The STAR data place stringent constraints on polarized parton distribution functions extracted at next-to-leading order from global analyses of inclusive deep inelastic scattering (DIS), semi-inclusive DIS, and RHIC $pp$ data. The measured asymmetries provide evidence for positive gluon polarization in the Bjorken-$x$ region $x>0.05$.

The $\jpsi$ $\pt$ spectrum and nuclear modification factor ($\raa$) are reported for $\pt < 5 \ \gevc$ and $|y|<1$ from 0\% to 60\% central Au+Au and Cu+Cu collisions at $\snn = 200 \ \gev$ at STAR. A significant suppression of $\pt$-integrated $\jpsi$ production is observed in central Au+Au events. The Cu+Cu data are consistent with no suppression, although the precision is limited by the available statistics. $\raa$ in Au+Au collisions exhibits a strong suppression at low transverse momentum and gradually increases with $\pt$. The data are compared to high-$\pt$ STAR results and previously published BNL Relativistic Heavy Ion Collider results. Comparing with model calculations, it is found that the invariant yields at low $\pt$ are significantly above hydrodynamic flow predictions but are consistent with models that include color screening and regeneration.

We report a measurement of high-p_T inclusive pi^0, eta, and direct photon production in p+p and d+Au collisions at sqrt(s_NN) = 200 GeV at midrapidity (0 < eta < 1). Photons from the decay pi^0 -> gamma gamma were detected in the Barrel Electromagnetic Calorimeter of the STAR experiment at the Relativistic Heavy Ion Collider. The eta -> gamma gamma decay was also observed and constituted the first eta measurement by STAR. The first direct photon cross section measurement by STAR is also presented, the signal was extracted statistically by subtracting the pi^0, eta, and omega(782) decay background from the inclusive photon distribution observed in the calorimeter. The analysis is described in detail, and the results are found to be in good agreement with earlier measurements and with next-to-leading order perturbative QCD calculations.

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.

Transverse momentum distributions and yields for $\pi^{\pm}$, $K^{\pm}$, $p$ and $\bar{p}$ in $p+p$ collisions at $\sqrt{s}$=200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter $T_{\rm inv}$, mean transverse momentum $<p_T>$ and yield per unit rapidity $dN/dy$ at each energy, and compare them to other measurements at different $\sqrt{s}$ in $p+p$ and $p+\bar{p}$ collisions. We also present the scaling properties such as $m_T$ scaling, $x_T$ scaling on the $p_T$ spectra between different energies. To discuss the mechanism of the particle production in $p+p$ collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.