Exclusion of heavy, broad resonances from precise measurements of WZ and VH final states at the LHC

TL;DR

This paper proposes a new search strategy for heavy vector resonances in specific final states at the LHC, using recent Higgs data to extend limits on resonance widths and explore discovery potential across various decay channels.

Contribution

It introduces a novel search method for heavy vector resonances involving multiple final states and uses a general vector-triplet model to evaluate discovery prospects and set new limits.

Findings

Excesses are expected in all final states, with the 1-lepton channel being most sensitive.

Current limits on resonance widths are extended up to 9% of the resonance mass.

The proposed approach enhances the discovery potential for heavy vector resonances at the LHC.

Abstract

A novel search for heavy vector resonances in the $H\rightarrow b\bar{b}$ and $Z\rightarrow b\bar{b}$ final states in association with a leptonically decaying $V$ ($Z$ or $W$) and $W$-only respectively, is proposed. It is argued that excesses with respect to the Standard Model prediction should be observed in all final states (0, 1 or 2 leptons), with the 1-lepton final state being the strongest. Since the relative strengths of these excesses depend on branching ratios and efficiencies, this is a clear signal for the presence of heavy resonances or their low mass tails. A general vector-triplet model is used to explore the discovery potential as a function of the resonance mass and width. Recent Higgs to $b\bar{b}$ observation data reported by the experiments ATLAS and CMS are used to test the model. Current limits are extended to resonance widths over mass as large as 9%.

Figures

Tables

Table 1. Table 1: Signal yield and resonant SM background yields for a 1.5 TeV resonance, g V = 5 subscript 𝑔 𝑉 5 g_{V}=5 and an integrated luminosity of 100 fb -1 , after baseline selection, and after an additional p ⟂ > 200 subscript 𝑝 perpendicular-to 200 p_{\perp}>200  GeV cut (last column). The resonance width over its mass is Γ / m ≃ similar-to-or-equals Γ 𝑚 absent \Gamma/m\simeq  10%.

Process	$\sigma \times BR$	A$\times ϵ$	Yield	Yield
$qq\to V^{\prime }$	[fb]	[%]	[$100$fb-1]	$p_{⟂}>$
				$200$ GeV
$Z^{\prime }\to ZH\to \nu \nu b\overline{b}$	1.58	2.45	3.78	3.51
$Z\to ZH\to \nu \nu b\overline{b}$	97.2	5.01	486	224
$ZZ\to \nu \nu jj$	2580	0.27	697	248
$W^{\prime }\to WH\to \ell \nu b\overline{b}$	3.87	5.5	21.3	19.4
$W^{\prime }\to WZ\to \ell \nu jj$	3.79	0.23	0.81	0.59
$W\to WH\to \ell \nu b\overline{b}$	225	1.44	324	148
$WZ\to \ell \nu jj$	4148	0.13	529	173
$Z^{\prime }\to ZH\to \ell \ell b\overline{b}$	0.553	4.76	2.6	2.2
$Z\to ZH\to \ell \ell b\overline{b}$	34.2	13.7	467	68.6
$ZZ\to \ell \ell jj$	910	0.96	875	72.6

Equations

$${\mathcal{L}}_V^{int}=$$

Full text

††journal: Eur. Phys. J. C

1 Introduction

Heavy vector resonances naturally appear in several extensions of the Standard Model (SM), such as GUT theories Pati ; Georgi ; Fritzcsh , composite Higgs Eichten ; Contino , little Higgs lHiggs1 ; lHiggs2 ; lHiggs3 , and models with vector $Z^{\prime}$ Zprime1 ; Zprime2 , and $W^{\prime}$ models Wprime . The LHC experiments, in most cases, have performed direct searches for heavy narrow resonances decaying to dibosons and $VH$ final states and have put limits to masses up to about 5.5 TeV ATLASVV_jjjj_19Dec17 ; ATLASVV_02lep_19Aug17 ; ATLASVV_1lep_19Oct17 ; ATLASVHbb17Dec17 ; ATLASVHbb_jj_21Jul17 ; CMSVHbb_leptons_8Jul18 ; CMSVV_llqq_6Dec17 ; CMSVV_lvqq_12Dec17 ; CMSVV_vvqq_11Jul17 ; CMS_VH_17 ; CMSVVandVH_17 ; ATLASHeavyCombination .

Broader resonances ($\Gamma/m>4\%$) appear for larger values of the resonance coupling strength to weak bosons and the Higgs, $g_{V}$, and depending on their mass, width, and production cross section, the LHC experiments can be sensitive to a large part of the resonance distribution. It is interesting to study to what extent the experiments can be sensitive to the full distribution or tails of such resonances and what the best observables are. In this work, we start addressing this question with a specific final state: a pair of $b$ quarks produced in association with a weak gauge boson ($Z$ or $W$) that subsequently decays leptonically (where leptons $\ell$ are electrons and muons). The interest in this final state is because in the large coupling regime, $g_{V}>3$, the heavy boson branching fraction is dominated by decays to $WZ$ and $VH$, while the decays to fermions are suppressed ourCHpaper1 . At the same time, both the Drell-Yan $q\bar{q}$ and the VBF production modes contribute, providing non-negligible cross sections. Here, we propose as an early indirect signal of the presence of vector resonances, the simultaneous excess in the 0, 1, and 2-lepton final states, for both $WZ\rightarrow\ell\nu b\bar{b}$ and $VH\rightarrow\ell\ell(\nu)b\bar{b}$, in a measurement of the $b\bar{b}$ quark invariant mass $M_{bb}$. The strength of these excesses in the different final states can be predicted by branching ratios and experimental efficiencies. We argue that the 0, 1, and 2-lepton final states should all be sensitive, albeit with very different sensitivity, to exotic $V^{\prime}\rightarrow VH$ decays. In particular, the 1-lepton state will be the most sensitive in Higgs excesses and in addition, it should also show a smaller level of excess in the $Z$-boson $b\bar{b}$ peak. The pattern of these excesses should be a clear sign of a heavy vector triplet (HVT), even before its observation as a (broad) mass resonance at the TeV scale. Using a general HVT framework that is also used by the experiments to model such heavy vector triplets hvt1 , we derive the sensitivity of a single LHC experiment as a function of the resonance width.

In this paper, we first introduce the HVT framework and summarize the recent LHC results. In Section 3 we present an analysis used to search for heavy vector triplets at the LHC, and we report on the search potential for the $b\bar{b}+$0, 1, and 2-lepton final states, for a range of masses and widths. Finally, the compatibility of our predictions with the recently published 13 TeV $H\rightarrow b\bar{b}$ data from ATLAS and CMS in terms of expected and observed limits, is discussed.

2 BSM Heavy Vector Resonances

A general HVT phenomenological Lagrangian can be used for the modelling of resonances predicted by a wide range of Beyond the Standard Model (BSM) scenarios hvt1 . The Lagrangian describing the interactions of these resonances $V^{a\prime}$, $a=1,2,3$ with quarks, leptons, vector bosons and the Higgs boson is shown below:

[TABLE]

where $q_{k}$ and $\ell_{k}$ are the quark and lepton weak doublets, $H$ is the Higgs doublet and $\sigma^{a}$ the three Pauli matrices. In this Lagrangian, the HVT triplet $V^{a\prime}=(W^{+\prime},W^{-\prime},Z^{\prime})$ interacts with the Higgs doublet, i.e. the longitudinal degrees of freedom of the SM $W$ and $Z$ bosons and the SM Higgs, with a coupling strength $g_{V}$. In order to allow for a broader class of models, this coupling strength can be varied by the parameter $c_{H}$, so in the Lagrangian the full coupling to the SM weak and Higgs bosons is $g_{V}c_{H}$. The HVT resonances also couple to the SM fermions, again through their coupling to the SM weak and Higgs bosons, $g^{2}/g_{V}$, where $g$ is the SM $SU(2)_{L}$ weak gauge coupling. This coupling between HVT resonances and fermions is also controlled by an additional parameter $c_{F}$ to allow for a broader range of models to be included, as follows: $g^{2}c_{F}/g_{V}$. As discussed in ATLASHeavyCombination , the experiments consider two Drell Yan (DY) production scenarios: model A is a scenario that reproduces the phenomenology of weakly coupled models based on an extended gauge symmetry Zprime1 . In this case, the couplings are $\frac{g^{2}c_{F}}{g_{V}}=-0.55$ and $g_{V}c_{H}=-0.56$, with the fermion coupling being universal. The second DY scenario, referred to as model B, implements a strongly coupled scenario as in composite Higgs models with $\frac{g^{2}c_{F}}{g_{V}}=0.14$ and $g_{V}c_{H}=-2.9$. In model B, the $V^{\prime}$ resonances are broader than in the weakly coupled scenario, model A, but for $|g_{V}c_{H}|\leq 3$ they remain narrow relative to the experimental resolution. For $|g_{V}c_{H}|>3$, the resonance intrinsic width becomes significant and cannot be neglected. In summary, the ATLAS and CMS benchmarks correspond for Model A to $c_{H}=-\frac{g^{2}}{g_{V}^{2}}$, $c_{F}=-\frac{1}{3}$ and $g_{V}=1$, and for model B to $c_{H}=-1$, $c_{F}=1$ and $g_{V}=3$. ATLAS and CMS experimental data from direct $m_{VH}$ and $m_{VV}$ searches, exclude part of the parameter phase space $(g_{V},c_{H},c_{F})$, for which the intrinsic width $\Gamma$ of the new bosons is dominated by the experimental resolution ($\sim 4\%$ of the mass). This is also the case, for example, for $|g_{V}c_{H}|\leq 3$, where negligible natural width is assumed (narrow width assumption). The goal of this work is to explore the part of the HVT parameter space with $|g_{V}c_{H}|>3$, where the heavy resonances have a significant natural width. Following the model B benchmark, we fix the two constant factors to $c_{H}=-1$, $c_{F}=1$ and allow $g_{V}$ to vary.

Results from recent experimental analyses that exclude part of the HVT phase space are shown in Fig. 1. The validity of these exclusions is actually limited due to the narrow width approximation. For larger, $|g_{V}c_{H}|>3$, couplings, the width grows and becomes comparable to or larger than the experimental resolution, leaving a large part of the HVT phase space unexcluded: $3<|g_{V}c_{H}|<4\pi$, where $g_{V}\simeq 4\pi$ is the perturbative limit. This is true even for resonance masses of order TeV.

In this work we argue that the presence of non-zero width resonances, could be observed as an excess of events in the $VH\rightarrow\ell\ell(\nu)b\bar{b}$ and $WZ\rightarrow\ell\nu b\bar{b}$ decays with 0, 1, or 2 leptons in the final state. Independent of the fact that due to the experimental resolution the close-by $Z\rightarrow b\bar{b}$ and $H\rightarrow b\bar{b}$ peaks have a partial overlap, in the 1-lepton case the excess must be present in both peaks, while in the 0 and 2-lepton final states it should appear only in the Higgs peak ($Z^{\prime}$ does not couple to $ZZ$). In addition, the excess in the 1-lepton category should be more significant, since in this case both $W^{+}$ and $W^{-}$ contribute. A combined analysis of the three final states can quantify the excess and correlate it to events with $b\bar{b}$ pairs of higher $p_{T}$ than in the SM Higgs production, since these pairs originate from heavy TeV-scale objects.

ATLAS and CMS recently published results in the search for the SM $H\rightarrow b\bar{b}$ in association with a gauge boson decaying leptonically ATLASSMvh ; CMSSMvh . This final state is identical with the one we propose here having sensitivity to HVT resonances. Both experiments report the measured $M_{bb}$ invariant mass as shown in Figures 2 and 3, for an integrated luminosity of $\sim 80$ fb*-1* at $13$ TeV.

Broad heavy resonances mainly decaying to dibosons and pairs of fermions, occur for larger values of the $g_{V}$ coupling ($g_{V}>3$ for $c_{H}=-1$ and $c_{F}=1$). The resonance width as a function of $g_{V}$ is shown in Fig. 4 and the branching ratios to a pair of fermions and diboson final states are shown in Fig. 5. In both figures the calculations shown were performed by implementing the HVT model in MadGraph5 madgraph5 . Fig. 4 shows that for $g_{V}>3$ the resonance natural width over the mass ratio exceeds $\sim 3\%$, and as it can be seen in Fig. 5, the resonance decay to dibosons ($VH$, $VV$) is dominant. Direct searches for resonances in the $VV$ and $VH$ channels assume narrow resonance width that corresponds to $g_{V}\leq 3$, leaving unexcluded a large part of the HVT model parameter space.

The magnitude of the impact of these resonances on SM measurements depends on their cross section which in turn, depends on the mass and the $g_{V}$ coupling. Heavy resonance Drell-Yan and Vector Boson Fusion (VBF) production cross sections for masses in the range 1-3 TeV are shown in Fig. 6. The DY mode is always dominant up to large values of $g_{V}$ close to the perturbative limit $g_{V}\simeq 4\pi$, and the total cross section drops with $g_{V}$. In Fig. 6 we can see that in the mass range $1-2$ TeV, the cross sections for $g_{V}>3$ range from $O(1)$ to $O(100)$ fb, leading to non-vanishing contributions at present and future measurements. For parts of the parameter space corresponding to low resonance masses and high $g_{V}$ that are theoretically excluded (parts for which the input electroweak parameters $\alpha_{EW}$, $G_{F}$ and $M_{Z}$, are not reproduced by the HVT model), the cross section is not provided in Fig. 6.

3 Resonance Search Potential at Colliders

To evaluate the expected signal yields for a broad range of masses and widths, we rely on Monte Carlo samples generated with MadGraph5, madgraph5 , interfaced to PYTHIA Pythia8:2008 . To simulate the response of an LHC-like experiment, realistic resolution and reconstruction efficiencies for electrons, muons, photons and jets were applied with the Delphes framework delphes . The specific cut-based event selection was based on the recent $VH$ with $H\rightarrow b\bar{b}$ results from ATLAS and CMS ATLASSMvh ; CMSSMvh . As discussed later, the efficiencies were normalized to the experimental ones. The selected events were split in three categories according to the final state: two $b$ jets with 0, 1 or 2 leptons. We follow more closely the ATLAS selection including the requirement of 2 or more jets out of which exactly two must be $b$ jets. We will call this SM selection, the baseline selection, in order to differentiate from additional discriminants which enhance the BSM Higgs signal. The BSM Higgs signal decays to a $b\bar{b}$ pair with a transverse momentum, $p_{\perp}$, significantly larger than the SM Higgs. For this reason, an additional requirement applied to the baseline analysis, is a cut on the transverse momentum of the $b\bar{b}$ system. As in the LHC analyses, we only consider resolved $b\bar{b}$ pairs although the search can be extended to a single, fat $b\bar{b}$ jet analysis. As a discriminant, the $b\bar{b}$ invariant mass $M_{bb}$ is used, which after full selection and subtraction of the $b\bar{b}$ continuum shows two peaks due to the presence of $Z\rightarrow b\bar{b}$ and $H\rightarrow b\bar{b}$.

The characteristic of the BSM signal is that the non-SM produced $H\rightarrow b\bar{b}$ has high $p_{\perp}$. This can be seen in Fig. 7, where the exotic Higgs has typical $p_{\perp}>150$ GeV while the SM Higgs is softer.

The expected signal and resonant SM background yields for the baseline selection are presented in Table 1 for the 0, 1 and 2-lepton final states, for a resonance of mass 1.5 TeV, $g_{V}=5$ and an integrated luminosity of 100 fb*-1*. For a $p_{\perp}>200$ GeV cut, we observe a significant increase in the fraction of BSM signal over the SM resonant yields, that allows to estimate the discovery potential as a function of the HVT model parameters.

In the calculations shown in Table 1, in order to realistically estimate the heavy resonance parameter space that can be probed at LHC as a function of luminosity, we have normalized the efficiencies to the published results by the LHC experiments. In this way, we maintain realistic non-resonant background levels and to some extent take systematic differences between the full simulation used by the experiments and the fast simulation used here partially into account. The signal efficiency $\times$ acceptance is at the level of few % and based on the signal cross-section expectations summarized in Fig. 6, observable excesses in $M_{bb}$ should be expected for broad resonance masses below 2 TeV, even for present LHC luminosities.

Example scenarios with parameters $M_{Z^{\prime}}=1.5$ TeV, $g_{V}=4$ and $g_{V}=5$ after baseline selection and a tight transverse momentum cut $p_{\perp}>300$ GeV, for which the bosonic decay branching ratio is dominant, are presented in Figures 8 and 9 respectively. In these Figures, the combined $M_{bb}$ distribution is shown for an integrated luminosity of 200 fb*-1*. Assuming that the new physics does not modify significantly the SM Higgs couplings, the excess over the SM expectation will indicate the presence of heavy vector resonances decaying into $VH$ and $WZ$ final states. As we have already seen from Table 1, most of the Higgs excess comes from the 1-lepton category, while all of the $ZW$ excess comes from the 1-lepton category but it is rather small.

The resonance search potential for a single LHC experiment and an integrated luminosity of 200 fb*-1*, is shown in Figures 10, 11, and  12, for the 0-lepton, 1-lepton and 2-lepton categories, respectively. As seen from these results, the search is sensitive to part of the phase space not currently probed by existing searches. For 200 fb*-1*, the sensitivity extends to masses up to about 2 TeV, and resonance widths $\Gamma/m\simeq$ 10%, although the higher the width the lower the resonance mass reach. For $\Gamma/m\simeq$ 10% only a mass up to 1.6 TeV can be reached. The reason for the search not being sensitive to larger widths is not only due to the reduction of cross section with $g_{V}$, but also due to a loss of efficiency because for $V^{\prime}$ masses greater than 1.5 GeV, the majority of $b$ quark pairs fall in the same jet. More precisely, for resonance masses of 2 TeV or higher, the fraction of resolved $b\bar{b}$ pairs is 5% or less, respectively. Here it should be stressed that an extension of the proposed search to include fat jets could push these limits at higher values.

In an attempt to test the prediction of the model against the $\sim$80 fb*-1* data recently reported by ATLAS and CMS, we estimate the observed limits as shown in Figures 13, 14, and 15. The observed limits are based on the signal strength values reported by the experiments in each category. The bands corresponding to the uncertainties in the Higgs signal strength $\mu$ are not shown. The observed limits should be compared to the expected limits with $p_{\perp}>200$ GeV. As seen in these results, the new LHC data already exclude new parts of the allowed resonance mass-width phase space. For instance from Fig. 14, HVT resonances with $\Gamma/m=9\%$ and mass up to 1.5 TeV, are excluded by the ATLAS experiment. The observed limits in the 1-lepton and 2-lepton categories exclude lower masses and widths than expected because of the corresponding excesses observed by the experiments. However, the uncertainties in each category are still large and more data are needed that will first settle if there is indeed an excess, and then allow to examine the compatibility of the pattern of the excess with the HVT model. It should be pointed out that the signal cross-sections were calculated at leading order, thus an additional theory uncertainty due to missing higher order corrections should be considered. These theoretical uncertainties are expected to be small with respect to the rest of the systematic uncertainties of this analysis.

4 Summary and Conclusions

In this work we explored the potential of a novel search of heavy vector resonances of non-zero natural width, decaying almost exclusively in the $H\rightarrow b\bar{b}$ and $Z\rightarrow b\bar{b}$ final states in association with a leptonically decaying $V$ ($Z$ or $W$) and $W$-only, respectively. For large $g_{V}$ coupling of the exotic resonance to the $W$ and $Z$ bosons and the SM Higgs, the branching ratio to $VH$ and to two gauge bosons $VV$ dominates, leading to a simultaneous excess to both Higgs $WH\rightarrow\ell\nu b\bar{b}$, $ZH\rightarrow\ell\ell b\bar{b}$ and non-Higgs $WZ\rightarrow\ell\nu b\bar{b}$ final states.

We showed that excesses of varying strengths should be observed in all final states (0, 1 or 2 leptons). The fact that the relative strengths of these excesses depend on branching ratios and efficiencies, provides a clear signature of the presence of heavy resonances or even their low mass tails. For a luminosity accessible by the LHC experiments of 200 fb*-1*, the search is sensitive up to resonance masses of 2 TeV and widths $\Gamma/m$ of 10%. A first test of a heavy vector triplet model against ATLAS and CMS data in terms of expected and observed limits was presented. Although experimental uncertainties on the $VH\rightarrow b\bar{b}$ signal strength $\mu_{VH}$ are still large, the LHC experiments are expected to accumulate significant amounts of data in 2018 and beyond, making the search proposed here a very useful tool in evaluating potential excesses in Higgs yield measurements when the Higgs is produced in association with a weak vector boson.

Acknowledgements

This work was supported by the Taiwanese Ministry of Science and Technology under grant number 106-2112-M-002-011-MY3.