A method and a system to produce, either in numerical simulations or in experiments, specified amounts of first, second and higher order PMD in a controlled manner, in particular large amounts. Parameters can be adjusted to obtain specific ranges of first, second and higher order PMD, and importance sampling can be used to determine the probability that the resulting PMD events can be obtained in realistic situations. Individual results obtained using specific parameter values can be combined to produce even larger ranges of PMD.