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

! Source position from phantom surface in cm.sposi=10.0nsebin=2! Number of source energy binsread(2,*) (esbin(i),i=1,nsebin)read(2,*) (espdf(i),i=1,nsebin)!---------------------------! Calculate CDF from pdf!---------------------------tnum=0.D0do ie=1,nsebintnum=tnum+espdf(ie)end doescdf(1)=espdf(1)/tnumdo ie=2,nsebinescdf(ie)=escdf(ie-1)+espdf(ie)/tnumend do Modify the maximum electron kinetic energy used.Changeekein=1.253! Kinetic energy of source photontoekein=esbin(nsebin) ! Maximum kinetic energy} Add sampling routines for source photon energy sampling.Changeekin=ekeintocall randomset(rnnow)do ie=1,nsebinif(rnnow.le.escdf(ie)) go to 1000end do1000 ekin=esbin(ie) Modify output statement concerning the source energy.Change300 FORMAT(/' Absorbed energy inside phantom for 1.253MeV photon'/to300 FORMAT(/' Absorbed energy inside phantom for Co-60 photon'/5. Run ucxyz phantom1.f by egs5run. In the case of Linux or CygwinEnter ucxyz phantom1 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 of DOSegs5run ucxyz phantom16. Check egs5job.out to confirm average source energy is nearly equal to 1.253MeV.Compare the obtained results with xyz phantom.out.4.2 Problem 21. cp ucxyz phantom1.f ucxyz phantom2.f2. cp ucxyz phantom1.data ucxyz phantom2.data3. cp ucxyz phantom1.inp ucxyz phantom2.inp41