Section \[sec:mcpred\] for a variety of models. Section \[sec:fitprod\] discusses the fit procedure. The results of fake data studies and the fit to external data are given in Section \[sec:results\]. In Section \[sec:discussion\] we interpret the results and discuss possible implications in cross section and neutrino oscillation analyses and Section \[sec:summary\] summarizes the results.

Interaction models in [<span style="font-variant:small-caps;">NEUT</span>]{} {#sec:neut}
============================================================================

The motivation for, and an overview of, new CCQE models has already been discussed. This section will briefly outline the important technical details of the models as implemented in [<span style="font-variant:small-caps;">NEUT</span>]{}, and highlight any caveats that should be borne in mind when fitting with them. The models used in the fits include the SF model, [multinucleon–neutrino]{}interactions, and RPA.

![The probability distribution for initial state protons within an oxygen nucleus for Benhar’s SF model [@sf] as a function of the removal energy ($E_{\mathrm{R}}$) and the magnitude of the nucleon momentum ($|\mathbf{p}|$). The SF is normalized such that the integral of this distribution is 1.[]{data-label="fig:SF_2D_distribution"}](figures/SF_oxygen_2D_new.pdf){width="0.9\columnwidth"}

The [<span style="font-variant:small-caps;">NEUT</span>]{}implementation of the SF model from Omar Benhar and collaborators [@sf] is described fully in Reference [@andy_thesis]. Although SF is a generic term, in this work it will specifically refer to the Benhar SF. 