Fluctuations of muon energy loss and simulation of ionization cooling

A Monte Carlo program to simulate muon ionization cooling is outlined. A Vavilov-type distribution to represent restricted ionization energy losses is derived. Above the restriction threshold μe scattering is treated event-by-event. Likewise μ-nucleus elastic scattering is simulated by a Gaussian below some angular threshold and treated individually above it. Other processes included are: incoherent μp scattering with nuclear protons, bremsstrahlung, pair production, and deep inelastic nuclear scattering. A small sample of results obtained with the code is included.