asyms.append(asym)
plot(coops, asyms, label='n_th=' + str(n_th))
gca().set_xscale('log')
xlabel('Cooperativity')
ylabel('Sideband asymmetry')
legend(loc='best')
show()
savefig(current_dir + '/figures/SB asymmetry on resonance.png')
savefig(current_dir + '/figures/SB asymmetry on resonance.pdf')
figure('SB asymmetry Painter style',figsize=(10,13))
subplot(211)
o.rsf = 10
o.n_th = 0
o.cooperativity = 0.01
title('Cooperativity = 0.01, n_th = 0')
o.freq = linspace(o.f_mech - 3000, o.f_mech+3000, 200)
o.delta_hz = o.f_mech
plot(o.freq, o.right.output.spectrum_sym(pi/2), '-b', label='drive on blue side (lower sb. resonant)')
o.delta_hz = -o.f_mech
plot(o.freq, o.right.output.spectrum_sym(pi/2), '-r', label='drive on red side (upper sb. resonant)')
xlabel('Frequency (Hz)')
ylabel('Phase quadrature spectrum (Vacua)')
legend(loc='best')
subplot(212)
ax=gca()
