;$Id: ppe.pro,v 1.9 2021/11/11 19:42:05 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  ;!p.charthick=1.5 & !p.thick=1.5 & !x.thick=1.5 & !y.thick=1.5
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
siz=1.4
!p.charsize=1.4
!x.margin=[7.3,0.4]
!y.margin=[3.1,0.4]
!x.title='!13E!6!s!dM!n!r!umax!n'
!y.title='!13E!6!s!dGW!n!r!usat!n'
print,"$mv idl.ps ~/tex/tina/BBN/fig/ppe.eps"
;
file='pe.txt'
file='pe_5em5.txt'  ;(kf600_red)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
kf=600.
xr=[1e-5,2e-3]*600
yr=[1e-9,2e-4]
;
xr=[1e-2,2e+0]
yr=[1e-9,2e-0]
;
plot_oo,EEM,EEGW,ps=-1,xr=xr,yr=yr,/nodata
oplot,EEM,EEGW,ps=-1,col=122
;
;  overplot slopes
;
;xx=[2e-5,1e-3]*600 & oplot,xx,11e3*(xx/600.)^2.70,li=3
;xx=[.05,.20] & oplot,xx,1e-3*xx^2,li=3
;xyouts,.07,1.7e-5,siz=siz,'!9A!6(!13E!6!s!dM!n!r!umax!n)!u2!n'
;xyouts,.12,6.6e-7,siz=siz,'!9A!6(!13E!6!s!dM!n!r!umax!n)!u2.7!n'
xyouts,.27,5.0e-6,siz=siz,'!7m!6=5!9X!610!u-5!n',col=122
;
file='pe_ramp0_5em5.txt' ;(kf600_red_dotted)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM,EEGW,ps=-1,li=1,col=122
;xyouts,.053,9.0e-7,siz=siz,'!7m!6=10!u-4!n',col=55
;
file='pe_ramp0.txt'
file='pe_2em7.txt'
file='pe_2em6.txt'
file='pe_f005.txt'
file='pe_ramp0_5em5.txt'
file='pe_1em4.txt' ;(kf600_orange)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM,EEGW,ps=-1,col=155
xyouts,.070,6.0e-6,siz=siz,'!7m!6=10!u-4!n',col=155
;
file='pe_2em6.txt' ;(kf600_blue)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM,EEGW,ps=-1,col=55
xyouts,.018,1.0e-7,siz=siz,'!7m!6=2!9X!610!u-6!n',col=55
;
file='pe_2em7.txt' ;(kf600_black)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM,EEGW,ps=-1
xyouts,.024,7.0e-9,siz=siz,'!7m!6=2!9X!610!u-7!n'
xyouts,.11,1.2e-6,siz=.8*siz,'!8k!6!df!n=600'
;
;-----------------------------------------------------------------------------
circ_sym,1.3,1
file='pe_2em2_kf6.txt'  ;(kf6_red)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM,EEGW,ps=-8,col=122
xyouts,.45,1.8e-1,siz=siz,'!7m!6=2!9X!610!u-2!n',col=122
;xx=[.45,1.98] & oplot,xx,10e-2*xx^2.2,col=155
;
circ_sym,1.3,1
file='pe_5em3_kf6.txt'  ;(kf6_blue)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM,EEGW,ps=8,col=55
xx=[.105,0.98] & oplot,xx,6e-2*xx^2.1,col=55
;xyouts,.13,6.6e-3,siz=siz,'!9A!6(!13E!6!s!dM!n!r!umax!n)!u2.1!n',col=55
xyouts,.11,4.0e-3,siz=siz,'!7m!6=5!9X!610!u-3!n',col=55
xyouts,.28,3.6e-2,siz=.8*siz,'!8k!6!df!n=6'
;
circ_sym,1.3,0
file='pe_line.txt'
file='pe_line_kf6.txt'  ;(not used)
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
;oplot,EEM,EEGW,ps=8,col=155
;xx=[.05,1.90] & oplot,xx,1e-1*xx^2.7,li=3
;
END