plot_phangs_ssps,against_uv_fast.pro
1.35 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
PRO plot_phangs_ssps_against_uv_fast
source_name='ngc0628'
data_dir='/Volumes/PILOT_FLIGHT1/PHANGS-JWST/DR1/'
restore,data_dir+'ngc0628_astrosat_minus_prediction_fast.sav',/verb
restore,data_dir+'ngc0628_astrosat_voronoi_prediction_fast.sav',/verb
;stop
window,0,xsize=800,ysize=800
!p.multi=[0,1,2]
xtit='(data-model) [MJy/sr]'
IF do_normalize THEN xtit='(data-model) normalized to 1e21 H/cm2 [MJy/sr]'
NHthreshold=0. & tit=''
;NHthreshold=0.1 & tit='NH > 1e20 H/cm2'
;NHthreshold=0.5 & tit='NH > 5e20 H/cm2'
;NHthreshold=1. & tit='NH > 1e21 H/cm2'
;NHthreshold=5. & tit='NH > 5e21 H/cm2'
;NHthreshold=10. & tit='NH > 1e22 H/cm2'
;NHthreshold=50. & tit='NH > 5e22 H/cm2'
;vmin=-0.1 & vmax=0.1 & Nbins=50
vmin=-1. & vmax=1. & Nbins=50
;vmin=-5. & vmax=5. & Nbins=100
;vmin=-10. & vmax=10. & Nbins=100
;im=diff_maps[*,*,0]
ims=actual_seds-predicted_seds
im=ims[*,k]
ind=where(im NE la_undef())
xtit='(data-model) normalized to 1e21 H/cm2 [MJy/sr]',
res=histogram(im[ind],locations=x,min=vmin,max=vmax,Nbins=Nbins)
yrange=[max([min(res),1.]),1.e4]
colors=['violet','blue','yellow','brown','red']
cgplot,x,res,/ylog,xtit=xtit,ytit='N',psym=10,yrange=yrange,/ysty,color='violet',tit=tit,/no_data
FOR k=0L,4 DO BEGIN
im=ims[*,k]
ind=where(im NE la_undef())
res=histogram(im[ind],locations=x,min=vmin,max=vmax,Nbins=Nbins)
cgoplot,x1,res1,color=colors[k],psym=10
ENDFOR
stop
END