Subtraction method for NLO corrections in Monte-Carlo event generators for Z boson production

We use a subtraction method to construct NLO corrections in a Monte-Carlo event generator for the case of vector boson production in Drell-Yan processes. Our calculations are carried out both for the Bengtsson-Sjostrand-van Zijl (BSZ) algorithm and for a modified algorithm proposed by Collins. In the case of the modified algorithm, we compute the relation between the parton distribution functions and the ones in the MSbar scheme; this relation is the same as the corresponding relation for DIS. For the BSZ algorithm, we show that there is no simple relation.