cg Version

tdet=7.62rdet=3.81tcov=0.1rtcov=0.1tgap=0.5rtgap=0.5write(6,330) tdet,rdet,tcov,rtcov,tgap,rtgap330 FORMAT(/' Detector length=',G15.5,' cm'/* ' Detector radius=',G15.5,' cm'/* ' Al cover thickness=',G10.2,' cm'/* ' Al cover side thickness=',G10.2,' cm'/* ' Front gap =',G10.2,' cm'/' Side gap =',G10.2,' cm'/)tdet=7.62rdet=3.81tcov=0.5rtcov=0.5write(6,330) tdet,rdet,tcov,rtcov330 FORMAT(/' Detector length=',G15.5,' cm'/* ' Detector radius=',G15.5,' cm'/* ' Al cover thickness=',G10.2,' cm'/* ' Al cover side thickness=',G10.2,' cm'/) ! --------------------------! Pulse height distribution! --------------------------write(6,370)! ----------------------! Absorbed energy in air! ----------------------avab = depes/ncountdepe2s = depe2s/ncountsigab = sqrt((depe2s - avab*avab)/ncount)write(6,362) avab,sigab362 FORMAT(' Absorbed energy in air =',G10.3,'+-',G9.2,' MeV/photon')avab = avab /33.97D-6 *1.602D-19sigab= sigab /33.97D-6 *1.602D-19write(6,364) avab,sigab364 FORMAT(' Output current =',G10.3,'+-',G9.2,' C/photon')! --------------------------! Pulse height distribution! --------------------------write(6,370) ausgab if (med(irl) .eq. 1) thendepe = depe + edepwtif (irl .eq. 1) thendepe = depe + edepwt23