<div dir="ltr">Hi Nicolas, <div><br></div><div>strange indeed. The two functions that compute the volume of the regions are:</div><div>- <a href="https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GateRegionDoseStat.cc#L67">https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GateRegionDoseStat.cc#L67</a></div><div>and </div><div>- <a href="https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GateRegionDoseStat.cc#L92">https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GateRegionDoseStat.cc#L92</a></div><div><br></div><div>Both called from here : <a href="https://github.com/OpenGATE/Gate/blob/86083caae3a4f41903b3d1e3c3d5571e807e041e/source/digits_hits/src/GateDoseActor.cc#L365">https://github.com/OpenGATE/Gate/blob/86083caae3a4f41903b3d1e3c3d5571e807e041e/source/digits_hits/src/GateDoseActor.cc#L365</a></div><div><br></div><div>I am unsure what is done there about the mMapLabelToSeveralRegions, but the behavior you observe may come from there. </div><div>Can you please provide us with a minimal working example that shows the bug ?</div><div><br></div><div>Cheer, </div><div>David</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, May 2, 2022 at 4:30 PM Nicolas Chouin <<a href="mailto:nicolas.chouin@oniris-nantes.fr">nicolas.chouin@oniris-nantes.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all,<br>
<br>
I have encountered one issue with the DoseByRegion option or maybe I did something wrong...<br>
My issue is related to the volumes displayed in the DoseByRegions output file. This volume did not correspond to the number of voxels (multiplied by the voxel size) I counted manually on my images. I can provide details of my macro if needed.<br>
Playing around with the DoseByRegions option, I observed that this issue only occurs for contiguous regions composed by the same material (but with a different ID; i.e. voxels have different values in each region). Even more striking, if I launched my simulations twice, only changing the material of the regions in the second simulation), I obtained different volumes in the DoseByRegions output files… whereas, obviously, the volume does not depend on the voxel composition.<br>
<br>
Has anyone encountered similar issues?<br>
<br>
Best regards,<br>
<br>
<br>
Nicolas Chouin, PhD<br>
Associate professor<br>
<br>
<br>
Research team: « Nuclear Oncology », CRCI2NA INSERM U1307, Nantes<br>
<br>
ONIRIS<br>
Site de la Chantrerie - CS 40706<br>
44307 Nantes, cedex 03, France<br>
Tel: +33-2-40-68-78-65<br>
<br>
_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><span style="color:rgb(0,0,0)">David Sarrut, Phd</span><br></div><div><font color="#000000">Directeur de recherche CNRS</font></div><div><font color="#000000">CREATIS, UMR CNRS 5220, Inserm U1294</font></div><div><font color="#000000">Centre de lutte contre le cancer Léon Bérard</font></div><div><font color="#000000">28 rue Laënnec, 69373 Lyon cedex 08</font></div><div><font color="#000000">Tel : 04 78 78 51 51 / 06 74 72 05 42</font></div><div><font color="#000000"><a href="http://dsarrut.github.io" target="_blank">http://dsarrut.github.io</a></font></div><div><font color="#000000">_________________________________</font></div><div><font color="#000000"> "2 + 2 = 5,  for extremely large values of 2"</font></div><div><font color="#000000">_________________________________</font></div></div></div>