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.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons with a Soeding-inspired analytic function --- statistical correlations only. Only one half of the (symmetric) matrix is stored.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic $\rho^0(770)$ photoproduction cross section off protons as a function of energy (various experiments)

Fit of elastic $\rho^0(770)$ photoproduction cross sections off protons as a function of energy (various experiments) --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single differential in the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross section off protons, single differential in the momentum transfer $\vert t\vert$ and in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Differential fiducial cross section of the Z boson $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the leading b-jet $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the Z boson $p_{\text{T}}$ in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $p_{\text{T}}$ of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the ratio between the $p_{\text{T}}$ and the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

The measured fiducial cross section vs $m(\ell\ell\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $p_{\mathrm{T}}^{\ell\ell\gamma}$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $\Delta\phi(ll,\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $p_{\mathrm{T}}^{\ell\ell\gamma}/m(\ell\ell\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

Measured differential cross section with associated uncertainties as a function of COS(THETA*). 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 DELTAYRAP(2GAMMA). 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 MULT(JET,PT>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 MULT(JET,PT>50 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 PT(JET1). 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 PT(JET2). 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 HT. 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 YRAP(JET1). 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 YRAP(JET2). 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 M(2JET). 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 DELTAYRAP(2JET). 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 ABSDPHI(2JET). 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 DPHI(2JET). 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 PT(2GAMMA2JET). 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 DPHI(2GAMMA,2JET). 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 TAUJET. 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 SUM(TAUJET). 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 PT(2GAMMA) [NJET=0,PT>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 PT(2GAMMA) [NJET=1,PT>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 PT(2GAMMA) [NJET=2,PT>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 PT(2GAMMA) [NJET>=3,PT>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.

The measured cross sections or cross section limits of the diphoton, VBF-enhanced, Nlepton $\geq$ 1, high $E_{T}^{miss}$, and ttH-enhanced fiducial regions are shown.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in bins of ABS(COS(THETA*)). 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.Each systematic uncertainty sources is fully uncorrelated with the other sources.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text. No factor are given for the Nlepton $\geq$ 1 and High-$E_{T}^{miss}$ fiducial regions as the gluon fusion contamination in both is negligible.

Isolation efficiencies in percent for gluon fusion $H\rightarrow\gamma\gamma$ 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 (as described in Section 9.1) for events that pass the diphoton kinematic criteria. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. 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 29) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis. No factors for the Nlepton $\geq$ 1 and High $E_{T}^{miss}$ fiducial regions are provided as the gluon fusion contamination is negligible.

Combined non-perturbative (Table 27) and particle-level isolation correction factors (Table 28) in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Regions of phase space where no reliable estimate could be obtained are listed as 100 without uncertainties. The uncertainties on the combined values properly take into account the correlations between both multiplicative factors.

Diphoton kinematic acceptances in percent for gluon-gluon fusion for the diphoton fiducial region and all differential variable bins studied in this paper, defined as the probability to fulfill the diphoton kinematic criteria: $p_{T}$/$m_{\gamma\gamma}$ < 0.35 (0.25) for the leading (subleading) photon and $|\eta_{\gamma\gamma}|$ < 2.37. The factors are evaluated using the Powheg NNLOPSevent generator. Uncertainties are taken from PDF variations. QCD scale variations have a negligible impact on these factors. The range of each bin is given in Table 26.

observed statistical correlations between pTyy, Njets, mjj, |DeltaPhijj| and pTj1

ggH default MC + XH predictions

XH ( = VBF + VH + ttH + bbH ) MC predictions

Best-fit values and uncertainties of the production-mode cross sections times branching ratio.

Best-fit values and uncertainties of the simplified template cross sections times branching ratio.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Best-fit values and uncertainties of the simplified template cross sections times branching ratio.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Observed correlations between the measured simplified template cross sections, including both the statistical and systematic uncertainties.

Measured cross sections for isolated-photon plus jet production as a function of $m^{\gamma-\rm jet}$.

Measured cross sections for isolated-photon plus jet production as a function of $|\cos\theta^{\star}|$.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $467 < m^{\gamma-\rm jet1} < 510$ GeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $510 < m^{\gamma-\rm jet1} < 590$ GeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $590 < m^{\gamma-\rm jet1} < 670$ GeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $670 < m^{\gamma-\rm jet1} < 760$ GeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $760 < m^{\gamma-\rm jet1} < 850$ GeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $850 < m^{\gamma-\rm jet1} < 950$ GeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $0.95 < m^{\gamma-\rm jet1} < 1.05$ TeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $1.05 < m^{\gamma-\rm jet1} < 1.25$ TeV.

Measured cross sections for isolated-photon plus 1jet production as a function of $|\cos\theta^{\star}|$ for the range $1.25 < m^{\gamma-\rm jet1} < 2.45$ TeV.

Measured cross sections for isolated-photon plus 2jet production as a function of $E_{\rm T}^{\gamma}$.

Measured cross sections for isolated-photon plus 2jet production as a function of $p_{\rm T}^{\rm jet2}$.

Measured cross sections for isolated-photon plus 2jet production as a function of $\Delta\phi^{\rm \gamma-jet2}$.

Measured cross sections for isolated-photon plus 2jet production as a function of $\Delta\phi^{\rm jet1-jet2}$.

Measured cross sections for isolated-photon plus 2jet production as a function of $\Delta\phi^{\rm \gamma-jet2}$ for the range $130< E_{\rm T}^{\gamma} <300$ GeV.

Measured cross sections for isolated-photon plus 2jet production as a function of $\Delta\phi^{\rm \gamma-jet2}$ for the range $E_{\rm T}^{\gamma} >300$ GeV.

Measured cross sections for isolated-photon plus 2jet production as a function of $\Delta\phi^{\rm jet1-jet2}$ for the range $130< E_{\rm T}^{\gamma} <300$ GeV.

Measured cross sections for isolated-photon plus 2jet production as a function of $\Delta\phi^{\rm jet1-jet2}$ for the range $E_{\rm T}^{\gamma} >300$ GeV.

Measured cross sections for isolated-photon plus 2jet production as a function of $\beta^{\rm jet1}$.

Measured cross sections for isolated-photon plus 2jet production as a function of $\beta^{\gamma}$.

Ratio of measured cross sections for isolated-photon plus 2jet production as a function of $\beta$.

Measured cross sections for isolated-photon plus 3jet production as a function of $E_{\rm T}^{\gamma}$.

Measured cross sections for isolated-photon plus 3jet production as a function of $p_{\rm T}^{\rm jet3}$.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm \gamma-jet3}$.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm jet1-jet3}$.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm jet2-jet3}$.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm \gamma-jet3}$ for the range $130< E_{\rm T}^{\gamma} <300$ GeV.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm \gamma-jet3}$ for the range $E_{\rm T}^{\gamma} >300$ GeV.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm jet1-jet3}$ for the range $130< E_{\rm T}^{\gamma} <300$ GeV.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm jet1-jet3}$ for the range $E_{\rm T}^{\gamma} >300$ GeV.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm jet2-jet3}$ for the range $130< E_{\rm T}^{\gamma} <300$ GeV.

Measured cross sections for isolated-photon plus 3jet production as a function of $\Delta\phi^{\rm jet2-jet3}$ for the range $E_{\rm T}^{\gamma} >300$ GeV.

Parton-level normalized $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 8 TeV. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level normalized $t\bar t$ differential cross-sections for the $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 8 TeV. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level normalized $t\bar t$ differential cross-sections for the $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 8 TeV. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level absolute $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 7 TeV. The cross-sections in the last bins include events (if any) beyond of the bin edges. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level absolute $t\bar t$ differential cross-sections for the $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 7 TeV. The cross-sections in the last bins include events (if any) beyond of the bin edges. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level absolute $t\bar t$ differential cross-sections for the $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 7 TeV. The cross-sections in the last bins include events (if any) beyond of the bin edges. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level absolute $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 8 TeV. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level absolute $t\bar t$ differential cross-sections for the $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 8 TeV. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Parton-level absolute $t\bar t$ differential cross-sections for the $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 8 TeV. The uncertainties quoted in the second column represent the statistical and systematic uncertainties added in quadrature.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 7 TeV. The elements of the covariance matrix are in units of 10$^{-6}$ GeV$^{-2}$.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 7 TeV. The elements of the covariance matrix are in units of 10$^{-6}$ GeV$^{-2}$.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 7 TeV. The elements of the covariance matrix are unit-less.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 8 TeV. The elements of the covariance matrix are in units of 10$^{-6}$ GeV$^{-2}$.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 8 TeV. The elements of the covariance matrix are in units of 10$^{-6}$ GeV$^{-2}$.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 8 TeV. The elements of the covariance matrix are unit-less.

Full covariance matrix of the absolute $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 7 TeV. The elements of the covariance matrix are in units of pb$^2$ GeV$^{-2}$.

Full covariance matrix of the absolute $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 7 TeV. The elements of the covariance matrix are in units of pb$^2$ GeV$^{-2}$.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 7 TeV. The elements of the covariance matrix are pb$^2$.

Full covariance matrix of the absolute $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 8 TeV. The elements of the covariance matrix are in units of pb$^2$ GeV$^{-2}$.

Full covariance matrix of the absolute $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 8 TeV. The elements of the covariance matrix are in units of pb$^2$ GeV$^{-2}$.

Full covariance matrix of the normalized $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 8 TeV. The elements of the covariance matrix are pb$^2$.

Statistical bin-to-bin correlations in the normalized $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 7 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the normalized $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 7 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the normalized $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 7 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the normalized $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 8 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the normalized $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 8 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the normalized $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 8 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the absolute $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 7 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the absolute $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 7 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the absolute $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 7 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the absolute $t\bar t$ differential cross-sections for $t\bar t$ system mass $m_{t\bar t}$ at $\sqrt{s}$ = 8 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the absolute $t\bar t$ differential cross-sections for $t\bar t$ system transverse momentum $p_{T, t\bar t}$ at $\sqrt{s}$ = 8 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.

Statistical bin-to-bin correlations in the absolute $t\bar t$ differential cross-sections for $t\bar t$ system absolute rapidity $|y_{t\bar t}|$ at $\sqrt{s}$ = 8 TeV. The off-diagonal correlations mainly due to bin migrations in unfolding and the normalization condition.