<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:-webkit-standard;
        panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
        {font-family:Menlo;
        panose-1:2 11 6 9 3 8 4 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:#0563C1;
        text-decoration:underline;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-CA" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi Dmitry,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">You can read about the calculation of dose and its uncertainty
<a href="https://opengate.readthedocs.io/en/latest/voxelized_source_and_phantom.html?highlight=dose%20actor#dose-collection">
here</a>. The number of events required to reach 1% uncertainty is a bit of a loaded question. It will depend on your dose actor resolution and the location of the volume of interest so it may take some trial and error to find a suitable number. I only have
 dose actor experience working with beta-emitters on the order of a few MeV for which I required 4x10^8 primaries to achieve reasonable uncertainties.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Larger volumes can have lower uncertainty due to containment of more energy-depositing events. Dose voxels near your proton beam will encounter more interactions and lower uncertainty than voxels far from the beam where
 few interactions occur. Also keep in mind that your voxel size will have an effect on the dose since absorbed dose is calculated as energy deposited per unit mass.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Hope this helps,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Matthew</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:-webkit-standard;color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:-webkit-standard;color:black">--</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:-webkit-standard;color:black">Matthew Strugari<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:-webkit-standard;color:black">PhD Candidate</span><span style="color:black"> at
</span><span lang="EN-US" style="font-size:10.5pt;font-family:-webkit-standard;color:black">Biomedical Translational Imaging Centre - BIOTIC,<br>
5890 University Ave, Halifax, NS, B3K 6R8</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Gate-users <gate-users-bounces@lists.opengatecollaboration.org> on behalf of Dmitry Hits <dmitry.hits@phys.ethz.ch><br>
<b>Date: </b>Tuesday, May 11, 2021 at 12:31<br>
<b>To: </b>Gate-Users <gate-users@lists.opengatecollaboration.org><br>
<b>Subject: </b>[Gate-users] Dose results interpretation<o:p></o:p></span></p>
</div>
<div align="center">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;border-radius:5px;overflow:hidden">
<tbody>
<tr>
<td style="padding:7.5pt 0cm 7.5pt 0cm">
<div align="center">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;background:#707372">
<tbody>
<tr>
<td valign="top" style="border:none;border-top:solid #FBE122 6.0pt;padding:3.0pt 6.0pt 3.0pt 6.0pt">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0">
<tbody>
<tr>
<td style="padding:0cm 0cm 0cm 0cm">
<p class="MsoNormal" style="line-height:12.0pt;mso-line-height-rule:exactly"><b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:white">CAUTION:</span></b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:white"> The Sender
 of this email is not from within Dalhousie.<o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<div>
<p class="MsoNormal"><span style="color:white">Hi Gaters, </span><o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I am trying to understand results output by the dose actor, which I enabled with the following lines:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">/gate/actor/addActor DoseActor             dose5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/save                   <span style="font-size:8.5pt;font-family:Menlo;color:black">results/dose_2021May11_171652_5.root</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/attachTo               box5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/stepHitType            random<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/setResolution          1 10 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/enableDose             true<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/enableUncertaintyDose  true<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">        /gate/actor/dose5/enableNumberOfHits     true<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<div>
<p class="MsoNormal">The source is a pencil beam of protons of 100 MeV. I simulated 10000 events.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Here is the histogram I got:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><img border="0" width="761" height="435" style="width:7.927in;height:4.5312in" id="FB1EE12D-028F-4D15-BF0D-945FF1B50B9B" src="cid:image001.png@01D74675.94466BC0"><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">My question is why is the error so large? How is it calculated?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Cheers,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Dmitry<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</body>
</html>