Trouble with simulation getting stuck in batch 1

[peytonlalain]

Hi! I am working on a simulation in EGSnrc of a planar Y-90 disc source in an extrapolation chamber where I am aiming to score the dose in the collector of my chamber's electrode. I have built the geometry and in egs_view everything looks correct, and I can't see any obvious issues with my source definition or my pegs4dat file, but when I try to run the simulation, it gets perpetually stuck in batch 1. I've tried running the simulation replacing my Y-90 source with a simple point source and it ran without issue, which leads me to think the geometry isn't an issue, but I'm not sure what else the problem could be. Below I have put the script I am running, which I currently have very simplified to identify where it is getting stuck, and I only have set to ncase of 1000 just to see if it will run. For reference, I am attempting to run in egs_chamber. I am also attaching my script as an egsinp file, as well as the pegs4dat file I am using. If anyone has any ideas for what could be wrong that I am missing, I am open to any suggestions. Thank you!

#############################################################################
## EGSnrc egs_chamber Input File: Y-90 Ophthalmic Brachytherapy
#############################################################################

:start geometry definition:

    ## 1. Y-90 Source Geometry (Stack: Taper -> Ring -> Disc -> Cap)
    :start geometry:
        name = Y90_disc_source
        library = egs_cones
        type = EGS_ConeStack
        axis = 0 0 -0.1 0 0 1

        :start layer: # Bottom taper
            thickness = 0.025 
            top radii = 0.27 3.5
            bottom radii = 0.30 3.5
            media = titanium aluminum
        :stop layer:

        :start layer: # Side wall
            thickness = 0.025
            top radii = 0.27 0.30 3.5
            bottom radii = 0.27 0.30 3.5
            media = air titanium aluminum
        :stop layer:

        :start layer: # Yttrium source disc
            thickness = 0.025
            top radii = 0.27 0.30 3.5
            bottom radii = 0.27 0.30 3.5
            media = yttrium titanium aluminum
        :stop layer:

        :start layer: # Titanium flat top cap
            thickness = 0.025 
            top radii = 0.30 3.5
            bottom radii = 0.30 3.5
            media = titanium aluminum
        :stop layer:
    :stop geometry:

    ## 2. Electrode Geometry (Collector + Base)
    :start geometry:
        name = electrode
        library = egs_cones
        type = EGS_ConeStack
        axis = 0 0 0.01 0 0 1

        :start layer: # Collector and Guard (C552)
            thickness = 0.635
            top radii = 0.20 0.508
            bottom radii = 0.20 0.508
            media = C552 C552
        :stop layer:

        :start layer: # Electrode base
            thickness = 2.177
            top radii = 1.50
            bottom radii = 1.50
            media = acrylic
        :stop layer:
    :stop geometry:
    
   ## 3. Source Holder
	:start geometry:
    	name = source_holder
    	library = egs_cones
    	type = EGS_ConeStack
    	axis = 0 0 -1.1001 0 0 1  # Shifted down by 0.0001 cm
    	:start layer:
        	thickness = 1.0
        	top radii = 3.50
        	bottom radii = 3.50
        	media = aluminum
    	:stop layer:
	:stop geometry:

	## 4. Air Box (Global Container)
    :start geometry:
        name = air_box
        library = egs_box
        box size = 10 10 10
        # SHIFT THE BOX SO (0,0,0) IS IN THE MIDDLE
        #origin = -5 -5 -5 
        :start media input:
            media = air
        :stop media input:
    :stop geometry:

    ## 5. Master Model
    :start geometry:
        name = full_ec_model
        library = egs_genvelope
        base geometry = air_box
        inscribed geometries = source_holder Y90_disc_source electrode
    :stop geometry:

    simulation geometry = full_ec_model

:stop geometry definition:

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

:start source definition:
	:start source:
        	name = y_90_base_source
        	library = egs_isotropic_source
			geometry = Y90_disc_source
        	region selection = IncludeSelected
        	selected regions = 5 
        	
        	:start shape:
            		type = box
            		box size = 0.6 0.6 0.03 
            		origin = -0.3 -0.3 -0.05125 
        	:stop shape:
       
        
        	:start spectrum:
            		type = monoenergetic
            		energy = 1
        	:stop spectrum:
    	:stop source:

    	:start source:
        	name = Y90_source
        	library = egs_radionuclide_source
        	base source = y_90_base_source
		
			charge = -1
        
        	:start spectrum:
            		type = radionuclide
            		nuclide = Y-90
        	:stop spectrum:
    	:stop source:

    	simulation source = Y90_source
:stop source definition:

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

:start run control: 
	ncase = 1000
:stop run control: 

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


:start scoring options:

    :start calculation geometry:
        geometry name = full_ec_model
        cavity regions = 25
        cavity mass = 0.3
        cavity geometry = full_ec_model
    :stop calculation geometry:

:stop scoring options:

EC_Y90_materials_2_pegs4dat.txt
EC_LV_Y90_geometry.txt

[rtownson]

See the box shape inputs:
https://nrc-cnrc.github.io/EGSnrc/doc/pirs898/classEGS__BoxShape.html

There's no origin parameter. To transform it you use a transformation input block.

Also, I don't actually see the global region number 5 in egs_view, are you sure that's your source region? selected regions takes a global region number, or a label.