with [<span style="font-variant:small-caps;">NEUT</span>]{}for each model using the default parameters given in Tables \[tab:RFG\_params\] and \[tab:SF\_params\] and the published flux for each dataset. The flux averaged cross section predictions were produced using the following method:

1.  For each event apply experiment-specific cuts and, if the event passes, calculate the relevant reconstructed quantity and fill the 1D or 2D event rate histogram.

2.  Calculate the event rate by integrating the MC event rate histogram (flux $\times$ cross section).

3.  Integrate the published flux histogram to get the average flux.

4.  Scale the filled histogram by the event rate divided by the average flux to get the flux averaged cross section per target nucleon.

5.  Divide the content of each bin by the bin width.

The default predictions for a variety of models available in [<span style="font-variant:small-caps;">NEUT</span>]{}, as well as the data, are shown in Figures \[fig:MIN\_20deg\_nominal\], \[fig:MB\_1D\_nominal\] and \[fig:MB\_2D\_nominal\] for the [MINER$\nu$A]{}, [MiniBooNE]{}single-differential and [MiniBooNE]{}double-differential samples, respectively. The Nieves 2p2h contribution is also shown on these plots for reference.

To produce a meaningful nominal $\chi^{2}$ for the [MiniBooNE]{}datasets, it is necessary to fit the [MiniBooNE]{}normalization parameters. The single and double-differential plots shown in Fig. \[fig:MB\_1D\_nominal\] and \[fig:MB\_2D\_nominal\] are scaled according to the 