The method to genrate the particles was simple - I used a von Nemann method to generate particle distributions conforming to the weight function formula (8) from the mentioned paper. I used the ROOT macro to do this.

I present the first results I obtained from the simulation:

This is a Trnsverse radius distribution for three different transverse momentum regions. (0 - 0.3, 0.3 - 0.5, 0.5 - 0.8 GeV)

I've made these plots for other ALFA values as well:

• ALFA=10.0
• ALFA=100.0

The last plot I've made was a p_T distribution with the dependence on ALFA value:

Further studies have to be made to fully understand these plots, and decide if they are correct.