Product of acceptance and efficiency for dimuon pairs as a function of generated mass in simulated events. The DY samples are used to represent spin-1 particles, and RS graviton samples are used for spin-2 particles.

The invariant mass spectra of dielectron events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass spectra of dielectron events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass spectra of dimuon events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass spectra of dimuon events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

Observed and expected background yields for different mass ranges in the dielectron channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dielectron channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the GSM model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the GSM model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the LR model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the LR model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the E6 model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the E6 model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

The observed local p-value for the dielectron channel as a function of the dilepton invariant mass.

The observed local p-value for the dielectron channel as a function of the dilepton invariant mass.

The observed local p-value for the dimuon channel as a function of the dilepton invariant mass.

The observed local p-value for the dimuon channel as a function of the dilepton invariant mass.

The observed local p-value for the combined channel as a function of the dilepton invariant mass.

The observed local p-value for the combined channel as a function of the dilepton invariant mass.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dielectron channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dielectron channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the combination of the dielectron and dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the combination of the dielectron and dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dielectron channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dielectron channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dimuon channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dimuon channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the combination of the dielectron and dimuon channels. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the combination of the dielectron and dimuon channels. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with central leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with central leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with forward leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with forward leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Invariant mass resolution for dielectron pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dielectron pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dimuon pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dimuon pairs as obtained from simulated DY samples as a function of generated mass.

Observed data and total background estimate in all bins used in the limit calculation for the CI limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 6 bins in dilepton invariant mass: [400,500,700,1100,1900,3500,10000] in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the CI limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 6 bins in dilepton invariant mass: [400,500,700,1100,1900,3500,10000] in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the ADD limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 5 bins in dilepton invariant mass: [1800,2200,2600,3000,3400,10000] ([1900,2200,2600,3000,3400,10000] for 2016) in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the ADD limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 5 bins in dilepton invariant mass: [1800,2200,2600,3000,3400,10000] ([1900,2200,2600,3000,3400,10000] for 2016) in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the CI limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the CI limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the ADD limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the ADD limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Dilepton exclusion limits at 95% CL on the CI scale (Lambda) for the six CI models considered for the muon channel. The limits are obtained for m > 400(2200) GeV in the case of constructive (destructive) interference.

Combined dilepton 95% CL exclusion limits on the cross section for the left-left constructive CI model. The red curve in the left plot shows the theoretical cross section as a function of the CI Scale Lambda.

Combined dilepton 95% CL exclusion limits on the CI scale (Lambda) for the six different CI models considered. Thelimits are obtained for m > 400(2200) GeV in the case of constructive (destructive) interference.

Exclusion limits at 95% CL on the UV cutoff for the electron channel with m> 1.8 TeV in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to leading order and a correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the UV cutoff for the muon channel with m> 1.8 TeV in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to leading order and a correction factor of 1.3 is applied.

Combined dilepton 95% CL exclusion limit on the cross section in the GRW convention with m > 1.8 TeV for the ADD model.

Combined dilepton 95% CL exclusion limits on the UV cutoff with m> 1.8 TeV in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to leading order and a correction factor of 1.3 is applied.

Individual and combined dilepton (this analysis) and diphoton 95% CL expected exclusion limits as a summary of all parameter conventions for the ADD model. Signal model cross sections are calculated up to leading order.

Individual and combined dilepton (this analysis) and diphoton 95% CL observed exclusion limits as a summary of all parameter conventions for the ADD model. Signal model cross sections are calculated up to leading order.