|Message: Cos-squared theta distribution using GPS||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I am trying to simulate cosmic ray neutrons but my scoring result does not match the angular spectrum I set up.
I start with assigning a cos-squared theta distribution to the neutron source in vacuum using macro file as following:
## Define the shape of the source
/gps/pos/rot1 -1 0 0
/gps/pos/rot2 0 1 0
/gps/pos/halfx 0.0001 m
/gps/pos/halfy 0.0001 m
/gps/pos/centre 0 0.35 0 m
## Set up the histogram for theta using "/gps/hist/point theta(radian) cos^2"
/gps/hist/point 0.000 1.0000
/gps/hist/point 0.017 0.9997
/gps/hist/point 0.035 0.9988
/gps/hist/point 0.052 0.9973
(omit the rest of the histogram, but I have it to go all the way to theta = 1.571 radian (90 degrees) )
And I set up a scorer following this tutorial http://geant4.slac.stanford.edu/tutorial/mcgill06/HandsOn3/HandsOn3.htm#_Installation
such that I have a hemisphere centered at the source with a radius of 0.5 meter. There are 90 rings scoring the incoming sphere surface current on the hemisphere, so each corresponds to the angle theta with an increment of 1 degree. However, the scoring result does not match the cosine-squared distribution. Here is part of the result of 1e4 neutrons on target:
Theta Scorer 1
Can anyone help me out? Any ideas and advice will be greatly appreciated. Thank you!
|Inline Depth:||Outline Depth:||Add message:|