<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>Dear Gate users,<br></div><div><br data-mce-bogus="1"></div><div>I followed David and Dimitrios advices and reduced the SetCutInRegion values in my code.<br data-mce-bogus="1"></div><div>David, thank you, I deleted the  range cut for positrons and I defined different range cuts for the volumes.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>I still obtain no significant Gy values differences as output <strong>with</strong> and <strong>without</strong> nanoparticles...<br data-mce-bogus="1"></div><div><br></div><div>I'm searching, but if someone as any idea, I would be very grateful for the help or any advice!<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Best regards,<br data-mce-bogus="1"></div><div>Sarah Blind<br data-mce-bogus="1"></div><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>De: </b>"Sarah Blind" <sarah.blind@univ-lorraine.fr><br><b>À: </b>"David Leibold" <D.Leibold@tudelft.nl><br><b>Cc: </b>"gate-users" <gate-users@lists.opengatecollaboration.org><br><b>Envoyé: </b>Lundi 7 Février 2022 17:27:59<br><b>Objet: </b>Re: [Gate-users] with and without nanoparticules -> surprising output<br></div><div><br></div><div data-marker="__QUOTED_TEXT__"><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>Dear David, dear Dimitrios,</div><br><div>Thank you very much for your quick responses!</div><div>I'm testing your advices, I will keep you informed of the result.</div><br><div>Regards,</div><div>Sarah Blind</div><br><div><span id="zwchr">----- Le 7 Fév 22, à 16:51, David Leibold <D.Leibold@tudelft.nl> a écrit :<br></span></div><div><blockquote style="border-left:2px solid #1010FF;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;">
Dear Sarah,
<div class=""><br class="">
</div>
<div class="">while I do not have the knowledge to tell you whether your result makes sense or not, I had the same thought as Dimitris. A colleague once told me to choose a range cut that is about a fourth or fifth of the smallest dimension involved in your
 simulation. I am not sure whether this also holds true at the very small length scales that you are using in your simulation; others might shed more light on this topic. </div>
<div class=""><br class="">
</div>
<div class="">You probably know that decreasing the range cut increases the simulation time, too. You might want to decrease the the range cut in small steps (factors of 10) and check both the simulation time and whether the results still change noticeably
 below a certain range cut.</div>
<div class=""><br class="">
</div>
<div class="">It might also make sense to set different range cuts for your various volumes, so the smaller your structure, the smaller the range cut. By doing so, you help keeping the computation time as low as possible.</div>
<div class=""><br class="">
</div>
<div class="">Last but not least, I noticed that you set a range cut for positrons. If you use X-ray energies as specified in your Gate macro (170 keV) then you don’t have to care about positrons, since the energy is too small to create them in the first place.</div>
<div class=""><br class="">
</div>
<div class="">Kind regards,</div>
<div class="">David Leibold</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">
<div><br class="">
<blockquote class="">
<div class="">On 7 Feb2022, at 10:15, Sarah Blind <<a href="mailto:sarah.blind@univ-lorraine.fr" class="" target="_blank">sarah.blind@univ-lorraine.fr</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">
<div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;" class="">
<div class="">Dear Gate users,<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I want to model a spheroid which contains 30 000 cells in each of which there are 10 clusters of 10 nanoparticles.<br class="">
</div>
<div class="">The composition of my nanoparticles is the following<br class="">
</div>
<div class=""><br class="">
</div>
<blockquote class="">
<div class="">    +el: name=Gadolinium  ;   f=0.15<br class="">
    +el: name=Silicon     ;   f=0.10<br class="">
    +el: name=Carbon      ;   f=0.25<br class="">
    +el: name=Nitrogen    ;   f=0.08<br class="">
    +el: name=Oxygen      ;   f=0.22<br class="">
    +el: name=Hydrogen    ;   f=0.20<br class="">
</div>
</blockquote>
<div class=""><br class="">
</div>
<div class="">My spheroid, the cells and the clusters are all spheres full with water.</div>
<div class="">For the moment, what I'm trying to observe is very simple : <strong class="">
with</strong> and <strong class="">without</strong> the nanoparticles (and thus only water) I should observe a difference, shouldn't I ?<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I use Geant4-DNA; with 10 000 000 photons sent ; here is what my .hdr outpout file looks like</div>
<div class=""><br class="">
</div>
<div class="">
<table style="width: 16.2672%; height: 54px;" cellpadding="3px" border="1px">
<tbody class="">
<tr style="height: 18px;" class="">
<td style="width: 1.18203%; height: 18px;" class="">Nanoparticles<br class="">
</td>
<td style="width: 4.63031%; height: 18px;" class="">Average in Gy<br class="">
</td>
</tr>
<tr style="height: 18px;" class="">
<td style="width: 1.18203%; height: 18px;" class="">No<br class="">
</td>
<td style="width: 4.63031%; height: 18px;" class="">0.002<br class="">
</td>
</tr>
<tr style="height: 18px;" class="">
<td style="width: 1.18203%; height: 18px;" class="">Yes<br class="">
</td>
<td style="width: 4.63031%; height: 18px;" class="">
<div class="">0.001</div>
</td>
</tr>
</tbody>
</table>
<p class="">There is almost no difference between with and without the nanoparticles, this surprises me.<br class="">
</p>
<p class="">I'm not a physicist, do you find this normal ?<br class="">
</p>
<p class="">Here is my Gate code :<br class="">
</p>
<p class=""><br class="">
</p>
<blockquote class="">
<p class="">======= world.mac =======<br class="">
<br class="">
<br class="">
/gate/world/geometry/setXLength 0.5 cm<br class="">
/gate/world/geometry/setYLength 0.5 cm<br class="">
/gate/world/geometry/setZLength 0.5 cm<br class="">
/gate/world/setMaterial Air<br class="">
<br class="">
<br class="">
/gate/world/daughters/name                       <strong class="">spheroide</strong><br class="">
/gate/world/daughters/insert                       sphere<br class="">
/gate/spheroide/geometry/setRmin             0.0 mm<br class="">
/gate/spheroide/geometry/setRmax            400 um<br class="">
/gate/spheroide/placement/setTranslation   0 0 0 mm<br class="">
/gate/spheroide/setMaterial                         G4_WATER<br class="">
/gate/spheroide/vis/setVisible                       1<br class="">
/gate/spheroide/vis/setColor                         blue<br class="">
/gate/spheroide/vis/forceSolid                      false<br class="">
<br class="">
<br class="">
/gate/spheroide/daughters/name                       <strong class="">cell</strong><br class="">
/gate/spheroide/daughters/insert                       sphere<br class="">
/gate/cell/geometry/setRmin                              0.0 mm<br class="">
/gate/cell/geometry/setRmax                             10 um<br class="">
/gate/cell/placement/setTranslation                   100 0 0 um<br class="">
/gate/cell/setMaterial                                         G4_WATER<br class="">
/gate/cell/vis/setVisible                                      1<br class="">
/gate/cell/vis/setColor                                         green<br class="">
/gate/cell/vis/forceSolid                                      false<br class="">
<br class="">
/gate/cell/daughters/name                           <strong class="">cluster</strong><br class="">
/gate/cell/daughters/insert                           sphere<br class="">
/gate/cluster/geometry/setRmin                   0 mm<br class="">
/gate/cluster/geometry/setRmax                  0.05 um<br class="">
/gate/cluster/placement/setTranslation        0 0 0 mm<br class="">
/gate/cluster/setMaterial                               G4_WATER<br class="">
/gate/cluster/vis/setVisible                            1<br class="">
/gate/cluster/vis/setColor                              yellow<br class="">
/gate/cluster/vis/forceSolid                           false<br class="">
<br class="">
/gate/cluster/daughters/name                          <strong class="">NP</strong><br class="">
/gate/cluster/daughters/insert                          sphere<br class="">
/gate/NP/geometry/setRmin                            0 mm<br class="">
/gate/NP/geometry/setRmax                           0.004 um<br class="">
/gate/NP/placement/setTranslation                 0 0 0 mm<br class="">
/gate/NP/setMaterial                                       AGuIX<br class="">
/gate/NP/vis/setVisible                                    1<br class="">
/gate/NP/vis/setColor                                      red<br class="">
/gate/NP/vis/forceSolid                                    true<br class="">
<br class="">
### repeaters ### <br class="">
<br class="">
/gate/NP/repeaters/insert                               genericRepeater<br class="">
/gate/NP/genericRepeater/setPlacementsFilename          ../data/NP.placements<br class="">
/gate/NP/genericRepeater/useRelativeTranslation         1<br class="">
<br class="">
/gate/cluster/repeaters/insert                               genericRepeater<br class="">
/gate/cluster/genericRepeater/setPlacementsFilename          ../data/cluster.placements<br class="">
/gate/cluster/genericRepeater/useRelativeTranslation         1<br class="">
<br class="">
<br class="">
/gate/cell/repeaters/insert cubicArray<br class="">
/gate/cell/cubicArray/setRepeatNumberX 32<br class="">
/gate/cell/cubicArray/setRepeatNumberY 32<br class="">
/gate/cell/cubicArray/setRepeatNumberZ 32<br class="">
/gate/cell/cubicArray/setRepeatVector 20 20 20 um<br class="">
<br class="">
</p>
<p class="">======= physics.mac =======<br class="">
</p>
<p class=""><br class="">
# Geant4-DNA<br class="">
/gate/physics/addPhysicsListMixed emstandard_opt3_mixed_emdna<br class="">
/gate/physics/SetDNAInRegion spheroide<br class="">
/gate/physics/ConstructProcessMixed<br class="">
<br class="">
/gate/physics/Gamma/SetCutInRegion      world 0.5 mm<br class="">
/gate/physics/Electron/SetCutInRegion   world 0.5 mm<br class="">
/gate/physics/Positron/SetCutInRegion   world 0.5 mm<br class="">
<br class="">
/gate/physics/displayCuts<br class="">
<br class="">
# Choice of pseudo-random number generation algorithm<br class="">
/gate/random/setEngineName MersenneTwister <br class="">
<br class="">
# Auger electron activation <br class="">
/gate/physics/addAtomDeexcitation<br class="">
<br class="">
<br class="">
/gate/physics/addProcess                    PhotoElectric gamma<br class="">
/gate/physics/processes/PhotoElectric/setModel        LivermoreModel<br class="">
/gate/physics/processes/PhotoElectric/setAugerElectron     true<br class="">
/gate/physics/processes/PhotoElectric/setDeltaRayCut     10 eV<br class="">
/gate/physics/processes/PhotoElectric/setXRayCut         10 eV<br class="">
/gate/physics/processes/PhotoElectric/setStepFunction     e- 0.2 0.1 nm<br class="">
<br class="">
/gate/physics/addProcess                    Compton gamma<br class="">
/gate/physics/processes/Compton/setModel             LivermoreModel<br class="">
/gate/physics/processes/Compton/setStepFunction         e- 0.2 0.1 nm<br class="">
<br class="">
/gate/physics/addProcess                    RayleighScattering<br class="">
/gate/physics/processes/RayleighScattering/setModel         LivermoreModel<br class="">
<br class="">
# transformation of a photon into an electron-positron pair<br class="">
<br class="">
/gate/physics/addProcess                     GammaConversion<br class="">
/gate/physics/processes/GammaConversion/setModel         LivermoreModel<br class="">
<br class="">
/gate/physics/addProcess                     ElectronIonisation<br class="">
/gate/physics/processes/ElectronIonisation/setModel         LivermoreModel e-<br class="">
/gate/physics/processes/ElectronIonisation/setModel         PenelopeModel e+<br class="">
/gate/physics/processes/ElectronIonisation/setStepFunction     e- 0.2 0.1 nm<br class="">
<br class="">
/gate/physics/addProcess                     Bremsstrahlung e-<br class="">
/gate/physics/addProcess                     Bremsstrahlung e+<br class="">
/gate/physics/processes/Bremsstrahlung/setModel         LivermoreModel e-<br class="">
/gate/physics/processes/Bremsstrahlung/setModel         PenelopeModel e+<br class="">
<br class="">
/gate/physics/addProcess                     eMultipleScattering e-<br class="">
/gate/physics/addProcess                     eMultipleScattering e+<br class="">
/gate/physics/processes/eMultipleScattering/setModel     UrbanModel e-<br class="">
/gate/physics/processes/eMultipleScattering/setModel     UrbanModel e+</p>
<p class=""><br class="">
</p>
<p class="">======= detectors.mac =======<br class="">
</p>
<p class=""><br class="">
</p>
<p class="">/gate/actor/addActor    DoseActor  dose<br class="">
/gate/actor/dose/save    output/dose.hdr<br class="">
/gate/actor/dose/attachTo    spheroide<br class="">
/gate/actor/dose/stepHitType    random<br class="">
/gate/actor/dose/setPosition    0.0 0.0 0.0 mm<br class="">
/gate/actor/dose/setSize    800 800 800 um<br class="">
/gate/actor/dose/setResolution    400 400 400<br class="">
/gate/actor/dose/enableDose true<br class="">
/gate/actor/dose/saveEveryNSeconds 30</p>
<p class=""><br class="">
</p>
<p class="">======= source.mac =======<br class="">
</p>
<p class=""><br class="">
</p>
<p class="">/gate/source/addSource mybeam gps<br class="">
/gate/source/mybeam/gps/particle gamma<br class="">
/gate/source/mybeam/gps/pos/type Beam<br class="">
/gate/source/mybeam/gps/ene/mono 170 keV<br class="">
/gate/source/mybeam/gps/ene/type Gauss<br class="">
/gate/source/mybeam/gps/ene/sigma 1.0 keV<br class="">
/gate/source/mybeam/gps/direction 0 0 -1<br class="">
<br class="">
/gate/source/mybeam/gps/position 0 0 0.2 cm<br class="">
/gate/source/mybeam/gps/pos/type Volume<br class="">
/gate/source/mybeam/gps/pos/shape Cylinder<br class="">
/gate/source/mybeam/gps/pos/inner_radius 0 mm<br class="">
/gate/source/mybeam/gps/pos/radius 500 um<br class="">
/gate/source/mybeam/gps/pos/halfz 0.01 cm<br class="">
/gate/source/mybeam/gps/ang/type      beam2d</p>
<p class=""><br class="">
</p>
</blockquote>
<p class="">Thank you for your help,</p>
<p class="">Sarah Blind<br class="">
             <br class="">
<br class="">
<br class="">
<br class="">
</p>
</div>
</div>
</div>
_______________________________________________<br class="">
Gate-users mailing list<br class="">
<a href="mailto:Gate-users@lists.opengatecollaboration.org" class="" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br class="">
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.opengatecollaboration.org_mailman_listinfo_gate-2Dusers&d=DwIGaQ&c=XYzUhXBD2cD-CornpT4QE19xOJBbRy-TBPLK0X9U2o8&r=lfXR69GFfS7NT-Wp5HZqbtBDrbYoilDdmNEG7fKW7aM&m=EyDsGXEZgNS4-rEZqgSrucregTGsuMrgNHirJC_C8KA&s=VXjPJlYs_w1iutWjCIWnPcjLp2pRXeazsqZbeDKMbY0&e=
</div>
</blockquote>
</div>
<br class="">
</div><br></blockquote></div><br><div>-- <br></div><div>Sarah Blind<br>Ingénieure de recherche - CNRS<br>Centre de Recherche en Automatique de Nancy (CRAN)   <br>Département Biologie, Signaux et Systèmes en Cancérologie et Neurosciences (BioSiS)<br>Site de la faculté de médecine <br>9, Avenue de la Forêt de Haye <br>54505 Vandoeuvre les Nancy <br>Tél : 03.72.74.61.13</div></div><br>_______________________________________________<br>Gate-users mailing list<br>Gate-users@lists.opengatecollaboration.org<br>http://lists.opengatecollaboration.org/mailman/listinfo/gate-users<br></div><div><br></div><div data-marker="__SIG_POST__">-- <br></div><div>Sarah Blind<br>Ingénieure de recherche - CNRS<br>Centre de Recherche en Automatique de Nancy (CRAN)   <br>Département Biologie, Signaux et Systèmes en Cancérologie et Neurosciences (BioSiS)<br>Site de la faculté de médecine <br>9, Avenue de la Forêt de Haye <br>54505 Vandoeuvre les Nancy <br>Tél : 03.72.74.61.13</div></div></body></html>