The observed and expected ($\pm2\sigma$) upper limits at 95% CL on the branching ratio for $H\rightarrow aa\rightarrow \gamma\gamma\tau\tau$ as a function of the resonance mass hypothesis $m_{a}$.

The polynomial parameterisation of the total signal selection efficiency as a function of $m_{a}$, defined as the fraction of the reconstructed and selected events out of the total number of events in the signal simulated samples.

Invariant mass $m_{J}$ of the resonance candidates in the region defined with forward photon $\eta_{\gamma} > 1.3$ and anti-tagged large-$R$ jet after the fit to data under the background-only hypothesis. The total systematic uncertainty is shown as the hatched band. Three representative $Z^{`}$ signal distributions are overlaid as red lines. The signal is shown for $g_q=0.2$ with production cross sections of 309 fb, 143 fb, and 34.2 fb for $m_{Z^{`}}=(20,~50,~\text{and}~125~\text{GeV}$), respectively.

Observed (expected) upper exclusion limits at 95$\%$ CL on the coupling strength $g_q$ of a new vector $Z^{`}$ particle decaying to a $q\bar q$ pair for a LHC DM WG benchmark signal model where decay rates for all decay modes except $Z^{`}\to q \bar q$ are set to 0, and the interference between the $Z^{`}$ and the SM $Z$ boson is neglected.

Observed (solid lines) and expected (dashed lines) upper exclusion limits at 95$\%$ CL on the product of the $Z^{`}$ production production cross section, branching ratio for $Z^{`}\to q\bar q$ and acceptance ($\sigma_{Z^{`}}\times B(Z^{`}\to q\bar q)\times A$) of a new $Z^{`}$ resonance for a LHC DM WG benchmark signal model where all decay modes except $Z^{`}\to q \bar q$ are set to 0, and the interference between the $Z^{`}$ and the SM $Z$ boson is neglected.

Correlation matrix of all the parameters including nuisance parameters used in the profile likelihood unfolding fit on $A_{C}$.

Measured inclusive fiducial $tt\gamma$ production cross section in the dilepton final state for the different dilepton-flavour channels and combined.

Absolute differential $tt\gamma$ production cross section as a function of $p_{T}(\gamma)$ . The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $|\eta |(\gamma)$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of min $\Delta R(\gamma, \ell)$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $\Delta R(\gamma, \ell_{1})$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $\Delta R(\gamma, \ell_{2})$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of min $\Delta R(\gamma, b)$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $|\Delta\eta(\ell\ell)|$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $\Delta \phi(\ell\ell)$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $p_{T}(\ell\ell) $. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $p_{T}(\ell_{1})+p_{T}(\ell_{2})$ . The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of min $\Delta R(\ell, j)$. The values provided in the table are not divided by the bin width.

Absolute differential $tt\gamma$ production cross section as a function of $p_{T}(j_{1})$ .

Normalized differential $tt\gamma$ production cross section as a function of $p_{T}(\gamma)$ . The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $|\eta |(\gamma)$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of min $\Delta R(\gamma, \ell)$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $\Delta R(\gamma, \ell_{1})$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $\Delta R(\gamma, \ell_{2})$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of min $\Delta R(\gamma, b)$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $|\Delta\eta(\ell\ell)|$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $\Delta \phi(\ell\ell)$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $p_{T}(\ell\ell) $. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $p_{T}(\ell_{1})+p_{T}(\ell_{2})$ . The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of min $\Delta R(\ell, j)$. The values provided in the table are not divided by the bin width.

Normalized differential $tt\gamma$ production cross section as a function of $p_{T}(j_{1})$ . The values provided in the table are not divided by the bin width.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $p_{T}(\gamma)$ .

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $p_{T}(\gamma)$ .

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $|\eta |(\gamma)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $|\eta |(\gamma)$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of min $\Delta R(\gamma, \ell)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of min $\Delta R(\gamma, \ell)$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $\Delta R(\gamma, \ell_{1})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $\Delta R(\gamma, \ell_{1})$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $\Delta R(\gamma, \ell_{2})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $\Delta R(\gamma, \ell_{2})$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of min $\Delta R(\gamma, b)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of min $\Delta R(\gamma, b)$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $|\Delta\eta(\ell\ell)|$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $|\Delta\eta(\ell\ell)|$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $\Delta \phi(\ell\ell)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $\Delta \phi(\ell\ell)$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $p_{T}(\ell\ell) $.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $p_{T}(\ell\ell) $.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $p_{T}(\ell_{1})+p_{T}(\ell_{2})$ .

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $p_{T}(\ell_{1})+p_{T}(\ell_{2})$ .

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of min $\Delta R(\ell, j)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of min $\Delta R(\ell, j)$.

Correlation matrix of the systematic uncertainty in the absolute differential cross section as a function of $p_{T}(j_{1})$ .

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of $p_{T}(j_{1})$ .

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c_{tZ}$, using the photon pT distribution from the dilepton analysis. The value of $c^{I}_{tZ}$ is fixed to zero in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c_{tZ}$, using the combination of photon pT distributions from the dilepton and lepton+jets analyses. The value of $c^{I}_{tZ}$ is fixed to zero in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c^{I}_{tZ}$, using the photon pT distribution from the dilepton analysis. The value of $c_{tZ}$ is fixed to zero in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c^{I}_{tZ}$, using the combination of photon pT distributions from the dilepton and lepton+jets analyses. The value of $c_{tZ}$ is fixed to zero in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c_{tZ}$, using the photon pT distribution from the dilepton analysis. The value of $c^{I}_{tZ}$ is profiled in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c_{tZ}$, using the combination of photon pT distributions from the dilepton and lepton+jets analyses. The value of $c^{I}_{tZ}$ is profiled in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c^{I}_{tZ}$, using the photon pT distribution from the dilepton analysis. The value of $c_{tZ}$ is profiled in the fit.

Negative log-likelihood difference from the best-fit value for the one-dimensional scans of the Wilson coefficient $c^{I}_{tZ}$, using the combination of photon pT distributions from the dilepton and lepton+jets analyses. The value of $c_{tZ}$ is profiled in the fit.

Negative log-likelihood difference from the best-fit value as a function of Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$ from the interpretation of the dilepton measurement.

Negative log-likelihood difference from the best-fit value as a function of Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$ from the interpretation of the dilepton and lepton+jets measurements combined.

One-dimensional 68 and 95% CL intervals obtained for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$, using the photon $p_{T}$ distribution from the dilepton analysis, or the combination of photon pT distributions from the dilepton and lepton+jets analyses.

Comparison of observed $95\%$ CL intervals for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$. Results are shown from a CMS ttZ measurement [JHEP 03 (2020) 056], from a CMS ttZ & tZq interpretation [arXiv:2107.13896], from a CMS ttG (lepton+jets) measurement [arXiv:2107.01508], from this measurement, and from a global fit by J. Ellis et al. [JHEP 04 (2021) 279].

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

Summary of the one-dimensional intervals at 68 and 95% CL.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.

Expected limit contour for GGM scenario

Observed limit contour for GGM scenario

Expected limit contour for GGM scenario

Observed limit contour for GGM scenario

Expected limit contour for Neutralino BF scenario

Observed limit contour for Neutralino BF scenario

Expected limit contour for Chargino BF scenario

Observed limit contour for Chargino BF scenario

Expected limit contour for T5Wg

Observed limit contour for T5Wg

Expected limit contour for Gluino BF scenario

Observed limit contour for Gluino BF scenario

Covariance matrix for all bins (additional material)

Correlation matrix for all bins (additional material)