pro flux_s_histo,help=help

if keyword_set(help) then begin
	print,'PURPOSE:'
	print,'	Plot histograms of flux and separation'
	print,'HELP:'
	print,'	pro flux_s_histo'
	print,'	no inputs, only outputs are ps files flux_hist.ps and s_hist.ps'
	print,'AUTHOR:'
	print,'	Rick Pernak, Goddard Space Flight Center'
	return
endif

openr,1,'/home/pernak/RHESSI/paper/spreadsheet/mem.csv'
mem = strarr(25)
readf,1,mem
close,1

openr,1,'/home/pernak/RHESSI/paper/spreadsheet/clean.csv'
clean = strarr(25)
readf,1,clean
close,1

openr,1,'/home/pernak/RHESSI/paper/spreadsheet/pixon.csv'
pixon = strarr(25)
readf,1,pixon
close,1

mfluxr = dblarr(25)
cfluxr = dblarr(25)
pfluxr = dblarr(25)
smem = dblarr(25)
sclean = dblarr(25)
spix = dblarr(25)
;m: mem, c:clean, p:pixon, r:ratio, s:separation angle

for ctr=0,24 do begin
	msplit = str_sep(mem[ctr],'	')
	mfluxr[ctr] = msplit[4]
	smem[ctr] = msplit[13]

	csplit = str_sep(clean[ctr],'	')
	cfluxr[ctr] = csplit[4]
	sclean[ctr] = csplit[13]

	psplit = str_sep(pixon[ctr],'	')
	pfluxr[ctr] = psplit[4]
	spix[ctr] = psplit[13]
	
endfor

set_plot,'ps'
device,filename='flux_hist.ps'
multiplot,/def
multiplot,[0,1,3,0,0]
plothist,mfluxr,bin=.1,yr=[0,15],/half
legend,['MEM'],box=0
multiplot
plothist,cfluxr,bin=.1,yr=[0,15],ytickformat='(A1)',/half
axis,yaxis=1,yticks=4,ytickv=[0,5,10,15]
legend,['Clean'],box=0
multiplot
plothist,pfluxr,bin=.1,yr=[0,15],xtitle='Flux1/(Flux1+Flux2)',/half
legend,['Pixon'],box=0
multiplot,/reset
multiplot,/def
device,/close

set_plot,'ps'
device,filename='s_hist.ps'
multiplot,/def
multiplot,[0,1,3,0,0]
plothist,smem,bin=5,xr=[0,50],yr=[0,15]
legend,['MEM'],box=0
multiplot
plothist,sclean,bin=5,xr=[0,50],yr=[0,15],ytickformat='(A1)'
axis,yaxis=1,yticks=4,ytickv=[0,5,10,15]
legend,['Clean'],box=0
multiplot
plothist,spix,bin=5,xr=[0,50],yr=[0,15],xtitle='s (Mm)'
legend,['Pixon'],box=0
multiplot,/reset
multiplot,/def
device,/close
set_plot,'x'

return
end