  We perform a simultaneous global analysis of hadron fragmentation functions
(FFs) to various charged hadrons at next-to-leading order in QCD. The world
data set includes results from electron-positron single-inclusive annihilation,
semi-inclusive deep inelastic scattering, as well as proton-proton collisions
including jet fragmentation measurements which lead to strong constraints on
the gluon fragmentations. By carefully selecting hadron kinematics to ensure
the validity of QCD factorization and the convergence of perturbative
calculations, we achieve a satisfying best fit with $\chi^2/$d.o.f.$=0.90$, in
the simultaneous extraction of FFs for light charged hadrons ($\pi^{\pm}$,
$K^{\pm}$ and $p/\bar{p}$). The total momentum of $u$, $d$ quarks and gluon
carried by light charged hadrons have been determined precisely. That urges
future precision measurements on fragmentation to neutral hadrons, which are
crucial for the test of fundamental sum rules in QCD fragmentation.
