  Merging the benefits of an NLO calculation with event generators is a topic
of current interest. Among other things, it is advantageous to be able to
generate unweighted events according to NLO matrix elements. I report on an
algorithm which generates a sequence of unweighted momentum configurations,
each configuration containing either n or n+1 four-vectors, such that for any
infrared safe observable the average over these configurations coincides with
the NLO calculation up to errors of a resolution variable.
