CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-HIG-21-001
Searches for additional Higgs bosons and vector leptoquarks in $\tau\tau$ final states in proton-proton collisions at $\sqrt{s}= $ 13 TeV
Abstract: Three searches are presented for signatures of physics beyond the standard model (SM) in $\tau\tau$ final states in proton-proton collisions at the CERN LHC, using a data sample collected with the CMS detector at a centre-of-mass energy of $\sqrt{s} = $ 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Upper limits at 95% confidence level (CL) are set on the products of the branching fraction for the decay into $\tau$ leptons and the cross sections for the production of a resonance $\phi$ in addition to the observed Higgs boson via gluon fusion ($\mathrm{gg}\phi$) or in association with b quarks, ranging from $\mathcal{O}$(10 pb) (at 60 GeV) to 0.3 fb (at 3.5 TeV) each. The data reveal two excesses for $\mathrm{gg}\phi$ production with local $p$-values equivalent to about three standard deviations at 0.1 and 1.2 TeV. In a search for $t$-channel exchange of a vector-like leptoquark $\mathrm{U}_1$, 95% CL upper limits are set on the $\mathrm{U}_1$ coupling to quarks and $\tau$ leptons ranging from 1 (at 1 TeV) to 6 (at 5 TeV), depending on the scenario. In the interpretation of minimal supersymmetric SM (MSSM) benchmark scenarios, additional Higgs bosons with masses below 350 GeV are excluded at 95% CL.
Figures & Tables Summary Additional Figures References CMS Publications
Figures

png pdf
Figure 1:
Diagrams for the production of neutral Higgs bosons (left) via gluon fusion (${\mathrm{g} \mathrm{g} \phi}$) and (middle and right) in association with b quarks (${\mathrm{b} \mathrm{b} \phi}$). In the middle panel a pair of b quarks is produced from two gluons. In the right panel ${\phi}$ is radiated from a b quark in the proton.

png pdf
Figure 1-a:
Diagrams for the production of neutral Higgs bosons (left) via gluon fusion (${\mathrm{g} \mathrm{g} \phi}$) and (middle and right) in association with b quarks (${\mathrm{b} \mathrm{b} \phi}$). In the middle panel a pair of b quarks is produced from two gluons. In the right panel ${\phi}$ is radiated from a b quark in the proton.

png pdf
Figure 1-b:
Diagrams for the production of neutral Higgs bosons (left) via gluon fusion (${\mathrm{g} \mathrm{g} \phi}$) and (middle and right) in association with b quarks (${\mathrm{b} \mathrm{b} \phi}$). In the middle panel a pair of b quarks is produced from two gluons. In the right panel ${\phi}$ is radiated from a b quark in the proton.

png pdf
Figure 1-c:
Diagrams for the production of neutral Higgs bosons (left) via gluon fusion (${\mathrm{g} \mathrm{g} \phi}$) and (middle and right) in association with b quarks (${\mathrm{b} \mathrm{b} \phi}$). In the middle panel a pair of b quarks is produced from two gluons. In the right panel ${\phi}$ is radiated from a b quark in the proton.

png pdf
Figure 2:
Diagram for the production of a pair of $\tau$ leptons via the $t$-channel exchange of a leptoquark ${\mathrm {U}_1}$.

png pdf
Figure 3:
Observed and expected distributions of (left) ${D_{\zeta}}$ in the e$\mu$ final state and (right) ${m_{\text {T}}^{\mu}}$ in the $\mu {\tau _\mathrm {h}}$ final state. The dashed vertical lines indicate the high-mass category definitions in each of the final states. A detailed discussion of the data modelling is given in Section 6. The distributions are shown in the "no b-tag'' category before any further event categorization and after a fit to the data in each corresponding variable. The grey shaded band represents the complete set of uncertainties used for signal extraction, after the fit.

png pdf
Figure 3-a:
Observed and expected distributions of (left) ${D_{\zeta}}$ in the e$\mu$ final state and (right) ${m_{\text {T}}^{\mu}}$ in the $\mu {\tau _\mathrm {h}}$ final state. The dashed vertical lines indicate the high-mass category definitions in each of the final states. A detailed discussion of the data modelling is given in Section 6. The distributions are shown in the "no b-tag'' category before any further event categorization and after a fit to the data in each corresponding variable. The grey shaded band represents the complete set of uncertainties used for signal extraction, after the fit.

png pdf
Figure 3-b:
Observed and expected distributions of (left) ${D_{\zeta}}$ in the e$\mu$ final state and (right) ${m_{\text {T}}^{\mu}}$ in the $\mu {\tau _\mathrm {h}}$ final state. The dashed vertical lines indicate the high-mass category definitions in each of the final states. A detailed discussion of the data modelling is given in Section 6. The distributions are shown in the "no b-tag'' category before any further event categorization and after a fit to the data in each corresponding variable. The grey shaded band represents the complete set of uncertainties used for signal extraction, after the fit.

png pdf
Figure 4:
Overview of the high-mass categories used for the extraction of the signal for the model-independent ${\phi}$ search for $ {m_{\phi}} \geq $ 250 GeV, and the vector-like leptoquark search and for the interpretation of the data in MSSM benchmark scenarios.

png pdf
Figure 5:
Overview of the low-mass categories used for the extraction of the signal for the model-independent ${\phi}$ search for 60 $\leq {m_{\phi}} < $ 250 GeV.

png pdf
Figure 6:
Composition of the differential signal for the MSSM interpretation of the data and the vector-like leptoquark search. On the left the A ${p_{\mathrm {T}}}$ density for the MSSM ${M_{\text {h}}^{125}}$ scenario for $ {m_{{\mathrm {A}}}} = $ 1.6 TeV and $ {\tan\beta} =30$ is shown, split by the contributions from the t quark only, the b quark only and the $\mathrm{t} \mathrm{b} $-interference term. On the right the ${m_{\text {T}}^{\text {tot}}}$ distribution in the ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final state, which is the most sensitive final state for the ${\mathrm {U}_1}$ search, is shown for ${\mathrm {U}_1}$ $t$-channel exchange with $ {m_{\text {U}}} = $ 1 TeV and $ {g_{\text {U}}} =1.5$, for the signal with and without the interference term.

png pdf
Figure 6-a:
Composition of the differential signal for the MSSM interpretation of the data and the vector-like leptoquark search. On the left the A ${p_{\mathrm {T}}}$ density for the MSSM ${M_{\text {h}}^{125}}$ scenario for $ {m_{{\mathrm {A}}}} = $ 1.6 TeV and $ {\tan\beta} =30$ is shown, split by the contributions from the t quark only, the b quark only and the $\mathrm{t} \mathrm{b} $-interference term. On the right the ${m_{\text {T}}^{\text {tot}}}$ distribution in the ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final state, which is the most sensitive final state for the ${\mathrm {U}_1}$ search, is shown for ${\mathrm {U}_1}$ $t$-channel exchange with $ {m_{\text {U}}} = $ 1 TeV and $ {g_{\text {U}}} =1.5$, for the signal with and without the interference term.

png pdf
Figure 6-b:
Composition of the differential signal for the MSSM interpretation of the data and the vector-like leptoquark search. On the left the A ${p_{\mathrm {T}}}$ density for the MSSM ${M_{\text {h}}^{125}}$ scenario for $ {m_{{\mathrm {A}}}} = $ 1.6 TeV and $ {\tan\beta} =30$ is shown, split by the contributions from the t quark only, the b quark only and the $\mathrm{t} \mathrm{b} $-interference term. On the right the ${m_{\text {T}}^{\text {tot}}}$ distribution in the ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final state, which is the most sensitive final state for the ${\mathrm {U}_1}$ search, is shown for ${\mathrm {U}_1}$ $t$-channel exchange with $ {m_{\text {U}}} = $ 1 TeV and $ {g_{\text {U}}} =1.5$, for the signal with and without the interference term.

png pdf
Figure 7:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 7-a:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 7-b:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 7-c:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 7-d:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 7-e:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 7-f:
Distributions of ${m_{\text {T}}^{\text {tot}}}$ in the global (left) no b-tag and (right) b-tag categories in the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states, for the most signal sensitive categories. For the e$\mu$ final state, the Medium-${D_{\zeta}}$ category is displayed, for the e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$ final states the Tight-$ {m_{\mathrm {T}}}$ categories are shown. The black horizontal line in the upper panel of each subfigure indicates the change from logarithmic to linear scale on the vertical axis. The distributions are shown for all data-taking years combined.

png pdf
Figure 8:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 8-a:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 8-b:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 8-c:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 8-d:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 8-e:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 8-f:
Distributions of ${m_{\tau \tau}}$ in the most signal sensitive categories: the $100\leq {{p_{\mathrm {T}}} ^{{\tau \tau}}} < $ 200 GeV (left) and $ {{p_{\mathrm {T}}} ^{{\tau \tau}}} \geq $ 200 GeV (right) categories of the global no b -tag category used for the model-independent ${\phi}$ search for 60 $ \leq {m_{\phi}} < $ 250 GeV for the (upper row) e$\mu$, (middle row) e${\tau _\mathrm {h}}$ and $\mu {\tau _\mathrm {h}}$, and (lower row) ${{\tau _\mathrm {h}} {\tau _\mathrm {h}}}$ final states. The distributions are shown for all data-taking years combined.

png pdf
Figure 9:
Expected and observed 95% CL upper limits on the product of the cross sections and branching fraction into $\tau$ leptons for ${\mathrm{g} \mathrm{g} \phi}$ and ${\mathrm{b} \mathrm{b} \phi}$ production in a mass range of 60 $\leq {m_{\phi}} \leq $ 3500 GeV, in addition to ${\mathrm {h_{obs}}}$. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark green and bright yellow bands indicate the 68 and 95% central intervals for the expected exclusion limit. The black dots correspond to the observed limits.

png pdf
Figure 9-a:
Expected and observed 95% CL upper limits on the product of the cross sections and branching fraction into $\tau$ leptons for ${\mathrm{g} \mathrm{g} \phi}$ and ${\mathrm{b} \mathrm{b} \phi}$ production in a mass range of 60 $\leq {m_{\phi}} \leq $ 3500 GeV, in addition to ${\mathrm {h_{obs}}}$. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark green and bright yellow bands indicate the 68 and 95% central intervals for the expected exclusion limit. The black dots correspond to the observed limits.

png pdf
Figure 9-b:
Expected and observed 95% CL upper limits on the product of the cross sections and branching fraction into $\tau$ leptons for ${\mathrm{g} \mathrm{g} \phi}$ and ${\mathrm{b} \mathrm{b} \phi}$ production in a mass range of 60 $\leq {m_{\phi}} \leq $ 3500 GeV, in addition to ${\mathrm {h_{obs}}}$. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark green and bright yellow bands indicate the 68 and 95% central intervals for the expected exclusion limit. The black dots correspond to the observed limits.

png pdf
Figure 10:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-a:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-b:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-c:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-d:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-e:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-f:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-g:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-h:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 10-i:
Maximum likelihood estimates, and 68 and 95% confidence level contours obtained from scans of the signal likelihood for the model-independent ${\phi}$ search. The scans are shown for selected values of ${m_{\phi}}$ between 60 GeV and 3.5 TeV.

png pdf
Figure 11:
Expected and observed 95% CL upper limits on ${g_{\text {U}}}$ in the VLQ BM 1 (left) and 2 (right) scenarios, in a mass range of $1\leq {m_{\text {U}}} \leq 5 TeV $. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark and bright grey bands indicate the central 68 and 95% intervals of the expected exclusion limit. The observed excluded parameter space is indicated by the coloured blue area. For both scenarios, the 95% confidence interval for the preferred region from the global fit of the low-energy observables presented in Ref. [73] is also shown by the green shaded area.

png pdf
Figure 11-a:
Expected and observed 95% CL upper limits on ${g_{\text {U}}}$ in the VLQ BM 1 (left) and 2 (right) scenarios, in a mass range of $1\leq {m_{\text {U}}} \leq 5 TeV $. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark and bright grey bands indicate the central 68 and 95% intervals of the expected exclusion limit. The observed excluded parameter space is indicated by the coloured blue area. For both scenarios, the 95% confidence interval for the preferred region from the global fit of the low-energy observables presented in Ref. [73] is also shown by the green shaded area.

png pdf
Figure 11-b:
Expected and observed 95% CL upper limits on ${g_{\text {U}}}$ in the VLQ BM 1 (left) and 2 (right) scenarios, in a mass range of $1\leq {m_{\text {U}}} \leq 5 TeV $. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark and bright grey bands indicate the central 68 and 95% intervals of the expected exclusion limit. The observed excluded parameter space is indicated by the coloured blue area. For both scenarios, the 95% confidence interval for the preferred region from the global fit of the low-energy observables presented in Ref. [73] is also shown by the green shaded area.

png pdf
Figure 12:
Expected 95% CL exclusion contours in the MSSM (left) ${M_{\text {h}}^{125}}$ and (right) ${M_{\text {h},\rm EFT}^{125}}$ scenarios. The expected median in the absence of a signal is shown as a dashed black line. The dark and bright grey bands indicate the associated 68 and 95% intervals of the expected exclusion. The observed exclusion contour is indicated by the coloured blue area. For both scenarios, those parts of the parameter space where ${m_{\mathrm{h}}}$ deviates by more then ${\pm}$3 GeV from the mass of ${\mathrm {h_{obs}}}$ are indicated by a red hatched area. For the ${M_{\text {h},\rm EFT}^{125}}$ scenario, the dashed blue line indicates the $ {m_{{\mathrm {A}}}} =2 {m_{\mathrm{t}}} $ threshold whereby the $ {\mathrm {A}} \to {\mathrm{t} {}\mathrm{\bar{t}}} $ decay starts to influence the ${{\mathrm {A}} \to \tau \tau}$ branching fraction. The ${\mathrm{H} \to \tau \tau}$ branching fraction is influenced somewhat more gradually close to this threshold since the A and H are not completely degenerate in mass.

png pdf
Figure 12-a:
Expected 95% CL exclusion contours in the MSSM (left) ${M_{\text {h}}^{125}}$ and (right) ${M_{\text {h},\rm EFT}^{125}}$ scenarios. The expected median in the absence of a signal is shown as a dashed black line. The dark and bright grey bands indicate the associated 68 and 95% intervals of the expected exclusion. The observed exclusion contour is indicated by the coloured blue area. For both scenarios, those parts of the parameter space where ${m_{\mathrm{h}}}$ deviates by more then ${\pm}$3 GeV from the mass of ${\mathrm {h_{obs}}}$ are indicated by a red hatched area. For the ${M_{\text {h},\rm EFT}^{125}}$ scenario, the dashed blue line indicates the $ {m_{{\mathrm {A}}}} =2 {m_{\mathrm{t}}} $ threshold whereby the $ {\mathrm {A}} \to {\mathrm{t} {}\mathrm{\bar{t}}} $ decay starts to influence the ${{\mathrm {A}} \to \tau \tau}$ branching fraction. The ${\mathrm{H} \to \tau \tau}$ branching fraction is influenced somewhat more gradually close to this threshold since the A and H are not completely degenerate in mass.

png pdf
Figure 12-b:
Expected 95% CL exclusion contours in the MSSM (left) ${M_{\text {h}}^{125}}$ and (right) ${M_{\text {h},\rm EFT}^{125}}$ scenarios. The expected median in the absence of a signal is shown as a dashed black line. The dark and bright grey bands indicate the associated 68 and 95% intervals of the expected exclusion. The observed exclusion contour is indicated by the coloured blue area. For both scenarios, those parts of the parameter space where ${m_{\mathrm{h}}}$ deviates by more then ${\pm}$3 GeV from the mass of ${\mathrm {h_{obs}}}$ are indicated by a red hatched area. For the ${M_{\text {h},\rm EFT}^{125}}$ scenario, the dashed blue line indicates the $ {m_{{\mathrm {A}}}} =2 {m_{\mathrm{t}}} $ threshold whereby the $ {\mathrm {A}} \to {\mathrm{t} {}\mathrm{\bar{t}}} $ decay starts to influence the ${{\mathrm {A}} \to \tau \tau}$ branching fraction. The ${\mathrm{H} \to \tau \tau}$ branching fraction is influenced somewhat more gradually close to this threshold since the A and H are not completely degenerate in mass.
Tables

png pdf
Table 1:
Summary of the best fit values and uncertainties of ${{\beta _{\mathrm {L}}} ^{\mathrm {s}\tau}}$ in the two considered ${\mathrm {U}_1}$ benchmark scenarios from Ref. [73].

png pdf
Table 2:
Selected working points of ${D_{\mathrm{e}}}$, ${D_{\mu}}$, and ${D_{\text {jet}}}$ depending on the ${\tau \tau}$ final state.

png pdf
Table 3:
Offline selection requirements applied to the electron, muon, and ${\tau _\mathrm {h}}$ candidates used for the selection of the $\tau$ pair. First and second lepton refer to the label of the final state in the first column. For the ${p_{\mathrm {T}}}$ requirements, the values in parentheses correspond to events that have been recorded based on different trigger paths in the online selection, depending on the data-taking year. A detailed discussion is given in the text.

png pdf
Table 4:
Event categories and discriminants used for the extraction of the signals, for the searches described in this note. We note that ${m_{\phi}}$ refers to the hypothesized mass of the model-independent ${\phi}$ search, while ${m_{\tau \tau}}$ refers to the reconstructed mass of the ${\tau \tau}$ system before the decays of the $\tau$ leptons, and thus to an estimate of ${m_{\phi}}$ in data. The variable $y_{l}$ refers to the output functions of the NNs used for signal extraction. [101].

png pdf
Table 5:
Background processes contributing to the event selection, as given in Section 5. The symbol $\ell $ corresponds to an electron or muon. The second column refers to the experimental signature in the analysis, the last four columns indicate the estimation methods used to model each corresponding signature, as described in Sections 6.1-6.4.

png pdf
Table 6:
Summary of systematic uncertainties discussed in the text. The first column indicates the source of uncertainty; the second the processes that it applies to; the third the variation; and the last how it is correlated with other uncertainties. A checkmark is given also for partial correlations. More details are given in the text.

png pdf
Table 7:
Contribution of MSSM signals to the ${m_{\text {T}}^{\text {tot}}}$ and NN output function template distributions used for signal extraction for the interpretation of the data in MSSM benchmark scenarios.
Summary
Three searches have been presented for signatures of physics beyond the standard model (SM) in ${\tau\tau}$ final states in proton-proton collisions at the LHC. The searches use a sample of proton-proton collisions collected with the CMS detector at a centre-of-mass energy of $\sqrt{s} = $ 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The data have been analysed in three different interpretations: in the context of a model-independent search for a (pseudo-)scalar resonance $\phi$ in addition to the observed Higgs boson, 95% confidence level (CL) upper limits have been set on the product of the branching fraction for the decay into $\tau$ leptons and the cross section for the production via gluon fusion or in association with b quarks, spanning four orders of magnitude from $\mathcal{O}$(10 pb) (for a mass of ${m_{\phi}} =$ 60 GeV) to 0.3 fb (for ${m_{\phi}} =$ 3.5 TeV) each. In the context of a search for non-resonant $t$-channel exchange of a vector-like leptoquark ${\mathrm{U}_1} $, 95% CL upper limits have been set on the coupling ${g_{\text{U}}}$ ranging from 1 (for a mass of ${m_{\text{U}}} =$ 1 TeV) to 6 (for ${m_{\text{U}}} =$ 5 TeV), depending on the considered scenario. In the interpretation of benchmark scenarios of the minimal supersymmetric SM (MSSM), additional Higgs bosons with masses below 350 GeV are excluded at 95% CL. Depending on the scenario, the sensitivity of the data to exclude values of the MSSM parameter ${\tan\beta}$ reaches up to 1.8 TeV. The data reveal two excesses for $\phi$ production via gluon fusion with local $p$-values equivalent to approximately three s.d. at 100 GeV and 1.2 TeV, which are found to be consistent across ${\tau\tau}$ final states and data-taking years.
Additional Figures

png pdf
Additional Figure 1:
Local p-value and significance as a function of $m_{\phi}$ for gg$ \phi $ production. The bb$ \phi $ production rate is profiled.

png pdf
Additional Figure 2:
Local p-value and significance as a function of $m_{\phi}$ for bb$ \phi $ production. The gg$ \phi $ production rate is profiled.

png pdf
Additional Figure 3:
Local p-value and significance as a function of $m_{\phi}$ for gg$ \phi $ production. The bb$ \phi $ production rate is set to zero.

png pdf
Additional Figure 4:
Local p-value and significance as a function of $m_{\phi}$ for bb$ \phi $ production. The gg$ \phi $ production rate is set to zero.

png pdf
Additional Figure 5:
Expected and observed 95% CL upper limits on $g_{\text {U}}$ in the benchmark scenario ("VLQ BM 3'') in which $\beta _{\mathrm {L}}^{\mathrm{b} \tau}$ is taken to be 1 and the other couplings are set to zero, in a mass range of 1 $\leq m_{\text {U}}\leq $ 5 TeV. The expected median of the exclusion limit in the absence of signal is shown by the dashed line. The dark and bright grey bands indicate the central 68 and 95% intervals of the expected exclusion limit. The observed excluded parameter space is indicated by the coloured blue area.

png pdf
Additional Figure 6:
Distribution of $m_{\tau \tau}$ in the e$\mu $ channel after a background-only fit to the data. The 100 $\leq {p_{\mathrm {T}}} ^{\tau \tau} < $ 200 GeV category of the global no b-tag category used for the model-independent $\phi $ search are displayed. The cross section for the gg$ \phi $ process with $m_{\phi} = $ 100 GeV is scaled to 5.8 pb for illustrative purposes. The distributions are shown for all data-taking years combined.

png pdf
Additional Figure 7:
Distribution of $m_{\tau \tau}$ in the e$\mu $ channel after a background-only fit to the data. The $ {p_{\mathrm {T}}} ^{\tau \tau}\geq $ 200 GeV category of the global no b-tag category used for the model-independent $\phi $ search are displayed. The cross section for the gg$ \phi $ process with $m_{\phi} = $ 100 GeV is scaled to 5.8 pb for illustrative purposes. The distributions are shown for all data-taking years combined.

png pdf
Additional Figure 8:
Distribution of $m_{\tau \tau}$ in the e$\tau _{\mathrm {h}}$ and $\mu \tau _{\mathrm {h}}$ channels after a background-only fit to the data. The 100 $\leq {p_{\mathrm {T}}} ^{\tau \tau} < $ 200 GeV category of the global no b-tag category used for the model-independent $\phi $ search are displayed. The cross section for the gg$ \phi $ process with $m_{\phi} = $ 100 GeV is scaled to 5.8 pb for illustrative purposes. The distributions are shown for all data-taking years combined.

png pdf
Additional Figure 9:
Distribution of $m_{\tau \tau}$ in the e$\tau _{\mathrm {h}}$ and $\mu \tau _{\mathrm {h}}$ channel after a background-only fit to the data. The $ {p_{\mathrm {T}}} ^{\tau \tau}\geq $ 200 GeV category of the global no b-tag category used for the model-independent $\phi $ search are displayed. The cross section for the gg$ \phi $ process with $m_{\phi} = $ 100 GeV is scaled to 5.8 pb for illustrative purposes. The distributions are shown for all data-taking years combined.

png pdf
Additional Figure 10:
Distribution of $m_{\tau \tau}$ in the $\tau _{\mathrm {h}}\tau _{\mathrm {h}}$ channel after a background-only fit to the data. The 100 $\leq {p_{\mathrm {T}}} ^{\tau \tau} < $ 200 GeV category of the global no b-tag category used for the model-independent $\phi $ search are displayed. The cross section for the gg$ \phi $ process with $m_{\phi} = $ 100 GeV is scaled to 5.8 pb for illustrative purposes. The distributions are shown for all data-taking years combined.

png pdf
Additional Figure 11:
Distribution of $m_{\tau \tau}$ in the $\tau _{\mathrm {h}}\tau _{\mathrm {h}}$ channel after a background-only fit to the data. The $ {p_{\mathrm {T}}} ^{\tau \tau}\geq $ 200 GeV category of the global no b-tag category used for the model-independent $\phi $ search are displayed. The cross section for the gg$ \phi $ process with $m_{\phi} = $ 100 GeV is scaled to 5.8 pb for illustrative purposes. The distributions are shown for all data-taking years combined.

png pdf
Additional Figure 12:
Distribution of $m_{\text {T}}^{\text {tot}}$ in the e$\mu $ final state after a background-only fit to the data. The Medium-$D_{\zeta}$ category of the global no b-tag category used for the model-independent $\phi $ search is displayed. The cross section for the gg$ \phi $ and PbPb$\phi $ processes with $m_{\phi} = $ 1200 GeV are scaled to 3.1 fb and 1.0 fb for illustrative purposes. The distribution is shown for all data-taking years combined.

png pdf
Additional Figure 13:
Distribution of $m_{\text {T}}^{\text {tot}}$ in the e$\mu $ final state after a background-only fit to the data. The Medium-$D_{\zeta}$ category of the global b-tag category used for the model-independent $\phi $ search is displayed. The cross section for the gg$ \phi $ and PbPb$\phi $ processes with $m_{\phi} = $ 1200 GeV are scaled to 3.1 fb and 1.0 fb for illustrative purposes. The distribution is shown for all data-taking years combined.

png pdf
Additional Figure 14:
Distribution of $m_{\text {T}}^{\text {tot}}$ in the e$\tau _{\text {h}}$ and $\mu \tau _{\text {h}}$ final states after a background-only fit to the data. The Tight-$m_{\text {T}}$ category of the global no b-tag category used for the model-independent $\phi $ search is displayed. The cross section for the gg$ \phi $ and PbPb$\phi $ processes with $m_{\phi} = $ 1200 GeV are scaled to 3.1 fb and 1.0 fb for illustrative purposes. The distribution is shown for all data-taking years combined.

png pdf
Additional Figure 15:
Distribution of $m_{\text {T}}^{\text {tot}}$ in the e$\tau _{\text {h}}$ and $\mu \tau _{\text {h}}$ final states after a background-only fit to the data. The Medium-$m_{\text {T}}$ category of the global b-tag category used for the model-independent $\phi $ search is displayed. The cross section for the gg$ \phi $ and PbPb$\phi $ processes with $m_{\phi} = $ 1200 GeV are scaled to 3.1 fb and 1.0 fb for illustrative purposes. The distribution is shown for all data-taking years combined.

png pdf
Additional Figure 16:
Distribution of $m_{\text {T}}^{\text {tot}}$ in the $\tau _{\text {h}}\tau _{\text {h}}$ final states after a background-only fit to the data. The global no b-tag category used for the model-independent $\phi $ search is displayed. The cross section for the gg$ \phi $ and PbPb$\phi $ processes with $m_{\phi} = $ 1200 GeV are scaled to 3.1 fb and 1.0 fb for illustrative purposes. The distribution is shown for all data-taking years combined.

png pdf
Additional Figure 17:
Distribution of $m_{\text {T}}^{\text {tot}}$ in the $\tau _{\text {h}}\tau _{\text {h}}$ final states after a background-only fit to the data. The global b-tag category used for the model-independent $\phi $ search is displayed. The cross section for the gg$ \phi $ and PbPb$\phi $ processes with $m_{\phi} = $ 1200 GeV are scaled to 3.1 fb and 1.0 fb for illustrative purposes. The distribution is shown for all data-taking years combined.
References
1 ATLAS Collaboration Observation of a new particle in the search for the standard model Higgs boson with the ATLAS detector at the LHC PLB 716 (2012) 1 1207.7214
2 CMS Collaboration Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC PLB 716 (2012) 30 CMS-HIG-12-028
1207.7235
3 CMS Collaboration Observation of a new boson with mass near 125 GeV in pp collisions at $ \sqrt{s} = $ 7 and 8 TeV JHEP 06 (2013) 081 CMS-HIG-12-036
1303.4571
4 CMS Collaboration A measurement of the Higgs boson mass in the diphoton decay channel PLB 805 (2020) 135425 CMS-HIG-19-004
2002.06398
5 ATLAS and CMS Collaborations Measurements of the Higgs boson production and decay rates and constraints on its couplings from a combined ATLAS and CMS analysis of the LHC pp collision data at $ \sqrt{s}= $ 7 and 8 TeV JHEP 08 (2016) 045 1606.02266
6 CMS Collaboration Combined measurements of Higgs boson couplings in proton-proton collisions at $ \sqrt{s}=$ 13 TeV EPJC 79 (2019) 421 CMS-HIG-17-031
1809.10733
7 ATLAS Collaboration Combined measurements of Higgs boson production and decay using up to 80 fb$ ^{-1} $ of proton-proton collision data at $ \sqrt{s}= $ 13 TeV collected with the ATLAS experiment PRD 101 (2020) 012002 1909.02845
8 CMS Collaboration Measurements of the Higgs boson width and anomalous HVV couplings from on-shell and off-shell production in the four-lepton final state PRD 99 (2019) 112003 CMS-HIG-18-002
1901.00174
9 \relax Yu. A. Golfand and E. P. Likhtman Extension of the algebra of Poincaré group generators and violation of p invariance JEPTL 13 (1971)323
10 J. Wess and B. Zumino Supergauge transformations in four dimensions NPB 70 (1974) 39
11 P. Fayet Supergauge invariant extension of the Higgs mechanism and a model for the electron and its neutrino NPB 90 (1975) 104
12 P. Fayet Spontaneously broken supersymmetric theories of weak, electromagnetic and strong interactions PLB 69 (1977) 489
13 DELPHI, OPAL, ALEPH, L3 and LEP Working Group for Higgs Boson Searches Collaboration Search for neutral MSSM Higgs bosons at LEP EPJC 47 (2006) 547 hep-ex/0602042
14 CDF Collaboration Search for Higgs bosons predicted in Two-Higgs-Doublet models via decays to tau lepton pairs in 1.96 TeV {\mathrm{p}}\bar{{\mathrm{p}}} $ collisions PRL 103 (2009) 201801 0906.1014
15 D0 Collaboration Search for neutral Higgs bosons in the multi-b-jet topology in 5.2 fb$^{-1}$ of $ {\mathrm{p}}\bar{{\mathrm{p}}} $ collisions at $ \sqrt{s} = $ 1.96 TeV PLB 698 (2011) 97 1011.1931
16 D0 Collaboration Search for Higgs bosons decaying to $ \tau^{+}\tau^{-} $ pairs in $ {\mathrm{p}}\bar{{\mathrm{p}}} $ collisions at $ \sqrt{s} = $ 1.96 TeV PLB 707 (2012) 323 1106.4555
17 CDF Collaboration Search for Higgs bosons produced in association with b-quarks PRD 85 (2012) 032005 1106.4782
18 CMS Collaboration Search for a Higgs boson decaying into a b-quark pair and produced in association with b quarks in proton-proton collisions at 7 TeV PLB 722 (2013) 207 CMS-HIG-12-033
1302.2892
19 CMS Collaboration Search for neutral MSSM Higgs bosons decaying into a pair of bottom quarks JHEP 11 (2015) 071 CMS-HIG-14-017
1506.08329
20 CMS Collaboration Search for beyond the standard model Higgs bosons decaying into a $ \mathrm{b\bar{b}} $ pair in pp collisions at $ \sqrt{s} = $ 13 Tev JHEP 08 (2018) 113 CMS-HIG-16-018
1805.12191
21 ATLAS Collaboration Search for heavy neutral Higgs bosons produced in association with b-quarks and decaying into b-quarks at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD 102 (2020) 032004 1907.02749
22 ATLAS Collaboration Search for the neutral Higgs bosons of the minimal supersymmetric standard model in pp collisions at $ \sqrt{s} = $ 7 TeV with the ATLAS detector JHEP 02 (2013) 095 1211.6956
23 CMS Collaboration Search for neutral MSSM Higgs bosons decaying to $ \mu^{+} \mu^{-} $ in pp collisions at $ \sqrt{s} = $ 7 and 8 TeV PLB 752 (2016) 221 CMS-HIG-13-024
1508.01437
24 CMS Collaboration Search for MSSM higgs bosons decaying to $ \mu^{+}\mu^{-}$ in proton-proton collisions at $ \sqrt{s} = $ 13 TeV PLB 798 (2019) 134992 CMS-HIG-18-010
1907.03152
25 ATLAS Collaboration Search for scalar resonances decaying into $ \mu^{+}\mu^{-} $ in events with and without b-tagged jets produced in proton-proton collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector JHEP 07 (2019) 117 1901.08144
26 ATLAS Collaboration Search for neutral Higgs bosons of the minimal supersymmetric standard model in pp collisions at $ \sqrt{s} = $ 8 TeV with the ATLAS detector JHEP 11 (2014) 056 1409.6064
27 ATLAS Collaboration Search for minimal supersymmetric standard model Higgs bosons H/A and for a Z' boson in the $ \tau \tau $ final state produced in pp collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector EPJC 76 (2016) 585 1608.00890
28 ATLAS Collaboration Search for additional heavy neutral Higgs and gauge bosons in the ditau final state produced in $ 36 fb$^{-1} of pp collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector JHEP 01 (2018) 055 1709.07242
29 CMS Collaboration Search for neutral minimal supersymmetric standard model Higgs bosons decaying to tau pairs in pp collisions at $ \sqrt{s} = $ 7 TeV PRL 106 (2011) 231801 CMS-HIG-10-002
1104.1619
30 CMS Collaboration Search for neutral Higgs bosons decaying to tau pairs in pp collisions at $ \sqrt{s} = $ 7 TeV PLB 713 (2012) 68 CMS-HIG-11-029
1202.4083
31 CMS Collaboration Search for neutral MSSM Higgs bosons decaying to a pair of tau leptons in pp collisions JHEP 10 (2014) 160 CMS-HIG-13-021
1408.3316
32 CMS Collaboration Search for additional neutral MSSM higgs bosons in the $ \tau\tau $ final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 09 (2018) 007 CMS-HIG-17-020
1803.06553
33 ATLAS Collaboration Search for heavy Higgs bosons decaying into two tau leptons with the ATLAS detector using pp collisions at $ \sqrt{s}= $ 13 TeV PRL 125 (2020) 051801 2002.12223
34 J. Steggemann Extended scalar sectors Ann. Rev. Nucl. Part. Sci. 70 (2020) 197
35 B. Diaz, M. Schmaltz, and Y.-M. Zhong The leptoquark hunter's guide: Pair production JHEP 10 (2017) 097 1706.05033
36 M. Schmaltz and Y.-M. Zhong The leptoquark hunter's guide: Large coupling JHEP 01 (2019) 132 1810.10017
37 CMS Collaboration Search for third-generation scalar leptoquarks and heavy right-handed neutrinos in final states with two tau leptons and two jets in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 07 (2017) 121 CMS-EXO-16-023
1703.03995
38 CMS Collaboration Search for a singly produced third-generation scalar leptoquark decaying to a $ \tau $ lepton and a bottom quark in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 115 CMS-EXO-17-029
1806.03472
39 CMS Collaboration Search for heavy neutrinos and third-generation leptoquarks in hadronic states of two $ \tau $ leptons and two jets in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 03 (2019) 170 CMS-EXO-17-016
1811.00806
40 ATLAS Collaboration Searches for third-generation scalar leptoquarks in $ \sqrt{s} = $ 13 TeV pp collisions with the ATLAS detector JHEP 06 (2019) 144 1902.08103
41 ATLAS Collaboration Search for a scalar partner of the top quark in the all-hadronic $ \mathrm{t{\bar{t}}} $ plus missing transverse momentum final state at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 80 (2020) 737 2004.14060
42 CMS Collaboration Search for singly and pair-produced leptoquarks coupling to third-generation fermions in proton-proton collisions at $\sqrt{s} =$ 13 TeV PLB 819 (2021) 136446 CMS-EXO-19-015
2012.04178
43 ATLAS Collaboration Search for new phenomena in pp collisions in final states with tau leptons, b-jets, and missing transverse momentum with the ATLAS detector PRD 104 (2021) 112005 2108.07665
44 ATLAS Collaboration Search for new phenomena in final states with b-jets and missing transverse momentum in $ \sqrt{s}= $ 13 TeV pp collisions with the ATLAS detector JHEP 05 (2021) 093 2101.12527
45 ATLAS Collaboration Search for pair production of third-generation scalar leptoquarks decaying into a top quark and a $ \tau $-lepton in pp collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector JHEP 06 (2021) 179 2101.11582
46 J. C. Pati and A. Salam Unified lepton-hadron symmetry and a gauge theory of the basic interactions PRD 8 (1973) 1240
47 J. C. Pati and A. Salam Lepton number as the fourth color PRD 10 (1974) 275
48 H. Georgi and S. L. Glashow Unity of all elementary particle forces PRL 32 (1974) 438
49 H. Fritzsch and P. Minkowski Unified interactions of leptons and hadrons Ann. Phys. 93 (1975) 193
50 S. Dimopoulos and L. Susskind Mass without scalars NPB 155 (1979) 237
51 S. Dimopoulos Technicolored signatures NPB 168 (1980) 69
52 E. Farhi and L. Susskind Technicolor PR 74 (1981) 277
53 K. D. Lane and M. V. Ramana Walking technicolor signatures at hadron colliders PRD 44 (1991) 2678
54 B. Schrempp and F. Schrempp Light leptoquarks PLB 153 (1985) 101
55 B. Gripaios Composite leptoquarks at the LHC JHEP 02 (2010) 045 0910.1789
56 G. R. Farrar and P. Fayet Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry PLB 76 (1978) 575
57 P. Ramond Dual theory for free fermions PRD 3 (1971) 2415
58 A. Neveu and J. H. Schwarz Factorizable dual model of pions NPB 31 (1971) 86
59 D. V. Volkov and V. P. Akulov Possible universal neutrino interaction JEPTL 16 (1972)438
60 J. Wess and B. Zumino A Lagrangian model invariant under supergauge transformations PLB 49 (1974) 52
61 H. P. Nilles Supersymmetry, supergravity and particle physics PR 110 (1984) 1
62 R. Barbier et al. R-parity violating supersymmetry PR 420 (2005) 1 hep-ph/0406039
63 M. Tanaka and R. Watanabe New physics in the weak interaction of $ \bar B\to D^{(*)}\tau\bar\nu $ PRD 87 (2013) 034028 1212.1878
64 R. Barbieri, G. Isidori, A. Pattori, and F. Senia Anomalies in b-decays and U(2) flavour symmetry EPJC 76 (2016) 67 1512.01560
65 D. A. Faroughy, A. Greljo, and J. F. Kamenik Confronting lepton flavor universality violation in B decays with high-$ {p_{\mathrm{T}}} $ tau lepton searches at LHC PLB 764 (2017) 126 1609.07138
66 M. Bordone, C. Cornella, J. Fuentes-Martin, and G. Isidori A three-site gauge model for flavor hierarchies and flavor anomalies PLB 779 (2018) 317 1712.01368
67 L. Di Luzio, A. Greljo, and M. Nardecchia Gauge leptoquark as the origin of B-physics anomalies PRD 96 (2017) 115011 1708.08450
68 A. Greljo and B. A. Stefanek Third family quark-lepton unification at the TeV scale PLB 782 (2018) 131 1802.04274
69 A. Angelescu et al. Single leptoquark solutions to the B-physics anomalies PRD 104 (2021) 055017 2103.12504
70 C. Cornella et al. Reading the footprints of the B-meson flavor anomalies JHEP 08 (2021) 050 2103.16558
71 LHCb Collaboration Test of lepton universality in beauty-quark decays Submitted to Nat. Phys 2103.11769
72 BaBar Collaboration Evidence for an excess of $ \bar{B} \to D^{(*)} \tau^-\bar{\nu}_\tau $ decays PRL 109 (2012) 101802 1205.5442
73 BaBar Collaboration Measurement of an excess of $ \bar{B} \to D^{(*)}\tau^- \bar{\nu}_\tau $ decays and implications for charged higgs bosons PRD 88 (2013) 072012 1303.0571
74 Belle Collaboration Measurement of the branching ratio of $ \bar{B} \to D^{(\ast)} \tau^- \bar{\nu}_\tau $ relative to $ \bar{B} \to D^{(\ast)} \ell^- \bar{\nu}_\ell $ decays with hadronic tagging at Belle PRD 92 (2015) 072014 1507.03233
75 LHCb Collaboration Measurement of the ratio of branching fractions $ \mathcal{B}(\bar{B}^0 \to D^{*+}\tau^{-}\bar{\nu}_{\tau})/\mathcal{B}(\bar{B}^0 \to D^{*+}\mu^{-}\bar{\nu}_{\mu}) $ PRL 115 (2015) 111803 1506.08614
76 Belle Collaboration Measurement of the $ \tau $ lepton polarization and $ R(D^*) $ in the decay $ \bar{B} \to D^* \tau^- \bar{\nu}_\tau $ PRL 118 (2017) 211801 1612.00529
77 LHCb Collaboration Test of lepton flavor universality by the measurement of the $ B^0 \to D^{*-} \tau^+ \nu_{\tau} $ branching fraction using three-prong $ \tau $ decays PRD 97 (2018) 072013 1711.02505
78 LHCb Collaboration Measurement of the ratio of the $ B^0 \to D^{*-} \tau^+ \nu_{\tau} $ and $ B^0 \to D^{*-} \mu^+ \nu_{\mu} $ branching fractions using three-prong $ \tau $-lepton decays PRL 120 (2018) 171802 1708.08856
79 T. D. Lee A theory of spontaneous T violation PRD (1973) 1226
80 G. C. Branco et al. Theory and phenomenology of two-Higgs-doublet models PR (2012) 1 1106.0034
81 E. Bagnaschi et al. MSSM Higgs boson searches at the LHC: Benchmark scenarios for Run 2 and beyond EPJC 79 (2019) 617 1808.07542
82 H. Bahl et al. HL-LHC and ILC sensitivities in the hunt for heavy Higgs bosons EPJC 80 (2020) 916 2005.14536
83 H. Bahl, S. Liebler, and T. Stefaniak MSSM Higgs benchmark scenarios for Run 2 and beyond: the low $ \tan\beta $ region EPJC 79 (2019) 279 1901.05933
84 E. A. Bagnaschi et al. Benchmark scenarios for MSSM Higgs boson searches at the LHC technical report, CERN, Geneva
85 LHC Higgs Cross Section Working Group -- MSSM subgroup LHCHWG MSSM ROOT files 2022
86 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
87 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
88 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
89 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
90 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
91 CMS Collaboration Performance of electron reconstruction and selection with the CMS detector in proton-proton collisions at $ \sqrt{s} = $ 8 TeV JINST 10 (2015) P06005 CMS-EGM-13-001
1502.02701
92 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC Submitted to JINST CMS-EGM-17-001
2012.06888
93 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
94 CMS Collaboration Identification of heavy-flavour jets with the CMS detector in pp collisions at 13 TeV JINST 13 (2018) P05011 CMS-BTV-16-002
1712.07158
95 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020) P12012 2008.10519
96 CMS Collaboration Performance of the DeepJet b tagging algorithm using 41.9/fb of data from proton-proton collisions at 13 TeV with Phase 1 CMS detector CDS
97 CMS Collaboration Performance of reconstruction and identification of $ \tau $ leptons decaying to hadrons and $ \nu_\tau $ in pp collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P10005 CMS-TAU-16-003
1809.02816
98 CMS Collaboration Identification of hadronic tau lepton decays using a deep neural network Submitted to JINST CMS-TAU-20-001
2201.08458
99 CMS Collaboration Performance of missing transverse momentum reconstruction in proton-proton collisions at $ \sqrt{s} = $ 13 TeV using the CMS detector JINST 14 (2019) P07004 CMS-JME-17-001
1903.06078
100 D. Bertolini, P. Harris, M. Low, and N. Tran Pileup per particle identification JHEP 10 (2014) 059 1407.6013
101 CDF Collaboration Search for neutral Higgs bosons of the minimal supersymmetric standard model decaying to $ \tau $ pairs in $ {\mathrm{p}}\bar{{\mathrm{p}}} $ collisions at $ \sqrt{s} = $ 1.96 TeV PRL 96 (2006) 011802 hep-ex/0508051
102 L. Bianchini, J. Conway, E. K. Friis, and C. Veelken Reconstruction of the Higgs mass in $ \mathrm{H}\to\tau\tau $ events by dynamical likelihood techniques J. Phys. Conf. Ser. 513 (2014) 022035
103 B. K. Bullock, K. Hagiwara, and A. D. Martin Tau polarization and its correlations as a probe of new physics NPB 395 (1993) 499
104 Particle Data Group Collaboration Review of particle physics PTEP 2020 (2020) 083C01
105 CMS Collaboration Measurement of Higgs boson production and decay to the $ \tau\tau $ final state CMS-PAS-HIG-18-032 CMS-PAS-HIG-18-032
106 CMS Collaboration An embedding technique to determine $ \tau\tau $ backgrounds in proton-proton collision data JINST 14 (2019) P06032 CMS-TAU-18-001
1903.01216
107 CMS Collaboration Measurement of the $ \mathrm{Z}/\gamma^{*} \to \tau\tau $ cross section in pp collisions at $ \sqrt{s} = $ 13 TeV and validation of $ \tau $ lepton analysis techniques EPJC 78 (2018) 708 CMS-HIG-15-007
1801.03535
108 J. Alwall et al. MadGraph 5: Going beyond JHEP 06 (2011) 128 1106.0522
109 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
110 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
111 J. Alwall et al. Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions EPJC 53 (2008) 473 0706.2569
112 S. Alioli, S.-O. Moch, and P. Uwer Hadronic top-quark pair-production with one jet and parton showering JHEP 01 (2012) 137 1110.5251
113 R. Frederix, E. Re, and P. Torrielli Single-top $ t $-channel hadroproduction in the four-flavour scheme with POWHEG and aMC@NLO JHEP 09 (2012) 130 1207.5391
114 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
115 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
116 S. Alioli, P. Nason, C. Oleari, and E. Re A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX JHEP 06 (2010) 043 1002.2581
117 T. Je\vzo and P. Nason On the treatment of resonances in next-to-leading order calculations matched to a parton shower JHEP 12 (2015) 065 1509.09071
118 E. Re Single-top Wt-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
119 K. Melnikov and F. Petriello Electroweak gauge boson production at hadron colliders through $ \mathcal{O}(\alpha_\text{s}^{2}) $ PRD 74 (2006) 114017 hep-ph/0609070
120 M. Czakon and A. Mitov Top++: A program for the calculation of the top-pair cross-section at hadron colliders CPC 185 (2014) 2930 1112.5675
121 N. Kidonakis Top quark production in Proceedings, Helmholtz International Summer School on Physics of Heavy Quarks and Hadrons (HQ 2013): JINR, Dubna, Russia, July 15-28, 2013, p. 139 2014 1311.0283
122 J. M. Campbell, R. K. Ellis, and C. Williams Vector boson pair production at the LHC JHEP 07 (2011) 018 1105.0020
123 T. Gehrmann et al. $ W^+W^- $ production at hadron colliders in next to next to leading order QCD PRL 113 (2014) 212001 1408.5243
124 E. Bagnaschi, G. Degrassi, P. Slavich, and A. Vicini Higgs production via gluon fusion in the POWHEG approach in the SM and in the MSSM JHEP 02 (2012) 088 1111.2854
125 P. Nason and C. Oleari NLO Higgs boson production via vector-boson fusion matched with shower in POWHEG JHEP 02 (2010) 037 0911.5299
126 G. Luisoni, P. Nason, C. Oleari, and F. Tramontano $ HW^{\pm} $/HZ + 0 and 1 jet at NLO with the POWHEG BOX interfaced to GoSam and their merging within MiNLO JHEP 10 (2013) 083 1306.2542
127 F. Granata, J. M. Lindert, C. Oleari, and S. Pozzorini NLO QCD+EW predictions for HV and HV+jet production including parton-shower effects JHEP 09 (2017) 012 1706.03522
128 K. Hamilton, P. Nason, E. Re, and G. Zanderighi NNLOPS simulation of Higgs boson production JHEP 10 (2013) 222 1309.0017
129 K. Hamilton, P. Nason, and G. Zanderighi Finite quark-mass effects in the NNLOPS POWHEG+MiNLO Higgs generator JHEP 05 (2015) 140 1501.04637
130 R. V. Harlander, H. Mantler, and M. Wiesemann Transverse momentum resummation for Higgs production via gluon fusion in the MSSM JHEP 11 (2014) 116 1409.0531
131 E. Bagnaschi et al. Resummation ambiguities in the Higgs transverse-momentum spectrum in the standard model and beyond JHEP 01 (2016) 090 1510.08850
132 E. Bagnaschi and A. Vicini The Higgs transverse momentum distribution in gluon fusion as a multiscale problem JHEP 01 (2016) 056 1505.00735
133 B. Jager, L. Reina, and D. Wackeroth Higgs boson production in association with b jets in the POWHEG BOX PRD 93 (2016) 014030 1509.05843
134 M. J. Baker, J. Fuentes-Mart\'\in, G. Isidori, and M. Konig High- $ {p_{\mathrm{T}}} $ signatures in vector-leptoquark models EPJC 79 (2019) 334 1901.10480
135 J. Butterworth et al. PDF4LHC recommendations for LHC Run II JPG 43 (2016) 023001 1510.03865
136 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
137 NNPDF Collaboration Parton distributions for the LHC Run II JHEP 04 (2015) 040 1410.8849
138 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
139 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
140 T. Sjostrand et al. An introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
141 S. Agostinelli et al. GEANT4--a simulation toolkit NIMA 506 (2003) 250
142 R. J. Barlow and C. Beeston Fitting using finite Monte Carlo samples CPC 77 (1993) 219
143 J. S. Conway Incorporating nuisance parameters in likelihoods for multisource spectra in PHYSTAT 2011, p. 115 2011 1103.0354
144 CMS Collaboration Measurements of inclusive W and Z cross sections in pp collisions at $ \sqrt{s} = $ 7 TeV JHEP 01 (2011) 080 CMS-EWK-10-002
1012.2466
145 CMS Collaboration Measurement of the differential cross section for top quark pair production in pp collisions at $ \sqrt{s} = $ 8 TeV EPJC 75 (2015) 542 CMS-TOP-12-028
1505.04480
146 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
147 CMS Collaboration CMS luminosity measurement for the 2017 data taking period at $ \sqrt{s} = $ 13 TeV CMS-PAS-LUM-17-004 CMS-PAS-LUM-17-004
148 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} = $ 13 TeV CMS-PAS-LUM-18-002 CMS-PAS-LUM-18-002
149 LHC Higgs Cross Section Working Group Handbook of LHC Higgs cross sections: 4. deciphering the nature of the Higgs sector technical report, CERN 1610.07922
150 S. Alioli, P. Nason, C. Oleari, and E. Re NLO Higgs boson production via gluon fusion matched with shower in POWHEG JHEP 04 (2009) 002 0812.0578
151 H. B. Hartanto, B. Jager, L. Reina, and D. Wackeroth Higgs boson production in association with top quarks in the POWHEG BOX PRD 91 (2015) 094003 1501.04498
152 ATLAS and CMS Collaborations Procedure for the LHC Higgs boson search combination in summer 2011 ATL-PHYS-PUB 2011-11, CMS NOTE 2011/005,
153 CMS Collaboration Combined results of searches for the standard model Higgs boson in pp collisions at $ \sqrt{s} = $ 7 TeV PLB 710 (2012) 26 CMS-HIG-11-032
1202.1488
154 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
155 T. Junk Confidence level computation for combining searches with small statistics NIMA 434 (1999) 435 hep-ex/9902006
156 A. L. Read Presentation of search results: The CL$ _{\text{s}} $ technique JPG 28 (2002) 2693
157 CMS Collaboration Search for a standard model-like Higgs boson in the mass range between 70 and 110 GeV in the diphoton final state in proton-proton collisions at $ \sqrt{s}= $ 8 and 13 TeV PLB (2019) 320 CMS-HIG-17-013
1811.08459
158 P. Slavich et al. Higgs-mass predictions in the MSSM and beyond EPJC 81 (2021) 450 2012.15629
159 A. L. Read Linear interpolation of histograms NIMA 425 (1999) 357
160 S. Heinemeyer, W. Hollik, and G. Weiglein FeynHiggs: A program for the calculation of the masses of the neutral CP-even Higgs bosons in the MSSM CPC 124 (2000) 76 hep-ph/9812320
161 S. Heinemeyer, W. Hollik, and G. Weiglein The masses of the neutral CP-even Higgs bosons in the MSSM: Accurate analysis at the two-loop level EPJC 9 (1999) 343 hep-ph/9812472
162 G. Degrassi et al. Towards high-precision predictions for the MSSM Higgs sector EPJC 28 (2003) 133 hep-ph/0212020
163 M. Frank et al. The Higgs boson masses and mixings of the complex MSSM in the Feynman-diagrammatic approach JHEP 02 (2007) 047 hep-ph/0611326
164 T. Hahn et al. High-precision predictions for the light CP-even Higgs boson mass of the minimal supersymmetric standard model PRL 112 (2014) 141801 1312.4937
165 H. Bahl and W. Hollik Precise prediction for the light MSSM Higgs boson mass combining effective field theory and fixed-order calculations EPJC 76 (2016), no. 9, 499 1608.01880
166 H. Bahl, S. Heinemeyer, W. Hollik, and G. Weiglein Reconciling EFT and hybrid calculations of the light MSSM Higgs-boson mass EPJC 78 (2018), no. 1, 57 1706.00346
167 H. Bahl et al. Precision calculations in the MSSM Higgs-boson sector with FeynHiggs 2.14 CPC 249 (2020) 107099 1811.09073
168 A. Djouadi, J. Kalinowski, and M. Spira HDECAY: A program for Higgs boson decays in the standard model and its supersymmetric extension CPC 108 (1998) 56 hep-ph/9704448
169 A. Djouadi, J. Kalinowski, M. Muehlleitner, and M. Spira HDECAY: Twenty$ _{++} $ years after CPC 238 (2019) 214 1801.09506
170 LHC Higgs Cross Section Working Group Handbook of LHC Higgs cross sections: 3. Higgs properties technical report, CERN 1307.1347
171 A. Denner et al. Standard model Higgs-boson branching ratios with uncertainties EPJC 71 (2011) 1753 1107.5909
172 R. V. Harlander, S. Liebler, and H. Mantler SusHi: A program for the calculation of Higgs production in gluon fusion and bottom-quark annihilation in the standard model and the MSSM CPC 184 (2013) 1605 1212.3249
173 R. V. Harlander, S. Liebler, and H. Mantler SusHi Bento: Beyond NNLO and the heavy-top limit CPC 212 (2017) 239 1605.03190
174 M. Spira, A. Djouadi, D. Graudenz, and P. M. Zerwas Higgs boson production at the LHC NPB 453 (1995) 17 hep-ph/9504378
175 R. Harlander and P. Kant Higgs production and decay: analytic results at next-to-leading order QCD JHEP 12 (2005) 015 hep-ph/0509189
176 R. V. Harlander and W. B. Kilgore Next-to-next-to-leading order Higgs production at hadron colliders PRL 88 (2002) 201801 hep-ph/0201206
177 C. Anastasiou and K. Melnikov Higgs boson production at hadron colliders in NNLO QCD NPB 646 (2002) 220 hep-ph/0207004
178 V. Ravindran, J. Smith, and W. L. van Neerven NNLO corrections to the total cross-section for Higgs boson production in hadron-hadron collisions NPB 665 (2003) 325 hep-ph/0302135
179 R. V. Harlander and W. B. Kilgore Production of a pseudo-scalar Higgs boson at hadron colliders at next-to-next-to leading order JHEP 10 (2002) 017 hep-ph/0208096
180 C. Anastasiou and K. Melnikov Pseudoscalar Higgs boson production at hadron colliders in next-to-next-to-leading order QCD PRD 67 (2003) 037501 hep-ph/0208115
181 C. Anastasiou et al. Higgs boson gluon-fusion production beyond threshold in N$ ^{3}LO $ QCD JHEP 03 (2015) 091 1411.3584
182 C. Anastasiou et al. Soft expansion of double-real-virtual corrections to Higgs production at N$ ^{3} $LO JHEP 08 (2015) 051 1505.04110
183 C. Anastasiou et al. High precision determination of the gluon fusion Higgs boson cross-section at the LHC JHEP 05 (2016) 058 1602.00695
184 U. Aglietti, R. Bonciani, G. Degrassi, and A. Vicini Two-loop light fermion contribution to Higgs production and decays PLB 595 (2004) 432 hep-ph/0404071
185 R. Bonciani, G. Degrassi, and A. Vicini On the generalized harmonic polylogarithms of one complex variable CPC 182 (2011) 1253 1007.1891
186 G. Degrassi and P. Slavich NLO QCD bottom corrections to Higgs boson production in the MSSM JHEP 11 (2010) 044 1007.3465
187 G. Degrassi, S. Di Vita, and P. Slavich NLO QCD corrections to pseudoscalar Higgs production in the MSSM JHEP 08 (2011) 128 1107.0914
188 G. Degrassi, S. Di Vita, and P. Slavich On the NLO QCD corrections to the production of the heaviest neutral Higgs scalar in the MSSM EPJC 72 (2012) 2032 1204.1016
189 L. Hofer, U. Nierste, and D. Scherer Resummation of tan-beta-enhanced supersymmetric loop corrections beyond the decoupling limit JHEP 10 (2009) 081 0907.5408
190 M. Bonvini, A. S. Papanastasiou, and F. J. Tackmann Resummation and matching of b-quark mass effects in bbH production JHEP 11 (2015) 196 1508.03288
191 M. Bonvini, A. S. Papanastasiou, and F. J. Tackmann Matched predictions for the bbH cross section at the 13 TeV LHC JHEP 10 (2016) 053 1605.01733
192 S. Forte, D. Napoletano, and M. Ubiali Higgs production in bottom-quark fusion in a matched scheme PLB 751 (2015) 331 1508.01529
193 S. Forte, D. Napoletano, and M. Ubiali Higgs production in bottom-quark fusion: matching beyond leading order PLB 763 (2016) 190 1607.00389
Compact Muon Solenoid
LHC, CERN