README:
This is a supplement to 'Small Protoplanetary Disks in the Orion Nebula Cluster and OMC1 with ALMA', Otter and Ginsburg et al 2021, ApJ. 

All the code for this work can be found in the github repo: https://github.com/jaotter/summer_research_2018
Though this repo is not well organized, so some of the important scripts have been included here, as well as the reduced data images.


description of files:

fits images:
Orion_SourceI_B3_continuum_r0.5.clean0.05mJy.allbaselines.huge.deepmask.image.tt0.pbcor.fits
Orion_SourceI_B6_continuum_r0.5.clean1mJy.150mplus.huge.image.tt0.pbcor.fits
Orion_SourceI_B7_continuum_r0.5.clean0.05mJy.250klplus.deepmask.image.tt0.pbcor.fits

These are the main reduced images. For the data reduction scripts, see https://github.com/keflavich/Orion_ALMA_2016.1.00165.S, and https://zenodo.org/record/2365835#.YXv_IjopCuU

other fits images:
B6_convolved_r0.5.clean1mJy.150mplus.huge.image.tt0.pbcor.fits
B7_convolved_r0.5.clean0.05mJy.250klplus.deepmask.image.tt0.pbcor.fits

These are the band 6 and 7 images convolved to the band 3 beam, with the script:
convolve_images.py

source measurements:
create_catalog.py - used to create first rough catalog from source positions with fluxes and size measurements.
multi_gauss_fit.py - code to fit multiple Gaussians for background fitting. Based on 'gaussfit_catalog'
fit_background.py - fit source and a background second Gaussian
upper_lim_fluxes.py - code to measure flux upper limits for non-detections
size_upper_lim.py - code to determine upper limit size that can be deconvolved as a function of S/N
calc_inf_vals.py - script to compute inferred values like dust masses

plotting:
multi_SED_plot.py - create Figure 9
spectral_indices.ipynb - create Figure 10
density_plot.py - create Figure 13
min_span_tree.py - create Figure 14