Quick Start#

Let’s get started MUFASA. Use the navigation bar on the right to skip ahead as needed.

Run Spectral Fits#

Note

Installation: If you haven’t installed MUFASA already, see Install.

Setting Up#

To set up spectral fits, initialize a Region object to manage the input data, model parameters, and output files in specified directories:

from mufasa import master_fitter as mf

# Define input and output paths
cubePath = "path/to/cube.fits"   # Path to the data cube
paraDir = "output_dir"           # Directory for the output files
paraNameRoot = "results"         # Prefix for the output file names

# Specify the spectral model type
fittype = 'nh3_multi_v'

# Initialize the Region object
region = mf.Region(cubePath, paraNameRoot, paraDir, fittype)

MUFASA currently supports the following spectral models to be specified with fittype:

  • NH₃ (1,1) (2-components; 'nh3_multi_v')

  • N₂H⁺ (1-0) (2-components; 'n2hp_multi_v')

Fitting#

To perform automated spectral fitting:

# Perform an automated two-components fit
region.master_2comp_fit(snr_min=3)

Note

SNR cutoff: The default snr_min=3 is recommended, which typically goes deeper than SNR=3 fairly efficiently. MUFASA supports snr_min=0, but the effort is much more exhaustive.

View Results#

Reading Saved Results#

If you have exited the fitting session already, you can quickly load the saved results by reinitializing region (see the Fitting Spectra section) and running:

# Load the saved fits for the 1 and 2 component models
region.load_fits(ncomp=[1, 2])

Plotting Spectral Fits#

To plot the fitted spectra for a specific position (x, y):

# Plot the fits for position (x, y) on a 3x3 grid for a 2-component model
region.ucube.plot_fits_grid(x, y, ncomp=2, size=3, xlim=None, ylim=None)

3D Scatter Plots#

To plot the fitted parameters in 3D, such as in position-position-velocity (PPV) space:

# Plot fitted results as an interactive 3D HTML file, saved to `savepath`
region.plot_ppv_scatter(savepath, vel_scale=0.5, showfig=True, auto_open_html=False)

Explore More#

  • For more examples and detailed guides, see Tutorials.

  • For a full reference of available modules, visit the API Reference.