  Next generation neutrino oscillation experiments utilize details of hadronic
final states to improve the precision of neutrino interaction measurements. The
hadronic system was often neglected or poorly modeled in the past, but they
have significant effects on high precision neutrino oscillation and
cross-section measurements. Among the physics of hadronic systems in neutrino
interactions, the hadronization model controls multiplicities and kinematics of
final state hadrons from the primary interaction vertex. For relatively high
invariant mass events, many neutrino experiments rely on the PYTHIA program.
Here, we show a possible improvement of this process in neutrino event
generators, by utilizing expertise from the HERMES experiment. Finally, we
estimate the impact on the systematics of hadronization models for neutrino
mass hierarchy analysis using atmospheric neutrinos such as the PINGU
experiment.
