spec_phot_2

import matplotlib.pyplot as plt
import numpy as np
from matplotlib.pyplot import figure
import matplotlib.pylab as pylab

#   y1 = LRIS
#   y2 = PFCam
#   x1, x2 = log(t) in seconds
y1 =  np.array([20.03,20.53,21.02,21.52,21.99,22.50,22.98,23.50,23.97,24.44,24.93,25.47,25.95,26.46,26.91,27.41])
x1 =  np.array([0.1513,0.3976,0.6180,0.8643,1.110,1.421,1.744,2.119,2.494,2.843,3.270,3.632,4.020,4.420,4.847,5.235])
F1 = 10**(-y1/2.5)*3640

y2 = np.array([17.98,18.48,19.00,19.49,19.98,20.49,21.00,21.53,21.99,22.46,22.98,23.48,23.97,24.42,24.96,25.44])
x2 = np.array([0.4280,0.6356,0.9077,1.256,1.593,1.942,2.291,2.705,3.131,3.493,3.907,4.308,4.709,5.109,5.497,5.885])
F2 = 10**(-y2/2.5)*3640

x1 = 10**x1; x2 = 10**x2    #    convert x1, x2 to seconds

figure(figsize=(15, 10), dpi=80)

params = {'legend.fontsize': 'xx-large',
         'axes.labelsize': 'xx-large',
         'axes.titlesize':'xx-large',
         'xtick.labelsize':'xx-large',
         'ytick.labelsize':'xx-large'}
pylab.rcParams.update(params)

ax1 = plt.subplot(1, 1, 1)
ax1.scatter(x1, F1, label = 'LRIS: m = -0.50', marker = '*', c = 'r', s = 50)
ax1.scatter(x2, F2, label = 'PFCam: m = -0.50', marker = 's', c = 'b', s = 50)

ax1.set_ylabel('Flux sensitivity, $\Delta S$ [Jy]')
ax1.set_xlabel('Integration time [s]')
ax1.set_xscale('log')
ax1.set_yscale('log')
ax1.grid(True)

ax2 = ax1.twiny()
ax3 = ax1.twinx()

hours = [0.03, 0.28, 2.78, 27.8, 278]
seconds = [1e2, 1e3, 1e4, 1e5, 1e6]
# new_pos = seconds
ax2.set_xscale('log')
ax2.set_xticks(seconds)
ax2.set_xticklabels(hours, rotation = 0)
ax2.set_xlabel('Integration time [hours]')
ax2.set_xlim(ax1.get_xlim())

WmHz   = [1e-30, 1e-31, 1e-32, 1e-33]
jansky = [1e-4, 1e-5, 1e-6, 1e-7]
ax3.set_yscale('log')
ax3.set_yticks(jansky)
ax3.set_yticklabels(WmHz, rotation = 0)
ax3.set_ylabel(r'Flux sensitivity [W m$^{-2}$ Hz$^{-1}$]')
ax3.set_ylim(ax1.get_xlim())