I_leakage = float(v_gateleakage.query('READ?'))*1e3 #leakage current [nA]
V_array = np.array([])
resistance_array = np.array([])
I_leakage_array = np.array([])
measured = R_measure('dev352',
offset = 0,
amplitude= measure_amplitude,
out_channel = measure_output_channnel,
in_channel = measure_input_channnel,
time_constant = demodulation_time_constant,
frequency = measure_frequency,
poll_length = deamodulation_duration,
device=device,
daq=daq,
out_mixer_channel=out_mixer_channel,
bias_resistor=R_bias,
in_range = in_range,
out_range = out_range,
diff = False,
add = True,
ac = ac,
calibration_factor = calibration_factor)
Resistance = measured[0]
Voltage = Resistance *i
resistance_array = np.append(resistance_array,Resistance)
V_array = np.append(V_array,Voltage)
I_leakage_array = np.append(I_leakage_array,I_leakage)
if event.type == QUIT: sys.exit()
elif event.type == KEYDOWN and event.dict['key'] == 101:
END = True
if END:
break
measured = R_measure(device_id='dev352',
amplitude=measure_amplitude,
out_channel=measure_output_channnel,
in_channel=measure_input_channnel,
time_constant=demodulation_time_constant,
frequency=freq,
poll_length=deamodulation_duration,
device=device,
daq=daq,
out_mixer_channel=out_mixer_channel,
bias_resistor=bias_resistor,
in_range=in_range,
out_range=out_range,
diff=diff,
add=add,
offset=offset,
ac=ac)
Vin = np.abs(measured[4] + 1j * measured[5])
# r = Vin * bias_resistor/(measure_amplitude - Vin) # thie relation is valid when the phase is close to zero.
measured[0] = Vin
line = [count, gate_voltage, leakage_current] + measured
break
print('GATE: {:6.4f}'.format(gate_voltage), 'V')
time.sleep(time_step)
measured = R_measure(device_id = 'dev352',
amplitude = measure_amplitude,
out_channel = measure_output_channnel,
in_channel = measure_input_channnel,
time_constant = demodulation_time_constant,
frequency = measure_frequency,
poll_length = deamodulation_duration,
device = device,
daq = daq,
out_mixer_channel = out_mixer_channel,
bias_resistor = bias_resistor,
in_range = 4e-3,
out_range = 100e-3,
diff = False,
calibration_factor = 1.45,
add = False,
offset = 0,
ac = False)
measured[0] = calibration_factor * measured[0] + shift
line = [count,gate_voltage, leakage_current] + measured
if save_data:
stlab.writeline(my_file_2,line)
colnames=colnames)
plt.rcParams["figure.figsize"] = [16, 9]
fig, (ax1, ax2, ax3) = subplots(3, 1, sharex=True)
for cnt in range(3):
R_measure(device_id='dev352',
amplitude=measure_amplitude,
out_channel=measure_output_channnel,
in_channel=measure_input_channnel,
BW=BW,
filter_order=filter_order,
demod_rate=demod_rate,
frequency=measure_frequency,
poll_length=poll_length,
device=device,
daq=daq,
out_mixer_channel=out_mixer_channel,
bias_resistor=R_bias,
in_range=in_range,
out_range=out_range,
diff=diff,
add=add,
offset=offset,
ac=ac,
initialize=True)
start_time = time.time()
for Vg_count, Vg in enumerate(Vg_list):
B2961A.RampVoltage(Vg, tt=deltaVg / gate_ramp_speed,