Beispiel #1
0
        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)))
Beispiel #2
0
'''
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