With a data sample containing 1.1×105 J/ψ→μ+μ− decays reconstructed with 16 MeV/c2 rms mass resolution, we have measured the differential cross sections versus Feynman-x, rapidity, and pT for the production of J/ψ and ψ’ in 800 GeV/c p-Au collisions. Our results are compared with leading-order QCD predictions and with previous measurements. While the shapes of the cross sections are in qualitative agreement with QCD predictions, the magnitudes disagree by factors of 7 (J/ψ) and 25 (ψ’). Assuming an appropriate form for the differential cross sections in regions not measured we derive a total J/ψ production cross section σ(p+N→J/ψ+X)=442±2±88 nb/nucleon and a (model-dependent) total ψ’ cross secton σ(p+N→ψ’+X)=75±5±22 nb/nucleon. For J/ψ produced at central rapidity, dσ(p+N→J/ψ+X)/dy‖y=0=230±5±46 nb/nucleon.
The differential cross sections dσ/dxF for J/ψ produced inclusively in 800 GeV/c p-Cu and p-Be collisions have been measured in the kinematic range 0.30≤xF≤0.95 through the decay mode J/ψ→μ+μ−. They are compared with the predictions of the semilocal duality model for several sets of parton density functions. No evidence for a suggested intrinsic charm contribution to the cross section is observed. The ratio of the differential cross sections for Cu and Be targets confirms the suppression of J/ψ production in heavy nuclei at large xF.