Two methods of folding a distribution of simulated muon bundle events underground with a finite size detector are compared. A general strategy to fully exploit the statistics produced by the event generator is described. The aim is to minimize the amount of computing time required in the analysis of data taken by large area detectors. The method presented here can also be used for other applications typical of cosmic ray physics.