We perform a simultaneous global analysis of hadron fragmentation functions (FFs) to various charged hadrons (π^{±}, K^{±}, and p/p[over ¯]) at next-to-leading order in QCD. The world data include results from electron-positron single-inclusive annihilation, semi-inclusive deep inelastic scattering, as well as proton-proton collisions including jet fragmentation measurements for the first time, 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 χ^{2}/d.o.f.=0.90. The total momentum of u, d quarks and gluon carried by light charged hadrons have been determined precisely, urging precision determinations of FFs to neutral hadrons for a test of fundamental sum rules in QCD fragmentation.