The genie neutrino Monte Carlo describes neutrino-induced hadronization with an effective model, known as Andreopoulos-Gallagher-Kehayias-Yang (agky), which is interfaced with pythia at high invariant mass. Only the low-mass agky model parameters were extracted from hadronic shower data from the FNAL 15 ft and BEBC experiments. In this paper, the first hadronization tune on averaged charged multiplicity data from deuterium and hydrogen bubble chamber experiments is presented, with a complete estimation of parameter uncertainties. A partial tune on deuterium data highlights the tensions between hydrogen and deuterium datasets
