Run haro29 GALEX images through to get SB HXZ: 02/25/10 Get exposure times: (Check exp FITS images) GI1_047067_UGCA281-fd-exp.fits 1647.62 GI1_047067_UGCA281-nd-exp.fits 1538.7 ======> fuv: 1647.62 s; nuv: 1538.7 s 1)Cut out smaller images: cutout subimage size: 800*800 centered around the galaxy center of cutout: imcopy GI1_047067_UGCA281-fd-int.fits[1516:2315,806:1605] haro29fcut.fit imcopy GI1_047067_UGCA281-nd-int.fits[1516:2315,806:1605] haro29ncut.fit -------------------------------------------- 2)Remove foreground stars and background galaxies: #Use the log files used by Deidre imedit haro29fcut.fit haro29fcl.fit cursor=haro29fcut.log imedit haro29ncut.fit haro29ncl.fit cursor=haro29ncut.log ========> haro29ncl.fit; haro29fcl.fit ---------------------------------------------------------- 3)Get the background now that the extraneous stuff has been removed: #Note: #The low counts rates in GALEX images leads to poisson, rather than gaussian noise. #which leads to the failure of using the usual sky-fitting routines in IRAF #Fit the sky of the original images usually failed. #So I selected many sky regions around the galaxy to get a mean sky level and the associated variance # #Here I first average the original images with 10*10 bin size ======= blkavg haro29fcl.fit haro29fblkav.fit 10 10 blkavg haro29ncl.fit haro29nblkav.fit 10 10 ========> haro29fblkav.fit; haro29nblkav.fit #Then use ELLIPSE to get mean pixel intensity and the associated rms in the elliptical annular sky regions around the galaxy. #I fixed the ELLIPSE parameters to that used to do surface photometry. # PACKAGE = isophote TASK = geompar (x0 = 40.536634) initial isophote center X (y0 = 40.066047) initial isophote center Y (ellip0 = 0.42) initial ellipticity (pa0 = 85.2) initial position angle (degrees) (sma0 = 5.) initial semi-major axis lenght (minsma = 1.) minimum semi-major axis lenght (maxsma = 20.) maximum semi-major axis lenght (step = 2.) sma step between successive ellipses (linear = yes) linear sma step ? (maxrit = INDEF) maximum sma lenght for iterative mode (recente= no) allows finding routine to re-center x0-y0 ? (xylearn= no) updates pset with new x0-y0 ? (physica= no) physical coordinate system ? (mode = al) ------- ellip haro29fblkav.fit fsky_ell ellip haro29nblkav.fit nsky_ell tprint fsky_ell.tab columns='sma,intens,pix_var,ndata' > fsky_ell.lst tprint nsky_ell.tab columns='sma,intens,pix_var,ndata' > nsky_ell.lst # I choose the surface photometry intensity on the 1.75/2.25 arcmin annulus (3/4th annul in fsky_ell.lst/nsky_ell.lst) # The sky variations around the galaxy are -- 5.37E-5 for FUV, 2.20E-4 for NUV imarith haro29fcl.fit - 2.83E-4 haro29fms.fit imarith haro29ncl.fit - 0.00215 haro29nms.fit =========> haro29nms.fit; haro29fms.fit #Add SKYMEAN keyword to these fits headers sethead haro29fms.fit skymean= 2.83E-4 sethead haro29nms.fit skymean= 0.00215 ------------------------------------------- #Geometrically transform the NUV and FUV images to match the pixel scale and orientation of # the original V image: #I used vn.geomap and vn.xy used by Deidre imheader haro29v.fits haro29v.fits[487,476][real]: TEMP55[1/1] =====> V image is 487 * 476. geotran haro29fms.fit haro29fms_geo.fit vn.geomap vn.xy geotran haro29nms.fit haro29nms_geo.fit vn.geomap vn.xy geotran haro29fcut.fit haro29fcut_geo.fit vn.geomap vn.xy geotran haro29ncut.fit haro29ncut_geo.fit vn.geomap vn.xy =====> haro29fms_geo.fit; haro29nms_geo.fit; haro29fcut.fit; haro29ncut_geo.fit ---------------------------------------------------------------------------------------- Ellipse photometry: ### #Use haro29n.tab as inellip to do photometry # copy haro2970fms_geo.pl/haro29nms_geo.pl (the original mask file) to the new GALEX directory # as haro29fms_geo.pl/haro29nms_geo.pl. ellipse haro29fms_geo.fit haro29f_ell inellip=haro29n.tab inter- region+ ellipse haro29nms_geo.fit haro29n_ell inellip=haro29n.tab inter- region+ ====> haro29f_ell.tab; haro29n_ell.tab isoexam haro29f_ell.tab isoexam haro29n_ell.tab ====> haro29f_ell.ps; haro29n_ell.ps tprint haro29f_ell.tab > haro29f.phot tprint haro29n_ell.tab > haro29n.phot ====> haro29f.phot; haro29n.phot ---------------------------------------------------------------------------------------- Create input_galex file: FUV photometry file: 0=no FUV NUV photometry file: 0=no NUV Final output file Number of ellipses Integration time --- FUV Integration time --- NUV Average sky (counts/s/pixel) --- FUV Average sky (counts/s/pixel) --- NUV b/a D(Mpc) E(B-V) total scale of V image in "/pixel R_25 (arcsec), 9999.9 if none R_D (arcsec), 9999.9 if none The original input_galex file should work fine, but change to whatever filenames you are using. haro29f.phot haro29n.phot haro29uv.final82 12 1647.62 1538.7 4.64955e-05 0.000353234 0.58 5.4 0.05 0.608 29.05 10.32 The sky is determined by scaling the sky to the V-band pixel size. scale as (V-band pixel size/1.5)**2. =====> input_galex ---------------------------------------------------------------------------------------- Calibrate the photometry and determine surface brightnesses: ../annuli_galex82.exe =====> haro29uv.final82 ---------------------------------------------------------------------------------------- Plot: uv_82.sm sm input uv_82.sm =====> haro290uv_82.eps -------------------------------------- Fit exponential disk: fields haro29uv.final82 2,5 lines="21-28" | polyfit stdin 1 > haro29n_.expfit This gives you mu_NUV = mu_NUV^0 + slope * R(arcmin) ../expfit.exe ### This calculates R_D for you. give it haro29n_.expfit mv expfit.out82 haro29n_.expout ===========> R_D= 0.132 +/- 0.010 arcmin R_D= 0.207 +/- 0.016 kpc mu_NUV_0= 21.591 +/- 0.300 y(100kpc)= 545.72 ----------------------------------- Color: ../avecolors.exe Final file=haro29uv.final82 Beg line for FUV-NUV(0=none)=21 End line for FUV-NUV(0=none)=25 =======>Ave FUV-NUV=0.061 +/- 0.153 fields haro29uv.final82 3,9 lines="21-25" | polyfit stdin 1 > haro29fmn_.expfit This calculates FUV-NUV = zeropt + slope * R(kpc) =====> zeropt = -0.16645 +/- 0.09684131; slope = 0.568125 +/- 0.2107239 y(100kpc)= 56.6460 ---------------- Revise the FUV-NUV lines in uv_82.sm update : line 107 , line 108 for NUV SB fitting plot. update : line 182 , line 183 for UV color plot with new Ave color Redraw haro29uv_82.eps ------------------------------------------ Combine with other photometry: Enter in galex2/photfiles_rev.list and cp to tempfiles.list and edit out all but galaxy that you are working on, or do them all at the end. ./collectv_lt.exe =====> haro29allphot82.dat