Normalized differential distributions of the $\Delta R (\gamma,\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading top quark $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

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

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

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

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

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

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

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

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

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Observed differential dijet spectrum using full Run-2 data.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for quark-quark type dijet resonances.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for quark-gluon type dijet resonances.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for gluon-gluon type dijet resonances.

The observed and expected 95% CL upper limits on the universal quark coupling $g_{q}'$ as a function of resonance mass for a leptophobic Z' resonance that only couples to quarks.

Predictions of a leading order (LO) QCD simulation, normalized to an integrated luminosity of 138 fb$^{-1}$. The number of events are examined within bins of the four-jet mass and the ratio $\alpha$, which is the average dijet mass divided by the four-jet mass.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\overline{m}_{\mathrm{2j}}$ plane from a signal simulation of a diquark with a width of 1.5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\overline{m}_{\mathrm{2j}}$ plane from a signal simulation of a diquark with a width of 5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\overline{m}_{\mathrm{2j}}$ plane from a signal simulation of a diquark with a width of 10% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\alpha$ plane from a signal simulation of a diquark with a width of 1.5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\alpha$ plane from a signal simulation of a diquark with a width of 5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\alpha$ plane from a signal simulation of a diquark with a width of 10% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

Signal differential distributions as a function of four-jet mass for $\alpha_{\mathrm{true}}$ = 0.25, diquark masses of 2, 5, 8.6 TeV and various widths, for all $\alpha$ bins inclusively. The integral of each distribution has been normalized to unity.

The product of acceptance, $A$, and efficiency, $\varepsilon$, of a resonant signal with $\alpha_{\mathrm{true}}$ = 0.25 vs. the diquark mass for various diquark widths, and for all $\alpha$ bins inclusively. The acceptance is defined as the fraction of generated events passing the kinematic selection criteria, while the efficiency is the fraction of signal events satisfying $m_{\mathrm{4j}} > 1.6$ TeV. We also show the signal acceptance alone in the curves where $\varepsilon = 1$.

The four-jet mass distribution in data for 0.22 < $\alpha$ < 0.24, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV, and $\Gamma/M_{\mathrm{S}}$ = 1.5%, 10% are also shown, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.24 < $\alpha$ < 0.26, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV, and $\Gamma/M_{\mathrm{S}}$ = 1.5%, 10% are also shown, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.26 < $\alpha$ < 0.28, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.28 < $\alpha$ < 0.30, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.30 < $\alpha$ < 0.32, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.32 < $\alpha$ < 0.34, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV, and $\Gamma/M_{\mathrm{S}}$ = 1.5%, 10% are also shown, with cross sections equal to the observed upper limits at 95% confidence level.

The inclusive four-jet mass distribution in data for $\alpha$ > 0.10, fitted with three background-only functions (Dijet-5p, PowExp-5p and ModDijet-5p), each with five free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.25$ and width of the initial resonance Y equal to 1.5%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 1.5%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.25$ and width of the initial resonance Y equal to 5%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 5%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.25$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.11$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.13$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.15$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.17$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.19$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.21$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.23$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.27$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.29$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.31$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.33$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.42$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.11$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.13$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.15$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.17$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.19$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.21$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.23$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.27$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.29$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.31$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.33$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.42$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

Observed local $p$-value for a four-jet resonance, Y, decaying to a pair of dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}}/M_{\mathrm{Y}}$ = 0.25, and various widths of Y superimposed.

Observed local $p$-value for a four-jet resonance, Y, decaying to a pair of dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}}/M_{\mathrm{Y}}$ = 0.29, and various widths of Y superimposed.

The table presents the cumulative cutflow for a signal with $M_{\mathrm{S}} = 2.0$ TeV, $M_{\chi} = 0.5$ TeV and different width hypotheses ($\Gamma/M_{\mathrm{S}} =$ 1.5, 5, and 10%). Each row corresponds to the number of signal events that survive all cuts up to and including the one listed. The percentage in parentheses shows the efficiency of the current cut alone, defined as the ratio of the number of events that survive all cuts up to and including this one to the number of events that survived all previous cuts.

The table presents the cumulative cutflow for a signal with $M_{\mathrm{S}} = 5.0$ TeV, $M_{\chi} = 1.25$ TeV and different width hypotheses ($\Gamma/M_{\mathrm{S}} =$ 1.5, 5, and 10%). Each row corresponds to the number of signal events that survive all cuts up to and including the one listed. The percentage in parentheses shows the efficiency of the current cut alone, defined as the ratio of the number of events that survive all cuts up to and including this one to the number of events that survived all previous cuts.

The table presents the cumulative cutflow for a signal with $M_{\mathrm{S}} = 8.6$ TeV, $M_{\chi} = 2.15$ TeV and different width hypotheses ($\Gamma/M_{\mathrm{S}} =$ 1.5, 5, and 10%). Each row corresponds to the number of signal events that survive all cuts up to and including the one listed. The percentage in parentheses shows the efficiency of the current cut alone, defined as the ratio of the number of events that survive all cuts up to and including this one to the number of events that survived all previous cuts.

Differential cross-section in bins of subleading muon $p_\mathrm{T]$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading muon $\eta$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading muon $\eta$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading muon $\phi$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading muon $\phi$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet $p_\mathrm{T]$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet $p_\mathrm{T]$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet rapidity. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet rapidity. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet azimuth. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet azimuth. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet mass. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet mass. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet constituent multiplicity. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet constituent multiplicity. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet $\tau_1$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet $\tau_1$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet $\tau_2$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet $\tau_2$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet $\tau_3$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of subleading charged particle jet $\tau_3$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of leading charged particle jet $\tau_{21}$. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

Differential cross-section in bins of $\Delta R$ between the leading charged particle jet and the dilepton system. The actual measurement is unbinned and available with examples at <a href="https://gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024">gitlab.cern.ch/atlas-physics/public/sm-z-jets-omnifold-2024</a>

The total uncertainty covariance matrix for measured differential cross sections in $m_{4\ell}$ bins in fiducial volume. The matrix elements are in unit of $[$fb/TeV$]^2$.

The total uncertainty covariance matrix for measured differential cross sections in $p_\mathrm{T}^{4\ell}$ bins in fiducial volume. The matrix elements are in unit of $[$fb/TeV$]^2$.

Measured differential four-jet cross section for R=0.4 jets, in bins of pT4, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of HT, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>500 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>1500 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>2000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar transverse momentum sum of all jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of leading two jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between diphoton and dijet systems. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between the two leading jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton thrust pT. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of the two photons. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of tau of highest-tau jet (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar sum of tau for all jets (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of third leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of dijet invariant mass. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the combined diphoton and dijet system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in jet multiplicity bins. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of leading jet transverse momentum for exclusive one-jet events. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Isolation efficiencies in percent for gluon fusion H -> 2gam for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see table~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Statistical bin-to-bin correlations between the 24 bins used in the analysis described in Phys.Lett. B753 (2016) 69-85 (arXiv:1508.02507). See Figure 3 of this paper for the bins used.

Breakdown of systematic uncertainties (in %) for the inclusive $b$-jet cross-section $\sigma(Zb)\times N_{b\text{-jet}}$ as a function of $b$-jet $p_T$.

The inclusive $b$-jet cross-section $\sigma(Zb)\times N_{b\text{-jet}}$ as a function of $b$-jet $|y|$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the inclusive $b$-jet cross-section $\sigma(Zb)\times N_{b\text{-jet}}$ as a function of $b$-jet $|y|$.

The inclusive $b$-jet cross-section $\sigma(Zb)\times N_{b\text{-jet}}$ as a function of $y_{\text{boost}}(Z,b)$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the inclusive $b$-jet cross-section $\sigma(Zb)\times N_{b\text{-jet}}$ as a function of $y_{\text{boost}}(Z,b)$.

The inclusive $b$-jet cross-section $\sigma^{*}(Zb)\times N_{b\text{-jet}}$ as a function of $\Delta y(Z,b)$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the inclusive $b$-jet cross-section $\sigma^{*}(Zb)\times N_{b\text{-jet}}$ as a function of $\Delta y(Z,b)$.

The inclusive $b$-jet cross-section $\sigma^{*}(Zb)\times N_{b\text{-jet}}$ as a function of $\Delta\phi(Z,b)$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the inclusive $b$-jet cross-section $\sigma^{*}(Zb)\times N_{b\text{-jet}}$ as a function of $\Delta\phi(Z,b)$.

The inclusive $b$-jet cross-section $\sigma^{*}(Zb)\times N_{b\text{-jet}}$ as a function of $\Delta R(Z,b)$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the inclusive $b$-jet cross-section $\sigma^{*}(Zb)\times N_{b\text{-jet}}$ as a function of $\Delta R(Z,b)$.

The cross-section $\sigma(Zb)$ as a function of $Z$ boson $p_T$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the cross-section $\sigma(Zb)$ as a function of $Z$ boson $p_T$.

The cross-section $\sigma(Zb)$ as a function of $Z$ boson $|y|$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the cross-section $\sigma(Zb)$ as a function of $Z$ boson $|y|$.

Integrated $Z+\ge 2$ $b$-jets cross-section.

Breakdown of systematic uncertainties (in %) for the integrated $Z+\ge 2$ $b$-jets cross-section.

The cross-section $\sigma(Zbb)$ as a function of $\Delta R(b,b)$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the cross-section $\sigma(Zbb)$ as a function of $\Delta R(b,b)$.

The cross-section $\sigma(Zbb)$ as a function of $\Delta m(b,b)$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the cross-section $\sigma(Zbb)$ as a function of $\Delta m(b,b)$.

The cross-section $\sigma(Zbb)$ as a function of $Z$ boson $p_T$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the cross-section $\sigma(Zbb)$ as a function of $Z$ boson $p_T$.

The cross-section $\sigma(Zbb)$ as a function of $Z$ boson $|y|$ together with the corresponding non-perturbative corrections.

Breakdown of systematic uncertainties (in %) for the cross-section $\sigma(Zbb)$ as a function of $Z$ boson $|y|$.

Normalized differential cross-sections for the absolute rapidity (|y|) of the ttbar system. The cross-section in each bin is given as the integral of the normalized differential cross-section over the bin width, divided by the bin width. The calculation of the cross-sections in the last bins includes events falling outside of the bin edges, and the normalization is done within the quoted bin width. The full covariance matrice is provided in Table 8 below.

Bin-wise full covariance matrices for the normalized differential cross-sections of the hadronically decaying top-quark PT. The elements of the covariance matrices are in units of TeV$^{-2}$.

Bin-wise full covariance matrices for the normalized differential cross-sections of the mass of the ttbar system. The elements of the covariance matrices are in units of TeV$^{-2}$.

Bin-wise full covariance matrices for the normalized differential cross-sections of the PT of the ttbar system. The elements of the covariance matrices are in units of TeV$^{-2}$.

Bin-wise full covariance matrices for the normalized differential cross-sections of the absolute rapidity (|y|) of the ttbar system.

The systematic uncertainties for the nominal muon-channel cross-section measurement. Some sources of uncertainty have both correlated and uncorrelated components. Correlated uncertainties arise from the uncertainty in the electroweak background contributions delta(e.w.)_cor, from corrections to the Monte Carlo modelling of the Z/gamma* pT spectra, and from the reconstruction, trigger and isolation efficiency corrections, given by delta(reco)_cor, delta(trig)_cor and delta(iso)_cor respectively. The uncertainty on the multijet and electroweak background cross sections, given by delta(multijet and delta(e.w.), respectively and muon momentum scale uncertainty, delta(pT_scale) are also correlated. Uncorrelated uncertainties are due to corrections for the trigger and isolation efficiencies, given by delta(trig)_unc and delta(iso)_unc respectively. The uncertainty from the muon resolution correction, delta(res), from the size of the signal Monte Carlo sample, delta(MC), and the uncertaintities due to corrections for the reconstruction, delta(reco)_unc, are also uncorrelated. The luminosity uncertainty 1.8% is not included.

The combined Born-level fiducial differential cross section dsigma/dm_ll, statistical delta(stat), total correlated delta(cor), uncorrelated delta(unc), and total delta(total) uncertainties, as well as individual correlated sources delta(cor)_i. The correlated uncertainties are a linear combination of the 13 correlated uncertainties in the nominal muon and electron channels. As the uncertainties on the combined result no longer originate from individual error sources they are numbered 1-13. Also shown is the correction factor used to derive the dressed cross section (D), and the NNLO extrapolation factor (A) used to derive the cross section for the full phase space, along with the uncertainties associated to variations in scale choice delta(A)_scale, and PDF uncertainty delta(A)_pdf_alpha_s. The luminosity uncertainty 1.8% is not included.

The extended muon channel Born-level fiducial differential cross section dsigma/dm_ll, with the statistical delta(stat), systematic delta(syst), and total delta(tot) uncertainties for each invariant mass bin. Also shown is the correction factor used to derive the dressed cross section (D), and the extrapolation factor (A) used to derive the cross section for the full phase space, along with the uncertainties associated to variations in scale delta(A)_scale, and PDF uncertainty delta(A)_pdf_alpha_s. The luminosity uncertainty 3.5% is not included.

The systematic uncertainties for the extended muon channel cross-section measurement in each invariant mass bin. Correlated uncertainties come from the reconstruction, trigger and isolation efficiency corrections, given by delta(reco), delta(trig) and delta(iso) respectively. The uncertainty on the multijet background cross section, delta(multijet) and the uncertainty on the muon momentum scale, delta(pT_scale), are also correlated across bins. Uncorrelated uncertainties are due to the uncertainty from the muon resolution correction, delta(res), and the sample size of the signal Monte Carlo sample, delta(MC). The luminosity uncertainty 3.5% is not included.

Theory predictions for NLO+LLPS and for fixed-order calculations at NLO and NNLO including higher-order electroweak corrections, for the nominal analysis of the differential cross section dsigma/dm_ll as a function of the invariant mass m_ll. The scale uncertainty is defined as the envelope of variations for 0.5< mu_R,mu_F<2.0 for Powheg. For Fewz the scale uncertainty is defined by the variation 0.5<mu_R=mu_F<2.0. The uncertainty labelled as "stat" in the table corresponds to the pdf+alpha_s uncertainty for Fewz or just the pdf uncertainty for Powheg.

Theory predictions for NLO+LLPS and for fixed-order calculations at NLO and NNLO including higher-order electroweak corrections, for the extended analysis of the differential cross section dsigma/dm_ll as a function of the invariant mass m_ll. The scale uncertainty is defined as the envelope of variations for 0.5<mu_R,mu_F<2.0 for Powheg. For Fewz the scale uncertainty is defined by the variation 0.5<mu_R=mu_F<2.0. The uncertainty labelled as "stat" in the table corresponds to the pdf+alpha_s uncertainty for Fewz or just the pdf uncertainty for Powheg.

Higher-order electroweak corrections in nominal analysis, Delta(HOEW), and the correction for the Photon Induced process, Delta(PI), together with its uncertainty derived from the uncertainty of the photon PDF as a function of the dilepton invariant mass m_ll. Also shown is the difference arising from the non-convergence of calculations derived with different electroweak schema, delta(scheme).

Higher-order electroweak corrections in extended analysis, Delta(HOEW), and the correction for the Photon Induced process, Delta(PI), together with its uncertainty derived from the uncertainty of the photon PDF as a function of the dilepton invariant mass m_ll. Also shown is the difference arising from the non-convergence of calculations derived with different electroweak schema, delta(scheme).

Values of chi^2 for the nominal and extended DY cross-section measurements for predictions based from Fewz at NLO and NNLO and Powheg using MSTW2008 PDFs accounting for the experimental correlated systematic uncertainties. The _pdf_scale columns show the chi^2 when the PDF and scale uncertainties on the theoretical predictions are taken into account.

The chi^2 values from NLO and NNLO QCD fits made using the ATLAS nominal and extended Drell--Yan cross-section measurements and the HERA-I dataset accounting for the experimental correlated systematic uncertainties.