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.

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology C.

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology D.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements.The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology C) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

B.R. times dN/dy of $^{3}_{\Lambda}$H vs y in 3 GeV 10-50% Au+Au collisions

B.R. times dN/dy of $^{4}_{\Lambda}$H vs y in 3 GeV 10-50% Au+Au collisions

B.R. times dN/dy at |y|<0.5 of $^{3}_{\Lambda}$H vs B.R in 3 GeV 0-10% Au+Au collisions

B.R. times dN/dy at |y|<0.5 of $^{4}_{\Lambda}$H vs B.R in 3 GeV 0-10% Au+Au collisions

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 0-10%

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 0-10%

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 10-50%

$^{3}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 10-50%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 0-10%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 0-10%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.75<y<-0.5, Au+Au 3 GeV, 0-10%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.25<y<0, Au+Au 3 GeV, 10-50%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.5<y<-0.25, Au+Au 3 GeV, 10-50%

$^{4}_{\Lambda}$H $p_T$ spectra times B.R., -0.75<y<-0.5, Au+Au 3 GeV, 10-50%

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF correlated with hard collision activity

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF without hard collision activity

Forward neutron single spin asymmetries in p+p collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Al collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Au collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

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.

$z_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

Fully corrected two subjet $z_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

$z_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

tracking efficiency in pp collisions as a function of track momentum

tracking efficiency in AuAu collisions as a function of track momentum

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.