Modify statements related geometry output.Changetdet=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'/)totdet=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'/) Add routines to calculated average absorbed energy of air and its statistical error.Change! --------------------------! Pulse height distribution! --------------------------write(6,370)to! ----------------------! Absorbed energy in air! ----------------------avab = depes/ncountdepe2s = depe2s/ncountsigab = sqrt((depe2s - avab*avab)/ncount)write(6,362) avab,sigab362 FORMAT(' Absorbed energy in air =',G15.5,'+-',G15.5,' MeV/photon')avab = avab /33.97D-6 *1.602D-19sigab= sigab /33.97D-6 *1.602D-19write(6,364) avab,sigab364 FORMAT(' Output current =',G15.5,'+-',G15.5,' C/photon')! --------------------------! Pulse height distribution! --------------------------write(6,370) Modify score condition at the detector region in ausgab to distinguish from air outsidedetector.toif (med(irl) .eq. 1) thendepe = depe + edepwtif (irl .eq. 1) thendepe = depe + edepwt5. Modify ucion<strong>cg</strong>v.data as follows:47
ROuter vacuum (region 4)Air region (region 3)Al Cover (region 2)Air (region 1)1.253MeVPhoton13.12Z-0.5 0.0 8.127.62-5.5Figure 5: Geometry of ucion<strong>cg</strong>v.fRCC 1 0.00 0.0 0.0 0.00 0.07.62 3.81RCC 2 0.00 0.0 -0.5 0.00 0.08.62 4.31RCC 3 0.00 0.0 -5.5 0.00 0.018.62 9.31RCC 4 0.00 0.0 -6.0 0.00 0.020.62 10.31ENDZ1 +1Z2 +2 -1Z3 +3 -2Z4 +4 -3END1 2 1 06. Check ucion<strong>cg</strong>v.data. cp ucion<strong>cg</strong>v.data ucion<strong>cg</strong>v.geo. Check ucion<strong>cg</strong>v.geo by using CGView as followsSelect "Making geometry data" of File option.Select Open File and assign ucion<strong>cg</strong>v.geo".Geometry is displayed when you select OK.Select "Geomtry Check" of Environment option.Select "Check Start".7. Modify ucion<strong>cg</strong>v.inp as follows.MIXT48
- Page 1 and 2: KEK Internal 2011-5December 2011R/D
- Page 3 and 4: Lecture Notes ofResponse calculatio
- Page 5 and 6: 3 Details of user code 313.1 Main p
- Page 7 and 8: 1 Combinatrial Geometry (CG)2 Combi
- Page 9 and 10: ody 3 2 body + ; OR B23A1Figure 2
- Page 11 and 12: Z4 +4Z5 +5 -3Z6 +6 -5END1 0 2 3 4 0
- Page 13 and 14: end if! ==============call block_se
- Page 15 and 16: ! pdf data for many sourcedeltae=0.
- Page 17 and 18: ! ---------------------------------
- Page 19 and 20: sigph=dsqrt((ph2s(ie)-avph*avph)/nc
- Page 21 and 22: 6 6.1 NaI 1. Cs-137 (0.662MeV) 2.
- Page 23 and 24: ekein=1.253ekein=1.333{ ! ---------
- Page 25 and 26: ROuter vacuum (region 6)Air region
- Page 27 and 28: nmed=2 medarr(1)='NAImedarr(2)='ALm
- Page 29 and 30: ROuter vacuum (region 4)Air region
- Page 31 and 32: EGS5 sample user code (ucnaicgv.f)R
- Page 33 and 34: operator, it means that the region
- Page 35 and 36: 2 Outlines of sample user code uccg
- Page 37 and 38: include 'auxcommons/edata.f'include
- Page 39 and 40: The material assignment is read in
- Page 41 and 42: ! Write batch numberwrite(39,fmt="(
- Page 43 and 44: ! ---------------------------------
- Page 45 and 46: a distance to the boundary of regio
- Page 47 and 48: { Add sampling routines for source
- Page 49 and 50: ROuter vacuum (region 6)Air region
- Page 51: to* xi0,yi0,zi0,avab,depes,depe2s,s
- Page 55 and 56: Appendix: Full listings of ucnaicgv
- Page 57 and 58: open(39,FILE='egs5job.pic',STATUS='
- Page 59 and 60: write(6,140)140 FORMAT(/,' HATCH-ca
- Page 61 and 62: ilines = ilines + 1write(6,280) eto
- Page 63 and 64: * 14X,' Positron')do ie=1,50elow=de
- Page 65 and 66: ! ---------------------------------
- Page 67 and 68: ccyidd=y_np+tval0*v_npzidd=z_np+tva
Change language