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) for index, delta_hz in enumerate(dets): subplot(len(dets),2,2*index+1)