[ T i t l e ] Input file for lecture about X-ray therapy [ P a r a m e t e r s ] icntl = 0 # (D=0) 3:ECH 5:NOR 6:SRC 7,8:GSH 11:DSH 12:DUMP maxcas = 1000000 # (D=10) number of particles per one batch maxbch = 100 # (D=10) number of batches emin(2) = 1.000000000E-10 # (D=1.0) cut-off energy of neutron (MeV) dmax(2) = 20.0000000 # (D=emin(2)) data max. energy of neutron (MeV) emin(12) = 1.000000000E-01 # (D=1.d9) cut-off energy of electron (MeV) emin(13) = 1.000000000E-01 # (D=1.d9) cut-off energy of positron (MeV) emin(14) = 1.000000000E-03 # (D=1.d9) cut-off energy of photon (MeV) dmax(12) = 1000.00000 # (D=emin(12)) data max. energy of electron (MeV) dmax(13) = 1000.00000 # (D=emin(13)) data max. energy of positron (MeV) dmax(14) = 1000.00000 # (D=emin(14)) data max. energy of photon (MeV) igamma = 2 # (D=0) 0:No, 1:Old, 2:EBITEM, 3:EBITEM+Isomer ipnint = 1 # (D=0) 0: no, 1: consider photo-nuclear reaction negs = 1 # (D=0) =1 EGS photon and electron # negs = -1 # (D=0) =1 EGS photon and electron file(6) = phits.out # (D=phits.out) general output file name file(7) = c:/phits/data/xsdir.jnd # (D=c:/phits/data/xsdir.jnd) nuclear data input file name file(20) = c:/phits/XS/egs # (D=c:/phits/XS/egs) EGS library data folder name [ S o u r c e ] s-type = 13 # mono-energetic axial source proj = electron # kind of incident nucleus e0 = 10.000 # energy of beam [MeV/u] r0 = 1.5000 # radius [cm] r1 = 1.0000 # radius [cm] for FWHM x0 = 0.0000 # (D=0.0) center position of x-axis [cm] y0 = 0.0000 # (D=0.0) center position of y-axis [cm] z0 = -1.00 # minimum position of z-axis [cm] z1 = -1.00 # maximum position of z-axis [cm] dir = 1.0000 # z-direction of beam [cosine] # totfact = 11.6E10 # (D=1.0) 250‚ŤA, 12micro second totfact = 0.25*12*1e-6/1.6e-19 [ M a t e r i a l ] MAT[ 1 ] # Water 1H 2.0000000E+00 16O 1.0000000E+00 MAT[ 2 ] # Air (1.20e-3 g/cm3) 14N 8.0000000E+00 16O 2.0000000E+00 MAT[ 3 ] # Tungsten target 180W 1.2000000E-01 182W 2.6498000E+01 183W 1.4314000E+01 184W 3.0642000E+01 186W 2.8426000E+01 MAT[ 4 ] # Copper target 63Cu 6.9170000E+01 65Cu 3.0830000E+01 MAT[ 5 ] # Collimator (17g/cm3, from AAPM report 72, TG-50) 184W -9.0500000E+01 58Ni -6.5000000E+00 56Fe -3.0000000E+00 MAT[ 6 ] # Flattening filter 63Cu 6.9170000E+01 65Cu 3.0830000E+01 [ Mat Name Color ] mat name size color # HSB color value 0 void 1.00 lightgray # { -0.800 1.000 1.000 } 1 Water 1.00 pastelblue # { 0.250 0.400 1.000 } 2 Air 1.00 pastelcyan # { 0.400 0.400 1.000 } 3 TungstenTarget 1.00 matblack # { -0.200 1.000 1.000 } 4 CopperTarget 1.00 orangeyellow # { 0.867 1.000 1.000 } 5 Collimator 1.00 magenta # { 0.067 1.000 1.000 } 6 FlatteningFilter 1.00 darkgreen # { 0.600 1.000 0.500 } [ C e l l ] 1 3 -19.25 1 -2 -4 $ W for X-ray generator 2 4 -8.94 2 -3 -4 $ Cu for X-ray generator 3 1 -1.0 -5 $ Water phantom 98 2 -1.20e-3 #1 #2 #3 -999 #13 #14 #15 #16 #17 #18 $ Air 99 -1 999 $ Outer region 13 5 -17.0 21 22 -23 -24 $ Primary collimator 14 5 -17.0 -31 $ Movable collimator (X-1) 15 5 -17.0 -32 $ Movable collimator (X-2) 16 5 -17.0 -33 $ Movable collimator (Y-1) 17 5 -17.0 -34 $ Movable collimator (Y-2) 18 6 -8.94 -35 $ Flattening filter [ S u r f a c e ] 1 pz 0.0 2 pz 0.15 3 pz 0.15+0.30 4 cz 1.5 5 rpp -15. 15. -15. 15. 30. 60. 999 so 100.0 21 trc 0.0 0.0 (0.45+6.0) 0.0 0.0 -6.0 2.0 0.5 22 pz 0.45 23 pz 0.45+6.0 24 cz 10.0 set:c1[2.5] 31 rpp (-10-c1) -c1 -5.0 5.0 11 16 32 rpp c1 (10+c1) -5.0 5.0 11 16 33 rpp -5.0 5.0 (-10-c1) -c1 17 22 34 rpp -5.0 5.0 c1 (10+c1) 17 22 35 trc 0.0 0.0 10 0.0 0.0 -1.5 4.0 1.0e-10 [ T - Deposit ] title = Dose in xyz mesh mesh = xyz # mesh type is xyz scoring mesh x-type = 2 xmin = -0.0000000 xmax = 0.2000000 nx = 600 y-type = 2 ymin = -0.0000000 ymax = 0.05000000 ny = 1 z-type = 2 zmin = -0.010000 zmax = 0.31000 nz = 960 unit = 5 # unit is [J/m^3/source] material = all # (D=all) number of specific material output = dose # total deposit energy 2D-type = 3 # 1:Cont, 2:Clust, 3:Color, 4:xyz, 5:mat, 6:Clust+Cont, 7:Col+Cont axis = xz # axis of output file = deposit_xz_Joule.out # file name of output for the above axis part = all gshow = 1 # 0: no 1:bnd, 2:bnd+mat, 3:bnd+reg 4:bnd+lat epsout = 1 [ T - Deposit ] title = Energy deposition in reg mesh mesh = reg # mesh type is xyz scoring mesh reg = 1 2 unit = 5 # unit is [J/m^3/source] material = all # (D=all) number of specific material output = dose # total deposit energy axis = reg # axis of output file = deposit.out # file name of output for the above axis part = all [ T - T r a c k ] title = Track in xyz mesh mesh = xyz # mesh type is xyz scoring mesh x-type = 2 # x-mesh is linear given by xmin, xmax and nx xmin = -25.00000 # minimum value of x-mesh points xmax = 25.00000 # maximum value of x-mesh points nx = 50 # number of x-mesh points y-type = 2 # y-mesh is linear given by ymin, ymax and ny ymin = -5.000000 # minimum value of y-mesh points ymax = 5.000000 # maximum value of y-mesh points ny = 1 # number of y-mesh points z-type = 2 # z-mesh is linear given by zmin, zmax and nz zmin = -30.0000 # minimum value of z-mesh points zmax = 60.00000 # maximum value of z-mesh points nz = 140 # number of z-mesh points e-type = 1 # e-mesh is given by the below data ne = 1 # number of e-mesh points 0.0 1000.0 unit = 1 # unit is [1/cm^2/source] 2D-type = 3 # 1:Cont, 2:Clust, 3:Color, 4:xyz, 5:mat, 6:Clust+Cont, 7:Col+Cont axis = xz # axis of output file = track_xz.out # file name of output for the above axis part = electron photon neutron gshow = 1 # 0: no 1:bnd, 2:bnd+mat, 3:bnd+reg 4:bnd+lat epsout = 1 # (D=0) generate eps file by ANGEL [ T - T r a c k ] title = Track in xyz mesh mesh = xyz # mesh type is xyz scoring mesh x-type = 2 # y-mesh is linear given by ymin, ymax and ny xmin = -5.000000 # minimum value of y-mesh points xmax = 5.000000 # maximum value of y-mesh points nx = 1 # number of y-mesh points y-type = 2 # x-mesh is linear given by xmin, xmax and nx ymin = -25.00000 # minimum value of x-mesh points ymax = 25.00000 # maximum value of x-mesh points ny = 50 # number of x-mesh points z-type = 2 # z-mesh is linear given by zmin, zmax and nz zmin = -30.0000 # minimum value of z-mesh points zmax = 60.00000 # maximum value of z-mesh points nz = 140 # number of z-mesh points e-type = 1 # e-mesh is given by the below data ne = 1 # number of e-mesh points 0.0 1000.0 unit = 1 # unit is [1/cm^2/source] 2D-type = 3 # 1:Cont, 2:Clust, 3:Color, 4:xyz, 5:mat, 6:Clust+Cont, 7:Col+Cont axis = yz # axis of output file = track_yz.out # file name of output for the above axis part = electron photon neutron gshow = 1 # 0: no 1:bnd, 2:bnd+mat, 3:bnd+reg 4:bnd+lat epsout = 1 # (D=0) generate eps file by ANGEL [ T - Deposit ] title = Energy deposition in reg mesh mesh = reg # mesh type is xyz scoring mesh reg = 3 unit = 0 # unit is [Gy/source] material = all # (D=all) number of specific material output = dose # total deposit energy axis = reg # axis of output file = deposit_org.out # file name of output for the above axis part = all [ T - C r o s s ] title = x-distribution of fluence mesh = xyz # mesh type is xyz scoring mesh x-type = 2 # x-mesh is linear given by xmin, xmax and nx xmin = -25.00000 # minimum value of x-mesh points xmax = 25.00000 # maximum value of x-mesh points nx = 50 # number of x-mesh points y-type = 2 # y-mesh is linear given by ymin, ymax and ny ymin = -5.000000 # minimum value of y-mesh points ymax = 5.000000 # maximum value of y-mesh points ny = 1 # number of y-mesh points z-type = 2 # z-mesh is linear given by zmin, zmax and nz zmin = -30.0000 # minimum value of z-mesh points zmax = 60.00000 # maximum value of z-mesh points nz = 9 # number of z-mesh points e-type = 1 # e-mesh is given by the below data ne = 1 # number of e-mesh points 0.0 50.0 unit = 1 # unit is [1/cm^2/source] axis = x # axis of output file = cross_x.out # file name of output for the above axis output = flux # surface crossing flux part = electron photon neutron angel = ymin[1.E05] ymax[1.E13] epsout = 1 # (D=0) generate eps file by ANGEL [ T - C r o s s ] title = energy spectrum mesh = xyz # mesh type is xyz scoring mesh x-type = 2 # x-mesh is linear given by xmin, xmax and nx xmin = -25.00000 # minimum value of x-mesh points xmax = 25.00000 # maximum value of x-mesh points nx = 1 # number of x-mesh points y-type = 2 # y-mesh is linear given by ymin, ymax and ny ymin = -5.000000 # minimum value of y-mesh points ymax = 5.000000 # maximum value of y-mesh points ny = 1 # number of y-mesh points z-type = 2 # z-mesh is linear given by zmin, zmax and nz zmin = -30.0000 # minimum value of z-mesh points zmax = 60.00000 # maximum value of z-mesh points nz = 9 # number of z-mesh points e-type = 3 # e-mesh is log given by emin, emax and ne emin = 0.01000 # minimum value of e-mesh points emax = 50.00000 # maximum value of e-mesh points ne = 50 # number of e-mesh points unit = 2 # unit is [1/cm^2/MeV/source] axis = eng # axis of output file = cross_eng.out # file name of output for the above axis output = flux # surface crossing flux part = electron photon neutron angel = ymin[1.E05] ymax[1.E13] epsout = 1 # (D=0) generate eps file by ANGEL [ T - Deposit ] title = z-distribution of dose mesh = xyz # mesh type is xyz scoring mesh x-type = 2 # x-mesh is linear given by xmin, xmax and nx xmin = -0.500000 # minimum value of x-mesh points xmax = 0.500000 # maximum value of x-mesh points nx = 1 # number of x-mesh points y-type = 2 # y-mesh is linear given by ymin, ymax and ny ymin = -0.500000 # minimum value of y-mesh points ymax = 0.500000 # maximum value of y-mesh points ny = 1 # number of y-mesh points z-type = 2 # z-mesh is linear given by zmin, zmax and nz zmin = 20.00000 # minimum value of z-mesh points zmax = 60.00000 # maximum value of z-mesh points nz = 120 # number of z-mesh points unit = 0 # unit is [Gy/source] material = all # (D=all) number of specific material output = dose # total deposit energy axis = z # axis of output file = deposit_z.out # file name of output for the above axis part = all angel = ymin[1.E-02] ymax[1.E-0] epsout = 1 # (D=0) generate eps file by ANGEL [END]