from pylab import *
import os.path as osp
from matplotlib.widgets import Slider, Button, RadioButtons
from omr import OmR
current_dir = osp.join(osp.dirname(__file__))
if current_dir=='':
current_dir = '.'
#close("pilier")
#figure("pilier")
ax=gca()
N_DET = 3
ax.set_color_cycle([cm.jet(k) for k in linspace(0,1,N_DET)])
pilier = OmR()#0.304)
pilier.set_params(m=25.e-9,
f_mech=4.e6,
q=2e6,
lambda_nm = 1064.,
losses=15.e-6,
transmission_input=30e-6,
transmission_output=1.e-6,
length=150e-6,
i_incident=1000e-6,
temp=1.,
delta_hz=10,
detection_losses=0.0)
pilier.freq = linspace(3.99e6,4.01e6, 1000)
"""
from pylab import *
import os.path as osp
from omr import OmR
current_dir = osp.join(osp.dirname(__file__))
if current_dir=='':
current_dir = '.'
#### 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)
from pylab import *
import os.path as osp
from omr import OmR
current_dir = osp.join(osp.dirname(__file__))
if current_dir=='':
current_dir = '.'
N_TEMPERATURES = 10
o = OmR()
o.set_params(m=6.75e-12,
f_mech=1e6,
q=1000.,
lambda_nm = 1064.,
losses=1.e-6,
transmission_input=300e-6,
length=0.00764,
i_incident=1e-3,
delta_hz=1.)
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
from pylab import *
import os.path as osp
from omr import OmR
current_dir = osp.join(osp.dirname(__file__))
if current_dir=='':
current_dir = '.'
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'),
from pylab import *
import os.path as osp
from omr import OmR
import curve
from pandas import Series
current_dir = osp.join(osp.dirname(__file__))
if current_dir=='':
current_dir = '.'
### Arcizet 2006
arcizet = OmR()
arcizet.set_params(m=190e-9,
f_mech=814e3,
q=1.e4,
lambda_nm = 1064.,
losses=1.e-6,
transmission_input=2*pi/3e4,
transmission_output=1e-6,
length=0.0024,
i_incident=2e-4,
temp=300,
delta_hz=1.)
arcizet.freq=linspace(813.e3, 816.e3, 1000)
kappa_hz = arcizet.kappa_hz
name = "Arcizet2006"
figure(name)
for power in array((1.4,2.5,4.5,6.5,9.5)):
shifts = []
