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

[Sarah Blind]

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



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