Defining a phasespace scoring volume rather than a plane

[matthew-efseaff]

Is it possible to define a geometry (e.g. a box) and score particles on entry? I'm familiar with:

          library = egs_planes
          type = EGS_Zplanes 

My .egsinp code:

:start geometry definition:

    :start geometry:
        library = egs_box
        box size = 10 10 0.007
        name = box_at_origin
    :stop geometry:

    :start geometry:
        library = egs_gtransformed
        my geometry = box_at_origin
        :start transformation:
            translation = 0 0 0  # for future flexibility
        :stop transformation:
        name = Detector2
    :stop geometry:

  simulation geometry = EGS_Mesh

:stop geometry definition:

:start ausgab object definition:

    :start ausgab object:
        library                  = egs_phsp_scoring
        name                     = Box_2
        output format            = EGSnrc
        particle type            = all  #, photons, or charged
        score multiple crossers  = yes  #no (default) or yes (EGSnrc format only)
        phase space geometry     = Detector2
        score particles on       = entry  #entry, exit, entry and exit [default]
    :stop ausgab object:

:stop ausgab object definition:

However, this fails on initialization (no error message, just stops). Thank you.

[rtownson]

This will crash because there is no material defined for the geometry. Add

        :start media input:
            media = air
        :stop media input:

to the box.

Also, for it to score on 'entry' the beam needs to start outside the geometry. You can do this by placing everything inside a bigger box. E.g.

:start geometry definition:

    :start geometry:
        library = egs_box
        box size = 100
        name = big_box
        :start media input:
            media = air
        :stop media input:
    :stop geometry:

    :start geometry:
        library = egs_box
        box size = 10 10 0.007
        name = box_at_origin
        :start media input:
            media = air
        :stop media input:
    :stop geometry:

    :start geometry:
        library = egs_gtransformed
        my geometry = box_at_origin
        :start transformation:
            translation = 0 0 0  # for future flexibility
        :stop transformation:
        name = Detector2
    :stop geometry:

    :start geometry:
        library = egs_genvelope
        name = world
        base geometry = big_box
        inscribed geometries = Detector2
     :stop geometry:

  simulation geometry = world

:stop geometry definition: