pro comp_alpha_plot,compare,help=help

if keyword_set(help) then begin
	print,'PURPOSE:'
	print,'	Plot old alpha (atan(y/x)) vs. new alpha (acos(...)) '
	print,' and old s vs. new s'
	print,'HELP:'
	print,'	pro comp_alpha_plot,compare'
	print,'	no inputs, only outputs are structure "compare," with all values '
	print,'	of alpha and s; and alpha_s_comparisons.ps file'
	print,'AUTHOR:'
	print,'	Rick Pernak, Goddard Space Flight Center'
	return
endif

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

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

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

;alpha: orientation angle, s:de-projected separation
;m: mem, c:clean, p:pixon

old_alpham=dblarr(25)
new_alpham=dblarr(25)
old_alphac=dblarr(25)
new_alphac=dblarr(25)
old_alphap=dblarr(25)
new_alphap=dblarr(25)

old_sm=dblarr(25)
new_sm=dblarr(25)
old_sc=dblarr(25)
new_sc=dblarr(25)
old_sp=dblarr(25)
new_sp=dblarr(25)

for ctr=0,24 do begin
	split = str_sep(mem_arr[ctr],'	')
	old = split[10] ;old individual alphas
	new = split[15] ;new individual alphas
	os = split[14]	;old individual s
	ns = split[13]	;new individual s
	old_alpham[ctr] = old
	new_alpham[ctr] = new
	old_sm[ctr] = os
	new_sm[ctr] = ns
endfor

for ctr=0,24 do begin
	split = str_sep(clean_arr[ctr],'	')
	old = split[10]		;old individual alphas
	new = split[15] 	;new individual alphas
	os = split[14]		;old individual s
	ns = split[13]		;new individual s
	old_alphac[ctr] = old
	new_alphac[ctr] = new
	old_sc[ctr] = os
	new_sc[ctr] = ns
endfor

for ctr=0,24 do begin
	split = str_sep(pixon_arr[ctr],'	')
	old = split[10]		;old individual alphas
	new = split[15] 	;new individual alphas
	os = split[14]		;old individual s
	ns = split[13]		;new individual s
	old_alphap[ctr] = old
	new_alphap[ctr] = new
	old_sp[ctr] = os
	new_sp[ctr] = ns
endfor

;comparison structure w/ all old and new values of alpha and s
compare= $
{proj_alpham:old_alpham,true_alpham:new_alpham,proj_smem:old_sm,true_smem:new_sm,$
proj_alphac:old_alphac,true_alphac:new_alphac,proj_sclean:old_sc,true_sclean:new_sc,$
proj_alphap:old_alphap,true_alphap:new_alphap,proj_spixon:old_sp,true_spixon:new_sp}

set_plot,'ps'
device,filename='alpha_s_comparisons.ps'

!p.multi = [0,3,2,0,0]

xr=[0,50]
yr=[0,50]
ytitle = 'Projected Separation (mM)'
xtitle = 'True Separation (mM)'

plot,old_sm,new_sm,xr=xr,yr=yr,psym=5,title='MEM s (Mm)',ytitle=ytitle,$
	xtitle=xtitle
oplot,xr,yr
for ctr=0,24 do begin
	xyouts,old_sm[ctr],new_sm[ctr],strcompress(ctr+1,/rem)
endfor

plot,old_sc,new_sc,xr=xr,yr=yr,psym=5,title='Clean s (Mm)',ytitle=ytitle,$
	xtitle=xtitle
oplot,xr,yr
for ctr=0,24 do begin
	xyouts,old_sc[ctr],new_sc[ctr],strcompress(ctr+1,/rem)
endfor

plot,old_sp,new_sp,xr=xr,yr=yr,psym=5,title='Pixon s (Mm)',ytitle=ytitle,$
	xtitle=xtitle
oplot,xr,yr
for ctr=0,24 do begin
	xyouts,old_sp[ctr],new_sp[ctr],strcompress(ctr+1,/rem)
endfor

ytitle = 'Projected Orientation (radians)'
xtitle = 'True Orientation (radians)'

plot,old_alpham,new_alpham,/iso,psym=5,title='MEM !7a!3',xtitle=xtitle,ytitle=ytitle
oplot,[0,2],[0,2]
for ctr=0,24 do begin
	xyouts,old_alpham[ctr],new_alpham[ctr],strcompress(ctr+1,/rem)
endfor

plot,old_alphac,new_alphac,/iso,psym=5,title='Clean !7a!3',xtitle=xtitle,ytitle=ytitle
oplot,[0,2],[0,2]
for ctr=0,24 do begin
	xyouts,old_alphac[ctr],new_alphac[ctr],strcompress(ctr+1,/rem)
endfor

plot,old_alphap,new_alphap,/iso,psym=5,title='Pixon !7a!3',xtitle=xtitle,ytitle=ytitle
oplot,[0,2],[0,2]
for ctr=0,24 do begin
	xyouts,old_alphap[ctr],new_alphap[ctr],strcompress(ctr+1,/rem)
endfor

!p.multi=0
device,/close
set_plot,'x'

return
end