Open CQUswl opened 1 year ago
This seems to only be an issue with tungsten from 521ICRU and howfarless turned on. If you turn off howfarless or change materials, it seems to be fine. This will be a tricky issue to investigate, so for now use the following, where we switch to pegsless
mode:
test #!GUI1.0
1
tungsten
0.521, 0.01, 0, 0
-1, -1, -1, 1
-0.1
0.004, 50
-0.1
0.004, 50
0
0.004, 50
0, 50, 0, 50, 0, 50, , 19.35
0, 0, 0, 0, 0, 0, 0, 0
0, 0, 0, 0, 0, 0, 0, 0
0, 100, 0, 100, 0, 100, 1, 0
0, 0, 0, 0, 0, 0, 0, 0
-1, 0, -0.01, 0.01, -0.01, 0.01, 90, 90, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0
0.2
20000000, 0, 99, 33, 97, 100.0, 0, 0, 1, 0, , 0, 0, 0, 1, 1, 0, 0
#########################
:Start MC Transport Parameter:
Global ECUT= 0.521
Global PCUT= 0.01
Global SMAX= 5
ESTEPE= 0.25
XIMAX= 0.5
Boundary crossing algorithm= PRESTA-I
Skin depth for BCA= 0
Electron-step algorithm= PRESTA-II
Spin effects= On
Brems angular sampling= Simple
Brems cross sections= BH
Bound Compton scattering= Off
Compton cross sections= default
Pair angular sampling= Simple
Pair cross sections= BH
Photoelectron angular sampling= Off
Rayleigh scattering= Off
Atomic relaxations= Off
Electron impact ionization= Off
Photon cross sections= xcom
Photon cross-sections output= Off
:Stop MC Transport Parameter:
#########################
:start media definition:
### particle production thresholds
ae = 0.521 # MeV
ap = 0.010 # MeV
ue = 50.511 # MeV
up = 50 # MeV
:start tungsten:
density correction file = tungsten
:stop tungsten:
### See $HEN_HOUSE/pegs4/density_corrections/ for more materials, or create
### you own: https://github.com/nrc-cnrc/EGSnrc/wiki/Generating-density-correction-files
:stop media definition:
This seems to only be an issue with tungsten from 521ICRU and howfarless turned on. If you turn off howfarless or change materials, it seems to be fine. This will be a tricky issue to investigate, so for now use the following, where we switch to
pegsless
mode:test #!GUI1.0 1 tungsten 0.521, 0.01, 0, 0 -1, -1, -1, 1 -0.1 0.004, 50 -0.1 0.004, 50 0 0.004, 50 0, 50, 0, 50, 0, 50, , 19.35 0, 0, 0, 0, 0, 0, 0, 0 0, 0, 0, 0, 0, 0, 0, 0 0, 100, 0, 100, 0, 100, 1, 0 0, 0, 0, 0, 0, 0, 0, 0 -1, 0, -0.01, 0.01, -0.01, 0.01, 90, 90, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 0.2 20000000, 0, 99, 33, 97, 100.0, 0, 0, 1, 0, , 0, 0, 0, 1, 1, 0, 0 ######################### :Start MC Transport Parameter: Global ECUT= 0.521 Global PCUT= 0.01 Global SMAX= 5 ESTEPE= 0.25 XIMAX= 0.5 Boundary crossing algorithm= PRESTA-I Skin depth for BCA= 0 Electron-step algorithm= PRESTA-II Spin effects= On Brems angular sampling= Simple Brems cross sections= BH Bound Compton scattering= Off Compton cross sections= default Pair angular sampling= Simple Pair cross sections= BH Photoelectron angular sampling= Off Rayleigh scattering= Off Atomic relaxations= Off Electron impact ionization= Off Photon cross sections= xcom Photon cross-sections output= Off :Stop MC Transport Parameter: ######################### :start media definition: ### particle production thresholds ae = 0.521 # MeV ap = 0.010 # MeV ue = 50.511 # MeV up = 50 # MeV :start tungsten: density correction file = tungsten :stop tungsten: ### See $HEN_HOUSE/pegs4/density_corrections/ for more materials, or create ### you own: https://github.com/nrc-cnrc/EGSnrc/wiki/Generating-density-correction-files :stop media definition:
This seems to only be an issue with tungsten from 521ICRU and howfarless turned on. If you turn off howfarless or change materials, it seems to be fine. This will be a tricky issue to investigate, so for now use the following, where we switch to
pegsless
mode:test #!GUI1.0 1 tungsten 0.521, 0.01, 0, 0 -1, -1, -1, 1 -0.1 0.004, 50 -0.1 0.004, 50 0 0.004, 50 0, 50, 0, 50, 0, 50, , 19.35 0, 0, 0, 0, 0, 0, 0, 0 0, 0, 0, 0, 0, 0, 0, 0 0, 100, 0, 100, 0, 100, 1, 0 0, 0, 0, 0, 0, 0, 0, 0 -1, 0, -0.01, 0.01, -0.01, 0.01, 90, 90, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 0.2 20000000, 0, 99, 33, 97, 100.0, 0, 0, 1, 0, , 0, 0, 0, 1, 1, 0, 0 ######################### :Start MC Transport Parameter: Global ECUT= 0.521 Global PCUT= 0.01 Global SMAX= 5 ESTEPE= 0.25 XIMAX= 0.5 Boundary crossing algorithm= PRESTA-I Skin depth for BCA= 0 Electron-step algorithm= PRESTA-II Spin effects= On Brems angular sampling= Simple Brems cross sections= BH Bound Compton scattering= Off Compton cross sections= default Pair angular sampling= Simple Pair cross sections= BH Photoelectron angular sampling= Off Rayleigh scattering= Off Atomic relaxations= Off Electron impact ionization= Off Photon cross sections= xcom Photon cross-sections output= Off :Stop MC Transport Parameter: ######################### :start media definition: ### particle production thresholds ae = 0.521 # MeV ap = 0.010 # MeV ue = 50.511 # MeV up = 50 # MeV :start tungsten: density correction file = tungsten :stop tungsten: ### See $HEN_HOUSE/pegs4/density_corrections/ for more materials, or create ### you own: https://github.com/nrc-cnrc/EGSnrc/wiki/Generating-density-correction-files :stop media definition:
I understand, thank you.
When running DOSXYZnrc GUI, if the number of histories exceeds 10000000, the error shown above will be displayed. The following are the input files: test #!GUI1.0 1 W521ICRU 0.521, 0.01, 0, 0 -1, -1, -1, 1 -0.1 0.004, 50 -0.1 0.004, 50 0 0.004, 50 0, 100, 0, 100, 0, 100, 1, 19.35 0, 0, 0, 0, 0, 0, 0, 0 0, 0, 0, 0, 0, 0, 0, 0 0, 100, 0, 100, 0, 100, 1, 0 0, 0, 0, 0, 0, 0, 0, 0 -1, 0, -0.01, 0.01, -0.01, 0.01, 90, 90, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 0.2 20000000, 0, 99, 33, 97, 100.0, 0, 0, 1, 0, , 0, 0, 0, 1, 1, 0 ######################### :Start MC Transport Parameter:
Global ECUT= 0.521 Global PCUT= 0.01 Global SMAX= 5 ESTEPE= 0.25 XIMAX= 0.5 Boundary crossing algorithm= PRESTA-I Skin depth for BCA= 0 Electron-step algorithm= PRESTA-II Spin effects= On Brems angular sampling= Simple Brems cross sections= BH Bound Compton scattering= Off Compton cross sections= default Pair angular sampling= Simple Pair cross sections= BH Photoelectron angular sampling= Off Rayleigh scattering= Off Atomic relaxations= Off Electron impact ionization= Off Photon cross sections= xcom Photon cross-sections output= Off
:Stop MC Transport Parameter: #########################