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

More documents

Recommendations

Info

Appendix: Full listings of ucxyz phantom.f48
!***********************************************************************!***************************** KEK, High Energy Accelerator Research *!***************************** Organization *!* u c x y z _ p h a n t o m* *!***************************** EGS5.0 USER CODE - 19 Nov 2012/1300 *!***********************************************************************!* This is a general User Code based on the cg geometry scheme. *!***********************************************************************! *! PROGRAMMERS: H. Hirayama *! Radiation Science Center *! Applied Science Laboratory *! KEK, High Energy Accelerator Research Organization *! 1-1, Oho, Tsukuba, Ibaraki, 305-0801 *! Japan *! *! E-mail: hideo.hirayama@kek.jp *! Telephone: +81-29-864-5489 *! Fax: +81-29-864-1993 *! Based on xyzdos.mor. *! *! 09AUG2004: CDF creation parts are corrected. H. Hirayama *!***********************************************************************!***********************************************************************! The ucxyz_phantom.f User Code requires a data-input file *! (e.g., ucxyz_phantom.data) that is read by subroutine getvoxel (with *! instructions in its header). The following shows the geometry for *! ucxyz_phantom.data. *! Trajectory data can be defined for PICT32 or CGview by setting *! npreci 1, 2 or 3, respectively. *! This user code corresponds to ucphantom_rec.mor for egs4. *! Use Ranlux random number generator. *!***********************************************************************! *! --------------------------------------------------- *! 3-Dimensional X-Y-Z Geometry (ucxyz_pahtom example) *! --------------------------------------------------- *! *! X (Y into page) *! ^ *! | *! | *! | Outer vacuum region *! +----+---+---+- -+---+----+ 15.0cm *! Vacuum | | | | | | | *! region +Air +H2O|Water (H2O) |Air | *! | | | | |H2O| | *! +----+---+---+- -+---+----+-- 0.5 *! | | | | | | | *! |Air |H2O|H2O| |H2O|Air | *! 1.253MeV photon| | | | | | | *! =========>+----+---+---+- -+---+----+-----> Z *! photons -5.0 0.0 1.0 19.0 20.0 *! 25.0 *! *! *!***********************************************************************!23456789|123456789|123456789|123456789|123456789|123456789|123456789|12!!-----------------------------------------------------------------------!------------------------------- main code -----------------------------!-----------------------------------------------------------------------!-----------------------------------------------------------------------! Step 1: Initialization!-----------------------------------------------------------------------implicit none! ------------! EGS5 COMMONs! ------------include 'include/egs5_h.f'include 'include/egs5_bounds.f'include 'include/egs5_edge.f'include 'include/egs5_elecin.f'include 'include/egs5_media.f'include 'include/egs5_misc.f'include 'include/egs5_stack.f'include 'include/egs5_thresh.f'include 'include/egs5_uphiot.f'! Main EGS "header" file49
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: (b) Add new material (Iron).ChangeA
Page 57 and 58: !----------------------------------
Page 59 and 60: faexp=0.D0faexps=0.D0faexp2s=0.D0fe
Page 61 and 62: ! =======================if (iwatch
Page 63 and 64: do k=kdosl(idose),kdosu(idose)do j=
Page 65 and 66: ! Record 10a ipeangsw, Switches for
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 ...

Create successful ePaper yourself

Delete template?

Save as template?