[Gate-users] Problem with 10^10 primaries send

[Sarah Blind]

Many thanks Hermann for your response! I will try this. 

Sarah 

----- Le 24 Fév 22, à 15:11, Hermann Fuchs <hermann.fuchs at meduniwien.ac.at> a écrit : 

> Hi Sarah,

> I believe this could be a problem with the implemented datatype. The number of
> particles is stored internally in GATE as an integer.

> The highest number an integer can store is~2.1e9

> Consequently simulating 1e9 particles works and 1e10 particles will not work any
> more.

> The easiest solution is to split your simulation in multiple parts, each having
> less particles than 2.1e9.

> All the best,

> Hermann
> On 23.02.22 16:37, Sarah Blind wrote:

>> Dear Gate users,

>> I have the following model : 1 cell inside of which there are 12 clusters of
>> 8000 nanoparticles each.
>> My source is composed of a gamma beam of 5 keV.
>> My simulation works fine up to 10^9 photons send. But with 10^10 primaries
>> nothing happens anymore (there is no error message).
>> Does someone have an idea what the problem might be ?

>> Thank you very much for your help,
>> Regards,
>> Sarah Blind

>> world.mac

>>> /gate/world/geometry/setXLength 30 um
>>> /gate/world/geometry/setYLength 30 um
>>> /gate/world/geometry/setZLength 30 um
>>> /gate/world/setMaterial Air

>>> /gate/world/daughters/name cell
>>> /gate/world/daughters/insert sphere
>>> /gate/cell/geometry/setRmin 0.0 mm
>>> /gate/cell/geometry/setRmax 10 um
>>> /gate/cell/placement/setTranslation 0 0 0 um
>>> /gate/cell/setMaterial G4_WATER
>>> /gate/cell/vis/setVisible 1
>>> /gate/cell/vis/setColor blue
>>> /gate/cell/vis/forceSolid false

>>> /gate/cell/daughters/name cluster
>>> /gate/cell/daughters/insert sphere
>>> /gate/cluster/geometry/setRmin 0.0 mm
>>> /gate/cluster/geometry/setRmax 500 nm
>>> /gate/cluster/placement/setTranslation 0 0 0 um
>>> /gate/cluster/setMaterial G4_WATER
>>> /gate/cluster/vis/setVisible 1
>>> /gate/cluster/vis/setColor yellow
>>> /gate/cluster/vis/forceSolid false

>>> /gate/cluster/daughters/name NP
>>> /gate/cluster/daughters/insert sphere
>>> /gate/NP/geometry/setRmin 0 mm
>>> /gate/NP/geometry/setRmax 2 nm
>>> /gate/NP/placement/setTranslation 0 0 0 mm
>>> /gate/NP/setMaterial Bismuth
>>> /gate/NP/vis/setVisible 1
>>> /gate/NP/vis/setColor red
>>> /gate/NP/vis/forceSolid true

>>> /gate/NP/repeaters/insert cubicArray
>>> /gate/NP/cubicArray/setRepeatNumberX 20
>>> /gate/NP/cubicArray/setRepeatNumberY 20
>>> /gate/NP/cubicArray/setRepeatNumberZ 20
>>> /gate/NP/cubicArray/setRepeatVector 5 5 5 nm

>>> /gate/cluster/repeaters/insert cubicArray
>>> /gate/cluster/cubicArray/setRepeatNumberX 3
>>> /gate/cluster/cubicArray/setRepeatNumberY 2
>>> /gate/cluster/cubicArray/setRepeatNumberZ 2
>>> /gate/cluster/cubicArray/setRepeatVector 5 5 5 um

>> physics.mac

>>> /gate/physics/addPhysicsList QGSP_BERT_EMV

>>> /gate/physics/Gamma/SetCutInRegion world 50 nm
>>> /gate/physics/Electron/SetCutInRegion world 50 nm

>>> /gate/physics/Gamma/SetCutInRegion cell 20 nm
>>> /gate/physics/Electron/SetCutInRegion cell 20 nm

>>> /gate/physics/Gamma/SetCutInRegion cluster 5 nm
>>> /gate/physics/Electron/SetCutInRegion cluster 5 nm

>>> /gate/physics/Gamma/SetCutInRegion NP 0.5 nm
>>> /gate/physics/Electron/SetCutInRegion NP 0.5 nm

>>> /gate/random/setEngineName MersenneTwister

>>> /gate/physics/addAtomDeexcitation

>>> /gate/physics/addProcess PhotoElectric gamma
>>> /gate/physics/processes/PhotoElectric/setModel LivermoreModel
>>> /gate/physics/processes/PhotoElectric/setAugerElectron true
>>> /gate/physics/processes/PhotoElectric/setDeltaRayCut 10 eV
>>> /gate/physics/processes/PhotoElectric/setXRayCut 10 eV
>>> /gate/physics/processes/PhotoElectric/setStepFunction e- 0.2 0.1 nm

>>> /gate/physics/addProcess Compton gamma
>>> /gate/physics/processes/Compton/setModel LivermoreModel
>>> /gate/physics/processes/Compton/setStepFunction e- 0.2 0.1 nm

>>> /gate/physics/addProcess RayleighScattering
>>> /gate/physics/processes/RayleighScattering/setModel LivermoreModel

>>> /gate/physics/addProcess GammaConversion
>>> /gate/physics/processes/GammaConversion/setModel LivermoreModel

>>> /gate/physics/addProcess ElectronIonisation
>>> /gate/physics/processes/ElectronIonisation/setModel LivermoreModel e-
>>> /gate/physics/processes/ElectronIonisation/setModel PenelopeModel e+
>>> /gate/physics/processes/ElectronIonisation/setStepFunction e- 0.2 0.1 nm

>>> /gate/physics/addProcess Bremsstrahlung e-
>>> /gate/physics/addProcess Bremsstrahlung e+
>>> /gate/physics/processes/Bremsstrahlung/setModel LivermoreModel e-
>>> /gate/physics/processes/Bremsstrahlung/setModel PenelopeModel e+

>>> /gate/physics/addProcess eMultipleScattering e-
>>> /gate/physics/addProcess eMultipleScattering e+
>>> /gate/physics/processes/eMultipleScattering/setModel UrbanModel e-
>>> /gate/physics/processes/eMultipleScattering/setModel UrbanModel e+

>>> /gate/physics/SetMinKineticEnergyInRegion world 1 eV
>>> /gate/physics/SetMinKineticEnergyInRegion cell 1 eV
>>> /gate/physics/ActivateSpecialCuts gamma

>> source.mac

>>> /gate/source/addSource mybeam gps
>>> /gate/source/mybeam/gps/particle gamma
>>> /gate/source/mybeam/gps/pos/type Beam
>>> /gate/source/mybeam/gps/ene/mono 5 keV
>>> /gate/source/mybeam/gps/ene/type Gauss
>>> /gate/source/mybeam/gps/ene/sigma 1.0 keV
>>> /gate/source/mybeam/gps/direction 0 0 -1

>>> /gate/source/mybeam/gps/position 0 0 15 um
>>> /gate/source/mybeam/gps/pos/type Volume
>>> /gate/source/mybeam/gps/pos/shape Cylinder
>>> /gate/source/mybeam/gps/pos/inner_radius 0 mm
>>> /gate/source/mybeam/gps/pos/radius 10 um
>>> /gate/source/mybeam/gps/pos/halfz 0.01 um
>>> /gate/source/mybeam/gps/ang/type beam2d

>> _______________________________________________
>> Gate-users mailing list [ mailto:Gate-users at lists.opengatecollaboration.org |
>> Gate-users at lists.opengatecollaboration.org ] [
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users |
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users ]

> --
> -------------
> 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. [ mailto:hermann.fuchs at meduniwien.ac.at | hermann.fuchs at meduniwien.ac.at ]

> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220224/182305eb/attachment-0001.html>