QED challenges at FCC-ee precision measurements

The expected experimental precision of the rates and asymmetries in the Future Circular Collider with electron positron beams (FCC-ee) in the centre of the mass energy range 88-365GeV considered for construction in CERN, will be better by a factor 5-200. This will be thanks to very high luminosity, factor up to $10^5$ higher than in the past LEP experiments. This poses the extraordinary challenge of improving the precision of the Standard Model predictions by a comparable factor. In particular the perturbative calculations of the trivial QED effects, which have to be removed from the experimental data, are considered to be a major challenge for almost all quantities to be measured at FCC-ee. The task of this paper is to summarize on the"state of the art"in this class of the calculations left from the LEP era and to examine what is to be done to match the precision of the FCC-ee experiments -- what kind of technical advancements are necessary. The above analysis will be done for most important observables of the FCC-ee like the total cross sections near $Z$ and $WW$ threshold, charge asymmetries, the invisible width of $Z$ boson, the spin asymmetry from $\tau$ lepton decay and the luminosity measurement.


Introduction
The high-energy high-luminosity future circular electron-positron collider FCC-ee [1][2][3], considered for construction at CERN, would feature luminosity up to a factor 10 5 higher than at LEP collider. Together with the improvements of the detector techniques it would allow to reduce an experimental error typically by a factor of 5-50, and in some cases even by factor 200. The uncertainty of the corrections due to trivial but large QED effects to be removed from the data would then become for many observables a dominant one and its reduction becomes a major challenge at the FCC-ee programme. For instance, some QED effects of order 0.1%, which at LEP could be neglected and accounted for in the error budget, will have to be calculated with two-digit precision, removed from data and their uncertainty hopefully below 0.001% will enter into a combined experimental and theoretical error. It is therefore important to review already now the present state of the art in the precise calculations of the SM for e + e − annihilation processes left over from the LEP era not only for multiloop pure electroweak corrections [4][5][6], but also for QED effects and to evaluate the prospects of the necessary future improvements in this area of the theoretical physics.
The modern techniques of calculating QED corrections in electron-positron colliders for arbitrary experimental cut-offs using Monte Carlo event generators were founded in 1980's for in PETRA and PEP experiments. Typically, these calculations were implementing complete O(α 1 ). Their precision tag (neglected higher orders) was of order 1-2%, see for instance ref. [7]. Analytic phase space integration of QED distributions was not playing at that time any major role, apart from being used in designing Monte Carlo algorithms. The LEP era has seen development of the entire new range of QED calculations for many new processes, mostly in the form of the MC event generators. However, for the analysing of data near the Z resonance, numerical programs based on analytic integration over the phase space were also playing an important role. Also first QED calculations at the O(α 2 ) have appeared [8], but the biggest boost in the precision came from resummation to infinite order of soft photon corrections [9] and from the calculations of the collinear logs due to small lepton masses to higher orders, up to O(α 3 ) [10]. Several of the MC programs were also armed with the complete O(α 1 ) electroweak (EW) corrections. The example of the most sophisticated MC event generator from the LEP era with soft photon resummation, complete O(α 2 ) QED and O(α 1 ) EW corrections is KKMC program [11]. Generally, in LEP data analysis near the Z peak, the use of MC generators was often limited to removing detector inefficiencies and partly removing experimental cut-off effects. Fitting EW parameters like masses of the top quark and the Higgs boson was done with non-MC programs like ZFITTER [12] and TOPAZ0 [13]. However, at the LEP2 above the W W threshold, MC event generators [14,15] were the only tools capable to calculate QED+EW Standard Model predictions for the total cross section and distributions of the e + e − → W + W − process and also were used to extract (fit) the mass of the W boson from data.
As it will be argued in the following, due to significant increase of the experimental and statistical precision at FCC-ee, the role of the MC event generators will increase. Near Z resonance a combination of MC programs and semianalytical non-MC calculations will be used in the data analysis [4] 1 . This will be due to many reasons, the rise of non-factorisable QED corrections above the level of the experimental precision, the need of further development of the resummation techniques of soft and collinear corrections which will be the only available within MC programs and the generally bigger role of the multiparticle final states, for which the MC technique is the only method of integration over the phase space. The future tools and techniques for practical calculations of the QED effects will gradually emerge in more detail in the following discussion of many observables to be measured at FCC-ee.
The content of the paper is the following: we start with a brief overview of perturbative QED techniques. In particular the soft and collinear factorisation is substantially different from that of QCD. It will be stressed that optimal strategy of truncation of the perturbative orders from the point of view of precision is not the simplistic order-by-order truncation but a more subtle approach taking into account mass logarithms at higher orders. Next we shall elaborate on the QED component of the systematic errors in the measurements of the mass and partial width of Z, the total cross section and asymmetries near the Z resonance and of the W mass measurement, according to the present state of art inherited from LEP era, describing in detail the main source of QED uncertainties in these experimental observables. Next, we shall compare present QED precision of the observables with the planned experimental systematic and statistical precision at FCC-ee, in order to estimate how much improvement is needed in the QED calculations. For each experimental observable we shall examine how difficult it will be to get sufficient improvement in terms of more perturbative orders, improvement in the resummation techniques and development of new software tools i.e. the entire new class of dedicated MC event generators. In the above discussion we shall briefly elaborate upon the question how to factorize in practice the resummed QED and the so called "pure electroweak" parts of the perturbative calculation such that they coexists in the MC event generators forming a complete perturbative expansion without any double counting nor under counting. A short summary will complete the paper. while the remaining second order correction O(α 2 L 1 f , α 2 L 0 f ) can be safely neglected. This is also visualized in Fig. 2(a).

General features of QED corrections
The next precision level of 3 · 10 −4 requires adding O(α 3 L 3 f , α 2 L 1 f ), see also Fig. 2(b) for the illustration. It is only at the precision level 1 · 10 −5 (relevant for the FCC-ee) that including O(α 2 L 0 f ) looks necessary, but it should to be accompanied in addition with O(α 4 L 4 f , α 3 L 2 f ) corrections in order to be useful in practice. The case of muon pair differs only slightly, as can be read from the RHS plot in Fig. 1.

Optimal strategy for advancing QED precision
In view of the above estimates of the numerical importance of the generic perturbative QED contributions it is easy to answer the question about the optimal strategy for advancing the precision of QED calculations. For example, performing the complete O(α 3 ) calculation may look justified because it will be automatically gauge invariant, all soft/collinear limits, mass effects will be correct. However, for precision at the ∼ 0.5% level, a much simpler (almost trivial) calculation taking only up to O(α 2 L 2 f ) into account is enough, see Fig. 1. Similarly, the inclusion of O(α 2 L 1 f ) in order to gain precision ∼ 0.05% level is in practice useless, unless O(α 3 L 3 f ) contribution enters into the game, see Fig. 1. It is also a common misconception that including O(α 2 L 2 f ) forces us into the use of the strictly collinear approximation 2 . It was demonstrated that in the EEX matrix element [16] with O(α 2 L 2 f , α 3 L 3 f ) corrections can be implemented at the exclusive (unintegrated) level, i.e. keeping angle distributions for two and more photons accessible for experimental cut-offs.
(a) 0.6% Figure 2: QED perturbative leading and subleading corrections. Rows represent corrections in consecutive perturbative orders -the first row is the Born contribution. The first column represents the leading logarithmic (LO) approximation and the second column depicts the next-to-leading (NLO) approximation. In the figure, terms selected for the same precision level are limited with the help of an additional line.

More on soft and collinear resummation in QED calculations
The other important question is that of the importance of the soft and collinear photon resummation. For totally inclusive observables in the final state (like FSR in the Z boson decay) mass logarithms do cancel completely among real and virtual ones according to the Kinoshita-Lee-Nauenberg theorem. For the cut-off on the total photon energy E max E beam , the O(α n ) YFS/Sudakov double logarithmic contributions for a strong cutoff may easily be huge, S (n) 1, and definitely require resummation. Only for a loose cut-off like E max /E beam ∼ 0.2 and very low precision requirements one could avoid soft photon resummation.
There are several techniques of soft photon resummation of the different levels of sophistication. In the most primitive version all photon energies are restricted to small values, the integration over photon angles and energies is done analytically, keeping only the sum of photon energies fixed. This we call inclusive soft photon exponentiation, IEX in short. In the Yennie-Frautschi-Suura (YFS) work [17] it was outlined how to match smoothly soft and hard photon distributions, covering the entire phase space, without any need of the cut-off on total photon energy. It was only in late 1980's that it was gradually worked out how to implement YFS soft photon factorisation and resummation within the Monte Carlo event generators [18,19]. This technique we call exclusive photon exponentiation, EEX in short. In the EEX methodology soft photon factorisation and resummation is done at the level of the multiphoton fully differential distributions, that is for spin amplitudes squared and spin summed/averaged.
Near the narrow resonances R another similar class of soft logarithms present in real and virtual corrections also requires resummation. In addition, due to the complicated pattern of the QED interferences between the initial and final state photons (cancellations due to short lifetime of the resonance) which operate at the amplitude level, one is forced to perform soft photon factorisation at the amplitude level, before squaring and spin summing. This technique was developed by the Frascati group [20]. Moreover, the need of proper implementation of spin effects in the τ pair production and decays at LEP also was enforcing the use of spin amplitudes. In order to meet the above requirements a new variant of YFS-inspired soft photon resummation was developed [21], in which soft photon factorisation was reformulated at the amplitude level and soft photon resummation is implemented numerically within the MC program. Narrow resonance effects were accommodated as in ref. [20,22]. This technique we call coherent exclusive exponentiation, CEEX in short. So far, the only implementation of CEEX technique is in the KKMC event generator [11]. The important message to theorists specialising in QED+EW multiloop calculations is the following: do not add soft real emissions to multiloop results in order to eliminate infrared singularitiesá la Bloch-Nordsieck, if you want these results to be used in the MC generators with IR resummation. Instead, you should subtract IR parts (YFS virtual formfactor) from the amplitudes, before squaring and spin summing 3 . Why? Because combining IR soft and real contributios and the differential cross section level is already done in the Monte Carlo.
The related important practical question is whether the use of collinear resummation of the mass logarithms is mandatory in QED to infinite order, like in QCD? Obviously in QED it is not mandatory and in practice it is usually enough to stop at some finite order, typically truncating infinite LO series at O(α 3 L 3 f ). In the FCC-ee environment it may be sometimes necessary to include O(α 4 L 4 f ). In case of the photon emission from leptons, employing the entire machinery of the collinear resummation technique up to LO+NLO or LO+NLO+NNLO level is not trivial due to finite lepton masses, handling properly the factorisation scale parameter, sectorisation of the phase space etc. In practice, it is more convenient and/or economical to perform soft photon resummation in first place and only then to include collinear resummation truncated to a convenient order 4 . This technique was used successfully in the many MC projects like KORALZ [23], KKMC [11], BHLUMI [24], YFSWW [25], KORALW [26], BHWIDE [27], YFSZZ [28], also in the semianalytic approaches like ZFITTER [12] and TOPAZ0 [13]. Moreover, soft photon resummation is of the highest possible priority for the resonant processes, obviously near the Z resonance, but also it will be mandatory in the W W production process, especially near the W W threshold.
Summarizing the above discussion: • The (complete) order-by-order perturbative calculation in QED is definitely not the economic way to to obtain predictions for cross sections or asymmetries with the precision below 0.1%.
• Soft photon resummation is an absolute necessity, especially for resonant processes. It exists in at least three different variants (IEX, EEX and CEEX).
• Do not follow Bloch-Nordsieck to eliminate IR singularities!
• Resummation of collinear mass logarithms ln(s/m 2 f ) is very useful, but in QED it is usually convenient to truncate it at some finite order.
• Approximation of small lepton mass m 2 f /s 1 should be exploited for electron and muons as much as possible, but for τ lepton ∝ m 2 τ /s, m 4 τ /s 2 terms may not be negligible at tree-level, while higher powers of this type in higher orders of α are probably irrelevant.

QED "deconvolution"
The so-called "deconvolution" of QED effects is the procedure of removing universal (process independent) QED effects from experimental data for the total cross section and angular differential cross sections in the process e + e − → ff (+nγ). It was a cornerstone of the final analysis of LEP data near the Z peak in ref. [29]. Cross sections, asymmetries, branching ratios, Z mass and partial widths derived from the "deconvoluted" integrated/differential experimental distributions were called EW pseudo-observables (EW-POs) and were (approximately) independent of the QED effects and experimental cutoffs, such that they could be combined among four LEP collaboration and SLD. The technique of EWPOs was defined and thoroughly elaborated in ref. [30].
N.B. The term "deconvolution" is a little bit misleading, because in practice it means fitting a certain theoretical formula to experimental data at one or several energies 5 .
It is tempting to assume that at FCC-ee one can use the same QED deconvolution as at LEP. However, this assumption may be wrong due to the much higher experimental precision at FCC-ee.
Let us comment briefly on some aspects of the factorisation of the soft and collinear QED corrections, which is the basis of QED deconvolution and an essential element in the construction of EWPOs. This is because one could worry that the entanglement of electromagnetic and weak interactions at multiloop corrections, the rise of non-factorisable interferences above the level of FCC-ee precision, may give rise to practical or principle problems with the clean separation of the QED universal corrections from the complete EW perturbative calculations, especially beyond the 1-st order.
In refs. [29,30] the simple version of QED ISR deconvolution relies on a well known simple convolution formula with the integration over a single variable (total ISR photon energy) over the product of two objects: the ISR radiator function 6 and Born-like differential or integrated cross section. In the case of a cut-off on the final fermion pair effective mass, the integration over the radiator function of FSR had to be included. The resulting double convolution formula for ISR×FSR was used in ref. [30] and also in ref. [16] for testing KKMC. In addition, in the ZFITTER and TOPAZ0 programs missing O(α 1 ) contributions were added into the game. In particular, missing initial-final state interference (IFI) and effects due to experimental cut-offs were also added in this way.
As it was discussed in ref. [30], the above simple treatment of IFI at LEP analysis was possible because near the Z peak, in the absence of strong cutoffs IFI is suppressed by an additional Γ Z /M Z factor, while away from the Z peak, where LEP data were of limited precision and IFI also could be either neglected or eliminated using an additive O(α 1 ) correction. The above treatment of IFI would be highly unsatisfactory at the FCC-ee precision.
The first mandatory thing on the way to improved treatment of IFI is soft photon resummation. The multiphoton convolution formula at the level of the matrix element in the soft photon approximation was constructed a long time before the LEP era, by the Frascati group [20,32], and it is a natural extension of the the exponentiation formula of Yennie-Frautschi-Suura [17] to a resonant process.
In the Frascati-type formula with soft photon resummation at the amplitude level the effective Born matrix element is clearly factorized out and could be exploited for constructing a better variant of the EWPO definition. However, as it was shown in ref. [33], after squaring, spin-summing, and integrating over photon angles the resulting convolution formula, has four convolution variables, two for ISR and FSR and two additional variables for IFI. It is not so handy as the traditional one, because the 5-dimensional (including cos θ) integration has to be done numerically using the MC method. In ref. [33] new MC code KKFoam implementing the above calculation was used to crosscheck the calculation of the KKMC program in the soft limit, where it should be by construction fully compatible with the Frascati approach.
In the final LEP data analysis [29] the coupling constants of Z, its mass and width inside the effective Born (differential) cross section were obtained from the fit to data taken typically only at 3-5 energies. The deconvoluted Born integrated cross section and asymmetries were not really coming directly from data, but were calculated from the (theoretical) fitted effective Born at s = M 2 Z . In other words, the resulting EW pseudoobservables (EPWOs) were encoded in the parameters inside the effective Born, EW ADLO [36] 80376 ± 33{6} 0.5 0.3 12 ±0.020{0.001} 1.0 · 10 −5 0.3 · 10 −5 100  [37]. They are all without theoretical uncertainty component. Last column shows improvement factor in QED theoretical calculations needed in order to be equal to experimental precision of FCC-ee measurements.
pseudo-parameters (EWPPs). The definition of the effective Born in ref. [29] was done at the spin amplitude level, see eq. (1.34) therein. Could the above scenario be repeated at FCC-ee using an effective Born defined at the amplitude level and factorizing QED at the amplitude level? The detailed numerical studies, as in ref. [30], at the FCC-ee precision level have to be done from the scratch in order to answer this question. In particular an additional uncertainty introduced by partial/incomplete inclusion of the SM effects in EWPOs extracted from data should be re-examined. In the recent ref. [4] it was argued that the combined use of more advanced versions of the ZFITTER/TOPAZ0 programs and of the MC programs of the KKMC class may provide solution. Summarizing: At the FCC-ee precision IFI requires resummation. Clean factorisation into a and a model-independent (Born-like) part still works at the spin amplitude level, even if it fails (due to non-factorisable contributions) at the amplitude squared level. Such a factorisation opens ways to a new more flexible definition of EW pseudo-observables, which would possibly cope with the FCC-ee precision, see subsection 5.7 for more details.

The most important QED-sensitive experimental observables at LEP and FCC-ee
The minimum improvement factor to be achieved in the precision of QED perturbative calculations for FCC-ee experiments, in order that these effects are controlled at the level of the pure experimental errors, are shown in the last column of Table 1. In the third column there we collect the values of the total error of the observables most sensitive to QED effects, as measured at LEP, including also explicitly QED theoretical uncertainty induced in the overall experimental systematic errors. These components are accompanied with the citations of the source papers from which they are taken. The projected much better experimental statistical and systematic errors at FCC-ee experiments are shown in 4th ad 5th column. They are taken from Table S.3 in ref. [2]. As a matter of fact, the actual improvement factor for FCC-ee should be even 2-3 times bigger than the one shown in Table 1, in order to be sure that QED effects are clearly a subdominant component in the corresponding overall systematic error for most of the observables measured at LEP experiments. The information in the Table 1 is the starting point for more detailed discussion in the following sections. The most sensitive to QED effects observables of the FCC-ee experiments listed in the table are total cross section, including low angle Bhabha for luminosity measurement, especially near the Z resonance (for Z mass and width, Z invisible width), cross section of Zγ final state (Z radiative return) above the Z peak, charge asymmetry for leptonic pair final states and spin asymmetry in the τ pair production. We are analysing a subset of experimental observables near the Z resonance and W W threshold which are most "vulnerable" to the QED effects -omitting many others.
Let us comment briefly on Table 1, before more detailed discussion in the following sections: QED uncertainty of M Z and Γ Z (derived from the lineshape) are taken from refs. [29,38]. Huge photonic corrections to the Z lineshape (30%) in all LEP experimental and theoretical studies are mastered using a formula derived in ref. [39]. However, the precision of QED corrections in ref. [29] is dominated by the uncertainty of the fermion pair correction of ref. [40], in spite of the fact that this kind of correction is small by itself. For more discussion on that see sect. 4.1. The QED uncertainty of R l Z is taken from [34] and is mainly due to t-channel 7 obscuring Γ ee . QED uncertainty of N ν from the Z radiative return is from the OPAL paper [35] and its error is enhanced due to the use of the outdated MC programs 8 with inferior precision δσ/σ = 2%. The factor 60 in last column is therefore a gross overestimate. On the other hand the precision of N ν from the peak cross section is limited by the luminosity error δL/L 0.6 · 10 −3 , mainly due to QED corrections for the low angle Bhabha luminometer 9 . Strictly speaking this luminosity error at LEP was mainly of the QCD origin because it was dominated by the vacuum polarisation contribution uncertainty due to the experimental error of low energy hadronic data.
The leptonic charge asymmetry of LEP experiments A lept. F B = 0.0171(10) from Table  2.13 in ref. [29] translates into EW mixing angle sin 2 θ ef f W = 0.23099(53), see Table 12.4 in ref. [34]). It is almost the same as for the muon pair final state alone. The theoretical error component was estimated in ref. [29] (Table 2.8) as δA 0,l F B = 0.01% from the difference between ZFITTER and TOPAZ0 see also ref. [30]. This difference provides for the technical precision of both programs and not what we need, that is uncertainty due to missing higher orders. Such an estimation of the QED perturbative uncertainty of A F B due to missing higher orders was provided in ref. [16]. From the comparison of the results of A F B (M Z ) from KKMC, O(α 1 ) KORALZ and ZFITTER (IFI included) in Fig. 28(b) for the loose cut-off on total photon energy the uncertainty of order δA F B 0.05% seems reasonable. It Table 1 from two spin asymmetries, P τ and A pol,τ FB , comes from the LEP summary of ref. [29] (Section 4.4). Both asymmetries have comparable experimental errors. QED effects in the MC simulation of the τ decays quoted there following ref. [42] are estimated to contribute δA τ = 0.0010 for P τ , which yields δ sin 2 θ ef f W = (1/8)δA τ = 0.00012. For A pol,τ FB the QED uncertainty according to ref. [29] is even smaller. The corresponding anticipated FCC-ee experimental statistical and systematic errors are taken from ref. [2], where it is shown for A pol,τ FB . We transformed it into the error of the EW mixing angle using relation δA pol,τ FB 6 · δ sin 2 θ ef f W . The combined LEP measurement of M W in ref. [36] comes mainly from the W mass reconstruction in the qqqq and lν l qq channels. The uncertainty of 6MeV due to higher order radiative corrections is estimated from comparison of the KandY [14,43] and RACOONWW [15] MC programs.
The LEP experimental error of A M Z ±3.5GeV being δA F B ∼ 2% is mainly statistical, and the QED error estimate δA F B ∼ 0.1% is taken from refs. [16,44].
Summarizing on the Table 1 one may say that the very minimum of the improvement on the QED calculations needed for FCC-ee near the Z resonance and near the W W threshold is typically a factor ∼ 3-60, with the exception of charge asymmetry, where a bigger improvement ∼ 100 will be necessary. It is, however, important to stress that in order to get to the same comfortable situation as at LEP, where QED uncertainties (except for the luminosity cross section) were typically at least factor 2-3 smaller than the statistical and systematic experimental error, one should really aim at the improvement factors being ∼ 3 times bigger than these of the last column in Table 1, that is a factor 10-180 and in the special case of A F B a factor 300. With values shown in the last column in Table 1 one will merely get into a rather uncomfortable situation like in the measurement of the luminosity and the total cross section at LEP1, for which theory uncertainty was comparable to the experimental error or even slightly bigger.

Present state of the art in calculating QED effects
In the following we are going to elaborate in a more detail on the present state of art in calculating of the QED effects, as in most cases inherited from the LEP era. It will be done for each observable listed in the tab. 1.

Cross section near the Z resonance -mass and width of Z
The most sizable QED corrections to the total cross section near the Z resonance i.e. the so called lineshape (∼ 30%), are due to multiple photon emissions from the initial e ± beams (ISR). The basic "effective radiator" formula for the ISR photonic QED corrections was defined in ref. [39], combining/summarizing results of several other works.
The master ISR formula of ref. [39] includes, in addition to the classic O(α 1 ) result [45], the relatively simple O(α 2 L 2 e ) contribution [46] and the more elaborate O(α 3 L 3 e ) corrections obtained in the analytical form 10 in ref. [48] and crosschecked using dedicated MC program [10]. The next correction O(α 4 L 4 e ) is also available, see ref. [49]. According to analysis of ref. [31], at the LEP precision it is negligible. ISR formula of ref. [39] also includes O(α 2 L 1 e ) photonic contribution from ref. [8], see alse refs. [50,51] for corrected O(α 2 L 0 e ) results. In addition to photonic corrections, it includes also the small contribution from fermion pair production of ref. [10]. It should be remembered that a similar but less precise ISR formula was presented earlier in the pioneering work of ref. [52]. Both ISR formulas of refs. [39,52] are employing soft photon resummation 11 .
The master analytical formula for ISR of ref. [39] was used in the numerical studies of [31,38] and in the analysis of all LEP experimental data and all related theoretical studies, as a basic tool for evaluating QED effects in the Z lineshape. In ref. [31] the total QED uncertainty due to photonic corrections was estimated to be δM Z , δΓ Z ≤ 0.1MeV.
The master convolution formula of ref. [39] with the ISR radiator function was neglecting contribution from the initial-final state interference (IFI). The uncertainty due to IFI to the lineshape observables was very carefully analyzed in ref. [38], exploiting older O(α 1 ) analytical work (without exponentiation) [53], and results from MC programs with [11] and without [23] soft photon resummation. The IFI correction was found to be of the order δσ/σ ∼ 0.2% for | √ s − M Z | ≤ 3 GeV to be linearly dependent on √ s and causing the shift of M Z of +0.17 MeV for lepton pair production and -0.17 MeV for hadronic final states, such that the net contribution of IFI to M Z is small due to the accidental cancellation. Possible uncertainties in the M Z measurement due to IFI-like missing O(α 2 ) and O(αα S ) corrections were discussed in ref. [38] quantitatively, concluding that the uncertainty of the IFI effect due to higher order corrections is δM Z 0.1 MeV and δΓ Z 0.1 MeV and of the total cross section at the top of the Z peak is merely δσ/σ < 10 −4 . From the analysis of the IFI effect in ref. [38] it was also clear that the size of the IFI effect and its uncertainty depend quite strongly on the centre of the mass (c.m.s) energy √ s and on experimental cut-offs, hence the role of the MC in controlling it was already quite important in the LEP data analysis. Similar analysis of the IFI effect in the lineshape observables was performed in ref. [30]. It was restricted to pure O(α 1 ) QED semi-analytical calculations using ZFITTER [12] and TOPAZ0 [13] programs, without soft photon resummation.
The mass of Z boson comes from fitting of the cross section across the Z resonance (lineshape). According to the final LEP1 and SLD data analysis of ref. [29], the dominant QED contribution to the lineshape, which translates into δM Z 0.2 MeV and δΓ Z 0.3 MeV, was not from photonic corrections, but rather from the light lepton and quark contributions calculated in ref. [40]. In the FCC-ee context, this contribution should be cross-checked and carefully re-analyzed.

Luminosity measurement
The overall error of the luminosity measurement quoted by all LEP experiments for LEP1 was below 0.05%, with the best experimental luminosity error 0.034% being that of the OPAL collaboration. The theoretical prediction of the QED dominated low angle Bhabha (LABH) cross section in all LEP collaborations was based on the calculation from BHLUMI 4.04 Monte Carlo event generator published in ref [24]. Its precision, following ref. [54], was quoted to be δσ/σ 0.061%. In fact, the luminosity cross section was the only observable in LEP experiments for which the theory uncertainty was bigger than the experimental statistical and systematic error. As seen in tab. 2 of ref. [55] its biggest component was in fact not of the QED origin but due to the vacuum polarisation effect, which was calculated using low energy experimental hadronic data. Since in the recent years the error of vacuum polarisation contribution was reduced significantly, the updated error of BHLUMI 4.04 predictions is now estimated to be 0.038% (see the same Table in ref. [55]) and is now dominated by uncertainties due to missing higher order perturbative QED corrections.
The principal luminosity measurement at FCC-ee will be done using the same type of the low angle Bhabha process as in LEP experiments, using a similar calorimetric detector [56]. Prospects of improving QED predictions for the FCC-ee luminometry will be discussed in the following.
In ref. [55]) ( Table 3) it was pointed out that due to twice wider angle of the FCC-ee Bhabha luminometer (64-86 mrads) than at LEP, the actual theoretical error according to the present state of the art (inherited from LEP) would be in fact 0.090% due to the bigger Z-exchange contribution. However, this error can be reduced to the negligible level already now using the BHWIDE Monte Carlo [27], see ref. [55]) for the details.

EW mixing angle from charge and spin asymmetries at LEP
All charge and spin asymmetries measured at LEP have their errors dominated by the statistical error [29,34,36,57,58]. Let us characterize briefly these (yet subdominant) QED corrections in the LEP measurements.
Charge and spin asymmetries depend on the ratios of the Z couplings. With the usual simplifications [34] and/or introduction of pseudo-observables (PO) [30], in which QED effects are "deconvoluted", it is convenient to express all charge and spin asymmetries (following notation of ref. [29]) in terms of In this work we consider leptonic charge asymmetry and two τ spin asymmetries The above defines convenient pseudo-observables (pseudo-parameters) sin 2 θ l ef f , l = e, µ, τ which allow us to compare information on Z couplings from various asymmetries. Using sin 2 θ l ef f = 0.2315 of tab. 1 (A l = 0.1472) and assuming lepton universality, the following simple relations relate uncertainties of asymmetries and of the EW mixing angle: How big were estimates of the theory uncertainties δA l F B , δA pol.τ FB and δ P τ due to QED effects in data analysis at LEP near the Z resonance and above?
The first systematic study of the QED uncertainties in A F B was attempted in the initial LEP workshop [59], in particular the value of the strongly suppressed IFI contribution A µ,IF I F B = 5 · 10 −4 (following ref. [60]) was established. As already said, at the Z peak ref. [29] cites the difference δA l F B = 10 −4 between ZFITTER and TOPAZ0 as a theory uncertainty estimate. A closer look into ref. [30] on which ref. [29] is based reveals that the story is more complicated. For IFI switched off, from the differences ZFITTER−TOPAZ0 and due to the change of the ISR radiation function (from the factorized to additive form) shown in ref. [30] one may indeed quote δA l F B (M Z ) 10 −4 and δA l F B (M Z ± 3GeV ) 3 · 10 −4 , especially for loose cut-offs 12 . On the other hand, for IFI switched on, the difference between ZFITTER and TOPAZ0 in ref. [30] for realistic cut-offs, within | √ s − M Z | < 3GeV , were quite sizable, of order of a few per mille 13 . Fortunately, in ref. [62] the source of these discrepancies between IFI implementation in ZFITTER and TOPAZ0 were identified and new precision estimates δA l F B (M Z ) = 0.0002 and δA l F B (M Z ± 3GeV ) = 0.0013 were provided. However, this estimate is not reliable for higher order effects in IFI, because both ZFITTER and TOPAZ0 implement essentially the same additive combination of the ISR and FSR collinear radiator functions with O(α 1 ) results, integrated analytically over a single real photon within some simple cut-offs. Hence the above estimates really represent the technical precision of both programs and not their physical precision.
As already noted in ref. [30] the necessary next step should be simultaneous exponentiation of ISR and IFI. This goal was achieved, almost in parallel with the above works, in the KKMC program [11]. In ref. [16] the comparison of KKMC with ZFITTER and O(α 1 ) KORALZ has shown that indeed for loose cut-offs one may conservatively estimate higher order QED effects to be δA µ F B (M Z ) = 0.0005. According to eq. (4.3) this estimate translates into δ sin 2 θ e ef f = 2.8 · 10 −4 . This QED uncertainty was about factor two below the LEP experimental precision at the Z peak, see Table 1.
Later on, at the LEP2 energies near and above the W W threshold the experimental LEP precision of A F B was worse due to smaller statistics, however, the IFI contribution was again quite important because it was not suppressed. For mild experimental cutoffs it could reach a few percent. The precision of the traditional O(α 1 ) calculations of IFI 12 Similar earlier study of ref [61] based on comparisons of analytical calculations versus of KORALZ Monte Carlo was quoting δA l F B (M Z ) = 0.0005 and δA l F B (M Z ± 3GeV ) = 0.005. 13 The authors of ref. [30] were recommending subtraction of IFI using MC programs like KORALZ, with complete O(α 1 ). was believed to be ∼ 1%, see ref. [62]. More systematic studies of the QED uncertainty, including IFI, at LEP2 energies were done in refs. [16,36,44]. In ref. [44] the overall QED uncertainty based mainly on comparisons of KKMC with ZFITTER and KORALZ was estimated to be δA µ F B (M Z ) = 0.004 − 0.005 and for the IFI component δA µ,IF I F B (M Z ) < 0.003. The more detailed study of ref. [16] has concluded that the overall QED uncertainty (including IFI) of the charge asymmetry prediction from KKMC is δA µ F B (M Z ) = 0.002 over the entire LEP2 energy range.

The invisible Z decay width from cross section and radiative return in LEP experiments
In LEP experiments the measurement of the Z invisible decay width was quantified as the deviation from the Standard Model expectation of the neutrino number and family generation number, N ν = 3. This parameter was measured to be N ν = 2.984 ± 0.008 [29]. It was deduced mainly from the total cross section at the Z peak σ(M Z ). The dominant contribution to its uncertainty was the luminosity error. According to ref. [29] the luminosity error contributed δN ν = 7.5 δL L . For δL L = 6.1·10 −4 it gives δN ν = 4.6·10 −3 . At the FCC-ee the above luminosity uncertainty will be again the dominant one. Getting δN ν = 1 · 10 −3 will require the luminosity error to be improved by a factor ∼ 5, down to δL L = 1 · 10 −4 , which looks feasible, see Sect. 5.3.
Another kind of measurement, the so called Z radiative return (ZRR) was also exploited in LEP experiments, see for instance ref. [35], where the less precise result N ν = 2.69 ± 0.15 was quoted. In the ZRR process one is examining the energy distribution of a photon emitted above some minimum angle from the beams, at the c.m.s. energy well above the Z resonance. Such a photon is also seen for Z decaying into neutrinos or other "invisible" particles. The Z resonance peak is seen in the energy distribution of the ZRR photon and N ν is obtained by means of comparing the integrated cross section of this process with the result of a reference Monte Carlo program. The error of the above LEP measurement δN ν = 0.15 was dominated by the statistical error. The QED component of this error was δN ν = 0.06, according to δN ν = 3δσ/σ, where δσ/σ = 2% was the QED/SM error attributed to the cross section calculations of KORALZ [23,63] and NUNUGPV98 [64] programs. Such a precision was then quite satisfactory in view of the large statistical error. The calculation of ZRR using the more advanced KKMC program [11] was not yet available at that time. Another theoretical study was done in ref. [65]. The improvement at FCC-ee down to δN ν = 1 · 10 −3 would require δσ/σ = 3 · 10 −4 precision for the MC programs calculating ZRR process (factor ∼ 60 improvement). However, with the advent of KKMC this factor is clearly an overestimate, see the following discussion.
With FCC-ee precision the error of the luminosity measurement will enter into the game and will have to be improved as well. On the other hand, thanks to high luminosity of FCC-ee it would be possible to eliminate the dependence on the luminosity error in the ZRR method by means of using the ratio of the photon distribution with invisible Z decay and with Z decaying in the muon pair. Such a method could not be exploited at LEP due to the limited statistics of the ZRR with muon pairs. At FCC-ee it makes sense and the prospects of its precision in terms of N ν are discussed in the following Sect. 5.5.

QED in W -pair production at LEP, M W measurement
The role of QED in the LEP2 data analysis near and above the W W threshold was somewhat different than near Z pole -there was no systematic attempt to "deconvolute" universal (process independent) QED effects and develop the technique of pseudoobservables -QED was usually kept together with pure electroweak corrections and effects due to semi-classical QED interaction of two massive W near the production threshold, the so-called Coulomb effect, had to be included in the game. Generally, there are four classes of QED effects in e + e − → W + W − or e + e − → 4f process: initial state corrections (ISR), final state corrections (FSR) in the decays of two W ± , final state Coulomb corrections (FSC) and the so-called non-factorisable interferences 14 between ISR and FSR in two W ± decays (IFF). The IFF corrections are suppressed due to relatively long lifetime of W 's. The effects due to ISR are numerically the biggest but easier to control, while the FSR effects can be also quite sizable for typical experimental cut-offs.
The experimental error of the measurement of the W boson mass from the total cross section of e + e − → W + W − near threshold at LEP2 experiments [36] was δM W = 200 MeV due to poor statistics of the data, quoting theoretical error in this measurement as negligible.
Much better experimental precision of δM W = 34 MeV was achieved in LEP2 experiments [36] from W mass reconstruction for lνqq and qqqq final states. For this method, the uncertainty of the O(α 1 ) theoretical calculations used for e + e − → W + W − process was estimated in ref. [36] from the difference between KandY [14,43] and RACOONWW [15,66] programs. In this way, for lνqq channel δM W = 8 MeV was obtained and for qqqq it was δM W = 5 MeV. In particular, the uncertainty due to ISR radiation alone δM W = 1 MeV was obtained in ref. [36] using KandY, RACOONWW and WPHACT [67] by means of switching on/off O(α 3 L 3 e ) ISR contribution. The general discussion of the theory issues in the precision SM calculations for the e + e − → W + W − process can be found in ref. [68]. However, there were many other works focusing on various specific issues. For instance the detailed analysis of the ISR effects can be found in ref. [69]. The ISR effect on the total cross section of e + e − → W + W − near the threshold, at 160 GeV is −28% and merely −7.5% at 205 GeV. However its uncertainty deduced from switching on/off non-IR O(L 3 e α 3 ) is not bigger than of 2.5 · 10 −4 at 160 GeV and 2.1 · 10 −4 at 205 GeV. The entire −28% ISR effect at 160GeV would translate into huge δM W ∼ 400 MeV, while its uncertainty is worth only δM W ∼ 0.2 MeV.
Dedicated analysis of the QED and non-QED effects in the M W reconstruction from the final states can be found in ref. [70] 15 . The uncertainty due to not included higher order and nonleading ISR effects was estimated in Table 1 of ref. [70] to be δM W < 1MeV, the FSR uncertainty was rated as δM W ∼ 2MeV and the IFF (following ref. [71]) at δM W < 2MeV. The above analysis is based on the results from KandY toolbox of the Monte Carlo and semianalytical programs [14].
SM calculations for e + e − → 4f process actually used in the LEP2 data analysis were all in form of the Monte Carlo event generators and they included both QED and the remaining EW corrections. In the following we shall provide brief descriptions of them commenting also on how the QED part was organized/implemented. There were two such MC codes with O(α) corrections to the signal doubly resonant CC03 graphs calculated respectively in the leading pole (LPA) or double pole (DPA) approximationand with the tree level matrix element for the remaining e + e − → 4f contributions: [25,72] generates the signal process e + e − → W + W − → 4f according to the LPA scheme, exploiting the O(α 1 ) calculations of Refs. [73][74][75][76] and with W ± decays simulated independently. Multiphotonic radiation for production part, e + e − → W + W − , is implemented in the YFS framework (EEX scheme). The hard photon ISR is corrected to the O(L 3 e α 3 ). The FSR in two separate W -decays is handled by PHOTOS. The KORALW 1.42 [26] code, which simulates the complete e + e − → 4f process at the tree level, has been combined on the event per event basis with the YFSWW3 code. This way a concurrent MC called KandY of ref. [14] has emerged, which simulates both the complete four-fermion final states and includes O(α) corrections to the W-pair production and decay. The semianalytical program Korwan for the improved Born approximation (IBA) is included for testing the main MC.
• RACOONWW [15] simulates the complete e + e − → 4f Born-level process, and the single real photon emission process e + e − → 4f + γ and implements O(α) electroweak virtual corrections in the DPA scheme [77,78] exploiting one-loop calculations for onshell W W production and decay. The ISR radiation is based on the QED collinear structure functions to second order with soft photon exponentiation. It includes also semianalytical program for the Improved Born Approximation (IBA) [79], based on CC03 graphs (doubly resonant) and universal corrections to these graphs, i.e. the Coulomb correction at the threshold, the running of effective couplings and collinear ISR.
Both programs include Coulomb corrections for off-shell W ± bosons. According to expert comparisons of the above two aproaches in refs. [68,80], DPA and LPA methods 16 of inserting O(α 1 ) EW corrections into doubly-resonant part of the matrix element in the e + e − → W + W − process are basicaly equvalent. The main differences are in the implementation of the QED part of the matrix element. The complete O(α) corrections to e + e − → W + W − cannot differ -it was checked in ref. [68] that the difference for the total cross section between implementations in YFSWW3 and RACOONWW of the virtual plus soft corrections in two programs is below 0.01%. The overall agreement of YFSWW3 and RACOONWW, including all physical effects, is of the order of 0.3% for the total cross section at 200 GeV [68]. The differences of the 16 Differences between two variants of LPA are explained in ref. [43].
O(αΓ W /M W ) are expected due to inefficiency of DPA/LPA. However, most of these differences arise from the different treatment of the QED corrections in both approaches, so it is interesting to look into them.
QED calculations in RACOONWW rely on the massless e + e − → 4f + γ matrix element, hence the complete standard full phase space cannot be used. Instead, two methods are used to deal with emission of photons collinear with fermions. In one of them the phase space is "sectorized", i.e. the real photon phase space is restricted with minimum angle of the photon to fermions (also minimum photon energy to exclude IR divergence) and the contribution below the minimum angle is integrated over analytically, recovering correct fermion mass dependence. In another method [81] a QED variant of the Catani-Seymour (CS) subtraction scheme [82] is used. In both these methods collinear ISR or FSR photon angular distributions are integrated over and effective longitudinal momentum distributions (effective radiator functions and/or PDFs) arise. Their original O(α 1 ) version is upgraded to include LO+NLO corrections up to O(L 3 f α 3 ). Soft photon resummation is also included in the effective radiator functions, but multiphotons are not present explicitly in the MC events. Non-factorisable interferences between the production and two decays of W s are reproduced up to O(α 1 ) (in the soft photon approximation) The methodology of QED treatment in YFSWW3 or KORALW is quite different from that in RACOONWW. Multiple photons are present explicitly in the MC events. The matrix element is constructed following EEX variant of the YFS soft photon factorisation and resummation. Non-soft higher order LO and NLO collinear photon universal contributions are added to matrix element in the exclusive (unintegrated) form, without collapsing to collinear photon distributions, to a δ(θ) function. The multiple photon emission is included in YFSWW3 for ISR and FSR out of W s, while in KORALW it is restricted to ISR. Single and double photon emissions in the decays of W s are added using PHOTOS [42,83] program 17 . Non-factorisable interferences between production and two decays of W s are not included. The YFS technique of factorizing and resuming the universal QED corrections employed in YFSWW3 provides for a clear separation of the universal QED corrections from the rest of the SM radiative corrections at any order, hence it is a very good candidate for the future inclusions of the non-QED EW corrections beyond the O(α 1 ), or in any attempt of better organisation (deconvolution) of the existing O(α 1 ) calculations.
After LEP2 data were analyzed, the complete O(α) corrections to e + e − → 4f have been also completed in Refs. [77,78] for the final states without repeated flavours and without CKM-suppressed states 18 . The comparison of the above new calculation with older LPA/DPA results gives new insight into their uncertainty beyond universal QED corrections. Generally these newer calculations confirm the precision estimates of the LPA/DPA approaches of the LEP2 era. Near the W W threshold the difference between the DPA/LPA and the complete 4-fermion O(α) is ∼ 2% of the 4f -Born. Far from the W W threshold the difference DPA vs. complete 4-fermion O(α) is smaller, drops to below 0.5%, in accordance with the stated precision of the relevant MC codes. (It rises back to 1-2% at 1-2 TeV energies.) Refs. [77,78] provide also estimation of the missing higher order EW corrections, with the important conclusion that higher order EW corrections are dominated by the QED contribution α 2 log(m 2 e /s) and are estimated at ≤ 0.1% for energies below 500 GeV. The other interesting finding is that the contribution of the higher order Coulomb effect is of the order of 0.2% at the threshold. 19 . The total precision of these four-fermion O(α) calculations has been estimated by the authors to be a few per mille.
Another calculation specialised to the near W W threshold energy, using the effective field theory (EFT) technique is reported in Refs. [85,86]. The dominant NNLO corrections to four-fermion process (µ −ν µ ud) were calculated using EFT for unstable particles -the best calculation was nick-named as N 3/2 LO EFT , because in EFT a different expansion parameter (relativistic velocity of W in the W W rest frame) is used to count the strength of particular corrections. The drawback of the EFT method is that it provides the inclusive results only. The effect of these pure N 3/2 LO EFT corrections on the W mass is estimated as 3 MeV.
The complete set of graphs for the e + e − → 4f process at the tree level was implemented in several MC codes [68]. Two of them were used in the actual data analysis of LEP2: standalone version of KORALW [26] and WPHACT [67]. The latter one implements the so-called "Fermion Loop" gauge restoring scheme of ref. [87].
The present state of the art (mostly inherited from LEP) can be summarized as 0.2% theoretical precision for the total cross section of the W -pair production process in the entire energy range. If statistics was not limited at LEP2, that would translate into δM W 3 MeV for the M W measurement from the threshold cross section. Looking at the future needs of the FCC-ee we can see that a factor of 10 improvement in precision is needed relative to the present state of the art. This is a relatively moderate goal as compared to other observables, see Table 1. We will discuss its feasibility in Sect. 5.6.

Prospects of the QED calculation improvements in the FCC-ee measurements
In contrast to the previous section, where we have elaborated on the present state of the art in the calculations of QED effects for observables selected in Table 1, in the following we shall list all possible developments in QED calculation needed for analysing future data at FCC-ee. Obviously, this will be to some extent speculative and it is quite probably that we are going to miss some new calculation fronts or methods which will really emerge in the future.

Cross section near the Z resonance and mass of Z
As seen in tab. 1, the hadronic cross section σ 0 had at the Z peak was measured with 0.09% total error with 0.06% component dominated by the theoretical uncertainty in the luminosity measurement. Anticipated factor 6 improvement at FCC-ee is expected due to reduction of luminosity error down to 0.01%, both in theory and experiment. According to ref. [31] the uncertainty due to QED ISR uncertainty in fitting σ 0 had to the experimental lineshape is 0.02%, while the uncertainty of the IFI contribution [38] is even smaller. The [40] quasi-photonic light fermion contributions (from electron, muon, τ pairs and light quarks) to the ISR radiator function also contribute in addition 0.02% to the uncertainty of σ 0 had . Most likely this is an overestimate because this uncertainty is approaching half of the effect itself.
The optimistic point of view might be that for the hadronic cross section σ 0 had , due to almost 100% detector acceptance, known QED calculations are sufficient to reach the 0.01% precision needed for FCC-ee -only better testing of the existing calculations and better estimates of the missing h.o. corrections are needed. In particular mixed corrections QED-QCD should be reexamined. The improvements on light fermion pair corrections require some effort but this looks feasible. However, in order to get back to the LEP situation where QED uncertainties were factor ∼ 2 − 3 smaller than experimental one then more work would be needed. This would also imply improvements in the luminosity measurement beyond what is described in Sect. 5.3, which will be hard to achieve.
On the other hand, R Z l = Γ h /Γ l is free of the luminosity problem but more sensitive to QED uncertainties due to relatively complicated acceptance, with the isolation cones around beams, more restrictive cut on the total photon energy (or on acollinearity), and the t-channel contribution for the e + e − final state. The factor 12 improvement in tab. 1 for R Z l from 0.60% down to 0.05%, or even more if the QED uncertainty is to be kept below the experimental error will require better calculations. The most demanding will be the FSR class of corrections, but ISR and IFI will have to be studied at the level close to the 0.01% precision level. Semianalytic calculations of the class like ZFITTER [12] and TOPAZ0 [13] used at LEP will not be sufficient for the task. The upgraded version of KKMC or a similarly accurate Monte Carlo program will have to be developed and used. Dedicated study is needed to determine what level of the improvement of the perturbative calculation is really needed to match FCC-ee precision of R Z l . Most likely it will be one order more accross precision boundaries in Fig. 2 at various perturbative orders in powers of α and mass logarithms for non-soft (IR-finite) QED corrections in the matrix element of the MC program.
The mass of the Z comes from fitting the cross section across the Z resonance (lineshape) such that most of the QED effects are removed. It is the hadronic cross section which matters mostly. In the final LEP1 and SLD analysis of ref. [29] the estimation of QED uncertainties in the measurement of M Z was taken from refs. [31] for ISR, from [38] for IFI and from [40] for light fermion pairs. According to ref. [31] the uncertainty due to photonic ISR corrections to light fermion pairs process is negligible because it is weakly dependent on the c.m.s. energy. This statement has to be cross-checked. The effect of IFI according to ref. [38] is sizable (∼ 0.17 MeV) and its uncertainty induces δM Z , δΓ Z ≤ 0.1MeV errors. Light fermion pairs also contribute δM Z , δΓ Z ≤ 0.1MeV. Both of these QED uncertainties are therefore of the size of the FCC-ee experimental errors, hence for the comfortable data analysis they should be reduced by at least a factor 2-3. The IFI effect in the M Z measurement has a similar size and the opposite sign for the leptonic and hadronic final states. It is strongly dependent on the experimental cut-offs. Calculations of O(α 1 ) without exponentiation used at LEP for evaluating IFI will not be sufficient at the FCC-ee precision. For extracting M Z and Γ Z from σ h (s) near the Z peak further progress in reducing uncertainty of IFI and light fermion pair contributions in M Z beyond the LEP state of the art will be necessary. It is possible that for this particular purpose hybrid approach of LEP era combining the use of MC event generators and semianalytical programs like ZFITTER/TOPAZ0 will still work. However, the alternative approach based entirely on the MC event generators will serve as crosscheck and it will also have an advantage to work for more difficult case of charge and spin asymmetries. For the moment KKMC is the leading candidate for further studies of the IFI effects. However, its CEEX matrix element should be upgraded to include O(α 3 L 3 e ) and non-IR parts of O(α 2 L 1 e ) pentaboxes. For improvements of the light fermion pair contributions one should exploit programs dedicated to four fermion final states, the same as for the production W W process.
Summarizing, reduction of QED uncertainties in σ 0 had below 0.1%, in R Z l below 0.05%, in M Z and Γ Z below 0.1 MeV is definitely feasible, but requires more work and improvements of the existing MC and semianalytical programs. Improvements on light fermion pair contributions seems to be the most urgent.

Charge and spin asymmetries at FCC-ee
Let us concentrate on symmetries at the Z peak, where they will be measured most precisely. As seen in Table 1, the charge and spin symmetries will be measured up to a factor ∼ 50 more precisely than at LEP, which will require MC tools for calculating differential distributions with realistic experimental cut-offs including O(α 2 ), possibly even O(α 3 ), electroweak, QCD and non-universal QED corrections, while resummed QED universal corrections will have to be included at even higher orders. Concerning feasibility of such higher order QED and EW calculations, let us cite statement in the summary of the recent workshop proceeding [4]: "The techniques for higher order SM corrections are basically understood, but not easily worked out or extended ... We anticipate that at the beginning of the FCC-ee campaign of precision measurements, the theory will be precise enough not to limit their physics interpretation." Presently, the only MC tool (event generator) which provides predictions for all charge and spin asymmetries for arbitrary experimental cut-offs, including O(α 1 ) EW corrections, complete O(α 2 ) QED corrections and soft-resummed universal QED corrections to infinite order, is the KKMC Monte Carlo [11], at the Z peak and far away from it. Its precision goes far beyond what was needed at LEP, as seen for example in the recent study of the IFI contribution to A F B (M Z ± 3.5GeV), where the precision δA F B ≤ 10 −4 was verified using auxiliary calculations. However, for the FCC-ee experimental precision of asymmetries quoted in Table 1, the QED part of matrix element has to be upgraded to include the next orders, up to O(α 3 L 3 f ) in the CEEX matrix element, and EW corrections should be upgraded in the MC matrix element to the level of known complete O(α 2 ) corrections [88][89][90][91]. Needless to say, TAUOLA MC used in KKMC for τ lepton decays will require additional testing and development. Of course, development from the scratch of another MC program of similar quality as KKMC would be of great help in the solid independent validation of the required precision.
On the other hand, let us point out to some important problems with the modelindependent representation of the data in a form of simple pseudo-observables like in Table 1, where the information on the fermion couplings to Z extracted from charge and spin asymmetries is quantified in terms of a single parameter, sin 2 θ ef f W . At LEP it was possible. The big question is whether for much higher precision at FCC-ee it will be still possible? In the methodology of ref. [29] a parameter sin 2 θ ef f W is just in one-to-one correspondence with the ratio of the (real) vector and axial Zff couplings in the effective Born, which was fit to the e + e − → ff data (taking into account factorisable QED corrections). 20 Whether this kind of "effective Born" will be still effective in parametrising FCC-ee data in the (SM-)independent way is an open question. Some known SM effects, could invalidate it, if they are not numerically small as compared to FCC-ee data precision. Most important among them is the s-channel non-resonant contribution, which at the √ s = M Z drops out because it is almost exactly real, while the Z contribution is purely imaginary. The reduction of the effective c.m.s. energy due to ISR would possibly invalidate that. The electroweak W W and γZ boxes also go beyond the simplistic effective Born ansatz. The S-matrix approach 21 , in which Zff couplings form the residue of the Z-pole, provide nice justification of the effective Born ansatz. These couplings in the S-matrix approach have small imaginary parts, which are already included in the LEP definition of the effective Born (see eq. 1.34 in ref. [29]). The effects due to EW boxes and imaginary parts of Z couplings were proven in ref. [30] to be small as compared to LEP data precision, but this might be not true in FCC-ee environment. Note also that they are unlikely to be affected by new physics.
In particular, although ISR is reducing c.m.s. energy by ∼ 100MeV, this effect is controlled to within 0.1 MeV [31], hence its effect in charge asymmetry of order 1% is probably controllable within the FCC-ee data precision. The effect of QED initial-final state interference is suppressed 22 at the Z peak, δA IF I F B (M Z ) 0.0005 [59], hence its calculation at the O(α 2 ) is probably adequate to take it into account -better quantitative study of the uncertainty of A IF I F B beyond O(α 1 ) at the Z peak is definitely needed. Neglecting imaginary parts was found in ref. [30] to induce δA F B 0.15%, hence it will be not negligible for FCC-ee precision at the δA F B 10 −5 level. On the other hand, EW boxes are a little bit less problematic, as they were found [30] to induce a δA F B 10 −4 effect only.
Summarizing, in the FCC-ee data analysis it is likely that the LEP-style definition of sin 2 θ ef f W will have to be either completely abandoned, or replaced by some variant in which all the above effects are not neglected but taken into account. So far, there was no detailed analysis concerning this issue. One may only guess that model-independent subtraction of the s-channel non-resonant contribution, using data far away from the Z resonance, may still work. Fitting the imaginary parts of the Z couplings in the effective Born to data, that is treating them as additional pseudo-observables, would be a brave decision. Including W W and γZ boxes and other O(α 1 ) EW (QCD) corrections in the effective Born, i.e. removing them from data on the way to pseudo-observables would really mean treating them the same way as QED and the major change in the meaning of the EW pseudo-observables (better say the EW pseudo-parameters). See also related discussion in subsection 5.7.
The τ lepton spin polarisation from the τ pair production at Z is less prone to the ISR QED effects, as it is weakly dependent on √ s. One specific source of the uncertainties for τ spin asymmetries is due to the limited quality of the τ decays τ → ν τ π, ν τ ρ used as a spin polarimeter. The effects due to nonperturbative QCD and QED effects in the τ decays should be controlled better than at LEP. However, it looks that achieving δ sin 2 θ ef f W 5 · 10 −5 using τ spin polarisation is within the reach of presently available MC tools like TAUOLA and PHOTOS [42,83], provided some extra testing is done [92]. In particular one should "calibrate" τ decay polarimetric features, implemented in τ decay simulation MCs, using high statistics τ decay samples from Belle experiments 23 . The influence of EW boxes or imaginary parts in the effective Z couplings on the pseudoobservables related to the tau spin asymmetries was not studied in the LEP era, because their precision was statistically limited. More quantitative studies are needed.
Altogether, it is expected that sin 2 θ ef f W from the τ lepton spin polarisation will be measured at FCC-ee as precisely as from charge asymmetries and will provide a powerful 22 Outside the Z resonance δA IF I F B 0.01. 23 At the c.m.s. energy of the Belle experiments single taus are not polarized, in the tau-pair production allows for this kind of testing. crosscheck on both measurements, due to very different experimental systematics. The QED induced uncertainties need to be re-examined both in the τ production and decay processes, but no serious problems are expected.
Just one example of possible problem: the decay τ ± → π ± π 0 ν τ (γ) is an efficient τ spin polarimeter, where the additional photon emission has to be taken into account very precisely. Hovewer, the above radiative process has to be distingushed from the cascade decay the cascade decay τ ± → π ± ων τ , ω → π 0 γ, which has a non-negligible combined branching ratio of 0.0015 [93]. The future high precision Belle II data should be used to analyse precisely this and other similar effects in the energy spectra of τ decay product used for measuring τ polarization at the FCC-ee precision level.
Summarizing, it looks that in general the main problem is not so much in the better QED and SM calculations of various asymmetries, but rather in the very survival of the methodology of pseudo-observables (pseudo-parameters) used in the model-independent representation of data for these asymmetries.

Luminosity measurement
The luminosity measurement at FCC-ee will be based again on the low angle Bhabha process [56], but one should also remember about the possible use of the e + e − → γγ process for the FCC-ee luminometry, which is statistically limited, but not vulnerable to uncertainty due to vacuum polarisation, see Sect. C.5 in ref. [4] for more details.
By the end of the LEP era any substantial improvement of the theoretical prediction for the low angle Bhabha (LABH) process used to determine the LEP luminosity was effectively blocked by the large uncertainty of the vacuum polarisation, which was in fact of the QCD origin, or more precisely due do experimental errors of the low energy hadronic data. Since then, this uncertainty was reduced by a factor four and by the time of FCC-ee experiments another factor two is probable. With the present vacuum polarisation error, the LEP luminosity error would reduce from 0.061% down to 0.038%, see tab. 2 in ref. [55]. In this way h.o. perturbative QED components in the uncertainty of the LABH cross section get dominant.
However, any further progress will not be possible without solid control of the so called technical precision i.e. any problems due to programming bugs, mistakes in the MC algorithm, numerical instabilities. The Monte Carlo event generator BHLUMI 4.04 of ref. [24] was the subject of many internal tests, in particular using elaborate comparisons with semi-analytical calculations in ref. [94]. These powerful crosschecks were unfortunately limited to not so realistic kinematic cut-offs. The FCC-ee luminometer detector similarly as at LEP will select/detect events in a way which cannot be dealt with using analytic methods. The Monte Carlo is the only way to implement them in the theoretical calculations. Due to all the above features of the luminosity detectors, the only reliable way of controlling technical precision for realistic event selection is to compare calculations from two different MC programs of the comparable quality. In the LEP era there was no other MC program of the same quality as BHLUMI 4.04, hence QED and technical uncertainties were lumped together. Luckily, it seems that there are two such candidates for the next generation improved MC's for LABH process, BHLUMI with the upgraded QED matrix element proposed in ref. [24] and another MC program BabaYaga [95], developed in recent years, provided its QED matrix element is upgraded.
Assuming that the problems of the low precision of the vacuum polarisation and of poor control of the technical precision will be solved, then yet another non-QED component in the budget of the expected theory uncertainty of the LABH process following the LEP era state of the art starts to dominate. As pointed out in ref. [55], due to the angular range of the FCC-ee luminosity detector being factor two wider, the contribution from the Z exchange near the Z peak will rise by factor of four. If one estimates its uncertainty the same way as at LEP, then the total luminosity error would jump to 0.097%, see Table 3 in ref. [55]. However, this uncertainty is rather easy to reduce using the MC program BHWIDE of ref. [27] developed for the wide angle Bhabha process. In ref. [55] it was shown that the uncertainty due to the Z exchange can be reduced substantially, perhaps down to 0.001%.
In this context it is worth to mention that vacuum polarization for negative Q 2 can be obtained directly using electron-muon scattering process. In the experiment proposed in ref. [96] with 150GeV muon beam scatering on fixed target electron it will be possible to measure directly QED effective coupling within the range 0 > Q 2 ≥ −0.140 GeV 2 . However, SM calculations for FCC-ee processes (low angle Bhabha) will need vacuum polarization for |Q 2 | ≥ 1 GeV 2 .
The path of reduction of uncertainties due to pure QED corrections was also described in ref. [55]. In particular it was pointed out that the dominant O(α 2 L e ) correction not included in BHLUMI 4.04 is already known since long time ago [97] and it can be added to the BHLUMI matrix element rather easily. However, the most elegant solution would be to implement in BHLUMI the QED matrix element with the same type of the CEEX resummation as in KKMC, including in addition the O(α 3 L 3 e ) corrections. According to ref. [55], this would reduce uncertainties due to all photonic corrections to the level below 10 −4 .
The quasi-photonic corrections due to emission of the light fermion pairs are small but difficult to master. Reduction of their uncertainty below 10 −4 level will require a special effort and most likely development of another dedicated MC program. The first step in this direction was done in ref. [101,102].
It is quite likely that at the FCC-ee times one may be in the situation in which the dominant contributions to the luminosity uncertainty ∼ 0.01% will be the experimental one from detector and indirectly the experimental one from the vacuum polarisation (QCD), while the pure QED contribution will be again subdominant.
Summarizing, according to ref. [27], the total uncertainty of the theoretical prediction for the luminosity cross section will be at the 10 −4 level, that is a factor ∼ 10 better that at LEP times. This will require an upgraded matrix element in BHLUMI to O(α 2 ) CEEX level, development of another MC program of comparable class to control technical precision and further reduction (a factor of 2) of the vacuum polarisation error contribution.

Measurement of α QED (M Z ) using A F B near Z resonance
The idea of using precision measurement of the charge asymmetry A F B near Z resonance in order to extract the QED coupling constant at the the electroweak scale α QED (M Z ) was proposed quite recently in ref. [103]. It is an alternative to the current method calculating α QED (M Z ) using as an input very precisely known α QED (0) and the low energy hadronic cross section as an input to dispersion relations. The new method requires charge asymmetry to be measured with the precision δA F B 10 −5 at √ s ± M Z ± 3.5 GeV. According to ref. [103] both statistical and experimental systematic error can be reduced to this fantastic low level, factor ∼ 200 better than what we have seen in LEP experimental results.
Concerning QED uncertainties, it was argued in ref. [103] that because both ISR and IFI corrections do not change sign across the Z peak, they will cancel in the difference However, in the analysis of ref. [33] based on the numerical results of KKMC it was shown that although such a cancellation is present, nevertheless it is not perfect and the remaining net effect in the difference is of order ∼ 1%. On the other hand, quite luckily, the QED effects near the the Z resonance, for a relatively sharp cut-off on the total photon energy E γ ≤ 0.02 E beam , are governed by soft photon emissions. Because of that, the resummation of higher order corrections using the soft photon approximation, which is the basis of the CEEX matrix element of KKMC, is very efficient, not only for the total cross section, but also for the charge asymmetry.
In the study of ref. [33] the results of KKMC for A F B (s ± ) and their difference were validated by means of comparing them with the newly developed MC program KKFoam with partial analytical integration and soft photon resummation. It was shown that predictions of KKMC can be trusted down to the δA F B (s ± ) ∼ 10 −4 precision. Also, some additional internal tests of KKMC indicate that the level of δA F B (s ± ) ∼ 10 −5 is attainable in the results of KKMC for the difference A F B (s + ) − A F B (s − ). The above results have to be consolidated, but they indicate that thanks to enormous power of the soft photon resummation at the amplitude level developed in ref. [16] the improvement of the QED prediction for A F B near the Z resonance by factor 200 needed for the FCC-ee experiments looks feasible. The above analysis concentrates on the pure photonic corrections to A F B and many other corrections like light pair emissions, QCD, electroweak corrections, in this extreme precision regime will have to be analyzed very carefully. In particular it should be checked to what extent all these corrections/effects cancel in the difference Summarizing, the present status of the QED uncertainties (mainly due to IFI) is that they can be controlled in A F B (s + ) − A F B (s − ) and in α QED (M Z ) down to the ∼ 10 −4 level. Another factor 10 improvement will require hard work but according to preliminary study in ref. [16] it looks feasible.

Invisible Z decay width from cross section and radiative return
In the determination of N ν from the peak cross section the error of the luminosity is the main obstacle in the precision improvement. Using the approximate rule δN ν = 7.5 δL L of ref. [29] we estimate that in order to get δN ν 10 −3 the luminosity relative error has to be improved from the present δL/L 5 · 10 −4 down to 1 · 10 −4 . As discussed in ref. [55] it is feasible, provided that the QED matrix element in BHLUMI is upgraded to the O(α 2 ) CEEX level, similarly as in KKMC, and the vacuum polarisation at t = −1GeV 2 is improved by another factor 2. The error of σ(M Z ) due to missing higher orders in the QED ISR calculations (including light fermion pairs emission) was conservatively estimated in ref. [31] to be δσ/σ(M Z ) 2 · 10 −4 . Factor 2 improvement is probably obtainable by means of the careful reexamination of the existing calculations. Reducing experimental systematic error to the same level is of course an independent important issue.
The prospects of reducing the QED uncertainty in the Z radiative return (ZRR) process are discussed in ref. [104]. Here we only summarize the main points.
First of all, in the direct determination of N ν from the ZRR neutrino-like cross section one has to take into account that the QED matrix element of KKMC is much better than that of KORALZ [105] and NUNUGPV98 [64] used at LEP. Nevertheless, due to selecting event with one real photon, one is effectively loosing one perturbative order and it is effectively "downgraded" from to the complete O(α 2 ) exp. to O(α 1 ) exp. . Validating precision of the KKMC for the ZRR process is not yet completed -an optimistic estimate of ref. [104] for the photonic uncertainty of the ZRR process is quite promising, δN ν /N ν 2.4 · 10 −5 for 105 GeV 24 and δN ν /N ν 2.2 · 10 −4 for 161 GeV. The present error due to luminosity δL/L 5 · 10 −4 [29] implies δN ν /N ν = 3 δL L 1.5 · 10 −3 and it would dominate, although its reduction by factor 2-5 according to ref. [55] is feasible. More conservatively, for full exploitation of the FCC-ee precision it would be recommended to develop a new dedicated MC program with the complete O(α 2 ) exp. matrix element for the ZRR neutrinolike process. In the mean time one should validate more precisely predictions of KKMC for ZRR by means of comparison with the semi-analytical programs based on collinear structure functions similar to KKsem [16] or NUNUGPV98 [64].
In the method of extracting N ν from the ratio of the ZRR photon spectrum for the invisible Z decay and Z decaying into muon pairs, the luminosity will cancel out. In addition the effect of the ISR will also cancel quite precisely as well, as shown in ref. [104]. On the other hand, due to the cutoff on the photon angle a sizable QED FSR effect in the muon pair process will remain in the ratio, see [104]. NB. The contribution from the IFI is also sizable, but does not contribute much to integrated cross section because it changes sign across the Z peak in the photon spectrum. Consequently, a high quality MC program for predicting ZRR with muon pairs will be mandatory. KKMC will be obviously helpful, provided its validation using auxiliary programs is pursued, but ultimately a genuine O(α 2 ) exp. for the muonic ZRR is highly desirable.
Summarizing, the optimistic estimate of the photonic QED correction in ref. [104], 24 The energy √ s = 105GeV is not included in the FCC-ee operational mode.
based on the internal crosschecks of KKMC, (dominated by the FSR effect in the muon-pair ZRR) is δN ν 0.9 · 10 −3 , for both c.m.s. energies, 105 GeV and 161 GeV. It would be necessary to develop a new MC program dedicated to the ZRR process in order to crosscheck the above result and gain another factor 3 needed for pushing the uncertainty of QED corrections to a subdominant level.

Cross section near the W W threshold (M W measurement)
Until recently there was no quantitative reliable estimate of the experimental precision of the M W measurement at FCC-ee from the W mass reconstruction using final state quarks and leptons 25 . According to recent studies presented at the FCC Week 2019 [106] both type of measurements, from final state mass reconstruction and from the threshold cross section will provide precision within a similar range of 0.28-0.45 MeV. In the following we shall refer mainly to the simpler measurement of M W at FCC-ee from the W -pair threshold cross section.
The mass of the W boson can be determined very precisely from the value of the total cross section near the threshold of the e + e − → W + W − process, at the c.m.s. energy where the cross section is most sensitive to M W , which is ∼ 162.5 GeV [3]. (Similarly, the best c.m.s. energy for Γ W determination is 157.5 GeV.) Knowing the accuracy of the cross section, the precision of the W mass determined form the threshold scan can be determined approximately from the following relation [107]:  [107] to be 0.91GeV/ √ pb at the W W threshold. The physics goals of FCC-ee are set to 10/ab and 3 × 10 7 events at the W W threshold [108]. For example, for the cross section σ W W 3 pb and the accuracy of the total cross section ∆σ W W σ W W = 1/ √ 3 × 10 7 = 0.02% the resulting accuracy of the W mass would be The ultimate goal of the measurement of the total cross section at the W W threshold at the FCC-ee is the precision of 0.02%. This requires at least a factor 10 improvement of the theory precision for the SM calculation of the W -pair production near the threshold. The first step in this direction is obvious -a MC program with the complete O(α) EW corrections to e + e − → 4f process must be available. This task is definitely doable since the analytical/numerical calculations already exist. That would consolidate the EW precision at the level of 0.2%. In order to reduce the precision further a O(α 2 ) calculation for the doubly resonant e + e − → W + W − process is needed, with some clever strategies inspired by EFT near the threshold.
A complete EW O(α 2 ) calculation for the e + e − → 4f and O(α 3 ) for e + e − → W + W − process seem to be very hard to do, but one may try to estimate uncertainties in case it is not available. Firstly, let us try to estimate an error due to the neglected O(α 3 ) EW corrections. There are no tools to determine their actual size. However, one can do naive scaling based on the lower order corrections. Namely, if at 161 GeV O(α) ∼ 2% of Born (ISR excluded) and O(α 2 ) ≤ 0.2% of Born, i.e. there is a factor of 10 suppression between them, then we can assume the same factor of 10 between O(α 2 ) and O(α 3 ) and we end up at O(α 3 ) ≤ 0.02% of Born, which shows that O(α 3 ) contribution should be just negligible.
Since O(α 2 ) complete calculation for the e + e − → 4f process is a difficult challenge, one should consider another, simpler scheme. This scheme would mimic the approach from LEP2 based on calculating separately the O(α 2 ) corrections to the production and decay of W -pair, i.e. to e + e − → W + W − and then W → 2f process 26 . Similarly as at LEP2, but at one order higher, the above calculations would be combined with the existing O(α 1 ) complete calculation for the e + e − → 4f process. What can we say about precision of such an approach? As earlier, we can only do a naive scaling from the lower order. Namely, analyzing the Single-pole (SP) and Double-pole (DP) contributions at O(α) we find from [77,78] that at 161 GeV we can estimate the O(α) SP as If we now assume, as argued earlier, that O(α 2 ) ≤ 0.2% of Born and we extrapolate from the first order that O(α 2 ) SP ∼ 8% O(α 2 ), then we arrive at O(α 2 ) SP ∼ 0.016%, within the targeted precision of 0.02%. Near the threshold, where the split into DP+SP+NP, using decomposition into powers of Γ W /M W becomes inefficient, one should try to exploit expansion in powers of β, β 2 = (s − 4M 2 W )/4M 2 W , as in EFT calculations of Refs. [85,86], see also the review in Ref. [6], but for the traditional diagrammatic calculation within the standard phase space, which can be implemented within the MC event generator. This may reduce significantly the number of diagrams in the O(α 2 ) for e + e − → 4f to a manageable subset.
In subsection 4.5 it was underlined that the implementation of the QED universal corrections in RACOONWW on one hand and in KORALW and YFSWW3 on the other hand was very different. The immediate question is whether these approaches can be easily extended to new Monte Carlo event generators including complete O(α 2 ) for e + e − → W + W − matrix element or any complete or partly complete O(α 2 ) for e + e − → 4f ? It looks that it will be very hard to extend the methodology of RACOONWW to O(α 2 ) -the subtraction technique of Catani-Seymour is limited to O(α 1 ) and the "sectorisation" technique becomes very complicated beyond O(α 1 ). On the other hand, an example of the MC, KKMC program, with complete O(α 2 ) QED for the e + e − → 2f process exists already, and YFS-inspired CEEX technique [11,16] of factorizing and resumming universal QED correction within the MC event generator will work at the practical level both for e + e − → W + W − and e + e − → 4f processes. In particular, the CEEX technique is very well suited for implementing "nonfactorisable" interferences between the W + W − production process and two W decays and between the two W decays as well. It extends the exponentiation of soft photons of ref. [17]. In the classical approach one exponentiates only interferences and emissions from external legs.
In the new CEEX-style scheme sketched in refs. [110,111] also photon radiation including interferences from the internal charged W bosons is exponentiated. One may ask whether it makes sense as there are no related infrared singularities present. The answer is yes, because W s can be treated as almost stable particles. The main features of the above scheme are the following: • The missing interferences for production-decay and decay-decay to the W W graphs are included to all orders in soft approximation.
• Shift in kinematics due to interferences/emissions (recoil) in W propagators is properly described.
• Higher order corrections can be easily added to the non-soft functions 27 resulting from the IR subtractions.
The above scheme would work for single-W or ZZ pair production as well. The implementation may start with simpler EEX-type scheme, already partly implemented for e + e − → W + W − in YFSWW3 and for W ± decays in WINHAC MC [84], without productiondecay and decay-decay interferences. The CEEX scheme would be implemented in the next step by means of reweighting MC events with multiplicative weight, similarly as it is done in KKMC. The above would take care of the doubly resonant part of e + e − → 4f . Implementation of the CEEX matrix element for the single-resonant part would be straightforward. For the non-resonant part one probably needs only CEEX for ISR and some crude FSR implementation, neglecting IFF interferences. Summarizing, the precision of the W W legacy MC codes from LEP2 is 2% at the threshold, as follows from the direct comparison with the complete calculation [77,78]. The inclusion of this complete O(α) calculation would improve EW precision to the level of 0.2%, as follows from the dominant O(α 2 ) calculation of [85,86] and estimates in [77,78]. To achieve 0.02% EW precision of the cross-section one has to compute and implement the O(α 2 ) EW corrections. We argued that it could be enough to calculate them in double/leading pole approximation, supported by the methods of near-thresholdimprovements of refs. [85,86]. The exponentiation of real and virtual radiation from intermediate W s proposed in ref. [110] can further improve the precision. The QCD corrections have to be analyzed and included as well. This way the overall SM precision tag ∼ 10 −4 for the e + e − → W + W − process seems feasible. The YFSWW+KORALW approach looks like a good starting point for the above developments. The same improvements of the theoretical calculations and MC programs are essential for M W measurements using final state mass reconstruction.

New ideas on pseudo-observables at FCC-ee
The system of electroweak pseudo-observable (EWPOs) used in final analysis of LEP data near Z resonance was forged in the interaction between theorists and experimentalists [29,30]. As already said, it may not work at the higher experimental precision of FCC-ee. The authors of Section C2 of ref. [4] are proposing to modify the EWPO scheme of LEP to the FCC-ee level. The main aim is to preserve its fundamental feature, that is to provide a flexible link between data and theory. EWPOs should provide model independent encapsulation (representation) of experimental data, where model independence means removing from data technical details of the detector (inefficiencies), kinematic cut-offs and large universal QED corrections.
Let us explain this new modified scheme of EWPOs following its graphical representation in Fig. 3. Similarly as in the LEP case, the EWPOs are encoded in the parameters of some effective Born (couplings, masses) such that EWPOs like cross sections, asymmetries, partial widths are in one-to-one correspondence with these parameters, which are called EW pseudo-parameters (EWPPs). The critical question is whether one may factorize-off and remove QED corrections at the FCC-ee precision level on the way from raw experimental data to EWPOs/EWPPs, (A) → (C) in Fig. 3?
As indicated in the previous section, the QED deconvolution of LEP, which was done at the amplitude squared level and ignoring non-factorisable IFI effects, will be unable to meet FCC-ee precision criteria. However, the QED factorisation at the amplitude level has the extra power for the task. The price to pay is that it cannot be done analytically, but has to be implemented using numerical sesummation within the Monte Carlo. The corresponding transformation of data (B) → (C) in Fig. 3 was done in the LEP scenario using ZFITTER/TOPAZ0 programs and in the proposed scheme the fitting parameters in the effective Born would be done by the Monte Carlo programs. An additional advantage is the full flexibility concerning the choice of experimental cut-offs at the stages (A) or (B) in Fig. 3. In addition, the MC program will have not only effective Born amplitudes convoluted with multiple photon emissions, but also O(α 1 ) or even O(α 2 ) non-IR electroweak corrections, thus it will be able to fit directly data to SM predictions, as indicated in (B) → (D) in Fig. 3. This is important, because what we want for new EWPOs/EWPPs is possibility of fitting data to SM, (C) → (D) in Fig. 3, as in case of LEP scenario. Direct fit of data to the SM predictions, (B) → (D), will provide a crosscheck on the precision loss in the desired two-step scenario (B) → (C) → (D). Another powerful crosscheck done in ref. [30] was to close the loop (B) → (C) → (D) → (B). This kind of cross-checks were also done in ref. [29] in the final LEP data analysis. Obviously, the role of the MC in the proposed scheme is more important than at LEP scenario. It will also require implementing in the MC programs certain provisions for fast evaluation of the changes of the SM predictions due to small variation of the SM input parameters, using the technique of weight differences.
It should be stressed that scenario proposed above is one of many possibilitities. It will be the aim of the future studies and practice of the FCC-ee data analysis to decide about its final shape. For instance, in the above LEP-like scenario the main objective in the (B) → (C) transition is to encapsulate experimental data in the form of Born like parametrization of the vector couplings of the Z exchange, while background nonresonant exchange is parametrized in a rather simple way as photon exchange. One may consider departing from the above by means of including in the simplified spin amplitudes used in (B) → (C) stage most of the O(α 1 ) EW corrections (like EW boxes), such that they are removed from data. In this way EWPOs/EWPPs would represent more cleanly possible signals of new physics. One may also think about extending list of EWPPs with parameters residing in the non-resonant part of spin amplitudes. This would help to describe them in a refined way, again with the aim of exposing more efficiently possible signals of new physics in experimental data.
The above analysis is only indicating certain possible avenues of reformulating the methodology of electroweak pseudo-observables and more detailed analysis is badly needed. In the end, experimentalists in close collaboration with theorists will decide on how FCCee data will be analyzed.  Table 2: Comparing experimental and theoretical errors at LEP and FCC-ee as in Table 1. 3rd column shows LEP experimental error together with uncertainty induced by QED and 4th column shows anticipated FCC-ee experimental statistical [systematic] errors. Additional factor ×3 in the 5-th column (4th in Table 1) reflects what is needed for QED effects to be subdominant. Rating from to marks whether the needed improvement is relatively straightforward, difficult or very difficult to achieve.

Summary
The main results of our study are indicated in Table 2, where we have indicated for selected observables, the same as in Table 1, the improvement factor needed in calculations of the QED effects in order to match the experimental precision anticipated in the FCC-ee experiment. We have also indicated explicitly the additional factor 3 necessary for these effects to become subdominant and the most important development to be done. We have also tried to rate how difficult it will be to achieve these targets. It is needless to say that many of the above estimates remain speculative and are not based on solid numerical results. However, this is always the case with estimates of the uncalculated higher order perturbative corrections, so they have to be always taken with a grain of salt. On the other hand, this kind of analysis is indispensable in planning directions and priorities of the future work.