;$Id: ppe_scl.pro,v 1.2 2021/09/21 02:59:42 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/!8k!6!df!n'
!y.title='!13E!6!s!dGW!n!r!usat!n'
print,"$mv idl.ps ~/tex/tina/BBN/fig/ppe_scl.eps"
;
file='pe.txt'
file='pe_5em5.txt'
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,6e-1]
yr=[1e-9,2e-0]
;
plot_oo,EEM/kf,EEGW,ps=-1,xr=xr,yr=yr,/nodata
oplot,EEM/kf,EEGW,ps=-1,col=122
xyouts,3e-4,5.0e-5,siz=siz,'!7m!6=5!9X!610!u-5!n',col=122
;
file='pe_ramp0_5em5.txt'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM/kf,EEGW,ps=-1,li=1,col=122
;xyouts,4e-5,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'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM/kf,EEGW,ps=-1,col=155
xyouts,5e-5,6.0e-6,siz=siz,'!7m!6=10!u-4!n',col=155
;
file='pe_2em6.txt'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM/kf,EEGW,ps=-1,col=55
xyouts,1.3e-5,1.0e-7,siz=siz,'!7m!6=2!9X!610!u-6!n',col=55
;
file='pe_2em7.txt'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM/kf,EEGW,ps=-1
xyouts,4.5e-5,7.0e-9,siz=siz*.9,'!7m!6=2!9X!610!u-7!n'
xyouts,2e-4,1.2e-6,siz=.8*siz,'!8k!6!df!n=600'
;
;-----------------------------------------------------------------------------
kf=6.
circ_sym,1.3,1
file='pe_2em2_kf6.txt'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM/kf,EEGW,ps=-8,col=122
xyouts,3e-2,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'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
oplot,EEM/kf,EEGW,ps=8,col=55
xyouts,5e-2,2.0e-3,siz=siz,'!7m!6=5!9X!610!u-3!n',col=55
;
circ_sym,1.3,0
file='pe_line.txt'
file='pe_line_kf6.txt'
a=rtable(file,7,head=1)
f0=reform(a[0,*])
nu=reform(a[1,*])
EEM=reform(a[2,*])
EEGW=reform(a[3,*])
xx=[5e-3,.40] & oplot,xx,2*xx^2,li=3
xyouts,2e-2,6e-3,siz=.8*siz,'!8k!6!df!n=6'
xyouts,6e-3,1.6e-5,siz=siz,'!9A!6(!13E!6!s!dM!n!r!umax!n)!u2!n'
;
END