Study of $Z \to ll\gamma$ decays at $\sqrt s~$= 8 TeV with the ATLAS detector

This paper presents a study of $Z \to ll\gamma~$decays with the ATLAS detector at the Large Hadron Collider. The analysis uses a proton-proton data sample corresponding to an integrated luminosity of 20.2 fb$^{-1}$ collected at a centre-of-mass energy $\sqrt{s}$ = 8 TeV. Integrated fiducial cross-sections together with normalised differential fiducial cross-sections, sensitive to the kinematics of final-state QED radiation, are obtained. The results are found to be in agreement with state-of-the-art predictions for final-state QED radiation. First measurements of $Z \to ll\gamma\gamma$ decays are also reported.


Introduction
The production of  bosons and their decay to lepton pairs at the Large Hadron Collider (LHC) through the Drell-Yan mechanism has been a topic of very fruitful and detailed studies in the LHC experiments, and the precision of the measurements has led over the past decade to impressive theoretical developments, mostly in the area of higher-accuracy quantum chromodynamic (QCD) predictions.
Precise predictions for quantum electrodynamic and electroweak (QED/EW) effects are also of prime importance for key measurements at the LHC, such as that of the -boson mass or of the weak mixing angle; this became very clear already at the end of the Run-1 data-taking period [1].The analysis described in this paper focuses on the final-state radiation of photons in Drell-Yan production of  bosons decaying to an electron or muon pair.This process, denoted in the following as QED FSR (Fig. 1), can be calculated separately from QED initial-state radiation or QED ISR (Fig. 2) and from initial-final-state interference (IFI) thanks to gauge invariance.Owing both to the kinematic selections applied to the final state particles, and to the narrow width of the  boson, the contributions from QED ISR and QED IFI in the fiducial phase space of the measurements described below are expected to be very small [2].
When considering QED radiation, fixed-order partitioning of the calculations has been shown at the time of the Large Electron-Positron Collider (LEP) to be sub-optimal in terms of obtaining the most precise predictions.Figure 3: Predicted ratios of expected differential cross-sections for the invariant mass,   +  , of the positively charged lepton and the photon (left) and of the angular separation, Δ  , between the photon and the closer of the two leptons (right) for electrons compared to muons.The ratios are shown for the bare and dressed lepton definitions described in the text.The predictions are obtained from a high-statistics inclusive -boson sample using Sherpa 2.2.4 and are shown for events containing two high- T leptons and one high- T photon satisfying the fiducial volume selections described in Table 2 in Section 6 (full symbols).Since most of the difference between electrons and muons observed in this case for the bulk of the distributions is due to the higher acceptance in  for the muons, the electron to muon ratios are also shown in the case of identical fiducial selections for both charged leptons (open symbols).
Photon resummation techniques using exponentiation were originally developed by Yennie-Frautschi-Suura (YFS) [3] for inclusive QED FSR in the case of  +  − collisions.They were then expanded to be used for exclusive QED ISR and FSR and for   collisions at the LHC [2].This analysis uses and compares several state-of-the-art implementations, namely KKMChh [4][5][6], Sherpa [7,8], and Photos [9].In the latter case, which applies only to QED FSR, Photos is used within Powheg+Pythia [10][11][12][13][14], as described in Section 3.These Monte-Carlo (MC) tools provide complete simulated events at LHC energies, in which the final state frequently contains several soft/collinear photons near the final-state leptons from the -boson decay, in addition to the single hard large-angle photon emission required by the analysis selection described in Section 4. The YFS approach resides on the idea of resumming the leading soft logarithms to all orders, rather than focusing on the leading collinear terms.These soft logarithms are largely independent of the inner process characteristics and can be calculated from the external particles and their four-momenta only.The big advantage of the YFS formalism is that in addition it allows for a systematic improvement of this eikonal approximation, order-by-order in the QED coupling constant [15,16].As discussed in Section 8, the measurement results presented in this paper cover normalised fiducial differential cross-sections and fiducial integrated cross-sections for  →  decays and a first measurement of  →  decays.
One often labels the leptons in such final-state configurations as "bare" leptons, a notation adopted throughout this paper, even though the meaning of such a bare lepton is by definition dependent on the tool used to describe multiple-photon emissions in QED FSR.One can also define a "dressed" lepton by combining the bare four-momenta of each lepton with that of QED FSR photons collected close to the lepton within a cone of size Δ = 0.1 around it.
Since soft/collinear QED emissions are regularised by the lepton mass, the predictions for the kinematic observables describing the final state will differ significantly in some regions of phase space between bare electrons and bare muons when considering the calculations using exponentiation described above.This is shown in the case of predictions from Sherpa 2.2.4 [8] in Fig. 3 for the shape of the normalised distribution for the invariant mass,   +  , of the positively charged lepton and the photon and for the angular separation, Δ  , between the photon and the closer of the two leptons.For large values of   +  , which are close to the edge of the phase space available for  →  decays, and for large values of Δ  , which correspond to specific configurations of  →  decays where the dilepton pair and the photon are back-to-back, Fig. 3 shows that one expects about 5% fewer  →  decays than  →  decays relative to the bulk of the distribution.It should be noted that most of the difference between electrons and muons observed in these plots for the bulk of the distributions is due to the higher acceptance in  for the muons (see Table 2).
Figure 3 also shows that the dressed lepton definition minimises the differences between electrons and muons.The use of dressed leptons also minimises the impact of the QED FSR radiation on the integrated cross-sections.It is often used for exclusive analyses of e.g./+jet production because it unambiguously assigns the QED FSR photons to the leptons rather than to hadronic jets when performing such exclusive measurements.However, it is not optimal for precise studies of the QED FSR radiation itself.
For the purpose of the analysis described in this paper, all results presented in Section 8 are shown separately for bare electrons and muons with the exception of the low-statistics measurement of  →  final states described in Section 9, in which bare electrons and muons are combined.However, all the dressed-lepton results are published in HEPDATA as well.The measured normalised differential distributions of the invariant mass of the negatively charged lepton and the photon are also published in HEPDATA, together with the charge asymmetry distribution and the full covariance matrix describing the uncertainties in all these measurements with their correlations.
detector surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, electromagnetic and hadron calorimeters, and a muon spectrometer.The inner tracking detector covers the pseudorapidity range || < 2.5.It consists of silicon pixel, silicon microstrip, and transition radiation tracking detectors.Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic (EM) energy measurements with high granularity.A steel/scintillator-tile hadronic calorimeter covers the central pseudorapidity range (|| < 1.7).The end-cap and forward regions are instrumented with LAr calorimeters for both the EM and hadronic energy measurements up to || = 4.9.The muon spectrometer surrounds the calorimeters and is based on three large superconducting air-core toroidal magnets with eight coils each.The field integral of the toroids ranges between 2.0 and 6.0 T m across most of the detector.The muon spectrometer includes a system of precision tracking chambers and fast detectors for triggering.A three-level trigger system is used to select events to be recorded.The first-level trigger is implemented in hardware and uses a subset of the detector information to accept events at a rate of at most 75 kHz.This is followed by two software-based trigger levels that together reduce the accepted event rate to 400 Hz on average depending on the data-taking conditions during 2012.An extensive software suite [18] is used in the reconstruction and analysis of real and simulated data, in detector operations, and in the trigger and data acquisition systems of the experiment.

Data and simulated event samples
The data used in this analysis was collected in 2012 using   collisions at a centre-of-mass energy of 8 TeV.The corresponding total integrated luminosity is 20.2±0.4 fb −1 .The average number of inelastic   interactions produced per bunch crossing for this dataset is 20.7.
The PowhegBox v1 MC programme [10][11][12][13] was used for the simulation of -boson decays at next-toleading order (NLO) in QCD using the CTEQ6L1 PDF set [19].It was interfaced to Pythia8 [14] for the modelling of the parton shower, hadronisation, and underlying event, with parameters set according to the AU2 CT10 tune [20].QED FSR is modelled using the exponentiated multi-photon emission provided by Photos [21].This combination of MC tools, denoted Powheg+Pythia8+Photos, is used as the main MC sample to estimate the effects of detector efficiency and resolution and to unfold the data presented in this paper.Initial-state radiation of photons (QED ISR) and initial-final-state interference (IFI) are not included in the simulation.For QED FSR photons, Photos is configured to include a matrix-element correction for -boson decay which brings it very close in accuracy (better than 10 −3 ) to NNLO QED [9] for the emission of up to two photons.
An alternative prediction for   →  production with up to three additional parton emissions at leading order in the strong coupling is provided by Sherpa 1.4.1 [7] using the CT10 PDF set [22].Sherpa 1.4.1 is configured such that the leading photon is generated at matrix-element level with leading-order accuracy.
Both FSR and ISR emissions, but not the small effects expected from IFI, are included in the simulation.
The YFS [3] approach is then used to exponentiate the radiation in a way similar to that done in Photos.The results obtained by using Sherpa 1.4.1 are compared to those obtained by using Powheg+Pythia8+Photos for a cross-check of the systematics related to the unfolding procedure.A specific version of Sherpa 1.4 (labelled Sherpa 1.4 ME  ) containing an exact description at tree level of the production of lepton pairs accompanied by two hard photons has also been used to compare its predictions to the unfolded  →  measurements.
The pseudorapidity is defined in terms of the polar angle  as  = − ln tan(/2).Angular distance is measured in units of For the comparisons between predictions and the unfolded data, additional generator-level predictions for   →  and   →  production have been produced using Sherpa 2.2.4 [8] and KKMChh [5,6].The QED treatment in Sherpa 2.2.4 is identical to that described above for Sherpa 1.4.1, while the QCD treatment is more advanced [23].KKMChh is an event generator for Z boson production and decays, that includes exponentiated multi-photon initial and final state radiation and virtual electroweak corrections.It separates the expected dominant contribution from QED FSR to the measurements presented in this paper from the negligible contributions from QED ISR and QED IFI.
To minimise the impact of differences in the QCD aspects on the measurements of QED FSR differential spectra, the MC predictions for which data and MC are not expected to agree to better than 0.5% should be reweighted to the data.This procedure was applied for the   T spectrum, which is only predicted at leading order in QCD for the two main MC samples used in the analysis.Other distributions have also been studied to assess the residual systematic uncertainties from the QCD aspects of the simulation in the measurements.This reweighting procedure is also applied to the MC predictions at generator level from Sherpa 2.2.4 and from KKMC-hh.
Background contributions from top-quark production ( t, , and -channel single top production) are generated using PowhegBox with the CT10 PDF in conjunction with Pythia6 [24] (v6.427 using the CTEQ6L1 PDF) with the Perugia 2011C tune for parton showering, hadronisation, and underlying event.Background contributions from +jet and  →  production have been simulated using Sherpa 1.4.3 [7] with its default tune and the CT10 PDFs.
The MC samples were processed through a full ATLAS detector simulation [25], based on Geant4 [26], and reconstructed with the same software as that used for the data.All MC samples are corrected with data-driven factors to account for small differences in photon and lepton trigger, reconstruction, identification and isolation efficiencies between data and simulation.Additional  −  interactions (pile-up), occurring in the same and neighbouring bunch crossings are simulated and overlaid at the detector hit level on top of the hard-scattering process from the MC simulation.

Selection of 𝒁 → 𝒍 𝒍𝜸 events
A  →  candidate is formed by selecting the opposite-sign same-flavour dilepton pair with mass closest to   and the highest  T photon in the event.No explicit requirements are made on the presence or absence of other activity in the event, such as additional photons or leptons or jets.Only prompt leptons/photons are considered for the signal, meaning that they should not originate from decays of hadrons or -leptons.Background events from processes producing fake or non-prompt photons/leptons are suppressed by the lepton/photon identification criteria described below, including in particular isolation requirements on both leptons and on the photon.
Event candidates in both data and MC simulation are required to have fired at least one unprescaled single lepton or dilepton trigger, with a lowest  T threshold of 24 GeV for single leptons and of 12-13 GeV for dilepton triggers.The single lepton triggers provide a high trigger efficiency with respect to the offline selection, which requires at least one lepton with  T >25 GeV.The trigger efficiency for  →  events satisfying all the selection criteria described below is about 99%.

Photon and lepton selection
Photon and electron candidates are reconstructed [27] from clusters of energy deposits in the electromagnetic calorimeter combined with information from charged tracks reconstructed in the inner detector.Electrons and photons are identified by shower shape and hadronic leakage variables.Photons are required to satisfy all the requirements on shower-shape variables corresponding to the Tight photon identification criteria of Ref. [28].Electron candidates are required to satisfy the Loose likelihood requirement of Ref. [29].
Muon candidates are reconstructed [30] from combined tracks in the inner detector and muon spectrometer with a transverse momentum above 10 GeV over the pseudorapidity range || < 2.5.Over the range 2.5 < || < 2.7, stand-alone muon candidates in the muon spectrometer are used, based on tracks with reconstructed segments in three spectrometer stations.
Electron and muon candidates are required to originate from the appropriate primary vertex (each primary vertex candidate is reconstructed from at least two associated tracks with  T > 0.4 GeV).The significance of the transverse impact parameter, defined as the absolute value of the track transverse impact parameter,  0 , measured relative to the beam trajectory, divided by its uncertainty,   0 , must satisfy | 0 |/  0 < 5 for electrons and | 0 |/  0 < 3 for muons.The longitudinal impact parameter is required to satisfy | 0 | < 10 mm.
The photon, electron, and muon candidates are required to be isolated from other particles using selections based on calorimeter isolation, which is defined as the sum of the calorimeter cell transverse energies,   , inside an isolation cone of Δ = 0.2 around the lepton or photon [31].In the case of the calorimeter isolation for photons/electrons, the cells corresponding to the calorimeter cluster itself are excluded from the sum.In addition, the calorimeter isolation is corrected on an event-by-event basis for the energy deposited by the photon or lepton candidate, and, using the method described in Refs.[32][33][34], for the contribution from the underlying event and pile-up.Each photon/electron (muon) is required to have a normalised calorimeter isolation smaller than 0.3 (0.2). Finally, photons close to leptons are excluded if Δ  < 0.4.

Signal region definition
Candidate  →  signal events are selected by requiring that they contain at least one pair of oppositecharge, same-flavour leptons and at least one photon.The highest- T lepton is required to have  T > 25 GeV, while the second lepton must have  T > 10 GeV.The selected dilepton pair in the event is that which has an invariant mass,   , closest to   .Only 2% of the events have more than one dilepton pair.
To further improve the purity of the QED FSR sample, the invariant mass of the dilepton pair,   , is required to be between 20 and 80 GeV, while the three-body invariant mass,   , is required to be between 80 and 100 GeV.These selections mainly suppress the background from -boson+jet production with a jet misidentified as a photon, which can be clearly seen in Fig. 4 outside the signal region.For similar reasons, these selections also reduce considerably the background from QED ISR.The integrated fiducial cross-section results are also extracted for a stricter selection on the invariant mass of the dilepton,   > 45 GeV.The criteria described above were applied to the data, leading to 30571 (resp.34948) events selected in the  →  (resp. → ) channel (see Table 1 and Section 5 for more details).The reconstructed differential distributions are shown in Fig. 5 for the three observables of interest, the invariant mass of the positively charged lepton and the photon,   +  , the angular separation between the photon and the closer of the two leptons, Δ  , and the transverse momentum of the photon,    .

Background estimation
The dominant sources of background to the  →  signal originate from top-quark pair production and decay with two high- T prompt leptons and one prompt photon in the final state and from +jet production with two high- T prompt leptons and one jet misidentified as a photon in the final state.Other smaller sources of background are also considered, such as   production,  production with two leptonic decays of the -lepton, but also backgrounds with one fake or non-prompt lepton such as  +jet production and backgrounds from   production with the electron from the -boson decay misidentified as a photon.The sum of all sources contributing to events containing two prompt leptons and one prompt photon is labelled "prompt background" and is predicted by simulation, while the sum of all other sources is labelled "fake  background" and is estimated from data as described below.The respective contributions of these prompt and fake  background sources to the data are listed separately with their total uncertainties for the dielectron and dimuon case in Table 1.
Table 1: Total observed rates of data selected for the measurements in the  →  and  →  channels.Also shown are the expected contributions with their total uncertainties from the two main sources of background, the prompt background containing two leptons and one photon, predominantly from top-quark pair production, and the fake  background containing two leptons and one jet misidentified as a photon, predominantly from +jet production (see text).The last row shows the expected signal with its total uncertainty (see Section 8.2).The overall background contribution to the signal is expected to be small, as can be seen from the small tails outside the pole region in the three-body invariant mass,   , distribution of Fig. 4. A simple extrapolation from the events with large   to the signal region leads to an expectation at the percent level.The nominal method used for the background estimation in this analysis is an unbinned maximum likelihood fit of the   distribution over a wide range, 60 <   < 120 GeV, using simulation for the prompt background sources (shape and normalisation) and also for the shape of the fake  background.This fit is performed for each bin of the differential cross-section measurements.

Channel
A cross-check of the fake  background estimate is performed with a data-driven two-dimensional sideband method using an inverted isolation selection for the photon candidate on one side and requiring that  background estimate.As shown in Table 1, the overall background is at the level of 2-3% of the expected signal with a total relative uncertainty of 13%.

Correction for detector effects
The reconstructed differential distributions presented in Fig. 5 are corrected for detector effects and bin-tobin migrations using an iterative Bayesian unfolding method applied to events that pass the detector-level selections.The method is applied with three iterations implemented in the RooUnfold framework [35].
The resulting unfolded distributions for the  →  and  →  signal are obtained at both bare-lepton and dressed-lepton level.
The particle-level selections shown in Table 2 are chosen to be as close as possible to the selections described in Section 4. The leptons and the photon are required to be prompt, meaning that they should not originate from decays of hadrons or -leptons.The small difference in pseudorapidity ranges for electrons and muons is retained at particle level since the electron and muon channels are not combined in this analysis (except for the low statistics  →  sample described in Section 9).
The response matrices, which connect the distributions at reconstruction and particle level, as well as the efficiency correction factors are derived using the Powheg+Pythia+Photos signal MC sample.The reconstruction, trigger and isolation efficiencies as well as the photon/lepton momentum scale and resolution in the MC simulation are corrected to match those determined in data.The transverse momentum distribution of the simulated  →  events was reweighted to that observed in data to minimise the only potentially significant bias from QCD modelling in the simulation (see Section 7 for more details).Since the charged leptons and the photon are precisely measured compared to the bin sizes chosen for the differential distributions, the systematic uncertainties related to the unfolding method were found to be below 0.1%.

Systematic uncertainties
This section describes the systematic uncertainties in the measurements of the normalised differential fiducial cross-sections for the observables presented in Fig. 5.These systematic uncertainties are grouped according to their source, and their typical relative values over most of the kinematic range of the measured observables are listed in Table 3, while the relative variations of the overall systematic uncertainty as a function of the measured observables are shown in Fig. 6, together with the variations of the statistical uncertainty in the data and of the total uncertainty in the measurements.The statistical uncertainty in the data is the dominant uncertainty over most of the measurement bins of all three observables.
The experimental systematic uncertainties in measured observables involving only leptons and photons are expected to be small, and particularly so for normalised distributions for which only the shape variations in these uncertainties remain.This can indeed be seen from Table 3 which shows contributions of at most a few per mille from the lepton/photon experimental systematic uncertainty sources.However, Fig. 6 shows that the overall systematic uncertainty is the dominant contribution at the edge of the phase space for the   +  and Δ  distributions.This arises from the very limited number of simulated events in these regions which inflates the overall systematic uncertainty even though it is statistical in nature.Over most of the phase space of interest, one expects that the shapes of the electron and muon channels are identical, so the  −  difference has been used as an additional cross-check of the experimental systematic uncertainties.The weighted mean of the  −  data/MC double ratio difference was found to be −0.0011± 0.0085 (stat) ± 0.00007 (syst).The systematic uncertainty is reduced here purely to the contribution of the experimental systematic uncertainties which are uncorrelated between electrons and muons.The  −  difference was thus found to be consistent with its statistical uncertainty, and no systematic trend was observed as a function of the measured differential distributions.
Table 3 shows that the theoretical uncertainties from QCD scale variations and PDFs, related to the acceptance of the photon/lepton fiducial cuts, are very small, at the per mille level, as expected.However, any QCD-related mis-modelling of -boson production not covered by these uncertainties may have a significant and direct impact on the QED-related observables of interest.As mentioned above, the simulated events were reweighted to reproduce the observed transverse momentum spectrum of  →  decays.This has potentially a direct impact on the photon transverse momentum spectrum, but it also affects at the percent level certain regions of phase space for the   +  and Δ  distributions (the most significant impact has been observed in the case of the Sherpa 1.4 predictions).After this reweighting procedure, the residual QCD modelling systematic uncertainties were found to be very small, at the few per mille level.This was further verified by also reweighting to data the other kinematic observables related to -boson production, namely the -boson rapidity and the angular coefficients [36] related to the lepton kinematics in the rest-frame of the -boson.As a final cross-check, the Sherpa 1.4 predictions for the reconstructed observables, after reweighting to the observed transverse momentum spectrum of  →  decays, were unfolded using the nominal Powheg+Pythia+Photos MC.The resulting systematic uncertainty, as derived from the difference between the generator-level Sherpa 1.4 prediction and that obtained after this unfolding cross-check, was found to be negligible, i.e. below the per mille level.

Results for 𝒁 → 𝒍 𝒍𝜸 process 8.1 Differential cross-sections
The measurement results are presented here as unfolded normalised differential cross-sections for the three observables of interest, namely the invariant mass of the positively charged lepton and the photon,   +  , the angular distance, Δ  , between the photon and the closer of the two leptons, and the transverse momentum,    , of the photon.These differential cross-sections are normalised to the total measured fiducial cross-section to minimise the experimental uncertainties and also to remove potentially dominant QCD-related theoretical uncertainties when comparing to predictions.The results below are shown separately for the electron and muon channels in the case of bare leptons.This maximises the sensitivity to expected small differences in the distributions between the two lepton flavours (see Figure 3).
Figure 7 shows in the case of bare leptons, separately for the electron and muon channels, the unfolded normalised differential distributions for the data with their uncertainties compared to the predictions with uncertainties from Powheg+Pythia8+Photos, Sherpa 2.2.4 (YFS), and KKMChh.As discussed in Section 3, Powheg+Pythia8+Photos considers only QED FSR and has an associated theoretical uncertainty of 0.2% obtained by varying the configuration of Photos from its default one to the slightly less precise prediction not including the -boson matrix-element correction discussed in Section 3.This small theoretical uncertainty assumes that the contributions from QED ISR and IFI are even smaller.The combined QED ISR/IFI contribution in the fiducial region of the analysis presented in this paper has been estimated with KKMChh to be 3 × 10 −3 for the integrated fiducial cross-sections, with an uncertainty below 10 −3 .These results have been confirmed using the CompHEP [37] code.In the case of the Sherpa 2.2.4 predictions which are also shown in Fig. 7, the uncertainty band is much larger and probably over-conservative, about 2%.This was estimated using a recipe similar to that described in [38], corresponding to higher-order QED corrections which are only known at the -boson pole and are not applied outside a very narrow region around   to assess this uncertainty.
The agreement between the measurements and the predictions from the MC simulations is reasonable for all three distributions over most of their range, both for electrons and muons.The agreement between the predictions themselves is better for   +  than their respective agreement with the data, suggesting perhaps that most of the fluctuations seen in the ratio plots for this observable are related to statistical fluctuations in the data since they are limited to a few bins.However, there is a noticeable discrepancy between the measurement in the electron channel and the Powheg+Pythia8+Photos prediction for values of Δ  larger than 2.5, at the edge of the phase space where the two leptons from -boson decay recoil against the photon. [GeV] ) )  The bottom panels represent the ratios of these predictions to the data, where the gray band around unity represents the total data uncertainties and the dashed lines around the predictions represent their respective total uncertainties (see text).

Integrated fiducial cross-sections
The normalised differential fiducial cross-section measurements presented above are only sensitive to the shape of the distributions, but they benefit from cancellations of experimental systematic uncertainties for the sources that do not vary much with lepton  T , such as the lepton efficiencies.This section presents the integrated fiducial cross-section measurements in the electron and muon channels and discusses the results for three different ranges of the dilepton mass.The uncertainties in the integrated fiducial cross-sections are summarised in Table 4.They include the uncertainty of 1.9% in the integrated luminosity and are somewhat larger for some of the experimental sources than are those listed in more detail in Table 3 for the reasons explained above.
The comparisons between the measured integrated cross-sections over three different dilepton mass ranges and the corresponding theory predictions, normalised to the QCD NNLO prediction used in [39] for the total inclusive cross-section for -boson production, are presented in Table 5.The agreement between the measurements and the predictions is reasonable within their comparable overall uncertainties.The uncertainties in the theory predictions are also taken from [39].As explained in Section 1 and shown in Table 5 for the theory predictions, the cross-sections for the  →  process are expected to be larger than those for the  →  process.This difference is enhanced significantly in this measurement by

First measurement of the 𝒁 → 𝒍 𝒍𝜸𝜸 process
This section presents a first measurement of the  →  process at the LHC (Fig. 8).Given the limited statistics expected for the dataset considered here, the candidate event samples for each charged lepton flavour were combined together, and the selection requirements on the sub-leading photon were loosened with respect to the leading one.The sample considered was the same one as for the  →  analysis in terms of the charged lepton requirements, with however a tighter selection requirement of  T > 15 GeV for the softer of the two leptons, but without any requirement on the dilepton mass.The two photons were required to be within the same pseudorapidity range as for the  →  analysis with thresholds in  T set to respectively 15 and 10 GeV and within an angular distance, Δ  of at least 0.4 from each charged lepton.In addition, the angular distance between the two photons, Δ  , was also required to be larger than 0.4.Figure 9 (top left) shows the invariant mass of the two leptons and two photons,   , versus the invariant mass of the two leptons and the highest  T photon,   .As illustrated in this figure, for values of   larger than 80 GeV, significant contributions of background from  →  +jet(s) events and from  →  +jets events are expected.
In order to reduce these contributions as well as possible contributions from  →  events with initial state QED radiation to the percent level or less, the invariant mass,   , was required to be below 80 GeV for both photons.These requirements are summarised as fiducial volume selections in Table 6.A total of 116 events, corresponding to 61 electron pairs and 55 muon pairs, were thus selected as  →  candidates.The dominant background from  → +jet events was estimated by applying the same two-dimensional sideband method as that described for the main analysis: the quality and isolation criteria were inverted for the second photon.The background was found to be less than one event and Fig. 9 (top right) illustrates this through the observed distribution of   which displays no events outside the pole region.
The data were compared to MC  →  samples produced using the same Powheg+Pythia8+Photos version used for the main analysis, Sherpa 1.4 ME  (see Section 3), and KKMChh.Figure 9 (top right, middle and bottom) shows that the predictions describe the data within their large statistical uncertainties for a selection of observables.
Finally, the data were corrected for detector efficiencies and resolution, and fiducial cross-sections were extracted together with unfolded distributions of several observables.The fiducial cross-section is measured to be 22.2 ± 2.1 (stat) ± 1.2 (exp) ± 4.1 (theo) fb for bare leptons.In contrast to the  →  case, the expected efficiencies from the Sherpa 1.4 ME  and Powheg+Pythia8+Photos MC samples are significantly different and this difference was included as an additional theoretical systematic uncertainty in the measurement result for the  →  cross-section quoted here.
The expected fiducial cross-section from Sherpa 1.4 ME  is 9.9±0.1 (stat) fb, that from Powheg+Pythia8+Photos is 9.4 ± 0.3 (stat) fb.Both values are significantly smaller than the measured one.However, the expected fiducial cross-section from Sherpa 2.2.4 is 20.9 ± 0.2 (stat) fb, in agreement with the data.Although the purely QED systematic uncertainty in the predictions cannot explain the observed discrepancies between the measurement and two of the predictions, the overall systematic uncertainty in the predictions, which should include QCD and mixed QCD/QED effects, is not known.This uncertainty is clearly of QCD origin, as shown by the two versions of Sherpa.Figure 10 shows the unfolded distributions for the data compared to the predictions.The uncertainties shown are the combined statistical and systematic uncertainties for the data and only the statistical uncertainties for the predictions.The agreement in shape between data and predictions is reasonable within the uncertainties of this first measurement.

Conclusions
This paper presents the first ATLAS measurement of unfolded normalised differential fiducial cross-sections for  →  and  →  final states together with a few integrated fiducial cross-sections.The overall accuracy of the differential measurements is at the 1-2% level over a large fraction of the overall phase space, with an average overall systematic uncertainty of 0.6% for the chosen observables,   +  , Δ  , and    .The results are in agreement with predictions of MC generators containing state-of-the-art QED FSR calculations such as Photos (in PowHeg+Pythia8+Photos) and YFS (in Sherpa 2.2.4).They are also in agreement with the recent calculations of KKMChh which provides a consistent and gauge-invariant breakdown of the different contributions to final states containing two charged leptons and a high- T photon, and thus demonstrates that the QED ISR and IFI contributions are very small for the fiducial selections of the analysis described in this paper.
Finally, the  →  process has been observed for the first time.The differential distributions for this process are described by the available calculations within present large statistical uncertainties of the measurement.The  →  fiducial cross-section is consistent with Sherpa 2.2.4 and is twice as large as predicted by PowHeg+Pythia8+Photos and Sherpa 1.4 ME  .

Figure 4 :Figure 5 :
Figure4: For  →  (left) and  →  (right) decay candidates in data, distribution of the three-body invariant mass,   , versus the dilepton mass,   .The dashed lines indicate the selection criteria applied in the analysis to these two reconstructed invariant masses.These criteria minimise in particular the background from  bosons produced in association with jets in which one jet is misidentified as a photon.Also shown are the projections of these figures onto the vertical (middle plot with only the selection on   applied) and horizontal (bottom plot with only the selection on   applied) axes with the corresponding vertical lines indicating the selection cuts described above.The curves show the signal predictions from the nominal Powheg+Pythia8+Photos MC normalised to the data.

Figure 6 :
Figure 6: Breakdown of relative uncertainties for the normalised unfolded differential cross-sections as a function of the observables of interest,   +  (top), Δ  (middle), and    (bottom), shown separately for the  →  (left) and  →  (right) channels.

Figure 7 :
Figure 7: Unfolded normalised differential fiducial cross-sections for the  →  (left) and  →  (right) channels.The results are shown in the case of bare leptons for the   +  (top), Δ  (middle), and    (bottom) distributions.The data are compared to the predictions from Powheg+Pythia8+Photos, Sherpa 2.2.4,and KKMChh.The bottom panels represent the ratios of these predictions to the data, where the gray band around unity represents the total data uncertainties and the dashed lines around the predictions represent their respective total uncertainties (see text).

Figure 9 :
Figure9: For all selected  →  candidates before applying any selection on the invariant mass of the two leptons and each photon described in the text, distribution of invariant mass of the two leptons and the two photons,   versus the invariant mass of the two leptons and of the highest  T photon,   (top left).The horizontal dashed (blue) rectangle indicates the region where the  →  signal is expected and the vertical dashed (green) rectangle indicates the region dominated by background from  →  events accompanied by one or more jets.Also shown after all selection requirements are the distributions of the total invariant mass (top right), of the invariant mass of the two leptons and the second photon, of the dilepton mass, of the angular distance between the two photons, and of the transverse momentum of the second photon.The data are compared to predictions from Powheg+Pythia8+Photos and Sherpa 1.4 ME  , normalised both to the data.The results are shown for bare leptons, and electrons and muons have been combined.The uncertainties shown are statistical.

Figure 10 :
Figure 10: For all selected  →  candidates, unfolded normalised distributions of   +  (top left), of Δ  (top right), of the angular distance between the first (middle left) and second (middle right) photon and the closest lepton, and of the transverse momentum spectrum of the two photons (bottom).The unfolded data are compared to predictions from Powheg+Pythia8+Photos, Sherpa 1.4 ME  , and KKMChh.The results are shown for bare leptons, and electrons and muons have been combined.The gray band represents the total data uncertainties and the dashed lines around the predictions represent their statistical uncertainties, which are dominant, for the predictions.

Table 3 :
Breakdown of systematic uncertainties for the normalised differential fiducial cross-sections.The values shown are typical over most of the kinematic range of the unfolded observables.

Table 4 :
Uncertainties in the integrated fiducial cross-section measurements.The values for bare leptons with   > 20 GeV are shown.

Table 5 :
Integrated fiducial cross-sections (pb) for three ranges of the dilepton mass   .The measurements are presented for bare leptons with the breakdown of their uncertainties together with the predictions from Powheg+Pythia8+Photos normalised as described in the text.