macro plot_fluka e=200 theta=90
* 
* Plot parameterizations of Bcal dynamic range
* 05/08/07  ES
*
* set options
*
option ndate
option nbox
*set stat 1111111
set stat 111
option stat
set fit 111
option fit
option grid
*
* plotting options
*
set *
set xmgl 4.
set ymgl 4.
set asiz 0.5
set xlab 2.
set ylab 1.
set xsiz 20.
set xmgl 3.
set ymgl 3.
set ysiz 20.
set gsiz 0.4
* set xwin 0.1
*
* set font definitions to bold roman 
*
set CFON -21
set GFON -21
set LFON -21
set TFON -21
set VFON -21
set txfp -21
set SMGU 0.02
set SMGR 0.02
set CSIZ 0.45
set VSIZ 0.3
set TSIZ 0.3
set YHTI 0.9
set HWID 3.0
set BWID 3.0
* set YHTI 2.
*
* open inputfile
*
* e = 2000
* theta = 20
infile = 'fluka_de_phot_' // [theta] // 'deg_' // [e] // 'mev.hbook'
hi/file 2 [infile] 0 -x
hi/create/title_global 'Energy in Bcal Segments E?[g]!=2000 MeV'
*
zone 2 2
*
* open metafile
*
* outfile = 'plot_fluka_' // [theta] // 'deg_' // [e] // 'mev.ps'
outfile = 'plot_fluka_' // [theta] // 'deg_' // [e] // 'mev.ps'
mess ' e=' [e] 'theta=' [theta] ' outfile=' [outfile]
for/file 66 [outfile]
meta 66 -111
* outfile = 'plot_fluka_' // [theta] // 'deg_' // [e] 'mev.eps'
*for/file 66 [outfile]
*meta 66 -113
*
npts1 = 101
*
xmin = 0 
xmax = 10
ymin = 10
ymax = 70
* option logx 
* option logy
*
npts = [npts1] - 1
sigma j=array([npts1],0#[npts1]-1)
sigma PDE = [xmin] + ([xmax]-[xmin])*j/[npts]
*
*
message 'wait'
* wait
*
label 1 1 ' '
csize = 0.20
igset chhe [csize]
set ndvx 505
set ndvy 505
* hplot/null [xmin] [xmax] [ymin] [ymax]
* hplot/atitle  'Energy(MeV)' 'Counts'
set pmci 2
set plci 2
*
jmax = 0
do j=1,[jmax]
nt/plot 77.Ed(1,[j])
nt/plot 77.Ed(2,[j])
nt/plot 77.Ed(3,[j])
nt/plot 77.Ed(4,[j])
* wait
enddo
*
do j=1,[jmax]
nt/plot 77.Edr(1,[j])
nt/plot 77.Edr(2,[j])
nt/plot 77.Edr(3,[j])
nt/plot 77.Edr(4,[j])
* wait
enddo
*
do j=1,[jmax]
nt/plot 77.Edl(1,[j])
nt/plot 77.Edl(2,[j])
nt/plot 77.Edl(3,[j])
nt/plot 77.Edl(4,[j])
* wait
enddo
*
do j=1,[jmax]
nt/plot 77.Per(1,[j])
nt/plot 77.Per(2,[j])
nt/plot 77.Per(3,[j])
nt/plot 77.Per(4,[j])
* wait
enddo
*
do j=1,[jmax]
nt/plot 77.Pel(1,[j])
nt/plot 77.Pel(2,[j])
nt/plot 77.Pel(3,[j])
nt/plot 77.Pel(4,[j])
* wait
enddo
*
do j=1,[jmax]
nt/plot 77.Ed(1,[j])+Ed(2,[j])+Ed(3,[j])+Ed(4,[j])
enddo
*
kmax=5
do k=1,[kmax]
jmax = 4
do j=1,[jmax]
nt/plot 77.plot_fluka.f([k],[j])
enddo
*
exe window#push
igset chhe 0.3
if ([k].eq.1) then 
   itx -0.2 -0.3 'Ed'
elseif ([k] .eq. 2) then
   itx -0.2 -0.3 'Edr'
elseif ([k] .eq. 3) then
   itx -0.2 -0.3 'Edl'
elseif ([k] .eq. 4) then
   itx -0.2 -0.3 'Per'
elseif ([k] .eq. 5) then
   itx -0.2 -0.3 'Pel'
endif
exe window#pop
*
jmax = 8
do j=5,[jmax]
nt/plot 77.plot_fluka.f([k],[j])
enddo
*
exe window#push
igset chhe 0.3
if ([k].eq.1) then 
   itx -0.2 -0.3 'Ed'
elseif ([k] .eq. 2) then
   itx -0.2 -0.3 'Edr'
elseif ([k] .eq. 3) then
   itx -0.2 -0.3 'Edl'
elseif ([k] .eq. 4) then
   itx -0.2 -0.3 'Per'
elseif ([k] .eq. 5) then
   itx -0.2 -0.3 'Pel'
endif
exe window#pop
*
jmax = 14
do j=11,[jmax]
nt/plot 77.plot_fluka.f([k],[j])
enddo
*
exe window#push
igset chhe 0.3
if ([k].eq.1) then 
   itx -0.2 -0.3 'Ed'
elseif ([k] .eq. 2) then
   itx -0.2 -0.3 'Edr'
elseif ([k] .eq. 3) then
   itx -0.2 -0.3 'Edl'
elseif ([k] .eq. 4) then
   itx -0.2 -0.3 'Per'
elseif ([k] .eq. 5) then
   itx -0.2 -0.3 'Pel'
endif
exe window#pop
*
*
enddo
*
do j=1,5
  nt/plot 77.plot_fluka.f([j],100)
  nt/plot 77.plot_fluka.f([j],100)%plot_fluka(1,100)
enddo
nt/plot 77.plot_fluka.f(4,100)%plot_fluka(5,100)
*
jmax = 0
do j=1,[jmax]
nt/plot 77.Ed([j],1)+Ed([j],2)+Ed([j],3)+Ed([j],4)+Ed([j],5)+Ed([j],6)+Ed([j],7)+Ed([j],8)
enddo
zone
* nt/plot 77.Ed(1,1)+Ed(1,2)+Ed(1,3)+Ed(1,4)+Ed(1,5)+Ed(1,6)+Ed(1,7)+Ed(1,8)+Ed(2,1)+Ed(2,2)+Ed(2,3)+Ed(2,4)+Ed(2,5)+Ed(2,6)+Ed(2,7)+Ed(2,8)+Ed(3,1)+Ed(3,2)+Ed(3,3)+Ed(3,4)+Ed(3,5)+Ed(3,6)+Ed(3,7)+Ed(3,8)+Ed(4,1)+Ed(4,2)+Ed(4,3)+Ed(4,4)+Ed(4,5)+Ed(4,6)+Ed(4,7)+Ed(4,8)
*
*
* project and plot histograms
*
hi/del 101
hi/create/1d 101 'Total Energy' 100 0 0.2*[e]
cd //lun2
nt/plot 101 77.plot_fluka.f(1,100)
hi/plot 101
*
close 66
*
exitm 
return