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)]: