Application of Beam-beam Interaction to a Particle Density Function

In conjunction with the Fokker-Planck propagator method, beam-beam 
interaction is applied to a particle density function.  Combined with
the propagators for the rest of the storage ring, we are able to
investigate the tail of the particle density function without unphysical run-to-run
fluctuation intrinsic to macro-particle Monte Carlo simulations.