OV-GSO-DC / dustem-wrapper_idl

Blame view

LabTools/IRAP/JPB/phangs_isrf_pipeline.pro 11.8 KB
Edit Raw Normal View History
2adbba75   Jean-Philippe Bernard   First commit 
1
2
 
PRO phangs_isrf_pipeline
eb7ea183   Jean-Philippe Bernard   included stellar ... 
3
 
;Note: faire un chmod a+rwx * apres transfert sur alma1
2adbba75   Jean-Philippe Bernard   First commit 
4
5
 

;PHANGS ISRF pipeline
9b37a060   Jean-Philippe Bernard   improved for smoo... 
6
7
 
save_data_dir=!phangs_data_dir+'/ISRF/WORK/'
str='ls '+save_data_dir+'ngc0628_*SPIRE3.sav'
f8a541d6   Jean-Philippe Bernard   updated to new fi... 
8
9
 
print,str
8c4defa5   Jean-Philippe Bernard   improved 
10
11
 
str='ls '+save_data_dir+'ngc0628_*images*.sav'
spawn,str
f8a541d6   Jean-Philippe Bernard   updated to new fi... 
12
 
str='ls '+save_data_dir+'ngc0628_*seds_muse_pixels*.sav'
8c4defa5   Jean-Philippe Bernard   improved 
13
 
spawn,str
2adbba75   Jean-Philippe Bernard   First commit 
14
e314ac19   Jean-Philippe Bernard   improved 
15
 
;==== make sky images
6fd06be5   Jean-Philippe Bernard   improved 
16
 
;produces _jwst_images.sav and _ref_header.sav
f8a541d6   Jean-Philippe Bernard   updated to new fi... 
17
 
;/data/projects/phangs/phangs_drive/Archive/PHANGS_JWST/DR1/ngc0628/f1000w_psf_matched
e314ac19   Jean-Philippe Bernard   improved 
18
 
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop
9b37a060   Jean-Philippe Bernard   improved for smoo... 
19
 
;produces _muse_filters_images.sav,_muse_data.sav, _muse_NH.sav
e314ac19   Jean-Philippe Bernard   improved 
20
 
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop
9b37a060   Jean-Philippe Bernard   improved for smoo... 
21
 
;produces _astrosat_images.sav
e314ac19   Jean-Philippe Bernard   improved 
22
 
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop
6fd06be5   Jean-Philippe Bernard   improved 
23
24
 
;produces _herschel_images.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop
8c4defa5   Jean-Philippe Bernard   improved 
25
 
;produces _HI_images.sav
14ae4f23   Jean-Philippe Bernard   improved 
26
 
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop
8c4defa5   Jean-Philippe Bernard   improved 
27
 
;produces _CO_images.sav
6fd06be5   Jean-Philippe Bernard   improved 
28
 
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop
e314ac19   Jean-Philippe Bernard   improved 
29
9b37a060   Jean-Philippe Bernard   improved for smoo... 
30
31
 
;produces _muse_data.sav,
phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop
3e3d714b   Jean-Philippe Bernard   removed hardcoded... 
32
e314ac19   Jean-Philippe Bernard   improved 
33
34
 
;==== extract SEDs on MUSE pixels
;produces _jwst_seds_muse_pixels.sav, _muse_seds_muse_pixels.sav', _astrosat_seds_muse_pixels.sav
9b37a060   Jean-Philippe Bernard   improved for smoo... 
35
 
;pointer for vid 40749, 41095, 41405, 41406, 41581 .... is not referenced in all_seds_indices
e314ac19   Jean-Philippe Bernard   improved 
36
 
phangs_extract_seds,source_name='ngc0628'
e314ac19   Jean-Philippe Bernard   improved 
37
6c2927ee   Jean-Philippe Bernard   improved 
38
 
;=== Compute ISRF prediction. This is done only at full resolution
de766928   Jean-Philippe Bernard   improved 
39
 
;=== produce _isrf_min_prediction.sav, _isrf_min_prediction.fits (product)
6c2927ee   Jean-Philippe Bernard   improved 
40
41
 
make_phangs_ssps_isrf_prediction,source_name='ngc0628',/save
de766928   Jean-Philippe Bernard   improved 
42
43
44
 
;==== Make isrf classes
;produces isrf_classes_one_ratio_on_*.fits, *_isrf_classes_map.fits, *_isrf_classes_voronoi.sav
make_phangs_isrf_classes,bidon,source_name='ngc0628',/save
e314ac19   Jean-Philippe Bernard   improved 
45
c4ea1c02   Jean-Philippe Bernard   improved 
46
47
48
49
50
 
;make_sed_phangs_tables,/test,/show_seds   ;test of plugin dustem_plugin_phangs_class_isrf.pro
;Generally very long. Do it using nohup (see make_phangs_grid.pro)
make_sed_phangs_tables,grid_type=2,isrf_class=15
make_sed_phangs_tables,grid_type=2,isrf_class=0
;...
de766928   Jean-Philippe Bernard   improved 
51
52
53
54
55
56
 

;=== Fit the SEDs with the ISRF Grids
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save

phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/from_restore
c4ea1c02   Jean-Philippe Bernard   improved 
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
 

;============================================================================================================
;===== NGC0628 at SPIRE3 resolution
;============================================================================================================
;resolution_filter='SPIRE3'
;ngc0628_jwst_images_SPIRE3.sav, ngc0628_ref_header_SPIRE3.sav
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_CO_images_SPIRE3.sav
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_HI_images_SPIRE3.sav
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_astrosat_data_SPIRE3.sav
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'

;ngc0628_muse_filters_data_SPIRE3.sav, ngc0628_muse_images_SPIRE3.sav (not produced _muse_NH_SPIRE3.sav)
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_herschel_images_SPIRE3.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'

phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'

phangs_extract_seds,source_name='ngc0628',resolution_filter='SPIRE3'
8cfdf0b6   Jean-Philippe Bernard   improved 
80
 
;ngc0628_isrf_min_prediction_SPIRE3.sav
eb7ea183   Jean-Philippe Bernard   included stellar ... 
81
 
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='SPIRE3',/save,/show,/nostop
c4ea1c02   Jean-Philippe Bernard   improved 
82
c4ea1c02   Jean-Philippe Bernard   improved 
83
84
 
make_phangs_isrf_classes,bidon,source_name='ngc0628',/save,resolution_filter='SPIRE3'
8cfdf0b6   Jean-Philippe Bernard   improved 
85
86
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save,resolution_filter='SPIRE3'
c7d22c08   Jean-Philippe Bernard   improved 
87
88
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save,resolution_filter='SPIRE3',/from_restore
966e3740   Jean-Philippe Bernard   improved 
89
90
91
92
93
94
95
 

;============================================================================================================
;===== NGC0628 at PACS3 resolution
;============================================================================================================
;resolution_filter='PACS3'
;ngc0628_jwst_images_PACS3.sav, ngc0628_ref_header_PACS3.sav
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
c7d22c08   Jean-Philippe Bernard   improved 
96
97
 
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
966e3740   Jean-Philippe Bernard   improved 
98
 
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
c7d22c08   Jean-Philippe Bernard   improved 
99
100
101
102
103
104
105
 
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'

phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'

phangs_extract_seds,source_name='ngc0628',resolution_filter='PACS3'
eb7ea183   Jean-Philippe Bernard   included stellar ... 
106
 
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='PACS3',/save,/show,/nostop
c7d22c08   Jean-Philippe Bernard   improved 
107
108
109
110
 

make_phangs_isrf_classes,bidon,source_name='ngc0628',/save,resolution_filter='PACS3'

phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save,resolution_filter='PACS3',/show_seds
966e3740   Jean-Philippe Bernard   improved 
111
524e900a   Jean-Philippe Bernard   fixed smoothing I... 
112
113
114
115
116
117
118
119
120
121
 
;============================================================================================================
;===== NGC0628 at PACS2 resolution
;============================================================================================================

phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'
c4ea1c02   Jean-Philippe Bernard   improved 
122
524e900a   Jean-Philippe Bernard   fixed smoothing I... 
123
124
125
126
127
128
 
phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS2'

phangs_extract_seds,source_name='ngc0628',resolution_filter='PACS2'
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='PACS2',/save

make_phangs_isrf_classes,bidon,source_name='ngc0628',/save,resolution_filter='PACS2'
eb7ea183   Jean-Philippe Bernard   included stellar ... 
129
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save,resolution_filter='PACS2',/subtract_star_light
524e900a   Jean-Philippe Bernard   fixed smoothing I... 
130
131
 

;===== j'en suis la sur mon mac
c4ea1c02   Jean-Philippe Bernard   improved 
132
133
134
eb7ea183   Jean-Philippe Bernard   included stellar ... 
135
136
137
138
139
140
524e900a   Jean-Philippe Bernard   fixed smoothing I... 
141
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/from_restore,resolution_filter='PACS3'
c4ea1c02   Jean-Philippe Bernard   improved 
142
dc862aaf   Jean-Philippe Bernard   improved 
143
144
de766928   Jean-Philippe Bernard   improved 
145
 
plot_phangs_ssps_against_uv_fast,bidon
e314ac19   Jean-Philippe Bernard   improved 
146
ccfa6760   Jean-Philippe Bernard   corrected ISRF cl... 
147
 
;==== Check the classes and do a plot
de766928   Jean-Philippe Bernard   improved 
148
149
 
;==== produce ngc0682_isrf_classes_one_ratio.sav
;==== requires : ngc0628_astrosat_voronoi_prediction_fast.sav
ccfa6760   Jean-Philippe Bernard   corrected ISRF cl... 
150
 
check_phangs_ssps_isrf_prediction,/show_map
e314ac19   Jean-Philippe Bernard   improved 
151
de766928   Jean-Philippe Bernard   improved 
152
153
154
 
;fit_phangs_ngc0628_nir_continuum,/from_restore,/show_plots     ;fit seds of some voronoi bins
;fit_phangs_ngc0628_nir_continuum,/from_restore,/brute_force
;fit_phangs_ngc0628_nir_continuum,/from_restore,/grid_brute_force
3b79f061   Jean-Philippe Bernard   improved 
155
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
156
 
;=== Compute ISRF prediction. This is done only at full resolution
2adbba75   Jean-Philippe Bernard   First commit 
157
 
make_phangs_ssps_isrf_prediction,source_name='ngc0628',/save
3b79f061   Jean-Philippe Bernard   improved 
158
159
 
check_phangs_ssps_isrf_prediction
ccfa6760   Jean-Philippe Bernard   corrected ISRF cl... 
160
 
;==== Make isrf classes
e2fb44a5   Jean-Philippe Bernard   improved 
161
 
make_phangs_isrf_classes,bidon,source_name='ngc0628',/fits
ccfa6760   Jean-Philippe Bernard   corrected ISRF cl... 
162
e314ac19   Jean-Philippe Bernard   improved 
163
ccfa6760   Jean-Philippe Bernard   corrected ISRF cl... 
164
165
166
167
168
169
 


;=== had to convert previous grids like this:
convert_all_grid_tables_from_v1_to_v2

;=== Fit the SEDs with the ISRF Grids
442151d0   Jean-Philippe Bernard   improved 
170
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save
2adbba75   Jean-Philippe Bernard   First commit 
171
442151d0   Jean-Philippe Bernard   improved 
172
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/from_restore
82beee16   Jean-Philippe Bernard   improved 
173
ee1fd9a9   Jean-Philippe Bernard   improved in the f... 
174
175
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
 
;============================================================================================================
;===== NGC0628 at SPIRE3 resolution
;============================================================================================================
;resolution_filter='SPIRE3'
;ngc0628_jwst_images_SPIRE3.sav, ngc0628_ref_header_SPIRE3.sav
phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_CO_images_SPIRE3.sav
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_HI_images_SPIRE3.sav
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_astrosat_data_SPIRE3.sav
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_muse_filters_data_SPIRE3.sav, ngc0628_muse_images_SPIRE3.sav (not produced _muse_NH_SPIRE3.sav)
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
;ngc0628_herschel_images_SPIRE3.sav
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
9b37a060   Jean-Philippe Bernard   improved for smoo... 
192
193
194
 

phangs_make_muse_filters_data,source_name='ngc0628',/save,/show,/nostop,resolution_filter='SPIRE3'
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
195
196
 
;ngc0628_jwst_seds_muse_pixels_SPIRE3.sav,ngc0628_seds_indices_SPIRE3.sav,ngc0628_muse_seds_muse_pixels_SPIRE3.sav
;ngc0628_astrosat_seds_muse_pixels_SPIRE3.sav,ngc0628_herschel_seds_muse_pixels_SPIRE3.sav,ngc0628_all_seds_muse_pixels_SPIRE3.sav
9b37a060   Jean-Philippe Bernard   improved for smoo... 
197
 
phangs_extract_seds,source_name='ngc0628',resolution_filter='SPIRE3',/reset_seds
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
 
;ngc0628_isrf_min_prediction_SPIRE3.sav
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='SPIRE3',/save
;isrf_classes_one_ratio_on_ngc0628_SPIRE3.fits,ngc0628_isrf_classes_map_SPIRE3.fits,ngc0628_isrf_classes_voronoi_SPIRE3.sav
make_phangs_isrf_classes,bidon,source_name='ngc0628',/fits,resolution_filter='SPIRE3'
;J'en SUIS LA
phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/include_herschel,/normalize,/fit_G0,/save,resolution_filter='SPIRE3'
;J'en SUIS LA






phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_co_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'
phangs_make_herschel_images,source_name='ngc0628',/save,/show,/nostop,resolution_filter='PACS3'

phangs_extract_seds,source_name='ngc0628',resolution_filter='PACS3'

;=== smooth MUSE isrfs to a given resolution
;phangs_smooth_muse_isrf,'ngc0628',resolution_filter='PACS3',/test
phangs_smooth_muse_isrf,'ngc0628',resolution_filter='PACS3',/save

;==== Make isrf classes
make_phangs_isrf_classes,bidon,source_name='ngc0628',/fits,resolution_filter='PACS3'
ee1fd9a9   Jean-Philippe Bernard   improved in the f... 
227
ee1fd9a9   Jean-Philippe Bernard   improved in the f... 
228
 
phangs_smooth_muse_isrf,'ngc0628',reso_filter='SPIRE1',/save
33d0e4eb   Jean-Philippe Bernard   improved 
229
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
230
231
232
233
82beee16   Jean-Philippe Bernard   improved 
234
235
236
237
 
;=== NGC3351
;Note: faire un chmod a+rwx * apres transfert sur alma1

source_name='ngc3351'
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
238
239
 
;J'en SUIS LA
82beee16   Jean-Philippe Bernard   improved 
240
 
;produces _jwst_images.sav and _ref_header.sav
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
241
 
phangs_make_jwst_images,source_name='ngc3351',/save,/show,/nostop
82beee16   Jean-Philippe Bernard   improved 
242
 
;produces _muse_filters_data.sav,_muse_images.sav, _muse_NH.sav
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
243
 
phangs_make_muse_filters_images,source_name='ngc3351',/save,/show,/nostop
82beee16   Jean-Philippe Bernard   improved 
244
 
;produces _astrosat_data.sav
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
245
 
phangs_make_astrosat_images,source_name='ngc3351',/save,/show,/nostop
82beee16   Jean-Philippe Bernard   improved 
246
 
;produces _herschel_images.sav
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
247
248
249
 
phangs_make_herschel_images,source_name='ngc3351',/save,/show,/nostop
phangs_make_hi_images,source_name='ngc3351',/save,/show,/nostop
phangs_make_co_images,source_name='ngc3351',/save,/show,/nostop
82beee16   Jean-Philippe Bernard   improved 
250
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
251
 
phangs_extract_seds,source_name='ngc3351',/save
82beee16   Jean-Philippe Bernard   improved 
252
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
253
 
make_phangs_ssps_isrf_prediction,source_name='ngc3351',/save
82beee16   Jean-Philippe Bernard   improved 
254
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
255
 
make_phangs_isrf_classes,source_name='ngc3351',/fits
82beee16   Jean-Philippe Bernard   improved 
256
9b0b6d7e   Jean-Philippe Bernard   finished implemen... 
257
 
phangs_brute_force_fit_with_isrf_grid,source_name='ngc3351',/include_herschel,/normalize,/fit_G0
f522d335   Jean-Philippe Bernard   improved 
258
2adbba75   Jean-Philippe Bernard   First commit 
259
 
END