#### PRX Regal
o = OmR() #\epsilon_p*\epsilon_d
o.set_params(m=6.75e-12,
f_mech=1.5243e6,
q=100.,
lambda_nm = 1064.,
losses=1.e-6,
transmission_input=30e-6,
transmission_output=10.e-6,
length=0.00764,
i_incident=1e-3,
temp=3.8e-4,
delta_hz=1.,
detection_losses=0.304)
o.gamma_m_hz = 2560.
tau = 1./(2*pi*1.7e6)
o.tau_input = tau/0.4
o.tau_output = tau/0.6
o.n_photons = 1.1e8
o.freq = linspace(1.5e6, 1.56e6, 1000)
#o.right.output._get_losses = lambda:
if True:
figname = 'PRX Regal fig1'
close(figname)
figure(figname, figsize=(15,10))
dets = (3e3, 6e3, 13e3, 20e3)
cav = OmR() #\epsilon_p*\epsilon_d
cav.set_params(m=6.75e-12,
f_mech=1.5243e6,
q=100.,
lambda_nm = 1064.,
losses=1.e-6,
transmission_input=30e-6,
transmission_output=10.e-6,
length=0.00764,
i_incident=1e-3,
temp=3.8e-4,
delta_hz=1.,
detection_losses=0.304)
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