Suppression of the J/psi nuclear modification factor has been seen as a trademark signature of final state effects in large collision systems for decades. In small systems, the nuclear modification was attributed to cold nuclear matter effects until the observation of strong differential suppression of the psi(2S) state in p/d+A collisions suggested the presence of final state effects. Results of J/psi and psi(2S) measurements in the dimuon decay channel are presented here for a collection of different small systems at the same collision energy. These small systems include: p+p, p+Al, p+Au and 3He+Au at collision energy square root s_NN = 200 GeV. The results are predominantly shown in the form of the nuclear modification factor, R_AB, the ratio of the J/psi or psi(2S) invariant yield per nucleon-nucleon collision in collisions of projectile A on target B to that in p+p collisions. Measurements of the J/psi and psi(2S) nuclear modification factor are then compared with shadowing and transport model predictions, as well as to complementary measurements at LHC energies.

ABSTRACT: Mesons composed of heavy quark-antiquark pairs, known as quarkonia, provide the only direct probe of the screening length in the deconfined state of quarks and gluons, known as the quark gluon plasma (QGP), which is believed to be produced in high energy heavy ion collisions. However, the observation of suppression of quarkonia production in heavy ion collisions at high energies is complicated by the modification of quarkonia production in normal nuclear matter. Measuring the modification of quarkonia production due to the effects of normal nuclear matter, often termed cold nuclear matter (CNM) effects, provide a critical baseline for understanding the properties of the QGP. Measurements of CNM effects on quarkonia production are also interesting in their own right, and can be measured independently in proton-nucleus (p+A) collisions. The modification of quarkonia production in p+A collisions provides insight into quarkonia production mechanisms unavailable through the study of proton-proton collisions alone. The study of quarkonia production in p+A collisions over a wide range of kinematic variables can also provide constraints on the modification of parton distribution functions in nuclei. In order to quantify the CNM effects present at the Relativistic Heavy Ion Collider (RHIC), the PHENIX experiment has recorded data on d+Au collisions at \sqsn=200 GeV. The analysis of J/psi→e+e- and psi'→e+e- production from that data set is presented here. Both J/psi and psi' production are found to be suppressed in d+Au relative to p+p collisions, with the suppression increasing for collisions with small impact parameters. The psi' production is found to be much more suppressed than J/psi production, a result which is unexpected based on measurements at lower collision energy and present theoretical pictures. A parametrization of the J/psi modification measured by PHENIX in terms of two CNM effects is also presented. One is the nuclear breakup of the forming quarkonium state through collisions with nucleons during the d+Au collision. The other is the modification of the gluon distribution in the Au nucleus. It is found that the two effects can be separated due to the very different impact parameter dependencies. A strongly non-linear geometric dependence on the modification of the gluon distribution function is observed, with the modification found to be concentrated near the center of the Au nucleus. This parametrization is also used to estimate the modification of J/psi production in Au+Au collisions due to CNM effects. This modification is compared to PHENIX measurements of J/psi production in Au+Au collisions. Suppression of J/psi production in Au+Au collisions beyond CNM effects is observed. This excess suppression is interpreted as suppression of J/psi production due to the formation of a QGP.