Measured $\frac{d\sigma}{d(-t^{\prime})}~[\mathrm{nb/GeV^{2}}]$ for reaction $\gamma p\to \{\Lambda \bar{\Lambda}\} p$ including data of $6.5 \leq E_{\gamma} \leq 11.5$ [GeV], splitted in 10 energy bins (each as a column in the table). The observable $(-t^{\prime})$ is in unit of $[\mathrm{nb/GeV^{2}}]$ and is divided into bins of width 0.03 $[\mathrm{GeV^{2}}]$ (each as a row in the table). The global systematic uncertainty is 19% (not included in the table), with contributions of 5% from kinematic fitting, 10% from data selection, 5% from flux normalization, 13% from tracking efficiency, 3% from model dependence, and 6% from run-period variations.

Measured $\frac{d\sigma}{d(-t^{\prime})}~[\mathrm{nb/GeV^{2}}]$ for reaction $\gamma p\to \{p \bar{\Lambda}\} \Lambda$ including data of $6.5 \leq E_{\gamma} \leq 11.5$ [GeV], splitted in 10 energy bins (each as a column in the table). The observable $(-t^{\prime})$ is in unit of $[\mathrm{nb/GeV^{2}}]$ and is divided into bins of width 0.075 $[\mathrm{GeV^{2}}]$ (each as a row in the table). The global systematic uncertainty is 22% (not included in the table), with contributions of 2% from kinematic fitting, 10% from data selection, 5% from flux normalization, 15% from tracking efficiency, 3% from model dependence, and 10% from run-period variations.

Measured $\frac{d\sigma}{d(-t^{\prime})}~[\mathrm{nb/GeV^{2}}]$ for reaction $\gamma p\to \{p \bar{p}\} p$ including data of $6.5 \leq E_{\gamma} \leq 11.5$ [GeV], splitted in 10 energy bins (each as a column in the table). The observable $(-t^{\prime})$ is in unit of $[\mathrm{nb/GeV^{2}}]$ and is divided into bins of width 0.01 $[\mathrm{GeV^{2}}]$ (each as a row in the table). The global systematic uncertainty is 13% (not included in the table), with contributions of 8% from kinematic fitting, 4% from data selection, 5% from flux normalization, 8% from tracking efficiency, 3% from model dependence, and 1% from run-period variations.

Measured $\frac{d\sigma}{d(-t^{\prime})}~[\mathrm{nb/GeV^{2}}]$ for reaction $\gamma p\to \{p \bar{p}\} p$ including data of $3.5 \leq E_{\gamma} \leq 6.0$ [GeV], splitted in 5 energy bins (each as a column in the table). The observable $(-t^{\prime})$ is in unit of $[\mathrm{nb/GeV^{2}}]$ and is divided into bins of width 0.1 $[\mathrm{GeV^{2}}]$ (each as a row in the table). The global systematic uncertainty is 13% (not included in the table), with contributions of 8% from kinematic fitting, 4% from data selection, 5% from flux normalization, 8% from tracking efficiency, 3% from model dependence, and 1% from run-period variations.

Measured total cross section $\sigma(\gamma\,p\to\{\Lambda\bar{\Lambda}\}p)$ in unit of [nb] with beam energy in the range $5.0 < E_{\gamma} < 11.5~[\mathrm{GeV}]$ with various bin width. The table only includes point-to-point statistical uncertainties. The global systematic uncertainty is 19% (not included in the table), with contributions of 5% from kinematic fitting, 10% from data selection, 5% from flux normalization, 13% from tracking efficiency, 3% from model dependence, and 6% from run-period variations.

Measured total cross section $\sigma(\gamma\,p\to\{p\bar{\Lambda}\}\Lambda)$ in unit of [nb] with beam energy in the range $5.0 < E_{\gamma} < 11.5~[\mathrm{GeV}]$ with various bin width. The table only includes point-to-point statistical uncertainties. The global systematic uncertainty is 22% (not included in the table), with contributions of 2% from kinematic fitting, 10% from data selection, 5% from flux normalization, 15% from tracking efficiency, 3% from model dependence, and 10% from run-period variations.

Measured total cross section $\sigma(\gamma\,p\to\{p\bar{p}\}p)$ in unit of [nb] with beam energy in the range $3.5 < E_{\gamma} < 11.5~[\mathrm{GeV}]$ with various bin width. The table only includes point-to-point statistical uncertainties. The global systematic uncertainty is 13% (not included in the table), with contributions of 8% from kinematic fitting, 4% from data selection, 5% from flux normalization, 8% from tracking efficiency, 3% from model dependence, and 1% from run-period variations.

Observed 95% $\text{CL}_\text{S}$ upper limits on the production cross-section times acceptance times branching ratio to jets, $\sigma \cdot A \cdot \text{BR}$, of Gaussian-shaped signals of 5%, 10%, and 15% width relative to their peak mass, $m_G$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J100 signal region.

Observed 95% $\text{CL}_\text{S}$ upper limits on the the value of the coupling to quarks, $g_q$, of the $Z^\prime$ model as a function of the resonance mass $m_{Z^\prime}$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J50 signal region.

Observed 95% $\text{CL}_\text{S}$ upper limits on the the value of the coupling to quarks, $g_q$, of the $Z^\prime$ model as a function of the resonance mass $m_{Z^\prime}$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J100 signal region.

Acceptance of simulated $Z^\prime$ signal events given the analysis event selection as a function of $m_{Z^\prime}$. The solid black line corresponds to the acceptance without a lower $m_{jj}$ requirement applied, the solid blue line includes the $m_{jj} > 344$ GeV requirement of the J50 signal region, and the solid red line includes the $m_{jj} > 481$ GeV requirement of the J100 signal region.

Numbers of events in each of the two considered datasets after consecutively applying the analysis selections.

Ratio between the prompt $\Xi_c^+$ baryons in a multiplicity class to the multiplicity-integrated (INEL $>$ 0) class.

The prompt production yield ratios between $\Xi_c^0$ and $\rm {D}^0$ measured in the same multiplicity classes.

The prompt production yield ratios between $\Xi_c^+$ and $\rm {D}^0$ measured in the same multiplicity classes.

The prompt production yield ratios between $\Xi_c^0$ and $\Lambda_c^+$ measured in the same multiplicity classes.

The prompt production yield ratios between $\Xi_c^+$ and $\Lambda_c^+$ measured in the same multiplicity classes.

The branching-fraction ratio of BR($\Xi_c^0 \rightarrow \Xi^-e^+\nu_e$)/BR($\Xi_c^0 \rightarrow \Xi^-\pi^+$)

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.

$p_\mathrm{{T}}$-differential production yield of $\overline{\Sigma}^{-}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $\overline{\Sigma}^{+}$/ $\overline{\Sigma}^{-}$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $\overline{\Sigma}^{+}$/ $\overline{\Sigma}^{-}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(p+\overline{p})/(\pi^+ + \pi^-)$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(\pi^+ + \pi^-)$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(\pi^+ + \pi^-)$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(p + \overline{p})$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(p + \overline{p})$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(p+\overline{p})/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.5<y_\mathrm{CMS}<0$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(\Lambda + \overline{\Lambda})/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.5<y_\mathrm{CMS}<0$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(p + \overline{p})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(p + \overline{p})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(\Lambda + \overline{\Lambda})/(p + \overline{p})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.5<y_\mathrm{CMS}<0$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(\Lambda + \overline{\Lambda})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(\Lambda + \overline{\Lambda})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

$p_\mathrm{{T}}$-differential nuclear modification factor $R_\mathrm{{pPb}}$ of $\overline{\Sigma}^{+}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

$p_\mathrm{{T}}$-differential nuclear modification factor $R_\mathrm{{pPb}}$ of $\overline{\Sigma}^{-}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Correlation matrix for the correlated systematic uncertainties of the forward rapidity $\eta$ meson cross section.

Correlation matrix for the correlated systematic uncertainties of the central rapidity $\eta$ meson cross section.

The invariant yields of $\phi$ and $\omega+\rho$ mesons as a function of $\langle N_{\rm part}\rangle$ in Au+Au collisions. The systematic uncertainties of type-A (uncorrelated) are combined with statistical uncertainties in quadrature and are labeled as stat. Type-B (correlated) systematic uncertainties are listed as sys.

$R_{AA}$ of $\phi$ meson as a function of $p_T$ for in Au+Au collisions at $\sqrt{s_{_{NN}}}=200$~GeV, compared with Cu+Au collisions at $1.2<|y|<2.2$. The systematic uncertainties of type-A (uncorrelated) are combined with statistical uncertainties in quadrature and are labeled as stat. Type-B (correlated) systematic uncertainties are listed as sys.

$R_{AA}$ of $\phi$ meson as a function of $\langle N_{\rm part}\rangle$ for $1.2<|y|<2.2$ in Au+Au collisions at $\sqrt{s_{_{NN}}}=200$~GeV. The systematic uncertainties of type-A (uncorrelated) are combined with statistical uncertainties in quadrature and are labeled as stat. Type-B (correlated) systematic uncertainties are listed as sys.

Observed profile likelihood scans of $\kappa_\lambda$. The scans are performed by varying only the coupling modifier of interest, while all other relevant coupling modifiers are fixed to unity.

Expected profile likelihood scans of $\kappa_\lambda$. The scans are performed by varying only the coupling modifier of interest, while all other relevant coupling modifiers are fixed to unity.

Observed profile likelihood scans of $\kappa_{2V}$. The scans are performed by varying only the coupling modifier of interest, while all other relevant coupling modifiers are fixed to unity.

Expected profile likelihood scans of $\kappa_{2V}$. The scans are performed by varying only the coupling modifier of interest, while all other relevant coupling modifiers are fixed to unity.

Confidence level contours at 68% (solid line) and 95% (dashed line) in the $(\kappa_\lambda, \kappa_{2V})$ parameter space, when all other coupling modifiers are fixed to their SM predictions. The corresponding expected contours are shown by the inner and outer shaded regions The SM prediction is indicated by the star, while the best-fit value is denoted by the cross.

Confidence level contours at 68% (solid line) and 95% (dashed line) in the $(\kappa_\lambda, \kappa_{2V})$ parameter space, when all other coupling modifiers are fixed to their SM predictions. The corresponding expected contours are shown by the inner and outer shaded regions The SM prediction is indicated by the star, while the best-fit value is denoted by the cross.

Observed $m_{\gamma\gamma}$ distributions for the high-mass categories in Run 2. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the high-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the high-mass categories in Run 2. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the high-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the high-mass categories in Run 2. 3he lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the high-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 2. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 2. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 2. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 2. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed $m_{\gamma\gamma}$ distributions for the low-mass categories in Run 3. The lines show the fit results for the continuum background only (light dotted), adding the single Higgs boson background (black dotted) and the full fit (solid).

Observed and expected profiled likelihood scans on $\kappa_\lambda$ for the low-mass (LM) and high-mass (HM) regions separately and for their combination. The scan is performed by floating the indicated parameter of interest while fixing all the others to their SM value.

Observed and expected profiled likelihood scans on $\kappa_\lambda$ for the low-mass (LM) and high-mass (HM) regions separately and for their combination. The scan is performed by floating the indicated parameter of interest while fixing all the others to their SM value.

Observed and expected profiled likelihood scans on $\kappa_\lambda$ for the low-mass (LM) and high-mass (HM) regions separately and for their combination. The scan is performed by floating the indicated parameter of interest while fixing all the others to their SM value.

Observed and expected profiled likelihood scans on $\kappa_\lambda$ for the low-mass (LM) and high-mass (HM) regions separately and for their combination. The scan is performed by floating the indicated parameter of interest while fixing all the others to their SM value.

Observed and expected profiled likelihood scans on $\kappa_\lambda$ for the low-mass (LM) and high-mass (HM) regions separately and for their combination. The scan is performed by floating the indicated parameter of interest while fixing all the others to their SM value.

Observed and expected profiled likelihood scans on $\kappa_\lambda$ for the low-mass (LM) and high-mass (HM) regions separately and for their combination. The scan is performed by floating the indicated parameter of interest while fixing all the others to their SM value.

The expected number of events (estimated by using simulation) from the SM HH signals and single Higgs boson production, and the expected number of events from the continuum background, evaluated in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window in Run 3 for the different low-mass (LM) and high-mass (HM) categories. For comparison, the number of data events is also shown. The uncertainties in the HH signals and single Higgs boson backgrounds include the systematic uncertainties discussed in Section 6. Asymmetric uncertainties arise primarily from the theory calculation of the SM ggF HH cross section and the large uncertainty in the yield of single Higgs bosons produced in ggF events in association with heavy-flavour jets. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious signal uncertainty.