phangsteam / megatable Goto Github PK

Tools for generating the PHANGS multi-wavelength data tables (a.k.a., mega-tables)

License: GNU General Public License v3.0

Python 100.00%

megatable's Introduction

MegaTable

This repository contains the source code for generating rich multiwavelength data tables (a.k.a. the "mega-tables") for the PHANGS team. The table construction and data aggregation schemes are described in the following papers:

Sun, Leroy, Rosolowsky, et al. (2022), "Molecular Cloud Populations in the Context of Their Host Galaxy Environments: A Multiwavelength Perspective"
Sun, Leroy, Ostriker, et al. (2023), "Star Formation Laws and Efficiencies across 80 Nearby Galaxies"

The current version of the repository corresponds to PHANGS mega-tables version 4.0.

Code Architecture

This repository consists of a python module named mega_table and a suite of python scripts and configuration files in the pipelines subdirectory.

The mega_table module provides the core infrastructure for table creation, manipulation, and input/output. This module relies heavily on the astropy.table subpackage. Most of the tools in this module are offered through three python classes:

mega_table.table.RadialMegaTable: Assemble measurements in radial bins with a given width
mega_table.table.TessellMegaTable: Assemble measurements according to a given tessellation pattern
mega_table.table.ApertureMegaTable: Assemble measurements in arbitrarily placed, fixed size apertures

The pipelines subdirectory includes the actual production code for the PHANGS mega-tables:

config_data_path.json: This file specifies the path to the underlying datasets on disk
config_tables.json: This file provides input parameters for table creation (e.g., table naming convention, bin width, FoV extent)
format_*.csv: These files controls the column-by-column content of the output tables (e.g., column names, physical units, descriptions)
make_*.py: These are the python scripts that define and incorporate the multiwavelength measurements in the PHANGS mega-tables

Contact

If you need help using the code in this repository for science applications, please reach out to Jiayi Sun. For bug reports and improvement suggestions, please open an Issue on Github.

megatable's People

Contributors

Stargazers

Watchers

Forkers

astrojysun ysong-astro ceibensteiner

megatable's Issues

Calculate average T_peak values in each aperture

Per Eva's request.

Adjust the "default" SFR prescription lists

Right now I am using FUV+W4 as default, or W4-ONLY if there is no FUV.

Should be using FUV+W4, then NUV+W4, then W4-ONLY.

Add alpha_vir estimates for the "cylinder" cloud model

This is mainly about the alpha_vir estimates for the CPROPS clouds. But of course the assumed l.o.s. depth should match.

-

Update input parameters

Latest galaxy masses and radii from Adam
R21 = 0.65

Incorporate team feedback

2021 team meeting breakout session minutes
2023 team meeting breakout session minutes

Add uncertainties on key cloud-scale/kpc-scale quantities

For quantities derived from PHANGS-ALMA data: use published uncertainty maps
For those derived from CPROPS catalogs: are there uncertainties in the catalogs?
For those derived from PHANGS-VLA: ask for Dyas's help to generate uncertainty maps?
For those derived from z0MGS: estimate statistical uncertainties from empty regions; what about systematic uncertainties?
For those derived from S4G: similar to z0MGS?
For those derived from narrow band Halpha: what to do?
For those derived from CO rotation curves: use uncertainties associated with the published rotation curves

Add rotation curve related quantities and refactor pipeline

Circular rotation speed V_circ and/or angular speed Omega_circ
Local shear parameter q = -d ln(Omega_circ) / d ln(r_gal)

Include Sigma_star estimates from a varying M/L

To incorporate the M/L maps from z0MGS and apply them on S4G ICA / IRAC1 maps.

Add fiducial estimate of averaged free-fall time on cloud scales

Assuming a fixed depth of H=100pc

To update Sigma_SFR estimates from narrow-band Halpha

Resample the convolved maps, and call them Sigma_SFR_HaONLY and Sigma_SFR_HaW4
Calculate the area-weighted average, and call them A<Sigma_SFR_HaONLY> and A<Sigma_SFR_HaW4>

Put together an ascii text file listing all columns and descriptions
Let the table generation code to put them into the

Best practice on the S4G stellar component maps?

Need to write a small piece of code to do the following selection:

Use the background subtracted and masked original 3.6um maps when ICA is not recommended (column Excluded? in "P5_table");
Choose between ICA1 and ICA2 when ICA is recommended (column ICA iter in "P5_table").

Also, need to think about the interactions between the ICA decomposition and the varying M/L ratio.

Channel width / spectral response curve correction

No correction is applied right now. Need to add this.

Add a warning message when coordinates get close to RA = 0/360

The code uses small angle approximation. It will break when a galaxy footprint cross the RA = 0/360 line...