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

Mon Mar 7 16:52:35 CET 2022

Dear Sarah,

there has been very recently a thread in the Gate user list on exactly this question, search for Running Gate on a Cluster (Mijat Paunovic) (14 and 15 February 2022). Basically, if your simulation setup is not time dependent (e.g. no rotating geometry), then the basic idea is to simply run a simulation multiple times and combine the output.

It depends a bit on whether you have multiple CPU cores available on your machine (which is usually the case these days), allowing you to execute runs in parallel, or whether you only have access to one CPU core, in which case you then repeat your simulation with this core multiple times.

Best regards,
David

On 7 Mar2022, at 14:47, Sarah Blind <sarah.blind at univ-lorraine.fr<mailto:sarah.blind at univ-lorraine.fr>> wrote:

Dear Dimitris Thanasas,

I come back to you about launching GATE with more than 10^9 primaries.
I have tried to set an activity to my source and to run the simulation for a certain time as you suggested.
It allowed me to launch my simulation with 10^10 photons, but not more.

Now I want to split my simulation, but I don't now how exactly and I can't find any information about it, would you give me some advices ?

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

----- Le 25 Fév 22, à 7:59, Dimitrios Thanasas <dthanas at uoa.gr<mailto:dthanas at uoa.gr>> a écrit :
Dear Sarah,

You could try to give an activity to your source, (e.g. 1 GBq - without
half life) and then to run the simulation for a certain time (e.g. 10s). I
think that this wiil work (you may need a statactor to know exactly the
number of the emitted particles).
If you try to split the simulation, you may need to change the engine seed
in order not to have the same results in each simulation.

Regards,
Dimitris

> 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
>
> BQ_BEGIN
>
>
> /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
>
> BQ_END
>
>
> source.mac
>
> BQ_BEGIN
>
>
> /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
>
>
>
>
>
> BQ_END
>
>
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users<https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.opengatecollaboration.org_mailman_listinfo_gate-2Dusers&d=DwQFaQ&c=XYzUhXBD2cD-CornpT4QE19xOJBbRy-TBPLK0X9U2o8&r=lfXR69GFfS7NT-Wp5HZqbtBDrbYoilDdmNEG7fKW7aM&m=6XIuSTRo-YlrPIDJvuacVBO_vPHHA9U1ntiFKhR5nOk&s=YrpArz4E5R17m2t-I5odtmrXiJM1ENqGmygg_6iw0Xg&e=>

--
Sarah Blind
Ingénieure de recherche - CNRS
Centre de Recherche en Automatique de Nancy (CRAN)
Département Biologie, Signaux et Systèmes en Cancérologie et Neurosciences (BioSiS)
Site de la faculté de médecine
9, Avenue de la Forêt de Haye
54505 Vandoeuvre les Nancy
Tél : 03.72.74.61.13
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
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=6XIuSTRo-YlrPIDJvuacVBO_vPHHA9U1ntiFKhR5nOk&s=YrpArz4E5R17m2t-I5odtmrXiJM1ENqGmygg_6iw0Xg&e=

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20220307/8bea59af/attachment-0001.html>