Toy MC for Z->mumu


I am trying to generate a toy model for plotting the pT distribution for Z->mu^{+}mu^{-} events. I see that there are tutorials on generating toy model datasets, but I am a little unsure on how to produce the aforementioned dataset. Does anyone have any ideas or references to direct me toward? Thanks.