<div dir="ltr"><div class="gmail_quote"><div class="gmail_attr" dir="ltr"><br></div><br><div dir="ltr"><div dir="ltr"><div class="gmail_attr" dir="ltr">Hi, all Gate users, </div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><br></div><div>I simulated mouse whole body dosimetry with 68Ga-labeled radiotracer for 10^8 number of events. and I got too many NA values(or value=0) in the 3D output ( (i.e., Dose and Edep<span><span>).</span></span></div><div><br></div><div>I discussed various causes affecting these results with my co-worker. <br></div><div>I guessed there is the wrong code in "geometry step", <span><span>because the black dots (NA values or value=0) are dispersed on the air materials. </span></span></div><div><span><span>But i couldn't find the exact ways to solve these problems. </span></span></div><div><br></div><div> The examples of the black dots are below.</div><div><div><img width="561" height="483" alt="image.png" src="cid:ii_kem973zd0"><br></div><br></div><div><div><br></div></div><div>My questions are obvious: </div><div><br></div><div>1. Why did i get these NA values in every 3D output? </div><div>2. How can i solve this problem? </div><div><br></div><div><br></div><div><div dir="auto">#=====================================================</div><div dir="auto"># VEROBOSITY</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">/control/execute mac/verbose.mac</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># VISUALISATION</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">/control/execute mac/visu.mac</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># GEOMETRY</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">/gate/geometry/setMaterialDatabase <span style="white-space:pre-wrap"> </span>data/GateMaterials.db</div><div dir="auto"><br></div><div dir="auto"># World</div><div dir="auto">/gate/world/geometry/setXLength <span style="white-space:pre-wrap"> </span>10 cm</div><div dir="auto">/gate/world/geometry/setYLength <span style="white-space:pre-wrap"> </span>10 cm</div><div dir="auto">/gate/world/geometry/setZLength <span style="white-space:pre-wrap"> </span>10 cm</div><div dir="auto">/gate/world/setMaterial <span style="white-space:pre-wrap"> </span>Air</div><div dir="auto">/gate/world/vis/setVisible <span style="white-space:pre-wrap"> </span>1</div><div dir="auto">/gate/world/vis/forceWireframe</div><div dir="auto"><br></div><div dir="auto"># Voxel phantom</div><div dir="auto">/gate/world/daughters/name ct</div><div dir="auto">/gate/world/daughters/insert ImageRegularParametrisedVolume</div><div dir="auto">/gate/ct/geometry/setRangeToMaterialFile <span style="white-space:pre-wrap"> </span>data/range_geo.dat</div><div dir="auto">/gate/ct/geometry/setImage <span style="white-space:pre-wrap"> </span>data/ct-gate.mhd</div><div dir="auto">/gate/ct/setSkipEqualMaterials <span style="white-space:pre-wrap"> </span>1</div><div dir="auto">/gate/ct/vis/setVisible <span style="white-space:pre-wrap"> </span>1</div><div dir="auto"># /gate/ct/vis/setColor <span style="white-space:pre-wrap"> </span>blue</div><div dir="auto">/gate/ct/vis/forceWireframe</div><div dir="auto"><br></div><div dir="auto"># /gate/geometry/rebuild</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># PHYSICS</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">/gate/physics/addPhysicsList <span style="white-space:pre-wrap"> </span>emstandard_opt3</div><div dir="auto">/gate/physics/addProcess <span style="white-space:pre-wrap"> </span>RadioactiveDecay #to generate decay process of F-18</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># DETECTORS</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">/gate/actor/addActor DoseActor doseDistribution</div><div dir="auto">/gate/actor/doseDistribution/attachTo <span style="white-space:pre-wrap"> </span> <span style="white-space:pre-wrap"> </span>ct</div><div dir="auto">/gate/actor/doseDistribution/stepHitType random</div><div dir="auto">/gate/actor/doseDistribution/setSize 30 28.8 99 mm #Size in mm</div><div dir="auto">/gate/actor/doseDistribution/setPosition 0 0 0 cm # Origin</div><div dir="auto">/gate/actor/doseDistribution/setResolution 58 48 165 #Size (pixel)</div><div dir="auto">/gate/actor/doseDistribution/enableEdep true</div><div dir="auto">/gate/actor/doseDistribution/enableUncertaintyEdep true</div><div dir="auto">/gate/actor/doseDistribution/enableSquaredEdep false</div><div dir="auto">/gate/actor/doseDistribution/enableDose true #from true</div><div dir="auto">/gate/actor/doseDistribution/setDoseAlgorithm <span style="white-space:pre-wrap"> </span>VolumeWeighting</div><div dir="auto"># /gate/actor/doseDistribution/enableDoseToWater true #new set</div><div dir="auto">/gate/actor/doseDistribution/enableUncertaintyDose true #from true</div><div dir="auto">/gate/actor/doseDistribution/enableSquaredDose false</div><div dir="auto">/gate/actor/doseDistribution/enableNumberOfHits true</div><div dir="auto">/gate/actor/doseDistribution/save output/output.mhd</div><div dir="auto"><br></div><div dir="auto">/gate/actor/addActor <span style="white-space:pre-wrap"> </span>SimulationStatisticActor stat</div><div dir="auto">/gate/actor/stat/save <span style="white-space:pre-wrap"> </span>output/stat.txt</div><div dir="auto">/gate/actor/stat/saveEveryNSeconds <span style="white-space:pre-wrap"> </span>60</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># INITIALISATION</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">/gate/run/initialize</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># SOURCE :</div><div dir="auto">#=====================================================</div><div dir="auto"><br></div><div dir="auto">#Description</div><div dir="auto"><br></div><div dir="auto">#Ion</div><div dir="auto">/gate/source/addSource<span style="white-space:pre-wrap"> </span>ion<span style="white-space:pre-wrap"> </span>voxel</div><div dir="auto">/gate/source/ion/reader/insert <span style="white-space:pre-wrap"> </span>image</div><div dir="auto">/gate/source/ion/imageReader/translator/insert <span style="white-space:pre-wrap"> </span>linear</div><div dir="auto">/gate/source/ion/imageReader/linearTranslator/setScale<span style="white-space:pre-wrap"> </span>1. Bq</div><div dir="auto">/gate/source/ion/imageReader/readFile <span style="white-space:pre-wrap"> </span>data/pt-gate.mhd </div><div dir="auto">/gate/source/ion/imageReader/verbose <span style="white-space:pre-wrap"> </span>1</div><div dir="auto">/gate/source/ion/setPosition <span style="white-space:pre-wrap"> </span>-15.0 -14.4 -49.5 mm</div><div dir="auto"><br></div><div dir="auto">/gate/source/ion/gps/particle <span style="white-space:pre-wrap"> </span>ion</div><div dir="auto">/gate/source/ion/gps/ion <span style="white-space:pre-wrap"> </span>9 18 0 0 #F-18</div><div dir="auto">/gate/source/ion/gps/energytype <span style="white-space:pre-wrap"> </span>Mono</div><div dir="auto">/gate/source/ion/gps/monoenergy <span style="white-space:pre-wrap"> </span>0. keV</div><div dir="auto">/gate/source/ion/gps/angtype <span style="white-space:pre-wrap"> </span>iso</div><div dir="auto">/gate/source/ion/setForcedUnstableFlag <span style="white-space:pre-wrap"> </span>true </div><div dir="auto">/gate/source/ion/useDefaultHalfLife<span style="white-space:pre-wrap"> </span>6586 s</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">/gate/source/ion/dump <span style="white-space:pre-wrap"> </span>1</div><div dir="auto"><br></div><div dir="auto"># /gate/source/list</div><div dir="auto"><br></div><div dir="auto">#=====================================================</div><div dir="auto"># RANDOM</div><div dir="auto">#=====================================================</div><div dir="auto">/gate/random/setEngineName MersenneTwister</div><div dir="auto">/gate/random/setEngineSeed auto</div><span><span><br></span></span></div><div> <div dir="auto">#=====================================================</div><div dir="auto"># START</div><div>#=====================================================</div><div><div dir="auto">/gate/application/noGlobalOutput</div><div>/gate/application/setTotalNumberOfPrimaries 100000000</div><div>/gate/application/setTimeStart 0s</div><div>/gate/application/setTimeSlice 1s</div><div>/gate/application/setTimeStop 60s</div><div>/gate/application/start</div><div><br></div><div><br></div><div><br></div><div><br></div><div>Thank you for your help</div><div><br></div><div>Best regards.</div><div><br></div><div>Subin Kim<span></span></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div>