Pencil beam in air

A simple simulation of a single pencil beam in air. The setup and the method presented in this example can be used for validation of lateral propagation of a single pencil beam in air. It is useful for validation and optimisation of beam envelope propagation. In this example an emittance model is used and the description can be found in Emittance model.

Input file

The beam is defined by means of emittance model with parameters corresponding to the beam sizes of 0.5, 0.3, 0.2, 0.4 and 0.8 cm (sigma) at distances from the isocentre of -20, -10, 0, 10 and 20 cm, respectively. The distance of the source to the isocentre is 40 cm. The dose is scored in an air phantom of size 5x5x40 cm³ with pixel size 1x1x10 cm³ . The material of the phantom inherits from air and changes the radiation length to large value to minimise the effect of Multiple Coulomb Scattering (MCS). This is required when the measurements of the lateral beam sizes are performed in air.

### Phantom ###
region<
    ID=Phantom
    O=[0, 0, -20]
    f=[ 0, 0, 1 ]
    u=[ 0, 1, 0 ]
    L=[ 5, 5, 40]
    voxels=[50,50,40]
    material=airNOMCS # material with no MCS effect
region>

### Beam ###
field: 1 ; O = [0,0,-40]; L=[10,10,20]; pivot = [0.5,0.5,0.2]
pb<
    ID = 1
    fieldID = 1
    particle = proton
    T = 150 # MeV
    Xsec = emittance

    # the Twiss parameters corresponding to the beam sizes of 0.5, 0.3, 0.2, 0.4 and 0.8 cm (sigma)
    # at distances from the ISO of -20, -10, 0, 10 and 20 cm, respectively
    twissAlphaX=-1.235900963624794
    twissBetaX=8.3916636857885
    emittanceX=0.0034387868648755714

    emittanceRefPlaneDistance=40; # Field_1 origin is at 40 cm from ISO
pb>

nprim=1e6

### Materials ###
material<
    ID=airNOMCS
    longName ='air without MCS'
    basedOn = air
    rho = 1.205e-3
    Lrad = 1e32 # high radiation length is set to minimise the impact of MCS
material>

Simple analysis

$ mhd_viewer.py out/Dose.mhd

Dose deposition in air.