def PlotFq_t(lit='wetzelsmooth'):
''' Plot the quiescent fraction evolution as a function of M* and t
'''
#prettyplot()
pretty_colors = prettycolors()
M_arr = np.arange(9.0, 12.5, 0.5)
t_arr = np.arange(5.0, 14.0, 0.2)
fig = plt.figure()
sub = fig.add_subplot(111)
for Mstar in M_arr:
fq_M_t = []
for tt in t_arr:
qf = gp.Fq()
fq_M_t.append(
qf.model(np.array([Mstar]), Util.get_zsnap(tt), lit=lit))
sub.plot(t_arr, fq_M_t, label=r"$M_* = " + str(round(Mstar, 2)) + "$")
sub.set_xlim([t_arr.min(), t_arr.max()])
sub.set_ylim([0., 1.])
#plt.show()
fig_file = ''.join(['figure/test/', 'fq_t.', lit, '.png'])
fig.savefig(fig_file, bbox_inches='tight')
return None
def PlotFq_WetzelComp():
''' Compare the quiescent function evolution for different analytic prescriptions
of Wetzel et al. (2013)
'''
lit = ['wetzel', 'wetzel_alternate']
M_arr = np.arange(9.0, 12.2, 0.2)
#z_arr = np.arange(0.0, 1.2, 0.2)
z_arr = np.array([0.0, 0.36, 0.66, 0.88])
prettyplot()
pretty_colors = prettycolors()
fig = plt.figure()
sub = fig.add_subplot(111)
qf = gp.Fq()
for il, l in enumerate(lit):
if il == 0:
ls = '-'
elif il == 1:
ls = '--'
for iz, z in enumerate(z_arr):
if iz == 0:
label = l.title()
else:
label = None
fq_M = qf.model(M_arr, z, lit=l)
sub.plot(M_arr,
fq_M,
c=pretty_colors[iz],
lw=3,
ls=ls,
label=label)
# Now plot Tinker et al. (2013)'s fQ^cen
for iz, z in enumerate(z_arr):
if z != 0.0:
fq_file = ''.join(
['dat/observations/cosmos_fq/', 'stats_z',
str(iz), '.fq_cen'])
mmm, fqcen = np.loadtxt(fq_file, unpack=True, usecols=[0, 1])
mmm = np.log10(mmm)
else:
fq_file = ''.join(
['dat/observations/cosmos_fq/', 'fcen_red_sdss_scatter.dat'])
mmm, fqcen = np.loadtxt(fq_file, unpack=True, usecols=[0, 1])
sub.scatter(mmm, fqcen, c=pretty_colors[iz], lw=0, s=10)
sub.set_xlim([9.0, 12.])
sub.set_xlabel(r'$\mathtt{M_{*}}$', fontsize=25)
sub.set_ylim([0., 1.])
sub.set_ylabel(r'$\mathtt{f_Q}$', fontsize=25)
sub.legend(loc='upper left')
fig_file = ''.join(['figure/test/', 'Fq_Wetzel_Comparison.png'])
fig.savefig(fig_file, bbox_inches='tight')
plt.close()
return None
def Plot_dFqdt():
''' Plot the evolution of the derivative of the quiescent fraction as a function of M* and t
for different parameterized models of f_Q
d f_Q / d t
'''
prettyplot()
pretty_colors = prettycolors()
dt = 0.01
M_arr = np.arange(9.0, 12.5, 0.5)
t_arr = np.arange(6.0, 13.0, 0.2)
fig = plt.figure()
sub = fig.add_subplot(111)
for iM, Mstar in enumerate(M_arr):
qf = gp.Fq()
lit = 'wetzelsmooth'
sub.plot(t_arr, [gp.dFqdt(Mstar, tt, lit=lit, dt=dt) for tt in t_arr],
ls='-',
lw=3,
c=pretty_colors[iM],
label=r"$\mathtt{M}_* = " + str(round(Mstar, 2)) + "$")
lit = 'wetzel'
sub.plot(t_arr, [gp.dFqdt(Mstar, tt, lit=lit, dt=dt) for tt in t_arr],
ls=':',
lw=3,
c=pretty_colors[iM])
lit = 'cosmosinterp'
t_blah = np.arange(8.0, 9.6, 0.1)
sub.plot(t_blah,
[gp.dFqdt(Mstar, tt, lit=lit, dt=dt) for tt in t_blah],
ls='--',
lw=3,
c=pretty_colors[iM])
lit = 'cosmosfit'
t_blah = np.arange(8.0, 9.6, 0.1)
sub.plot(t_blah,
[gp.dFqdt(Mstar, tt, lit=lit, dt=dt) for tt in t_blah],
ls='-.',
lw=3,
c=pretty_colors[iM])
sub.set_xlim([t_arr.min(), t_arr.max()])
sub.legend(loc='upper left')
sub.set_ylim([0., 0.2])
#plt.show()
fig_file = ''.join(['figure/test/', 'dfqdt.png'])
fig.savefig(fig_file, bbox_inches='tight')
return None
def PlotFq_WetzelTinker_comparison():
''' Compare the quiescent function evolution for different analytic prescriptions
'''
lit = ['wetzel', 'wetzel_alternate', 'cosmos_tinker']
prettyplot()
pretty_colors = prettycolors()
M_arr = np.arange(9.0, 12.2, 0.2)
z_arr = np.arange(0.0, 1.2, 0.2)
prettyplot()
pretty_colors = prettycolors()
fig = plt.figure(figsize=(20, 5))
qf = gp.Fq()
for il, l in enumerate(lit):
sub = fig.add_subplot(1, 4, il + 1)
for iz, z in enumerate(z_arr):
if iz == 0:
label = (l.replace('_', ' ')).title()
else:
label = None
fq_M = qf.model(M_arr, z, lit=l)
sub.plot(M_arr, fq_M, c=pretty_colors[iz], lw=3, label=label)
sub.set_xlim([9.0, 12.])
sub.set_xlabel(r'$\mathtt{M_{*}}$', fontsize=25)
sub.set_ylim([0., 1.])
if il == 0:
sub.set_ylabel(r'$\mathtt{f_Q}$', fontsize=25)
else:
sub.set_yticklabels([])
sub.legend(loc='upper left')
fig_file = ''.join(
['figure/test/', 'Fq_Comparison.', '_'.join(lit), '.png'])
fig.savefig(fig_file, bbox_inches='tight')
plt.close()
return None
def Plot_rhoSFR(source='li-march',
sfms_prop=None,
fq_lit='wetzel',
zslope=None):
''' Calculate the SFR density by integrating
rho_SFR = int n(M*) SFR_MS(M*) (1 - fQ(M*)) dM*
where n(M*) is given by the stellar mass function, SFR_MS is the SFMS SFR
fQ is the quiescent fraction
'''
if sfms_prop is None:
if zslope is None:
sfms_prop = {'name': 'linear', 'zslope': 1.5}
else:
sfms_prop = {'name': 'linear', 'zslope': zslope}
# SMF
smf_obj = gp.SMF()
n_smf = smf_obj.analytic
# SFMS Star Formation Rate
sfr_ms = sfr_evol.AverageLogSFR_sfms
# Quiescent Fraction
qf_obj = gp.Fq()
fq_model = qf_obj.model
def rhoSFR_integrand(Mstar, z):
M_smf, phi_smf = n_smf(z, source=source)
n_smf_M = interpolate.interp1d(M_smf, phi_smf)
return n_smf_M(Mstar) * 10.**sfr_ms(Mstar, z, sfms_prop=sfms_prop) * (
1. - fq_model(Mstar, z, lit=fq_lit))
z_range = np.array([0.1, 0.5, 1.0]) #np.arange(0.1, 1.0, 0.2)
rhoSFR_z = []
for zz in z_range:
print zz
rhoSFR_int_wrap = lambda mm: rhoSFR_integrand(mm, zz)
rhoSFR_int, rhoSFR_int_err = quad(rhoSFR_int_wrap,
7.5,
12.0,
epsrel=0.1)
rhoSFR_z.append(rhoSFR_int)
z_behroozi = [
0.036209892874743854, 0.40889165564288144, 0.6761600913028816,
0.9616372266749669, 1.170544445309055, 1.4822192825461054,
1.6049075298939122
]
sfr_cosmic = [
0.010063810047184728, 0.02228775362744163, 0.03802861399461939,
0.06050175265988725, 0.0805521892279678, 0.10063810047184728,
0.10862045113189879
]
prettyplot()
pretty_colors = prettycolors()
fig = plt.figure()
sub = fig.add_subplot(111)
sub.scatter(z_range, rhoSFR_z, s=10)
sub.plot(z_range, rhoSFR_z, c='k', lw=2)
sub.plot(z_behroozi, sfr_cosmic, lw=3, c='r', label='Behroozi Best-fit')
sub.set_xlabel(r'Redshift $\mathtt{(z)}$', fontsize=25)
sub.set_xlim([0.0, z_range.max()])
sub.set_ylim([0.001, 0.1])
sub.set_ylabel(r'$\rho_\mathtt{SFR}$', fontsize=25)
sub.set_yscale('log')
sub.legend(loc='upper left')
plt.show()
fig_name = ''.join([
'figure/test/', 'rhoSFR', '.SMF_', source, '.fq_', fq_lit,
'.SFMS_zslope',
str(round(sfms_prop['zslope'], 1)), '.png'
])
fig.savefig(fig_name, bbox_inches='tight')
plt.close()
return None