cg Version

include 'include/egs5_h.f'include 'include/egs5_epcont.f'include 'include/egs5_misc.f'include 'include/egs5_stack.f'include 'include/egs5_useful.f'include 'auxcommons/aux_h.f'include 'auxcommons/etaly1.f'include 'auxcommons/lines.f'include 'auxcommons/ntaly1.f'include 'auxcommons/watch.f'common/totals/* depe,deltae,spec(3,50),maxpictreal*8 depe,deltae,specinteger maxpictinteger* iargreal*8* edepwtinteger* ie,iql,irl,ntype! ------------------------! Set some local variables! ------------------------irl = ir(np)iql = iq(np)edepwt = edep*wt(np)! Main EGS "header" file! COMMONs required by EGS5 code! Auxiliary-code "header" file! Auxiliary-code COMMONs! Variables to score! Arguments! Local variables! -----------------------------------------------------------! Keep track of energy deposition (for conservation purposes)! -----------------------------------------------------------if (iarg .lt. 5) thenesum(iql+2,irl,iarg+1) = esum(iql+2,irl,iarg+1) + edepwtnsum(iql+2,irl,iarg+1) = nsum(iql+2,irl,iarg+1) + 1end if! -----------------------------------------------------------------! Print out particle transport information (if switch is turned on)! -----------------------------------------------------------------! ========================if (iwatch .gt. 0) call swatch(iarg,iwatch)! ========================if(iarg .ge. 5) return! ----------------------------------------------! Score energy deposition inside NaI detector! ----------------------------------------------if (med(irl) .eq. 1) thendepe = depe + edepwt! ----------------------------------------------------! Score particle information if it enters from outside! ----------------------------------------------------if (irl .ne. irold .and. iarg .eq. 0) thenif (iql .eq. 0) then! photonntype=1ie = e(np)/deltae +1if(ie .gt. 50) ie = 50elseif (iql .eq. -1) then ! electronntype=2ie = (e(np) - RM)/deltae +1if(ie .gt. 50) ie = 50else! positronntype=3ie = (e(np) - RM)/deltae +1if(ie .gt. 50) ie = 50end ifspec(ntype,ie) = spec(ntype,ie) + wt(np)end ifend if59