Is there a similar code like STATDOSE to analyze 3-dimensional dose distributions generated by C++ class library?
Thanks.
Note that the post authors are incorrect. Find our current active community here: https://www.reddit.com/r/EGSnrc
Is there a similar code like STATDOSE to analyze 3-dimensional dose distributions generated by C++ class library?
Thanks.
Hi everyone,
Can i combine two input files written in different platforms like using C++ and Fortran(Mortran) during BEAM simulation?Any suggestions??Thanks.

Dear EGSnrc users
I start to compile Dosxyznrc exemple using mac os yossimite
I have dis problem :
please could you help me?
I like to start a gui like dosxyznrc from the desktop. For this I installed on Ubuntu 16.04 gnome-panel and created a file Starter.desktop in the folder home/.local/share/applications:
#!/usr/bin/env xdg-open
[Desktop Entry]
Version=1.0
Type=Application
Terminal=false
Icon[de_DE]=gnome-panel-launcher
Name[de_DE]=Starter
Exec=gnome-desktop-item-edit .local/share/applications –create-new
Comment[de_DE]=Desktop-Symbol erstellen
Name=Starter
Comment=Desktop-Symbol erstellen
Icon=gnome-panel-launcher
Using this starter I tested it by creating a desktop-symbol for starting a second firefox as was shown on youtube and this was successful:
#!/usr/bin/env xdg-open
[Desktop Entry]
Version=1.0
Type=Application
Terminal=false
Icon[de_DE]=firefox
Name[de_DE]=My2Fx
Exec=firefox %u
Comment[de_DE]=second Firefox
Name=My2Fx
Comment=seconf Firefox
Icon=firefox
But using the starter for creating a desktop symbol for starting the dosxyznrc_gui resulted in the error message:
Beim Starten der Anwendung ist ein Fehler aufgetreten
(Running the application resulted in an error (my translation))
#!/usr/bin/env xdg-open
[Desktop Entry]
Version=1.0
Type=Application
Terminal=false
Icon[de_DE]=gnome-panel-launcher
Name[de_DE]=dosxyznrc
Exec=dosxyznrc_gui
Name=dosxyznrc
Icon=gnome-panel-launcher
The file dosxyznrc.desktop is executable but doesn´t work. I tried also with the option %u as is used for starting the second firefox but this didn´t help.
Thank you in advance for any help
Gunter

Hello every member
can you help me to modelate this geometry in dosxyznrc
GS_CDGeometry::howfar: ireg<0, but position appears inside
name=pinPointPTW base name=my_cd_planes inscribed name=baseCapaPMMA
x=(0.30156,-0.560007,0.0963128) u=(0.546865,-0.523682,0.653219) ibase=3 icd=1
distance to boundaries: base=(-1,2.22804) inscribed=(-1,0.153666)
**** Last 2 geometry steps:
old region=1, position=(0.267206,-0.527109,-33.4447), direction=(0.546865,-0.523682,0.653219) intended step=0.0628194, new region=1, step=0.0628194
old region=1, position=(0.30156,-0.560007,-33.4037), direction=(0.546865,-0.523682,0.653219) intended step=0.00688051, new region=1, step=0.00688051
The above error occurs frequently when I run the following geometry in egs_chamber.
*******************************************************************************
:start geometry definition:
################################################################################
################################################################################
###MUNDO
##Testada e Funcional
:start geometry:
library = egs_box
name = world
box size = 50 50 58.1996
:start media input:
media = AIR700ICRU
:stop media input:
:start transformation:
translation = 0 0 -28.1996
:stop transformation:
:stop geometry:
################################################################################
################################################################################
###FANTOMA ÁGUA
##Testada e Funcional
:start geometry:
library = egs_box
name = phantom
box size = 30 30 30
:start media input:
media = H2O700ICRU
:stop media input:
:start transformation:
translation = 0 0 -43.1996
:stop transformation:
:stop geometry:
################################################################################
################################################################################
###Definição da geometria para XCSE
##Testada e funcional
:start geometry:
library = egs_box
name = boxXCSE
box size = 5 1 29.5
:start media input:
media = H2O700ICRU
:stop media input:
:start transformation:
translation = 0 0 -43.1996
:stop transformation:
:stop geometry:
######################################################################
###Definição da geometria para IPSS
##Testada e funcional
:start geometry:
library = egs_box
name = boxIPSS
box size = 4.9 0.9 29.
:start media input:
media = H2O700ICRU
:stop media input:
:start transformation:
translation = 0 0 -43.1997
:stop transformation:
:stop geometry:
######################################################################
###Envelopamento XCSE + IPSS
##Testada e funcional
:start geometry:
library = egs_genvelope
name = GeoIPSSXCSE #(Vol = 0.00418879 cm³)
base geometry = boxXCSE
inscribed geometries = boxIPSS
:stop geometry:
######################################################################
######################################################################
### DETECTOR PADRÃO
##Testada e Funcional
:start geometry:
name = my_cd_planes
library = egs_planes
type = EGS_Xplanes
positions = -0.53 0 0.13 0.298 1.52 #1.8
# No media required
:stop geometry:
# Esferas concêntricas envolvendo o volume sensível de ar, o
# parede de grafite e a capa de PMMA
:start geometry:
name = my_cd_sphere
library = egs_spheres
midpoint = 0 0 0
radii = 0.145 0.2 0.53
:start media input:
media = AIR700ICRU 226C700ICRU PMMA700ICRU
set medium = 1 1
set medium = 2 2
set medium = 3 3
:stop media input:
:stop geometry:
:start geometry:
name = base_calota #bc
library = egs_cylinders
type = EGS_XCylinders
radii = 0.03 0.13 0.2 0.53
:start media input:
media = AL700ICRU AIR700ICRU 226C700ICRU PMMA700ICRU
set medium = 1 1
set medium = 2 2
set medium = 3 3
set medium = 4 4
:stop media input:
:stop geometry:
:start geometry:
name = base_bc #bbc
library = egs_cylinders
type = EGS_XCylinders
radii = 0.03 0.145 0.53
:start media input:
media = AL700ICRU 226C700ICRU PMMA700ICRU
set medium = 1 1
set medium = 2 2
set medium = 3 3
:stop media input:
:stop geometry:
:start geometry:
name = baseCapaPMMA
library = egs_cylinders
type = EGS_XCylinders
radii = 0.03 0.67
:start media input:
media = AL700ICRU PMMA700ICRU
set medium = 1 1
set medium = 2 2
:stop media input:
:stop geometry:
# The composite geometry
:start geometry:
name = pinPointPTW
library = egs_cdgeometry
base geometry = my_cd_planes
# set geometry = 1 geom means:
# “in region 1 of the basegeometry, use geometry “geom”
set geometry = 0 my_cd_sphere
set geometry = 1 base_calota
set geometry = 2 base_bc
set geometry = 3 baseCapaPMMA
#set geometry = 4 base_geral
# The final region numbers are attributed by the cd geometry object;
# Use the viewer to determine region numbers
:stop geometry:
######################################################################
######################################################################
### World + Phantom + GeoIPSSXCSE
##Testada e Funcional
:start geometry:
library = egs_genvelope
name = world_phantom #(Vol = 0.00418879 cm³)
base geometry = world
inscribed geometries = phantom
:stop geometry:
:start geometry:
library = egs_genvelope
name = completWorld #(Vol = 0.00418879 cm³)
base geometry = world_phantom
inscribed geometries = GeoIPSSXCSE
:stop geometry:
################################################################################
################################################################################
###Detectores Esféricos + completWorld
##Testada e funcional
:start input loop:
loop count = 40
loop variable = 1 posicionadorZ -29.0 -0.5 #-28.2496 -0.1
loop variable = 0 numerador 1 1
:start geometry:
library = egs_gtransformed
my geometry = pinPointPTW
name = Z$(numerador)
:start transformation:
translation = 0 0 $(posicionadorZ)
:stop transformation:
:stop geometry:
:stop input loop:
:start input loop:
loop count = 40
loop variable = 0 numerador 1 1
:start geometry:
library = egs_genvelope
name = dose_Z$(numerador)
base geometry = world_phantom
inscribed geometries = Z$(numerador)
:stop geometry:
:stop input loop:
simulation geometry = dose_Z1
:stop geometry definition:
######################################################################
#DEFINIÇÃO DA FONTE
:start source definition:
:start source:
name = HEAD
library = egs_phsp_source
phase space file = C2100Cone21mm_0p13.egsphsp2
particle type = all
#cutout = -1 1 -2 2
reuse photons = 10
reuse electrons = 10
:stop source:
:start source:
library = egs_transformed_source
name = HEADSOURCE
source name = HEAD
nrcycl = 0
:start transformation:
translation = 0 0 0.01
rotation = 0 3.14159
:stop transformation:
:stop source:
simulation source = HEADSOURCE
:stop source definition:
################################################################################
# CONTROLE DA SIMULAÇÃO
:start run control:
ncase = 1000 #100000000
nbatch = 10
calculation = first
geometry error limit = 100
:stop run control:
################################################################################
#SCORING PARÂMETROS
:start scoring options: #ok
calculation type = dose
silent = 0
:start calculation geometry: #ok
geometry name = world_phantom #completWorld
cavity regions = 1
cavity mass = 127.89
cavity geometry = phantom #GeoIPSSXCSE
enhance regions = 1
enhancement = 128
:stop calculation geometry: #ok
:start input loop: #ok
loop count = 40
loop variable = 0 numerador 1 1
:start calculation geometry: #ok
geometry name = dose_Z$(numerador)
cavity regions = 12 13
cavity mass = 0.000020685
cavity geometry = Z$(numerador)
enhance regions = 1
enhancement = 128
:stop calculation geometry: #ok
:stop input loop: #ok
:stop scoring options: #ok
################################################################################
#TRACK SCORING
:start ausgab object definition:
:start ausgab object:
library = egs_track_scoring
score photons = yes
score electrons = yes
score positrons = yes
:stop ausgab object:
:stop ausgab object definition:
:start variance reduction:
cs enhancement = 1
TmpPhsp = 1 # i.e., score phase space and use it once in each specified
# geometry
:start range rejection:
rejection = 512
Esave = 0.512
cavity geometry = GeoIPSSXCSE #this is a dummy
rejection range medium = H2O700ICRU
:stop range rejection:
:stop variance reduction:
:start MC transport parameter:
Global ECUT = 0.521 # Global electron transport cutoff
Global PCUT = 0.001 # Global photon transport cutoff
Set ECUT 0 0 = 6.2
Set PCUT 0 0 = 6.2
:stop MC transport parameter:
******************************************************************************
Could someone help me?
Hello everyone,
I’m new to EGSnrc and wonder if I miss something. I recently buildt a beowulf cluster for running EGSnrc simulations. I managed to successfully run some of the example files but I run into a problem when calculating egs_chamber inputs in parallel. The exb command starts several caclulations, but no output whatsosever appears, not even an error file is generated. When I try to look them up they don’t appear in the queue anymore, so I guess they end very quick or never even get started right. so far ai onyl had that Problemin egs_chamber. Does anyone by chance have an idea what I should look into now to solve that?
Hi,
Can someone please share the process of simulating a VMAT plan with DOSXYZnrc using the DICOM plan from the treatment planning system?
Or share a sample input file?
Thank you
Sotiri Stathakis
Just running a test DOSXYZnrc simulation using (an IAEA phase space file) in Windows 10. Getting the following error (even for 1500 histories):
* Warning * Used all particles fromsource file
Restarting from first particle infile
Strangely the same input file on windows 7 using an earlier version of EGS is running fine. Any one has encountered any such issue. Otherwise both BEAMnrc and DOSXYZnrc are running fine for other simulations
Thanks,
Basim
Hello everyone,
If i modify the pre-existing CMs during beam simulation, i edit/modify CMNAME_macros.mortran and CMNAME_cm.mortran. After that what else i need to do? is this enough to run the simulation with modified CM?Is there any better way to create new CM?Help will be much appreciated!! Thanks.