In [2]:

    

import numpy as np 
import matplotlib.pyplot as plt
import glob
plt.ion()
# Loads nicard and scope
manager.startModule('logic','cavitylogic')
filenames = glob.glob('./*data.dat')

print(filenames)
outlier_cutoff = 1.5


    

    




[]

In [ ]:

    




    

In [2]:

    

fileNum=2
cavitylogic._load_full_sweep(, filename=filenames[fileNum])
cavitylogic._trim_data()
cavitylogic._data_split_up()
cavitylogic._polyfit_SG(order = 3,plot=False)
resonances = cavitylogic._peak_search()
corrected_resonances = cavitylogic._find_missing_resonances(resonances)


    

    




Traceback (most recent call last):
  File "C:\Users\ChildressLab\Software\Qudi\logic\jupyterkernel\qzmqkernel.py", line 795, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-2-c427b2e95357>", line 2, in <module>
    cavitylogic._load_full_sweep(filepath='', filename=filenames[fileNum])
IndexError: list index out of range

In [61]:

    

delta_peaks = cavitylogic._update_peaks_distance(corrected_resonances)
plt.plot(delta_peaks,'x-')
plt.show()
    
plt.plot(cavitylogic.RampUp_signalSG_polyfit, cavitylogic.RampUp_signalR)
plt.plot(cavitylogic.RampUp_signalSG_polyfit[corrected_resonances],  cavitylogic.RampUp_signalR[corrected_resonances],'r',marker ='o')
plt.xlabel('Absolute Voltage_from_NI_card (V)')
plt.ylabel('Reflection (mircon)')

plt.grid('on')
#plt.xlim(8.0,9.2)
plt.show()


    

In [62]:

    

cavitylogic.time_trim[-1]-cavitylogic.time_trim[0]


    

    
Out[62]:





1.5

In [63]:

    

cavitylogic._fit_ramp()


    

In [64]:

    

A = cavitylogic._ni.sweep_function(cavitylogic.RampUp_time[corrected_resonances]-0.5e-3,*cavitylogic.popt)
AD = cavitylogic.RampUp_signalSG[corrected_resonances]
B = cavitylogic.RampUp_signalNI[corrected_resonances]
C = cavitylogic.volts_trim[1]
D = cavitylogic.volts_trim[3]
E = cavitylogic.volts[3]
plt.plot(cavitylogic.time,E, '-*', linewidth =0.3)
plt.plot(cavitylogic.time_trim,C, '--', linewidth =0.3)
plt.plot(cavitylogic.time_trim,D, '--', linewidth =0.3)
plt.plot(cavitylogic.RampUp_time[corrected_resonances],AD,'o')
plt.plot(cavitylogic.RampUp_time[corrected_resonances],A,'*')
plt.plot(cavitylogic.RampUp_time[corrected_resonances],B,'x')
plt.grid()
plt.xlim(0.1,0.95)
plt.ylim(-3.5,1)


    

    
Out[64]:





(-3.5, 1)






    

In [ ]:

    

# Go to position
Rescale NI:
cavitylogic._ni.cavity_set_position(20e-6) # 0 V
cavitylogic.run_single()
cavitylogic._scope.aquire_data()
# Measure position