Multiplicity dependent J=ψ and ψ(2S) production at forward and backward rapidity in p + p collisions at (Formula Presented) GeV

Recent measurements of J=ψ production as a function of event charged-particle multiplicity at the collision energies of both the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) show enhanced J=ψ production yields with increasing multiplicity. One potential explanation for this type of dependence is multiparton interactions (MPI). We present the first study of potential autocorrelations at RHIC energies and forward and backward rapidity of self-normalized J=ψ yields and ψ(2S) to J=ψ ratio, as a function of self-normalized multiplicity in p + p collisions. In addition, detailed PYTHIA studies tuned to RHIC energies were performed to investigate the MPI impacts. We find that the PHENIX data at RHIC are consistent with recent LHC measurements and can only be described by PYTHIA calculations that include MPI effects. The forward and backward ψ(2S) to J=ψ ratio is found to be less dependent on the charged-particle multiplicity.