<div dir="ltr">Hi Niels<div><br></div><div>One way, or at least the way that I used to learn, was: create a voxel and attach it to an specific voxel of your phantom, the active the option for energy loss:</div><div><h3 style="color:black;background-image:none;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;margin:0.3em 0px 0px;overflow:hidden;padding-top:0.5em;padding-bottom:0px;border-bottom:none;font-size:1.2em;line-height:1.6;font-family:sans-serif"><span class="gmail-mw-headline" id="gmail-Energy_spectrum">Energy spectrum</span></h3><p style="margin:0.5em 0px;line-height:inherit;color:rgb(37,37,37);font-family:sans-serif;font-size:14px">This actor builds four histograms: the initial kinetic energy of each track (energySpectrum), the energy deposition per event (edepHisto), the energy deposition per track (edepTrackHisto), the energy loss per track (eLossHisto). These histograms are stored in a root file. They takes into account the weight of particles. They are three commands to define the boundaries and the binning of the energy spectrum and three commands to define the boundaries and the binning of the energy loss histograms (edepHisto, edepTrackHisto, eLossHisto).</p><pre style="font-family:monospace,Courier;color:black;background-color:rgb(249,249,249);border:1px solid rgb(221,221,221);padding:1em;white-space:pre-wrap;line-height:1.3em;font-size:14px">/gate/actor/addActor  EnergySpectrumActor  <b>MyActor</b>
/gate/actor/<b>MyActor</b>/save  MyOutputFile.root
/gate/actor/<b>MyActor</b>/energySpectrum/setEmin 0 eV
/gate/actor/<b>MyActor</b>/energySpectrum/setEmax 10 GeV
/gate/actor/<b>MyActor</b>/energySpectrum/setNumberOfBins 200
/gate/actor/<b>MyActor</b>/energyLossHisto/setEmin 0 eV
/gate/actor/<b>MyActor</b>/energyLossHisto/setEmax 15 MeV
/gate/actor/<b>MyActor</b>/energyLossHisto/setNumberOfBins 120</pre></div><div><br></div><div>You will get an energy deposit spectrum in this specif voxel, I did the energy conversion and considering the voxel volume and material density, manually I calculated the dose, then I get a factor for my dose image.</div><div>The reason for the factor is due to emission probabilities of the source and other issue like time and activity. If you are already fixing the activity you will get in your spectrum of deposit energy in one line emission probability and in other energy deposited for time unit.</div><div><br></div><div><b>For example for this par data:</b></div><div><b>Energy (MeV)    Probability</b></div><div><b>0,001                  0,550919</b></div><div><b><br></b></div><div><b>To get the Dose:</b></div><div><b><br></b></div><div><b>Dose[Gy-MBq/s] =  0,550919*0,001[MeV]*10E6[eV/MeV]*1,602176487E-19[J/eV]*(1/m[kg])*1E6[MBq/Bq]*3,826[particles/Bq-s]</b></div><div><br></div><div>This last number correspond to the emission probability of all the emission of the radionuclide (gammas+betas+IC electrons+Auger Electron).<br></div><div>In your case if you use only photon annihilation the probability is 1,935, but if your source include also a positron espectrum the factor is 2,908. <b>Please check here in table 5.2:</b></div><div><b><a href="http://www.nucleide.org/DDEP_WG/Nuclides/F-18_tables.pdf">http://www.nucleide.org/DDEP_WG/Nuclides/F-18_tables.pdf</a></b><br></div><div><br></div><div><i><b><u><font size="4">BUT, if are using GATE source take into account:</font></u></b></i></div><div><pre style="font-family:monospace,Courier;color:black;background-color:rgb(249,249,249);border:1px solid rgb(221,221,221);padding:1em;white-space:pre-wrap;line-height:1.3em;font-size:14px"><i>/gate/source/NAME/gps/particle e+
/gate/source/NAME/gps/energytype Fluor18
/gate/source/NAME/setForcedUnstableFlag true
/gate/source/NAME/setForcedHalfLife 6586 s 
</i></pre><p style="margin:0.5em 0px;line-height:inherit;color:rgb(37,37,37);font-family:sans-serif;font-size:14px"><i>In the above example, more properties of fluorine-18 have been added by using the helper keyword Fluor18 (see below): half-life and energy distribution.<b> Note that the branching ratios are not respected with this kind of simulation because no decay is simulated, i.e. emission is 100% positron. </b>You may want to lower the activity by an appropriate factor to take this better point into account.</i></p></div><div><br></div><div><b>And sincerely I don´t how it work, I construct my sources complete with MIRD book data (<a href="https://books.google.de/books/about/MIRD.html?id=zEyo06YS5NMC&redir_esc=y">https://books.google.de/books/about/MIRD.html?id=zEyo06YS5NMC&redir_esc=y</a>)</b></div><div><br></div><div>A last point if you already are fixing the activity, I am not sure but I would expect that your calculation will give you Dose in [Gy/s], then if you include a desintegration curve in your calculation you will get directly Dose in [Gy].</div><div><br></div><div>Sorry, it was pretty long but it is how I deal with that.<br></div><div><br></div><div>Eventually you will fine an easy way, but it is the larger but with better explanation.</div><div><br></div><div>Regards</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-08-09 10:24 GMT+02:00 Niels van der Werf <span dir="ltr"><<a href="mailto:nrvdwerf@gmail.com" target="_blank">nrvdwerf@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear all,<br>
<br>
Last week I asked a question about my simulation, which is now running with the help of many reactions on my previous questions.<br>
<br>
However, the results which I get from this calculation are higher than what I expected, with approximately a factor 2. Last week, therefore, I send out this question to the GATE-mailing list. Unfortunately, I got no response yet. I would therefore like to ask you if you maybe have the time to look at my question, and see if you possible have any idea where I should look for to acquire the correct dose results.<br>
<br>
For my simulation, I used the fluor-18 file from the dosimetry example from the gate-repository as my source. The simulation ran with 1 Bq in the brain-tissue, and 6 000 000 primary events. With a mean energy of 250 keV for the beta-particle of F-18 I would have expected a total deposited energy of 1.500.000 MeV for the 6 million primary events. However, when I sum over all pixel values from the resulting Edep.mhd / Edep.raw file I get 3.237.062 MeV.<br>
<br>
Do you have any idea where this difference comes from?<br>
<br>
Thanks in advance for your help!<br>
<br>
Kind regards,<br>
Niels<br>
______________________________<wbr>_________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.<wbr>opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.<wbr>opengatecollaboration.org/<wbr>mailman/listinfo/gate-users</a></blockquote></div><br></div>