1 gaskelld 1.1 !!!SF is 40 Em bins of 5 MeV (0-200 MeV) by 40 Pm bins of 20 MeV (0-800 MeV/c)
 2              
 3              !SF is 80 Em bins of 5 MeV             by 40 Pm bins of 20 MeV.
 4              
 5              destroy *
 6              read\scalar benharsf_?1.dat pbins ebins
 7              read benharsf_?1.dat\(2) pmbin embin sfpbin sfnbin dp de
 8              
 9              stat\nomess sfpbin normfacp\sum
10              stat\nomess sfnbin normfacn\sum
11              sfpbin=sfpbin/normfacp
12              sfnbin=sfnbin/normfacn
13              
14              sort embin pmbin sfpbin sfnbin
15              em = embin[1:#:pbins]
16              rebin sfpbin sfp_em pbins
17              rebin sfnbin sfn_em pbins
18              
19              sort pmbin embin sfpbin sfnbin
20              pm = pmbin[1:#:ebins]
21              rebin sfpbin sfp_pm ebins
22 gaskelld 1.1 rebin sfnbin sfn_pm ebins
23              
24              clear
25              portrait
26              window 3
27              set ylog 10
28              xlabel `Pm'
29              pch 2 .8 1
30              auto on
31              scales 0 1000 5 -4 0 4
32              gr pm sfp_pm
33              pch 12 .8 2
34              cur pm sfn_pm
35              
36              copy sfp_pm sfpsum_pm
37              copy sfn_pm sfnsum_pm
38              do ind=[2:len(pm):1]
39               stat\nomess sfp_pm[1:ind] sfpsum_pm[ind]\sum
40               stat\nomess sfn_pm[1:ind] sfnsum_pm[ind]\sum
41              enddo
42              sfpsum_pm=sfpsum_pm/sfpsum_pm[len(pm)]
43 gaskelld 1.1 sfnsum_pm=sfnsum_pm/sfnsum_pm[len(pm)]
44              
45              set %xloc 40
46              set %yloc 80
47              set cursor -7
48              text `Raw Proton Norm. = '//rchar(nint(10000*normfacp)/10000)
49              set %yloc 75
50              text `Raw Neutron Norm. = '//rchar(nint(10000*normfacn)/10000)
51              
52              window 4
53              xlabel `Em'
54              set ylog 10
55              scales 0 250 5 -4 0 0
56              pch 2 .8 1
57              gr em sfp_em
58              pch 12 .8 2
59              cur em sfn_em
60              
61              !terminal
62              
63              !cl
64 gaskelld 1.1 !set ylog 10
65              !landscape
66              !window 0
67              !nbins=10		!must be mult of 40
68              !scales 0 400 8 -8 0 8
69              !graph\axesonly
70              !do ind=[1:nbins:1]
71              ! pmlo=nint(800/nbins)*(ind-1)+1
72              ! pmhi=nint(800/nbins)*ind
73              ! destroy x y
74              ! copy embin sfpbin xtmp ytmp IFF (pmbin>=pmlo & pmbin<=pmhi)
75              ! sort xtmp ytmp
76              ! x=xtmp[1:#:nint(40/nbins)]
77              ! rebin ytmp y nint(40/nbins)
78              ! pch ind 1 ind
79              ! destroy x y IFF (y<1.d-30)
80              
81              ! stat\nomess y ysum\sum
82              ! y=y/ysum		!fixed normalization
83              
84              ! cur x y
85 gaskelld 1.1 !enddo