A computer simulation of four charm production and fragmentation in e+e- collisions

Theoretical predictions of the rate of exclusive production of two charmonia in e(+)e(-) annihilation differ significantly from measurements by the Belle collaboration. Up to now, there has been no conclusive explanation for this substantial discrepancy. However, recent theoretical works have suggested that it may be possible to identify the origin of the discrepancy from a study of four charm inclusive events. We present necessary computational steps required for such an experimental study at B factories. Event generation of the reaction e(+)e(-) -> c (c) over barc (c) over bar and the fragmentation of charm quarks are simulated. These two important computational steps provide all the information necessary to measure the total cross section of the four charm hadron final state at B factories.