a serious concern for neutrino oscillation experiments. This paper shows the effect of fitting current CCQE and multinucleon models to the [MiniBooNE]{} [@mb-ccqe-2010; @mb-ccqe-antinu-2013] and [MINER$\nu$A]{} [@minerva-nu-ccqe; @minerva-antinu-ccqe] datasets to a variety of models implemented in [<span style="font-variant:small-caps;">NEUT</span>]{}by members of T2K’s Neutrino Interactions Working Group (NIWG). Previous constraints on the CCQE model produced by the NIWG and used in T2K oscillation analyses only considered an RFG model, and recommended the [<span style="font-variant:small-caps;">NEUT</span>]{}default central value for the axial mass [$M_{\textrm{A}}$]{}= 1.21 GeV$/c^{2}$ based on the value found by the K2K experiment [@k2kCCQE], with an error large enough to cover fits to the [MiniBooNE]{}neutrino mode CCQE dataset [@mb-ccqe-2010], as is fully described in Reference [@Abe:2015awa]. This work improves on the previous situation by including more sophisticated effects proposed to explain the large axial mass anomaly, and by using all of the newly available CCQE data to constrain all model parameters without reference to the default [<span style="font-variant:small-caps;">NEUT</span>]{}model.

The models which have been implemented in the [<span style="font-variant:small-caps;">NEUT</span>]{}generator are discussed in Section \[sec:neut\] and Section \[sec:nuwro\] discusses cross-generator validation. Section \[sec:extdata\] gives a brief overview of the [MiniBooNE]{}and [MINER$\nu$A]{}data used in the fit. The nominal [<span style="font-variant:small-caps;">NEUT</span>]{}predictions for these datasets are shown in 