<div dir="ltr"><div>Hello,</div><div><br></div><div>First of all, concerning the last email, you put in your code <br></div><div>/gate/source/Cu64Source/gps/monoenergy 0. MeV</div><div>which generally means that you declare that your source emits 0 energy. <br></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><br></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><br></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><br><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">tienphong.pham@laposte.net</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>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:rgb(0,0,0)"><br>Hi Zhengzhi and Tien-Phong,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><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:rgb(0,0,0)"><span>/gate/physics/addPhysicsList emstandard_opt4<br></span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><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:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:11pt" color="#000000" face="Calibri, sans-serif"><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>
<div></div>
# P H Y S I C A L P R O C E S S E S #<br>
<div></div>
# #<br>
<div></div>
#-------------------oooooOOOOO00000OOOOOooooo---------------------#<br>
<div></div>
<br>/gate/physics/addProcess PhotoElectric<br>
<div></div>
/gate/physics/processes/PhotoElectric/setModel StandardModel<br>
<div></div>
<br>/gate/physics/addProcess Compton<br>
<div></div>
/gate/physics/processes/Compton/setModel StandardModel<br>
<div></div>
<br>/gate/physics/addProcess RayleighScattering<br>
<div></div>
/gate/physics/processes/RayleighScattering/setModel PenelopeModel<br>
<div></div>
<br>/gate/physics/addProcess ElectronIonisation<br>
<div></div>
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-<br>
<div></div>
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+<br>
<div></div>
<br>/gate/physics/addProcess Bremsstrahlung<br>
<div></div>
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-<br>
<div></div>
/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>
/gate/physics/addProcess eMultipleScattering e-<br>
<div></div>
<br>/gate/physics/processList Enabled<br>
<div></div>
/gate/physics/processList Initialized<br>
<div></div>
<br>/gate/physics/addProcess IonIonisation<br>
<div></div>
/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 rgb(204,204,204);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="http:///service/home/~/?id=27694&part=2&auth=co"></p>
</div>
<br>_______________________________________________<br>
<div></div>
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>
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>