GEMS: galaxy fitting catalogs and testing parametric galaxy fitting codes: GALFIT and GIM2D

In the context of measuring the structures of intermediate-redshift galaxies with HST ACS surveys, we tune, test, and compare two widely used fitting codes (GALFIT and GIM2D) for fitting single-component Sersic models to both simulated and real galaxy data. Our study focuses on the GEMS survey with the sensitivity of typical HST survey data, and we include our final catalog of fit results for all 41,495 objects detected in GEMS. Using simulations, we find that fitting accuracy depends sensitively on galaxy profile shape. Exponential disks are well fit and have small measurement errors, whereas fits to de Vaucouleurs profiles show larger uncertainties owing to the large amount of light at large radii. Both codes provide reliable fits with little systematic error for galaxies with effective surface brightnesses brighter than that of the sky; the formal uncertainties returned by these codes significantly underestimate the true uncertainties (as estimated using the simulations). We find that GIM2D suffers significant systematic errors for spheroids with close companions owing to the difficulty of effectively masking out neighboring galaxy light; there appears to be no work-around to this important systematic in GIM2D's current implementation. While this crowding error affects only a small fraction of galaxies in GEMS, it must be accounted for in the analysis of deeper cosmological images or of more crowded fields with GIM2D. In contrast, GALFIT results are robust to the presence of neighbors because it can simultaneously fit the profiles of multiple companions as well as the galaxy of interest. We find GALFIT's robustness to nearby companions and factor of greater than or similar to 20 faster runtime speed are important advantages over GIM2D for analyzing large HST ACS data sets.