We present our results on neutral pion (π0) 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.