Sources

[Emmanuel Farhi (71 83)]

Hello,

Yes, you're right, creating a rectangular source is very easy (just change angle and radius directly into dimentions for the two first random numbers). For a Maxwellian distribution, I hesitate between using a flat energy distribution, but weightening with p=cte.exp(-E/T), or using E = -T log(rand01/cte), in order to produce more neutrons of certain energies. The first solution is simpler, but some neutrons with low weight are not very useful. The second solution appears to me more efficient in computation.
What do you think about it ?

Cheers. E.

--------------------------------------------
Emmanuel FARHI, http://www.ldv.univ-montp2.fr:7082/~manuf   \|/ ____ \|/
TAS-Group, Institut Laue-Langevin (ILL) Grenoble            ~@-/ oO \-@~
Avenue des Martyrs, BP 156, 38042 Grenoble Cedex 9,France   /_( \__/ )_\
Work :Tel (33/0) 4 76 20 71 83. Fax (33/0) 4 76 48 39 06       \__U_/
La Grande Arche, Chateau d'Uriage, 38410 Saint Martin d'Uriage 04 76 59 73 94