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

Hermann Fuchs hermann.fuchs at meduniwien.ac.at
Thu Feb 24 15:11:56 CET 2022


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
> Gate-users at lists.opengatecollaboration.org
> 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.hermann.fuchs at meduniwien.ac.at
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220224/328ac292/attachment.html>


More information about the Gate-users mailing list