mpl.rcParams['xtick.labelsize'] = 10
mpl.rcParams['ytick.labelsize'] = 10
mpl.rcParams['font.family'] = 'serif'
mpl.rcParams['font.serif'] = 'Times New Roman'
if __name__ == '__main__':
import h5py
from CALIFAUtils.scripts import get_h5_data_masked
baseCode = 'Bzca6e'
#baseCode = 'Bgsd6e'
h5 = h5py.File('/Users/lacerda/dev/astro/SFRRadialProfiles/SFRSDRadialProfiles_minpopx0.05_%s.h5' % baseCode)
N_T = h5['HEADER/N_T'].value.astype(np.int)
tSF__T = h5['HEADER/tSF__T'].value.astype(np.float_)
SFR__Tg = [ get_h5_data_masked(h5, 'SFR__Tg/%d'%iT, **dict(dtype = np.float_)) for iT in xrange(N_T) ]
SFRSD__Tg = [ get_h5_data_masked(h5, 'SFRSD__Tg/%d'%iT, **dict(dtype = np.float_)) for iT in xrange(N_T) ]
SFR_Ha__Tg = [ get_h5_data_masked(h5, 'SFR_Ha__Tg/%d'%iT, **dict(dtype = np.float_)) for iT in xrange(N_T) ]
SFRSD_Ha__Tg = [ get_h5_data_masked(h5, 'SFRSD_Ha__Tg/%d'%iT, **dict(dtype = np.float_)) for iT in xrange(N_T) ]
f = plt.figure()
page_size_inches = [10, 8]
f.set_size_inches(page_size_inches)
f.suptitle(r'%s %d gals minpopx 5%%' % (baseCode, h5['HEADER/N_gals'].value.astype(np.int)))
default_sc_kwargs = dict(marker = 'o', s = 10, alpha = 0.8, edgecolor = 'none', label = '')
default_ols_plot_kwargs = dict(c = 'r', ls = '--', lw = 2, label = 'OLS')
default_ols_kwargs = dict(c = 'r', pos_x = 0.98, pos_y = 0.01, fs = 12, rms = True, text = True, kwargs_plot = default_ols_plot_kwargs)
grid_shape = (2, 2)
xlabel = r'$\log\ SFR^\star(t_\star)\ [M_\odot yr^{-1}]$'
ylabel = r'$\log\ SFR^{neb}\ [M_\odot yr^{-1}]$'
from matplotlib import pyplot as plt from CALIFAUtils.scripts import get_h5_data_masked from CALIFAUtils.plots import plotOLSbisectorAxis mpl.rcParams['font.size'] = 12 mpl.rcParams['axes.labelsize'] = 10 mpl.rcParams['axes.titlesize'] = 12 mpl.rcParams['xtick.labelsize'] = 10 mpl.rcParams['ytick.labelsize'] = 10 mpl.rcParams['font.family'] = 'serif' mpl.rcParams['font.serif'] = 'Times New Roman' h5file = sys.argv[1] h5 = h5py.File(h5file, 'r') SFR_int = get_h5_data_masked(h5, 'SFR_int__Tg', **dict(dtype = np.float_)) SFR_Ha_int = get_h5_data_masked(h5, 'SFR_Ha_int__Tg', **dict(dtype = np.float_)) SFR_Ha_int_masked = get_h5_data_masked(h5, 'SFR_Ha_int_masked__Tg', **dict(dtype = np.float_)) O3N2_int_masked = get_h5_data_masked(h5, 'O3N2M13_int_masked__Tg', **dict(dtype = np.float_)) O3N2_int = np.array(h5['data/O3N2M13_int__g'].value, dtype = np.float_) O3N2_int = np.ma.masked_array(O3N2_int, mask = np.isnan(O3N2_int)) sc_kwargs = dict(marker = 'o', s = 50, edgecolor = 'none', label = '') kwargs_ols_plot = dict(c = 'r', ls = '--', lw = 2, label = 'OLS') kwargs_ols = dict(c = 'k', pos_x = 0.98, pos_y = 0.01, fs = 12, rms = True, text = True, kwargs_plot = kwargs_ols_plot) f = plt.figure() grid_shape = (1, 2) ax1 = plt.subplot2grid(grid_shape, loc = (0, 0)) ax2 = plt.subplot2grid(grid_shape, loc = (0, 1)) xm1, ym1 = C.ma_mask_xyz(x = O3N2_int, y = np.ma.log10(SFR_int[1]) - np.ma.log10(SFR_Ha_int[1]))
try:
h5file = sys.argv[1]
except IndexError:
print 'usage: %s HDF5FILE' % (sys.argv[0])
exit(1)
h5 = h5py.File(h5file, 'r')
RbinCenter__r = h5['HEADER/RbinCenter__r'].value
rbinstep = h5['HEADER/RbinStep'].value
tSF__T = h5[ 'HEADER/tSF__T'].value
NRbins = h5['HEADER/N_rbins'].value
N_gals = h5['HEADER/N_gals']
N_T = h5['HEADER/N_T'].value
gals = h5['HEADER/gals'].value
morf__g = get_h5_data_masked(h5, 'morf__g', **dict(dtype = np.float_))
aSFRSD__Trg = get_h5_data_masked(h5, 'aSFRSD__Trg', **dict(dtype = np.float_))
aSFRSD_Ha__Trg = get_h5_data_masked(h5, 'aSFRSD_Ha__Trg', **dict(dtype = np.float_))
N_zones_notmasked__Tg = get_h5_data_masked(h5, 'N_zones_notmasked__Tg', **dict(dtype = np.float_))
colortipo = ['brown', 'red', 'orange', 'green', '#00D0C9', '#0076C9', 'blue']
Ntype = len(colortipo)
mtypes = [ -2, -1, 0, 1, 2, 3, 4 ]
#mtypes = [ 0, 1, 2, 3, 4, 5, 6 ]
mtype_labels = [ 'E', 'S0', 'Sa', 'Sb', 'Sbc', 'Sc', 'Sd' ]
halfbinstep = np.diff(mtypes)[0]/2.
tickpos = np.linspace(mtypes[0] + halfbinstep, mtypes[-1] - halfbinstep, Ntype)
aSFRSD__Ttr = np.empty((N_T, Ntype, NRbins), dtype = np.float_)
aSFRSD_Ha__Ttr = np.empty((N_T, Ntype, NRbins), dtype = np.float_)
nbinelem_aSFRSD__Ttr = np.empty((N_T, Ntype, NRbins), dtype = np.int_)
