Search for triboson $W^{\pm}W^{\pm}W^{\mp}$ production in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector

This paper reports a search for triboson $W^{\pm}W^{\pm}W^{\mp}$ production in two decay channels ($W^{\pm}W^{\pm}W^{\mp}\rightarrow \ell^{\pm}\nu\ell^{\pm}\nu\ell^{\mp}\nu$ and $W^{\pm}W^{\pm}W^{\mp}\rightarrow \ell^{\pm}\nu\ell^{\pm}\nu{}jj$ with $\ell=e, \mu$) in proton-proton collision data corresponding to an integrated luminosity of 20.3 fb$^{-1}$ at a centre-of-mass energy of 8 TeV with the ATLAS detector at the Large Hadron Collider. Events with exactly three charged leptons, or two leptons with the same electric charge in association with two jets, are selected. The total number of events observed in data is consistent with the Standard Model (SM) predictions. The observed 95 % confidence level upper limit on the SM $W^{\pm}W^{\pm}W^{\mp}$ production cross section is found to be 730 fb with an expected limit of 560 fb in the absence of SM $W^{\pm}W^{\pm}W^{\mp}$ production. Limits are also set on $WWWW$ anomalous quartic gauge couplings.


Introduction
The triple gauge couplings (TGCs) and quartic gauge couplings (QGCs) that describe the strengths of the triple and quartic gauge boson self-interactions are completely determined by the non-Abelian nature of the electroweak SU(2) L × U(1) Y gauge structure in the Standard Model (SM). These interactions contribute directly to diboson and triboson production at colliders. Studies of triboson production can test these interactions and any possible observed deviation from the theoretical prediction would provide hints of new physics at a higher energy scale. Compared with TGCs, QGCs are usually harder to study due to the, in general, smaller production cross sections of the relevant processes.
In the SM, charged QGC interactions (W W W W , W W Z Z, W W Zγ and W W γ γ ) are allowed whereas neutral QGC interactions (Z Z Z Z, Z Z Zγ , Z Zγ γ , Z γ γ γ and γ γ γ γ ) are forbidden. Searches  ] used W ± W ± vector-boson scattering events, whereas this paper presents the first search for W W W W QGC interactions via triboson W ± W ± W ∓ production and sets the first limit on the total SM W ± W ± W ∓ production crosssection using proton-proton ( pp) collision data collected with the ATLAS detector and corresponding to an integrated luminosity of 20.3 fb −1 [18] at a centre-of-mass energy of 8 TeV. Two decay channels, W ± W ± W ∓ → ± ν ± ν ∓ ν and W ± W ± W ∓ → ± ν ± ν j j, with = e or μ, are considered and are hereafter referred to simply as ν ν ν and ν ν j j channels, respectively.

The ATLAS detector
The ATLAS detector [19] is composed of an inner tracking detector (ID) surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, electromagnetic and hadronic calorimeters, and a muon spectrometer (MS). The ID consists of three subsystems: the pixel and silicon microstrip detectors that cover |η| < 2.5 in pseudorapidity, 1 and the outer transition radiation tracker that has an acceptance range of |η| < 2.0. The finely-segmented electromagnetic calorimeter is composed of lead absorbers with liquid argon (LAr) as the active material, spanning |η| < 3.2. In the region |η| < 1.8, a pre-sampler detector using a thin layer of LAr is used to correct for the energy loss by electrons and photons upstream of the calorimeter. The hadronic 1 The ATLAS experiment uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector. The x-axis points from the IP to the centre of the LHC ring, the y-axis points upward, and the z-axis is along one of the proton beam directions. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). Transverse momentum ( p T ) is defined relative to the beam axis and is calculated as p T = p sin θ where p is the momentum. tile calorimeter (|η| < 1.7) consists of steel absorbers and scintillating tiles and is located directly outside the envelope of the barrel electromagnetic calorimeter. The endcap hadronic calorimeters use LAr as active material, with copper as absorber material, while the forward calorimeters use LAr as active material, with copper absorber for the first layer, dedicated to electromagnetic measurements, and tungsten for other layers, dedicated to hadronic measurements. The MS is composed of three large superconducting air-core toroidal magnets, a system of three stations of tracking chambers in the range |η| < 2.7, and a muon trigger system in the range |η| < 2.4. The precision muon momentum measurement is performed by monitored drift tubes everywhere except in the innermost layer for the range |η| > 2.0 where cathode strip chambers are used instead. The muon trigger system is composed of resistive plate chambers in the barrel region (|η| < 1.05) and thin gap chambers in the endcap region (1.05 < |η| < 2.4).
The ATLAS trigger system has three distinct levels referred to as L1, L2, and the event filter. Each trigger level refines the decisions made at the previous level. The L1 trigger is implemented in hardware and uses a subset of detector information to reduce the event rate to a design value of at most 75 kHz. The L2 and event filter are software-based trigger levels and together reduce the event rate to about 400 Hz.
Events used were selected by single-lepton triggers with a transverse momentum, p T , threshold of 24 GeV for both muons and electrons, along with an isolation requirement. The single-lepton triggers are complemented with triggers having a higher p T threshold (60 GeV for electrons and 36 GeV for muons) and no isolation requirement in order to increase the acceptance at high p T .

Object reconstruction and event selection
Each event is required to have at least one primary vertex reconstructed from at least three tracks with p T > 400 MeV. If there are multiple primary vertices reconstructed in the event due to additional pp interactions (pile-up) in the same or a neighbouring bunch crossing, the vertex with the highest p 2 T , calculated using all associated tracks, is taken as the primary collision vertex. The mean number of interactions per bunch crossing in this data set is 20.7.
Electron candidates [20] are required to have p T > 20 GeV and |η| < 2.47. Candidates within the transition region between the barrel and endcap calorimeters (1.37 < |η| < 1.52) are rejected. In addition, they must satisfy the tight quality definition described in Ref. [21]. Muon candidates are reconstructed by combining tracks in the ID with tracks in the MS and have p T > 20 GeV and |η| < 2.5. The ID tracks associated with these muons must pass a number of quality requirements [22].
To ensure that lepton candidates originate from the primary vertex, a requirement is placed on the longitudinal impact parameter, z 0 , multiplied by the sine of the track polar angle, θ , such that the absolute value is smaller than 0.5 mm (|z 0 × sin θ | < 0.5 mm). A requirement is also placed on the transverse impact parameter, d 0 , divided by its resolution (σ d 0 ), such that |d 0 /σ d 0 | < 3. To suppress the contribution from hadronic jets which are misidentified as leptons, signal leptons are required to be isolated in both the ID and the calorimeter. The calorimeter isolation is defined as E ConeX T /E T whereas the ID isolation is defined as p ConeX is the transverse energy (momentum) deposited in the calorimeter (the scalar sum of the p T of tracks with p T > 1 GeV) within a cone of size R = ( η) 2 + ( φ) 2 = X around the lepton. The transverse momentum from the lepton itself is excluded in the calculations of E ConeX T and p ConeX T . Different lepton isolation criteria are applied in the two channels to maximize the signal efficiency while suppressing the backgrounds. In the ν ν ν channel, Jets are reconstructed from clusters of energy in the calorimeter using the anti-k t algorithm [23] with radius parameter R = 0.4. Jet energies are calibrated using energyand η-dependent correction factors derived using Monte Carlo (MC) simulation and validated by studies of collision data [24]. For jets with p T < 50 GeV and |η| < 2.4, at least 50% of the summed scalar p T of the tracks within a cone of size R = 0.4 around the jet axis must originate from the primary vertex. This requirement reduces the number of jet candidates originating from pile-up vertices. Jets containing b-hadrons ("b-jets") with |η| < 2.5 and p T > 25 GeV are identified using the impact parameter significance of tracks in the jet and secondary vertices reconstructed from these tracks [25,26]. In the ν ν ν and ν ν j j channels, the efficiency of the b-tagging algorithm used is 85 and 70%, respectively.
The measurement of the two-dimensional missing transverse momentum vector, p miss At most one jet with p T > 25 GeV and |η| < 4.5 b-jet veto No identified b−jets with p T > 25 GeV and |η| < 2.5 The experimental signature of the ν ν ν channel is the presence of three charged leptons and E miss T . The signature of the ν ν j j channel is the presence of two same-charge leptons, E miss T , and two jets with an invariant mass close to 80 GeV. The selection requirements used to define the signal regions described in the following are obtained from a multidimensional optimization to maximize the sensitivity to the W ± W ± W ∓ process and to reduce the contributions from SM background processes.
To select ν ν ν candidates, events are required to have exactly three charged leptons with p T > 20 GeV, at most one jet with p T > 25 GeV and |η| < 4.5, and no identified b-jets. In addition, the absolute value of the azimuthal angle between the trilepton system and the p miss T , |φ 3 − φ p miss T |, is required to be above 2.5. Eight different final states with equal production probability are considered based on the flavour and the charge of the leptons, namely e ± e ± e ∓ , e ± e ∓ μ ± , e ± e ∓ μ ∓ , e ± e ± μ ∓ , μ ± μ ∓ e ± , μ ± μ ∓ e ∓ , μ ± μ ± e ∓ , and μ ± μ ± μ ∓ . Three separate signal regions are defined based on the number of same-flavour opposite-sign (SFOS) lepton pairs in the event: 0 SFOS (e ± e ± μ ∓ and μ ± μ ± e ∓ ), 1 SFOS (e ± e ∓ μ ± , e ± e ∓ μ ∓ , μ ± μ ∓ e ± , and μ ± μ ∓ e ∓ ), and 2 SFOS (e ± e ± e ∓ and μ ± μ ± μ ∓ ). In the 0-SFOS case, the invariant mass of the same-flavour lepton pair, m , is required to be greater than 20 GeV. If there are at least two electrons in the event, the di-electron invariant mass, m ee , is required to have |m ee − m Z | > 15 GeV, where m Z is the pole mass of the Z boson [28]. No requirement is applied on the E miss T variable, as it was found to not discriminate between signal and backgrounds. In the 1-SFOS case, the SFOS dilepton invariant mass, m SFOS , is required to be outside of the region m Z − 35 GeV < m SFOS < m Z + 20 GeV. In addition, events are required to satisfy E miss T > 45 GeV. Finally, in the 2-SFOS case, the SFOS dilepton invariant masses are required to have |m SFOS − m Z | > 20 GeV while the E miss T must be greater than 55 GeV. The selection criteria for m SFOS and E miss T are mainly used to reduce the contributions from the Z +jets and W Z+jets processes. Table 1 shows the kinematic selection criteria used for the ν ν ν channel.
To select ν ν j j candidates, events are required to have exactly two leptons with the same electric charge, at least two jets, and no identified b-jets. Three different final states are considered based on the lepton flavour, namely e ± e ± , e ± μ ± , and μ ± μ ± . The lepton p T (E miss T ) threshold is set to 30 (55) GeV to reduce the SM background contributions, though the E miss T criterion is not applied for the μ ± μ ± final state due to the smaller Z +jets background expected in this channel. The leading (sub-leading) p T jet must have p T > 30 (20) GeV and |η| < 2.5. The two jets are required to have 65 GeV < m j j < 105 GeV and | η j j | < 1.5 in order to distinguish the signal from the W ± W ± backgrounds, where m j j is the dijet invariant mass and η j j is the pseudorapidity separation between the two jets. The dilepton system is required to have m > 40 GeV and in the case of the e ± e ± final state, m ee must have m ee < 80 GeV or m ee > 100 GeV in order to suppress events with two opposite-sign prompt leptons where the charge of one of the electrons is misidentified. To reduce the contributions from W Z+jets and Z Z+jets production, events are removed if they contain additional leptons reconstructed with p T > 6 GeV passing looser identification quality requirements, with a medium identification requirement for electrons as defined in Ref.
[21] and the minimum identification required for muon reconstruction. Table 2 shows the kinematic selection criteria used for the ν ν j j channel. Table 2 Selection criteria for the ν ν j j channel, split based on the lepton flavour: e ± e ± , e ± μ ± , and μ ± μ ± ν ν j j e ± e ± e ± μ ± μ ± μ ±  [30]. However, they are not considered in this analysis. In order to determine W ± W ± W ∓ production cross sections, events are generated at NLO in QCD using Mad-Graph5_aMC@NLO [31] including on-shell diagrams as well as Higgs associated diagrams. The CT10 NLO parton distribution function (PDF) [32] is used. Subsequent decays of unstable particles and parton showers are handled by pythia8 [33]. Fiducial cross sections are calculated using the generator-level lepton, jet, and E miss T definitions as described in Ref. [34]. Generator-level prompt leptons (those not originating from hadron and τ lepton decays) are dressed with prompt photons within a cone of size R = 0.1. Generatorlevel jets are reconstructed by applying the anti-k t algorithm with radius parameter R = 0.4 on all final-state particles after parton showering and hadronisation. The E miss T variable is calculated using all generator-level neutrinos. The same kinematic selection criteria as listed in Tables 1 and 2 are applied on these objects, with the exception of the b-jet veto requirements in the ν ν ν channel and the lepton quality requirements. To take into account the effect of the lepton isolation in the fiducial region, any lepton pairs must satisfy R( , ) > 0.1, and in the ν ν j j channel any lepton-jet pairs must satisfy R( j, ) > 0.3. Electrons or muons from τ decays are not included.
The fiducial cross section is predicted to be 309 ± 7 (stat.) ± 15 (PDF) ± 8 (scale) ab in the ν ν ν channel and 286 ± 6 (stat.) ± 15 (PDF) ± 10 (scale) ab in the ν ν j j channel. Uncertainties due to the PDFs are computed using an envelope of the CT10, NNPDF3.0 [35], and MSTW2008 [36] NLO PDF 68 or 90% (for CT10) confidence level (CL) uncertainties, following the recommendation of Ref. [37]. The renormalization and factorization scales are set to the invariant mass of the W W W system. Scale uncertainties are estimated by varying the two scales independently up and down by a factor of two and taking the largest variation from the nominal cross-section values.
In order to combine the measurements from the two decay channels, a common phase space is defined where each W boson can decay either leptonically (including τ leptons) or hadronically, pp → W ± W ± W ∓ + X , with no kinematic requirements placed on the final-state leptons but with jets restricted to have p T > 10 GeV. The extrapolation factor from the fiducial phase space to the total phase space is large, but it is mainly due to the well-known W boson decay branching ratios. The total cross section in this common phase space is 241.5 ± 0.1 (stat.) ±10.3 (PDF) ±6.3 (scale) fb.
In order to determine the detector reconstruction effects on the signal selection, W ± W ± W ∓ signal samples are generated with vbfnlo [29, [38][39][40] at LO. The parton shower and hadronisation are performed by pythia8. The fiducial cross sections are seen to be consistent between vbfnlo and MadGraph5_aMC@NLO when computed at the same order. The vbfnlo LO fiducial cross sections are normalized to the NLO fiducial cross section predicted by Mad-Graph5_aMC@NLO for the signal yield calculations. These events are processed through the full ATLAS detector simulation [41] based on Geant 4 [42]. To simulate the effect of multiple pp interactions occurring during the same or a neighbouring bunch crossing, minimum-bias interactions are generated and overlaid on the hard-scattering process. These events are then processed through the same object reconstruction and identification algorithms as used on data. MC events are reweighted so that the pile-up conditions in the simulation match the data. Additional corrections are made to the simulated samples to account for small differences between the simulation and the data for the object identification and reconstruction efficiencies, the trigger efficiencies, and the energy and momentum scales and resolutions. While excluded in the fiducial cross-section definition, the contribution from events with W → τ ν → ννν decays are counted as signal in the vbfnlo signal sample used in the final event selection. These events contribute up to 20% of the predicted signal yield. This approach is used to ease comparisons of the obtained cross-section limits with alternative cross-section predictions that may not simulate tau decays.

Background estimation
The SM processes that mimic the W ± W ± W ∓ signal signature can be grouped into five categories: • The W Z/γ * +jets process that produces three prompt leptons or two prompt leptons with the same electric charge (referred to as "W Z background"); • The W γ +jets or Z γ +jets processes where the photon is misreconstructed as a lepton (referred to as "V γ background", where V = W, Z ); • Processes other than W Z/γ * +jets that produce three prompt leptons or two prompt leptons with the same electric charge (referred to as "other prompt background"); • Processes that produce two or three prompt charged leptons, but the charge of one lepton is misidentified (referred to as "charge-flip background"); • Processes that have one or two non-prompt leptons originating either from misidentified jets or from hadronic decays (referred to as "fake-lepton background").
The dominant irreducible background originates from the W Z(→ ± ν ± ∓ )+jets process and is estimated using simulated events. In the ν ν ν channel these events are generated with Powheg-BOX [43][44][45][46] and hadronised with pythia8 and in the ν ν j j channel they are generated with Sherpa [47]. In the ν ν ν channel, the inclusive W Z+jets cross section is normalized using a scale factor (1.08 ± 0.10) derived from a W Z-enriched region in data. This region is obtained by requiring exactly one SFOS lepton pair with |m SFOS − m Z | < 15 GeV. In the ν ν j j channel, the cross section is normalized to the NLO calculation in QCD from vbfnlo [48] in the specified fiducial phase space with a normalization factor of 1.04 ± 0.09.
The V γ background contributes when the photon is misidentified as an electron. In the ν ν ν channel, this originates primarily from the Z γ process and its contribution is estimated using events generated with Sherpa. In the ν ν j j channel, this comes primarily from electroweak and strong production of W γ j j events. Strong production of W γ j j [49] is estimated using Alpgen [50] interfaced to Herwig [51] and Jimmy [52] for simulation of the parton shower, fragmentation, hadronisation and the underlying event. The electroweak production of W γ j j [53] is modelled using Sherpa.
Other SM processes that produce multiple prompt leptons include Z Z, tt V , Z W W , Z Z Z, W ± W ± j j production, and double parton scattering processes. The production of Z Z is modelled with Powheg-BOX [46] and hadronised with pythia8 in the ν ν ν channel and is modelled with Sherpa in the ν ν j j channel. The tt V [54], Z W W [55], and Z Z Z [55] processes are modelled using MadGraph5_aMC@NLO together with pythia8 for both channels. The non-resonant W ± W ± j j background [56] is only important for the ν ν j j channel and its contribution is estimated using Sherpa. Contributions from double parton scattering processes are found to be negligible in both channels.
The charge-flip background originates from processes where the charge of at least one prompt lepton is misidentified. This occurs primarily when a lepton from a hard bremsstrahlung photon conversion is recorded instead of the signal lepton. It mainly contributes to the 0-SFOS signal region in the ν ν ν channel and the e ± e ± / e ± μ ± signal regions in the ν ν j j channel. The electron charge misidentification rate is measured using Z → e + e − events. In the ν ν ν channel, the charge-flip background is estimated by using these rates to re-weight the MC estimate of W Z and Z Z events based on the probability for opposite-sign events of this kind to migrate into the 0-SFOS category. In the ν ν j j channel, the background is estimated by applying these rates on data events satisfying all signal selection criteria except the two leptons are required to have opposite-sign.
Contributions from fake-lepton backgrounds are estimated in data, using different approaches in the two chan- nels. In the ν ν ν channel, the probabilities of prompt leptons or non-prompt leptons to satisfy the signal lepton criteria are computed using a tag-and-probe method whereby a well-reconstructed "tag" lepton is used to identify the event and a second "probe" lepton is used to study the probabilities without bias. A tag lepton must satisfy the signal lepton requirements while a looser lepton selection criterion is defined for probe leptons with the lepton isolation requirements removed and the electron quality requirement loosened to medium as defined in Ref.
[21]. The probability for a prompt lepton to satisfy the signal lepton criteria is estimated using Z → + − events with the tag-and-probe lepton pair required to have the same-flavour, opposite-sign and an invariant mass within 10 GeV of the pole mass of the Z boson.
The probability for a non-prompt lepton from hadronic activity to satisfy the signal lepton requirement is estimated using the tag-and-probe method in a W +jets-enriched region with E miss T > 10 GeV, the tag lepton is a muon with p T > 40 GeV, and the tag and probe leptons have the same electric charge. The probabilities are calculated separately for electrons and muons. A loosely identified set of data is also selected by requiring at least three loose leptons as defined above. This set of data, along with these probabilities are then used to estimate the background in the signal region with the matrix method [57].
In the ν ν j j channel, events that contain one signal lepton and one "lepton-like" jet are selected. A "lepton-like" jet satisfies all signal lepton selection criteria except that the isolation requirements are 0.14 < E Cone0.3 T / p T < 2 and 0.06 < p Cone0.3 T / p T < 2 for electrons, and 0.07 < E Cone0.3 T / p T < 2 and 0.07 < p Cone0.3 T / p T < 2 for muons. In addition, the |d 0 /σ d 0 | and |z 0 × sin θ | selection criteria are loosened to 10 mm and 5 mm, respectively. These events are dominated by non-prompt leptons and are scaled by a fake factor to estimate the non-prompt background. The fake factor is the ratio of the number of jets satisfying the signal lepton identification criteria to the number of jets satisfying the "lepton-like" jet criteria. It is measured as a function of the jet p T and η from a dijet-enriched sample selected by requiring a lepton back-to-back with a jet ( φ j > 2.8) and E miss T < 40 GeV.

Validation of background estimates
The background predictions are tested in several validation regions (VRs). These VRs are defined to be close to the signal region with a few selection criteria removed or inverted. They generally have dominant contributions from one or two background sources and a negligible contribution from the signal process. The signal and background predictions are compared to data for each VR in Table 3.
In the ν ν ν channel, three VRs are considered. The first VR, called the pre-selection region, tests the modelling of the W Z+jets background by requiring exactly three signal leptons. The distribution of the trilepton transverse mass, m 3 where p 3 T is the p T of the trilepton system, is shown at the top left of Fig. 2. This VR includes the three signal regions (0, 1, and 2 SFOS), but the effect of the signal is considered negligible at this stage of the selection, as shown in Table 3. The W Z+jets purity is estimated to be around 70% in this region. The second region, called the fake-lepton region, tests the modelling of the fakelepton background by requiring exactly three signal leptons with no SFOS lepton pairs and at least one b−jet. The distribution of the jet multiplicity, N jet , is shown at the top right of Fig. 2. The purity of the fake-lepton background is estimated to be around 80% in this region. The third region, called Z γ region, tests the modelling of the Z γ background by requiring the presence of only μ + μ − e ± events where the trilepton invariant mass is close to the Z resonance peak. This restricts the main contributions to originate from the Z γ → μ + μ − γ and Z → μ + μ − → μ + μ − γ processes. The Z γ purity is estimated to be around 70% in this region. The data are seen to be well described by the background in all three VRs.
In the ν ν j j channel, five VRs are considered. The modelling of the charge-flip background is tested using e ± e ± events with 80 GeV < m < 100 GeV. The purity of the  Fig. 2. The purity of the W Z+jets is estimated to be around 60% in this region. The modelling of backgrounds from non-prompt leptons is tested in a b−tagged region that requires at least one b−jet. The purity of the nonprompt lepton background is estimated to be around 80% in this region. The m j j modelling is checked by examining events with masses m j j in the regions m j j < 65 GeV or m j j > 105 GeV. The distribution of m j j in this region is shown at the bottom right of Fig. 2. Finally, conversion and prompt backgrounds are tested in a region with at most one jet, called the ≤ 1 jet region. The purity of the conversion and prompt backgrounds is estimated to be around 70% in this region. As for the ν ν ν channel, good agreement is observed between the data and the prediction in all five VRs.

Systematic uncertainties
Systematic uncertainties in the signal and background predictions arise from the measurement of the integrated lumi- nosity, from the experimental and theoretical modelling of the signal acceptance and detection efficiency, and from the background estimation. The effect of the systematic uncertainties on the overall signal and background yields are evaluated separately for the ν ν ν and ν ν j j channels. The results are summarised in Table 4. The systematic uncertainties are included as nuisance parameters in the profile likelihood described in Sect. 7. Correlations of systematic uncertainties arising from common sources are maintained across signal and background processes and channels. The experimental uncertainties include the uncertainties on the lepton and jet energy and momentum scales and resolutions, on the efficiencies of the lepton and jet reconstruction and identification, and on the modelling of E miss T and b-jets. They are evaluated separately for both the signal and background estimations. For the expected signal yield, the major contributions in the ν ν ν channel come from uncertainties in the lepton reconstruction and identification efficiencies as well as lepton energy/momentum resolution and scale modelling (±1.6%), E miss T modelling (±1.1%), and jet energy scale and resolution (±2.3%). The contributions in the ν ν j j channel come from uncertainties in the lepton efficiencies and energy/momentum modelling (±2.1%), E miss T modelling (±0.7%), b-jet identification (±2.2%), and jet energy resolution and scale modelling (±21%). Larger systematic uncertainties due to the jet energy scale and resolution are expected in the ν ν j j channel due to the dijet requirements, in particular in the dijet invariant mass. For the background yields estimated from MC simulation, the major contributions in the ν ν ν channel come from uncertainties in lepton reconstruction and identification efficiencies (±1.8%), E miss T modelling (±1.4%), and jet energy modelling (±2.8%). The major contributions in the ν ν j j channel come from uncertainties in lepton efficiencies and energy modelling (±3.3%), E miss T modelling (±1.8%), b-jet identification (±2.2%), and jet energy modelling (±15%).
The estimates of the data-driven fake-lepton background also have uncertainties specific to each channel. In the ν ν ν channel, the systematic uncertainty results from the uncertainties on the probabilities of candidate leptons that satisfy the looser lepton selection criteria to also satisfy the signal lepton selection criteria. For prompt leptons this uncertainty is ±(5 to 10)% while for fake leptons and misidentified/non-prompt leptons this uncertainty is ±(80 to 90)%. The latter uncertainty is a conservative estimate which accounts for differences in the heavy-flavour and lightflavour composition between the signal region and the control region where the fake-lepton efficiency is determined for these leptons. In the ν ν j j channel, the systematic uncertainty results from the uncertainties in the measurement of the fake factors, which is estimated to be ±(20 to 30)%. Statistical uncertainties in the samples used for the matrix method and the fake-factor method also contribute to the overall uncertainty of the estimation of the fake-lepton background. The total uncertainty in the overall fake-lepton background yield is ±13% in the ν ν ν channel and ±8% in the ν ν j j channel.
The charge-flip background is only relevant for the e ± e ± /e ± μ ± final state in the ν ν j j channel and for the 0-SFOS region in the ν ν ν channel. Its uncertainty is dominated by the statistical precision with which the electron charge misidentification rate is determined from the available data. Since the charge-flip background estimation uses the number of Z (→ e + e − ) + 2 jets events, the number of events in the data also contributes to the overall systematic uncertainty. In the ν ν ν channel, the uncertainty on the charge-flip background estimate is ±0.5% in the 0-SFOS region but is ±0.04% for the total background estimate in all three signal regions. In the ν ν j j channel, the total system- atic uncertainty in the overall background yield due to the uncertainty in the charge-flip background estimate is found to be ±2.2%.
There are also uncertainties in the overall normalization of the signal and MC background cross sections. Uncertainties in the signal cross section are those described in Sect. 4. These are not, however, included as uncertainties in the model and merely serve as a comparison for the final measurement in Sect. 7. The normalizations of the SM background cross sections described in Sect. 5.1 have their own associated uncertainties. The uncertainty in the predicted W Z+jets background cross section is the most important one since it is the largest irreducible background. The size of the uncertainty relative to the predicted W Z+jets background is ±10% in the ν ν ν channel and ±(16 to 23)% in the ν ν j j channel depending on the production mechanism. This uncertainty is based on the measurement performed in the control region, for the ν ν ν channel as described in Sect. 5, while it is a combination of the scale, PDF and parton shower uncertainties estimated as in Ref.
[12], for the ν ν j j channel. The remaining uncertainties are mostly negligible in the overall background prediction. The normalization uncertainty in the total background prediction is around ±8% in the ν ν ν channel and ±13% in the ν ν j j channel.
The uncertainty on the integrated luminosity is ±1.9%, affecting the overall normalization of both the signal and background processes estimated from MC simulation. It is derived following the methodology detailed in Ref. [18]. The uncertainties associated with the pile-up reweighting of the events are estimated to be no more than ±0.1% for the signal and the backgrounds.

Cross-section measurement
The signal and background predictions together with their uncertainties are compared to the data for six signal regions in Table 5. The expected signal yields are calculated using the SM W ± W ± W ∓ cross sections listed in Sect. 4. The expected numbers of signal plus background events are consistent with the numbers of events observed in data in all regions. Figure 3 shows the m 3 T distribution for the ν ν ν channel and the distribution of the sum of the scalar p T for all selected objects, T + E miss T , for the ν ν j j channel, after summing over the three signal regions in each channel. Good agreement between data and the signal-plusbackground model is observed for both distributions.
The amount of W ± W ± W ∓ signal in the selected data set is determined using the numbers of expected signal and background events as well the numbers of observed events in the data. The signal strength, μ, is the parameter of interest, defined as a scale factor multiplying the cross section times branching ratio predicted by the SM. A test statistic based on the profile-likelihood ratio [58] is used to extract μ from a maximum-likelihood fit of the signal-plusbackground model to the data. The likelihood, L, is given by where the index c represents one of the two analysis channels, i represents one of the three signal regions in each channel, n obs is number of observed events, n sig, SM is the expected number of signal events based on the SM calculations, and n bkg is the expected number of background events. The effect of a systematic uncertainty k on the likelihood is modelled with a nuisance parameter, θ k , constrained with a corresponding Gaussian probability density function g(θ k ). The test statistic, t μ , is defined as whereμ is the unconditional maximum-likelihood (ML) estimators of the independent signal strength in the categories being compared,θ are the unconditional ML estimators of the nuisance parameters, andθ(μ) are the conditional ML estimators of θ for a given value of μ. The significance of μ is obtained with the above test statistic, and is estimated using 100,000 MC pseudo-experiments to determine how well the fit result agrees with the background-only hypothesis. The observed (expected) significance of a positive signal cross section is 0.96 σ (1.05 σ ) for the combination of the two channels. Most of the sensitivity comes from the 0-SFOS category in the ν ν ν channel and the μ ± μ ± category in the ν ν j j channel. The most significant deviation from the signal-plus-background hypothesis occurs in the e ± e ± region, where zero events are observed and 4.0 background and 0.46 signal events are expected. The probability that the background fluctuates down to zero events is 2.3%. The central value of μ corresponds to the minimum of the negative log-likelihood distribution. The measured fiducial cross section in each channel is obtained by multiplying μ by the expected value of the fiducial cross section in that channel. The measured total cross section is obtained by combining the results for the two channels and then extrapolating to the total phase space using the signal acceptance. The log-likelihood scans for the total crosssection measurement are evaluated with and without systematic uncertainties and are shown in Fig. 4. The expected and observed fiducial and total cross sections are summarized in Table 6. The presence of the W ± W ± W ∓ signal is also assessed using a one-sided 95% CL upper limit on the production cross section using the CL s method of Ref. [59]. The limits are evaluated using 2000 MC pseudo-experiments. The observed (expected) upper limit on the fiducial cross section in the absence of W ± W ± W ∓ production is found to be 1.3 fb (1.1 fb) in the ν ν ν channel and 1.1 fb (0.9 fb) in the ν ν j j channel. The observed (expected) upper limit in the absence of W ± W ± W ∓ production on the total cross section is 730 fb (560 fb) when the two channels are combined. If the SM W ± W ± W ∓ signal is also considered, the expected upper limit on the total cross section is 850 fb.

Limits on anomalous quartic gauge couplings (aQGCs)
Contributions from sources beyond the SM to the W ± W ± W ∓ production process can be expressed in a model-independent way using higher-dimensional operators leading to W W W W aQGCs. The parameterization of aQGCs is based on Ref. [60] in a linear representation [61] considering only dimensioneight operators involving four gauge bosons. There are 18 dimension-eight operators built from the covariant derivative of the Higgs field D μ , the SU(2) L field strength W i μν , and U(1) Y field strength B μν . Only the two terms built exclusively from D μ and with aQGC parameters f S,0 / 4 and f S,1 / 4 are considered in this analysis: where is the energy scale of the new physics. These two operators only affect massive bosons and do not depend on the gauge boson momenta since no SU(2) L or U(1) Y field strengths are included. As a result, they are important for the study of longitudinal vector-boson scattering. Similar parameters were studied before by the ATLAS and CMS Collaborations in Refs. [8, 10,16]. The effective Lagrangian approach leads to tree-level unitarity violation. This can be avoided by introducing a form factor [62] as where α corresponds to one of the two couplings, α 0 is the value of the aQGC at low energy,ŝ is the square of the partonic centre-of-mass energy, and FF is the form-factor cutoff scale. However, there is no theoretical algorithm to predict for which form-factor cutoff scale the cross section would violate unitarity. Therefore different values of FF are considered with FF = 0.5, 1, 2, and 3 TeV as well as FF = ∞, which corresponds to the non-unitarized case. Events with aQGCs are generated with vbfnlo at LO and passed through the ATLAS detector simulation. A grid of samples is obtained using different parameters of f S,0 / 4 and f S,1 / 4 values. The interpolation between these points is performed with a 2-dimensional quadratic function in the ( f S,0 / 4 , f S,1 / 4 ) space. The LO samples are scaled using Table 6 The predicted and observed fiducial cross sections for the ν ν ν and ν ν j j channels and the predicted and observed total cross section for the combination of the two channels   T distribution in the ν ν ν channel and the p T distribution in the ν ν j j channel, summed over the three signal regions in each channel. Even though aQGC events tend to have leptons or jets with larger momenta, the detection efficiency for events in the fiducial region is found to be consistent with the one obtained for the SM sample within 20%. The efficiencies of the aQGC samples are used with their statistical and systematic uncertainties to derive the 95% confidence intervals (CI) on aQGC, while the largest observed deviation of the aQGC efficiencies from the SM one is used as an extra systematic uncertainty. Frequentist CI on the anomalous coupling are computed by forming a profile-likelihoodratio test that incorporates the observed and expected numbers of signal events for different values of the anomalous couplings. Table 7 shows the expected and observed 95% CI on f S,0 / 4 ( f S,1 / 4 ) with different FF values, assuming f S,1 / 4 ( f S,0 / 4 ) to be zero. Figure 5 shows the twodimensional 95% CL contour limits of f S,0 / 4 vs f S,1 / 4 in the cases where FF = 1 TeV and FF = ∞. For FF = ∞, the limits can be compared to the stronger limits obtained by the CMS Collaboration in Ref.
[16] in a different production channel. Other parameterization (α 4 , α 5 ) of new physics have been introduced in Refs. [63][64][65]. The limits presented in this paper can be converted into limits on α 4 and α 5 following the formalism defined in the Appendix of Ref. [60] and using Equations (60) and (61) 10]. The latter were obtained using a different unitarization scheme. Since that scheme is not applicable to triboson production, a combination of the limits is not possible.

Summary
A search for triboson W ± W ± W ∓ production in two decay channels (W ± W ± W ∓ → ± ν ± ν ∓ ν and W ± W ± W ∓ → ± ν ± ν j j with = e, μ) is reported, using proton-proton collision data corresponding to an integrated luminosity of 20.3 fb −1 at a centre-of-mass energy of 8 TeV collected by the ATLAS detector at the LHC. Events with exactly three charged leptons or two same-charge leptons in association with two jets are selected. The data are found to be in good agreement with the SM predictions in all signal regions. The observed 95% CL upper limit on the SM W ± W ± W ∓ production cross section is found to be 730 fb with an expected limit of 560 fb in the absence of W ± W ± W ∓ production. Limits are also set on the aQGC parameters f S,0 / 4 and f S,1 / 4 .