<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Hi Sarah,</p>
<p>I believe this could be a problem with the implemented datatype.
The number of particles is stored internally in GATE as an
integer.<br>
</p>
<p>The highest number an integer can store is~2.1e9</p>
<p>Consequently simulating 1e9 particles works and 1e10 particles
will not work any more.</p>
<p>The easiest solution is to split your simulation in multiple
parts, each having less particles than 2.1e9.</p>
<p>All the best,</p>
<p>Hermann<br>
</p>
<div class="moz-cite-prefix">On 23.02.22 16:37, Sarah Blind wrote:<br>
</div>
<blockquote type="cite"
cite="mid:1614461897.1525339.1645630677551.JavaMail.zimbra@univ-lorraine.fr">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div style="font-family: arial, helvetica, sans-serif; font-size:
12pt; color: #000000">
<div>Dear Gate users, <br data-mce-bogus="1">
</div>
<div><br data-mce-bogus="1">
</div>
<div>I have the following model : 1 cell inside of which there
are 12 clusters of 8000 nanoparticles each.<br
data-mce-bogus="1">
</div>
<div>My source is composed of a gamma beam of 5 keV.<br
data-mce-bogus="1">
</div>
<div>My simulation works fine up to 10^9 photons send. But with
10^10 primaries nothing happens anymore (there is no error
message).<br data-mce-bogus="1">
</div>
<div>Does someone have an idea what the problem might be ?<br
data-mce-bogus="1">
</div>
<div><br data-mce-bogus="1">
</div>
<div>Thank you very much for your help,<br data-mce-bogus="1">
</div>
<div>Regards,<br data-mce-bogus="1">
</div>
<div>Sarah Blind<br data-mce-bogus="1">
</div>
<div><br data-mce-bogus="1">
</div>
<div>world.mac<br data-mce-bogus="1">
</div>
<div><br data-mce-bogus="1">
</div>
<blockquote>
<div>/gate/world/geometry/setXLength 30 um<br>
/gate/world/geometry/setYLength 30 um<br>
/gate/world/geometry/setZLength 30 um<br>
/gate/world/setMaterial Air<br>
<br>
/gate/world/daughters/name cell<br>
/gate/world/daughters/insert sphere<br>
/gate/cell/geometry/setRmin 0.0 mm<br>
/gate/cell/geometry/setRmax 10 um<br>
/gate/cell/placement/setTranslation 0 0 0 um<br>
/gate/cell/setMaterial G4_WATER<br>
/gate/cell/vis/setVisible 1<br>
/gate/cell/vis/setColor blue<br>
/gate/cell/vis/forceSolid false<br>
<br>
/gate/cell/daughters/name cluster<br>
/gate/cell/daughters/insert sphere<br>
/gate/cluster/geometry/setRmin 0.0 mm<br>
/gate/cluster/geometry/setRmax 500 nm<br>
/gate/cluster/placement/setTranslation 0 0 0
um<br>
/gate/cluster/setMaterial G4_WATER<br>
/gate/cluster/vis/setVisible 1<br>
/gate/cluster/vis/setColor yellow<br>
/gate/cluster/vis/forceSolid false<br>
<br>
/gate/cluster/daughters/name NP<br>
/gate/cluster/daughters/insert
sphere<br>
/gate/NP/geometry/setRmin 0 mm<br>
/gate/NP/geometry/setRmax 2 nm<br>
/gate/NP/placement/setTranslation 0 0
0 mm<br>
/gate/NP/setMaterial
Bismuth<br>
/gate/NP/vis/setVisible 1<br>
/gate/NP/vis/setColor red<br>
/gate/NP/vis/forceSolid true<br>
<br>
<br>
/gate/NP/repeaters/insert cubicArray<br>
/gate/NP/cubicArray/setRepeatNumberX 20<br>
/gate/NP/cubicArray/setRepeatNumberY 20<br>
/gate/NP/cubicArray/setRepeatNumberZ 20<br>
/gate/NP/cubicArray/setRepeatVector 5 5 5 nm<br>
<br>
/gate/cluster/repeaters/insert cubicArray<br>
/gate/cluster/cubicArray/setRepeatNumberX 3<br>
/gate/cluster/cubicArray/setRepeatNumberY 2<br>
/gate/cluster/cubicArray/setRepeatNumberZ 2<br>
/gate/cluster/cubicArray/setRepeatVector 5 5 5 um<br>
</div>
<div><br data-mce-bogus="1">
</div>
</blockquote>
<div>physics.mac<br data-mce-bogus="1">
</div>
<blockquote>
<div><br>
/gate/physics/addPhysicsList QGSP_BERT_EMV<br>
<br>
/gate/physics/Gamma/SetCutInRegion world 50 nm<br>
/gate/physics/Electron/SetCutInRegion world 50 nm<br>
<br>
/gate/physics/Gamma/SetCutInRegion cell 20 nm<br>
/gate/physics/Electron/SetCutInRegion cell 20 nm<br>
<br>
/gate/physics/Gamma/SetCutInRegion cluster 5 nm<br>
/gate/physics/Electron/SetCutInRegion cluster 5 nm<br>
<br>
/gate/physics/Gamma/SetCutInRegion NP 0.5 nm<br>
/gate/physics/Electron/SetCutInRegion NP 0.5 nm</div>
<div><br>
/gate/random/setEngineName MersenneTwister <br>
<br>
/gate/physics/addAtomDeexcitation<br>
<br>
/gate/physics/addProcess PhotoElectric
gamma<br>
/gate/physics/processes/PhotoElectric/setModel
LivermoreModel<br>
/gate/physics/processes/PhotoElectric/setAugerElectron
true<br>
/gate/physics/processes/PhotoElectric/setDeltaRayCut 10
eV<br>
/gate/physics/processes/PhotoElectric/setXRayCut 10
eV<br>
/gate/physics/processes/PhotoElectric/setStepFunction e-
0.2 0.1 nm<br>
<br>
/gate/physics/addProcess Compton gamma<br>
/gate/physics/processes/Compton/setModel
LivermoreModel<br>
/gate/physics/processes/Compton/setStepFunction e-
0.2 0.1 nm<br>
<br>
/gate/physics/addProcess
RayleighScattering<br>
/gate/physics/processes/RayleighScattering/setModel
LivermoreModel<br>
<br>
/gate/physics/addProcess GammaConversion<br>
/gate/physics/processes/GammaConversion/setModel
LivermoreModel<br>
<br>
/gate/physics/addProcess
ElectronIonisation<br>
/gate/physics/processes/ElectronIonisation/setModel
LivermoreModel e-<br>
/gate/physics/processes/ElectronIonisation/setModel
PenelopeModel e+<br>
/gate/physics/processes/ElectronIonisation/setStepFunction
e- 0.2 0.1 nm<br>
<br>
/gate/physics/addProcess Bremsstrahlung
e-<br>
/gate/physics/addProcess Bremsstrahlung
e+<br>
/gate/physics/processes/Bremsstrahlung/setModel
LivermoreModel e-<br>
/gate/physics/processes/Bremsstrahlung/setModel
PenelopeModel e+<br>
<br>
/gate/physics/addProcess
eMultipleScattering e-<br>
/gate/physics/addProcess
eMultipleScattering e+<br>
/gate/physics/processes/eMultipleScattering/setModel
UrbanModel e-<br>
/gate/physics/processes/eMultipleScattering/setModel
UrbanModel e+<br>
<br>
/gate/physics/SetMinKineticEnergyInRegion world
1 eV<br>
/gate/physics/SetMinKineticEnergyInRegion cell 1
eV<br>
/gate/physics/ActivateSpecialCuts gamma</div>
</blockquote>
<div><br data-mce-bogus="1">
</div>
<div>source.mac</div>
<blockquote>
<div><br>
/gate/source/addSource mybeam gps<br>
/gate/source/mybeam/gps/particle gamma<br>
/gate/source/mybeam/gps/pos/type Beam<br>
/gate/source/mybeam/gps/ene/mono 5 keV<br>
/gate/source/mybeam/gps/ene/type Gauss<br>
/gate/source/mybeam/gps/ene/sigma 1.0 keV<br>
/gate/source/mybeam/gps/direction 0 0 -1<br>
<br>
/gate/source/mybeam/gps/position 0 0 15 um<br>
/gate/source/mybeam/gps/pos/type Volume<br>
/gate/source/mybeam/gps/pos/shape Cylinder<br>
/gate/source/mybeam/gps/pos/inner_radius 0 mm<br>
/gate/source/mybeam/gps/pos/radius 10 um<br>
/gate/source/mybeam/gps/pos/halfz 0.01 um<br>
/gate/source/mybeam/gps/ang/type beam2d<br>
<br>
<br>
<br>
<br data-mce-bogus="1">
</div>
</blockquote>
<div><br data-mce-bogus="1">
</div>
<div><br data-mce-bogus="1">
</div>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-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>
<pre class="moz-signature" cols="72">--
-------------
DI Hermann Fuchs, PhD
Department of Radiation Oncology
Medical University Vienna
Währinger Gürtel 18-20
A-1090 Wien
Tel. + 43 / 2622 / 26100 / 628
Mail. <a class="moz-txt-link-abbreviated" href="mailto:hermann.fuchs@meduniwien.ac.at">hermann.fuchs@meduniwien.ac.at</a></pre>
</body>
</html>