cav.gamma_m_hz = 2560.
tau = 1./(2*pi*1.7e6)
cav.tau_input = tau/0.4
cav.tau_output = tau/0.6
cav.n_photons = 1.1e8
cav.freq = linspace(1.5e6, 1.56e6, 1000)
figure('mean_fields_omr', figsize=(12.9, 10.9))
last=None
for index, ((input_trans, out_trans, loss), title_) in enumerate((((1e-6,5e-6, 0), 'undercoupled'),
((5e-6, 5e-6,0), 'critically coupled'),
((5e-6, 0.5e-6,0), 'over coupled'),
((1e-6, 1e-6, 5e-6), 'loss dominated'))):
last = subplot(2, 2 ,index + 1, sharex=last)
cav.transmission_input = input_trans
cav.transmission_output = out_trans
cav.losses = loss
title(title_)
reflected = []
transmitted = []
deltas = linspace(-25*cav.bandwidth_hz, 25*cav.bandwidth_hz, 300)
for delta_hz in deltas:
cav.delta_hz = delta_hz
reflected.append(cav.left.output.mean_field)#/cav.intra.mean_field))
transmitted.append(cav.right.output.mean_field)#/cav.intra.mean_field))
reflected = array(reflected)
transmitted = array(transmitted)
plot(real(reflected), imag(reflected),'o', label='reflected')
plot(real(transmitted), imag(transmitted),'o', label='transmitted')
xlabel('Real')
