2026-07-072026-07-07http://salesiana.dossiersoluciones.com/handle/123456789/200518In this talk, we present the results from our recent global reanalysis of nuclear parton distribution functions (nPDFs), where the DGLAP-evolving nPDFs are constrained by nuclear hard process data from deep inelastic $l+A$ scattering (DIS) and the Drell-Yan (DY) process in $p+A$ collisions, and by sum rules. The main improvements over our earlier work {\em EKS98} are the automated $χ^2$ minimization, better controllable fit functions and possibility for error estimates. The obtained 16-parameter fit to N=514 datapoints is good, $χ^2/{\rm d.o.f}=0.82$. Fit quality comparison and the error estimates obtained show that the old {\em EKS98} parametrization is fully consistent with the present automated reanalysis. Comparison with other global nPDF analyses is presented as well. Within the DGLAP framework we also discuss the possibility of incorporating a clearly stronger gluon shadowing, which is suggested by the RHIC BRAHMS data from d+Au collisions.10 pages, talk given at the workshop "High-pT physics at LHC", March 23-27, 2007, Department of Physics, University of Jyvaskyla, FinlandHigh Energy Physics - Phenomenologytext