EVALUATION: show felix

This commit is contained in:
Valentin Boettcher 2024-05-21 17:03:55 -04:00
parent 393610e0cc
commit 645a92d9e6
2 changed files with 28 additions and 3 deletions

View file

@ -9,7 +9,7 @@ from plot_utils import wrap_plot
def transient_rabi(): def transient_rabi():
"""A transient rabi oscillation without noise.""" """A transient rabi oscillation without noise."""
params = Params(η=0.0001, δ=1 / 4, d=0.1, laser_detuning=0.01, Δ=0.005, N=2) params = Params(η=0.0001, δ=1 / 4, d=0.1, laser_detuning=0.1, Δ=0.005, N=2)
t = time_axis(params, 3, 0.1) t = time_axis(params, 3, 0.1)
solution = solve(t, params) solution = solve(t, params)
signal = output_signal(t, solution.y, params.laser_detuning) signal = output_signal(t, solution.y, params.laser_detuning)
@ -59,8 +59,7 @@ def ringdown_after_rabi():
"""Demonstrates the nonstationary ringdown of the resonator after turning off the EOM and laser drive.""" """Demonstrates the nonstationary ringdown of the resonator after turning off the EOM and laser drive."""
off_lifetime = 4 off_lifetime = 4
laser_detuning = 0.1 laser_detuning = 0.1
params = Params(η=0.0001, d=0.01, laser_detuning=laser_detuning, Δ=0.00, N=2)
params = Params(η=0.0001, d=0.01, laser_detuning=laser_detuning, Δ=0.00, N=4)
params.laser_off_time = params.lifetimes(off_lifetime) params.laser_off_time = params.lifetimes(off_lifetime)
params.drive_off_time = params.lifetimes(off_lifetime) params.drive_off_time = params.lifetimes(off_lifetime)

View file

@ -0,0 +1,26 @@
from ringfit import data
import matplotlib.pyplot as plt
from ringfit.data import *
from ringfit.plotting import *
from ringfit.fit import *
from rabifun.analysis import *
import numpy as np
# %% load data
path = "../../data/08_05_24/characterization_first"
scan = ScanData.from_dir(path)
# %% Fourier
freq, fft = fourier_transform(
scan.time, scan.output, low_cutoff=1000, high_cutoff=10**8
)
fig = plt.figure("interactive")
fig.clf()
ax, ax2 = fig.subplots(1, 2)
ax.set_yscale("log")
ax.plot(freq * 10 ** (-6), np.abs(fft) ** 2)
plot_scan(scan, ax=ax2, linewidth=0.1)
ax2.set_xlim(0.02, 0.020001)