|Message: Dose Normalization||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'm having some difficulties when trying to normalize my outputs, like the Total Dose.
First I'm running with a simple geometry: a spherical target sorrounded by a spherical shell with 2mm thickness.
To compute the dose, I'm using the primitive scorer G4PSDoseDeposit and using GPS as the source.
I took a look at a ppt document of G. Santin about normalization, and as I understood, the "real world" Dose is calculated as:
Dr = Nr*Ds/Ns,
where Nr is the number of particles in the real world, Ds is the dose calculated in the simulation and Ns is the number of events simulated, right?
Now, the "real world" number of particles Nr is calculated by the flux. That's where I'm a little confused. I do have the fluxes as functions of the energy, in a table calculated by Spenvis and the respective macro created. As I understad, the number Nr is calculated by:
Nr = flux * 4*pi^2*r^2.
I don't know which flux to use, since I have a table based on energies. And the radius is defined by me on the source, right? I defined the source like this:
/gps/hist/point 1.000000E-01 2.407800E+01
/gps/hist/point 1.500000E-01 2.397600E+01
/gps/hist/point 2.000000E-01 1.974300E+01
/gps/hist/point 3.000000E-01 1.135400E+01
/gps/hist/point 4.000000E-01 1.059100E+01
/gps/hist/point 5.000000E-01 9.840100E+00
/gps/hist/point 6.000000E-01 1.154300E+01
/gps/hist/point 7.000000E-01 1.277300E+01
/gps/hist/point 1.000000E+00 8.546000E+00
/gps/hist/point 1.500000E+00 3.924800E+00
/gps/hist/point 2.000000E+00 3.838800E+00
/gps/hist/point 3.000000E+00 3.667500E+00
/gps/hist/point 4.000000E+00 3.831800E+00
/gps/hist/point 5.000000E+00 3.993600E+00
/gps/hist/point 6.000000E+00 4.341400E+00
/gps/hist/point 7.000000E+00 4.587900E+00
/gps/hist/point 1.000000E+01 3.820900E+00
/gps/hist/point 1.500000E+01 2.932800E+00
/gps/hist/point 2.000000E+01 2.538000E+00
/gps/hist/point 3.000000E+01 2.066400E+00
/gps/hist/point 4.000000E+01 1.757500E+00
/gps/hist/point 5.000000E+01 1.592800E+00
/gps/hist/point 6.000000E+01 1.462600E+00
/gps/hist/point 7.000000E+01 1.370200E+00
/gps/hist/point 1.000000E+02 1.118000E+00
/gps/hist/point 1.500000E+02 7.635400E-01
/gps/hist/point 2.000000E+02 4.962900E-01
/gps/hist/point 3.000000E+02 2.196200E-01
/gps/hist/point 4.000000E+02 2.586700E-02
/gps/pos/centre 0. 0. 0. mm
/gps/pos/radius 9 mm
/gps/ang/mintheta 0.000E+00 deg
/gps/ang/maxtheta 9.000E+01 deg
So, my radius will be 9mm.
Finally, my normalized dose changes values when I change the number of primary events...this sounds incorret, since I'm normalizing it.
I know I'm doing it wrong, just don't know where.
Please help? I'm really having a hard time on this.
Thank you in advance.
Best Regards, Leonardo.