<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear all,<div class=""><br class=""></div><div class="">I am computing Lu177 dosimetry using as input a CT and a PET. The conversion of the CT is done using the Schneider2000 table and the conversion of the PET is linear (1 Bq). The isotope used is Lu177 defined with a discrete spectrum and the number of primaries is 10^8. Overall the result looks ok, but there is unexpected high dose in some areas where there is Air. In the attachment you can find some examples and below the .mac file. Is that a known issue with GATE or is it due to an error in the .mac? Also, any other comment on how to improve the .mac file is welcome!</div><div class=""><br class=""></div><div class="">Best regards,</div><div class=""><br class=""></div><div class="">Francesca De Benetti</div><div class=""><br class=""></div><div class=""><img apple-inline="yes" id="5FA4986A-22E0-4DF1-93F7-04BE65372203" width="640" height="220" src="cid:278A0028-63F0-4F8B-84EF-0F0D1BE2378A" class=""><img apple-inline="yes" id="89846282-B7D0-4361-81C9-198A1F082D0A" width="640" height="345" src="cid:A2740FBC-DF0D-4ACF-BD14-2BD8E16668C2" class=""><img apple-inline="yes" id="AE7DA21A-2494-4722-BE01-30F20B56ED8C" width="640" height="345" src="cid:30790DF3-9735-4245-9819-03CE4FF30C85" class=""></div><div class=""><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># GEOMETRY</div><div class="">#=========================================================================================================</div><div class="">/gate/geometry/setMaterialDatabase <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/GateMaterials.db</div><div class=""><br class=""></div><div class=""># World</div><div class="">/gate/world/geometry/setXLength <span class="Apple-tab-span" style="white-space:pre"> </span>3 m</div><div class="">/gate/world/geometry/setYLength <span class="Apple-tab-span" style="white-space:pre"> </span>3 m</div><div class="">/gate/world/geometry/setZLength <span class="Apple-tab-span" style="white-space:pre"> </span>3 m</div><div class="">/gate/world/setMaterial <span class="Apple-tab-span" style="white-space:pre"> </span>{worldMaterial}</div><div class="">/gate/world/vis/setVisible <span class="Apple-tab-span" style="white-space:pre"> </span>0</div><div class=""><br class=""></div><div class=""># Automatic generation of the conversion values</div><div class="">/gate/HounsfieldMaterialGenerator/SetMaterialTable <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{conversion}MaterialsTable.txt</div><div class="">/gate/HounsfieldMaterialGenerator/SetDensityTable <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{conversion}DensitiesTable.txt</div><div class="">/gate/HounsfieldMaterialGenerator/SetDensityTolerance <span class="Apple-tab-span" style="white-space:pre"> </span><span class="Apple-tab-span" style="white-space:pre"> </span>0.1 g/cm3</div><div class="">/gate/HounsfieldMaterialGenerator/SetOutputMaterialDatabaseFilename <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{conversion}_CT-HUmaterials.db</div><div class="">/gate/HounsfieldMaterialGenerator/SetOutputHUMaterialFilename <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{conversion}_CT-HU2mat.txt </div><div class="">/gate/HounsfieldMaterialGenerator/Generate</div><div class=""><br class=""></div><div class=""># Daughter</div><div class="">/gate/world/daughters/name <span class="Apple-tab-span" style="white-space:pre"> </span>patient</div><div class="">/gate/world/daughters/insert <span class="Apple-tab-span" style="white-space:pre"> </span>ImageNestedParametrisedVolume</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span></div><div class="">/gate/geometry/setMaterialDatabase <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{conversion}_CT-HUmaterials.db</div><div class="">/gate/patient/geometry/setHUToMaterialFile <span class="Apple-tab-span" style="white-space:pre"> </span> <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{conversion}_CT-HU2mat.txt</div><div class=""><br class=""></div><div class="">/gate/patient/geometry/setImage <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/paddedCT.mhd</div><div class=""><br class=""></div><div class="">/gate/patient/placement/setTranslation <span class="Apple-tab-span" style="white-space:pre"> </span>0 0 0 mm</div><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># PHYSICS</div><div class=""><div class="">#=========================================================================================================</div></div><div class=""><br class=""></div><div class=""># Coulomb (do we need it?) scattering, multiple scattering, ionisation, Bremsstrahlung (change the version to add the pair production for electrons)</div><div class="">/gate/physics/addPhysicsList <span class="Apple-tab-span" style="white-space:pre"> </span>emlivermore</div><div class=""><br class=""></div><div class="">/gate/physics/addProcess <span class="Apple-tab-span" style="white-space:pre"> </span>Compton</div><div class="">/gate/physics/processes/Compton/setModel <span class="Apple-tab-span" style="white-space:pre"> </span>LivermoreModel</div><div class=""><br class=""></div><div class="">/gate/physics/addProcess <span class="Apple-tab-span" style="white-space:pre"> </span>RayleighScattering</div><div class="">/gate/physics/processes/RayleighScattering/setModel <span class="Apple-tab-span" style="white-space:pre"> </span>LivermoreModel</div><div class=""><br class=""></div><div class="">/gate/physics/addProcess<span class="Apple-tab-span" style="white-space:pre"> </span>RadioactiveDecay</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span></div><div class="">/gate/physics/addProcess <span class="Apple-tab-span" style="white-space:pre"> </span>HadronElastic</div><div class="">/gate/physics/processes/HadronElastic/setModel <span class="Apple-tab-span" style="white-space:pre"> </span>G4HadronElastic</div><div class="">/gate/physics/processes/HadronElastic/setDataSet <span class="Apple-tab-span" style="white-space:pre"> </span>G4HadronElasticDataSet</div><div class=""><br class=""></div><div class="">/gate/physics/processList <span class="Apple-tab-span" style="white-space:pre"> </span>Enabled </div><div class="">/gate/physics/processList <span class="Apple-tab-span" style="white-space:pre"> </span>Initialized</div><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># CUTS</div><div class=""><div class="">#=========================================================================================================</div></div><div class=""><br class=""></div><div class="">#/gate/physics/Electron/SetCutInRegion <span class="Apple-tab-span" style="white-space:pre"> </span>world 4 mm</div><div class="">#/gate/physics/Positron/SetCutInRegion <span class="Apple-tab-span" style="white-space:pre"> </span>world 4 mm</div><div class="">#/gate/physics/Gamma/SetCutInRegion <span class="Apple-tab-span" style="white-space:pre"> </span>world 4 mm</div><div class="">#/gate/physics/Proton/SetCutInRegion <span class="Apple-tab-span" style="white-space:pre"> </span>world 4 mm</div><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># ACTOR</div><div class="">#=========================================================================================================</div><div class=""><br class=""></div><div class="">/gate/actor/addActor <span class="Apple-tab-span" style="white-space:pre"> </span>DoseActor doseDistribution</div><div class="">/gate/actor/doseDistribution/save <span class="Apple-tab-span" style="white-space:pre"> </span>{outputFolder}/{outputDose}.mhd</div><div class="">/gate/actor/doseDistribution/attachTo <span class="Apple-tab-span" style="white-space:pre"> </span>patient</div><div class="">/gate/actor/doseDistribution/stepHitType <span class="Apple-tab-span" style="white-space:pre"> </span>random</div><div class="">/gate/actor/doseDistribution/setSize <span class="Apple-tab-span" style="white-space:pre"> </span>416.88 348.82 1104 mm</div><div class="">/gate/actor/doseDistribution/setPosition <span class="Apple-tab-span" style="white-space:pre"> </span>0 0 0 cm</div><div class="">/gate/actor/doseDistribution/setResolution <span class="Apple-tab-span" style="white-space:pre"> </span>588 492 368 </div><div class="">/gate/actor/doseDistribution/enableEdep <span class="Apple-tab-span" style="white-space:pre"> </span>false</div><div class="">/gate/actor/doseDistribution/enableUncertaintyEdep <span class="Apple-tab-span" style="white-space:pre"> </span>false</div><div class="">/gate/actor/doseDistribution/enableDose <span class="Apple-tab-span" style="white-space:pre"> </span>true</div><div class="">/gate/actor/doseDistribution/enableSquaredDose <span class="Apple-tab-span" style="white-space:pre"> </span>false</div><div class="">/gate/actor/doseDistribution/enableUncertaintyDose <span class="Apple-tab-span" style="white-space:pre"> </span>false</div><div class="">/gate/actor/doseDistribution/enableNumberOfHits <span class="Apple-tab-span" style="white-space:pre"> </span>false</div><div class="">/gate/actor/doseDistribution/saveEveryNSeconds <span class="Apple-tab-span" style="white-space:pre"> </span>60</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># INITIALIZATION</div><div class="">#=========================================================================================================</div><div class=""><br class=""></div><div class="">/gate/run/initialize</div><div class="">/gate/physics/displayCuts</div><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># SOURCE</div><div class="">#=========================================================================================================</div><div class=""><br class=""></div><div class="">/gate/source/addSource <span class="Apple-tab-span" style="white-space:pre"> </span>spect voxel</div><div class="">/gate/source/spect/reader/insert <span class="Apple-tab-span" style="white-space:pre"> </span>image</div><div class=""><br class=""></div><div class="">/gate/source/spect/imageReader/translator/insert <span class="Apple-tab-span" style="white-space:pre"> </span>linear</div><div class="">/gate/source/spect/imageReader/linearTranslator/setScale <span class="Apple-tab-span" style="white-space:pre"> </span>1 Bq</div><div class="">/gate/source/spect/imageReader/readFile <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/paddedPET.mhd</div><div class="">/gate/source/spect/imageReader/verbose <span class="Apple-tab-span" style="white-space:pre"> </span>1</div><div class=""><br class=""></div><div class="">/gate/source/spect/setPosition <span class="Apple-tab-span" style="white-space:pre"> </span>-208.44 -174.41 -552 mm</div><div class="">/gate/source/spect/dump <span class="Apple-tab-span" style="white-space: pre;"> </span>1</div><div class=""><br class=""></div><div class="">/gate/source/spect/gps/particle <span class="Apple-tab-span" style="white-space:pre"> </span>e-</div><div class="">/gate/source/spect/gps/energytype <span class="Apple-tab-span" style="white-space:pre"> </span>UserSpectrum</div><div class="">/gate/source/spect/gps/setSpectrumFile <span class="Apple-tab-span" style="white-space:pre"> </span>{dataDir}/{spectrumFile}.txt</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span></div><div class="">/gate/source/spect/setForcedUnstableFlag <span class="Apple-tab-span" style="white-space:pre"> </span>true</div><div class="">/gate/source/spect/setForcedHalfLife <span class="Apple-tab-span" style="white-space:pre"> </span>581817 s</div><div class=""><br class=""></div><div class="">/gate/source/spect/gps/angtype <span class="Apple-tab-span" style="white-space:pre"> </span>iso</div><div class="">/gate/source/spect/gps/centre <span class="Apple-tab-span" style="white-space:pre"> </span>0 0 0 mm</div><div class="">/gate/source/list</div><div class=""><br class=""></div><div class="">#=========================================================================================================</div><div class=""># START</div><div class="">#=========================================================================================================</div><div class=""><br class=""></div><div class="">/gate/random/setEngineName <span class="Apple-tab-span" style="white-space:pre"> </span>MersenneTwister</div><div class="">/gate/random/setEngineSeed <span class="Apple-tab-span" style="white-space:pre"> </span>default</div><div class="">/gate/application/setTotalNumberOfPrimaries <span class="Apple-tab-span" style="white-space:pre"> </span>{nPrimaries}</div><div class="">/gate/application/start</div><div class=""><br class=""></div></div></body></html>