| Geant4 Cross Reference |
1 macro egs4.kumac 1 macro egs4.kumac
2 * 2 *
3 set * ; opt * 3 set * ; opt *
4 * 4 *
5 filecase keep 5 filecase keep
6 6
7 set *fon -40 7 set *fon -40
8 set mscf 1.5 8 set mscf 1.5
9 9
10 set ygti 0.5 ; set gsiz 0.35 ; set ymgu 1.5 10 set ygti 0.5 ; set gsiz 0.35 ; set ymgu 1.5
11 set vsiz 0.25 ; set yval 0.3 ; set ylab 0.55 11 set vsiz 0.25 ; set yval 0.3 ; set ylab 0.55
12 set asiz 0.25 ; set yhti 0.9 ; set xlab 1.2 12 set asiz 0.25 ; set yhti 0.9 ; set xlab 1.2
13 13
14 title '30 GeV e- on 20 X0 iron' 14 title '30 GeV e- on 20 X0 iron'
15 h/file 1 9.0p1.hbook 15 h/file 1 9.0p1.hbook
16 16
17 set dmod 1 ; set hcol 4 ; set plci 4 17 set dmod 1 ; set hcol 4 ; set plci 4
18 h/pl //lun1/4 18 h/pl //lun1/4
19 key 14 9.0 1 'Geant4-9.0p1' ! l 19 key 14 9.0 1 'Geant4-9.0p1' ! l
20 20
21 exec data 21 exec data
22 22
23 set dmod 2 ; set hcol 2 ; set plci 2 23 set dmod 2 ; set hcol 2 ; set plci 2
24 h/pl 1 s 24 h/pl 1 s
25 key 14 8.0 2 'EGS4' ! l 25 key 14 8.0 2 'EGS4' ! l
26 26
27 atitle t '(100/E0) (dE/dt)' 27 atitle t '(100/E0) (dE/dt)'
28 exec Redraw-Axis 28 exec Redraw-Axis
29 29
30 close 1 30 close 1
31 31
32 return 32 return
33 33
34 macro data 34 macro data
35 * 35 *
36 * e- 30 GeV on 20 X0 Iron 36 * e- 30 GeV on 20 X0 Iron
37 * Particle Data Group. Physics Review D 50-3 ( 37 * Particle Data Group. Physics Review D 50-3 (1994)
38 * longitudinal energy deposition 38 * longitudinal energy deposition
39 * (from 0 to 20 X0 by bin of 0.5 X0 --> 40 bin 39 * (from 0 to 20 X0 by bin of 0.5 X0 --> 40 bins)
40 * 40 *
41 * *** read data 41 * *** read data
42 nbin=40 42 nbin=40
43 vector/create vect(2,[nbin]) 43 vector/create vect(2,[nbin])
44 vector/read vect pdg.ascii 44 vector/read vect pdg.ascii
45 **vector/write vect ! (1x,f6.1,3x,e12.6) 45 **vector/write vect ! (1x,f6.1,3x,e12.6)
46 * 46 *
47 * *** create and fill histogram 47 * *** create and fill histogram
48 Emin= 0. 48 Emin= 0.
49 Emax= 20. 49 Emax= 20.
50 create/1dhisto 1 '20 X0 iron. e- GeV' [nbin] 50 create/1dhisto 1 '20 X0 iron. e- GeV' [nbin] [Emin] [Emax]
51 put/contents 1 vect(2) 51 put/contents 1 vect(2)
52 52
53 **h/pl 1 53 **h/pl 1
54 54
55 * *** save histo 55 * *** save histo
56 * 56 *
57 **h/file 111 pdg.hbook ! n 57 **h/file 111 pdg.hbook ! n
58 **hrout 0 58 **hrout 0
59 59
60 return 60 return
61 61
62 Macro Redraw-Axis 62 Macro Redraw-Axis
63 set plci 1 63 set plci 1
64 X1 = $GRAFINFO('WNXMIN') 64 X1 = $GRAFINFO('WNXMIN')
65 X2 = $GRAFINFO('WNXMAX') 65 X2 = $GRAFINFO('WNXMAX')
66 Y1 = $GRAFINFO('WNYMIN') 66 Y1 = $GRAFINFO('WNYMIN')
67 Y2 = $GRAFINFO('WNYMAX') 67 Y2 = $GRAFINFO('WNYMAX')
68 IS = $GRAFINFO('?FAIS') 68 IS = $GRAFINFO('?FAIS')
69 IC = $GRAFINFO('?FACI') 69 IC = $GRAFINFO('?FACI')
70 Set FAIS 0 ; Set FACI 1 70 Set FAIS 0 ; Set FACI 1
71 Box [X1] [X2] [Y1] [Y2] 71 Box [X1] [X2] [Y1] [Y2]
72 Set FAIS [IS] ; Set FACI [IC] 72 Set FAIS [IS] ; Set FACI [IC]
73 Tick XL [X1] ! 73 Tick XL [X1] !
74 * Tick YB ! [Y1] 74 * Tick YB ! [Y1]
75 Return 75 Return
76 76