diff --git a/LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro b/LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro
index 0952c6e..8a7a2c9 100644
--- a/LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro
+++ b/LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro
@@ -109,7 +109,11 @@ chi2s=fltarr(Nvor)+la_undef()
rchi2s=fltarr(Nvor)+la_undef()
dG0s=fltarr(Nvor)+la_undef()
dYpahs=fltarr(Nvor)+la_undef()
-dYVSGs=fltarr(Nvor)+la_undef()
+dYVSGs=fltarr(Nvor)+la_undef()
+G0_hit=fltarr(Nvor)+la_undef()
+Ypah_hit=fltarr(Nvor)+la_undef()
+Yvsg_hit=fltarr(Nvor)+la_undef()
+
use_NHmap=NHCO ;used NH map in 1.e21 (from CO)
;use_NHmap=NH_map ;used NH map in 1.e21 (from MUSE)
@@ -158,10 +162,10 @@ FOR isrf_class=class_min,class_max DO BEGIN
;=== Get the G0 value for that vid
vid=ind_class[vvid]
;print,vid
- IF vid EQ 12641 THEN BEGIN
- message,'stopped at vid =12641',/continue
- stop
- ENDIF
+ ;IF vid EQ 12641 THEN BEGIN
+ ; message,'stopped at vid =12641',/continue
+ ; stop
+ ;ENDIF
G0=interpol(ISRFs[*,vid],isrf_wavelengths,1.0)/Mathis_1mic
fixed_parameters_values=[G0]
sed=*seds_ptr[vid]
@@ -183,11 +187,16 @@ FOR isrf_class=class_min,class_max DO BEGIN
;Ypahs[vid]=params[1]/fact ;extenssive quantities must be divided by normalization factor
;Yvsgs[vid]=params[2]/fact
G0s[vid]=G0
- Ypahs[vid]=params[0]/fact ;extenssive quantities must be divided by normalization factor
- Yvsgs[vid]=params[1]/fact
+ ;Ypahs[vid]=params[1]/fact ;extenssive quantities must be divided by normalization factor
+ ;Yvsgs[vid]=params[2]/fact
+ Ypahs[vid]=weighted_params[1]/fact ;extenssive quantities must be divided by normalization factor
+ Yvsgs[vid]=weighted_params[2]/fact
facts[vid]=fact
chi2s[vid]=chi2
rchi2s[vid]=rchi2
+ G0_hit[vid]=params_hit[0]
+ Ypah_hit[vid]=params_hit[1]
+ Yvsg_hit[vid]=params_hit[2]
;This would be if there was no fixed_parameters
;dG0s[vid]=params_uncertainties[0]
;dYpahs[vid]=params_uncertainties[1]/fact ;extenssive quantities must be devided by normalization factor
@@ -225,6 +234,9 @@ dG0_map=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'))+la_undef()
dYpah_map=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'))+la_undef()
dYVSG_map=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'))+la_undef()
fact_map=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'))+la_undef()
+Ypah_hit_map=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'))+la_undef()
+YVSG_hit_map=fltarr(sxpar(href,'NAXIS1'),sxpar(href,'NAXIS2'))+la_undef()
+
FOR vid=0LL,Nvor-1 DO BEGIN
IF vid mod 100 EQ 0 THEN print,(1.*vid)/Nvor*100.
;ind=where(voronoi_id EQ vid)
@@ -237,14 +249,10 @@ FOR vid=0LL,Nvor-1 DO BEGIN
chi2_map[*all_seds_indices[vid]]=chi2s[vid]
rchi2_map[*all_seds_indices[vid]]=rchi2s[vid]
fact_map[*all_seds_indices[vid]]=facts[vid]
+ Ypah_hit_map[*all_seds_indices[vid]]=Ypah_hit[vid]
+ YVSG_hit_map[*all_seds_indices[vid]]=Yvsg_hit[vid]
ENDFOR
-;FOR vid=0LL,Nvor-1 DO GO_map[*all_seds_indices[vid]]=GOs[vid]
-;FOR vid=0LL,Nvor-1 DO Ypah_map[*all_seds_indices[vid]]=Ypahs[vid]
-;FOR vid=0LL,Nvor-1 DO YVSG_map[*all_seds_indices[vid]]=Yvsgs[vid]
-;FOR vid=0LL,Nvor-1 DO chi2_map[*all_seds_indices[vid]]=chi2s[vid]
-;stop
-
;G0 histogram
win=0L
window,win,xsize=800,ysize=800 & win=win+1
@@ -269,7 +277,7 @@ imrange=[-2.,6]
image_cont20,la_log10(Ypah_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Ypah)'
window,win,xsize=800,ysize=900 & win=win+1
-imrange=[-2.,2]
+imrange=[-2.,6]
;imrange=[-3.,0.]-3. ;why ??
;image_cont20,la_log10(YVSG_map),Href,/square,imrange=[-2.,6],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Yvsg)'
image_cont20,la_log10(YVSG_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Yvsg)'
@@ -280,9 +288,17 @@ imrange=[3,7]
image_cont20,la_log10(chi2_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='log10(Chi2)',off_bar=obp
window,win,xsize=800,ysize=900 & win=win+1
-imrange=[-1,2]
+imrange=[-3,2]
image_cont20,la_log10(fact_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='log10(fact)',off_bar=obp
+window,win,xsize=800,ysize=900 & win=win+1
+imrange=[-2,2]
+image_cont20,Ypah_hit_map,Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='Ypah hit',off_bar=obp
+
+window,win,xsize=800,ysize=900 & win=win+1
+imrange=[-2,2]
+image_cont20,Yvsg_hit_map,Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='Yvsg hit',off_bar=obp
+
stop
jwst_4_coutours=fltarr((size(jwst_images))[1],(size(jwst_images))[2],2)
jwst_4_coutours[*,*,0]=jwst_images[*,*,0,5]
--
libgit2 0.21.2