Measurement of $WW/WZ \to \ell \nu q q^{\prime}$ production with the hadronically decaying boson reconstructed as one or two jets in $pp$ collisions at $\sqrt{s}=8$ TeV with ATLAS, and constraints on anomalous gauge couplings
ATLAS Collaboration

TL;DR
This study measures $WW/WZ$ production cross-sections in proton-proton collisions at 8 TeV using ATLAS data, reconstructs hadronic decays as jets, and sets limits on anomalous gauge couplings, confirming Standard Model predictions.
Contribution
It provides the first measurements of $WW/WZ$ production with hadronic decays reconstructed as one or two jets at 8 TeV and constrains anomalous gauge couplings.
Findings
Measured cross-sections consistent with Standard Model.
Reconstructed hadronic bosons as single or double jets.
Set 95% confidence limits on anomalous couplings.
Abstract
This paper presents a study of the production of or boson pairs, with one boson decaying to or and one or boson decaying hadronically. The analysis uses 20.2 fb of TeV collision data, collected by the ATLAS detector at the Large Hadron Collider. Cross-sections for production are measured in high- fiducial regions defined close to the experimental event selection. The cross-section is measured for the case where the hadronically decaying boson is reconstructed as two resolved jets, and the case where it is reconstructed as a single jet. The transverse momentum distribution of the hadronically decaying boson is used to search for new physics. Observations are consistent with the Standard Model predictions, and confidence intervals are calculated for parameters describing anomalous triple gauge-boson…
| Signal | ||
|---|---|---|
| Total Expected Signal | ||
| Background | ||
| Single top-quark | ||
| Multijet | ||
| - | ||
| Total Expected Background | ||
| Total SM Expected | ||
| Observed | ||
| () | ||
| () |
| Lepton | with and | |
|---|---|---|
| Jet | with , , | with , , |
| No small- jets with , , | ||
| , | ||
| Global | ||
| Source of uncertainty | Relative uncertainty for |
|---|---|
| Top-quark background modelling | 13% |
| Signal modelling | 12% |
| modelling | 4% |
| Multijet background modelling | 1% |
| Small- jet energy/resolution | 9% |
| Other experimental (leptons, pile-up) | 4% |
| Luminosity | 2% |
| MC statistics | 9% |
| Data statistics | 14% |
| Source of uncertainty | Relative uncertainty for |
|---|---|
| modelling | 60% |
| Top-quark background modelling | 32% |
| Signal modelling | 15% |
| Multijet background modelling | 13% |
| Large- jet energy/resolution | 45% |
| Small- jet energy/resolution | 16% |
| Other experimental (leptons, pile-up) | 3% |
| Luminosity | 2% |
| MC statistics | 19% |
| Data statistics | 33% |
| Form factor | Parameter | Observed | Expected | Observed | Expected |
|---|---|---|---|---|---|
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| [ ] | [ ] | [ ] | [ ] | ||
| Parameter | Observed | Expected | Observed | Expected |
|---|---|---|---|---|
| [ ] | [ ] | [ ] | [ ] | |
| [ ] | [ ] | [ ] | [ ] | |
| = | [ ] | [ ] | [ ] | [ ] |
| Parameter | Observed [] | Expected [] | Observed [] | Expected [] |
|---|---|---|---|---|
| [ ] | [ ] | [ ] | [ ] | |
| [ ] | [ ] | [ ] | [ ] | |
| [ ] | [ ] | [ ] | [ ] | |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
\AtlasTitle
Measurement of production with the hadronically decaying boson reconstructed as one or two jets in collisions at with ATLAS, and constraints on anomalous gauge couplings
\AtlasRefCodeSTDM-2015-23 \PreprintIdNumberCERN-EP-2017-060 \arXivId1706.01702 \AtlasJournalRefEur. Phys. J. C77 (2017) 563 \AtlasDOIdoi:10.1140/epjc/s10052-017-5084-2 \AtlasAbstractThis paper presents a study of the production of or boson pairs, with one boson decaying to or and one or boson decaying hadronically. The analysis uses 20.2\leavevmode\nobreak\ \mbox{fb{}^{-1}} of collision data, collected by the ATLAS detector at the Large Hadron Collider. Cross-sections for production are measured in high- fiducial regions defined close to the experimental event selection. The cross-section is measured for the case where the hadronically decaying boson is reconstructed as two resolved jets, and the case where it is reconstructed as a single jet. The transverse momentum distribution of the hadronically decaying boson is used to search for new physics. Observations are consistent with the Standard Model predictions, and confidence intervals are calculated for parameters describing anomalous triple gauge-boson couplings.
\size@chapter\sectfont
Contents
@afterheading@starttoc
toc
1 Introduction
Measurements of the production of two massive vector gauge bosons (hereafter, “diboson” production) represent an important test of the Standard Model (SM) of particle physics. Diboson measurements are powerful probes of the electroweak theory of the SM, in particular the structure of the triple gauge-boson couplings (TGCs) [1, 2]. In addition, precise diboson measurements are a valuable test of higher-order calculations in quantum chromodynamics (QCD).
Measurements of and production in the leptonic channels and () have been performed by the ATLAS and CMS collaborations in collisions at and [3, 4, 5, 6, 7, 8, 9], and by the Tevatron experiments in collisions [10, 11, 12, 13]. Measurements in the semileptonic channel () have been performed by ATLAS [14] and CMS [15] at , and by the Tevatron experiments in collisions [16, 17]. The semileptonic channel offers features complementary to the leptonic channels. On the one hand, the presence of jets and the large background from and production limit the experimental precision. On the other hand, the semileptonic channel has an approximately six times higher branching fraction than the fully leptonic channels. Also, for , the original diboson kinematics can be better reconstructed in an final state than in an final state, since the latter has two invisible particles, rather than only one in . Both of these advantages are particularly beneficial for searching for beyond-the-Standard-Model (BSM) enhancements of diboson production due to heavy new particles, which could modify the diboson spectrum at high transverse momentum () of the bosons [18].
It is possible to reconstruct the decay as two small-radius jets (“small-” jets, denoted by ) or as a single large-radius jet (“large-” jet, denoted by ). Reconstructing the decay as a large- jet enables an increased reconstruction efficiency at high , thus improving the sensitivity to BSM signals. In addition, by applying grooming [19] techniques such as trimming [20] to the large- jets, it is possible to better distinguish events containing decays from background events [21].
In this paper, measurements of fiducial cross-sections are presented in phase spaces containing a candidate with high . Two fiducial cross-sections are measured, in phase spaces chosen to closely match the two experimental selections used in this paper. The first event selection, denoted , reconstructs the decay as two small- jets, while the second one, denoted , reconstructs the as a single large- jet. Previous cross-section measurements of have not exploited large- jets.
A search for anomalous triple gauge-boson couplings (aTGCs) is also presented in this paper, using both the and channels. Previous searches for charged aTGC contributions to production have been conducted by the ATLAS Collaboration [14] using collisions, by the CMS Collaboration [15, 22] using and collisions, and by the D0 [23] and CDF [24] collaborations using collisions. Most published aTGC searches in the channel have reconstructed the as two small- jets, with the exception of Ref. [22], which reconstructed the as a single large- jet.
2 Analysis overview
As mentioned above, measurements of production are performed using either two small- jets or a single large- jet to reconstruct the hadronically decaying boson. For both channels, the leptonically decaying boson is reconstructed by requiring the presence of a lepton (electron or muon) and missing transverse momentum.
After applying stringent event selection requirements, the signal-to-background ratio remains quite low at 5–10, because of the large W + jets background. In order to distinguish the SM signal from the background, the dijet mass distribution (in the channel) or the mass distribution of the large- jet (in the channel) is used as a discriminating variable. The signal events peak near the mass in these distributions, while the shape of the dominant background is smoothly falling. In both channels, the signal is extracted from a fit to the discriminating variable. Wide fitting ranges are used, in order to allow the backgrounds to be constrained by the data.
A fiducial cross-section is measured separately in the and the channel; the fiducial phase spaces for the measurements are defined to be close to the experimental event selections. The fiducial cross-section in each channel is extracted from the previously mentioned fits. The events in the two channels partially overlap, because there are some events for which the decay can be reconstructed both as two small- jets and as one large- jet. In order to simplify the interpretation of the results and allow easier comparison with theoretical predictions, the overlap events are not removed, and both measurements are presented separately. No combination of the and cross-section measurements is performed. The electron and muon channels are combined when performing the measurements, since little improvement in sensitivity is expected from separating by lepton flavour. Event kinematics and the signal-to-background ratio are similar in the electron and muon channels, and the dominant sources of uncertainty are unrelated to lepton flavour.
A search for aTGC contributions is also performed in the and channels. The event selection is the same as for the cross-section measurements, except that a tighter requirement is made on the dijet mass or on the mass of the large- jet. The search is performed by fitting the distribution of the dijet system ( channel) or of the large- jet ( channel). These distributions are sensitive to aTGCs, which are expected to lead to deviations from the SM prediction at high .
3 ATLAS detector
The ATLAS detector [25], which surrounds one of the interaction points of the Large Hadron Collider (LHC) [26], is built of several subdetectors. The first subdetector layer consists of the inner detector (ID), which provides charged-particle tracking for .111ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the -axis along the beam pipe. The -axis points from the IP to the centre of the LHC ring, and the -axis points upward. Cylindrical coordinates are used in the transverse plane, being the azimuthal angle around the -axis. The pseudorapidity is defined in terms of the polar angle as . The ID is further subdivided into (ordered from innermost to outermost) a pixel detector, a silicon-microstrip tracker, and a transition radiation tracker. Surrounding the ID there is a superconducting solenoid that provides a 2 T magnetic field. Outside of the solenoid, there is an electromagnetic (EM) calorimeter based on liquid-argon technology, which provides coverage up to . Additionally, a scintillator-tile calorimeter provides hadronic energy measurements in the range , and liquid-argon-based endcap and forward calorimeters extend the EM and hadronic measurements up to . A muon spectrometer, consisting of tracking and triggering detectors and three toroidal magnets, surrounds the calorimeters; it provides muon tracking and identification up to and triggering capability up to .
A three-level trigger system is used to select the most interesting events for data storage [27]. An initial hardware-based trigger stage is followed by two software-based triggers, which reduce the final event rate to about 400 Hz.
4 Data and Monte Carlo samples
This analysis is based on an integrated luminosity of 20.2\pm 0.4\leavevmode\nobreak\ \mbox{fb{}^{-1}} of collisions recorded by the ATLAS detector in 2012. Events are required to pass one of several single-lepton triggers. The triggers require either an isolated electron or muon with , or an electron (muon) having without an isolation requirement.
The nominal signal Monte Carlo (MC) samples consist of events generated at next-to-leading order (NLO) in QCD using MC@NLO v4.07 [28] interfaced with Herwig v6.520 [29] and Jimmy v4.31 [30] for the simulation of parton showering, hadronization, and the underlying event. The CT10 parton distribution function (PDF) set [31] and parameter values from the AUET2 tune [32] are used for these samples. The and bosons are generated on-shell by MC@NLO and decayed subsequently by Herwig. The same MC configuration is also used to model aTGC contributions to production, using an event reweighting feature built into MC@NLO.
In order to study systematic uncertainties, alternative samples are generated at NLO in QCD with Powheg-Box [33, 34, 35] using the CT10 PDF set. The parton showering and hadronization is modelled with Pythia 8.175 [36] using the AU2 tune [37]. Off-shell and decays are included; the decays have a requirement of and .
Another set of alternative samples are generated with Sherpa v1.4.1 [38, 39, 40, 41]. These samples are generated at leading order (LO) in QCD, but include up to three additional partons in the matrix element. Off-shell and decays are included; the decays have a requirement of and .
Contributions from are only at the level after applying the full event selection and are thus neglected. Signal MC samples for non-resonant production are not used in the analysis, but the contribution from this process is estimated as described in Section 10, and included in the final cross-section predictions.
The and backgrounds (collectively referred to as ) are modelled at LO in QCD with Sherpa v1.4.1, with up to four additional final-state partons. The CT10 PDF set is used for these samples, and they are normalized using inclusive cross-sections that are next-to-next-to-leading order (NNLO) in QCD, obtained using FEWZ [42]. For studies of systematic uncertainties, alternative samples are generated with Alpgen [43] interfaced with Pythia 6.426 [44], modelling the process at LO in QCD with up to five final-state partons. These additional samples use the Perugia 2011C tune [45] and the CTEQ6L1 PDF set [46].
The MC samples for the and single-top-quark (-channel, -channel, and ) processes (collectively referred to as top-quark processes) are generated with Powheg-Box [47, 48, 49] interfaced with Pythia 6.426 [44] (or Pythia 6.427 for the -channel single-top-quark process). All of these samples use the CT10 PDF set for the matrix element, the CTEQ6L1 PDF set for the parton shower, and the Perugia 2011C tune.
The background process is modelled with Powheg interfaced with Pythia 8. The sample is normalized using the NLO prediction from MCFM [50, 51].
The MC samples are passed through a GEANT4-based [52] simulation of the ATLAS detector [53]. For some of the MC samples, a fast simulation is used that makes use of a parameterization of the showers in the calorimeter. The hard-scattering processes in the MC samples are overlaid with simulated minimum-bias events in order to model additional collisions in the same or neighbouring bunch crossings (“pile-up”). The MC samples are reweighted so that their pile-up profile matches that observed in the data.
5 Event reconstruction
This analysis considers events with exactly one lepton (electron or muon), missing transverse momentum, and either two small- jets or one large- jet.
In each event, primary vertices are reconstructed, which must be formed from at least three tracks with . In case an event has multiple primary vertices (due to pile-up), the primary vertex with the highest of the associated tracks is defined as the hard-scatter vertex.
Electron candidates are formed from energy clusters in the EM calorimeter matched to ID tracks. They are required to have and . Candidates in the transition region between the barrel and endcaps of the EM calorimeter, , are excluded. In order to ensure that the electron candidates are consistent with having been produced at the hard-scatter vertex, the transverse impact parameter and longitudinal impact parameter are required to satisfy and mm, respectively, where is the uncertainty in the measured . Both and are measured with respect to the hard-scatter vertex. Electron candidates must also satisfy the “tight” cut-based identification criteria from Ref. [54], based on track parameters and on the shower shapes in the calorimeter. Candidates must also pass isolation requirements based on calorimeter and track measurements. The calorimeter isolation requires , where is defined as the scalar transverse energy sum of the calorimeter energy deposits within a cone centred on the electron candidate (excluding transverse energy from the candidate itself), divided by the of the electron candidate. Similarly, the track isolation requires , where is the scalar sum of the of the tracks within a cone centred on the electron candidate (excluding the of the candidate’s track itself), divided by the electron candidate’s .
Muon candidates are formed from the combination of a track in the muon spectrometer and one in the ID. They are required to have and . Their impact parameters must satisfy and mm. The candidates must also satisfy the isolation criteria and , where and are defined analogously to the electron case.
Small- jets are reconstructed from topological energy clusters [55] in the calorimeter using the anti- algorithm [56] with radius parameter . The jet energies are calibrated as described in Ref. [57] and are corrected for pile-up. They are required to have and for the channel. Small- jets with are used in the channel as part of a jet veto (see Section 6). In order to remove jets originating from pile-up, small- jets having and are required to have an absolute value of the “jet vertex fraction” variable (JVF) [58] greater than 0.5.
In the channel, large- jets are reconstructed using the anti- algorithm with radius parameter , and are trimmed [20] using a subjet radius of and a momentum-fraction parameter ; the trimming procedure discards soft subjets from the large- jets and reduces their sensitivity to pile-up [21]. They are required to have and . The energies of the small- and large- jets and the masses of the large- jets are calibrated using - and -dependent scale factors [57, 59].
If an electron and a muon candidate share the same ID track, the electron candidate is rejected. If a small- jet is within of a selected electron candidate, the jet is rejected; if the jet is within of a selected electron, the electron candidate is rejected. Muon candidates are rejected if they are within of a small- jet. Finally, large- jets are rejected if they are within of a selected lepton candidate. In the object selection stage, small- jets and large- jets are allowed to overlap; however, in the event selection stage a requirement is applied between the small- and large- jets, as explained in Section 6.
The missing transverse momentum is computed as the negative vector sum of the transverse momentum of all the detected objects in the event, including reconstructed jets, photons, electrons, and muons. An additional “soft term” is included that accounts for the of clusters in the calorimeter which are not associated with any specific reconstructed object [60]. The magnitude of is denoted .
6 Event selection
Two independent sets of event selection criteria are developed that target different event topologies: the selection, described in Section 6.1, and the selection, described in Section 6.2. The channel and channel differ significantly from one another in their kinematics, expected signal yields, and signal-to-background ratios. Therefore, the event selection criteria are optimized separately for the two channels.
For both the and selections, all events are required to contain at least one primary vertex. Events must have exactly one good electron or muon candidate. Events are vetoed if they contain any additional lepton candidates that have and satisfy a looser set of selection criteria.
6.1 channel
Events must have and a transverse mass222The transverse mass is defined as , where is the transverse momentum of the lepton candidate. . Events must contain exactly two small- jets. The requirement of exactly two jets substantially reduces the background from top-quark decays. The pseudorapidity separation of the selected jets is required to satisfy , in order to improve the signal-to-background ratio.
In order to reduce the multijet background not removed by the requirement, an azimuthal-angle difference between the direction and the direction of the leading- jet of is required. Also, both the and candidates must pass requirements on their transverse momenta: and , where . These requirements enhance the separation between the signal and background distributions in the dijet mass.
As described in Section 8, the signal is extracted using a maximum-likelihood (ML) fit to the dijet mass () distribution. In the dijet mass calculation, the mass of each individual jet is set to zero, which makes the variable easier to model in the MC simulation. Since the signal is extracted from a fit to , only a loose requirement is made on this variable: .
6.2 channel
Events must contain exactly one large- jet with and . The backgrounds from top-quark decays are suppressed by rejecting events containing any small- jets with and that are separated from the large- jet by . In order to suppress the multijet background, a requirement of is applied. The trimmed mass of the large- jet, , must be , and the signal is measured from the ML fit to .
Since the and event selections are done independently, some events pass both selections. About of the signal MC events that pass the selection also pass the selection, while about of the signal MC events that pass the selection also pass the selection.
7 Background estimation
The methods for estimating the expected background yields and kinematic distributions are described in this section. The estimates from this section are used as inputs to the ML fit in which the signal is measured while the backgrounds are allowed to vary within their systematic uncertainties. In that ML fit, the normalization is allowed to vary without constraint, so the estimates given in this section are pre-fit estimates.
Most of the backgrounds (, , , single top-quark, and ) are estimated using MC simulation, with data-driven corrections applied in some cases, as described later in this section. By far the largest background in the analysis is from , followed by top-quark production. Despite the latter background’s subdominant contribution, it plays an important role because it contains contributions from real decays, which make it more difficult to distinguish from the signal. About of the top-quark background is due to production, and the remainder comes from single-top-quark production.
Multijet processes form another source of background. Multijet events can pass the event selection if they contain non-prompt leptons (produced from semileptonic decays of - and -hadrons) or “fake” leptons (resulting from misidentified jets). The multijet backgrounds are estimated using data-driven techniques, as described in Sections 7.1 and 7.2.
7.1 channel
The background prediction is MC-based, but data-driven corrections are applied to the MC prediction in order to improve the description of the jet kinematics. A control region (CR) is defined identically to the signal region, except that the region is vetoed, in order to remove most of the signal events. One-dimensional reweighting functions of the variables and are derived from this CR. These reweighting functions have approximately effects on the shapes of the and distributions. Data–MC comparisons in the CR are shown in Figure 1, before and after application of the reweighting functions. All further results in this paper are shown with these two reweighting functions applied to the MC samples. The same reweighting functions are used for both the and processes. It was checked that the reweighting functions obtained from the low- and high- portions of the control region are compatible.
The top-quark background is modelled with MC simulation, and is cross-checked in a validation region containing three small- jets, one of which is -tagged using the MV1 algorithm [61, 62]. Good agreement is observed between the data and the MC simulation, so no corrections are applied to the prediction. The background from events is also modelled with MC simulation.
The data-driven multijet background estimate makes use of a multijet CR. The multijet CR is formed by selecting events in data that pass the same selection requirements as for the signal region, except that the lepton quality criteria are modified in order to produce a CR enriched in non-prompt and fake leptons. Lepton candidates satisfying these modified criteria are called “anti-identified” lepton candidates. Anti-identified muon candidates must have a non-negligible impact parameter, , and satisfy looser isolation criteria than the signal muon candidates. Anti-identified electrons must fail the “tight” but satisfy the “medium” cut-based identification criteria from Ref. [54], and are also required to contain a hit in the innermost layer of the pixel detector. In addition, the isolation criteria are modified for anti-identified electron candidates, in order to enrich the sample in non-prompt and fake electrons.
The shapes of the kinematic distributions (such as , , and ) of the multijet background are estimated from events in the multijet CR, after subtracting the MC predictions of the non-multijet contributions to the CR. These non-multijet contributions are about 20% (50%) of the total in the electron (muon) channel. The overall multijet background event yield is estimated from a fit to the distribution of events that pass the full signal region selection, except that the requirements on and (and also and for the muon channel) are removed in order to enhance the number of multijet events. This selection is referred to as the extended signal region. In this fit, the multijet shape is estimated from an extended multijet CR, defined analogously to the extended signal region, but requiring the lepton to pass the anti-identified-lepton selection. The shapes of the other backgrounds are estimated using MC samples. The multijet event yield obtained from this fit is then extrapolated to the signal region, using the ratio of events in the multijet CR and the extended multijet CR, corrected for non-multijet contributions. The multijet background estimates are performed separately for the electron and muon channels. Only about of the total multijet background is in the muon channel.
The expected signal and background yields in the signal region are given in Table 1, and compared to the number of events observed in data. The predictions for the distribution shapes of the signal and backgrounds are shown in Figure 2.
7.2 channel
In the channel, the , , and top-quark backgrounds are estimated using MC samples. The MC predictions for the two largest backgrounds ( and top-quark production) are corrected by scale factors obtained from dedicated control regions.
The top-quark control region (top CR) is formed by events satisfying the signal region selection, except that the presence of at least one small- -tagged jet with and is required instead of applying the nominal veto on small-R jets. The jets are -tagged using the MV1 algorithm [61, 62], using a working point with a -tagging efficiency of about 70% and a gluon/light-quark jet rejection factor of over 100 in events. About of the events in this top CR originate from top-quark backgrounds. There is a deficit in data in the top CR relative to the MC prediction, which is attributed to a mismodelling of the top-quark backgrounds. A global scale factor of for the top-quark backgrounds is obtained from this CR, after subtracting the prediction for non-top-quark backgrounds. The data in the top CR is shown in Figure 3, compared to the SM prediction after application of the top-quark scale factor. This scale factor is applied to the top-quark background predictions in the signal region.
The control region for the background ( CR) is obtained by applying the standard signal region selection, but adding the requirement that or . This additional requirement removes almost all of the signal events and also a large fraction of the top-quark events. About of the events in this CR originate from backgrounds. The top-quark background prediction in the CR is scaled by the top-quark scale factor obtained above. A data deficit is observed in the CR relative to the prediction. A global scale factor of is obtained for the background, after subtracting the expected contributions from the other signal/background processes. A comparison between the data and the prediction in the CR is shown in Figure 3, after application of the scale factor. The scale factor is applied to the prediction in the signal region.
The method for estimating the multijet background is similar to that used in the channel. As in the channel, a multijet CR is defined by requiring an “anti-identified” lepton candidate. The shapes of the kinematic distributions are estimated from this CR using the same method as in the channel. The non-multijet background contributions to the CR are about 6% of the total. The multijet event yield is estimated from a fit to the distribution, as in the channel, but the only requirement that is removed for the definiton of the extended signal region/multijet CR is the requirement. The multijet background is found to be negligible for the muon channel, so only the contribution in the electron channel is considered for the final results.
The numbers of expected and observed events in the signal region are summarized in Table 1. The previously mentioned top-quark and scale factors are applied to the predictions. The contribution from events is expected to be very small in the channel, so it is neglected. The nominal predictions for the distribution shapes of the signal and backgrounds are shown in Figure 2.
8 Cross-section extraction
The fiducial cross-section for production is measured independently for the and phase spaces, in both cases using the formula:
[TABLE]
where is the measured signal yield, is the integrated luminosity, and is a factor that corrects for experimental acceptance and efficiencies. Since this analysis measures as the sum of the and processes, which can each have different acceptances and efficiencies, is given by:
[TABLE]
where the are reconstruction correction factors and the variable is the predicted ratio of the fiducial cross-section to the fiducial cross-section. The and values are estimated from MC simulation. The factors are defined as the predicted number of signal events passing the reconstruction-level event selection divided by the number of events in the fiducial phase space defined with generator-level particles. The factors account for reconstruction inefficiencies, resolution effects, and for contributions to the signal region from events that do not decay to (such as or ); the latter are included in the numerator and not in the denominator. The cross-section is measured for the sum of the electron and muon channels, so is computed as a weighted average over the electron and muon channels. The fiducial cross-section measurement therefore assumes that the signal MC simulation correctly predicts the ratio of to and of electrons to muons. The value of is in the channel and in the channel, including systematic uncertainties (see Section 9).
The fiducial phase spaces for the and channels are defined in Sections 8.1 and 8.2, respectively. These fiducial phase spaces partially overlap. In order to cope with the small signal-to-background ratios in this analysis (5–10), the cross-section is extracted using a binned ML fit to the distribution (in the analysis) or the distribution (in the analysis). The ML fits are performed on the sum of the electron and muon channels. It was cross-checked that the electron and muon channels are compatible, in both the and channels.
In the ML fits, the value of and the background yield are both free to vary without constraint. Systematic uncertainties in the signal and backgrounds are incorporated in the fit by including nuisance parameters that are allowed to vary within prior constraints. The nuisance parameters allow the luminosity, , the non- background yields, and the and shapes of the signal and background distributions to vary within their systematic uncertainties. The correlations between the uncertainty in and the uncertainty in the signal shapes are accounted for in the fit. The sources of systematic uncertainty and the methods to assess these uncertainties are described in detail in Section 9.
8.1 fiducial phase space
The fiducial phase space is defined to closely match the experimental event selection. The phase-space definition requires a pair with the bosons decaying as and , where . Events containing other kinds of decay channels (such as events or with the decaying to ), are not included in the fiducial phase-space definition. Such events can still pass the experimental event selection (where they are included in the signal category), and they are accounted for in the definition.
Leptons selected in the fiducial region must have and . The four-momentum of the lepton is modified by adding to it the four-momenta of all the photons within , excluding photons produced by hadron decays. Particle-level anti- jets are constructed using as constituents all stable particles, excluding muons and neutrinos. Stable particles are defined as those having a mean lifetime of ps. The particle-level jets must have and . Jets within of a selected electron are rejected, and then leptons within of a remaining jet are rejected. The true in the event is defined as the magnitude of the vector sum of all the neutrinos.
The event must have exactly one lepton and two jets matching the above definitions. The remaining requirements for the fiducial phase space are summarized in Table 2, and are analogous to the experimental event selection, but are defined using the lepton, , and particle-level jets described in this section.
8.2 fiducial phase space
As in the channel, the fiducial phase-space definition requires a pair with and . Leptons, , and particle-level jets are defined in the same way as in the channel, except that two sets of leptons and small- jets are considered: central leptons (small- jets) are required to have (), and extended leptons and small- jets are required to have . Particle-level large- jets are defined by applying the anti- algorithm with radius parameter to all stable particles, excluding muons and neutrinos. No trimming is applied to these jets. The large- jets are required to have and . Central (extended) small- jets that are within of a central (extended) electron are rejected. Then, central leptons are rejected if they are within of a remaining central small- jet. Large- jets are rejected if they are within of any remaining central leptons. Events are required to contain exactly one central lepton and one large- jet with the above definitions, and events are discarded if they contain any extended small- jets with . The event must also have , and the large- jet must have a mass greater than . The fiducial phase-space definition is summarized in Table 2.
9 Systematic uncertainties
Systematic uncertainties in the measured can be due to uncertainties in , , and/or . Uncertainties in the measured can in turn be due to uncertainties in the background yields or in the shapes of the kinematic distributions (, ) of the signal and backgrounds (hereafter called “shape uncertainties”). The dominant systematic uncertainties in the measurement are those affecting the measured .
A wide variety of detector-related experimental uncertainties are considered, which affect , the predicted background yields, and the signal and background shapes. The most important of these uncertainties are those related to the jet reconstruction. Uncertainties in the small- jet energy scale and resolution are accounted for [57, 63]. In the channel, uncertainties in the large- jet energy and jet mass scales are also taken into account. The scale uncertaities of the large- jets are estimated using a double-ratio method that compares calorimeter- and track-jets in data and MC simulation [21]. The energy and mass resolution uncertainties of large- jets are estimated by smearing the jet energies/masses so as to degrade the resolutions by 20%; this approach is based on prior studies of large- jets [64, 65]. The systematic uncertainty due to the JVF requirement is also included [66]. In addition to the jet-related uncertainties, there are also systematic uncertainties in the electron and muon reconstruction (including triggering, object reconstruction, identification, and the energy scale and resolution) [54, 67, 68, 69, 70]. The effects of the jet and lepton uncertainties are propagated to the calculation, and an additional systematic uncertainty in the soft terms entering the calculation is also included [60].
In the cross-section fits, the yield is taken to be a free parameter, while several uncertainties in the modelling of its shape are accounted for (in addition to the shape uncertainties from the previously mentioned detector effects). Systematic uncertainties in the shape are estimated by varying the MC event generator used (Sherpa compared to Alpgen+Pythia). The differences between the predictions of the two generators are taken as additional systematic uncertainties. Additional uncertainties in the shape are estimated by varying the renormalization and factorization scales by factors of and , and by varying the scale used in Sherpa for matching the matrix elements to the parton showers [39] from its nominal value of to alternative values of and . In the channel, the uncertainty in the shapes of the predictions due to the two kinematic reweighting functions (see Section 7.1) is estimated by including the full difference between applying and not applying each reweighting function as additional systematic uncertainties. In the channel, an uncertainty of in the cross-section ratio is also included; this uncertainty is ignored in the channel as it has a negligible effect.
For the background, uncertainties due to the matrix-element event generator, parton shower/hadronization model, and amount of initial- and final-state radiation are all included. The theoretical uncertainties in the top-quark background cross-sections are also taken into account. In the channel, instead of using the theoretical cross-section uncertainty, the top-quark background is assigned a normalization uncertainty of to account for the uncertainty in the data-driven scale factor. Systematic uncertainties in the multijet background estimate are also included, which affect both its normalization and its shape. These uncertainties are derived from studies of variations of the data-driven estimate, such as changing the control region definitions and varying the non-multijet background subtraction. The uncertainty in the multijet yield amounts to () for the electron (muon) channel in the channel. In the channel, an uncertainty of is assigned to the multijet yield in the electron channel, while the multijet background is neglected in the muon channel. A uncertainty is assigned to the event yield in the channel, to account for uncertainties in the cross-section and the extrapolation to the fiducial phase space.
Additionally, the uncertainty in the modelling of pile-up interactions is accounted for [71]. The uncertainty in the integrated luminosity is also included, computed as described in Ref. [72]. The statistical uncertainty of the MC samples is taken into account, which affects each bin in the ML fits in an uncorrelated way.
Uncertainties in the signal shapes and in the parameter due to variations of the signal model are computed by varying the renormalization and factorization scales by factors of and , and by comparing the nominal MC@NLO signal samples to alternative samples generated with Sherpa and Powheg +Pythia 8. The effect on from the uncertainties in the CT10 PDF set is also taken into account; the PDF uncertainty has a negligible impact on the signal shapes.
The measured values are compared to theoretical predictions from MC@NLO. The uncertainty in the theoretical prediction is calculated including the uncertainties due to renormalization and factorization scales. Since the fiducial phase spaces contain a veto on additional jets, the Stewart–Tackmann procedure [73] is used to estimate the scale uncertainties. These uncertainties are also propagated to the theoretical value which enters into the calculation, although the effect of this on the measured is very small (). PDF-induced uncertainties in the theoretical prediction are also taken into account.
10 Cross-section results
The result of the ML fit to the distribution for the channel is shown in Figure 4. The fit is performed on the sum of events in the electron and muon channels. The observed significance is , including statistical and systematic uncertainties,333The significance is calculated based on the profile-likelihood ratio of the background-only and signal-and-background hypotheses. This ratio is converted to a significance using the asymptotic approximation [74]. while the expected significance, calculated using the Asimov data set [74], is . The fitted background normalization is times its pre-fit value, while the fitted top-quark background normalization is times its pre-fit value.
The fiducial cross-section for the signal process is extracted from the fit as described in Section 8, and the result is
[TABLE]
The impacts of the various systematic uncertainties on the cross-section measurement are shown in Table 3. The measurement can be compared to the theoretical prediction of
[TABLE]
The theoretical prediction is obtained using MC@NLO for the prediction. The prediction is also included, and is calculated using the total NLO cross-section prediction [75] multiplied by the acceptance from MC@NLO. The contribution increases the fiducial cross-section prediction by in both the and channels. Given the relatively small contribution, the possible differences in acceptance between the and processes are neglected. The uncertainty in the MC@NLO prediction is described in Section 9.
The result of the fit for the channel is shown in Figure 5. Although the signal-to-background ratio is better in this case than in the channel, the total number of signal events is much smaller. The observed significance of the result is (including statistical and systematic uncertainties), compared to an expected significance of . The fitted (top-quark) background normalization is () times its pre-fit value.
The extracted fiducial cross-section for the signal process is
[TABLE]
which is compatible with the theoretical prediction of
[TABLE]
The breakdown of the uncertainties contributing to the fiducial cross-section measurement is shown in Table 4.
The cross-section measurements are summarized in Figure 6. As mentioned in Section 8, the two cross-section measurements are performed in partially overlapping phase spaces. The uncertainty in the theory prediction is significantly larger in the channel than in the channel. The theoretical uncertainty in the channel is dominated by the scale uncertainties, which are particularly large because of the aggressive jet veto in this channel (only about 30% of signal MC events pass the jet veto in the channel, compared to about 80% in the channel).
11 Constraints on anomalous gauge couplings
In many extensions of the SM, diboson production can be modified, such as through new resonances that couple to bosons. If the scale of new physics is sufficiently high, new resonances may not be visible in the current data; however, diboson production could still be affected below the new-physics scale, in the form of modified couplings. One common framework for parameterizing new physics in diboson production is an effective Lagrangian [1] of the form:
[TABLE]
where or , , and . The six parameters , , and (hereafter called “aTGC parameters”) are all zero in the SM. The parameter is zero because of EM gauge invariance, leaving five free aTGC parameters, which describe deviations of the triple gauge-boson couplings from their SM predictions. It is common to apply the so-called LEP constraint [76], which imposes gauge invariance, and which reduces the number of independent aTGC parameters to three, by introducing the following constraints: and , where is the weak mixing angle. Since aTGC parameters lead to violation of unitarity at high energies, form factors are often applied to them in order to ensure unitarity:
[TABLE]
where is one of the aTGC parameters, is the square of the diboson invariant mass, and is the form factor’s energy scale.
An alternative framework for describing modifications of diboson production is an effective field theory (EFT) [77, 78] that is assumed to be valid below an energy scale , and which introduces three CP-conserving dimension-six operators:
[TABLE]
Here, is the Higgs doublet field, is the covariant derivative, and and are the field strength tensors of the and gauge boson fields. The coefficients of these operators (EFT parameters), , , and , are zero in the SM and can be related to the LEP-constraint aTGC parameters as follows:
[TABLE]
This relation only holds if no form factor is applied to the aTGCs. The effect of aTGC/EFT parameters on the process is neglected.
The aTGC and EFT parameters both tend to increase the diboson cross-section at high and high invariant mass of the diboson system. Both the channel and the channel can be used to search for these BSM enhancements. The channel, although currently less sensitive as a SM measurement, is expected to provide a higher sensitivity to the aTGC/EFT models, because of the better efficiency at high . On the other hand, the channel, where the SM measurement is clearly established, is useful as a complementary search channel that probes a different energy range.
In this analysis, the new-physics search uses signal regions with exactly the same event selection as the cross-section measurements, except that the requirement is tightened to in the channel and the requirement is tightened to in the channel. These tighter requirements lead to an increase in the signal-to-background ratio. In the channel, events which fail the requirement (i.e. or ) are put into a sideband control region. The background is neglected in the new-physics search, due to its very small expected contribution.
The search makes use of the ( channel) or ( channel) distribution. Hereafter, is used to refer to both and . The distributions of the events in the signal regions are shown in Figure 7. This figure also shows the expected enhancement at high in the presence of different EFT parameter values. As can be seen from the figure, no significant deviation from the SM prediction is observed; therefore, confidence intervals are computed for the aTGC and EFT parameters.
The confidence intervals are computed from binned ML fits to the distributions. The intervals are calculated using a frequentist Feldman–Cousins approach [79]. In the channel, simultaneous fits to the distributions in the signal region and sideband CR are used, while in the channel, only the distribution in the signal region is used. Since the and selections overlap, the confidence intervals are calculated separately for the and selections. In the fits, the SM and background predictions are allowed to vary within their uncertainties. The measured cross sections of Section 10 are consistent with theoretical SM predictions, but have large associated uncertainties; for this reason the theoretical prediction is used here. The systematic uncertainties in the normalizations and shapes of the signal and backgrounds are accounted for through nuisance parameters. The systematic uncertainties that have the largest impact on the results are the jet-related uncertainties (in both channels) and the uncertainty from the limited size of the MC samples (in the channel).
The observed 95% confidence intervals for the aTGC parameters are shown in Table 5, without applying the LEP constraint. The confidence intervals for a given aTGC parameter are computed while fixing the other aTGC parameters to zero. The confidence intervals are shown separately for the and selections, and the expected confidence intervals under the SM hypothesis are also shown for comparison. Confidence intervals for the aTGC parameters are shown for and for the case of no form factor (equivalent to ). The value of is chosen in order to ensure unitarity over the range of aTGC parameter values to which this analysis is sensitive [80].
The selection has significantly better sensitivity to the aTGC parameters. No combination of the and constraints is performed, since it is expected that the channel would dominate the combination. The sensitivity to the aTGC parameters in the channel mainly comes from the bins, whereas the sensitivity in the channel mainly comes from the 300–600 GeV bins. Since the channel probes a lower range, its sensitivity is less degraded by the form factors (which have a larger effect at higher ) than the channel.
In addition, the observed and expected confidence intervals for the aTGC parameters in the LEP-constraint scenario are given in Table 6 for .
The observed and expected confidence intervals for the EFT parameters are shown in Table 7, separately for the and selections. Confidence regions for combinations of two EFT parameters are shown in Figure 8; for each combination the third EFT parameter is held fixed to zero. Although the constraints from the channel are less stringent than those from the channel, they probe a complementary phase space. The sensitivity of the channel is similar to the most sensitive previous analyses to publish constraints on these parameters [3, 6, 5, 22]. The channel probes a similar phase space to Ref. [22]; these analyses benefit from their ability to reconstruct high- decays.
12 Conclusion
The production of , with being a or boson, is measured using 20.2\leavevmode\nobreak\ \mbox{fb{}^{-1}} of collisions at at the LHC with the ATLAS detector. The measurements focus on production where the bosons have large transverse momentum. Fiducial cross-sections for the process are measured in two different, but partially overlapping, phase spaces.
The first phase space, denoted , targets a hadronically decaying boson whose decay products can be distinguished as two jets. In this phase space, the process is measured with a significance of , and the fiducial cross-section is measured to be , in agreement with the MC@NLO prediction of .
The second phase space, denoted , contains a single jet consistent with the collimated decay products of a high- boson. The process is measured with a significance of in this phase space. The fiducial cross-section for this phase space is measured to be , consistent with the MC@NLO prediction of .
The events are also used to search for new physics modifying triple gauge-boson vertices, which could lead to enhancements of the cross-section at high of the bosons. No evidence is found for new physics, and confidence intervals are computed for anomalous coupling parameters. The constraints on new physics are also interpreted in terms of an effective field theory. The channel is found to be significantly more sensitive to the new-physics parameters than the channel, which demonstrates the power of large-radius jet substructure techniques. The constraints from this analysis on the new physics parameters are comparable to the previous most stringent constraints from other diboson analyses.
Acknowledgements
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently.
We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DSM/IRFU, France; SRNSF, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, the Canada Council, CANARIE, CRC, Compute Canada, FQRNT, and the Ontario Innovation Trust, Canada; EPLANET, ERC, ERDF, FP7, Horizon 2020 and Marie Skłodowska-Curie Actions, European Union; Investissements d’Avenir Labex and Idex, ANR, Région Auvergne and Fondation Partager le Savoir, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF; BSF, GIF and Minerva, Israel; BRF, Norway; CERCA Programme Generalitat de Catalunya, Generalitat Valenciana, Spain; the Royal Society and Leverhulme Trust, United Kingdom.
The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [81].
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Kaoru Hagiwara, R.D. Peccei, D. Zeppenfeld and K. Hikasa “Probing the weak boson sector in e + e − → W + W − → superscript 𝑒 superscript 𝑒 superscript 𝑊 superscript 𝑊 e^{+}e^{-}\to W^{+}W^{-} ” In Nucl. Phys. B 282 , 1987, pp. 253 DOI: 10.1016/0550-3213(87)90685-7 · doi ↗
- 2[2] Anja Butter et al. “The Gauge-Higgs Legacy of the LHC Run I” In JHEP 07 , 2016, pp. 152 DOI: 10.1007/JHEP 07(2016)152 · doi ↗
- 3[3] ATLAS Collaboration “Measurement of total and differential W + W − superscript 𝑊 superscript 𝑊 W^{+}W^{-} production cross sections in proton–proton collisions at s = 8 Te V 𝑠 8 Te V \sqrt{s}=8\;\mbox{Te V} with the ATLAS detector and limits on anomalous triple-gauge-boson couplings” In JHEP 09 , 2016, pp. 029 DOI: 10.1007/JHEP 09(2016)029 · doi ↗
- 4[4] ATLAS Collaboration “Measurement of the W + W − superscript 𝑊 superscript 𝑊 W^{+}W^{-} production cross section in p p 𝑝 𝑝 pp collisions at a centre-of-mass energy of s 𝑠 \sqrt{s} = 13 Te V with the ATLAS experiment”, 2017 ar Xiv: 1702.04519 [hep-ex]
- 5[5] CMS Collaboration “Measurement of the W + W − superscript 𝑊 superscript 𝑊 W^{+}W^{-} cross section in p p 𝑝 𝑝 pp collisions at s = 8 Te V 𝑠 8 Te V \sqrt{s}=8\;\mbox{Te V} and limits on anomalous gauge couplings” In Eur. Phys. J. C 76 , 2016, pp. 401 DOI: 10.1140/epjc/s 10052-016-4219-1 · doi ↗
- 6[6] ATLAS Collaboration “Measurements of W ± Z superscript 𝑊 plus-or-minus 𝑍 W^{\pm}Z production cross sections in p p 𝑝 𝑝 pp collisions at s = 8 Te V 𝑠 8 Te V \sqrt{s}=8\;\mbox{Te V} with the ATLAS detector and limits on anomalous gauge boson self-couplings” In Phys. Rev. D 93 , 2016, pp. 092004 DOI: 10.1103/Phys Rev D.93.092004 · doi ↗
- 7[7] ATLAS Collaboration “Measurement of the W ± Z superscript 𝑊 plus-or-minus 𝑍 W^{\pm}Z boson pair-production cross section in p p 𝑝 𝑝 pp collisions at s = 13 Te V 𝑠 13 Te V \sqrt{s}=13\;\mbox{Te V} with the ATLAS Detector” In Phys. Lett. B 762 , 2016, pp. 1 DOI: 10.1016/j.physletb.2016.08.052 · doi ↗
- 8[8] CMS Collaboration “Measurement of the WZ production cross section in pp collisions at s 𝑠 \sqrt{s} = 7 and 8 Te V and search for anomalous triple gauge couplings at s 𝑠 \sqrt{s} = 8 Te V” In Eur. Phys. J. C 77.4 , 2017, pp. 236 DOI: 10.1140/epjc/s 10052-017-4730-z · doi ↗
