Level-1 muon trigger efficiency in cosmic-ray muon data (blue) and signal simulation (red) as a function of $d_0$, for the Level-1 trigger $p_T$ threshold used in the 2018 analysis triggers. The denominator in the efficiency calculation is the number of STA muons with $|\eta| < 1.2$ and $p_T > 28$ GeV.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA-to-TMS muon association procedure, as a function of true $L_{xy}$, in all simulated $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal samples combined. The fractions are computed relative to the number of signal events passing the trigger and containing two STA muons with more than 12 muon detector hits and $p_T > 10$ GeV matched to generated muons from $X \rightarrow \mu \mu$ decays.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA-to-TMS muon association procedure, as a function of true $L_{xy}$, in all simulated $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal samples combined. The fractions are computed relative to the number of signal events passing the trigger and containing two STA muons with more than 12 muon detector hits and $p_T > 10$ GeV matched to generated muons from $X \rightarrow \mu \mu$ decays.

Comparison of the number of events observed in 2016 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-STA dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2018 data in the STA-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 30 and 60 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legends also include the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{\mu \mu}$ intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $6 < min(d_0 / \sigma_{d_0}) \leq 10$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $10 < min(d_0 / \sigma_{d_0}) \leq 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, in representative $m_{\mu \mu}$ intervals in the $min(d_0 / \sigma_{d_0}) > 20$ bin. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility. The legend also includes the total number of observed events as well as the number of expected background events obtained inclusively, by applying the background evaluation method to the events in all $m_{Z_D}$ and min($d_0 / \sigma_{d_0}$) intervals combined.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

Comparison of the number of events observed in 2016 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

Comparison of the number of events observed in 2018 data in the TMS-TMS dimuon category with the expected number of background events, as a function of the smaller of the two $d_0 / \sigma_{d_0}$ values in the TMS-TMS dimuon. The black points with crosses show the number of observed events; the green and yellow components of the stacked histograms represent the estimated numbers of DY and QCD events, respectively. The last bin includes events in the overflow. The uncertainties in the total expected background (shaded area) are statistical only. Signal contributions expected from simulated $H \rightarrow Z_D Z_D$ with $m_{Z_D}$ of 20 and 50 GeV are shown in red and blue, respectively. Their yields are set to the corresponding combined median expected exclusion limits at 95% CL, scaled up as indicated in the legend to improve visibility.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(\Phi \rightarrow XX)B(X \rightarrow \mu \mu)$ as a function of $c\tau(X)$ in the heavy-scalar model, for $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 10\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 10\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 20\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 30\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 30\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 40\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 40\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 50\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 60\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

The 95% CL upper limits on $\sigma(H \rightarrow Z_DZ_D)B(Z_D \rightarrow \mu \mu)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m_{Z_D} = 60\ GeV$. The median expected limits obtained from the STA-STA, STA-TMS, and TMS-TMS dimuon categories are shown as dashed green, blue, and red curves, respectively; the combined median expected limits are shown as dashed black curves; the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits. The horizontal lines in gray correspond to the theoretical predictions for values of $B(H \rightarrow Z_DZ_D)$ indicated next to the lines.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $c\tau(Z_D)$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $c\tau(Z_D)$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $\epsilon$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Observed 95% CL exclusion contours in the HAHM model, in the ($m(Z_D)$, $\epsilon$) plane. The contours correspond to several representative values of $B(H \rightarrow Z_DZ_D$) ranging from 0.005 to 1%.

Background estimation and observed number of events in the STA-STA dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown. The mass interval is followed by the estimated and observed counts for the given year. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the STA-STA dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown, followed by the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ for the given year. The quoted uncertainties are statistical only.

Background estimation and observed number of events in the TMS-TMS dimuon category in 2016 data. The mass interval is followed by the estimated and observed counts within each $min(d_0 / \sigma_{d_0})$ bin in this mass interval. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the TMS-TMS dimuon category in 2016 data. For each mass interval, the table shows the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ in each of the three $\text{min}(d_0 / \sigma_{d_0})$ bins. The quoted uncertainties are statistical only

Background estimation and observed number of events in the TMS-TMS dimuon category in 2018 data. The mass interval is followed by the estimated and observed counts within each $min(d_0 / \sigma_{d_0})$ bin in this mass interval. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the TMS-TMS dimuon category in 2016 data. For each mass interval, the table shows the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ in each of the three $\text{min}(d_0 / \sigma_{d_0})$ bins. The quoted uncertainties are statistical only

Correspondence between the mass intervals in the TMS-TMS category and the parameters of the simulated signal samples.

Correspondence between the probed LLP masses and the chosen mass intervals in the TMS-TMS category.

Background estimation and observed number of events in the STA-TMS dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown. The mass interval is followed by the estimated and observed counts for the given year. The quoted uncertainties are statistical only.

Background estimations and observed numbers of events in the STA-TMS dimuon category in 2016 and 2018 data. For each probed LLP mass, the chosen mass interval is shown, followed by the predicted background yield $N^\text{est}_\text{bkg}$ and the observed number of events $N^\text{obs}$ for the given year. The quoted uncertainties are statistical only.

Number of events passing consecutive sets of selection criteria for 2018 collision data and the signal process $\Phi(125) \rightarrow XX(20\ GeV, c\tau = 13\ cm) \rightarrow \mu\mu$. Each row introduces a new criterion that is applied in addition to the selection of the previous row. In addition to the total number of events, N(events), the event yields of the individual dimuon vertex categories, STA-STA, TMS-TMS, and STA-TMS, are shown in separate columns for each data set. In these columns, events containing selected dimuons of different categories are independently counted for each category.

Number of events passing consecutive sets of selection criteria, in 2018 data and in a sample of simulated $\Phi \rightarrow XX \rightarrow \mu\mu$ signal events with $m(H) = 125\ GeV$, $m(X) = 20\ GeV$, and $c\tau = 13\ cm$. Each row introduces a new criterion that is applied in addition to the selection of the previous row. In addition to the total number of events $N(\text{total})$, the event yields in the individual dimuon categories, STA-STA, TMS-TMS, and STA-TMS, are shown in separate columns for each data set. In these columns, events containing selected dimuons of different categories are counted independently for each category.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 350\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 125\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 200\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 400\ GeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 150\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 150\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1000\ GeV$ and $m(X) = 350\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $\Phi \rightarrow XX \rightarrow 4\mu$ signal process with $m(\Phi) = 1\ TeV$ and $m(X) = 350\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 10\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 10\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 20\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 20\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 30\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 30\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 40\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 40\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 50\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 50\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 60\ GeV$. The figure shows efficiencies in the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well as the combined efficiency (black) calculated as the sum of the efficiencies of the individual categories. The signal efficiencies for the 2016 and 2018 datasets are shown as dashed and solid lines, respectively.

Overall signal efficiencies as a function of $c\tau$ for the $H \rightarrow Z_DZ_D \rightarrow \mu\mu + anything$ signal process with $m(H) = 125\ GeV$ and $m(Z_D) = 60\ GeV$. The plot shows efficiencies of the three dimuon categories, STA-STA (green), TMS-TMS (red), and STA-TMS (blue), as well the combined efficiency (black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper. The efficiencies in the 2016 and 2018 data sets are shown as dashed and solid curves, respectively.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true} < 20\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-STA dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true} < 20\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-STA dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true} < 20\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true} < 20\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true} < 20\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the TMS-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true} < 20\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the TMS-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-STA dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-STA dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the TMS-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $20\ cm < L_{xy}^\mathrm{true} < 70\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the TMS-TMS dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2016 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-STA dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper. Efficiencies for dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the STA-TMS and TMS-TMS dimuon categories are equal to zero.

Signal efficiencies as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the $\Phi \rightarrow XX \rightarrow \mu\mu + anything$ signal model, in 2018 samples. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria applied in the STA-STA dimuon category to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper. Efficiencies for dimuons with $70\ cm < L_{xy}^\mathrm{true} < 320\ cm$ in the STA-TMS and TMS-TMS dimuon categories are equal to zero.

Measured fiducial differential cross sections for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Systematic uncertainties for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Measured fiducial differential cross sections for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Systematic uncertainties for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Final state radiation correction used in the $\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-p_{T}^{Z}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Correlation matrix of statistical uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Measured total $Z$-boson cross-section for different datasets. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $y^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Systematic uncertainties in the single differential cross-sections in interval regions of $y^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Upper limits on the production cross section times branching fraction of the b* RH hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the b* VL hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+b hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+t hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for prompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Fraction of nonprompt $J/\psi$ mesons (in \%) in ($p_\text{T},y$) intervals. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-0.2$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-1$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=+1$ rather than zero, in ($p_\text{T},y$) intervals.

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.