  The nuclear parton distribution functions (nPDFs) of gluons are known to be
difficult to determine with fits of deep inelastic scattering (DIS) and
Drell-Yan (DY) data alone. Therefore, the nCTEQ15 analysis of nuclear PDFs
added inclusive neutral pion production data from RHIC to help in constraining
the gluon. In this analysis, we present a new global analysis of nuclear PDFs
based on a much larger set of single inclusive light hadron data from RHIC and
the LHC. Using our new nCTEQ code (nCTEQ++) with an optimized version of INCNLO
we study systematically the limitations of the theory and the impact of the
fragmentation function uncertainty.
