Search for WW / W Z resonance production in (cid:96)ν qq ﬁnal states in pp collisions at √ s = 13 TeV with the ATLAS detector The

A search is conducted for new resonances decaying into a WW or WZ boson pair, where one W boson decays leptonically and the other W or Z boson decays hadronically. It is based on proton–proton collision data with an integrated luminosity of 36.1 fb − 1 collected with the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of √ s = 13 TeV in 2015 and 2016. The search is sensitive to diboson resonance production via vector– boson fusion as well as quark–antiquark annihilation and gluon–gluon fusion mechanisms. No signiﬁcant excess of events is observed with respect to the Standard Model backgrounds. Several benchmark models are used to interpret the results. Limits on the production cross section are set for a new narrow scalar resonance, a new heavy vector-boson and a spin-2 Kaluza–Klein graviton. A search is conducted for resonant WW and WZ production decaying into semileptonic ( (cid:96)ν qq ) ﬁnal states using 36.1 fb − 1 of pp collision data collected at a centre-of-mass energy of √ s = 13 TeV by the ATLAS detector at the LHC during 2015 and 2016. The analysis is carried out in two di ﬀ erent kinematic topologies of the hadronically decaying W / Z boson, which can be reconstructed either as two small- R jets or one large- R jet. The data are compatible with the Standard Model background hypothesis and the largest local excess observed is approximately 2.7 σ , which is not signiﬁcant. Limits on the production cross section are obtained as a function of the resonance mass for models predicting a narrow scalar boson, a heavy spin-1 vector boson and a spin-2 KK graviton. Two di ﬀ erent production modes are considered, the vector-boson fusion and the gluon–gluon fusion or quark–antiquark annihilation, and independent limits are set. Masses below 2730 GeV and 3000 GeV are excluded at 95% CL for the Z (cid:48) in models A and B of the HVT parametrization, respectively. For the W (cid:48) resonance, the corresponding limits obtained exclude masses below 2800 GeV and 2990 GeV. Additionally, RS G KK signals with k / ¯ M Pl = 1.0 produced via gluon–gluon fusion are excluded at 95% CL below 1750 GeV. This search has signiﬁcantly extended previous ATLAS high-mass limits [11], by 390-660 GeV, depending on the model.


Introduction
Diboson resonances are predicted in a number of extensions to the Standard Model (SM), such as composite Higgs models [1,2], warped extra dimensions [3][4][5], models with an extended Higgs sector [6,7] and grand unified theories [8][9][10]. Searches for diboson resonances in various decay channels have been carried out by the ATLAS and CMS collaborations at the Large Hadron Collider (LHC), but no evidence of such resonances has been observed [11][12][13][14][15][16][17][18]. The most recent ATLAS searches using data collected in 2015 and 2016 have been performed in the ZZ/ZW final state [17], with one Z decaying to leptons, and the fully hadronic final state with boson-tagged jets [18].
This paper reports on a search for a charged or neutral resonance, in a mass range from 300 GeV to 5000 GeV, that decays into a WZ or WW boson pair. The semileptonic final state where one W boson decays leptonically (W → ν with = e, µ) and the other W/Z boson (denoted by V) decays hadronically (V → qq /qq with q, q quarks) is considered. The search uses pp collision data at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 36.1 fb −1 , collected by the ATLAS experiment in 2015 and 2016. The strategy for identification of resonances depends on the ability to resolve the quarks from the hadronically decaying V boson. For high-mass resonances, the opening angles between the quarks from V boson decays are small and both quarks can be identified as a single jet. This case is referred to as the merged analysis and is denoted by νJ. In contrast, separate identification of the two quarks from low-mass resonances is referred to as the resolved analysis and is denoted by ν j j.
In addition to a larger data sample, the search makes use of several improvements to the methodology compared to the previous ATLAS result [11]. The resolved analysis has been included, and in addition the event selections are optimized for two different production modes: the vector-boson fusion (VBF) and the gluon-gluon fusion (ggF) or quark-antiquark (qq) annihilation. In addition, a new mass reconstruction algorithm is implemented for hadronically decaying W/Z bosons that are highly Lorentz boosted. It is based on both the calorimeter energy deposits and the charged tracks instead of calorimeter information alone, as used in the previous publication [11].
The VBF process (pp → VV j j) is characterized by the presence of two jets with a large rapidity gap resulting from quarks from which a vector boson is radiated. The absence of this topology is interpreted as ggF or qq production, collectively referred to as ggF/qq in this paper. Results are provided for the VBF and ggF/qq categories separately and possible signal leakage between categories is neglected.
The spectrum of the reconstructed invariant mass of the WV resonance candidates, m(WV), is examined for localized excesses over the expected SM background. Three signal models are used to optimize the event selection, assess the sensitivity of the search and interpret the data: an additional heavy Higgs boson predicted by many theories beyond the SM, a heavy vector triplet (HVT) parameterization based on a simplified phenomenological Lagrangian [19, 20] and a bulk Randall-Sundrum (RS) model [21]. transverse energy of the electron itself and correcting for the expected pile-up contributions. The isolation requirement for electrons is chosen to ensure approximately 95% and 99% selection efficiency, for signal and veto electrons, respectively.
Muons are reconstructed by combining an ID track with an MS track that has compatible trajectory [63]. Based on the quality of their reconstruction and identification, signal muons are required to pass the "medium" selection with p T > 27 GeV and |η| < 2.5, while veto muons are required to pass the "loose" selection, p T > 7 GeV and |η| < 2.7. In addition, a similar isolation requirement to that used for electron candidates, only using tracks within a cone of ∆R = 0.3, is applied to signal and veto muon candidates with an efficiency of 99%.
To reject non-prompt leptons, requirements of |d 0 |/σ d 0 < 5 (3) and |z 0 sin θ| < 0.5 mm are imposed on the tracks associated with the electrons (muons), where d 0 is the transverse impact parameter with respect to the measured beam-line position, σ d 0 is the corresponding uncertainty, z 0 is the longitudinal impact parameter with respect to the primary vertex and θ is the polar angle of the track. 2 Jets are reconstructed using the anti-k t algorithm [64] implemented in the FastJet package [65] from three-dimensional topological clusters of energy deposits in the calorimeter [66], with two different radius parameters: R = 1.0 for large-R jets (denoted by J) and R = 0.4 for small-R jets.
Small-R jets [67] are required to have p T > 20 GeV and |η| < 2.4, while jets considered for the tagging of VBF events are required to have p T > 30 GeV and |η| < 4.5. For jets with p T < 60 GeV and |η| < 2.4 a jet-vertex-tagger multivariate discriminant [68], based on tracking and vertexing information, is applied to select jets that originate from the primary vertex. The selected working point provides at least 92% efficiency.
An overlap removal procedure is applied to prevent using the same energy deposits in more than one electron, muon or jet. Small-R jets are discarded if they are within a cone of size ∆R = 0.2 around the direction of an electron candidate. However, if the distance between a jet and an electron candidate is within 0.2 < ∆R < min(0.4, 0.04 + 10/p e T ), the jet is retained but the nearby electron is rejected from the analysis. A muon candidate lying within ∆R < min(0.4, 0.04 + 10/p µ T ) from a small-R jet is discarded unless it is within ∆R < 0.2 and satisfies one of the two following requirements: (a) the small-R jet has fewer than three tracks; (b) p µ T /p j T > 0.5 and p µ T / p T > 0.7, where p T is the sum of the transverse momenta of tracks associated with the small-R jet. In this case, the muon is retained but the nearby small-R jet is rejected.
Small-R jets containing b-hadrons are identified using the MV2c10 b-tagging algorithm [69, 70] with an efficiency of 85%, determined with tt simulated events. The corresponding misidentification rates are approximately 3% and 30% for selecting jets originating from light quark and charm quark, respectively. For simulated samples the b-tagging efficiencies are corrected to match those measured in data [69].
Large-R jets [71, 72] are formed from constituent energy deposits and are trimmed to mitigate pile-up effects and soft radiation. The jet constituents are reclustered into subjets using the k t algorithm with R = 0.2 [73], removing those which carry less than 5% of the p T of the original jet [74]. To overcome the limited angular resolution of the calorimeter, the mass of a large-R jet is computed using a combination of calorimeter and tracking information [75]. The mass is defined as: In the merged analysis, a baseline selection on the D 2 variable providing 80% efficiency for V signals is applied to all large-R jets. To further distinguish hadronically decaying V bosons from jets originating from non-top quarks or gluons, boson tagging algorithms (V-tagging) based on the combined large-R jet mass and the D 2 variable are constructed. The requirements on the D 2 variable and the mass window depend on the jet p T and are defined separately for the W and Z bosons. In this paper, working points resulting in 50% and 80% signal selection efficiency are used, as defined in Sec. 5.
The missing transverse momentum (E miss T ) is the absolute value of the negative vectorial sum of the transverse momenta of electrons, muons, and small-R jets. Reconstructed charged-particle tracks originating from the primary vertex and not matched to any electron, muon, or jet are also included in the E miss T reconstruction [80].
The neutrino momentum four-vector is reconstructed by imposing a W boson mass constraint on the charged-lepton-neutrino system. The neutrino transverse momentum components are set equal to the missing transverse momentum of the event and the unknown z-component of the momentum (p z ) is obtained from the resulting quadratic equation. The p z is chosen as either the smaller, in absolute value, of the two real solutions or, if the solution is complex, its real part.
The selection criteria outlined above are the result of a signature-dependent optimisation using the asymptotic significance.

Trigger and event selection
Events are selected that contain exactly one charged signal lepton and no additional veto electrons or muons. Single-electron triggers with minimum transverse energy (E T ) thresholds of 24 GeV and 26 GeV in 2015 and 2016, as well as 60 GeV are applied to record events in the electron final state. The low threshold triggers require electron candidates to pass isolation requirements resulting in at least 90% efficiency, depending on the lepton p T . As for the muon final state, the events are recorded either by a single-muon trigger or an E miss T triggers. The single-muon trigger, with p T > 20 (26) GeV in 2015 (2016), 3 The angular exponent β, defined in Ref. [77], is set to unity. 6 is subject to a large inefficiency due to limited trigger hardware coverage. The E miss T trigger has an online threshold of 70 GeV for the 2015 data and of 90-110 GeV for the 2016 data, where the muon track p T is not used to compute E miss T in the trigger algorithm. Therefore, it is fully efficient for W → µν with p T (W) > 200 GeV and it is used in the merged analysis, where a high-p T lepton is expected, to recover the single-muon trigger inefficiency. Events recorded by single-lepton triggers, where the signal lepton matches the trigger lepton, and E miss T triggers are selected.
The sensitivity to resonances of different masses is optimized by classifying the events according to the topology, production mechanism and amount of background. The event selection criteria are summarized in Tables 1 and 2 for the merged and resolved analyses respectively. Figure 1 illustrates the jet selections used to reconstruct the hadronically decaying V boson candidates in the signal and control regions of the analysis. The mass of either the large-R jet (m J ) or the system of two small-R jets (m j j ) is used to define "mass windows".  Figure 1: (a) Illustration of the merged WW (shaded area) and WZ (dashed lines) signal regions (SR) according to the large-R jets selection. The 50% and 80% V-tagging efficiency (ε V ) working points, based on the combined cut of the D 2 and m J , are used to form the high-purity (HP) and low-purity (LP) regions respectively. For each working point, a jet mass requirement is imposed and an upper bound on the substructure variable is set. Since both requirements depend on the p T of the large-R jet, an absolute definition is not given in the figure. (b) Definitions of the resolved WW and WZ SR based on the dijet mass selection. In both channels, the SR mass sidebands are used to define the W+jets control region (CR).
The unique kinematic signature of the VBF process is used to define event categories enriched in this production mechanism and maximize the sensitivity by reducing the SM backgrounds. Events with two small-R ("tag") jets with invariant mass m tag ( j, j) > 770 GeV and pseudorapidity gap between them |∆η tag ( j, j)| > 4.7 are classified as VBF candidates. In case there are more than two tag-jets, the pair with the largest invariant mass is chosen. Events that fail the VBF selection are assigned to the ggF/qq category.
Events belonging to the VBF or ggF/qq categories are further assigned to the merged or resolved regions as follows: • Merged signal region: the large-R jet with the highest p T is selected as the candidate for the hadronically decaying V boson, requiring no overlap with either of the tag-jets in the VBF category (∆R( j tag , J) > 1.0). Furthermore, the event is required to have E miss T > 100 GeV to suppress the multijet contamination. The leptonically decaying W candidate is required to have a leptonneutrino system with transverse momentum p T ( ν) > 200 GeV. A threshold of 0.2 is set on the ratio E miss T /p T (eν) in the electron channel in order to further suppress the multijet background. In the desired signal topology, the two bosons are produced from a heavy resonance decay and their transverse momenta are expected to be close to half the reconstructed resonance mass. As a result, a threshold of 0.4 (0.3) is applied to p T (J)/m(WV) and p T ( ν)/m(WV) in the ggF/qq (VBF) category. Furthermore, events are rejected if there is a b-tagged jet present with a separation of ∆R > 1.0 from the hadronically decaying V candidate. The latter requirement rejects more than 70% of background events from tt production while keeping more than 95% of signal events, independently of the resonance mass. The remaining events are assigned to the high-purity (HP) region if the large-R jet satisfies the V-tagging 50% efficiency working point, for both the mass window and the D 2 variable, as defined in Section 4. Otherwise, events are assigned to the low-purity (LP) region if the 80% efficiency working point is satisfied for the large-R jet. The improvement in cross-section sensitivity resulting from combining the HP and LP regions reaches up to 36% for resonances with 5.0 TeV mass. The selected HP and LP events can simultaneously pass both the WW and the WZ selections if the large-R jet passes both the W and Z selections.
• Resolved signal region: events not satisfying the selection criteria of the merged signal region and with E miss T > 60 GeV and p T ( ν) > 75 GeV are considered. The hadronically decaying V candidate is formed by combining the two small-R jets, excluding VBF tag-jets, with the highest p T and requiring their invariant mass to be between 66 and 94 (82 and 106) GeV in order to be consistent with the W (Z) boson mass. The two selected small-R jets are required to have p T > 45 GeV (60 GeV for the highest p T jet) and the azimuthal angle separation between jets, lepton and E miss T directions must satisfy ∆φ( j, ) > 1.0, ∆φ( j, E miss T ) > 1.0, ∆φ( j, j) < 1.5 and ∆φ( , E miss T ) < 1.5. In the calculation of the WV invariant mass, a V mass constraint is imposed on the two small-R jets by Table 2: Summary of the selection criteria in the resolved analysis for the WW and WZ signal regions (SR), W+jets control region (W CR) and tt control region (tt CR). The events are also categorized according to their production mechanism, the VBF selection is prioritized and the remaining events are assigned to the ggF/qq category. [ rescaling the p T of the dijet system to be p j j T and m( j j) are the transverse momentum and the invariant mass of the dijet system respectively, and m(V) is the known value of the V boson mass. Studies using MC simulated events show that the mass constraint reduces the uncertainties due to the jet energy scale and results in an approximately 20% improvement of the resolution of the reconstructed diboson resonance mass, which ranges between 20 GeV and 120 GeV across the mass spectrum. In addition, selected events in the ggF/qq (VBF) category are required to satisfy p T ( j j)/m(WV) > 0.35 (0.3) and p T ( ν)/m(WV) > 0.35 (0.3). Events are rejected from the WW selection if both jets from the V boson decay are tagged as b-tagged jets. Furthermore, events with one or more b-tagged jets, not compatible with the V boson decay, are also removed. As in the merged signal region, a threshold of 0.2 is set on the ratio E miss T /p T (eν) to suppress the multijet background in the eνqq channel.
The signal efficiency times acceptance ( × A), defined as the ratio of the number of signal events in the signal region to the number of generated signal events, is presented as a function of the WV → νqq resonance mass in Figures 2 and 3 for all the generated benchmark signals. Experimental factors, such as the detector coverage and the pile-up activity, lead to low tagging efficiency of the VBF jets resulting in small × A. Priority is given to the VBF category, using the selection outlined previously, aiming to increase the sensitivity to genuine VBF signal events that have a small signal × A in the VBF category and the high fraction of the VBF signal that leaks in the ggF/qq category. The leakage occurs due to inefficiencies related to the reconstruction and identification of the "tag" jets, and results in a small deterioration in sensitivity after accounting for the background. Concerning the × A of the various analyses, the resolved analysis is more sensitive in the low mass region, while the merged analysis is more efficient in the high mass region with a relatively constant × A. In the ggF/qq category, the × A values are generally lower for the scalar signal because the two bosons are produced less centrally than for the spin-1 and spin-2 signals, and the p T (V)/m(WV) requirements reject more signal.

Background estimation
Simulation studies indicate that the dominant background sources are W+jets and tt events. The W+jets contribution is found to be approximately 50%, 70% and 60%-65% in the high-purity, low-purity and resolved ggF/qq signal regions, respectively, while the corresponding fractions in the VBF category are 40%, 60% and 40%-55%. In the resolved analysis, the W+jets contribution is higher in the WZ channel than the WW channel because of the different selections on b-jets. The tt contamination in the ggF/qq category is estimated to be 30% (20%) in the high-purity (low-purity) and 25% (30%) in the resolved WW (WZ) signal regions. The contribution from tt production in the VBF category is 50%, 30% and 35%-50% in the high-purity, low-purity and the resolved signal regions, respectively. Smaller background contributions arise from Z+jets, single-top and SM diboson production. Control regions for the high-and low-purity categories as well as the resolved category are defined for events that fail the selection criteria of the signal regions in order to estimate the dominant background contributions: • The W+jets control regions are formed from events satisfying the signal region selection except for the invariant mass requirement of the hadronically decaying V candidate. The mass is required to be in the sideband region which is defined as m( j j) < 66 GeV or 106 < m( j j) < 200 GeV for the resolved analysis. In the merged analysis, the sideband regions are formed by events satisfying the respective D 2 selections but not the mass window requirement for the 80% efficiency working point. Approximately 65% and 77% of the selected events are from W+jets production in the ggF/qq category of the merged and resolved analyses, respectively. The remaining events are primarily from tt production. The contribution from W+jets processes is 50% and 65% for the merged and resolved analyses, respectively, in the VBF category.
• The tt control regions are formed from events satisfying the signal region selection except for the b-jet requirement, which is inverted. Studies using simulated events show that 77%-87% of the selected ggF/qq and VBF category events are from tt production and the rest are from single-top, V+jets or diboson production, for both the merged and the resolved event topologies.
The shapes of the mass distributions for events from production of W+jets and tt are modelled using simulated events. Their normalizations are determined from a combined fit to the events in the signal and control regions, as detailed in Sec. 8. Concerning the subdominant background contributions from Z+jets, single-top and SM diboson production, simulation is used to obtain the shapes and normalizations, which are subsequently constrained within statistical, experimental and theoretical uncertainties.
The contribution from multijet production primarily consists of events with jets or photon conversions misidentified as leptons or real but non-prompt leptons from decays of heavy-flavour hadrons. The multijet background in the merged event topology is estimated by a fit to the E miss T distribution of events that satisfy all the signal selection criteria but without any E miss T requirement. The shape of multijet events is obtained from an independent data control sample that satisfies the signal selection criteria except for the E miss T requirement and the lepton requirement: the leptons are required to satisfy the veto lepton selection, defined in Section 4, but not the signal lepton selection. Contributions from other processes with prompt leptons to the control sample are subtracted from the data using samples of simulated events in the extraction of the multijet background shape. In the fit, the normalizations of the W+jets and multijet components are allowed to float, with all the other backgrounds fixed to their predicted cross sections. Following this procedure, the multijet background in the merged event topology is found to be negligible.
A fake-factor method is implemented to estimate the multijet background contribution in the resolved topology. The "signal lepton" control region is formed by events that have exactly one signal lepton and exactly one small-R jet. The same event selection criteria are applied to the events in the "inverted lepton" control region except for the lepton requirement: the selected electron candidate is required to pass the "medium" but fail the "tight" requirements, and the selected muon candidate is required to fail the nominal but pass a looser isolation requirement. The fake-factor is defined as the ratio of the number of events in the signal lepton control region to the number of events in the inverted lepton control region, after subtracting contributions from prompt leptons as estimated by the simulation. The fake-factor is calculated as a function of the lepton p T and η, and E miss T . It is subsequently used to reweight a sample of events selected with the inverted lepton selection, as previously described, that satisfy the rest of the signal region selection.

Systematic uncertainties
Systematic uncertainty sources impacting the search can be divided into four categories: experimental uncertainties related to the detector or reconstruction algorithms, uncertainties in the estimations of background contributions, uncertainties in modelling the signal and statistical uncertainties in the MC predictions. Two kinds of background uncertainties are provided, normalization and shape uncertainties. Normalization uncertainties are extracted from data and MC simulation comparisons, while shape uncertainties are accounted for by varying MC parameters.
Modelling uncertainties affecting the shape of the final mass discriminant are estimated for the W+jets background. These include uncertainties in the renormalization and factorization scales, the CKKW [81,82] matching scales, the resummation scale, the PDF and α s . The scale uncertainties are obtained by doubling and halving the corresponding parameters in the nominal generator. Potential systematic uncertainties due to choices of parton shower and matrix element implementations are estimated by comparing the nominal MC samples to the alternative samples generated using MadGraph.
The uncertainty in the shape of the m(WV) distribution from the tt background is estimated by comparing the nominal Powheg+Pythia sample to the alternative samples described in Section 3. The factorization and renormalization scales of the nominal generator are varied, in a similar manner as the W+jets parameters, and their difference from the nominal sample is also applied as a systematic uncertainty.
The SM diboson production cross section is fixed to the inclusive next-to-leading-order calculation with a 30% systematic uncertainty in the normalization. The m(WV) distribution shape uncertainty of the diboson background is estimated by comparing the predictions based on the alternative Powheg-Box MC samples to those of the nominal Sherpa MC samples.
Systematic uncertainties in the multijet background estimate are only considered in the resolved analysis, as this background contribution in the merged analysis is negligible. These are obtained by varying the lepton or isolation selection used in the fake-factor calculation. In addition, the statistical uncertainties of the measured fake-factors and the systematic uncertainties in the prompt lepton contribution in the measurement of the fake-factors, are taken into account in the estimation of systematic uncertainties of the multijet background modelling. The effect of this uncertainty is found to be marginal in the fit.
Experimental uncertainties related to leptons, jets and E miss T are considered, affecting the shape and normalization of both the background and the signal distributions. These are estimated for the trigger efficiencies, the energy scale and resolution of small-R jets [67] and large-R jets [71], lepton identification, reconstruction and isolation efficiencies, lepton momentum scales and resolutions [61][62][63], b-tagging efficiency and misidentification rates [69, 70], and missing transverse momentum resolution [80].
For central small-R jets (|η| < 2.0), the total relative uncertainty in the jet energy scale [67] ranges from about 6% for jets with p T of 25 GeV to about 2% for p T of 1000 GeV. The uncertainty in the small-R jet energy resolution ranges from 10%-20% for jets with p T of 20 GeV to less than 5% for jets with p T > 200 GeV.
The uncertainties in the scale of the D 2 variable and in the large-R jet energy and mass are estimated by comparing the ratio of calorimeter-based to track-based energy and mass measurements in dijet data and simulation [71]. These uncertainties range between 2% and 5%. An uncertainty of 2% is assigned to the large-R jet energy resolution and uncertainties of 20% and 15% are assigned to the resolution of the large-R jet mass and D 2 , respectively.

13
The dominant uncertainties in the signal acceptance arise from the choice of PDF and the uncertainty in the amount of initial-and final-state radiation (ISR and FSR, respectively) in simulated events. The cross section obtained with the nominal PDF set is compared to those of the MMHT 2014 PDF [83] and CT14 PDF [84] to derive the uncertainties in the acceptance. The prescription in Ref. [85] is followed and the envelope of the uncertainties associated to the three PDF sets is used. The ISR/FSR contributions are computed by varying the parton shower and multi-parton interaction parameters following the prescription in Ref. [41].
The uncertainty in the combined 2015+2016 integrated luminosity is 3.2%. It is derived, following a methodology similar to that detailed in Ref. [86], from a preliminary calibration of the luminosity scale using x-y beam-separation scans performed in August 2015 and May 2016. This uncertainty is applied to the yields predicted by the simulation.

Results
The results are extracted by performing a simultaneous binned maximum-likelihood fit to the m(WV) distributions in the signal regions and the W+jets and tt control regions. The WW and WZ channels are treated individually, without combining their respective regions. A test statistic based on the profile likelihood ratio [87] is used to test hypothesized values of the global signal-strength factor (µ), separately for each model considered. The likelihood is defined as the product of the Poisson likelihoods for all signal and control regions for a given production mechanism category and channel (WW or WZ), simultaneously for the electron and muon channels. The fit includes six contributions, corresponding to W+jets, tt, single-top, Z+jets, diboson and multijet events. The main background sources, namely W+jets and tt, are constrained by the corresponding control regions and are treated as uncorrelated among the resolved and merged signal regions. For each of these backgrounds, a normalization factor, defined as the ratio of the number of simulated events after the fit to the number of simulated events before the fit, is derived and the results are collectively presented in Table 3. In all regions and categories, the normalization factors are found to be compatible with 1.0.
Systematic uncertainties are taken into account as constrained nuisance parameters with Gaussian or lognormal distributions. For each source of systematic uncertainty, the correlations across bins of m(WV) distributions and between different kinematic regions, as well as those between signal and background, are taken into account. The number of bins and the bin widths in each signal region are optimized according to the expected background event distribution and detector resolution. In the merged region, the diboson invariant mass range extends from 500 GeV to 5000 GeV divided into twenty (eleven) bins in the ggF/qq (VBF) category. The resolved region is covered by ten (nine) bins of varying width in the ggF/qq (VBF) category, beginning at 300 GeV and ending at 1500 GeV, due to the selection efficiency. In all regions, the overflow events are included in the last bin. The m(WV) distributions are presented in Figures 4 and  5 after the VBF and ggF/qq categorizations, respectively, for the merged and the resolved regions.
The list of leading sources of uncertainty in the best-fit µ value is given in Table 4 together with their relative importance (∆µ/µ). The values are quoted separately for the VBF and ggF/qq categories, and for the case of high and low mass signal samples, for which the merged and resolved topologies reach the highest sensitivity respectively. The largest systematic uncertainties are related to the background modelling and jet measurements and these are most important at lower masses. Exclusion limits are calculated using the CL s method [88], in the asymptotic approximation, at the 95% confidence level (CL) for resonance masses below 1.0 (1.6) TeV in the VBF (ggF/qq) category. For higher masses, the small number of expected events makes the asymptotic approximation imprecise and the limits are calculated using pseudo-experiments. The limits are calculated by fitting the merged highand low-purity signal regions simultaneously with the corresponding resolved region. The calculation is performed separately in each final state, WW or WZ, and the largest local excess observed is approximately 2.7 σ, which is not significant. The observed and expected upper limits on the cross sections for all generated benchmark signal models are shown in Figures 6 and 7 for the VBF and ggF/qq categories respectively. Because of the small deterioration in sensitivity after accounting for the background and the unknown ratio of the various production mechanisms in the models that are considered, the interpretation in the VBF (ggF/qq) category assumes there is no signal leakage from ggF/qq (VBF) processes. Table 5 summarizes exclusion limits on the mass for the various signal hypotheses as extracted from the ggF/qq category. For signal produced via the VBF mechanism and all scalar signals, only upper limits on the cross sections are set.

WW Selection Excluded
HVT

20
A search is conducted for resonant WW and WZ production decaying into semileptonic ( νqq) final states using 36.1 fb −1 of pp collision data collected at a centre-of-mass energy of √ s = 13 TeV by the ATLAS detector at the LHC during 2015 and 2016. The analysis is carried out in two different kinematic topologies of the hadronically decaying W/Z boson, which can be reconstructed either as two small-R jets or one large-R jet. The data are compatible with the Standard Model background hypothesis and the largest local excess observed is approximately 2.7 σ, which is not significant. Limits on the production cross section are obtained as a function of the resonance mass for models predicting a narrow scalar boson, a heavy spin-1 vector boson and a spin-2 KK graviton. Two different production modes are considered, the vector-boson fusion and the gluon-gluon fusion or quark-antiquark annihilation, and independent limits are set. Masses below 2730 GeV and 3000 GeV are excluded at 95% CL for the Z in models A and B of the HVT parametrization, respectively. For the W resonance, the corresponding limits obtained exclude masses below 2800 GeV and 2990 GeV. Additionally, RS G KK signals with k/M Pl = 1.0 produced via gluon-gluon fusion are excluded at 95% CL below 1750 GeV.