<div>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN">Hi Sofia,<br /></span></span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN">I’ve done exactly as you said but it doesn’t seem to work. I still have the difference between the spectral ion source and spectral definition. </span></span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN">This is my spectrum that I obtained ( please see the attached file)<br /></span></span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">Furthermore, I did test with the source ion F18 and then compared with spectral definition. Unfortunately, I’ve got the same problem as for the Cu64 source ion.</span></span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">There are two sources that I tried: </span></span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">----------------------------------------------------------------------------------</span></span><span style="font-family: 'Times New Roman',serif;" lang="EN-US"></span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/addSource sourceF18</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/setActivity 1. MBq</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;">/gate/source/sourceF18/gps/particle e+</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/setForcedUnstableFlag true</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/energytype Fluor18</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/setForcedHalfLife 6586 s</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/type Volume</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/shape Sphere</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/radius 1. mm</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/centre 0. 0. 0. cm</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/confine NULL </span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/sourceF18/gps/angtype iso</span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">-------------------------------------------------------------------------------------</span></span><span style="font-family: 'Times New Roman',serif;" lang="EN-US"></span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/addSource fluor18</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/setActivity 1. MBq</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/particle ion</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/ion 9 18 0 0 </span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/energytype Mono</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/monoenergy 0. MeV</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/angtype iso</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/number 1</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/centre 0. 0. 0. cm</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/type Volume</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/shape Sphere</span></p>
<p style="margin-bottom: .0001pt; text-align: justify; line-height: normal;"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">/gate/source/fluor18/gps/radius 1. mm</span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">As far as I know, the difference between energy spectra of electron and positron can be explained by the fact that electron is attracted by the nucleus once it is created while the positron is repelled. So, I think I’m missing a model concerning coulomb interactions for the radioactive decay. I’m trying to look deeper into it, but I’ve yet to succeed. Maybe someone else has any other ideas to fix this problem?</span></span></p>
<p style="text-align: justify;"><span class="hps"><span style="font-family: 'Times New Roman',serif;" lang="EN-US">Best regards,</span></span></p>
</div>
<div></div>
<div>----- Mail d'origine -----<br />De: Sofia Kolovi <sofiaphys@gmail.com><br />À: tienphong pham <tienphong.pham@laposte.net><br />Cc: Matthew Strugari <matthew.strugari@dal.ca>, gate-users <gate-users@lists.opengatecollaboration.org><br />Envoyé: Mon, 02 Mar 2020 21:47:56 +0100 (CET)<br />Objet: Re: [Gate-users] Re : Re: 64 Cu positron spectra problem</div>
<div></div>
<div dir="ltr">
<div>Hello,</div>
<div></div>
<div>First of all, concerning the last email, you put in your code</div>
<div>/gate/source/Cu64Source/gps/monoenergy 0. MeV</div>
<div>which generally means that you declare that your source emits 0 energy.</div>
<div>The colleague above suggested to insert</div>
<div><span>/gate/source/pointSource/gps/energy 0 MeV</span></div>
<div><span>in order to declare that your source has 0 kinetic energy.<br /></span></div>
<div><span> </span></div>
<div><span>Apart from that, you may try to see if there is a bug with the radioactive decay for this ion.<br /></span></div>
<div><span>You can simply use the UserSpectrum (described in Source section of the documentation) where you describe on your own the energy spectrum emitted by your source and then compare with what you got using Radioactive Decay.<br /></span></div>
<div><span> </span></div>
<div><span>I hope you find this helpful,</span></div>
<div><span>Best regards,</span></div>
<div><span>Sofia<br /></span></div>
</div>
<div></div>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Mon, 2 Mar 2020 at 20:40, <<a href="mailto:tienphong.pham@laposte.net" target="_blank">tienphong.pham@laposte.net</a>> wrote:</div>
<blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid #cccccc; padding-left: 1ex;">
<div>Hi Matthew and Zhengzhi Liu,</div>
<div></div>
<div>Thank you very much for your quick response. I tried doing as you guys suggested but it still doesn't work. I don't know what's wrong yet.</div>
<div></div>
<div>Here is my code: </div>
<div></div>
<div>
<div>#=====================================================<br /># GEOMETRY<br />#=====================================================</div>
<div>/gate/geometry/setMaterialDatabase data/GateMaterials.db</div>
<div># World<br />/gate/world/geometry/setXLength 1. m<br />/gate/world/geometry/setYLength 1. m<br />/gate/world/geometry/setZLength 1. m<br />/gate/world/setMaterial Water</div>
<div># Water Box<br />/gate/world/daughters/name phantom<br />/gate/world/daughters/insert box<br />/gate/phantom/geometry/setXLength 20. cm<br />/gate/phantom/geometry/setYLength 20. cm<br />/gate/phantom/geometry/setZLength 20. cm<br />/gate/phantom/placement/setTranslation 0. 0. 0. cm<br />/gate/phantom/setMaterial Water<br />/gate/phantom/vis/setVisible 1<br />/gate/phantom/vis/setColor blue</div>
<div>#=====================================================<br /># PHYSIC <br />#=====================================================</div>
<div>/gate/physics/addPhysicsList emstandard_opt4 <br />/gate/physics/addProcess RadioactiveDecay</div>
<div>#=====================================================<br /># ACTOR<br />#=====================================================</div>
<div>#===========Positron spectra=========================================================================</div>
<div>/gate/actor/addActor EnergySpectrumActor Pos_Spectrum<br />/gate/actor/Pos_Spectrum/save ./output/Pos_Spectrum.root<br />/gate/actor/Pos_Spectrum/attachTo phantom<br />/gate/actor/Pos_Spectrum/addFilter particleFilter<br />/gate/actor/Pos_Spectrum/particleFilter/addParticle e+<br />/gate/actor/Pos_Spectrum/energySpectrum/setEmin 10 keV<br />/gate/actor/Pos_Spectrum/energySpectrum/setEmax 1 MeV<br />/gate/actor/Pos_Spectrum/energySpectrum/setNumberOfBins 1000<br />/gate/actor/Pos_Spectrum/enableNbPartSpectrum true</div>
<div>#===========Positron spectra==========================================================================</div>
<div>#===========Beta Spectra==============================================================================</div>
<div><br />/gate/actor/addActor EnergySpectrumActor Beta_Spectrum<br />/gate/actor/Beta_Spectrum/save ./output/Beta_Spectrum.root<br />/gate/actor/Beta_Spectrum/attachTo phantom</div>
<div>/gate/actor/Beta_Spectrum/addFilter particleFilter<br />/gate/actor/Beta_Spectrum/particleFilter/addParticle e-</div>
<div>/gate/actor/Beta_Spectrum/energySpectrum/setEmin 10 keV<br />/gate/actor/Beta_Spectrum/energySpectrum/setEmax 1 MeV<br />/gate/actor/Beta_Spectrum/energySpectrum/setNumberOfBins 1000<br />/gate/actor/Beta_Spectrum/enableNbPartSpectrum true</div>
<div>#===========Beta spectra=================================================================================</div>
<div><br />/gate/run/initialize</div>
<div>#=====================================================<br /># Source<br />#=====================================================</div>
<div>/gate/source/addSource Cu64Source gps<br />/gate/source/Cu64Source/gps/particle ion<br />/gate/source/Cu64Source/gps/ion 29 64 0 0<br />/gate/source/Cu64Source/useDefaultHalfLife<br />/gate/source/Cu64Source/setActivity 1. MBq<br />/gate/source/Cu64Source/gps/energytype Mono<br />/gate/source/Cu64Source/gps/monoenergy 0. MeV<br />/gate/source/Cu64Source/gps/angtype iso<br />/gate/source/Cu64Source/gps/centre 0. 0. 0. cm<br />/gate/source/Cu64Source/gps/type Volume<br />/gate/source/Cu64Source/gps/shape Sphere<br />/gate/source/Cu64Source/gps/radius 1. mm</div>
<div>#=====================================================<br /># MAIN<br />#=====================================================</div>
<div>/gate/random/setEngineName MersenneTwister<br />/gate/random/setEngineSeed 10</div>
<div># PARTICLE<br />/gate/application/setTotalNumberOfPrimaries 1000000<br />/gate/application/start</div>
<div></div>
<div></div>
<div>Tell me, please, what i'm missing here?</div>
<div></div>
<div></div>
<div>Thank you again for your help.</div>
<div></div>
<div>Tien-Phong PHAM</div>
</div>
<div></div>
<div></div>
<div>----- Mail d'origine -----<br />De: Matthew Strugari <<a href="mailto:matthew.strugari@dal.ca" target="_blank">matthew.strugari@dal.ca</a>><br />À: Zhengzhi Liu <<a href="mailto:zliu36@stanford.edu" target="_blank">zliu36@stanford.edu</a>>, Tien-Phong PHAM <<a href="mailto:tienphong.pham@laposte.net" target="_blank">tienphong.pham@laposte.net</a>><br />Cc: gate-users <<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>><br />Envoyé: Mon, 02 Mar 2020 19:17:43 +0100 (CET)<br />Objet: Re: [Gate-users] 64 Cu positron spectra problem</div>
<div></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: #000000;"><br />Hi Zhengzhi and Tien-Phong,</div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: #000000;"></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: #000000;"><br />Please refer to <a id="gmail-m_5133377077665207965LPlnk903858" href="https://opengate.readthedocs.io/en/latest/setting_up_the_physics.html#" target="_blank"><br />https://opengate.readthedocs.io/en/latest/setting_up_the_physics.html#</a> for information on physics list. The implementation suggested by Zhengzhi is being deprecated and should not be used. Rather, use the physics list builder mechanism:</div>
<blockquote style="margin-top: 0px; margin-bottom: 0px;">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: #000000;"><span>/gate/physics/addPhysicsList emstandard_opt4<br /></span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: #000000;"><span>/gate/physics/addProcess RadioactiveDecay</span></div>
</blockquote>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt;"><span>You may replace emstandard_opt4 with the physics list builder most suitable for your needs (see<br /><a id="gmail-m_5133377077665207965LPlnk759400" href="http://geant4-userdoc.web.cern.ch/geant4-userdoc/UsersGuides/PhysicsListGuide/html/index.html" target="_blank"><br />http://geant4-userdoc.web.cern.ch/geant4-userdoc/UsersGuides/PhysicsListGuide/html/index.html</a> for more info on physics lists.) Also, t<span>he following line may be required with your source definition to ensure that your Cu-64 source has no kinetic energy:</span><br /></span></div>
<blockquote style="margin-top: 0px; margin-bottom: 0px;">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt;"><span>/gate/source/pointSource/gps/energy 0 MeV</span></div>
</blockquote>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt;"></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt;"><span>Regards,</span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt;"><span>Matthew</span></div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt;"><span> </span></div>
<div><hr style="display: inline-block; width: 98%;" />
<div id="gmail-m_5133377077665207965divRplyFwdMsg" dir="ltr"><span style="color: #000000; font-family: Calibri,sans-serif; font-size: 11pt;"><b>From:</b> Gate-users <<a href="mailto:gate-users-bounces@lists.opengatecollaboration.org" target="_blank">gate-users-bounces@lists.opengatecollaboration.org</a>> on behalf of Zhengzhi Liu <<a href="mailto:zliu36@stanford.edu" target="_blank">zliu36@stanford.edu</a>><br /><b>Sent:</b> March 2, 2020 1:22 PM<br /><b>To:</b> Tien-Phong PHAM <<a href="mailto:tienphong.pham@laposte.net" target="_blank">tienphong.pham@laposte.net</a>><br /><b>Cc:</b> gate-users <<a href="mailto:gate-users@lists.opengatecollaboration.org" target="_blank">gate-users@lists.opengatecollaboration.org</a>><br /><b>Subject:</b> Re: [Gate-users] 64 Cu positron spectra problem</span>
<div> </div>
</div>
<div>
<div dir="ltr">Hi Tien-Phong, <br />
<div></div>
<div>Below is my physics.mac Could you compare this with yours? You might have missed something in your physics.mac. </div>
<div></div>
<div>Good luck,</div>
<div>Zhengzhi </div>
<div></div>
<div></div>
<div></div>
<div>#-------------------oooooOOOOO00000OOOOOooooo---------------------#<br />
<div></div>
<br /># #<br />
<div></div>
<br /># P H Y S I C A L P R O C E S S E S #<br />
<div></div>
<br /># #<br />
<div></div>
<br />#-------------------oooooOOOOO00000OOOOOooooo---------------------#<br />
<div></div>
<br />/gate/physics/addProcess PhotoElectric<br />
<div></div>
<br />/gate/physics/processes/PhotoElectric/setModel StandardModel<br />
<div></div>
<br />/gate/physics/addProcess Compton<br />
<div></div>
<br />/gate/physics/processes/Compton/setModel StandardModel<br />
<div></div>
<br />/gate/physics/addProcess RayleighScattering<br />
<div></div>
<br />/gate/physics/processes/RayleighScattering/setModel PenelopeModel<br />
<div></div>
<br />/gate/physics/addProcess ElectronIonisation<br />
<div></div>
<br />/gate/physics/processes/ElectronIonisation/setModel StandardModel e-<br />
<div></div>
<br />/gate/physics/processes/ElectronIonisation/setModel StandardModel e+<br />
<div></div>
<br />/gate/physics/addProcess Bremsstrahlung<br />
<div></div>
<br />/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-<br />
<div></div>
<br />/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+<br />
<div></div>
<br />/gate/physics/addProcess PositronAnnihilation<br />
<div></div>
<br />/gate/physics/addProcess eMultipleScattering e+<br />
<div></div>
<br />/gate/physics/addProcess eMultipleScattering e-<br />
<div></div>
<br />/gate/physics/processList Enabled<br />
<div></div>
<br />/gate/physics/processList Initialized<br />
<div></div>
<br />/gate/physics/addProcess IonIonisation<br />
<div></div>
<br />/gate/physics/addProcess RadioactiveDecay</div>
</div>
<br />
<div>
<div dir="ltr">On Mon, Mar 2, 2020 at 6:15 AM Tien-Phong PHAM <<a href="mailto:tienphong.pham@laposte.net" target="_blank">tienphong.pham@laposte.net</a>> wrote:</div>
<blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid #cccccc; padding-left: 1ex;">
<div>
<p>Dear Gate-users,</p>
<p>I’m trying to simulate the energy spectrum of a Cu-64 source using the ion definition in GATE v8.2, but i have a problem with the positron emission spectra. In fact, the energy distribution with my ion source is very different with spectral definition in<br /> the low energy region. I mean the energy spectrum of positron is similar to the energy spectrum of electron.</p>
<p>In this case, i just used the « RadioactiveDecay » as physics list</p>
<p>I would like to ask if someone has experience with ion type source can help me.</p>
<p>Many thanks for your help !<img id="gmail-m_5133377077665207965x_gmail-m_7204315926884882512EDA6169A-608B-48A1-9684-15705A3482A1" style="max-width: 100%; height: auto;" src="#" /></p>
</div>
<br />_______________________________________________<br />
<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>
</div>
</div>
<div></div>
<div></div>
<br />_______________________________________________<br />
<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>
<div></div>
<div></div>