  We present our results on neutral pion pi0 production in pp and dAu
collisions at RHIC energy. Pion spectra are calculated in a next-to-leading
order (NLO) perturbative QCD-based model. The model includes the transverse
component of the initial parton distribution (``intrinsic kT''). We compare our
results to the available experimental data from RHIC, and fit the data with
high precision. The calculation tuned this way is repeated for the dAu
collision, and used to investigate the interplay of shadowing and multiple
scattering at RHIC. The centrality dependence of the nuclear modification
factor shows a measurable difference between different shadowing
parameterizations.
