<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 <matthew.strugari@dal.ca><br />À: Zhengzhi Liu <zliu36@stanford.edu>, Tien-Phong PHAM <tienphong.pham@laposte.net><br />Cc: gate-users <gate-users@lists.opengatecollaboration.org><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="LPlnk903858" 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="LPlnk759400" 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="divRplyFwdMsg" dir="ltr"><span style="color: #000000; font-family: Calibri,sans-serif; font-size: 11pt;" color="#000000" face="Calibri, sans-serif"><b>From:</b> Gate-users <gate-users-bounces@lists.opengatecollaboration.org> on behalf of Zhengzhi Liu <zliu36@stanford.edu><br /><b>Sent:</b> March 2, 2020 1:22 PM<br /><b>To:</b> Tien-Phong PHAM <tienphong.pham@laposte.net><br /><b>Cc:</b> gate-users <gate-users@lists.opengatecollaboration.org><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 class="x_gmail_quote">
<div class="x_gmail_attr" 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 class="x_gmail_quote" 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="x_gmail-m_7204315926884882512EDA6169A-608B-48A1-9684-15705A3482A1" style="max-width: 100%; height: auto;" src="/service/home/~/?id=27694&part=2&auth=co" dfsrc="cid:1709c433b969093c11c1" /></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" target="_blank" rel="noreferrer">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></blockquote>
</div>
</div>
</div>
<div></div>
<div></div>