In [1]:
%pylab inline
from __future__ import division
import matplotlib.pyplot as plt
import numpy as np
# Matplotlib related
# Matplotlib default settings
rcdef = plt.rcParams.copy()
pylab.rcParams['figure.figsize'] = 12, 10
pylab.rcParams['xtick.major.size'] = 8.0
pylab.rcParams['xtick.major.width'] = 2.5
pylab.rcParams['xtick.minor.size'] = 4.0
pylab.rcParams['xtick.minor.width'] = 2.5
pylab.rcParams['ytick.major.size'] = 8.0
pylab.rcParams['ytick.major.width'] = 2.5
pylab.rcParams['ytick.minor.size'] = 4.0
pylab.rcParams['ytick.minor.width'] = 2.5
import matplotlib.pyplot as plt
from matplotlib.patches import Ellipse
from matplotlib.ticker import NullFormatter
from matplotlib.ticker import MaxNLocator
In [2]:
mega_u_pivot = 3798.8
hsc_g_pivot = 4782.2
hsc_r_pivot = 6101.7
hsc_i_pivot = 7648.0
hsc_z_pivot = 8883.0
hsc_y_pivot = 9750.8
wfc3_f125w_pivot = 12486.1
wfc3_f160w_pivot = 15370.8
sdss_u = 3543.0
sdss_g = 4770.0
sdss_r = 6231.0
sdss_i = 7625.0
sdss_z = 9134.0
In [3]:
redshift = np.linspace(0.0, 1.0, 99)
mega_u_obs = mega_u_pivot / (1.0 + redshift)
hsc_g_obs = hsc_g_pivot / (1.0 + redshift)
hsc_r_obs = hsc_r_pivot / (1.0 + redshift)
hsc_i_obs = hsc_i_pivot / (1.0 + redshift)
hsc_z_obs = hsc_z_pivot / (1.0 + redshift)
hsc_y_obs = hsc_y_pivot / (1.0 + redshift)
wfc3_f125w_obs = wfc3_f125w_pivot / (1.0 + redshift)
wfc3_f160w_obs = wfc3_f160w_pivot / (1.0 + redshift)
sdss_u_obs = sdss_u / (1.0 + redshift)
sdss_g_obs = sdss_g / (1.0 + redshift)
sdss_r_obs = sdss_r / (1.0 + redshift)
sdss_i_obs = sdss_i / (1.0 + redshift)
sdss_z_obs = sdss_z / (1.0 + redshift)
In [9]:
z1_range = np.linspace(0.20, 0.35, 59)
z2_range = np.linspace(0.35, 0.50, 59)
z3_range = np.linspace(0.4, 0.5, 59)
z4_range = np.linspace(0.5, 0.6, 59)
In [6]:
fig = plt.figure(figsize=(14, 12))
fig.subplots_adjust(left=0.1, right=0.95, bottom=0.15)
ax1 = fig.add_subplot(111)
ax1.plot(redshift, mega_u_obs, 'b', linestyle='--', linewidth=3.5, alpha=0.8, label='MegaCam-u')
ax1.plot(redshift, hsc_g_obs, 'g', linestyle='-', linewidth=3.5, alpha=0.8, label='HSC-g')
ax1.plot(redshift, hsc_r_obs, 'y', linestyle='-', linewidth=3.5, alpha=0.8, label='HSC-r')
ax1.plot(redshift, hsc_i_obs, 'r', linestyle='-', linewidth=3.5, alpha=0.8, label='HSC-i')
ax1.plot(redshift, hsc_z_obs, 'm', linestyle='-', linewidth=3.5, alpha=0.8, label='HSC-z')
ax1.plot(redshift, hsc_y_obs, 'k', linestyle='-', linewidth=3.5, alpha=0.8, label='HSC-y')
ax1.plot(redshift, wfc3_f125w_obs, 'r', linestyle='-.', linewidth=3.5, alpha=0.6, label='WFC3/F125W')
ax1.plot(redshift, wfc3_f160w_obs, 'k', linestyle='-.', linewidth=3.5, alpha=0.6, label='WFC3/F160W')
ax1.axhline(sdss_g, color='g', linewidth=16.0, alpha=0.25)
ax1.axhline(sdss_r, color='y', linewidth=16.0, alpha=0.25)
ax1.axhline(sdss_i, color='r', linewidth=16.0, alpha=0.25)
ax1.set_xlim(0.0, 1.0)
ax1.set_ylim(1600.0, 15999.0)
ax1.set_xlabel('Redshift', size=28.0)
ax1.set_ylabel('Restframe Pivot Wavelength ($\AA$)', size=28.0)
for tick in ax1.xaxis.get_major_ticks():
tick.label.set_fontsize(20)
for tick in ax1.yaxis.get_major_ticks():
tick.label.set_fontsize(20)
ax1.spines['top'].set_linewidth(3.5)
ax1.spines['right'].set_linewidth(3.5)
ax1.spines['bottom'].set_linewidth(3.5)
ax1.spines['left'].set_linewidth(3.5)
ax1.legend(fontsize=20)
#ax1.text(-23.45, -23.0, '$0.3<z<0.4$', size=32.0)
Out[6]:
In [14]:
fig = plt.figure(figsize=(12, 12))
fig.subplots_adjust(left=0.14, right=0.95, bottom=0.12, top=0.95)
ax1 = fig.add_subplot(111)
ax1.plot(redshift, mega_u_obs, 'b', linestyle='--', linewidth=4.0, alpha=0.9, label='MegaCam-u')
ax1.plot(redshift, hsc_g_obs, 'g', linestyle='-', linewidth=4.0, alpha=0.9, label='HSC-g')
ax1.plot(redshift, hsc_r_obs, 'y', linestyle='-', linewidth=4.0, alpha=0.9, label='HSC-r')
ax1.plot(redshift, hsc_i_obs, 'r', linestyle='-', linewidth=4.0, alpha=0.9, label='HSC-i')
ax1.plot(redshift, hsc_z_obs, 'm', linestyle='-', linewidth=4.0, alpha=0.9, label='HSC-z')
ax1.plot(redshift, hsc_y_obs, 'k', linestyle='--', linewidth=4.0, alpha=0.9, label='HSC-y')
# Choice
ax1.plot(z1_range, hsc_i_pivot/(1.0+z1_range), 'r', linestyle='-', linewidth=10)
ax1.plot(z2_range, hsc_z_pivot/(1.0+z2_range), 'm', linestyle='-', linewidth=10)
#ax1.plot(z3_range, hsc_z_pivot/(1.0+z3_range), 'm', linestyle='-', linewidth=10)
ax1.axhline(sdss_g, color='g', linewidth=16.0, alpha=0.25)
ax1.axhline(sdss_r, color='y', linewidth=16.0, alpha=0.25)
ax1.axhline(sdss_i, color='r', linewidth=16.0, alpha=0.25)
ax1.set_xlim(0.09, 0.68)
ax1.set_ylim(1600.0, 9999.0)
ax1.set_xlabel('Redshift', size=28.0)
ax1.set_ylabel('Restframe Pivot Wavelength ($\AA$)', size=28.0)
for tick in ax1.xaxis.get_major_ticks():
tick.label.set_fontsize(20)
for tick in ax1.yaxis.get_major_ticks():
tick.label.set_fontsize(20)
ax1.spines['top'].set_linewidth(3.5)
ax1.spines['right'].set_linewidth(3.5)
ax1.spines['bottom'].set_linewidth(3.5)
ax1.spines['left'].set_linewidth(3.5)
ax1.legend(fontsize=20)
#ax1.text(-23.45, -23.0, '$0.3<z<0.4$', size=32.0)
fig.savefig('hsc_filter_redshift.png', dpi=200.0)
In [ ]: