Single Event Sensitivity (SES) for the $K^{+}\rightarrow\pi^{+}X$ search as a function of X mass.

Model-independent constraints on the branching ratio of the $K^{+}\rightarrow\pi^{+}X$ decay

Model-independent constraints on the branching ratio of the $K^{+}\rightarrow\pi^{+}X$ decay

Observed model-independent upper limits at 90 % CL of $\mathcal{B}(K^{+}\rightarrow\pi^{+}X)$ as function of $X$ mass, for several $X$ lifetime hypotheses, assuming $X$ decays to visible SM particles.

Observed model-independent upper limits at 90 % CL of $\mathcal{B}(K^{+}\rightarrow\pi^{+}X)$ as function of $X$ mass, for several $X$ lifetime hypotheses, assuming $X$ decays to visible SM particles.

Di-muon mass spectrum of selected $K^{+}\rightarrow\pi^{+}\mu^{+}\mu^{-}$ events from 2017–2018 data.

Di-muon mass spectrum of selected $K^{+}\rightarrow\pi^{+}\mu^{+}\mu^{-}$ events from 2017–2018 data.

Expected and observed model-independent upper limits at 90 % CL for $\mathcal{B}(K^{+}\rightarrow\pi^{+}X)\times\mathcal{B}(X\rightarrow\mu^{+}\mu^{-})$ as a function of $m_{X}$ for $\tau_{X}=0$.

Expected and observed model-independent upper limits at 90 % CL for $\mathcal{B}(K^{+}\rightarrow\pi^{+}X)\times\mathcal{B}(X\rightarrow\mu^{+}\mu^{-})$ as a function of $m_{X}$ for $\tau_{X}=0$.

Observed model-independent upper limits at 90 % CL for $\mathcal{B}(K^{+}\rightarrow\pi^{+}X)\times\mathcal{B}(X\rightarrow\mu^{+}\mu^{-})$ for several $\tau_{X}$ values.

Observed model-independent upper limits at 90 % CL for $\mathcal{B}(K^{+}\rightarrow\pi^{+}X)\times\mathcal{B}(X\rightarrow\mu^{+}\mu^{-})$ for several $\tau_{X}$ values.

Branching ratio of the $K^{+}\rightarrow\pi^{+}A^{\prime}$ decay divided by the kinetic mixing coupling squared, $\varepsilon^{2}$, as a function of $m_{A^{\prime}}$ for the BC2 model with dark photon $A^{\prime}$.

Branching ratio of the $K^{+}\rightarrow\pi^{+}A^{\prime}$ decay divided by the kinetic mixing coupling squared, $\varepsilon^{2}$, as a function of $m_{A^{\prime}}$ for the BC2 model with dark photon $A^{\prime}$.

Excluded region, at 90 % CL, of the parameter space $(m_{A^{\prime}},\varepsilon)$ for a dark photon $A^{\prime}$, decaying invisibly, in the BC2 model for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded region, at 90 % CL, of the parameter space $(m_{A^{\prime}},\varepsilon)$ for a dark photon $A^{\prime}$, decaying invisibly, in the BC2 model for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded region, at 90 % CL, of the parameter space $(m_{A^{\prime}},\varepsilon)$ for a dark photon $A^{\prime}$, decaying invisibly, in the BC2 model for the NA62 $\pi^{0}\rightarrow{\rm inv}$ search.

Excluded region, at 90 % CL, of the parameter space $(m_{A^{\prime}},\varepsilon)$ for a dark photon $A^{\prime}$, decaying invisibly, in the BC2 model for the NA62 $\pi^{0}\rightarrow{\rm inv}$ search.

Branching ratio of the $K^{+}\rightarrow\pi^{+}S$ decay divided by $\text{sin}^{2}\theta$, as a function of $m_{S}$.

Branching ratio of the $K^{+}\rightarrow\pi^{+}S$ decay divided by $\text{sin}^{2}\theta$, as a function of $m_{S}$.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta)$ for a dark scalar $S$, in the BC4-inv model for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta)$ for a dark scalar $S$, in the BC4-inv model for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta)$ for a dark scalar $S$, in the BC4-inv model for the NA62 $\pi^{0}\rightarrow{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta)$ for a dark scalar $S$, in the BC4-inv model for the NA62 $\pi^{0}\rightarrow{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta$) for a dar scalar S, in the BC4 model, from the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta$) for a dar scalar S, in the BC4 model, from the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta$) for a dar scalar S, in the BC4 model, from the NA62 $K^{+}\rightarrow\pi^{+}\mu^{+}\mu^{-}$ study.

Excluded regions, at 90 % CL, of the parameter space $(m_{S},\rm{sin}^{2}\theta$) for a dar scalar S, in the BC4 model, from the NA62 $K^{+}\rightarrow\pi^{+}\mu^{+}\mu^{-}$ study.

Branching ratio of the $K^{+}\rightarrow\pi^{+}a$ decay divided by $(C_{ff}/\Lambda)^{2}$, as a function of $m_{a}$, assuming Λ = 1 TeV.

Branching ratio of the $K^{+}\rightarrow\pi^{+}a$ decay divided by $(C_{ff}/\Lambda)^{2}$, as a function of $m_{a}$, assuming Λ = 1 TeV.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\Lambda)$ for an ALP $a$, in the BC10-inv model, evaluated assuming $\Lambda = 1$ TeV, for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\Lambda)$ for an ALP $a$, in the BC10-inv model, evaluated assuming $\Lambda = 1$ TeV, for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\Lambda)$ for an ALP $a$, in the BC10-inv model, evaluated assuming $\Lambda = 1$ TeV, for the NA62 $\pi^{0}\rightarrow\rm{inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\Lambda)$ for an ALP $a$, in the BC10-inv model, evaluated assuming $\Lambda = 1$ TeV, for the NA62 $\pi^{0}\rightarrow\rm{inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\lambda)$ for an ALP $a$, in the BC10 model, evaluated assuming $\Lambda = 1$ TeV for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\lambda)$ for an ALP $a$, in the BC10 model, evaluated assuming $\Lambda = 1$ TeV for the NA62 $K^{+}\rightarrow\pi^{+}X_{\rm inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\lambda)$ for an ALP $a$, in the BC10 model, evaluated assuming $\Lambda = 1$ TeV for the NA62 $\pi^{0}\rightarrow\rm{inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\lambda)$ for an ALP $a$, in the BC10 model, evaluated assuming $\Lambda = 1$ TeV for the NA62 $\pi^{0}\rightarrow\rm{inv}$ search.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\lambda)$ for an ALP $a$, in the BC10 model, evaluated assuming $\Lambda = 1$ TeV for the NA62 $K^{+}\rightarrow\pi^{+}\mu^{+}\mu^{-}$ study.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{ff}/\lambda)$ for an ALP $a$, in the BC10 model, evaluated assuming $\Lambda = 1$ TeV for the NA62 $K^{+}\rightarrow\pi^{+}\mu^{+}\mu^{-}$ study.

Branching ratio of the $K^{+}\rightarrow\pi^{+}a$ decay divided by $(C_{GG}/\Lambda)^{2}$, as a function of $m_{a}$, assuming $\Lambda = 1$ TeV.

Branching ratio of the $K^{+}\rightarrow\pi^{+}a$ decay divided by $(C_{GG}/\Lambda)^{2}$, as a function of $m_{a}$, assuming $\Lambda = 1$ TeV.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{GG}/\Lambda)$ for an ALP $a$, of the BC11 model, evaluated assuming $\Lambda = 1$ TeV, from the NA62 search for $K^{+}\rightarrow\pi^{+}X_{\rm inv}$.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{GG}/\Lambda)$ for an ALP $a$, of the BC11 model, evaluated assuming $\Lambda = 1$ TeV, from the NA62 search for $K^{+}\rightarrow\pi^{+}X_{\rm inv}$.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{GG}/\Lambda)$ for an ALP $a$, of the BC11 model, evaluated assuming $\Lambda = 1$ TeV, from the NA62 search for $\pi^{0}\rightarrow\rm{inv}$.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{GG}/\Lambda)$ for an ALP $a$, of the BC11 model, evaluated assuming $\Lambda = 1$ TeV, from the NA62 search for $\pi^{0}\rightarrow\rm{inv}$.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{GG}/\Lambda)$ for an ALP $a$, of the BC11 model, evaluated assuming $\Lambda = 1$ TeV, from the NA62 $K^{+}\rightarrow\pi^{+}\gamma\gamma$ study.

Excluded regions, at 90 % CL, of the parameter space $(m_{a}, C_{GG}/\Lambda)$ for an ALP $a$, of the BC11 model, evaluated assuming $\Lambda = 1$ TeV, from the NA62 $K^{+}\rightarrow\pi^{+}\gamma\gamma$ study.

Data from Figure 3, panel a, U+U

Data from Figure 3, panel b, Au+Au

Data from Figure 3, panel b, U+U

Data from Figure 6, panel a, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 6, panel b, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 6, panel c, 0-0.5% Centrality, 0.2<p_{T}<3 GeV/c

Data from Figure 22, panel a, U+U/Au+Au

Data from Figure 22, panel a, U+U/Au+Au

Data from Figure 30, upper row panel a, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, upper row panel b, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, upper row panel c, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, upper row panel d, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, middle row panel a, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, middle row panel b, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, middle row panel c, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, middle row panel d, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, lower row panel a, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, lower row panel b, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, lower row panel c, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 30, lower row panel d, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel a, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel a, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel b, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel b, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel c, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel c, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel d, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel d, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel e, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel e, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel f, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel f, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel g, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel g, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel g, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel g, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel i, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel i, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel g, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel g, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel k, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel k, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel l, Au+Au, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 31, panel l, U+U, 2subevent method, 0.2<p_{T}<3 GeV/c

Data from Figure 32, upper row panel a, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, upper row panel b, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, upper row panel c, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, upper row panel d, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, middle row panel a, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, middle row panel b, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, middle row panel c, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, middle row panel d, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, lower row panel a, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, lower row panel b, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, lower row panel c, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 32, lower row panel d, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel a, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel a, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel b, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel b, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel c, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel c, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel d, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel d, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel e, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel e, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel f, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel f, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel g, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel g, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel g, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel g, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel i, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel i, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel g, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel g, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel k, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel k, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel l, Au+Au, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 33, panel l, U+U, 2subevent method, 0.5<p_{T}<3 GeV/c

Data from Figure 34, upper row panel a, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, upper row panel b, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, upper row panel c, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, upper row panel d, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, middle row panel a, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, middle row panel b, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, middle row panel c, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, middle row panel d, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, lower row panel a, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, lower row panel b, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, lower row panel c, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 34, lower row panel d, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel a, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel a, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel b, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel b, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel c, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel c, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel d, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel d, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel e, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel e, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel f, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel f, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel g, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel g, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel g, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel g, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel i, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel i, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel g, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel g, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel k, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel k, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel l, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 35, panel l, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 36, upper row panel a, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, upper row panel b, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, upper row panel c, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, upper row panel d, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, middle row panel a, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, middle row panel b, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, middle row panel c, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, middle row panel d, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, lower row panel a, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, lower row panel b, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, lower row panel c, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 36, lower row panel d, 2subevent method, 0.5<p_{T}<2 GeV/c

Data from Figure 37, panel a, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel a, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel b, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel b, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel c, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel c, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel d, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel d, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel e, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel e, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel f, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel f, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel g, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel g, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel g, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel g, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel i, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel i, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel g, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel g, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel k, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel k, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel l, Au+Au, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 37, panel l, U+U, 2subevent method, 0.2<p_{T}<2 GeV/c

Data from Figure 38, upper row panel a, 2subevent method, four p_{T} ranges

Data from Figure 38, upper row panel b, 2subevent method, four p_{T} ranges

Data from Figure 38, upper row panel c, 2subevent method, four p_{T} ranges

Data from Figure 38, upper row panel d, 2subevent method, four p_{T} ranges

Data from Figure 38, middle row panel a, 2subevent method, four p_{T} ranges

Data from Figure 38, middle row panel b, 2subevent method, four p_{T} ranges

Data from Figure 38, middle row panel c, 2subevent method, 0.2< p_{T} <2 Gev/c, 0.5<p_{T}<2 GeV/c

Data from Figure 38, middle row panel c, 2subevent method, 0.2< p_{T} <3 Gev/c, 0.5<p_{T}<3 GeV/c

Data from Figure 38, middle row panel d, 2subevent method, four p_{T} ranges

Data from Figure 38, lower row panel a, 2subevent method, four p_{T} ranges

Data from Figure 38, lower row panel b, 2subevent method, four p_{T} ranges

Data from Figure 38, lower row panel c, 2subevent method, 0.2< p_{T} <2 Gev/c, 0.5<p_{T}<2 GeV/c

Data from Figure 38, lower row panel c, 2subevent method, 0.2< p_{T} <3 Gev/c, 0.5<p_{T}<3 GeV/c

Data from Figure 38, lower row panel d, 2subevent method, four p_{T} ranges

Corrected Yield R=0.5 pi0+jet 15-20 pp at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 gamma+jet 10-15 pp at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 gamma+jet 10-15 pp at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 pi0+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 pi0+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 pi0+jet 15-20 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 pi0+jet 15-20 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 gamma+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 gamma+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

IAA(Δφ) R=0.5 gamma+jet 10-15 GeV/c

IAA(Δφ) R=0.2 gamma+jet 10-15 GeV/c

IAA(Δφ) R=0.5 pi0+jet 10-15 GeV/c

IAA(Δφ) R=0.2 pi0+jet 10-15 GeV/c

IAA(Δφ) R=0.5 pi0+jet 15-20 GeV/c

IAA(Δφ) R=0.2 pi0+jet 15-20 GeV/c

Observed best fit differential cross section for the $|y_{H}|$ observable

Observed best fit differential cross section for the $m_{jj}$ (GeV) observable

Observed best fit differential cross section for the $\Delta\eta_{jj}$ observable.

Observed best fit differential cross section for the $\tau_{C}^{j}$ (GeV) observable

Bin-to-bin correlation matrix for the $p_{T}^{H}$ spectrum

Bin-to-bin correlation matrix for the $N_{jets}$ spectrum

Bin-to-bin correlation matrix for the $p_{T}^{j0}$ spectrum

Bin-to-bin correlation matrix for the $|y_{H}|$ spectrum

Bin-to-bin correlation matrix for the $m_{jj}$ spectrum

Bin-to-bin correlation matrix for the $\Delta\eta_{jj}$ spectrum

Bin-to-bin correlation matrix for the $\tau_{C}^{j}$ spectrum

Rotation matrix for the PCA of the SMEFT Wilson coefficients

Correlation matrix of the linear combinations of Wilson coefficients obtained from the PCA, obtained by fitting the observed data to the $p_{T}^{H}$ spectra of all decay channels.

Summary of observed and expected confidence intervals for the first eigenvectors obtained from the PCA of the SMEFT Wilson coefficients.

90% CL upper limit in axion-like particle coupling to SM fermions at UV scale Lambda=1TeV vs mass parameter space for di-lepton final states.

90% CL upper limit in axion-like particle coupling to SM gluons at UV scale Lambda=1TeV vs mass parameter space for hadronic final states.

Number of expected $a \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $C_{bs}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\eta$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $C_{bs}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\gamma$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $C_{bs}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $C_{bs}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $C_{bs}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\eta$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\gamma$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\pi^0$ assuming production coupling $\theta_{\pi^0} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\eta$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\gamma$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta$ assuming production coupling $\theta_{\eta} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\eta$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\gamma$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in mixing with $\eta^\prime$ assuming production coupling $\theta_{\eta^\prime} = 1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\eta$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\gamma$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $a \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Primakoff production assuming production coupling $C_{\gamma\gamma}/\Lambda=1\,\mathrm{GeV}^{-1}$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to K^+K^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to K^+K^-$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to K^+K^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for mixing production assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in light meson decays assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in light meson decays assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in light meson decays assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in light meson decays assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $A^\prime \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for production in light meson decays assuming production coupling $\varepsilon=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to K^+K^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to \pi^+\pi^-$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to \pi^+\pi^-\pi^0\pi^0$ events for BR = 1 in the plane (mass, width) after full selection, for Bremsstrahlung production assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to e^+e^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to \mu^+\mu^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to K^+K^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

Number of expected $S \to \pi^+\pi^-$ events for BR = 1 in the plane (mass, width) after full selection, for production in B meson decays assuming production coupling $\mathrm{sin}^2\theta=1$. Acceptance for particles that reach the FV and decay therein and mass resolution.

The relative systematic uncertainties of the fitted number of neutrino interactions.

The relative systematic uncertainties of the fitted number of anti-neutrino interactions.

The relative systematic uncertainties of the fitted number of neutrino and anti-neutrino interactions.

Number of observed CC $\nu_{\mu}$ neutrino candidate events in data and simulation.

Number of observed CC $\nu_{\mu}$ neutrino candidate events per inv. bin width.

Number of CC $\nu_{\mu$} neutrino interactions in the fiducial volume.

Number of CC $\nu_{\mu$} neutrino interactions in the fiducial volume per bin width.

Number of CC $\bar{nu}_{\mu$} neutrino interactions in the fiducial volume per bin width.

Number of CC $\nu_{\mu}$ and $\bar{nu}_{\mu$} neutrino interactions in the fiducial volume per bin width.

Neutrino cross section

Anti-neutrino cross section

Neutrino and anti-neutrino cross section

Neutrino cross section divided by energy.

Anti-neutrino cross section divided by energy.

Neutrino and anti-neutrino cross section divided by energy.

Neutrino flux divided by the cross sectional area of the fiducial volume.

Anti-Neutrino flux divided by the cross sectional area of the fiducial volume.

Neutrino and anti-Neutrino flux divided by the cross sectional area of the fiducial volume.

Neutrino flux per bin width divided by the cross sectional area of the fiducial volume.

Anti-neutrino flux per bin width divided by the cross sectional area of the fiducial volume.

Neutrino and anti-neutrino flux per bin width divided by the cross sectional area of the fiducial volume.

Number of neutrino interactions from pion and kaon decays.

Covariance matrix of the fitted number of neutrinos from pion and kaon decays.

Response matrix to unfold the reconstructed, calibrated muon momentum $q/p^{'}_{\mu}$ to the true neutrino energy $-L/E_{\nu}$.

Distribution of the number of clusters in the IFT for neutrino-like data, muon-like data, and simulation.

Distribution of the muon angle for neutrino-like data, muon-like data, and simulation.

Distribution of the muon charge over momentum for neutrino-like data, muon-like data, and simulation.

Efficiency in bins of the neutrino energy $-L/E_{\nu}$.

Covariance matrix from the likelihood of the number of neutrino interactions in bins of the neutrino energy $-L/E_{\nu}$.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 and ISR starting scale SpaceShower:pT0Ref = 1 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 and ISR starting scale SudakovCommon:pTmin=0.7 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

The charged multiplicity normalzied excess yield of 54.4 GeV in 0-80% centrality.

T vs. uB from 27 GeV and 54.4 GeV