Observable phi-eta-star at LHC and second-order QED matrix element in Z/gamma to l+ l- decays

In a recent publication by ATLAS collaboration a new observable, the so-called phi-star angle, was used for precise measurement of transverse Z momentum. One of the dominant systematic errors for this measurement originates from the theoretical control of QED final-state bremsstrahlung. At present, it is estimated at the 0.3% level for the shape of the phi-star distribution. In this paper we discuss the possible effects of the second-order QED matrix element for that quantity. For that purpose, results from simulations based on the Yennie--Frautchi--Suura (YFS) exponentiation and featuring the second-order matrix elements are used and compared with the case when the matrix element is restricted to the first order. From this study we conclude that in order to reach the precision below 0.3% for the phi-star distribution at the LHC, inclusion of the second-order QED matrix element in a respective Monte Carlo event generator is necessary.


Introduction
The main purpose of LHC experiments is to search for new elementary particles and interactions. The Higgs-like particle observability papers [1,2] document a fundamental achievement for this goal. Another important class of the LHC measurements aims on high-precision consistency tests of the Standard Model. In this respect, a precise measurement of the W mass plays a particularly important role [3]. To reduce systematic errors, improvements in the measurement techniques are desirable. It is generally believed that higher precision can be achieved thanks to the use, whenever possible, of leptonic degrees of freedom rather than hadronic ones. In this way the precision better by even an order of magnitude can be achieved. At present, the best measurements of W mass have been performed in pp collisions at Tevatron [4,5,6]. In the CDF collaboration the precision of 19 MeV on the W mass was achieved. They have estimated that the largest contribution from theory to the systematic error originates from parton distribution functions (PDF's) and initial-state hadronic interactions in general.
An approach [7,8,9] based on measuring the φ * η angle, instead of Z transverse momentum (p Z T ), may offer a significant improvement. In Ref. [10] the first results of φ * η measurements with ATLAS experiment were presented already. In that experimental publication the systematic error of 0.3% to φ * η due to implementation of QED final state radiation (FSR) in the Monte Carlo generators used, was estimated in proportion of differences observed between PHOTOS [11,12,13] and SHERPA [14] predictions. If one understands the pattern of these differences, one can hopefully reduce systematic errors further. In this sense presented here work is continuation of recent efforts in this direction documented in Ref. [15].
One of the possible sources of the above differences may be the second-order matrix element for QED final-state radiation (FSR) in Z decays. Such a matrix element is missing in PHOTOS and in SHERPA, but its dominant contribution is taken into account thanks to iterative and/or multiphoton nature of the algorithms. Multitude of tests for PHOTOS were devoted to this point over time [12,15,16]. In particular, it was shown [11] that despite the fact that PHOTOS is not using a matrix-element-based kernel, it agrees substantially better with distributions obtained from the KKMC Monte Carlo [17] if exclusive exponentiation featuring the second-order matrix element is used rather than if the matrix element is restricted to the first-order only. The KKMC program is particularly suitable as a source of precise numerical benchmarks, because its precision was studied at the LEP time in a great detail, see e.g. Ref. [18], and with precision requirements surpassing the present-day requirements of the LHC experiments. It is known since more than 20 years now [19] that second-order terms are necessary for an exponentiation-based Monte Carlo generator if it is expected to assure a high precision in the case of precision observables of e + e − colliders. That is why such terms were implemented in the e + e − → l + l − generators KORALZ [20] and in the KKMC program as well. The precision of 0.06% was evaluated in [21] for small-angle Bhabha scattering simulated with the Monte Carlo generator BHLUMI [22].
For the LHC applications, the second-order QED FSR matrix element was found to be important for simulations of background for Higgs-boson searches in the 2γ channel [16]. With time and improving experimental precision, contribution from the second-order QED FSR matrix element may become important for more inclusive observables at the LHC as well. First attempts to implement the second-order matrix element (ME) to the LHC Monte Carlo generators are on-going; a good example is [23].
Before such projects are completed, let us provide some quantitative results using the LEP-era Monte Carlo program KKMC. For emulation of the initial-state hadronic interactions, parton distribution functions (PDF's) combined with the transverse momentum smearing based on the RESBOS program [24] are used only. Such an approach doesn't take into account the detector simulations, nonetheless it should be sufficient to quantify the size of the effect of the second-order ME.
Our note is organized as follows. In Section 2 we describe our Monte Carlo set-up. Section 3 is devoted to a careful definition of the φ * η distribution, including selections criteria, etc. In Section 4 we present our numerical results. Section 5 summarizes the paper with conclusions.

Monte Carlo set-up
The KKMC event generator is constructed for e + e − → µ + µ − (τ + τ − , qq) processes and for centre-of-mass system energies from the production threshold up to 200 GeV. For our purposes it was adapted to work for quark-anti-quark annihilation process qq → l + l − . However, in this case, only QED FSR can be generated. By default, the second-order matrix element is active, but the program can be downgraded to the first-order FSR as well. We use this program as a building block in our simulation chain. With the help of WINHAC, the MC event generator for Drell-Yan processes [25,26], using the MSTW2008NLO [27] parametrisation of PDF's we generate a series of the Bjorken variables x 1 and x 2 for pp collisions at √ s = 7 TeV corresponding to production of Z-bosons in a very narrow invariant mass range: M Z ± 1 MeV. Then we generate monochromatic Z decays including QED FSR bremsstrahlung. Let us stress that we aim at the estimation of a small secondorder QED FSR effect and not at the complete predictions. That is why we can limit ourselves to the dominant QCD initial-state effects only.
The generated sample has to be convoluted with the initial-state QCD effects. For the sake of convolution, each generated qq → Z → e + e − n(γ) event is boosted to the laboratory frame using x 1 and x 2 of the incoming quarks generated by WINHAC Monte Carlo [25,26]. To emulate transverse momentum of Z we use smearing on the basis of a p Z T distribution (histogram) obtained from the RESBOS program [24], which was found to model data within the 4% accuracy, as concluded in Ref. [10]. From this distribution we randomly generate p Z T of the Z boson, to be used in the boost to the laboratory frame. This is a crude approximation which, again, is acceptable in evaluation of the considered effect, only because it is expected to be small. In Fig. 1, distributions of p Z T and y Z (for uū annihilation) are reproduced for the reference.
Two versions of the results will be compared: the one when the exclusive exponentiation, as embedded in KKMC, is featuring the second-order matrix element (option CEEX2) and the another one when it is limited to the first-order only (option CEEX1), as defined in Ref. [18]. The technical advantage of correlated samples will help to control differences for the two options in less populated regions of the phase space. The CEEX1 variant will be implemented with an appropriate weight. This is advantageous, as we expect the difference of CEEX2 and CEEX1 to be rather small, but at the same time limitations due to restricted treatment of the initial-state hadronic interactions are present. We have to limit ourselves (for the particular run) to a fixed flavour of quarks entering the hard process, moreover the virtuality of the intermediate state has to be fixed as well. For the presented plots the incoming u-quarks were chosen and the intermediate-state virtuality was fixed in KKMC generation to Z boson mass. The results for the incoming d-quarks essentially coincide with the presented ones. 3 Definition of observable φ * η As proposed and studied in Refs. [7,8,9], an alternative observable to study the Z transverse momentum (p Z T ) is φ * η , defined [8] as: where φ acop ≡ π − ∆φ, ∆φ is the azimuthal opening angle between the two leptons, and the angle θ * η is a measure of the scattering angle of the leptons with respect to the proton beam direction in the rest frame of the dilepton system. The angle θ * η is defined [8] by where η − and η + are the pseudorapidities of the negatively and positively charged lepton, respectively.
The observable φ * η is expected to be less sensitive to experimental resolution and it can probe largely the same physics as p Z T for small p Z T or φ * η . The theoretical calculations for the φ * η are documented in Refs. [28,29]. The first experimental measurement of the φ * η by the D0 Collaboration [30]. demonstrated that the order of magnitude improvements in experimental precision could be achieved with the φ * η technique. The φ * η was then employed by the ATLAS and LHCb experiments in recent publications [10,31]. The first measurement of the normalized φ * η distribution at √ s = 7 TeV pp collisions performed by the ATLAS experiment is very likely one of the most precise measurement at the LHC, with the total uncertainty at the level of 0.5-0.8%. The normalized differential cross section was defined [10] as 1/σ fid · dσ fid /dφ * η , where σ fid is measured within the fiducial lepton acceptance defined by the lepton (ℓ = e, µ) transverse momentum p ℓ T > 20 GeV, the lepton pseudorapidity |η ℓ | < 2.4 and the invariant mass of the lepton pair: 66 GeV < m ℓℓ < 116 GeV. In [10] three lepton definitions with respect to the QED FSR were used. A "bare" level defined directly by the final-state leptons after QED FSR will be followed in our paper. A "dressed" level is defined if lepton four-momenta are recombined with photons radiated within a cone of ∆R = (∆η) 2 + (∆φ) 2 = 0.1. Finally, a "Born" level is defined when four-momenta of the final-state leptons before Monte Carlo simulation of QED FSR are used. The measured φ * η distribution in data after correcting only for all detector effects is at bare level. To account for QED FSR effects and to obtain the φ * η distribution at Born and dressed level a correction, computed using PHOTOS was applied to the data.
The dominant systematic uncertainty of this measurement is due to φ * η -dependent modelling of QED FSR in this correction, assigned to be 0.3% by comparing predictions from PHOTOS (interfaced in POWHEG+PYTHIA6) [32,33] and SHERPA, as studied in Ref. [34]. This uncertainty is extracted by looking at the ratio between the normalized φ * η distribution predictions at the dressed or bare levels and the one at the Born level in PHOTOS and in SHERPA.
The uncertainties on this correction are extracted by looking at the ratio between the normalized φ * η distribution at the dressed or bare levels and the one at the Born level in PHOTOS and in SHERPA (in all cases no detector smearing effects are introduced).
The ratio of the predictions of PHOTOS and SHERPA in the case of the electron channel is shown 1 in Fig. 6.5 of Ref. [34]. Statistical fluctuations are too large to draw firm conclusions, nevertheless a hint of difference (wavy structure) seems present which appear to be more profound when comparing the two programs for the bare case. One can only partly profit from the statistical correlation of the left and right hand side plots of Fig. 6.5 which are constructed from the same sample of events.
We can investigate this difference by means of comparing the results from KKMC: at the first-order mode and at the second order mode. This comparison is expected to be sensitive to a pure QED effect. The difference as a function of φ * η can be used as an correction or as estimation of the systematic error for the missing second-order matrix element of QED. We present details of the comparison in the following section.

Numerical results
We combine histograms and plot superimposed distributions for φ * η in the case when the first and second-order QED FSR matrix elements are used in KKMC. KKMC is not suitable for simulations of events in pp collisions. We overcome this obstacle by switching off all the QED initial-state radiation in KKMC, the final-state radiation is then retained only. We replace the incoming electrons with quarks. Generation of qq → l + l − nγ events at the fixed centre-of-mass energy can be then performed. We use several fixed quark-pair virtualities close to the Z peak 2 .
For the KKMC-based plots the following cuts were used. These cuts are similar but not identical with the one used in experimental studies described earlier. We request p l ± T > 20 GeV, |η l ± | < 2.4, in our case the region 1.37 < |η l ± | < 1.52 was rejected. Finally, acoplanarity for leptons was requested to be smaller than 3. We have used bare electrons only. We have investigated several options for the cut on lepton pair mass m ee and not only the default one used by the ATLAS experiment: 66 GeV < m ee < 116 GeV. The cut removes hard-photon configurations and reduces higher-order QED effects, as can be seen in Fig. 3: from 0.2-0.3% to below 0.1%. The sample used for construction of these plots was with the Z virtuality equal to M Z , but we have checked the contributions from higher virtualities of an intermediate state as well. Then, see Fig. 3, the correction was higher, even if the m ee cut was present -it approached 2%.
In Fig. 2 we show the results with realistic cuts for φ * η defined in the previous section, where the effects of the longitudinal momentum of Z due to PDF's are taken into account. The transverse momentum, p Z T , of Z is, however, set to zero. As a consequence, φ * η is non-zero only because of QED FSR. At least one photon is needed for events to populate other bins than the first one. That is why the effects of the second-order ME corrections are numerically large. They represent leading correction: nearly all events are placed in the first bin of the histogram, corresponding to φ * η = 0 or below the detector resolution, whereas contents of all other bins are at least two order of magnitude smaller. The secondorder matrix element represent the correction which grows from 1% to 15% at the end of the spectrum. As expected, the CEEX1 prediction is larger for large φ * η . The β 1 correction of the YFS exponentiation is negative for hard photons [18].
For Fig. 3 no restrictions on p Z T spectrum were introduced and one can see that the correction is small over the whole range of φ * η , nonetheless grows up to 0.2% at the end of the spectrum. This is consistent with what we observe in Fig. 6.5 of Ref. [34] at the end of the spectrum for bare leptons. The distributions of the Z transverse momentum p Z T and its rapidity y Z without any cuts, used in our MC simulations, are shown in Fig. 1.
It is technically simple to prepare such kind of plots for the comparisons of PHOTOS and KKMC, but similar results have recently been presented in Ref. [15]. Also the comparisons with other variants of the matrix element implementation in KKMC, such as EEX1, EEX2, EEX3, may be instructive for detailed study aiming at the precision implementation of the exclusive exponentiation algorithms into pp collisions Monte Carlo generators. This is, however, out of the scope of this work.  Figure 2: The φ * η distribution: the comparisons of the CEEX2 and CEEX1 results. The sample of uū → e + e − (nγ) events is generated. The virtuality of the Z boson equal to its mass is used. The boost to the laboratory frame of Z is performed prior to histogramming with cuts. The longitudinal momentum of Z is generated according to WINHAC in pp collision at √ s = 7 TeV, while p Z T = 0. In the first bin, the configurations with φ * η smaller than 0.001 are collected.
Let us now comment on the reliability of presented results. In left hand-side plot of Fig. 6.5 from Ref. [34], the corrections for the φ * η calculated with dressed electrons, are compared for the SHERPA and PHOTOS cases; they coincide down to about 0.1%. The distinct treatment of hadronic initial-state interactions in PYTHIA and SHERPA can not therefore be the prime source of the wavy structure in φ * η corrections calculated for the similar plot but with bare electrons, see right-hand side of Fig. 6.5. Note that the two plots of Fig. 6.5 are in part statistically correlated as they are constructed from the same sample of events.
The difference between dressed and bare cases likely originates from the different treatment of collinear photons. We may deduce from Fig. 6.5 and our Fig. 3 that absent in SHERPA second-order QED FSR matrix element can in part be responsible for the observed hint of the difference in corrections calculated for bare electrons with the two programs. In PHOTOS second-, and higher-order QED FSR terms, important in collinear limits, are taken into account thanks to iterative algorithm.  Figure 3: The φ * η distribution: the comparisons of the CEEX2 and CEEX1 results. The sample of uū → e + e − (nγ) events is generated. The virtuality of the Z boson equal to its mass is used. The longitudinal momentum of Z is generated according to WINHAC in pp collision at √ s = 7 TeV, while its transverse momentum is generated according to the prescription given in the text. The boost to the laboratory frame of Z is performed prior to histogramming with cuts. In the first bin, the configurations with φ * η smaller than 0.001 are collected. In the upper row of plots, cuts as explained in the text except cut on m ee , are used. The plot with 66 GeV < m ee < 116 GeV is on the bottom-left side. On the bottomright side, where in addition to this cut, exceptionally the events with the Z/γ * virtuality of 115 GeV are taken, a wavy structure of the correction appears. For the virtuality above 116 GeV the correction becomes larger, 2% or more, because the no-bremsstrahlung events do not contribute. This demonstrates the limitations of our approach. Contributions from all Z/γ * virtualities have to be combined for realistic prediction, rather than of fixed ones, we have used for the estimate of the effect size.

Summary
In this paper we have studied the effect of the second-order matrix element with respect to the first-order one. We have concentrated mainly on the case when bare electrons were used, because it was contributing to the discussion of the sytematic error in measurement of Ref. [10]. The Monte Carlo simulation based on the exponentiation was used for the QED final-state radiation. The φ * η distribution in Z/γ → e + e − n(γ) production and decay at LHC was studied. We have started with the case when p T of the decaying Z was set to zero, see Fig. 2, the effect of the introduction of the second-order matrix element was large, at the level of 1-15%. However, the bulk of the distribution resided in the first bin of the histogram of φ * η < 0.001. When the initial-state Z-boson p Z T is taken into account, then the φ * η distribution is representing convolution of the dominant effect of initial-state quark transverse momenta and the QED FSR correction. The effect of p Z T of the incoming Z dominates. The corrections due to the second order matrix element are smaller in the case of dressed leptons. This is consistent if the effect is due to consecutive collinear photons. Already the second-hardest photon requires the CEEX2 matrix element to be present in the exponentiation-based MC simulations. The ratio of the parameters defining the range of the lepton-pair invariant mass m ll to the virtuality of the intermediate Z/γ * used in our generation seems to be decisive for the details of the shape of the CEEX2-CEEX1 differences.
We can conclude that starting from the precision better than 0.3% for the φ * η distribution, the dominant contribution of the second-order QED FSR matrix element has to be included in discussion of systematic errors for Monte Carlo simulations used for the LHC experiments.
The above effects are slightly larger than presently aimed precision for the fermion-pair emission or interference, as discussed in Ref. [15] for example. The φ * η is an example of relatively inclusive observable, where at the least, dominant parts of the second-order QED FSR effects have to be taken into account, even if the exclusive exponentiation is in use.
We could observe that the effect of the second-order matrix element was more significant in the case of the φ * η observable defined with the bare leptons than with the dressed ones, again confirming that the dominant effect is of the leading-logarithmic nature. That is why, one may expect that for many applications, the systematic error concluded in Ref. [10] is an overestimation. To understand when it is actually the case, more detailed studies of cut-off dependence than presented in this paper are necessary.