##Note this galaxy has very deep UV obsess, so need large aperture to do surface photometry Run D63 GALEX images through to get SB HXZ: 02/08/10 Get exposure times: (Check exp FITS files) GI3_050003_HolmbergII-fd-exp.fits 22591.6 GI3_050003_HolmbergII-nd-exp.fits 21083.5 ======> fuv: 22591.6 s; nuv: 21083.5 s 1)Cut out smaller images: cutout subimage size: 800*800 centered around the galaxy center of cutout: 1932; 1909 imcopy GI3_050004_HolmbergI-fd-int.fits[1530:2329,1510:2309] d63fcut.fit imcopy GI3_050004_HolmbergI-nd-int.fits[1530:2329,1510:2309] d63ncut.fit -------------------------------------------- 2)Remove foreground stars and background galaxies: a)use Daofind to find point sources on d50fut.fit and d216ncut.fit respectively. and exclude all sources on top of the galaxy ======> f_tgmsk.fit.coo.1; n_tgmsk.fit.coo.1 txdump n_tgmsk.fit.coo.1 xcenter,ycenter > ncut.reg txdump f_tgmsk.fit.coo.1 xcenter,ycenter > fcut.reg imedit d63ncut.fit ntmp.fit cursor=ncut.reg logfile=n1.log radius=6 display- imedit d63fcut.fit ftmp.fit cursor=fcut.reg logfile=f1.log radius=6 display- b)check V-band image to further mark some point sources on top of the galaxy. Then display FUV/NUV color images to include only extreme red sources marked on V image, In the meanwhile, further imedit out any residuals on the off-galaxy regions. imedit ntmp.fit d63ncl.fit imedit ftmp.fit d63fcl.fit ========> d63ncl.fit; d63fcl.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 d63fcl.fit d63fblkav.fit 10 10 blkavg d63ncl.fit d63nblkav.fit 10 10 ========> d63fblkav.fit; d63nblkav.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. #The sky regions are located as where the intensities fall into the background noise PACKAGE = isophote TASK = geompar (x0 = 41.414307) initial isophote center X (y0 = 40.834485) initial isophote center Y (ellip0 = 0.05) initial ellipticity (pa0 = -0.2) initial position angle (degrees) (sma0 = 10.) initial semi-major axis lenght (minsma = 10.) minimum semi-major axis lenght (maxsma = 40.) maximum semi-major axis lenght (step = 3.) 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-y (xylearn= no) updates pset with new x0-y0 ? (physica= no) physical coordinate system ? (mode = al) ------- ellip d63fblkav.fit fsk_ell ellipse d63nblkav.fit nsk_ell tprint fsk_ell.tab columns='sma,intens,pix_var,ndata' > fsky_ell.lst tprint nsk_ell.tab columns='sma,intens,pix_var,ndata' > nsky_ell.lst # I choose the surface photometry intensity on the 5.5 arcmin annulus (5th annul in *_ell.lst) # The variations of the sky mean around the galaxy are -- 2.88E-5 for FUV, 1.00E-4 for NUV imarith d63fcl.fit - 3.06E-4 d63fms.fit imarith d63ncl.fit - 0.00232 d63nms.fit =========> d63nms.fit; d63fms.fit #Add SKYMEAN keyword to fits header sethead d63fms.fit skymean=3.06E-4 sethead d63nms.fit skymean=0.00232 ------------------------------------------- #Geometrically transform the NUV and FUV images to match the pixel scale and orientation of # the original V image: #I use 9 field stars in v and nuv images center d63ncut.fit output=d63ncut.fit.ctr.1 center d63v.fits output=d63v.fits.ctr.1 txdump d63v.fits.ctr.1 xc,yc yes > v.xy txdump d63ncut.fit.ctr.1 xc,yc yes > n.xy joinlines v.xy n.xy > vn.xy imheader d63v.fits d63v.fits[956,932][real]: d63 geomap vn.xy vn.geomap 1 956 1 932 fitgeo=general functio=polyn xxorder=2 xyorder=2 xxterms=half =====> V image is 956 * 932. rms in geomap is 0.093 and 0.093 for x and y geotran d63fms.fit d63fms_geo.fit vn.geomap vn.xy geotran d63nms.fit d63nms_geo.fit vn.geomap vn.xy geotran d63fcut.fit d63fcut_geo.fit vn.geomap vn.xy geotran d63ncut.fit d63ncut_geo.fit vn.geomap vn.xy =====> d63fms_geo.fit; d63nms_geo.fit; d63fcut.fit; d63ncut_geo.fit ---------------------------------------------------------------------------------------- Ellipse photometry: #Use d63ncl.tab generated previously to do photometry ### # copy ddo63fms_geo.pl/ddo63nms_geo.pl (the original mask file) to the new GALEX directory # as d63fms_geo.pl/d63nms_geo.pl. ellip d63fms_geo.fit d63f_ell inellip=d63ncl.tab inter- region+ ellip d63nms_geo.fit d63n_ell inellip=d63ncl.tab inter- region+ ====> d63f_ell.tab; d63n_ell.tab isoexam d63f_ell.tab isoexam d63n_ell.tab ====> d63f_ell.ps; d63n_ell.ps tprint d63f_ell.tab > d63f.phot tprint d63n_ell.tab > d63n.phot ====> d63f.phot; d63n.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. d63f.phot d63n.phot d63uv.final82 14 22591.6 21083.5 0.000174890 0.00132596 0.95 3.8 0.0625 1.134 77.49 163.2 The sky is determined by scaling the sky to the V-band pixel size. Then scale as (V-band pixel size/1.5)**2. =====> input_galex ---------------------------------------------------------------------------------------- Calibrate the photometry and determine surface brightnesses: ../annuli_galex82.exe =====> d63uv.final82 ---------------------------------------------------------------------------------------- Plot: uv_82.sm sm input uv_82.sm =====> d75uv_82.eps -------------------------------------- Fit exponential disk: fields d63uv.final82 2,5 lines="23-29" | polyfit stdin 1 > d63n_inner.expfit This gives you mu_NUV = mu_NUV^0 + slope * R(arcmin) ../expfit.exe ### This calculates R_D for you. give it d63n_inner.expfit mv expfit.out82 d63n_inner.expout fields d63uv.final82 2,5 lines="29-36" | polyfit stdin 1 > d63n_outer.expfit ../expfit.exe ### This calculates R_D for you. give it d63n_outer.expfit mv expfit.out82 d63n_outer.expout # Calculate the radius R_Br and surface brightness where they intersect. ../intersect.exe 1st file returned by POLYFIT: d63n_inner.expfit 2nd file returned by POLYFIT: d63n_outer.expfit ======> intersection x0 : 1.167 intersection x0 in Kpc: 1.290 intersection y0: 25.655 # =====> R_Br= 1.157 arcmin = 1.279 kpc at mu_NUV = 25.629 R_D inner = 55.022 +/- 23.106 arcmin = 60.844 +/- 25.551 kpc mu_NUV^0 inner = 25.606 +/- 0.035 R_D outer = 0.422 +/- 0.027 arcmin = 0.466 +/- 0.030 kpc kpc y(100 kpc) outer = 255.49 ----------------------------------- Color: ../avecolors.exe Final file=d63uv.final82 Beg line for FUV-NUV(0=none)=23 End line for FUV-NUV(0=none)=36 =======>Ave FUV-NUV= 0.168 +/- 0.133 fields d63uv.final82 3,9 lines="23-29" | polyfit stdin 1 > d63fmn_inner.expfit This calculates FUV-NUV = zeropt + slope * R(kpc) =====> zeropt = 0.1877857 +/- 0.01030301; slope = -0.1642857 +/- 0.01282723 fields d63uv.final82 3,9 lines="29-36" | polyfit stdin 1 > d63fmn_outer.expfit This calculates FUV-NUV = zeropt + slope * R(kpc) =====> zeropt = -0.2441789 +/- 0.1681947; slope = 0.2263416 +/- 0.07845998 y (100 kpc) = 22.3900 # ../intersect.exe Intersect at R_colBR= 1.106 arcmin = 1.223 kpc at FUV-NUV=0.006 # fields d63uv.final82 3,9 lines="23-36" | polyfit stdin 1 > d63fmn_.expfit ======> zeropt= 0.02050806 +/- 0.05889754; slope= 0.1010692 +/- 0.03489703 y (100 kpc) = 10.1274 ---------------- Revise the FUV-NUV lines in uv_82.sm update : line 107 , line 108 , line 109 for NUV SB fitting plot. update : line 182 , line 183 for UV color plot with new Ave color ------------------------------------------ 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 =====> d63allphot82.dat