Search for production of $WW/WZ$ resonances decaying to a lepton, neutrino and jets in $pp$ collisions at $\sqrt{s}$ = 8 TeV with the ATLAS detector

A search is presented for narrow diboson resonances decaying to $WW$ or $WZ$ in the final state where one $W$ boson decays leptonically (to an electron or a muon plus a neutrino) and the other $W/Z$ boson decays hadronically. The analysis is performed using an integrated luminosity of 20.3 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 8 TeV collected by the ATLAS detector at the LHC. No evidence for resonant diboson production is observed, and resonance masses below 700 GeV and 1490 GeV are excluded at 95% confidence level for the spin-2 Randall-Sundrum bulk graviton $G^*$ with coupling constant of 1.0 and the extended gauge model $W'$ boson respectively.


Introduction
Several new physics scenarios beyond the Standard Model (SM), such as technicolour [1][2][3], warped extra dimensions [4][5][6], and grand unified theories [7], predict new particles that predominantly decay to a pair of on-shell gauge bosons.In this paper, a search for such particles in the form of WW/WZ resonances where one W boson decays leptonically (W → ν with = e, µ) and the other W/Z boson decays hadronically (W/Z → q q /q q, with q, q = u, c, d, s or b) is presented.This search makes use of jet-substructure techniques for highly boosted W/Z bosons decaying hadronically and is optimized to significantly improve the sensitivity to high mass resonances compared to previous searches.
Two benchmark signal models are used to optimize the analysis strategy and interpret the search results.A spin-2 Kaluza-Klein (KK) graviton (G * ) is used to model a narrow resonance decaying to a WW final state.The KK graviton interpretation is based on an extended Randall-Sundrum model of a warped extra dimension (RS1) [8] where the SM fields can propagate into the bulk of the extra dimension.This extended "bulk" RS model, referred to as bulk RS hereafter, avoids constraints on the original RS1 model from limits on flavour-changing neutral currents and electroweak precision tests, and has a dimensionless coupling constant k/ MPl ∼ 1, where k is the curvature of the warped extra dimension and MPl = M Pl / √ 8π is the reduced Planck mass.A spin-1 gauge boson (W ) of the Sequential Standard Model with modified coupling to WZ, also referred to as the extended gauge model (EGM) [7], is used to model a narrow resonance that decays to a WZ final state.The EGM introduces W and Z bosons with SM couplings to fermions and with the coupling strength of the heavy W to WZ modified by a mixing factor ξ = c × (m W /m W ) 2 relative to the SM couplings, where m W and m W are the pole masses of the W and W bosons respectively, and c is a coupling scaling factor.In this scenario the partial width of the W boson scales linearly with m W , leading to a narrow resonance over the accessible mass range.The width of the W resonance at 1 TeV is approximately 35 GeV.
Searches for these particles in several decay channels have been performed at the Tevatron and the Large Hadron Collider (LHC) and are reported elsewhere [9][10][11][12][13].Previous results from the ATLAS experiment in the q q channel excluded EGM W bosons with masses up to 1.59 TeV for WZ final states and RS1 gravitons with masses up to 740 GeV for ZZ final states [13].The CMS experiment set limits on the production cross sections of bulk RS gravitons as well as excluded RS1 gravitons with k/ MPl = 0.1 for masses up to 1.2 TeV and W bosons for masses up to 1.7 TeV [9].This analysis is based on pp collision data at a centre-of-mass energy √ s = 8 TeV corresponding to an integrated luminosity of 20.3 fb −1 collected by the ATLAS experiment at the LHC.The uncertainty on the integrated luminosity is ±2.8%.It is determined, following the same methodology as that detailed in Ref. [14], from a calibration of the luminosity scale derived from beam-separation scans performed in November 2012.

The ATLAS detector
The ATLAS detector [15] is a general-purpose particle detector used to investigate a broad range of physics processes.It includes inner tracking devices surrounded by a superconducting solenoid, electromagnetic and hadronic calorimeters and a muon spectrometer with a toroidal magnetic field.The inner Simulated event samples are used to model the shape and normalization of most SM background processes.The main background sources in the analysis arise from W bosons produced in association with jets (W+jets), followed by top-quark and multijet production, with smaller contributions from dibosons and Z+jets.Production of W and Z bosons in association with up to five jets is simulated using Sherpa 1.4.1 [21] with the CT10 PDFs [22], where b-and c-quarks are treated as massive particles.Samples generated with MC@NLO [23] and interfaced to Herwig [24] for hadronization and to Jimmy [25] for the underlying event are used for t t production as well as for single top-quark production in the schannel and the Wt process.The t t cross section is normalized to the calculation at NNLO in QCD including resummation of next-to-next-to-leading logarithmic soft gluon terms with Top++2.0 [26][27][28][29][30][31][32].Single top-quark production in the t-channel is simulated with AcerMC [33] interfaced to Pythia6 [34].Diboson samples (WW, WZ and ZZ) are generated with Herwig and Jimmy.
The effect of multiple pp interactions in the same and neighbouring bunch crossings (pile-up) is included by overlaying minimum-bias events simulated with Pythia8 on each generated signal and background event.The number of overlaid events is such that the distribution of the average number of interactions per pp bunch crossing in the simulation matches that observed in the data (on average 21 interactions per bunch crossing).The generated samples are processed through the Geant4-based detector simulation [35,36] or a fast simulation using a parameterization of the performance of the calorimeters and Geant4 for the other parts of the detector [37], and the standard ATLAS reconstruction software used for collision data.

Event selection
Events are required to have a vertex with at least three associated tracks, each with transverse momentum p T > 400 MeV.The primary vertex is chosen to be the reconstructed vertex with the largest track p 2 T .The main physics objects used in this analysis are electrons, muons, jets and missing transverse momentum.Electrons are selected from clusters of energy depositions in the calorimeter that match a track reconstructed in the ID and satisfy "tight" identification criteria defined in Ref. [38].The electrons are required to have transverse momentum p T > 25 GeV and |η| < 2.47, excluding the transition region between the barrel and endcaps in the LAr calorimeter (1.37 < |η| < 1.52).Muons are reconstructed by combining ID and MS tracks that have consistent trajectories and curvatures [39].The muon tracks are required to have p T > 25 GeV and |η| < 2.5.In addition, leptons are required to be isolated from other tracks and calorimetric activity.The scalar sum of transverse momenta of tracks with p T > 1 GeV within ∆R = (∆η) 2 + (∆φ) 2 = 0.2 around the lepton track is required to be less than 15% of the lepton p T .Similarly, the sum of transverse energy deposits in the calorimeter within a cone of ∆R = 0.2, excluding the transverse energy from the lepton and corrected for the expected pile-up contribution, is required to be less than 14% of the lepton p T .In order to ensure that leptons originate from the interaction point, a requirement of |d 0 |/σ d 0 < 6 (3.5) and |z 0 sin θ| < 0.5 mm is imposed on the electrons (muons), where d 0 (z 0 ) is the transverse (longitudinal) impact parameter of the lepton with respect to the reconstructed primary vertex and σ d 0 is the uncertainty on the measured d 0 .
In this analysis, jets are reconstructed from three-dimensional clusters of energy depositions in the calorimeter using two different algorithms.The jet constituents are considered massless.The low-p T hadronically decaying W/Z candidates are selected by combining the two highest-p T jets which are constructed by the anti-k t algorithm [40] with a distance parameter of R = 0.4.These jets are referred to as small-R jets and denoted by " j" hereafter.The energy of small-R jets is corrected for losses in passive material, the non-compensating response of the calorimeter, and extra energy due to multiple pp interactions [41].The small-R jets are required to have p T > 30 GeV and |η| < 2.8.For jets with p T < 50 GeV, the summed scalar p T of associated tracks from the reconstructed primary vertex is required to be at least 50% of the summed scalar p T of all associated tracks.In the pseudorapidity range |η| < 2.5, jets containing hadrons from b-quarks are identified using the MV1 b-tagging algorithm [42] with an efficiency of 70%, determined from t t simulated events, and with a misidentification rate for selecting light-quark or gluon jets of less than 1%.
For high-p T W/Z bosons, such as the ones from a resonance with mass above 1 TeV, the hadronically decaying W/Z candidates are identified using a single large-R jet, referred to as "J" hereafter.The Cambridge-Aachen jet clustering algorithm [43] with a distance parameter of R = 1.2 is used.This jet algorithm offers the advantage of allowing the usage of a splitting and filtering algorithm similar to that described in Ref. [44] but optimized for the identification of highly boosted boson decays.To exploit the characteristics of the decay of massive bosons into a light-quark pair, the splitting and filtering algorithm used here does not impose a mass relation between the large-R jet and its subjets [45].The momentum balance is defined as T )∆R 12 /m 12 , where p j1 T and p j2 T are the transverse momenta of the two leading subjets, ∆R 12 is their separation and m 12 is their invariant mass.To suppress jets from gluon radiation and splitting, √ y f is required to be greater than 0.45.Furthermore, the large-R jets are required to have p T > 400 GeV and |η| < 2.0.
The missing transverse momentum (with magnitude E miss T ) is calculated as the negative of the vectorial sum of the transverse momenta of all electrons, muons, and jets, as well as calibrated calorimeter energy clusters within |η| < 4.9 that are not associated with any other objects [46].
The data used were recorded by single-electron and single-muon triggers, which are fully efficient for leptons with p T > 25 GeV.The analysis selects events that contain exactly one reconstructed electron or muon matching a lepton trigger candidate, E miss T > 30 GeV and no b-tagged small-R jets.The transverse momentum of the neutrino from the leptonically decaying W boson is assumed to be equal to the missing transverse momentum.The momentum of the neutrino in the z-direction, p z , is obtained by imposing the W boson mass constraint on the lepton and neutrino system, which leads to a quadratic equation.The p z is defined as either the real component of the complex solution or the smaller in absolute value of the two real solutions.
In order to maximize the sensitivity to resonances with different masses, three different optimized sets of selection criteria are used to classify the events according to the p T of the leptonically decaying W candidate (p ν T ) and hadronically decaying W/Z candidate (p T or p J T ), namely the "low-p T resolved region" (LRR), "high-p T resolved region" (HRR) and "merged region" (MR), where the highly boosted W/Z decay products are observed as a single merged jet in the final state.To ensure the orthogonality of the signal regions, events are assigned exclusively to the first region for which the criteria are fulfilled, applying sequentially the MR, HRR, and LRR event selection.The hadronically decaying W/Z candidate is formed by combining the two small-R jets with highest p T in the resolved regions and its invariant mass m j j is required to be between 65 and 105 GeV.In the LRR (HRR), the event is required to have p ν T > 100 (300) GeV, p j j T > 100 (300) GeV and ∆φ( j, E miss T ) > 1, where ∆φ( j, E miss T ) is the azimuthal angle between the leading jet and the missing transverse momentum.The HRR additionally requires the two leading jets to have p T > 80 GeV.In the MR region, the large-R jet with the highest p T is selected as the hadronically decaying W/Z candidate and p ν T > 400 GeV is also imposed.The jet mass of the selected large-R jet (m J ) is required to be consistent with a W/Z boson mass (65 < m J < 105 GeV) and the azimuthal angle between the jet and the missing transverse momentum, ∆φ(J, E miss T ), is required to satisfy ∆φ(J, E miss T ) > 1.The signal acceptance times efficiency after all selection requirements increases from about 5% at m W = 300 GeV to a plateau of around 25% for m W > 500 GeV for W → WZ → νq q with = e, µ, τ.

Background estimation
The reconstructed WW/WZ mass, m ν j j (m νJ ), defined as the invariant mass of the ν j j ( νJ) system, is used to distinguish the signal from the background.The background distributions from W/Z+jets, t t, single top-quark and diboson processes are modelled using simulated events.The background shape from multijet production is obtained from an independent data sample that satisfies the signal selection criteria except for the lepton requirement: the electrons are required to satisfy a looser identification criterion ("medium" in Ref. [38]) but not meet the "tight" selection criteria; the selected muons are required to satisfy all the selection criteria after inverting the transverse impact parameter significance cut.The contribution from other processes is subtracted from data in the extraction of the multijet background shape.
The background contributions from t t, single top-quark and diboson production are normalized to the number of background events predicted by simulation.The p T (W) distribution in the W+jets simulated sample is corrected by comparing it to data in the LRR sidebands defined as 40 < m j j < 65 GeV or 105 < m j j < 200 GeV.The normalizations of the W/Z+jets and multijet background contributions are derived in a control data sample which is obtained by requiring the mass of the hadronic W/Z candidate to be within the m J (m j j ) sidebands.They are determined from binned minimum χ 2 fits to the E miss T distributions in the control data samples corresponding to each signal region and channel separately.The fitted parameters are the normalizations of these two processes.The difference of the W/Z+jets normalization from the expected background from simulation ranges between 1% and 18%.
The multijet background templates were validated in the electron channel using samples enriched in multijet events, obtained by inverting the E miss T requirement.The description of the t t background in simulation was validated in a sample dominated by top-pair events by requiring at least one b-tagged small-R jet.Good agreement within uncertainties is observed between data and expectation in these validation regions.

Systematic uncertainties
The main systematic uncertainty on the background estimation is the uncertainty on the normalization of W/Z+jets background obtained from the fit described above.This uncertainty is 3-4% in the LRR and HRR, and 13-19% in the MR.An uncertainty on the shape of the W/Z+jets background is obtained in the LRR by comparing data and simulation in the m j j sidebands, leading to an approximately 5% uncertainty for m ν j j < 600 GeV.Due to the low numbers of data events in the sidebands for the HRR and MR, the W+jets shape uncertainty in these regions is evaluated by comparing a sample of simulated events from Sherpa with a sample of simulated events from Alpgen [47] interfaced to Pythia6.The uncertainty in the shape of the t t mass distribution is estimated by comparing a sample from MC@NLO interfaced to Herwig with a sample from Powheg [48][49][50] interfaced to Pythia6.The uncertainty on the shape of the multijet background is evaluated by using alternative templates obtained by removing the calorimeterbased lepton isolation cuts.For the remaining background processes, detector-related uncertainties from the small-R jet energy scale and resolution, large-R jet energy and mass scale, lepton reconstruction and identification efficiencies, lepton momentum scales and resolutions, and missing transverse momentum were considered when evaluating possible systematic effects on the shape or normalization of the background estimation and are found to have a minor impact.The large-R jet energy and mass scale uncertainties are evaluated by comparing the ratio of calorimeter-based to track-based measurements in dijet data and simulation, and are validated by in-situ data of high-p T W production in association with jets.
The dominant uncertainty on the signal arises from initial-and final-state radiation modelling in Pythia and is less than 12% (6%) for G * (W ).Uncertainties due to the choice of PDFs are below 1%.

Results and interpretation
Table 1 shows the number of events predicted and observed in each signal region.The reconstructed m ν j j (m νJ ) distributions for data and predicted background events as well as selected benchmark signal models in the three signal regions are shown in Fig. 1 for the combined electron and muon channels.Good agreement is observed between the data and the background prediction.In the absence of a significant excess, the result is interpreted as 95% confidence level (CL) upper limits on the production cross section times branching ratio for the G * and W models.These upper limits are determined with the CL s modified frequentist formalism [51] with a profile-likelihood test statistic [52].The test statistic is evaluated with a maximum-likelihood fit of signal models and background predictions to the reconstructed m ν j j (m νJ )   spectra.Systematic uncertainties are taken into account as nuisance parameters with Gaussian sampling distributions.For each source of systematic uncertainty, the correlations across bins and between different kinematic regions, as well as those between signal and background, are taken into account.The likelihood fit is performed for signal pole masses between 300 and 800 GeV for the LRR, 600-1000 GeV for the HRR and 800-2000 GeV for the MR.Overlapping regions are fit simultaneously.Figure 2 shows 95% CL upper limits on the production cross section multiplied by the branching fraction into WW (WZ) for the bulk RS G * (EGM W ) as a function of the resonance pole mass.The theoretical predictions for the EGM W with a scale factor c = 1 and the bulk RS G * with coupling constant k/ MPl = 1, shown in the figure, allow observed lower mass limits of 1490 GeV for the W and 700 GeV for the G * to be extracted.The search is performed using an integrated luminosity of 20.3 fb −1 of pp collisions at √ s = 8 TeV collected by the ATLAS detector at the LHC.A set of event selections for bulk RS G * and EGM W boson signal is derived using simulated events and applied to the data.No evidence for resonant diboson production is observed and 95% CL upper limits on the production cross section times branching fraction of G * and W are determined.Resonance masses below 700 GeV are excluded for the spin-2 RS graviton G * and masses below 1490 GeV are excluded for the spin-1 EGM W boson at 95% CL.The analysis also sets the most stringent limits to date on the production cross section for W -like resonances decaying to WZ with masses around 2 TeV.The results represent a significant improvement over previously reported limits [11] in the same final state due to an increased data set size and the development of new techniques to analyse highly boosted bosons that decay hadronically.NCN, Poland; GRICES and FCT, Portugal; MNE/IFA, Romania; MES of Russia and ROSATOM, Russian Federation; JINR; MSTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SER, SNSF and Cantons of Bern and Geneva, Switzerland; NSC, Taiwan; TAEK, Turkey; STFC, the Royal Society and Leverhulme Trust, United Kingdom; DOE and NSF, United States of America.
[ p Also at Department of Physics, National Tsing Hua University, Taiwan q Also at Department of Physics, The University of Texas at Austin, Austin TX, United States of America r Also at Institute of Theoretical Physics, Ilia State University, Tbilisi, Georgia s Also at CERN, Geneva, Switzerland t Also at Georgian Technical University (GTU),Tbilisi, Georgia u Also at Ochadai Academic Production, Ochanomizu University, Tokyo, Japan v Also at Manhattan College, New York NY, United States of America w Also at Institute of Physics, Academia Sinica, Taipei, Taiwan x Also at LAL, Université Paris-Sud and CNRS/IN2P3, Orsay, France y Also at Academia Sinica Grid Computing, Institute of Physics, Academia Sinica, Taipei, Taiwan z Also at Dipartimento di Fisica, Sapienza Università di Roma, Roma, Italy aa Also at Moscow Institute of Physics and Technology State University, Dolgoprudny, Russia ab Also at Section de Physique, Université de Genève, Geneva, Switzerland ac Also at International School for Advanced Studies (SISSA), Trieste, Italy ad Also at Department of Physics and Astronomy, University of South Carolina, Columbia SC, United States of America ae Also at School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China a f Also at Faculty of Physics, M.V.

Figure 1 :
Figure1: Reconstructed mass distributions in data and the predicted backgrounds in the three kinematic regions referred to in the text as the low-p T resolved region (top left), high-p T resolved region (top right) and merged region (bottom).G * and W signal hypotheses of masses 400, 800 and 1200 GeV are also shown.The band denotes the statistical and systematic uncertainty on the background before the fit to the data.The lower panels show the ratio of data to the SM background estimate.

Figure 2 :
Figure 2: Observed and expected 95% CL upper limits on the cross section times branching fraction as a function of the resonance pole mass for the G * (top) and EGM W (bottom).The LO (NNLO) theoretical cross section for the G * (EGM W ) production is also shown.The inner and outer bands around the expected limits represent ±1σ and ±2σ variations respectively.The band around the W cross section corresponds to the NNLO theory uncertainty.

Table 1 :
Event yields in signal regions for data, predicted background contributions, and G * and W signals. Errors are shown before the fit to the data.The errors on the total background and total signal correspond to the full statistical and systematic uncertainty, while the errors on each background component include the full systematic uncertainty only.The signal points correspond to a resonance mass of 400 (LRR), 800 (HRR) and 1200 (MR) GeV.search for WW and WZ resonances decaying to a lepton, neutrino and jets is presented in this paper. A 9] CMS Collaboration, Search for massive resonances in dijet systems containing jets tagged as W or Z boson decays in pp collisions at √ s = 8 TeV, J. High Energy Phys. 08 (2014) 173, arXiv:1405.1994[hep-ex].Search for WZ resonances in the fully leptonic channel using pp collisions at √ s = 8 TeV with the ATLAS detector, Phys.Lett.B 737 (2014) 223, arXiv:1406.4456[hep-ex].TeV using the ATLAS detector at the LHC, Eur.Phys.J. C 73 (2013) 2518, arXiv:1302.4393[hep-ex].Also at Novosibirsk State University, Novosibirsk, Russia d Also at TRIUMF, Vancouver BC, Canada e Also at Department of Physics, California State University, Fresno CA, United States of America f Also at Department of Physics, University of Fribourg, Fribourg, Switzerland g Also at Departamento de Fisica e Astronomia, Faculdade de Ciencias, Universidade do Porto, Portugal h Also at Tomsk State University, Tomsk, Russia i Also at CPPM, Aix-Marseille Université and CNRS/IN2P3, Marseille, France j Also at Università di Napoli Parthenope, Napoli, Italy k Also at Institute of Particle Physics (IPP), Canada l Also at Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United Kingdom m Also at Department of Physics, St. Petersburg State Polytechnical University, St. Petersburg, Russia n Also at Louisiana Tech University, Ruston LA, United States of America o Also at Institucio Catalana de Recerca i Estudis Avancats, ICREA, Barcelona, Spain c Lomonosov Moscow State University, Moscow, Russia ag Also at National Research Nuclear University MEPhI, Moscow, Russia ah Also at Department of Physics, Stanford University, Stanford CA, United States of America ai Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest, Hungary a j Also at Department of Physics, Oxford University, Oxford, United Kingdom ak Also at Department of Physics, The University of Michigan, Ann Arbor MI, United States of America al Also at Discipline of Physics, University of KwaZulu-Natal, Durban, South Africa am Also at University of Malaya, Department of Physics, Kuala Lumpur, Malaysia * Deceased