95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of pT, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

95% CL interval on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Differential cross section (multiplied by the dimuon branching fraction) of prompt J/$\psi$ mesons in pPb collisions at backward $y_{\mathrm{CM}}$.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for prompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pp collisions.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for prompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pPb collisions.

Transverse momentum dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in seven $y_{\mathrm{CM}}$ ranges.

Rapidity dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 6.5<$p_{\mathrm{T}}$<10 GeV.

Rapidity dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 10<$p_{\mathrm{T}}$<30 GeV.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 0<$y_{\mathrm{CM}}$<0.9 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 0.9<$y_{\mathrm{CM}}$<1.5 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 1.5<$y_{\mathrm{CM}}$<1.93 range.

Dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons on the hadronic activity in the event, given by the transverse energy deposited in the CMS detector at large pseudorapidities $E_{\mathrm{T}}^{\mathrm{HF}|\eta|>4}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pp collisions at forward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pp collisions at backward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pPb collisions at forward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pPb collisions at backward $y_{\mathrm{CM}}$.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for nonprompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pp collisions.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for nonprompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pPb collisions.

Transverse momentum dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in seven $y_{\mathrm{CM}}$ ranges.

Rapidity dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 6.5<$p_{\mathrm{T}}$<10 GeV.

Rapidity dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 10<$p_{\mathrm{T}}$<30 GeV.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 0<$y_{\mathrm{CM}}$<0.9 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 0.9<$y_{\mathrm{CM}}$<1.5 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 1.5<$y_{\mathrm{CM}}$<1.93 range.

Dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons on the hadronic activity in the event, given by the transverse energy deposited in the CMS detector at large pseudorapidities $E_{\mathrm{T}}^{\mathrm{HF}|\eta|>4}$.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 4 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Observed signal strength ($\mu$) for a SM Higgs boson at 125.38 GeV. The labels untagged combined,dijet combined and combined represent the results obtained from simultaneous fits of the untagged categories, dijet categories, and full set of categories, respectively. The black solid line shows $\mu=1$, and the red dashed line shows the best fit value $\hat{\mu}=2.4$ of all categories combined.

Upper limit (95% CL) on the signal strength ($\mu$) relative to the SM prediction, as a function of the assumed value of the Higgs boson mass used in the fit.

Jet shape ratios ($\rho(\Delta r)_{\mathrm{PbPb}}/\rho(\Delta r)_{\mathrm{pp}}$) for inclusive and b jets as function of $\Delta r$ at $\sqrt{s_{NN}} = 5.02$ TeV.

Jet shape ratios ($\rho(\Delta r)_{\mathrm{PbPb}}/\rho(\Delta r)_{\mathrm{pp}}$) for inclusive and b jets as function of $\Delta r$ at $\sqrt{s_{NN}} = 5.02$ TeV.

Jet shape ratios ($\rho(\Delta r)_{\mathrm{PbPb}}/\rho(\Delta r)_{\mathrm{pp}}$) for inclusive and b jets as function of $\Delta r$ at $\sqrt{s_{NN}} = 5.02$ TeV.

Transverse momentum profile difference $P(\Delta r)_{\mathrm{PbPb}}-P(\Delta r)_{\mathrm{pp}}$ vs. $\Delta r$ for inclusive and b jets at $\sqrt{s_{NN}} = 5.02$ TeV.

Transverse momentum profile difference $P(\Delta r)_{\mathrm{PbPb}}-P(\Delta r)_{\mathrm{pp}}$ vs. $\Delta r$ for inclusive and b jets at $\sqrt{s_{NN}} = 5.02$ TeV.

Transverse momentum profile difference $P(\Delta r)_{\mathrm{PbPb}}-P(\Delta r)_{\mathrm{pp}}$ vs. $\Delta r$ for inclusive and b jets at $\sqrt{s_{NN}} = 5.02$ TeV.

Jet shape ratios ($\rho(\Delta r)_{\mathrm{b}}/\rho(\Delta r)_{\mathrm{incl.}}$) as function of $\Delta r$ from pp and PbPb collisions at $\sqrt{s_{NN}}= 5.02$ TeV.

Jet shape ratios ($\rho(\Delta r)_{\mathrm{b}}/\rho(\Delta r)_{\mathrm{incl.}}$) as function of $\Delta r$ from pp and PbPb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

Jet shape ratios ($\rho(\Delta r)_{\mathrm{b}}/\rho(\Delta r)_{\mathrm{incl.}}$) as function of $\Delta r$ from pp and PbPb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

Self-normalized zg distributions in PbPb and smeared pp collisions in the 30-50 centrality event class for jets with PTJET 160-180 GeV.

Self-normalized zg distributions in PbPb and smeared pp collisions in the 10-30 centrality event class for jets with PTJET 160-180 GeV.

Self-normalized zg distributions in PbPb and smeared pp collisions in the 0-10 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 50-80 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 30-50 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10-30 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 0-10 centrality event class for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 140-160 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 160-180 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 180-200 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 200-250 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 250-300 GeV.

Ratio of zg distributions in PbPb and smeared pp collisions in the 10% most central events for jets with PTJET 300-500 GeV.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 10-20% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 20-30% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 30-40% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 40-50% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ in 50-60% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 0-0.2% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 0-5% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 0-10% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 10-20% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 20-30% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 30-40% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 40-50% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ in 50-60% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ for multiplicity class $220 \leq N^{offline}_{trk}<260$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ for multiplicity class $185 \leq N^{offline}_{trk}<220$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ for multiplicity class $150 \leq N^{offline}_{trk}<185$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) elliptic flow, $v^{(\alpha)}_2$, as a function of $p_T$ $120 \leq N^{offline}_{trk}<150$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $220 \leq N^{offline}_{trk}<260$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $185 \leq N^{offline}_{trk}<220$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $150 \leq N^{offline}_{trk}<185$ in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) triangular flow, $v^{(\alpha)}_3$, as a function of $p_T$ for multiplicity class $120 \leq N^{offline}_{trk}<150$ in pPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 0-0.2% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 0-5% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 10-20% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 20-30% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 30-40% in PbPb collisions.

Pearson correlation coefficients $r_2$ and $r_3$ reconstructed using the leading and subleading flows as a function of $p^{a}_{T} - p^{b}_{T}$ in bin $p^{a}_{T}$ for centrality class 40-50% in PbPb collisions.

Elliptic and triangular ratios of subleading and leading flow as a function of multiplicity for the highest bin 2.5 < $p_T$ < 3.0 GeV in PbPb collisions.

Elliptic and triangular ratios of subleading and leading flow as a function of multiplicity for the highest bin 2.5 < $p_T$ < 3.0 GeV in pPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 0-0.2% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 0-5% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 0-10% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 10-20% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 20-30% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 30-40% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 40-50% centrality PbPb collisions.

Leading ($\alpha$ = 1) and subleading ($\alpha$ = 2) multiplicity modes, $v^{(\alpha)}_0$, as a function of $p_T$ in 50-60% centrality PbPb collisions.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to b quarks and $\tau$ leptons.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to light-flavor quarks and $\tau$ leptons.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 2.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 3.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 0.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 2.0$.

Observed and expected distributions of the jet multiplicity in the e+jet btag control region.

Observed and expected distributions of the jet multiplicity in the $\mu$+jet btag control region.

Measured min-$d_{xy}$ distribution in the 2-match category, after requiring the $d_{xy}$ muon to pass the isolation requirement $I_{\mathrm{rel}}^{\mathrm{PF}} <0.25$ (i. e., the B and D bins of the ABCD plane). Overlaid with a red histogram is the background predicted from the region of the ABCD plane failing the same requirement (the A and C bins), as well as three signal benchmark hypotheses (as defined in the legends), assuming $\alpha_D = \alpha_{\mathrm{EM}}$ (the fine-structure constant). The red hatched bands correspond to the background prediction uncertainty. The last bin includes the overflow.

Two-dimensional exclusion surfaces for $\Delta = 0.1 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Two-dimensional exclusion surfaces for $\Delta = 0.4 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).