In the case <strong>of</strong> DOSrun5again ucxyz <strong>phantom</strong>35. Confirm the density at the lung regions. Compare the absorbed dose <strong>distribution</strong><strong>with</strong> xyz pantom.out.6. Check the trajectories using CGview.4.4 Problem 41. cp ucxyz <strong>phantom</strong>3.data ucxyz <strong>phantom</strong>4.data2. cp ucxyz <strong>phantom</strong>3.inp ucxyz <strong>phantom</strong>4.inp3. Modify ucxyz <strong>phantom</strong>4.data as follows.(a) Modify the density at tumor regions.Change1,3,1,3, 5,14, 1, 0.300, 0.00, 0.00 lung1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impacto1,3,1,3, 5, 7, 1, 0.300, 0.00, 0.00 lung1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3, 8, 9, 1, 0.000, 0.00, 0.00 tumor1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3,10,14, 1, 0.300, 0.00, 0.00 lung1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac4. Execute run5again.Enter ucxyz <strong>phantom</strong>4 as the file name for unit 4 and simply enter "return" asthe file name for unit 25.5. Enter simply return for "Enter name <strong>of</strong> the executable".6. Confirm the density at the tumor regions. Compare the absorbed dose <strong>distribution</strong><strong>with</strong> xyz pantom.out.4.5 Ploblem 51. cp ucxyz <strong>phantom</strong>2.f ucxyz <strong>phantom</strong>5.f2. cp ucxyz <strong>phantom</strong>4.data ucxyz <strong>phantom</strong>5.data3. cp ucxyz <strong>phantom</strong>4.inp ucxyz <strong>phantom</strong>5.inp4. Modify ucxyz <strong>phantom</strong>5.f as follows: Increase number <strong>of</strong> medium.ChangetoMXMED=2MXMED=35. Modify ucxyz <strong>phantom</strong>4.data as follows.(a) Modify the number <strong>of</strong> materials from '2' to '3'.Changeto2 nmed (I10)3 nmed (I10)46
(b) Add new material (Iron).ChangeAIR-AT-NTPtoAIR-AT-NTPFE(c) Add characteristic dimension <strong>of</strong> iron.Changeto1.0, 1.0 chard1.0, 1.0, 1.0 chard(d) Modify material assignment etc.Changemedia(j,2) (24A1)media(j,2) (24A1)media(j,3) (24A1)related the "tumor" regions.1,3,1,3, 2, 4, 1, 0.000, 0.00, 0.00 tissue1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3, 5, 7, 1, 0.300, 0.00, 0.00 lung1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3, 8, 9, 1, 0.000, 0.00, 0.00 tumor1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3,10,14, 1, 0.300, 0.00, 0.00 lung1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3,15,17, 1, 0.000, 0.00, 0.00 tissue1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impacto1,3,1,3, 2, 6, 1, 0.000, 0.00, 0.00 tissue1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3, 7, 7, 3, 0.000, 0.00, 0.00 iron1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac1,3,1,3, 8,21, 1, 0.000, 0.00, 0.00 tissue1 1 0 0 0 0 0 0 peang,edgefl,auger,ray,pola,incoh,pr<strong>of</strong>,impac6. Modify ucxyz <strong>phantom</strong>5.inp as follows.(a) Add input data for iron.ELEM&INP IRAYL=1 /ENDFEFEFEENER&INP AE=0.521,AP=0.010,UE=0.711,UP=0.2 /ENDPWLF&INP /ENDDECK&INP /END7. Run ucxyz <strong>phantom</strong>5.f by egs5run. In the case <strong>of</strong> Linux or CygwinEnter ucxyz <strong>phantom</strong>5 as the user code.Simply enter "return" as the file name for unit 4 and 25.Enter 1 for "Does this user code read from the terminal?". In the case <strong>of</strong> DOSegs5run ucxyz <strong>phantom</strong>58. Check egs5job.out to confirm proper setting <strong>of</strong> iron region. Compare the absorbeddose <strong>distribution</strong> <strong>with</strong> xyz <strong>phantom</strong>.out.9. Check the trajectories using CGview to confirm almost all photons stopping atthe iron region.47
- 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: &INP AE=0.521,AP=0.0100,UE=2.011,UP
- 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 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