Final state radiation correction used in the $\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-p_{T}^{Z}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Correlation matrix of statistical uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Measured total $Z$-boson cross-section for different datasets. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $y^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Systematic uncertainties in the single differential cross-sections in interval regions of $y^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

NuE background+LEE fractional covariance matrix after the 1mu1p constraint from arXiv:2110.14080

NuE background fractional covariance matrix before the 1mu1p constraint from arXiv:2110.14080

NuE background+LEE fractional covariance matrix before the 1mu1p constraint from arXiv:2110.14080

NuE simulation from arXiv:2110.14080

Covariance matrix of the $\nu_\mu$CC inclusive total cross section per nucleon in neutrino energy bins.

Covariance matrix of the $\nu_\mu$CC inclusive differential cross section per nucleon in muon energy bins.

Covariance matrix of the $\nu_\mu$CC inclusive differential cross section per nucleon in hadronic energy transfer bins.

Additional smearing matrix of the $\nu_\mu$CC inclusive total cross section per nucleon in neutrino energy bins.

Additional smearing matrix of the $\nu_\mu$CC inclusive differential cross section per nucleon in muon energy bins.

Additional smearing matrix of the $\nu_\mu$CC inclusive total cross section per nucleon in neutrino energy bins.

Partially contained $\nu_e$CC signal efficiency as a function of true neutrino energy. Each bin shows the fraction of Monte-Carlo $\nu_e$ CC events with true neutrino interaction vertex in the fiducial volume (3 cm inside the TPC active volume) which are selected in this channel.

Fully contained $\nu_e$CC true neutrino energy vs reconstructed neutrino energy. Each bin shows the relative number of selected fully contained Monte-Carlo $\nu_e$ CC events with true neutrino interaction vertex in the fiducial volume (3 cm inside the TPC active volume) which have the corresponding true neutrino energy and reconstructed neutrino energy values. Each axis has 60 bins covering an energy range from 0 to 3 GeV, corresponding to 0.05 GeV per bin.

Partially contained $\nu_e$CC true neutrino energy vs reconstructed neutrino energy. Each bin shows the relative number of selected partially contained Monte-Carlo $\nu_e$ CC events with true neutrino interaction vertex in the fiducial volume (3 cm inside the TPC active volume) which have the corresponding true neutrino energy and reconstructed neutrino energy values. Each axis has 60 bins covering an energy range from 0 to 3 GeV, corresponding to 0.05 GeV per bin.

Fully contained $\nu_\mu$CC true neutrino energy vs reconstructed neutrino energy. Each bin shows the relative number of selected fully contained Monte-Carlo $\nu_\mu$ CC events with true neutrino interaction vertex in the fiducial volume (3 cm inside the TPC active volume) which have the corresponding true neutrino energy and reconstructed neutrino energy values. Each axis has 60 bins covering an energy range from 0 to 3 GeV, corresponding to 0.05 GeV per bin.

Partially contained $\nu_\mu$CC true neutrino energy vs reconstructed neutrino energy. Each bin shows the relative number of selected partially contained Monte-Carlo $\nu_\mu$ CC events with true neutrino interaction vertex in the fiducial volume (3 cm inside the TPC active volume) which have the corresponding true neutrino energy and reconstructed neutrino energy values. Each axis has 60 bins covering an energy range from 0 to 3 GeV, corresponding to 0.05 GeV per bin.

Fully contained $\nu_e$CC data, signal, background, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The background includes neutral current events, $\nu_\mu$CC events, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining intrinsic $\nu_e$CC events. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Partially contained $\nu_e$CC data, signal, background, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The background includes neutral current events, $\nu_\mu$CC events, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining intrinsic $\nu_e$CC events. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Fully contained $\nu_\mu$CC data, signal, and background predictions. Events in the $\nu_e$CC or $\pi^0$ selections have been removed. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The background includes neutral current events, $\nu_e$CC events, $\nu_\mu$CC events with greater than or equal to one $\pi^0$, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining $\nu_\mu$CC events.

Partially contained $\nu_\mu$CC data, signal, and background predictions. Events in the $\nu_e$CC or $\pi^0$ selections have been removed. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The background includes neutral current events, $\nu_e$CC events, $\nu_\mu$CC events with greater than or equal to one $\pi^0$, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining $\nu_\mu$CC events.

Fully contained $\nu_\mu$CC$\pi^0$ data, signal, background predictions. Events in the $\nu_e$CC selection have been removed. Note that the rightmost bin is an overflow bin, containing all events with reconstructed $\pi^0$ kinetic energy greater than 1 GeV. The background includes neutral current events, $\nu_e$CC events, $\nu_\mu$CC events without a $\pi^0$, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining $\nu_\mu$CC$\pi^0$ events.

Partially contained $\nu_\mu$CC$\pi^0$ data, signal, and background predictions. Events in the $\nu_e$CC selection have been removed. Note that the rightmost bin is an overflow bin, containing all events with reconstructed $\pi^0$ kinetic energy greater than 1 GeV. The background includes neutral current events, $\nu_e$CC events, $\nu_\mu$CC events without a $\pi^0$, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining $\nu_\mu$CC$\pi^0$ events.

NC$\pi^0$ data, signal, and background predictions. Events in the $\nu_e$CC selection have been removed. Note that the rightmost bin is an overflow bin, containing all events with reconstructed $\pi^0$ kinetic energy greater than 1 GeV. The background includes $\nu_e$CC events, $\nu_\mu$CC events, neutral current events without a $\pi^0$, events with a true neutrino interaction vertex outside the fiducial volume (3 cm inside the TPC active volume), and cosmic ray backgrounds. The signal includes the remaining NC$\pi^0$ events.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.0), the standard prediction with no low energy excess. No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.1). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.2). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.3). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.4). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.5). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.6). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.7). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.8). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.9). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.0), the median unfolded MiniBooNE LEE under a $\nu_e$CC hypothesis. No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.1). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.2). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.3). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.4). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.5). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.6). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.7). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.8). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.9). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

7 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=2.0). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC PC channel. The 53-78th bins/rows/columns correspond to the $\nu_\mu$CC FC channel. The 79-104th bins/rows/columns correspond to the $\nu_\mu$CC PC channel. The 105-115th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 116-126th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 127-137th bins/rows/columns correspond to the NC$\pi^0$ channel.

Fully contained $\nu_e$CC $0pX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Partially contained $\nu_e$CC $0pX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Fully contained $\nu_e$CC $NpX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Partially contained $\nu_e$CC $NpX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Fully contained $\nu_\mu$CC $0pX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Partially contained $\nu_\mu$CC $0pX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Fully contained $\nu_\mu$CC $NpX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

Partially contained $\nu_\mu$CC $NpX\pi$ data, prediction, and LEE(x=1) predictions. Note that the rightmost bin is an overflow bin, containing all events with reconstructed neutrino energy greater than 2.5 GeV. The LEE(x=1) includes the predicted excess from an unfolding of the MiniBooNE LEE under a $\nu_e$CC hypothesis.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.0), the standard prediction with no low energy excess. No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.1). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.2). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.3). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.4). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.5). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.6). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.7). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.8). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=0.9). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.0), the median unfolded MiniBooNE LEE under a $\nu_e$CC hypothesis. No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.1). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.2). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.3). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.4). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.5). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.6). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.7). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.8). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=1.9). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

11 channel covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. For the data statistical uncertainty covariance matrix, (only diagonal elements, not included here), the Neyman, Pearson, or combined Neyman and Pearson (CNP) techniques can be used. This corresponds to LEE(x=2.0). No constraints have been applied at this stage. The 1-26th bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ FC channel. The 27-52nd bins/rows/columns correspond to the $\nu_e$CC $0pX\pi$ PC channel. The 53-78th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ FC channel. The 79-104th bins/rows/columns correspond to the $\nu_e$CC $NpX\pi$ PC channel. The 105-130th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ FC channel. The 131-156th bins/rows/columns correspond to the $\nu_\mu$CC $0pX\pi$ PC channel. The 157-182nd bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ FC channel. The 183-208th bins/rows/columns correspond to the $\nu_\mu$CC $NpX\pi$ PC channel. The 209-219th bins/rows/columns correspond to the CC$\pi^0$ FC channel. The 220-230th bins/rows/columns correspond to the CC$\pi^0$ PC channel. The 231-241st bins/rows/columns correspond to the NC$\pi^0$ channel.

Data/MC ratio as a function of reconstructed $\pi^0$ momentum for the 2$\gamma$0p selection

Energy spectra for the 1$\gamma$1p selected events. The figure shows the unconstrained background prediction and breakdowns as a function of reconstructed shower energy.

Energy spectra for the 1$\gamma$1p selected events. This figure shows the total background prediction with systematic uncertainty both before and after the 2$\gamma$ constraint.

Energy spectra for the 1$\gamma$0p selected events. The figure shows the unconstrained background prediction and breakdowns as a function of reconstructed shower energy.

Energy spectra for the 1$\gamma$0p selected events. This figure shows the total background prediction with systematic uncertainty both before and after the 2$\gamma$ constraint.

The observed event rate for the 1$\gamma$1p event sample, and comparisons to unconstrained (left) and constrained (right) background and LEE model predictions.

The observed event rate for the 1$\gamma$0p event sample, and comparisons to unconstrained (left) and constrained (right) background and LEE model predictions.

The fractional systematic covariance matrix for the final selected 1$\gamma$1p, 1$\gamma$0p, 2$\gamma$1p, and 2$\gamma$0p event samples, in bins of reconstructed shower energy for the 1$\gamma$ distributions, and reconstructed $\pi^{0}$ momentum for the 2$\gamma$ distributions. The bin boundaries are defined in Sec. II in supplementary material. This covariance matrix is evaluated at central value prediction and all variances have been rounded to three significant digits.

The fractional systematic covariance matrix for the final selected 1$\gamma$1p, 1$\gamma$0p, 2$\gamma$1p, and 2$\gamma$0p event samples, in bins of reconstructed shower energy for the 1$\gamma$1p signal, 1$\gamma$1p background, 1$\gamma$0p signal and 1$\gamma$0p background distributions, and reconstructed $\pi^{0}$ momentum for the 2$\gamma$ distributions. The bin boundaries are defined in Sec. II in supplementary material. This covariance matrix is evaluated at central value prediction and variances have been rounded to three significant digits.

Fourth $q^2$ moment in GeV$^8$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

First $q^2$ moment in GeV$^2$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Second $q^2$ moment in GeV$^4$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Third $q^2$ moment in GeV$^6$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Fourth $q^2$ moment in GeV$^8$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Second central $q^2$ moment in GeV$^4$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Third central $q^2$ moment in GeV$^6$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Fourth central $q^2$ moment in GeV$^8$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Second central $q^2$ moment in GeV$^4$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Third central $q^2$ moment in GeV$^6$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Fourth central $q^2$ moment in GeV$^8$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the first $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the first $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second central central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second central central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for prompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Fraction of nonprompt $J/\psi$ mesons (in \%) in ($p_\text{T},y$) intervals. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-0.2$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-1$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=+1$ rather than zero, in ($p_\text{T},y$) intervals.

The measured differential branching fractions as a function of the invariant hadronic mass squared of the $X_u$ system ($M_X^2$).

The measured differential branching fractions as a function of the light-cone momenta of the hadronic $X_u$ system [$P_+ = (E_X^B - |old{p}_X^B|)$].

The measured differential branching fractions as a function of the light-cone momenta of the hadronic $X_u$ system [$P_- = (E_X^B + |old{p}_X^B|)$].

The background-subtracted yields as a function of $E_\ell^B$.

The background-subtracted yields as a function of $q^2$.

The background-subtracted yields as a function of $M_X$.

The background-subtracted yields as a function of $M_X^2$.

The background-subtracted yields as a function of $P_+$. [$P_+ = (E_X^B - |P_X^B|)$]

The background-subtracted yields as a function of $P_-$. [$P_- = (E_X^B + |P_X^B|)$].

The full experimental correlation matrix for the background-subtracted spectrum of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The full experimental correlation matrix for the background-subtracted spectrum of $q^2$).

The full experimental correlation matrix for the background-subtracted spectrum of $M_X$.

The full experimental correlation matrix for the background-subtracted spectrum of $M_X^2$.

The full experimental correlation matrix for the background-subtracted spectrum of $P_+$.

The full experimental correlation matrix for the background-subtracted spectrum of $P_-$.

Migration matrix of $E_\ell^B$ [in %].

Migration matrix of $q^2$ [in %].

Migration matrix of $M_X$ [in %].

Migration matrix of $M_X^2$ [in %].

Migration matrix of $P_+$ [in %].

Migration matrix of $P_-$ [in %].

The total correction factors and phase space acceptance as a function of the invariant hadronic mass of the $X_u$ system ($M_X$).

The total correction factors and phase space acceptance as a function of the invariant hadronic mass squared of the $X_u$ system ($M_X^2$).

The total correction factors and phase space acceptance as a function of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The total correction factors and phase space acceptance as a function of the four-momentum-transfer squared of the $B$ to the $X_u$ system ($q^2$).

The total correction factors and phase space acceptance as a function of $P_+$.

The total correction factors and phase space acceptance as a function of $P_-$.

The statistical correlations of the differential branching fractions are shown. The matrix elements are ordered for $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$.

The full experimental (statistical and systematical) correlations of the differential branching fractions are shown. The matrix elements are ordered for $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$.

The total partial branching fraction with $E_{\ell}^{B} > 1$ GeV. The order of entries is: the 2D fit result of Phys. Rev. D 104, 012008 (2021), the results calcuated with the differential branching fractions of $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$

The full experimental correlations between the total partial branching fractions from summing the individual bins are shown. The matrix elements are saved in the order of $M_{X}$, $M_{X}^{2}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$.

The first moment of the measured differential branching fraction of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The second moment of the measured differential branching fraction of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The third moment of the measured differential branching fraction of the lepton energy in the $B$ rest frame ($E_\ell^B$).

The first moment of the measured differential branching fraction of $q^2$.

The second moment of the measured differential branching fraction of $q^2$.

The third moment of the measured differential branching fraction of $q^2$.

The first moment of the measured differential branching fraction of $M_{X}$.

The second moment of the measured differential branching fraction of $M_{X}$.

The third moment of the measured differential branching fraction of $M_{X}$.

The first moment of the measured differential branching fraction of $P_{+}$.

The second moment of the measured differential branching fraction of $P_{+}$.

The third moment of the measured differential branching fraction of $P_{+}$.

The first moment of the measured differential branching fraction of $P_{-}$.

The second moment of the measured differential branching fraction of $P_{-}$.

The third moment of the measured differential branching fraction of $P_{-}$.

The matrix elements are saved in the order of $M_{X}$, $q^{2}$, $E_{\ell}^{B}$, $P_{+}$ and $P_{-}$. For each variable, the entries is ordered for the first, second and third moments.

CL$_{s}$ value as a function of the branching fraction of $B^{+} \rightarrow K^{+} \nu \bar{\nu}$ for expected and observed signal yields and the corresponding upper limits at 90% confidence level (CL). The expected limit is derived for the background-only hypothesis. The observed limit is derived from a simultaneous fit to the on-resonance and off-resonance data, corresponding to an integrated luminosity of 63 fb$^{−1}$ and 9 fb$^{−1}$, respectively.

Signal efficiency as a function of the dineutrino invariant mass squared $q^{2}$ for events in the signal region (SR) ($\rm{BDT}_{1} > 0.9$ and $\rm{BDT}_{2} > 0.95$). The error bars indicate the statistical uncertainty.