Date - 2004.03.09 Tester - J.E. Hibbard (CV) Platform: Merger Version: Daily AIPS++ Version 1.9 Build 473 Files used to test imageprofilefitter fitting results (see Defect AOCso4745; see also 2004-03-11-hibbard.txt for getting coordsystem right on AIPS datacube). ------------------------------------------------- pro fitgau, infile, firstguess ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; program to fit gaussians to data in infile using IDL fitting ;; routines GAUSSFIT (reports no errors) and CURVEFIT ;; infile=datafile (x,y,*,*,*,mask) ;; firstguess=array of first guesses for [Io,Xmean,FWHM] ;; JEH 3/2004 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; GAUSSFIT and gaussian take sigma, not FWHM, as input: firstguess[2]=firstguess[2]/2.355 ;; readcol,infile,x,data,est,fit,resid,mask, f='F,F,F,F,F,I' yfit1=GAUSSFIT(x,data,answer1,estimates=firstguess,nterms=3) answer2=firstguess weights=1/data yfit2=CURVEFIT(x,data,weights,answer2,sig2,FUNCTION_NAME='gaussian') ;; ;; Convert fluxes to mJy, widths back to FWHM: firstguess[0]=firstguess[0]*1000. firstguess[2]=firstguess[2]*2.355 answer1[0]=answer1[0]*1000. answer1[2]=answer1[2]*2.355 answer2[0]=answer2[0]*1000. answer2[2]=answer2[2]*2.355 sig2[0]=sig2[0]*1000. sig2[2]=sig2[2]*2.355 print,'results for file ',infile,':' print,";; Peak Mean FWHM" print,";;ESTIMATE: ",STRING(firstguess,FORMAT='(f9.2,2f15.2)') print,";;GAUSSFIT: ",STRING(answer1,FORMAT='(f9.2,2f15.2)') print,";;CURVEFIT: ",STRING(answer2[0],sig2[0], $ answer2[1],sig2[1],answer2[2],sig2[2], $ FORMAT='(F9.2,"(",F5.2,")",F9.2,"(",F5.2,")",F9.2,"(",F5.2,")")') ; FORMAT='(F9.2,"(",F4.2,")",F9.2,"(",F4.2,")",F9.2,"(",F4.2,")")') ; loadct,30 plot,x,data oplot,x,yfit1,color=150,thick=2,linestyle=3 oplot,x,yfit2,color=150,thick=3,linestyle=2 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro gaussian,x,params,f,pder ;; function to compute gaussian at x, given params=[peak,mean,sigma] sigma=params[2] num=(x-params[1])/sigma f=(params[0]*exp(-num*num/2.)) IF N_PARAMS() GE 4 THEN $ pder = [ [f/params[0]], [f*num/sigma], [f*num*num/sigma] ] end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;HERE ARE MY NOTES ON RESULTS OF FITS IN IDL/AIPS++/AIPS: JEH 3/9/04 ;;In AIPS: take LMV in, trans 312, slice, ;;blc 1 238 297; trc 83 238 297;go slice; ;;blc 1 242 282; trc 83 242 282;go slice; ;;blc 1 245 275; trc 83 245 275;go slice; ;;blc 1 249 263; trc 83 249 263;go slice; ;;blc 1 255 249; trc 83 255 249;go slice; ;;task 'slfit';ngauss 1;dopos 1;domax 1;dowid 1; ;;inver 1;tvsl;gpos 870;gmax 0.018;gwid 11.6;go slfit ;;merger> SLFIT1: Results in SLice units ;;merger> SLFIT1: Cmp Parameter Value Error ;;merger> SLFIT1: 1 Peak 1.953942E+01 2.2886E-01 MilliJY/BEAM ;;merger> SLFIT1: 1 Position 8.797791E+02 1.0949E-01 Kilo FELO-HEL ;;merger> SLFIT1: 1 Fullwhm 1.906065E+01 2.5792E-01 Kilo FELO-HEL ;;inver 2;tvsl;gpos 850;gmax 0.0184;gwid 11.0;go slfit ;;merger> SLFIT1: Cmp Parameter Value Error ;;merger> SLFIT1: 1 Peak 1.737267E+01 2.3312E-01 MilliJY/BEAM ;;merger> SLFIT1: 1 Position 8.489680E+02 1.6514E-01 Kilo FELO-HEL ;;merger> SLFIT1: 1 Fullwhm 2.508781E+01 3.8882E-01 Kilo FELO-HEL ;;inver 3;tvsl;gpos 840;gmax 0.009;gwid 15.5;go slfit ;;merger> SLFIT1: Cmp Parameter Value Error ;;merger> SLFIT1: 1 Peak 1.107632E+01 1.9270E-01 MilliJY/BEAM ;;merger> SLFIT1: 1 Position 8.376895E+02 2.6514E-01 Kilo FELO-HEL ;;merger> SLFIT1: 1 Fullwhm 3.105991E+01 6.2429E-01 Kilo FELO-HEL ;;inver 4;tvsl;gpos 810;gmax 0.012;gwid 14;go slfit ;;merger> SLFIT1: Cmp Parameter Value Error ;;merger> SLFIT1: 1 Peak 1.269988E+01 2.2969E-01 MilliJY/BEAM ;;merger> SLFIT1: 1 Position 8.140298E+02 2.2264E-01 Kilo FELO-HEL ;;merger> SLFIT1: 1 Fullwhm 2.508513E+01 5.2422E-01 Kilo FELO-HEL ;;inver 5;tvsl;gpos 785;gmax 0.010;gwid 11;go slfit ;;merger> SLFIT1: Cmp Parameter Value Error ;;merger> SLFIT1: 1 Peak 1.040294E+01 2.2683E-01 MilliJY/BEAM ;;merger> SLFIT1: 1 Position 7.894862E+02 2.2441E-01 Kilo FELO-HEL ;;merger> SLFIT1: 1 Fullwhm 2.095891E+01 5.2840E-01 Kilo FELO-HEL ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; RESULTS FROM DRAWSPEC (curtosey of Harvey Liszt, 3/10/04): ;; ;;1 GAUSSIANS Fitting Source=*238 297, Scan Number=20718, Line=H I ;;CHI**2 = 37.15 with 22 deg/freedom; channel spacing = -2.590077637 ;;RMS of fit = 0.0013 over 25 chans: RMS of spectrum = 0.0010717 ;; ;; ---CENTER---------PEAK----------FWHM----------AREA---- ;; 1 879.79164899 0.0195364 19.072511573 0.3966304976 ;;± 0.2669751121 0.0005578 0.6293347580 0.0111329938 ;; ;;1 GAUSSIANS Fitting Source=*242 282, Scan Number=13042, Line=H I ;;CHI**2 = 25.42 with 26 deg/freedom; channel spacing = -2.590077637 ;;RMS of fit = 0.0012 over 29 chans: RMS of spectrum = 0.0012704 ;; ;; ---CENTER---------PEAK----------FWHM----------AREA---- ;; 1 848.97646412 0.0173293 25.200970019 0.4648706149 ;;± 0.4100559187 0.0005750 0.9657002943 0.0151590829 ;; ;;1 GAUSSIANS Fitting Source=*245 275, Scan Number=9461, Line=H I ;;CHI**2 = 21.39 with 23 deg/freedom; channel spacing = -2.590077637 ;;RMS of fit = 0.0012 over 26 chans: RMS of spectrum = 0.0013095 ;; ;; ---CENTER---------PEAK----------FWHM----------AREA---- ;; 1 837.67857091 0.0110968 30.917417925 0.3652024664 ;;± 0.7327359152 0.0005371 1.7477876165 0.0174869142 ;; ;;1 GAUSSIANS Fitting Source=*249 263, Scan Number=3321, Line=H I ;;CHI**2 = 67.53 with 23 deg/freedom; channel spacing = -2.590077637 ;;RMS of fit = 0.0015 over 26 chans: RMS of spectrum = 0.0009575 ;; ;; ---CENTER---------PEAK----------FWHM----------AREA---- ;; 1 814.03489483 0.0126808 25.168825215 0.3397362263 ;;± 0.4221052939 0.0004338 0.9949496191 0.0114251740 ;; ;;1 GAUSSIANS Fitting Source=*255 249, Scan Number=-3841, Line=H I ;;CHI**2 = 39.32 with 23 deg/freedom; channel spacing = -2.590077637 ;;RMS of fit = 0.0014 over 26 chans: RMS of spectrum = 0.0011161 ;; ;; ---CENTER---------PEAK----------FWHM----------AREA---- ;; 1 789.46824755 0.0104205 20.893526831 0.2317574311 ;;± 0.5455047501 0.0005548 1.2849445359 0.0121309495 ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;RESULTS FOR EVERYTHING: ;;IDL> fitgau,'fit1.txt',[0.018,794.8,23.1] ;; Peak Mean FWHM ;;ESTIMATE: 18.00 794.80 23.10 ;;aips++: 19.54(0.61) 793.69(0.29) 18.97(0.68) ;;DRAWSPEC: 19.54(0.56) 879.79(0.27) 19.07(0.63) ;;IDL GAUSSFIT: 19.54 793.63 18.97 ;;AIPS SLFIT: 19.54(0.23) 879.78(0.11) 19.06(0.26) ;; ;;IDL> fitgau,'fit2.txt',[0.0184,762.,22.05] ;; Peak Mean FWHM ;;ESTIMATE: 18.40 762.00 22.05 ;;aips++: 17.37(0.61) 762.98(0.43) 24.96(1.02) ;;DRAWSPEC: 17.33(0.57) 848.98(0.41) 25.20(0.97) ;;IDL GAUSSFIT: 17.38 762.98 24.96 ;;AIPS SLFIT: 17.37(0.23) 848.97(0.17) 25.09(0.39) ;; ;;IDL> fitgau,'fit3.txt',[0.009,751.8,30.9] ;; Peak Mean FWHM ;;ESTIMATE: 9.00 751.80 30.90 ;;aips++: 11.08(0.52) 751.76(0.71) 30.91(1.67) ;;DRAWSPEC: 11.10(0.54) 837.68(0.73) 30.92(1.75) ;;IDL GAUSSFIT: 11.08 751.76 30.91 ;;AIPS SLFIT: 11.08(0.19) 837.69(0.27) 31.06(0.62) ;; ;;IDL> fitgau,'fit4.txt',[0.012,728.4,21.9] ;; Peak Mean FWHM ;;ESTIMATE: 12.00 728.40 21.90 ;;aips++: 12.70(0.60) 728.20(0.58) 24.96(1.36) ;;DRAWSPEC: 12.68(0.43) 814.03(0.42) 25.17(0.99) ;;IDL GAUSSFIT: 12.71 728.21 24.92 ;;AIPS SLFIT: 12.70(0.23) 814.03(0.22) 25.09(0.52) ;; ;;IDL> fitgau,'fit5.txt',[0.0099,702.7,21.9] ;; Peak Mean FWHM ;;ESTIMATE: 9.90 702.70 21.90 ;;aips++: 10.40(0.63) 703.80(0.61) 20.86(1.45) ;;DRAWSPEC: 10.42(0.55) 789.47(0.55) 20.89(1.28) ;;IDL GAUSSFIT: 10.38 703.80 20.96 ;;AIPS SLFIT: 10.40(0.23) 789.49(0.22) 20.96(0.53) ;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;