In a Glauber Monte Carlo framework, taking account of the transverse spatial distribution of hard partons in the nucleon, we analyse the nuclear modification factor R dAu for J/ψ in d+Au collisions with the EPS09 shadowing parametrization. After the influence of nucleon hard partons distribution is considered, a clearly upward correction is revealed for the dependence of R dAu on N coll in peripheral d+Au collisions, however, an unconspicuous correction is shown for the results versus p T . The theoretical results are in good agreement with the experimental data from PHENIX. 24.85.+p, 12.38.Bx For J/ψ production is sensitive to both cold nuclear matter (CNM) and Quark-Gluon Plasma (hot-dense matter) [1] , the "normal suppression"in the CNM should be subtracted in any interpretation of charmonium production in heavy ion collisions. The Glauber Monte Carlo (GMC) approach [2−6] , which can simulate experimentally observable quantities and analyse real data, is an ideal tool to study the CNM effects provided in deuteron-gold (d+Au) collisions [7] . In the framework of the GMC, only the essential J/ψ production process [8] , g + g → cc, is considered. Thus, the subprocess cross section can be obtained by extracting from the proton-proton (p+p) experimental data [9] and given as the Monte Carlo inputs of rapidity (y) and transverse momentum (p T ) [3,4] . The shadowing and "normal nuclear absorption"effects, which are two dominant CNM effects at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) energies, can also be considered properly in the GMC framework. Since shadowing should depend on the spatial position of the interaction parton within the nucleus, the inhomogeneous shadowing effect [10,11] should also be taken into account. In this paper, the latest shadowing parametrization EPS09 (Eskola, Paukkunen and Salgado) [12] is used.
PACSAt ultra-high energy domain, the mechanism of inelastic hadronic collisions is dominated by the contribution from small-x gluons and the influence of the transverse spatial distribution of hard partons in the nucleon become important [13,14] . Unfortunately, comparing with the longitudinal momentum distribution of partons in the nucleon, the measurements of the transverse spatial distribution of partons are rather limited [15] . In this paper, two kinds of hard partons distributions in the nucleon are used. One is assumed that the hard partons are uniformly distributed in a hard sphere nucleon, the other is derived from fits to J/ψ photo-production data at HERA and FNAL [13,15] . The corresponding results with both of them will be given in the following part of the paper.