45keV_XCOM_HighRes.egsinp 17 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368
  1. TITLE= Energy deposition kernel for 045 keV photons
  2. ##########################
  3. :start I/O control:
  4. IRESTART= first # first (0) First run for this data set
  5. # restart (1) Restart a previous run
  6. # analyze (3) Just read in the raw data and do the statistical analysis
  7. # parallel (5) Combine results from previous parallel runs
  8. STORE DATA ARRAYS= yes # yes: (0) Store data arrays for re-use
  9. # no : (1) don't store them
  10. PRINT OUT EDK FILE= yes # produce energy deposition file
  11. # in old kernels' format ?
  12. # yes: (0) EDK stored in old format files
  13. # no : (1) don't produce EDK files in old format
  14. :stop I/O control:
  15. ##########################
  16. ##########################
  17. :start Monte Carlo inputs:
  18. NUMBER OF HISTORIES= 999999999 #
  19. INITIAL RANDOM NO. SEEDS= 1, 5858
  20. IFULL= ENERGY DEPOSITION KERNEL # cavity calculation (0) calculate dose in cavity regions
  21. # energy deposition kernels (1)
  22. # dose calculation (2)
  23. # dose and edk (3)
  24. DOPPLER BROADENING= On # On or Off
  25. # On: Default EGSnrc implementation
  26. # Off: Neglects Doppler broadening
  27. :stop Monte Carlo inputs:
  28. #########################
  29. ##########################
  30. :start geometrical inputs:
  31. NUMBER OF CONES= 48 # number of cones (individual or by group)
  32. # If omitted or ZERO, pure spherical geometry
  33. # assumed.
  34. ANGLES= 3.75 # ANGLES defining the geometry (reals)
  35. # No needed in pure spherical geometries.
  36. #
  37. # For group input there must be as many entries
  38. # as for the NUMBER OF CONES, i.e. :
  39. # NCON1,NCON2,...,NCONn
  40. # DANG1,DANG2,...,DANGn
  41. #
  42. # For individual input, ncones must be equal
  43. # to the number of entries, i.e.:
  44. # ncones
  45. # DANG1, DANG2,...,DANGncones
  46. # NUMBER OF SPHERES = 24 # number of spheres (individual or by group)
  47. # For individual inputs as below, number of spheres can be omitted
  48. RADII= 0.01,0.03,0.05,0.08,0.11,0.15,0.2,0.3,0.4,0.6,0.8,1.1,
  49. 1.5,2.0,3.0,4.0,6.0,8.0,10.0,15.0,20.0,30.0,45.0,60.0
  50. # RADII= 60.
  51. # Below is radii used by Mackie and Rogers which was scaled above to include lower radii
  52. # RADII= 0.05,0.1,0.15,0.2,0.3,0.4,0.5,0.6,0.8,1.0, # radii of spheres defining the geometry (reals)
  53. # 1.5,2.0,3.0,4.0,5.0,6.0,8.0,10.0,15.0,20.0, #
  54. # 30.0,40.0,50.0,60.0 # For group input there must be as many entries
  55. # as for the NUMBER OF SPHERES, i.e. :
  56. # NSPH1,NSPH2,...,NSPHn
  57. # DRAD1,DRAD2,...,DRADn
  58. MEDIA= WaterLiq; #Media in the problem:
  59. #These must match exactly, including case, one
  60. #of the media names in the pegs4 data set being
  61. #used in the problem.
  62. #The maximum length of name is 24 characters
  63. #They are automatically left justified on input.
  64. MEDNUM= 1 # define what medium goes where
  65. # use region numbers to define this
  66. # (region numbers start at 2 and increase
  67. # number from smallest angle to largest angle
  68. # of the conical intervals and innermost radius to
  69. # outermost radius)
  70. #Next we specify which media are in
  71. #which geometric regions
  72. #note that by default all regions contain
  73. #medium 1 and which medium to input as 1 should
  74. #be selected with this in mind.
  75. START REGION= 2 #This puts water everywhere
  76. STOP REGION= 1153
  77. :stop geometrical inputs:
  78. #########################
  79. ##########################
  80. :start source inputs:
  81. INCIDENT PARTICLE= photon # electron,photon,positron
  82. INCIDENT ENERGY= monoenergetic # monoenergetic, spectrum;
  83. INCIDENT KINETIC ENERGY(MEV)= 0.045 # only use for "monoenergetic"
  84. #If INCIDENT ENERGY= spectrum:
  85. # SPEC FILENAME= full name of file containing energy spectrum
  86. # SPEC IOUTSP= include # none,include;
  87. # none: no spectrum data in .egslst file
  88. # include: output spectrum data to .egslst file
  89. SOURCE NUMBER= 0 # 0,1,2
  90. # 0: point source AT origin, emission along Z-axis
  91. # 1: point source AT origin, isotropically radiating in 4Pi
  92. # 2: point source NEAR origin, emission along Z-axis
  93. # ZIN= 0.000001 # only for source number 2
  94. # source offset on Z-axis
  95. # Option used to emulate old way of calculating EDK
  96. :stop source inputs:
  97. #########################
  98. ##########################
  99. :start MC transport parameter:
  100. Global ECUT= 0.512 # Electron cutoff for transport
  101. Global PCUT= 0.0010 # Photon cutoff for transport
  102. Global SMAX= 0.0 # Maximum step size in cm (not needed
  103. # unless old PRESTA algorithm used)
  104. ESTEPE= 0.25 # Max fractional continuous energy loss
  105. # per step. Use 0.25 unless using
  106. # PRESTA-I
  107. XImax= 0.5 # Max first elastic scattering moment
  108. # per step. Using default.
  109. Skin depth for BCA= 3 # Distance from a boundary (in elastic
  110. # MFP) at which the algorithm will go
  111. # into single scattering mode (using
  112. # default here)
  113. Boundary crossing algorithm= EXACT # exact,PRESTA-I;
  114. # exact: cross boundaries in single scattering
  115. # mode (distance at which to go into
  116. # single scattering mode determined by
  117. # "Skin depth for BCA"
  118. # PRESTA-I: cross boundaries with lateral
  119. # correlations off and force multiple
  120. # scattering mode
  121. Electron-step algorithm= PRESTA-II # PRESTA-II,PRESTA-I;
  122. # Determines the algorithm used to take
  123. # into account lateral and longitudinal
  124. # correlations in a condensed history
  125. # step
  126. Spin effects= on # Off (default),On;
  127. # Turns off/on spin effects for electron
  128. # elastic scattering. Spin On is
  129. # ABSOLUTELY necessary for good
  130. # backscattering calculations. Will
  131. # make a difference even in `well
  132. # conditioned' situations (e.g. depth
  133. # dose curves).
  134. Brems angular sampling= KM # Simple,KM (default);
  135. # Simple: leading term of Koch-Motz
  136. # dist'n used to determine angle
  137. # of bremsstrahlung photons
  138. # KM: Koch-Motz distribution used to
  139. # determine angle
  140. Triplet production= Off # On or Off (default).
  141. # Turns on/off simulation of triplet production.
  142. # On: Borsellino's first Born approximation is
  143. # used to sample triplet events based on the
  144. # triplet cross-section data.
  145. Brems cross sections= nrc # BH (default),NIST;
  146. # BH: Bethe-Heitler cross-sections used
  147. # NIST: NIST cross-sections used
  148. Bound Compton scattering= On # Off (default),On, simple or norej;
  149. # Off: Klein-Nishina used for Compton
  150. # scattering
  151. # On: Impulse approximation used for
  152. # Compton scattering
  153. # simple: impulse approximation incoherent
  154. # scattering function used (i.e., no
  155. # Doppler broadenning).
  156. # norej: the actual total bound Compton cross
  157. # section is used and there are no
  158. # rejections at run time.
  159. Radiative Compton corrections= On # On or Off (default).
  160. # On: include radiative corrections for Compton
  161. # scattering.
  162. Electron Impact Ionization= On # Off (default), On, ik, casnati, kolbenstvedt,
  163. # gryzinski, penelope.
  164. # On or ik: use Kawrakow's theory to derive
  165. # EII cross-sections.
  166. # casnati: use the cross-sections of Casnati from
  167. # $HEN_HOUSE/data/eii_casnati.data.
  168. # Similarly for kolbenstvedt, gryzinski and
  169. # penelope. Case-sensitive except for Off, On or
  170. # ik options.
  171. Pair angular sampling= Simple # Off, Simple (default),KM);
  172. # Simple: use leading term of K-M
  173. # dist'n
  174. # KM: use complete Koch and Motz dist'n
  175. # Off: angle of pairs is m/E--like old EGS4
  176. Photoelectron angular sampling= On # Off (default),On;
  177. # Off: Photoelectrons get direction of
  178. # photon that creates them
  179. # On: Sauter's formula is used
  180. Pair cross sections= NRC # BH (default) or NRC.
  181. # BH: use Bethe-Heitler pair production
  182. # cross-sections.
  183. # NRC: use NRC pair production cross-sections
  184. # (in file $HEN_HOUSE/data/pair_nrc1.data).
  185. Photon cross sections= xcom # si (Storm-Israel--the default),
  186. # epdl (Evaluated Photon Data Library),
  187. # xcom (NIST Photon Cross Sections Database)
  188. # pegs4 (PEGS4 file photon data)
  189. # User can supply their own cross-section
  190. # data as well.
  191. # Hence this entry is case-sensitive.
  192. Photon cross-sections output= Off # Off (default) or On.
  193. #On: file $EGS_HOME/user_code/inputfile.xsections
  194. # is created with the photon cross-section
  195. # data used.
  196. Compton cross sections= compton_sigma.data # Bound Compton cross-section data.
  197. # User-supplied bound Compton
  198. # cross-sections in the file
  199. # comp_xsections_compton.data in
  200. # directory $HEN_HOUSE/data/, where
  201. # comp_xsections is the name supplied
  202. # for this input. Uses compton_sigma.data
  203. # by default.
  204. Rayleigh scattering= On # Off (default),On, custom;
  205. # Off: no coherent scattering
  206. # On: simulates coherent scattering
  207. # custom: user must provide media names
  208. # for wich form factor (FF) files
  209. # will be provided. For the rest
  210. # of the media, default atomic FF
  211. # used.
  212. #
  213. # IF 'custom' Rayleigh option then:
  214. #
  215. #ff media names = A list of media names (must match media found in
  216. # PEGS4 data file) for which the user is going to
  217. # provide custom Rayleigh form factor data.
  218. #ff file names = A list of names of files containing the Rayleigh
  219. # form factor data for the media specified by
  220. # the ff media names = input above. Full directory
  221. # paths must be given for all files, and for each medium
  222. # specified, iray_ff_media(i), there must be a
  223. # corresponding file name, iray_ff_file(i). For
  224. # example files, see the directory
  225. # $HEN_HOUSE/data/molecular_form_factors.
  226. Atomic relaxations= On # Off (default),On;
  227. # On: use correct cross section
  228. # for p.e. events and shell vacancies
  229. # for Compton & p.e. events are relaxed
  230. # via emission of fluorescent X-Rays,
  231. # Auger and Koster-Cronig electrons
  232. # electrons
  233. # Atomic relaxations, Rayleigh scattering, Photoelectron angular sampling and
  234. # Bound compton scattering can be turned on/off on a region by region basis.
  235. # Instead of simply "On" or "Off" for these cases put:
  236. # Atomic relaxations= On (or Off) in Regions
  237. # Relaxations start region= 1, 40 #turns relaxations on in regions 1-10 and
  238. # Relaxations stop region= 10, 99 #40-99
  239. #
  240. # Rayleigh scattering= On (or Off) in Regions
  241. # Rayleigh start region= 1, 40
  242. # Rayleigh stop region= 10, 99
  243. #
  244. # Photoelectron angular sampling= On (or Off) in Regions
  245. # PE sampling start region= 1, 40
  246. # PE sampling stop region= 10, 99
  247. #
  248. # Bound Compton scattering= On (or Off) in Regions
  249. # Bound Compton start region= 1, 40
  250. # Bound Compton stop region= 10, 99
  251. #ECUT, PCUT and SMAX can also be set on a region-by-region basis.
  252. #
  253. #Set XXXX= f_value1, f_value2, ...
  254. #Set XXXX start region= i_value1, i_value2, ...
  255. #Set XXXX stop region= j_value1, j_value2, ...
  256. #
  257. #where XXXX is ECUT, PCUT or SMAX ,
  258. #f_value1, f_value2,... are the desired values for XXXX
  259. #and i_value_i and j_value_i are the start and
  260. #stop regions.
  261. Set PCUT= 0
  262. Set PCUT start region= 1
  263. Set PCUT stop region= 1
  264. Set ECUT= 0
  265. Set ECUT start region= 1
  266. Set ECUT stop region= 1
  267. Set SMAX= 0
  268. Set SMAX start region= 1
  269. Set SMAX stop region= 1
  270. :stop MC transport parameter:
  271. #########################
  272. ##########################
  273. :start variance reduction:
  274. ELECTRON RANGE REJECTION= on #Off,On;
  275. #On: if charged particle energy is below ESAVEIN
  276. # and it cannot get out of current region
  277. # with energy > ECUT, the particle is
  278. # terminated
  279. ESAVEIN= 0.005 #Energy below which range rejection is
  280. #considered
  281. EXPONENTIAL TRANSFORM C= 0.0000 # parameter for pathlength biasing
  282. # <0 for path length shortening
  283. # (useful for surface problems)
  284. # >0 for path length increase
  285. # (useful for shielding problems)
  286. # if 0.0, no biasing done
  287. PHOTON FORCING= On #Off (default),On;
  288. #On: force photons to interact in geometry
  289. START FORCING= 1 #Start forcing at this interaction number
  290. STOP FORCING AFTER= 1 #Number of photon interactions after which
  291. #to stop forcing photon interactions
  292. :stop variance reduction:
  293. #########################
  294. #########################
  295. :start plot control:
  296. PLOTTING= histogram #Off: do not create plot files
  297. #histogram: create histogram plots
  298. #point: create xy plots
  299. PLOT RADIAL REGION IX= 1,12,24 #Indices of spheres for which to plot depth-
  300. #dose data (0 for no depth-dose plots)
  301. PLOT CONICAL REGION IC= 1,24,48 #Indices of cones for which to plot dose vs
  302. #radius data (0 for no dose vs radius plots)
  303. :stop plot control:
  304. ########################