Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] acej_v1_0.gz(13 Kbytes)
Manuscript Title: Monte Carlo simulation of radiative corrections to the processes e+e- - > mu+mu- and e+e- -> q(bar)q in the Z0 region.
Authors: F.A. Berends, R. Kleiss, S. Jadach
Program title: MUSTRAAL
Catalogue identifier: ACEJ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 29(1983)185
Programming language: Fortran.
Computer: CDC CYBER 173.
Operating system: NOS/BE.
RAM: 12K words
Word size: 60
Keywords: Elementary, Particle physics, Quantum electrodynamics, Radiative corrections, Monte carlo simulation, Muon, Quark, Photon, Electroweak theory, Bremsstrahlung.
Classification: 11.4.

Nature of problem:
Radiative corrections to a process are experiment dependent. A Monte Carlo simulation of both mu+mu- and mu+mu-gamma events is ideally suited for imposing any experimental cuts. The MC program presented here improves on F.A. Berends and R. Kleiss, Nucl. Phys. B177(1981)237 in two respects: it is faster and it includes Z0 production in the GWS model.

Solution method:
The program decides to generate either a soft or a hard bremsstrahlung event. In both cases events are first generated according to an approximate distribution. By weighing procedure this is transformed into the exact distribution.

The virtual corrections take only the QED part into account. Moreover, terms of order mmu/E (muon mass over beam energy) are neglected, which restricts the applicability at low energies. With this restriction the program can also be used for heavy charged leptons and after some changes in coupling constants for quarks.

Running time:
160 events per CPU s.