Uniform particle generation in solid angle [message #11678] |
Mon, 18 April 2011 09:07 |
|
Dear all,
I was wondering if there's an easy way to generate particles along a unit 3-vector. As I noticed, the box generator in emc_complete.c doesn't do this. Apparently, it draws a random number within the specified ranges of theta and phi. What I need is something like G4RandomDirection function. I hit the problem when I roughly calculated the solid angle coverage of the full EMC to be about 92%. Whereas, the efficiency calculation using the box generator gives me something around 85% which indicates a bigger probability for particle generation at the poles (entrance and exit holes of the EMC)! Am I missing something here or there is other alternative to the box generator?
Best regards
|
|
|