#Ttheo = np.linspace(T[0],T[-1], 1000)
Xtd, Ytd, Xt, Yt = fk.Curve(l_trans, l_ref, harmonic, phi,
a * max(abs(Z)), offset, T[0], T[-1], 1000)
delay = T
X = Z.real
X_s = Z_s.real
Y = Z.imag
Y_s = Z_s.imag
T_d, X_d, Y_d = fk.CutDataSet(T, X, Y, data_window)
Z = X_d + 1j * Y_d
scaling = 1E-4 #max(Z.real)/max(Ttheo)
# Fourier traffo
Zg = fk.GaussWindow(T_d, Z, suscept)
Td = T[1] - T[0]
if gauss:
wn, dft = fk.DFT(T_d,
Zg,
Td,
l_ref,
harmonic,
zeroPaddingFactor=zeroPaddingFactor)
else:
wn, dft = fk.DFT(T_d,
Z,
Td,
l_ref,
harmonic,
zeroPaddingFactor=zeroPaddingFactor)
delay = T
X = Z.real
X_s = Z_s.real
Y = Z.imag
Y_s = Z_s.imag
#print T
T_d = T[:]
X_d = X[:]
Y_d = Y[:]
T_d, X_d, Y_d = fk.CutDataSet(T, X, Y, data_window)
Z = X_d + 1j * Y_d
# Fourier traffo
Zg = fk.GaussWindow(T_d, Z, False)
Td = T[1] - T[0]
wn, dft = fk.DFT(T_d,
Zg,
Td,
l_ref,
harmonic,
zeroPaddingFactor=zeroPaddingFactor)
FWHM, peakFullWidth = fk.PeakWidth(T, 0.1, False,
zeroPaddingFactor)
FWHM = FWHM
#print(FWHM)
# Plot time domain
figTD = plt.figure('scan_{0}_{1}_d{2}_TD_{3}'.format(
run, dev, d, i),
Ttheo = np.linspace(delay[0],delay[-1], 30000)
Xtd,Ytd,Xt,Yt = fk.Curve(spc.h*spc.c/(spc.e*l_trans)*1e9, l_ref, harmonic, phi, a*max(abs(Z)), offset, delay[0], delay[-1], 30000) # theo curve for transition
Xtd,Ytd,Xt,Yt = fk.Curve(l_fel/harmonic, l_ref, harmonic, phi, a*max(abs(Z)), offset, delay[0], delay[-1], 30000) # theo curve for laser CWL
T_d,X_d,Y_d = fk.CutDataSet(delay, Z.real, Z.imag, data_window)
Z_cut = X_d + 1j*Y_d
Z_plot = Z_cut
Z_cut = fk.Phasing_TD(Z_cut,T_d,E_r) #time domain phasing of data set to take into account position of data window
scaling = 1 #max(Z.real)/max(Ttheo)
# Fourier traffo
Td = np.mean(np.unique(np.diff(delay)))
if gauss:
Zg = fk.GaussWindow(T_d, Z_cut, suscept)
wn, dft = fk.DFT(T_d, Zg, Td, l_ref , harmonic, zeroPaddingFactor = zeroPaddingFactor)
else:
wn, dft = fk.DFT(T_d, Z_cut, Td, l_ref , harmonic, zeroPaddingFactor = zeroPaddingFactor)
# dft = dft[::-1]
FWHM , peakFullWidth = fk.PeakWidth(T_d, 0.05, False, zeroPaddingFactor)
#Fouriert Trafo of the theoretical dataset
wn_theo, dft_theo = fk.DFT(t_fano_theo, I_fano_theo, t_fano_theo[1]-t_fano_theo[0], l_ref , harmonic, zeroPaddingFactor = zeroPaddingFactor)
''' Fitting the dephasing for 6th harmonic R'''
def decay(t,amp,tau):
return amp*np.exp(-t/tau) + 0.001828
def decay_theo(t,amp,tau):
# spectrogram
# =============================================================================
#analysing 5H
Z5 = Z[4]
Z5 *= np.exp(1j * np.pi / 180. * phas_corr) #phasing data
T_d, X_d, Y_d = fk.CutDataSet(T, Z5.real, Z5.imag, data_window)
Z5_cut = X_d + 1j * Y_d
Z5_cut = fk.Phasing_TD(Z5_cut, T_d, E_r - E_undersamp)
scaling = 1E-4 #max(Z.real)/max(Ttheo)
# Fourier traffo
Td = np.mean(np.unique(np.diff(T)))
if gauss:
Z5g = fk.GaussWindow(T_d, Z5_cut, suscept)
wn, dft = fk.DFT(T_d,
Z5g,
Td,
l_ref,
harmonic,
zeroPaddingFactor=zeroPaddingFactor)
else:
wn, dft = fk.DFT(T_d,
Z5_cut,
Td,
l_ref,
harmonic,
zeroPaddingFactor=zeroPaddingFactor)
FWHM, peakFullWidth = fk.PeakWidth(T_d, peakMargin, suscept, zeroPaddingFactor)