arXiv:hep-ph/9407308v3 22 Jul 1994 b Shimizu Y. Kurihara, Y. a Kawabata, S. Ishikawa, T. processe Fujimoto, J. four-fermion radiative in diagrams Non-resonant ( e oe ut ag.Tennrsnn diagrams non-resonant The by dominated are large. be- contribution quite of relative comes their off-shellness ones, the resonant with the only di- competing non-resonant agrams the threshold, Below of subset diagrams. with the exper- in cope appearing to the of violations applied gauge was on set electron cut a outgoing the a to of constraints, addition loose in imental results, these In in states final four-fermion two Introduction 1. the for diagrams non-resonant the of contribution relative h rs etosfor sections cross the ffc thge nry(%at (7% energy higher at effect atrl eo h artrsod(pt 27% to (up threshold pair impor- W shown at an the was below play It role diagrams tant [1]. non-resonant meeting the Sotchi that the at ported E Ieeg ag n o uuelna oldr ( colliders contribu linear Their future for studies. and range other energy in Specia II far LEP presented. so are overlooked case diagrams radiative and non-radiative the lrt ihtewl-nw w-htn( two-photon well-known the with ilarity o h aitv aecnb eno i.b2 As Fig.1b,2. diagrams on seen of be type can case these produc- radiative of W the for no examples or Some W tion. single essentially dedicated are a pro- cesses non-resonant on The based [2]. states generator four-fermion final selected anal- the of by ysis workshop This this at energy. confirmed the was effect with res- falls of which that diagrams to onant non-resonant contrast in of growing contribution keeps processes the with and increases energy section the cross the Hence processes. e + API23CR,France LAPP-IN2P3/CNRS, iaiTty olbrto,KK Japan KEK, Collaboration, Minami-Tateya − h opeete ee rs eto for section cross level tree complete The w er g,acmlt aclto of calculation complete a ago, years Two e √ ν ¯ − e s u d → ¯ 5 e)adehbtanon-negligible a exhibit and GeV) 150 = n 7 (¯ 27% and ) e − ν ¯ e u d ¯ t udu -channel and e + d ¯ e r eotd iia eaiu o h aitv aei s is case radiative the for behaviour Similar reported. are ) e − + → e − γ e − − → ν ¯ W √ e e u s + udu ¯ d rpsi sim- in graphs ¯ e 9 GeV). 190 = and − d ¯ a re- was , collisions, e θ + e γ γ + e angle − − e − − → W γ → √ ) udu ¯ s e 50GV.Eet,a h artrsod fodr3% order of threshold, pair W the at Effects, GeV). =500 − d ¯ ν ¯ vn eeaosbsdo h RC akg for package GRACE the on based generators Event . aitv hne ece 42.5%. reaches channel radiative e neeti ruh oteeeto h non-resonant the of effect the to brought is interest l u o etn h aiiyo ipie n therefore and generators. simplified faster benchmark of validity a as the used testing be for can four-fermion work radiative This the level processes. we of tree one [4][5], complete of the including time, computation years, first the last for these present, in out ried uncer- tainties. additional introduce state will final W, as bremsstrahlung, recognized intermediate mistakingly be the How- can from which algorithm. coming accounted reconstruction photons mass properly ever, the and in detected for is unless photon reconstruction, mass the additional the to on contribute uncertainties ra- will state corrections contri- Final diative threshold, relative diagrams. pair non-resonant the of W enhance bution the therefore, to will, Close ISR energy. center mass available the of reduce (ISR) radiations state inhv ob nlddfrapeiecomparison precise the a at for measurements included with be to have in tion corrections radiative the that uncertainties. overall measurement the width to and contribute mass will concerned, background they as and is counted processes, reconstruction be should mass diagrams W these the as far dγ int h oa rs eto speetdfrthe for presented is section cross total the to tion ¯ lhuhanme fsuishsbe car- been has studies of number a Although thsbe eonzd ic oetm [3], time some since recognized, been has It scmue n icse ncmaio with comparison in discussed and computed is a n .Perret-Gallix D. and on At hown. s e + e − √ oldr.Initial colliders. b s W 0 e,the GeV, 500 = arproduc- pair 1 2

γ e γ u e u e W e d W d u u ν W γ e W νe d d γ W ν ν e ν e e e ν Z W e Z e e Z e γ e e e e e W (a) Figure 2. (Fig.1 cont.) Some examples of + − γ ν diagrams from the radiative process e e e e − ν ¯ → e e u e ν¯eudγ: Other non-resonant diagrams. d W W d W the denominator can vanishes for positive squared u u γ γ γ momentum transfer. A gauge invariance checking e e e e program is automatically built by the system. Only the kinematics for complex processes has to be written by the user. The issue, here, is to ν e νe e e figure out what is the best transform to apply to the basic integration variables to regularize the u u W W divergencies appearing in the integration of the d γ d γ d matrix element. The best kinematics will give e e the best accuracy for a given computation time. e e γ γ (b) 3. Four-fermion production The results presented hereafter have been ob- Figure 1. Some examples of diagrams from the + − − tained using the following set of parameters: radiative process e e e ν¯eudγ¯ : a) Resonant → diagrams, b)Non-resonant γ W diagrams. MZ = 91.1 GeV − ΓZ = 2.534 GeV 2. The GRACE system α = 1/137 2 2 2 sin θW = 1 (M /M ) The GRACE system [6] has been used to − W Z perform the very lengthy computations involved MW = 80 GeV in the study of most radiative processes. The mu = md = 0.1 GeV. GRACE package is a complete set of tools for computing tree level processes. All the phases The width of the W is taken from the Particle involved in a given computation are covered: Data Table; ΓW = 2.25 GeV. The gauge boson from the process definition (specification of the (W,Z) widths are assumed to be constant in the initial and final particles) to the event genera- calculation. Furthermore some realistic experi- tor. It is composed essentially of three compo- mental cuts have been introduced: nents: the diagram generator, the construction of ◦ θe > 8, 20, 30, 40 the matrix element based on helicity amplitude ◦ ◦ function from the CHANEL [7] library and the 172 > θγ, θq,q¯ > 8 multi-dimensional phase space integration pack- Eγ > 1 GeV, Ee,q,q¯ > 1 GeV age BASES [8] associated with the event genera- tor SPRING [8]. Fermion masses are properly introduced in the where θe,θγ and θq,q¯ are angles measured from − helicity amplitudes. The particle width is intro- the incident e beam and Ee,q,q¯ are the energies duced into the gauge boson propagators when of final e−, q andq ¯, respectively. The electron 3 polar angle cut is introduced to avoid an undesir- able gauge violation for extremely forward direc- tions [9] caused inevitably by the finite width of W boson when the subset of γ W diagrams is considered. − The gauge invariance of the amplitude is checked numerically for a random selection of the boson gauge parameters at several points in the phase space. The errors are within the preci- sion of numerical calculation (typically less than O(10−12) in double precision) when the W width is turned off. Further gauge cancellations are checked to confirm that the obtained cross sec- tions are stable against finite width. The pro- cedure is as follows: first extract the product of gauge boson propagators without width as an overall factor of the whole amplitude. Then re- place them by those with finite width. As a result, one is left with rigorously gauge invariant ampli- tude[10]. For the resonant diagrams, this method gives the same results as ordinally method. While for the γ W diagrams, this gives 10% smaller + − − ¯ − Figure 3. Cross sections for a) e e e ν¯eud results. and b) e+e− udu¯ d¯ versus C.M. energy→ with ◦ → θe 8 . The solid curves are for all diagrams 3.1. The non-radiative case and≥ dashed ones for three resonant graphs. The The contribution of the non-resonant diagrams dotted line in a) is the contribution from the γ − ¯ has been considered for two final states, e ν¯eud W and Z W processes. The estimates at 500− ¯ andudu ¯ d. Fig.3 shows the cross sections ver- GeV are indicated− by arrows. sus the center of mass energy around the W pair threshold. The solid line represents the case where all diagrams are included, the dashed line polar angle cut. The non-resonant contribution the contribution of the resonant diagrams only. vanishes as the cut is raised. In theudu ¯ d¯ channel, one can observe the onset of the Z resonant diagrams as a small shoulder 3.2. The radiative case starting around 2MZ . The behaviour below the In the covariant gauge, the GRACE package W pair threshold can be better seen in fig.4 where generates 752 diagrams when all the interactions the ratio of the three main resonant diagrams over are included, but actually 142 diagrams are used the complete calculation is plotted for the same in the unitary gauge with Higgs-fermion coupling − range of center of mass energy. In the e ν¯eud¯ being set to 0. The full tree level computation channel, the effect is of the order of 3% around requires, for a decent precision, some 100 hours √s = 160 GeV, while inudu ¯ d¯ it reaches± 28%. on a HP 735 computer. Thanks to the vectoriza- At the energy of future LC, 500 GeV, the− contri- tion of the integration package BASES, the most bution of non-resonant diagrams is quite signifi- extensive computation could be performed on the − cant in e ν¯eud¯ but is negligible inudu ¯ d¯. KEK vector processor within 1 hour. Fig.6 shows Similarly to the γ γ events, the non-resonant the contribution of the non-resonant terms in the processes are predominantly− peaked in the for- total cross section versus the C.M. energy. The ward/backward direction. Fig.5 shows the be- behaviour is similar to the non-radiative case. In haviour of the cross section versus the electron fig.7, one can see more clearly the contribution 4

Figure 4. Energy dependence of the relative con- Figure 5. Electron angular cut dependence for + − + − − tribution of resonant diagrams in a) e e e e e ν¯eud¯ (C.M. energy = 180 GeV). Solid − + − ◦ → → e ν¯eud¯ and b) e e udu¯ d¯ with θe 8 to line: all diagrams, dashed line: three resonant the full computation. → ≥ diagrams and dotted line: γ W and Z W processes. − − of non-resonant diagrams similar to fig.4 for non- − radiative case. two non-radiative computations e ν¯eud¯andudu ¯ d¯ The dependence of the cross section on the elec- has been carried out. The contribution of the tron polar angle cut θe is presented in fig.8. Fig.9 non-resonant processes is large below threshold shows the distribution of the invariant mass of the (at least forudu ¯ d¯ channel) which makes the W pairs of final state fermions. The general shape is, mass determination from the excitation curve as expected, a Breit-Wigner distribution. How- more difficult. At high energy, the contribu- − ever some background coming from the onset of tion increases, reaching 46% (e ν¯eud¯) and 42.5% − non-resonant processes can be seen. In addition, (e ν¯eudγ¯ ) at √s = 500 GeV. Determination of for the radiative process (fig.10), a mass shift of the W mass from direct reconstruction will have − the order of 1 GeV is observed on the e ν¯e in- to deal with the intrinsic background from the variant mass due to final state radiation. non-resonant processes prominently in the for- ward direction. Furthermore the final state radia- 4. Conclusions tive correction will introduce a shift in the mass of the W unless the photon is properly taken into A complete calculation of the radiative four- account in the algorithm of mass reconstruction. + − − fermion final state e e e ν¯eudγ¯ has been However, this last point may prove to be diffi- performed, for the first time,→ using automatic cult as the photon may be produced either from Feynman diagram calculation techniques imple- the initial state, the intermediate charged boson mented in the GRACE system. Comparison with or the final state fermions. Only this last class 5

+ − − Figure 6. Total cross section for e e e ν¯eudγ¯ Figure 7. Energy dependence of the contribution ◦→ + − − versus the C.M. energy with θe 8 . The solid of non-resonant diagrams in e e e ν¯eudγ¯ ≥ ◦ → line represents the complete cross section while with θe 8 . The ratio of contribution of reso- the dotted line the resonant term only. nant diagrams≥ versus all diagrams is shown. plays a role in the W mass determination. More la Recherche Scientifique (CNRS) and l’Institut detailed studies are needed to check if additional de Physique Nucleaire et Physique des Particules kinematical cuts can improve the mass determi- (IN2P3) in France. nation. However to quantify the actual relative contri- REFERENCES bution of radiative events in the total cross sec- tion, soft photon and one-loop contributions must 1. T. Ishikawa, T. Kaneko, S. Kawabata, Y. be taken into account. Large negative contribu- Kurihara, Y. Shimizu and H. Tanaka, KEK tions are expected, reducing, consequently, the Preprint 92-210, 1992 and in the Proceedings number of radiative events and softening the mass of the 7th workshop on high energy physics reconstruction problem. and quantum field theory, Sotchi, Russia, 1992. 5. Acknowledgements 2. M. Pittau, in these proceedings. 3. M. Lemoine and M. Veltman, Nucl. Phys. This work has taken place in the framework B164 (1980) 445 of the KEK-LAPP collaboration supported in 4. H. Tanaka, T. Kaneko and Y. Shimizu, Com- part by the Ministry of Education, Science and put. Phys. Commun. 64(1991) 149. Culture (Monbusho) under the Grant-in-Aid for 5. See related contributions in these proceed- International Scientific Research Program No. ings. 04044158 in Japan and by le Centre National de 6. T. Ishikawa, T. Kaneko, K. Kato, S. Kawa- 6

+ − − Figure 8. Total cross section for e e e ν¯eudγ¯ versus the electron polar angle cuts→ in degree. The solid line: all diagrams, the dashed line: reso- nant diagrams, and the dotted line: resonant and γ W . − bata, Y. Shimizu and H. Tanaka, KEK Re- port 92-19, 1993, The Grace manual Ver. 1 7. H. Tanaka, Comput. Phys. Commun. 58 (1990) 153 and see also [4]. 8. S. Kawabata, Comput. Phys. Commun. 41(1986) 127. 9. E.E. Boos et al., Phys. Lett. B326(1994) 190. Figure 9. Invariant mass for non-radiative process 10. A.E. Pukhov and V.A. Ilyin, private commu- in solid line for the ud¯ system and in dotted line nication. − ◦ for the e ν¯e system with θe 8 . ≥ 7

Figure 10. Invariant mass for radiative process in solid line for the ud¯ system and in dotted line for − ◦ the e ν¯e system with θe 8 . ≥