<div dir="ltr">Thanks Jenny, David and Nicolas for your replies.<div><br></div><div>I deleted everything in the physics section and replaced it with the below lines. I tried a few different libraries.</div><div><br></div><div><div style="font-size:13px"><b><span style="font-size:12.8000001907349px">/gate/physics/addPhysicsList QGSP_BIC_EMY</span><br></b></div><div style="font-size:13px"><span style="font-size:12.8000001907349px"><b>/gate/physics/addProcess RadioactiveDecay</b></span></div></div><div style="font-size:13px"><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">And got the below error </span>every time<span style="font-size:12.8000001907349px">.</span></div><div style="font-size:13px"><span style="font-size:12.8000001907349px"><br></span></div><div style><div style><b>### === Deexcitation model UAtomDeexcitation is activated for 3 regions:</b></div><div style><b> sphere</b></div><div style><b><br></b></div><div style><b> *** Break *** segmentation violation</b></div></div><div style="font-size:13px"><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">After hunting through the mailing list and a bit of google-fu it looks like a memory issue. Im currently running vGate with 4G of RAM. Can someone </span>running<span style="font-size:12.8000001907349px"> more RAM please confirm the macro will run fine. Im in the process of installing Gate on a server (16G) to try it myself but am running into proxy issues :(</span></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">Cheers</span></div><div><span style="font-size:12.8000001907349px">Jye</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 2, 2015 at 9:45 AM, Jye Smith <span dir="ltr"><<a 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 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>
</blockquote></div><br></div>