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

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.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL including the BDT cut as a function of the signal mass.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL without the BDT cut as a function of the signal mass.

The number of events when applying each cut in the analysis for $m_a$ = 16 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 20 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 30 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 40 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 50 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 60 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

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

Distribution of the missing transverse energy in CRTop. The Z+HF and $t\bar{t}$ scale factors, described in the text, have been applied to the simulated samples. The distribution labeled "Signal" is for the model with ($m_a$, $m_{\tilde{\chi}_{1}^{0}}$, $m_{\tilde{\chi}_{2}^{0}}$) = (45 GeV, 10 GeV, 80 GeV).

Distribution of the dijet invariant mass in the signal region, shown together with the parameterized background model (labelled "Bkg Model"). For reference, the MC prediction for the SM background is also shown (labelled "SM MC"). The Z+HF and $t\bar{t}$ scale factors, described in the text, have been applied to the simulated samples. The signal region is defined to have dijet invariant mass > 20 GeV. The distribution labeled "Signal" is for the model with ($m_a$, $m_{\tilde{\chi}_{1}^{0}}$, $m_{\tilde{\chi}_{2}^{0}}$) = (45 GeV, 10 GeV, 80 GeV).

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 65$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 95$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 110$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 20$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 30$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 65$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 95$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 110$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 20$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 30$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Unweighted and weighted number of events after each stage of selection for the NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}}) = (45,10,80)$ GeV and the Z boson decaying to $e^{+}e^{-}, \;\mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The weighted number of events corresponds to an integrated luminosity of $139 \;\mathrm{fb}^{-1}$. The "skimming selection" required either a single electron or muon with $p_{T} > 100\;$ GeV or a pair of electrons or muons each with $p_{T} > 20\;$ GeV. The preselection requirements include the trigger, absence of a bad jet or bad muon, and two leptons, where a lepton is either an electron or a muon.

Acceptance and efficiency of this analysis for the signal models considered in this paper. The signal is an NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where the values of $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}})$ are varied. The Z boson decays to $e^{+}e^{-}, \mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The product of acceptance times efficiency is defined by the fraction of simulated events that pass all the selection criteria of this analysis. The acceptance is defined by the fraction of simulated events that pass all the selection criteria as applied to Monte Carlo truth-level quantities. The efficiency is then defined as the acceptance times efficiency, divided by the acceptance.

Acceptance and efficiency of this analysis for the signal models considered in this paper. The signal is an NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where the values of $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}})$ are varied. The Z boson decays to $e^{+}e^{-}, \mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The product of acceptance times efficiency is defined by the fraction of simulated events that pass all the selection criteria of this analysis. The acceptance is defined by the fraction of simulated events that pass all the selection criteria as applied to Monte Carlo truth-level quantities. The efficiency is then defined as the acceptance times efficiency, divided by the acceptance.

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.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.