tau = 1./(2*pi*1.7e6)
o.tau_in = tau/0.4
o.tau_ex = tau/0.6
o.n_th = 0
o.freq = linspace(0.99e6, 1.01e6, 400)
o.delta_hz = 0
figure('SB asymmetry on resonance', figsize=(10,13))
subplot(211)
o.n_th = 0
o.cooperativity = 0.1
title('Cooperativity = 0.1, n_th = 0')
plot(o.freq, o.right.output.spectrum_ss('upper'), '-r', label='Upper sideband')
plot(o.freq, o.right.output.spectrum_ss('lower'), '-b', label='Lower sideband')
xlabel('Frequency (Hz)')
ylabel('Single sided spectrum (Vacua)')
legend(loc='best')
subplot(212)
ax=gca()
ax.set_color_cycle([cm.bwr(k) for k in linspace(0,1,N_TEMPERATURES)])
coops = logspace(-2, 2, 100)
o.freq = array([1e6])
for n_th in [0.01*2**n for n in range(N_TEMPERATURES)]:
