We introduce the pQCD framework based on collinear factorization for inclusive hadron production in pp and pA collisions and review the current status of theoretical calculations, their uncertainties, and how they compare to data. Particular emphasis is put on what we currently know about parton-to-hadron fragmentation functions, their range of applicability, and how they are determined from global QCD analyses.