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