Lecture Notes of Dose distribution calculation inside phantom with ...

More documents

Recommendations

Info

! Unit ifto Output file.!! ----------! INPUT FILE! ----------! Record 1 title (80A1) Title line.! --------! Record 2 nmed (I10) Number of media in problem.! --------! Record 3 media(j,i) (24A1) Media names (j=1,24, I=1,nmed lines).! -------- Note that entire volume is initially! set to medium.! ---------! Record 4 chard characteristic distance for each medium! ---------! --------! Record 5 maxx, maxy, maxz Number of voxels in the X,Y,Z directions! -------- If 0! input, one per line, the maxx + 1 x boundaries!! if maxx < 0! input smallest x boundary, followed by abs(maxx) pairs! one per line: voxel width, # voxels with this width!! for example: starting at record 5! -1,-1,-1! 0.0! 1.0,16! 0.0! 1.0,16! 0.0! 1.0,16! defines a 16x16x16 cube of 1cm**3 voxels with a total of 4097 reg! or! -1,-1,3! 0.0! 1.0,16! 0.0! 1.0,16! 0.0! 5.0! 10.0! defines a 16x16x10 cube with 1x1x5 cm voxels stacked 2 deep! --------! Record 7 ybound! --------! Record 8 zbound! --------! ---------! Record 9 ecutin,pcutin, For medium 1 at first! ---------ipeangsw, Switches for PE-angle sampling,! iedgeflsw, K & L-edge fluorescence,! iaugersw, K & L-Auger! iraysw, Rayleigh scattering,! ipolarsw, Linearly-polarized photon scattering,! incohrsw, S/Z rejection,! iprofrsw, Doppler broadening,! impacrsw electron impact ionization (0=off, 1=on).! --------! Record 10 il,iu, jl,ju, kl,ku, medtmp, rhotmp,ecutin,pcutin! -------- Line is repeated until a blank line found! All regions default to medium 1 with its! default density unless changed here.! For all voxels with! IL
! Record 10a ipeangsw, Switches for PE-angle sampling,! --------- iedgeflsw, K & L-edge fluorescence,! iaugersw, K & L-Auger! iraysw, Rayleigh scattering,! ipolarsw, Linearly-polarized photon scattering,! incohrsw, S/Z rejection,! iprofrsw, Doppler broadening,! impacrsw electron impact ionization (0=off, 1=on).!! --------! Record 11 il,iu, jl,ju, kl,ku,izscan! --------- Regions for which the dose will be output.! IZSCAN non-zero to get z-scan per page,! otherwise output is an x-scan per page.!! ---------! Record 12 iwatch Switch for tracking events with swatch:! --------- (0=No, 1=each interaction,! 2=each step)! ---------! Record 13 ibrdst,iprdst, Switches for bremsstrahlung and pair! --------- ibrspl,nbrspl production ANGLE SAMPLING, and brems-! strahlung SPLITTING:! ibrdst=0 No (use default: theta=m/E)! 1 Yes (recommended)! iprdst=0 No (use default: theta=m/E)! 1 Yes (low-order distribution)! 2 Yes (recommended)! ibrspl=0 No! 1 Yes (NBRSPL=splitting factor)! ----------------------------------------------------------------------subroutine getvoxel(ifto)implicit noneinclude 'include/egs5_h.f'include 'include/egs5_bounds.f'include 'include/egs5_brempr.f'include 'include/egs5_edge.f'include 'include/egs5_eiicom.f'include 'include/egs5_elecin.f'include 'include/egs5_media.f'include 'include/egs5_misc.f'include 'include/egs5_thresh.f'include 'include/egs5_uphiot.f'include 'include/egs5_useful.f'include 'include/egs5_userpr.f'include 'include/egs5_usersc.f'include 'include/egs5_uservr.f'include 'include/egs5_userxt.f'include 'pegscommons/mscom.f'include 'auxcommons/aux_h.f'include 'auxcommons/edata.f'include 'auxcommons/geoxyz.f'include 'auxcommons/instuf.f'include 'auxcommons/voxel.f'include 'auxcommons/watch.f'include 'include/randomm.f'integer iftoreal*8* ecutin,ecutmn,ek0,pcutin,rhotmp,totphi,* thetax,thetay,thetaz,xlower,* xupper,ylower,yupper,width! Main EGS "header" file! COMMONs required by EGS5 code! PEGS common! Auxiliary-code "header" file! Auxiliary-code COMMONs! Additional (non-EGS5) COMMON! Argumentinteger i,igroup,ii,iiz,il,in,irl,iu,ixinu,* ixx,izn,j,jl,ju,jxx,jyinu,k,kl,ku,maxbd,maxx,maxy,* maxz,medtmp,moreOutput,n,ner,ngroup,nn,nnn! Local variablesinteger ipeangsw,iedgeflsw,iaugersw,iraysw,ipolarsw,incohrsw,* iprofrsw,impacrsw59
Page 1 and 2:
KEK Internal 2011-8December 2011R/D
Page 3 and 4:
Lecture Notes ofDose distribution c
Page 5 and 6:
2.8.1 Statistical uncertainty . . .
Page 7 and 8:
1 ucxyz phantom.f ucxyz phantom.f
Page 9 and 10:
common/score/! Variables to score*
Page 11 and 12:
emaxe = 0.D0 ! dummy value to extra
Page 13 and 14: vin=dsin(phai)*sinthdis=sposi/winxp
Page 15 and 16: 2.9 Step 9: Output of results ()wr
Page 17 and 18: ipeangswiedgeswiaugerswirayswipolar
Page 19 and 20: ! ---------------------------------
Page 21 and 22: end ifend ifend ifdnear(np)=dnearlr
Page 23 and 24: 4 ucxyz phantom.f (egs5job.out, e
Page 25 and 26: 4.2 1. cp ucxyz phantom1.f ucxyz ph
Page 27 and 28: call randomset(rnnow)do ie=1,nsebin
Page 29 and 30: Linux 4 ucxyz phantom3 25 Enter
Page 31 and 32: Linux ucxyz phantom5 , 4 25 "Doe
Page 33 and 34: 1 Outlines of sample user code ucxy
Page 35 and 36: common/score/! Variables to score*
Page 37 and 38: Output the material data and parame
Page 39 and 40: totke = totke + availkelatchi=0! Ke
Page 41 and 42: * rhor(irl),dose(i,j,kkk),doseun(i,
Page 43 and 44: Ifahistorynumberislessthanmaxpict,
Page 45 and 46: 3 Exercise problems3.1 Problem 1 :
Page 47 and 48: ! Source position from phantom surf
Page 49 and 50: At the above data, a first 201 is t
Page 51 and 52: &INP AE=0.521,AP=0.0100,UE=2.011,UP
Page 53 and 54: (b) Add new material (Iron).ChangeA
Page 55 and 56: !**********************************
Page 57 and 58: !----------------------------------
Page 59 and 60: faexp=0.D0faexps=0.D0faexp2s=0.D0fe
Page 61 and 62: ! =======================if (iwatch
Page 63: do k=kdosl(idose),kdosu(idose)do j=
Page 67 and 68: if (maxz.gt.LKMAX) thenwrite(ifto,'
Page 69 and 70: zbound(in+1)=zbound(in)+widthend do
Page 71 and 72: if(idosu(idgrp).le.0 .or. idosu(idg
Page 73 and 74: !--------------------------------if
Page 75 and 76: end ifend ifelse if (u(np) .lt. 0.0
show all

Lecture Notes of Dose distribution calculation inside phantom with ...

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?