The paper proposes a method to match next-to-leading order (NLO) QCD calculations with parton shower Monte Carlo (MC) simulations. The method ensures that fully exclusive events are generated, with hadronization following the MC model, and total exclusive rates are accurate to NLO. It handles hard emissions as in NLO computations while soft/collinear emissions are managed by the MC simulation with the same logarithmic accuracy. The method also addresses double-counting issues and allows for unweighted event generation. The approach is illustrated using a toy model and applied to the hadroproduction of W+W− pairs in QCD. The results show that the MC@NLO method reduces numerical dependence on subleading logarithms and reproduces the shapes of MC curves in soft-dominated regions.The paper proposes a method to match next-to-leading order (NLO) QCD calculations with parton shower Monte Carlo (MC) simulations. The method ensures that fully exclusive events are generated, with hadronization following the MC model, and total exclusive rates are accurate to NLO. It handles hard emissions as in NLO computations while soft/collinear emissions are managed by the MC simulation with the same logarithmic accuracy. The method also addresses double-counting issues and allows for unweighted event generation. The approach is illustrated using a toy model and applied to the hadroproduction of W+W− pairs in QCD. The results show that the MC@NLO method reduces numerical dependence on subleading logarithms and reproduces the shapes of MC curves in soft-dominated regions.