We present the matching between a next-to-leading order computation in QCD and the PYTHIA parton shower Monte Carlo, according to the MC@NLO formalism. We study the case of initial-state radiation, and consider in particular single vector boson hadroproduction.