Hello

Hello,

I have a question about the IAEA phase space files. I could not find information about different files for same energy, field size and LINAC. There are three name formats:

1. Ending with part1, part2, part3, part 4

2. Ending with w1,w2, w3

3. Ending with 01, 02, 03, 04…

I guess part1, part 2… and 01,02,…. these two formats show the parts of same phase space file and need to be combined by addphsp before using for simulation. But for w1,w2,w3…I got no idea.

I have checked the header files but could not figure out…

Please give some hint about use of these phase space files if there is anything specific?

I am adding screen shorts from the IAEA website page for Phase Space Data-base….

Thanks!

I try to run sprrznrc and got error code below.

I try to run sprrznrc and got error code below.

If you are not trying to reset transport parameters,

ignore all the output until the message

****************** end input transport parameter *********************

**************ERROR**************

>>Radiative Compton corrections<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Pair cross sections<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Triplet production<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>scale L EII cross-sections<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Compton cross sections<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Electric Field<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Magnetic Field<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>EM ESTEPE<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Photonuclear attenuation<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

**************ERROR**************

>>Photonuclear cross sections<< NOT FOUND

END OF DELIMETER: MC TRANSPORT PARAMETER

****************** end input transport parameter *********************

Please recommend me what is the error? and introduce me how to fix this problem.

and This below, is my code

TITLE= sprrznrc_19mm_graphiteCavityChamber

##########################

:start I/O control:

IWATCH= off

STORE INITIAL RANDOM NUMBERS= no

IRESTART= first

STORE DATA ARRAYS= yes

SPR OUTPUT= slabs/cylinders

SPR IN CYLINDER IX= 1, 2

SPR IN SLAB IZ= 0

:stop I/O control:

#########################

##########################

:start Monte Carlo inputs:

NUMBER OF HISTORIES= 1000000

INITIAL RANDOM NO. SEEDS= 1, 3

MAX CPU HOURS ALLOWED= 90.000

IFULL= dose and stoppers

STATISTICAL ACCURACY SOUGHT= 0.0000

PHOTON REGENERATION= yes

:stop Monte Carlo inputs:

#########################

##########################

:start geometrical inputs:

METHOD OF INPUT= individual

Z OF FRONT FACE= 0.

DEPTH BOUNDARIES= 0.3, 0.64, 1.6, 1.9

RADII= 0.3, 0.75, 0.85, 0.95

MEDIA= C87ICRU521,

AIR521ICRU;

DESCRIPTION BY= planes

MEDNUM= 1, 2, 2, 2

START ZSLAB= 1, 2, 2, 2

STOP ZSLAB= 4, 3, 4, 2

START RING= 1, 2, 3, 4

STOP RING= 1, 2, 3, 4

:stop geometrical inputs:

#########################

##########################

:start source inputs:

INCIDENT PARTICLE= photon

SOURCE NUMBER= 0

SOURCE OPTIONS= 5, 0, 0, 5

INCIDENT ENERGY= spectrum

SPEC FILENAME= D:EGSnrcHEN_HOUSEspectraegsnrccs137.spectrum

SPEC IOUTSP= include

:stop source inputs:

#########################

##########################

:start MC transport parameter:

Global ECUT= 0.521

Global PCUT= 0.01

Global SMAX= 1e10

ESTEPE= 0.25

XImax= 0.5

Skin depth for BCA= 3

Boundary crossing algorithm= EXACT

Electron-step algorithm= PRESTA-II

Spin effects= on

Brems angular sampling= KM

Brems cross sections= BH

Photon cross sections= PEGS4

Electron Impact Ionization= Off

Photon cross-sections output= on

Bound Compton scattering= On

Pair angular sampling= Simple

Photoelectron angular sampling= On

Rayleigh scattering= On

Atomic relaxations= On

Set PCUT= 0

Set PCUT start region= 1

Set PCUT stop region= 1

Set ECUT= 0

Set ECUT start region= 1

Set ECUT stop region= 1

Set SMAX= 0

Set SMAX start region= 1

Set SMAX stop region= 1

:stop MC transport parameter:

#########################

##########################

:start media definition:

AE= 0.521

UE= 50.511

AP= 0.01

UP= 50

material data file= D:/EGSnrc/HEN_HOUSE/pegs4/data/material.dat

:stop media definition:

#########################

##########################

:start variance reduction:

ELECTRON RANGE REJECTION= off

ESAVEIN= 0.0

RUSSIAN ROULETTE DEPTH= 0.0000

RUSSIAN ROULETTE FRACTION= 0.0000

EXPONENTIAL TRANSFORM C= 0.0000

PHOTON FORCING= on

START FORCING= 1

STOP FORCING AFTER= 1

:stop variance reduction:

#########################

Hello all

Hello all,

Just wondering if anyone had success running an IAEA-phase space file? I’m using source 20 in DOSXYZnrc and trying to run a VMAT (so dynamic mlc’s and jaws throughout the run) simulation. I’ve successfully run simulations using an .egsphsp phase space file generated by my accelerator, but when I try running a simulation using an IAEA-phase space file, the process quits after stating “Performing BEAM calibration run…Running 1000000 particles to calibrate BEAM accelerator”. Any tips on where I might have went wrong?

Hi, can someone please help me with this error in egs_brachy. Marc Chamberland

Hi, can someone please help me with this error in egs_brachy. Marc Chamberland

EGS_CDGeometry::howfar: ireg<0, but position appears inside

name=seed base name=planes inscribed name=seed_center

x=(-0.0363666,0.0162328,-0.0391393) u=(0.305875,-0.166924,-0.937324) ibase=1 icd=1

distance to boundaries: base=(0,0.158281) inscribed=(0,0.0196619)

I am using an egsphant phantom and the OncoSeed_6711. Is_convex is set to false.

Thank you

Hello, I,m getting this error message :

Hello, I,m getting this error message :

****************************************

In EGS_Application::initSimulation()

VHP_OrganData: wrong endianess.

This machine has endianess 1, data was written on 0

This requires byte swaping, which has not been implemented yet

EGS_VHPGeometry: failed to construct organ data

createGeometry(EGS_VHPGeometry): failed to construct geometrycreateSingleGeometry: got null geometry

library: egs_vhp_geometry

input:

:START geometry:

library = egs_vhp_geometry

name = fax06

slice range = 0 180

phantom data = /home/ismail/EGSnrc/egs_home

media data = /home/ismail/EGSnrc/egs_home

:STOP geometry:

********************************************

Any help?

Hello

Hello,

I just started working with EGSnrc. I want to simulate a photon source with a gussian shape and I need PDD und profile curves in water phantom of this source. Is it possible to simulate a gaussian shaped source in DOSXYZ? Or is their another way to get the curves I need?

I would be very appreciated if someone can help me.

Best Regards.