<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="-1">Hello,<br>
I had some problems using Radioactive Decay with some physics list
in the past.<br>
I'm quite sure that RadioactiveDecay didn't work with all physics
list (may be it's OK now with GATE v7).<br>
With emstandard it's OK.... but with QGSP_BIC_EMY I think it
didn't work for me.<br>
Regards.<br>
Nicolas<br>
</font>
<pre class="moz-signature" cols="72">***********************************************************************
Nicolas ARBOR
Maître de Conférences
Université de Strasbourg - UFR Physique et Ingénierie
Groupe de recherche RaMsEs - Institut Pluridisciplinaire Hubert Curien
23, rue du Loess BP28 67037 Strasbourg Cedex 2
mail : <a class="moz-txt-link-abbreviated" href="mailto:nicolas.arbor@iphc.cnrs.fr">nicolas.arbor@iphc.cnrs.fr</a> / tel : 03.88.10.64.27
***********************************************************************</pre>
<div class="moz-cite-prefix">Le 02/04/2015 10:52, David Sarrut a
écrit :<br>
</div>
<blockquote
cite="mid:CAKTJC4Dna056p0i_fbTeTFM92CjOBkGDUyiOXU=OJCVuiv2D9A@mail.gmail.com"
type="cite">
<div dir="ltr">Hello Jye,
<div><br>
</div>
<div>maybe try </div>
<div><span style="font-size:12.8000001907349px">/gate/physics/addPhysicsList
QGSP_BIC_EMY</span><br>
</div>
<div>
<div style=""><span style="font-size:12.8000001907349px">/gate/physics/addProcess
RadioactiveDecay</span><br>
</div>
<div style="font-size:12.8000001907349px"><br>
</div>
</div>
<div style="font-size:12.8000001907349px">or </div>
<div style="">
<div style="font-size:small"><span
style="font-size:12.8000001907349px">/gate/physics/addPhysicsList
emstandard_opt3</span></div>
<div style="font-size:12.8000001907349px">
<div style="font-size:small"><span
style="font-size:12.8000001907349px">/gate/physics/addProcess
RadioactiveDecay</span><br>
</div>
</div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px"><br>
</span></div>
<div style=""><span style="font-size:12.8000001907349px">instead
of all the lines in your # PHYSICS section</span></div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px"><br>
</span></div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px">(I think decay is not
included in the phys list, to be checked)</span></div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px"><br>
</span></div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px">David</span></div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px"><br>
</span></div>
<div style="font-size:12.8000001907349px"><span
style="font-size:12.8000001907349px"><br>
</span></div>
</div>
<div style="font-size:12.8000001907349px"><br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Apr 2, 2015 at 1:45 AM, Jye
Smith <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:jye.smith@gmail.com" target="_blank">jye.smith@gmail.com</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>Hello GATE,</div>
<div><br>
</div>
<div>Im currently using GATE for therapy dosimetry
applications and initially ran some test macros with
131I, 177Lu, 90Y to validate the results against
OLINDA. Everything looked great, even when running the
macro with mono energetic particles and comparing the
results against published literature.</div>
<div><br>
</div>
<div>Recently I decided to have a look at other diagnostic
isotopes 99mTc, 67Ga etc and some alpha emitters e.g.
223Ra. I found the marco I was using (originally based
on the mouse dosimetry example) did not give results
close OLINDA when I simply changed the ion source. Ive
got a feeling Im not including all of the required
parameters in the Physics and tried the PhysicsList
option with out luck.</div>
<div><br>
</div>
<div>Ive uploaded a working example to <a
moz-do-not-send="true"
href="http://dropcanvas.com/53wei" target="_blank">http://dropcanvas.com/53wei</a>
which is a simple ~100gm sphere I can compare to OLINDAs
sphere results. The .xls file is a summary of the
results. You can see that for 18F and 177Lu the
difference in dose to the sphere is within a few
percent. However for 67Ga, 99mTc and 223Ra the
difference is large. It appear all positron/beta
emitters I model the results or good. For gamma and
alpha emitters the results are very different.</div>
<div><br>
</div>
<div>I think the correct PhysicsList I need is
QGSP_BIC_EMY (e.g. /gate/physics/addPhysicsList
QGSP_BIC_EMY) and it should allow me to simulate gamma,
beta and alphas by just changing the ion source. Ive
been doing this from the command line like below;</div>
<div><br>
</div>
<div>Gate -a
[sourcefile,source.dat][numparticles,1000000][outfilename,1.0M_18F.hdr][AtomicNumber,9][AtomicMass,18][ExcitationEnergy,0]
main.mac<br>
</div>
<div><br>
</div>
<div>Any help in correcting the Physics list (or what
every else is wrong) is greatly appreciated.</div>
<div><br>
</div>
<div>Cheers</div>
<div>Jye</div>
<div><br>
</div>
<div>main.mac copied here for convenience and to help
anyone searching for solutions in the future.</div>
<div><br>
</div>
<div>
<div>#=====================================================</div>
<div># How to run macro</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>#Gate -a
[sourcefile,source.dat][numparticles,1000000][outfilename,1.0M_18F.hdr][AtomicNumber,9][AtomicMass,18][ExcitationEnergy,0]
main.mac</div>
<div>#Gate -a
[sourcefile,source.dat][numparticles,1000000][outfilename,1.0M_177Lu.hdr][AtomicNumber,71][AtomicMass,177][ExcitationEnergy,0]
main.mac</div>
<div>#Gate -a
[sourcefile,source.dat][numparticles,1000000][outfilename,1.0M_67Ga.hdr][AtomicNumber,31][AtomicMass,67][ExcitationEnergy,0]
main.mac</div>
<div>#Gate -a
[sourcefile,source.dat][numparticles,1000000][outfilename,1.0M_99mTc.hdr][AtomicNumber,43][AtomicMass,99][ExcitationEnergy,142.6833]
main.mac</div>
<div>#Gate -a
[sourcefile,source.dat][numparticles,1000000][outfilename,1.0M_223Ra.hdr][AtomicNumber,88][AtomicMass,223][ExcitationEnergy,0]
main.mac</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># GEOMETRY</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/geometry/setMaterialDatabase GateMaterials.db</div>
<div><br>
</div>
<div># World</div>
<div>/gate/world/geometry/setXLength 1 m</div>
<div>/gate/world/geometry/setYLength 1 m</div>
<div>/gate/world/geometry/setZLength 1 m</div>
<div>/gate/world/setMaterial Air</div>
<div><br>
</div>
<div># Global Box</div>
<div>/gate/world/daughters/name sphere</div>
<div>/gate/world/daughters/insert compressedMatrix</div>
<div><br>
</div>
<div>/gate/sphere/geometry/insertReader interfile</div>
<div>/gate/sphere/interfileReader/insertTranslator range</div>
<div>/gate/sphere/interfileReader/rangeTranslator/readTable
range_geo.dat</div>
<div><br>
</div>
<div>/gate/sphere/interfileReader/readFile sphere.iff</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># PHYSICS</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/physics/addProcess RadioactiveDecay</div>
<div><br>
</div>
<div>/gate/physics/addProcess PhotoElectric</div>
<div>/gate/physics/processes/PhotoElectric/setModel
StandardModel</div>
<div><br>
</div>
<div>/gate/physics/addProcess Compton</div>
<div>/gate/physics/processes/Compton/setModel
StandardModel</div>
<div><br>
</div>
<div>/gate/physics/addProcess RayleighScattering</div>
<div>/gate/physics/processes/RayleighScattering/setModel
PenelopeModel</div>
<div><br>
</div>
<div>/gate/physics/addProcess ElectronIonisation</div>
<div>/gate/physics/processes/ElectronIonisation/setModel
StandardModel e-</div>
<div>/gate/physics/processes/ElectronIonisation/setModel
StandardModel e+</div>
<div><br>
</div>
<div>/gate/physics/addProcess Bremsstrahlung e-</div>
<div>/gate/physics/addProcess Bremsstrahlung e+</div>
<div>/gate/physics/processes/Bremsstrahlung/setModel
StandardModel e-</div>
<div>/gate/physics/processes/Bremsstrahlung/setModel
StandardModel e+</div>
<div><br>
</div>
<div>/gate/physics/addProcess PositronAnnihilation e+</div>
<div><br>
</div>
<div>/gate/physics/addProcess MultipleScattering e+</div>
<div>/gate/physics/addProcess MultipleScattering e-</div>
<div><br>
</div>
<div>/gate/physics/addProcess GammaConversion</div>
<div>/gate/physics/processes/GammaConversion/setModel
StandardModel</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># DETECTORS</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/actor/addActor
DoseActor doseDistribution</div>
<div>/gate/actor/doseDistribution/attachTo <span
style="white-space:pre-wrap"> </span> <span
style="white-space:pre-wrap"> </span> sphere</div>
<div>/gate/actor/doseDistribution/stepHitType
random</div>
<div>/gate/actor/doseDistribution/setPosition
0 0 0 cm</div>
<div><br>
</div>
<div>/gate/actor/doseDistribution/setResolution
128 128 128</div>
<div><br>
</div>
<div>/gate/actor/doseDistribution/enableEdep
false</div>
<div>/gate/actor/doseDistribution/enableUncertaintyEdep
false</div>
<div>/gate/actor/doseDistribution/enableSquaredEdep
false</div>
<div>/gate/actor/doseDistribution/enableDose
true</div>
<div>/gate/actor/doseDistribution/enableUncertaintyDose
false</div>
<div>/gate/actor/doseDistribution/enableSquaredDose
false</div>
<div>/gate/actor/doseDistribution/enableNumberOfHits
false</div>
<div>/gate/actor/doseDistribution/save <span
style="white-space:pre-wrap"> </span>output/{outfilename}</div>
<div><br>
</div>
<div>/gate/actor/addActor <span
style="white-space:pre-wrap"> </span>SimulationStatisticActor
stat</div>
<div>/gate/actor/stat/save <span
style="white-space:pre-wrap"> </span>output/{outfilename}-output-stat.txt</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># INITIALISATION</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/run/initialize</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># SOURCE :</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/source/addSource isotope voxel</div>
<div><br>
</div>
<div>/gate/source/isotope/reader/insert interfile</div>
<div>/gate/source/isotope/interfileReader/translator/insert
range</div>
<div>/gate/source/isotope/interfileReader/rangeTranslator/readTable
{sourcefile}</div>
<div><br>
</div>
<div>/gate/source/isotope/interfileReader/readFile
sphere.iff</div>
<div><br>
</div>
<div>/gate/source/isotope/setPosition -64. -64.
-64. mm </div>
<div><br>
</div>
<div>/gate/source/isotope/gps/particle ion</div>
<div>/gate/source/isotope/gps/ion {AtomicNumber}
{AtomicMass} 0 {ExcitationEnergy}</div>
<div><br>
</div>
<div>/gate/source/isotope/gps/energytype Mono</div>
<div>/gate/source/isotope/gps/monoenergy 0. keV</div>
<div><br>
</div>
<div>/gate/source/isotope/gps/angtype iso </div>
<div><br>
</div>
<div>#=====================================================</div>
<div># VISUALISATION</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/vis/disable</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># RANDOM</div>
<div>#=====================================================</div>
<div>/gate/random/setEngineName MersenneTwister</div>
<div>/gate/random/setEngineSeed auto</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># START</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/application/noGlobalOutput</div>
<div>/gate/application/setTotalNumberOfPrimaries
{numparticles}</div>
<div>/gate/application/start</div>
<div><br>
</div>
<div>exit</div>
</div>
<div><br>
</div>
</div>
<br>
_______________________________________________<br>
Gate-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a><br>
<a moz-do-not-send="true"
href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users"
target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="gmail_signature">
<div dir="ltr">David Sarrut, Phd<br>
Directeur de recherche CNRS<br>
CREATIS, UMR CNRS 5220, Inserm U 1044
<div>Centre de lutte contre le cancer Léon Bérard<br>
28 rue Laënnec, 69373 Lyon cedex 08<br>
Tel : 04 78 78 51 51 / 06 74 72 05 42<br>
<a moz-do-not-send="true"
href="http://www.creatis.insa-lyon.fr/%7Edsarrut"
target="_blank">http://www.creatis.insa-lyon.fr/~dsarrut</a><br>
_________________________________</div>
<div> "2 + 2 = 5, for extremely large values of 2"<br>
_________________________________</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Gate-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a>
<a class="moz-txt-link-freetext" href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></pre>
</blockquote>
<br>
</body>
</html>