the initial state nucleon distribution. As a result, changing ${\ensuremath{p_{\mathrm{F}}}\xspace}^{\mbox{\scriptsize{RFG}}}$ affects a wide range of [$Q^{2}$]{}, whereas changing ${\ensuremath{p_{\mathrm{F}}}\xspace}^{\mbox{\scriptsize{SF}}}$ only affects very low [$Q^{2}$]{}events.

The [multinucleon–neutrino]{}(2p2h) model from Nieves [*et al.*]{} [@nieves; @nievesExtension] has been implemented in [<span style="font-variant:small-caps;">NEUT</span>]{}as described in Reference [@peterthesis]. The cross section as a function of neutrino energy and the outgoing lepton kinematics was made available by the authors of the model and is implemented as a series of lookup tables for various nuclear targets and neutrino species. The tables provided had hadronic variables integrated out, so a generic model based on Reference [@multinucleon] for simulating the initial and final hadronic states was used for generating [<span style="font-variant:small-caps;">NEUT</span>]{}events[^1]. The discrepancy between the leptonic and hadronic simulation makes the current [<span style="font-variant:small-caps;">NEUT</span>]{}implementation of the Nieves model inadequate for comparisons with experimental measurements of the final state hadrons from [CCQE]{}events (such as can be found in Reference [@minerva_ccqe_p_2014]). For this reason, only leptonic measurements are used in this analysis. As the Nieves model is very complex, the current [<span style="font-variant:small-caps;">NEUT</span>]{}implementation does not allow fundamental model parameters to be changed. For simplicity, only a simple scaling parameter which changes the normalization of 2p2h events has been considered in this 