p.errorbar(phys_radius, data["Vt"], yerr=data['eVt'],
fmt='-', color='b', label="This work",
drawstyle='steps-mid')
p.plot(phys_radius, vcirc(phys_radius, *plot_pars), 'r-')
p.ylabel(r"V$_{\mathrm{circ}}$ (km / s)")
p.xlabel(r"Radius (kpc)")
p.grid()
p.draw()
p.savefig(paper1_figures_path("M33_vrot_{}_wfit.pdf".format(params)))
p.savefig(paper1_figures_path("M33_vrot_{}_wfit.png".format(params)))
# raw_input("Next plot?")
p.clf()
# load in the Corbelli curve for comparison
corbelli = Table.read(c_hi_analysispath("rotation_curves/corbelli_rotation_curve.csv"))
p.errorbar(phys_radius, data["Vt"], yerr=data['eVt'],
fmt='-', color='b', label="This work", drawstyle='steps-mid')
p.errorbar(corbelli["R"][corbelli["R"] <= 8.0],
corbelli["Vr"][corbelli["R"] <= 8.0],
yerr=corbelli["dVr"][corbelli["R"] <= 8.0],
fmt='--', color='r', label="Corbelli et al. (2014)",
drawstyle='steps-mid')
p.ylabel(r"V$_{\mathrm{circ}}$ (km / s)")
p.xlabel(r"Radius (kpc)")
p.legend(loc='lower right')
p.grid()
p.draw()
p.savefig(paper1_figures_path("M33_vrot_{}_wCorbelli.pdf".format(params)))
'''
Run entire analysis. The cube should already have been exported into a FITS
file from CASA.
'''
import os
import sys
from subprocess import call
from analysis.paths import c_hi_analysispath, fourteenB_HI_data_path
from analysis.constants import moment1_name
os.sys.path.insert(0, c_hi_analysispath(""))
# os.chdir(c_hi_analysispath(""))
#################
# Initial masking.
#################
# Mask w/ pbcov limit.
execfile(c_hi_analysispath("pbcov_masking.py"))
# Create the source mask
execfile(c_hi_analysispath("make_signal_mask.py"))
# Create moment arrays
execfile(c_hi_analysispath("HI_make_moments.py"))
# Make some moment figures