Warning, /geant4/examples/extended/electromagnetic/TestEm2/egs4/sum.kumac is written in an unsupported language. File is not indexed.
0001 macro sum.kumac
0002 *
0003 set * ; opt *
0004 *
0005 filecase keep
0006
0007 set *fon -40
0008 set mscf 1.5
0009
0010 set ygti 0.5 ; set gsiz 0.35 ; set ymgu 1.5
0011 set vsiz 0.25 ; set yval 0.3 ; set ylab 0.55
0012 set asiz 0.25 ; set yhti 0.9 ; set xlab 1.2
0013
0014 title '30 GeV e- on 20 X0 iron'
0015 h/file 1 9.0p1.hbook
0016
0017 set dmod 1 ; set hcol 4 ; set plci 4
0018 h/pl //lun1/6
0019 key 14 60.0 1 'Geant4-9.0p1' ! l
0020
0021 exec data
0022
0023 set dmod 2 ; set hcol 2 ; set plci 2
0024 h/pl 2 s
0025 key 14 55.0 2 'EGS4' ! l
0026
0027 atitle t '100 Edep/E0'
0028 exec Redraw-Axis
0029
0030 close 1
0031
0032 return
0033
0034 macro data
0035 *
0036 * e- 30 GeV on 20 X0 Iron
0037 * Particle Data Group. Physics Review D 50-3 (1994)
0038 * cumulative longitudinal energy deposition
0039 * (from 0 to 20 X0 by bin of 0.5 X0 --> 40 bins)
0040 *
0041 * *** read data
0042 nbin=40
0043 vector/create vect(2,[nbin])
0044 vector/read vect sum.ascii
0045 **vector/write vect ! (1x,f6.1,3x,e12.6)
0046 *
0047 * *** create and fill histogram
0048 Emin= 0.25
0049 Emax= 20.25
0050 create/1dhisto 2 '20 X0 iron. e- GeV' [nbin] [Emin] [Emax]
0051 put/contents 2 vect(2)
0052
0053 **h/pl 2
0054
0055 * *** save histo
0056 *
0057 **h/file 222 pdgsum.hbook ! n
0058 **hrout 0
0059
0060 return
0061
0062 Macro Redraw-Axis
0063 set plci 1
0064 X1 = $GRAFINFO('WNXMIN')
0065 X2 = $GRAFINFO('WNXMAX')
0066 Y1 = $GRAFINFO('WNYMIN')
0067 Y2 = $GRAFINFO('WNYMAX')
0068 IS = $GRAFINFO('?FAIS')
0069 IC = $GRAFINFO('?FACI')
0070 Set FAIS 0 ; Set FACI 1
0071 Box [X1] [X2] [Y1] [Y2]
0072 Set FAIS [IS] ; Set FACI [IC]
0073 Tick XL [X1] !
0074 * Tick YB ! [Y1]
0075 Return
0076
