ELLIPSECT: A surface brightness analysis tool for GALFIT 3

This paper introduces EllipSect, a Python-based tool that complements GALFIT 3 by generating publication-quality visualizations, exportable data, and non-parametric measurements to facilitate the quantitative assessment of galaxy surface brightness models.

The Gist

Imagine you are an astronomer trying to understand the shape and structure of a distant galaxy. You have a powerful tool called GALFIT, which is like a master chef that can "cook up" a mathematical recipe to describe how bright a galaxy is in every single spot. It breaks the galaxy down into ingredients like a central bulge, a swirling disk, or a bar.

However, GALFIT has a problem: it's a bit of a "black box." It gives you a list of numbers and a complex data file, but it doesn't give you a nice, easy-to-read picture of how well that recipe actually matches the real galaxy. It's like the chef handing you a list of ingredients and saying, "Here's the recipe," but refusing to let you taste the dish to see if it's good.

Enter EllipSect.

What is EllipSect?

Think of EllipSect as the tasting spoon and the food critic for the GALFIT chef. It is a new, user-friendly tool (written in the Python programming language) that takes GALFIT's raw data and turns it into clear, beautiful pictures and easy-to-understand measurements.

Here is how it works, using some everyday analogies:

1. The "Slice and Dice" Method (The Algorithm)

Imagine the galaxy is a giant, slightly squashed pizza.

• GALFIT tries to describe the whole pizza mathematically.
• EllipSect takes that pizza and cuts it into many thin, pie-shaped slices (sectors) radiating from the center.
• It then measures the "cheese" (light) on each slice at different distances from the center.
• By averaging these slices, it creates a smooth, easy-to-read graph that shows exactly how the brightness changes from the center of the galaxy to its edge. This helps astronomers see if the mathematical recipe (the model) actually fits the real pizza (the galaxy).

2. The "Quality Control" Check

When you bake a cake, you want to know if it's too dry or if the frosting is uneven.

• EllipSect creates a "Residual Image." Imagine taking a photo of the real galaxy and subtracting the "recipe" photo from it. What's left? The mistakes!
• If the recipe was perfect, the residual image would be blank (just a black sky).
• If there are bright spots or dark holes in the residual image, it means the recipe is missing something. Maybe the galaxy has a hidden dust lane or an extra arm that the model didn't catch. EllipSect highlights these errors so astronomers can fix their models.

3. The "Total Size" Calculator

GALFIT is great at measuring the size of individual parts (like just the bulge or just the disk). But what if you want to know the size of the entire galaxy?

• Imagine a galaxy made of a bright core, a fuzzy halo, and a faint outer ring. GALFIT might give you the size of the core and the size of the ring separately.
• EllipSect acts like a smart calculator that adds all those pieces together to tell you the Total Effective Radius. It answers the question: "How big is the galaxy if we include 90% of its light?" This is crucial for comparing galaxies to each other.

4. The "Sky Background" Detective

One of the hardest things in astronomy is knowing how bright the "empty" space around the galaxy actually is. If you guess the background brightness wrong, your measurement of the galaxy will be wrong.

• EllipSect has two clever tricks to figure this out:
  • The Gradient Method: It looks at rings of sky moving outward from the galaxy, like ripples in a pond, to find where the galaxy's light fades away and the true sky begins.
  • The Random Box Method: It throws random "boxes" of pixels around the galaxy (like throwing darts at a board) to sample the sky and calculate an average.
• This ensures the final measurements aren't skewed by a bad guess about the background.

Why is this a big deal?

Before EllipSect, astronomers had to use old, clunky software (like IRAF) that was hard to use and often required converting files back and forth, which was slow and prone to errors.

EllipSect is like upgrading from a typewriter to a modern word processor. It:

• Saves Time: It automates the boring math.
• Saves Sanity: It produces "publication-quality" graphs that look great in scientific papers.
• Adds Smarts: It calculates special numbers (like the "Bumpiness" of a galaxy or how "tidal" its shape is) that GALFIT doesn't provide on its own.

The Bottom Line

EllipSect is the bridge between complex mathematical models and human understanding. It takes the raw, messy output of a galaxy-fitting program and turns it into a clear story, helping astronomers see the true shape, size, and structure of the universe's building blocks. It's the tool that turns a spreadsheet of numbers into a picture that tells a story.

Technical Summary

1. Problem Statement

While GALFIT 3 is the standard tool for 2D surface brightness (SB) modeling of galaxies, it lacks a dedicated, user-friendly graphical interface to visualize and quantitatively assess its own output.

• Visualization Gap: Users traditionally rely on external tools (like IRAF's ellipse) to compare input images with GALFIT models. However, IRAF is deprecated, and converting GALFIT's complex output formats (FITS cubes and parameter files) for use with other codes is time-consuming and error-prone.
• Limitations of 1D Approaches: Traditional 1D isophotal fitting (e.g., IRAF ellipse) struggles with crowded fields (clusters, compact groups) where galaxies overlap. It often requires masking, which fails to account for the extended light of neighboring galaxies. Furthermore, 1D methods assume varying ellipticity and position angles with radius, which can be problematic when analyzing simultaneous multi-component fits where components are independent.
• Missing Metrics: GALFIT 3 provides parameters for individual components (e.g., effective radius $r_e$ for a bulge) but does not calculate global, non-parametric metrics for the entire multi-component system, such as the total effective radius ($r^t_e$), Petrosian radius, or specific goodness-of-fit statistics like AIC and BIC tailored to the model.

2. Methodology

EllipSect is a Python-based tool designed as a post-fit analysis utility that ingests GALFIT 3 output files directly without requiring manual format conversion.

• Input Data: It reads the GALFIT 3 output FITS cube (original image, model image, residual image) and the parameter file (galfit.nn).
• Sector-Based Photometry: Instead of fitting ellipses to isophotes, EllipSect divides the galaxy and its model into angular sectors around the center using the sectors_photometry routine from the mgefit library.
  • It uses the geometric parameters (position angle, ellipticity) from the last component of the GALFIT model (typically the outer disk or exponential component) as fixed defaults for the entire analysis, though these are user-adjustable.
  • It calculates surface brightness (SB) at multiple radii within these sectors.
• Sky Background Estimation: EllipSect offers two methods to compute the sky background, which GALFIT 3 lacks:
  1. Sky Gradient Method: Constructs concentric elliptical rings and analyzes the gradient of mean sky values to find the stable region (following Barden et al., 2012).
  2. Random Box Method: Places random boxes in the field outside the galaxy mask to calculate the median sky value (following Gao & Ho, 2017).
• Non-Parametric Calculations: It computes metrics that GALFIT 3 does not provide, including:
  • Total Effective Radius ($r^t_e$): Calculated by summing the flux of all components and finding the radius enclosing 50% of the total light via the bisection method.
  • Petrosian and Kron Radii.
  • Goodness-of-Fit: Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC), including an alternative BIC based on resolution elements.
  • Structural Parameters: Bumpiness (irregularity), Tidal parameter (deviation from smooth model), and Signal-to-Noise (S/N).
• Corrections: It automatically retrieves distance moduli and galactic extinction from NED to calculate absolute magnitudes and luminosities. It applies corrections for surface brightness dimming but excludes $k$-corrections (which require color data).

3. Key Contributions

• Integrated Workflow: Eliminates the need for manual data conversion between GALFIT 3 and plotting tools, providing publication-quality figures (SB profiles, residuals, and component breakdowns) directly from GALFIT output.
• Multi-Component Analysis: Enables the visualization and quantitative assessment of individual components (e.g., bulge, disk, bar) within a single fit, allowing users to identify redundant or irrelevant components.
• Global Metrics for Complex Systems: Introduces the calculation of the total effective radius ($r^t_e$) for multi-component fits, a critical parameter for comparing galaxies with numerical simulations or assessing metallicity gradients, which cannot be derived simply from individual component $r_e$ values.
• Robust Sky Determination: Provides independent sky background estimation methods to validate or refine the sky value used in the GALFIT fit, reducing parameter degeneracy.
• Automation Ready: Designed to be scriptable and compatible with automated fitting pipelines (e.g., the DGCG code for cluster galaxies).

4. Results and Examples

The paper demonstrates EllipSect's capabilities using three distinct datasets:

1. Holm 15A (Abell 85 BCG): A cD galaxy fitted with a 7-Gaussian Multi-Gaussian Expansion (MGE) model. EllipSect successfully plotted the SB profile, individual Gaussian components, and residuals, showing the model's ability to capture the galaxy's extended core.
2. M61 (2MASS Ks-band): A spiral galaxy fitted with Sérsic and Exponential components. The tool generated multi-plot visualizations showing SB variations at different azimuthal angles, highlighting regions where the model deviates from the data.
3. NGC 4261 (HST): A single Sérsic fit to an HST image. The tool produced high-resolution plots of the galaxy, model, and residual, demonstrating its utility with space-based data.

In all cases, EllipSect generated publication-ready plots and exported data files containing photometric variables (magnitudes, radii, AIC/BIC) for further analysis.

5. Significance

• Enhanced Reliability: By providing a direct link between 2D modeling and quantitative assessment, EllipSect helps astronomers verify that GALFIT models are physically appropriate for the data, particularly in complex, crowded environments where 1D methods fail.
• Standardization: It standardizes the calculation of global parameters (like $r^t_e$) for multi-component galaxies, facilitating more accurate comparisons between observational data and theoretical models.
• Accessibility: As a Python-based, open-source tool, it lowers the barrier to entry for analyzing GALFIT outputs, replacing deprecated IRAF workflows with modern, reproducible code.
• Future Scalability: The tool is designed to handle large datasets (hundreds of galaxies in clusters), making it a vital component for upcoming large-scale surveys and automated galaxy morphology studies.

In summary, EllipSect bridges the gap between 2D parametric modeling and detailed photometric analysis, offering a comprehensive suite of tools to validate, visualize, and quantify galaxy structures derived from GALFIT 3.