  To calculate dilepton rates in a Monte Carlo simulation of ultrarelativistic
heavy ion collisions, one usually scales the number of similar QCD processes by
a ratio of the corresponding differential probabilities. We derive the formula
for such a ratio especially for dilepton bremsstrahlung processes. We also
discuss the non-triviality of including higher order corrections to direct
Drell-Yan process. The resultant mass spectra from our Monte Carlo simulation
are consistent with the semi-analytical calculation using dilepton
fragmentation functions.
