THE EGS5 CODE SYSTEM

pnorm(1,2) =0.0pnorm(2,2) =0.0pnorm(3,2)= 1.0! Normal vectors!-----------------------------------------------------------------------! Step 7: Initialization-for-ausgab!-----------------------------------------------------------------------do i=1,50ebin(i) = 0.0! Zero scoring array before startingend dobwidth = 0.002!-----------------------------------------------------------------------! Step 8: Shower-call!-----------------------------------------------------------------------! Initiate the shower ncase timesncase=10000do i=1,NCASEcall shower(iqin,ein,xin,yin,zin,uin,vin,win,irin,wtin)end do!-----------------------------------------------------------------------! Step 8: Output-of-results!-----------------------------------------------------------------------! Use log10(10000.0) as maximum valuebinmax=dlog10(10000.d0)if (iedgfl(2).eq.1) thenwrite(6,160) ein,pcoord(3,2)160 format(/’ Reflected photon spectrum’/’ for a’,F8.2,* ’ MeV pencil beam of photons on a’,F7.2,* ’ cm thick slab of lead’/’ with fluorescence photon’//T6,* ’Energy counts/incident photon’/* 25X,’ log(counts for 10^4 incident photons)’)elsewrite(6,170) ein,pcoord(3,2)170 format(’ Reflected photon spectrum’/’ for a’,F8.2,* ’ MeV pencil beam of photons on a’,F7.2,* ’ cm thick slab of lead’/’ without fluorescence photon’//T6,* ’Energy counts/incident photon’/* 25X,’ log(counts for 10^4 incident photons)’)end ifdo j=1,48line(j)=’ ’end do! Blank entire output arraydo j=1,50if(ebin(j).gt.0) then182