<div dir="ltr">Hi Ayeleth, <div><br></div><div>unsure why you cannot get information in the sinograms, but maybe look at this exercice that write and read a phase-space (it is for radiotherapy, but the idea should be similar) : </div><div><a href="https://gitlab.in2p3.fr/davidsarrut/gate-exercices/-/tree/master/linac">https://gitlab.in2p3.fr/davidsarrut/gate-exercices/-/tree/master/linac</a><br></div><div><br></div><div>Using actor frames or world frames is up to your needs. Maybe start with world frames it may be simpler. </div><div><br></div><div>I recommend to check your simulation with very few particle (1e4) to see if everything goes well </div><div><br></div><div>cheers, </div><div>David</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Apr 25, 2022 at 10:59 AM Devir-Wolfman, Ayeleth (GE Healthcare) <<a href="mailto:Ayeleth.Devir-Wolfman@ge.com">Ayeleth.Devir-Wolfman@ge.com</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">
<div lang="EN-US" style="overflow-wrap: break-word;">
<div class="gmail-m_6861192307080730571WordSection1">
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Good day to all!<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Sending my question again, hopefully with more responses this time.<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I'm try to implement phase space actor into my simulation, to be later used as the source of my simulation.<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I add the actor as following to the main macro:<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/addActor PhaseSpaceActor phaseSpace<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/save output/phaseSpace.root<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/attachTo Phantom<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"># Store not only entering particles, but also outgoing particles<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/storeOutgoingParticles true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/storeSecondaries false<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/setMaxFileSize 50 GB<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableProductionVolume false<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableProductionProcess false<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"># Store coordinates in the actor coordinates<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/useVolumeFrame true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableEkine true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableXPosition true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableYPosition true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableZPosition true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableXDirection true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableYDirection true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableZDirection true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableParticleName true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">/gate/actor/phaseSpace/enableWeight true<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Then after I have the phaseSpace.root file, I try to replace my source with this actor, by:<u></u><u></u></p>
<ol style="margin-top:0cm" start="1" type="1">
<li class="gmail-m_6861192307080730571MsoListParagraph" style="margin-right:0cm;text-align:left;direction:ltr;unicode-bidi:embed">
commenting out the above actor commands<u></u><u></u></li><li class="gmail-m_6861192307080730571MsoListParagraph" style="margin-right:0cm;text-align:left;direction:ltr;unicode-bidi:embed">
commenting out the previous source<u></u><u></u></li><li class="gmail-m_6861192307080730571MsoListParagraph" style="margin-right:0cm;text-align:left;direction:ltr;unicode-bidi:embed">
adding the following commands instead of the source, at the same place of the main macro:<u></u><u></u></li></ol>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"># Phase-Space actor as source<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">gate/source/addSource pS_Source phaseSpace<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">gate/source/pS_Source/addPhaseSpaceFile data/phaseSpace.root<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">gate/source/pS_Source/attachTo Phantom<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<ol style="margin-top:0cm" start="1" type="1">
<li class="gmail-m_6861192307080730571MsoListParagraph" style="margin-right:0cm;text-align:left;direction:ltr;unicode-bidi:embed">
My problem is that the output sinograms are empty, and that if I check the stat actor it shows that the simulation time has reached 1.79769e+299 seconds, where the simulation was run for 5 seconds.<u></u><u></u></li></ol>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<ol style="margin-top:0cm" start="2" type="1">
<li class="gmail-m_6861192307080730571MsoListParagraph" style="margin-right:0cm;text-align:left;direction:ltr;unicode-bidi:embed">
Also, I was wondering what is more recommended for this case, to save the particles coordinates in the actor frame or in the world frame (i.e. setting the command /gate/actor/phaseSpace/useVolumeFrame to be either true or false).<u></u><u></u></li></ol>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<ol style="margin-top:0cm" start="3" type="1">
<li class="gmail-m_6861192307080730571MsoListParagraph" style="margin-right:0cm;text-align:left;direction:ltr;unicode-bidi:embed">
And one general question: what is the difference between using the phase space actor and using the tracker/detector mode?<u></u><u></u></li></ol>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">I would greatly appreciate any help on the matter.<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed"><u></u> <u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Thank you,<u></u><u></u></p>
<p class="MsoNormal" style="text-align:left;direction:ltr;unicode-bidi:embed">Ayeleth.<span style="font-family:Arial,sans-serif"><u></u><u></u></span></p>
</div>
</div>
_______________________________________________<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></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>