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 = []
cooling = []
phis = linspace(-4,4,100)
arcizet.delta_hz = 0
arcizet.i_intra = power
for index, ph in enumerate(phis): #(0.03, 0.06, 0.09, 0.11, 0.13)):#
arcizet.delta_hz = ph*kappa_hz/2
y = arcizet.right.output.spectrum_sym(pi/2)#dummy_spec()
c = curve.Curve()
c.set_data(Series(abs(y), index=arcizet.freq))
f,a = c.fit('lorentz')
print f.params
shifts.append(f.params['x0'])#arcizet.freq[y.argmax()])
cooling.append(2*f.params['bandwidth'])#y.max())
#gca().set_yscale('log')
subplot(211)
plot(phis, array(shifts) - arcizet.omega_m/(2*pi),label=str(power) + ' W')
xlabel('phi')
ylabel('freq shift (Hz)')
subplot(212)
