def get_trace(self, mode='lin'):
'''
Send trigger to device, wait for aquiring the data,
and read back the data from the device.
'''
qt.mstart()
startfreq = self.get_start_freq(query=False)
stopfreq = self.get_stop_freq(query=False)
numpoints = self.get_numpoints(query=False)
if mode == 'lin':
freqs = numpy.linspace(startfreq, stopfreq, numpoints)
elif mode == 'log':
freqs = numpy.logspace(numpy.log10(startfreq),
numpy.log10(stopfreq), numpoints)
else:
print 'mode needs to be either "lin" or "log"!'
return False
sweep_time = self.get_sweep_time(query=False)
print 'sending trigger to network analyzer, and wait to finish'
print 'estimated waiting time: %.2f s' % sweep_time
self.send_trigger()
qt.msleep(sweep_time)
print 'readout network analyzer'
reply = self.read()
reply = numpy.array(reply)
qt.mend()
return (freqs, reply)
def get_trace(self, mode='lin'):
'''
Send trigger to device, wait for aquiring the data,
and read back the data from the device.
'''
qt.mstart()
startfreq = self.get_start_freq(query=False)
stopfreq = self.get_stop_freq(query=False)
numpoints = self.get_numpoints(query=False)
if mode=='lin':
freqs = numpy.linspace(startfreq,stopfreq,numpoints)
elif mode=='log':
freqs = numpy.logspace(numpy.log10(startfreq),numpy.log10(stopfreq),numpoints)
else:
print 'mode needs to be either "lin" or "log"!'
return False
sweep_time = self.get_sweep_time(query=False)
print 'sending trigger to network analyzer, and wait to finish'
print 'estimated waiting time: %.2f s' % sweep_time
self.send_trigger()
qt.msleep(sweep_time)
print 'readout network analyzer'
reply = self.read()
reply = numpy.array(reply)
qt.mend()
return (freqs, reply)
def run(self):
qt.mstart()
self.Fire()
qt.msleep(0.1)
while not self.buttons[8]:
print "=============================="
self.Fire()
qt.msleep(0.1)
qt.mend()
print "So long, and thanks for all the fish!"
def run(self):
qt.mstart()
self.Fire()
qt.msleep(0.1)
while not self.buttons[8]:
print "=============================="
self.Fire()
qt.msleep(0.1)
qt.mend()
print "So long, and thanks for all the fish!"
def get_trace(self):
'''
This function performs a full measurement.
First a trigger is sent to initiate a sweep.
An estimate is made of the time the sweep takes.
After the estimated time the data is queried from
the device. Usually the estimated time is a bit
lower then the actual time, the device will
respond as soon as it is finished.
It is assumed that the instrument is already on
'trigger hold'-mode.
Input:
None
Ouptput:
freqs (array of floats): The frequencies at which the
reflection / transmission was
measured
reply (arrray of foats): Measured data
'''
qt.mstart()
startfreq = self.get_start_freq(query=False)
stopfreq = self.get_stop_freq(query=False)
numpoints = self.get_numpoints(query=False)
IF_Bandwidth = self.get_IF_Bandwidth(query=False)
freqs = numpy.linspace(startfreq,stopfreq,numpoints)
sweep_time = numpoints / IF_Bandwidth
print 'sending trigger to network analyzer, and wait to finish'
print 'estimated waiting time: %.2f s' % sweep_time
self.send_trigger()
qt.msleep(sweep_time)
print 'reading out network analyzer'
reply = self.read()
reply = numpy.array(reply)
qt.mend()
return (freqs, reply)
def get_trace(self):
'''
This function performs a full measurement.
First a trigger is sent to initiate a sweep.
An estimate is made of the time the sweep takes.
After the estimated time the data is queried from
the device. Usually the estimated time is a bit
lower then the actual time, the device will
respond as soon as it is finished.
It is assumed that the instrument is already on
'trigger hold'-mode.
Input:
None
Ouptput:
freqs (array of floats): The frequencies at which the
reflection / transmission was
measured
reply (arrray of foats): Measured data
'''
qt.mstart()
startfreq = self.get_start_freq(query=False)
stopfreq = self.get_stop_freq(query=False)
numpoints = self.get_numpoints(query=False)
IF_Bandwidth = self.get_IF_Bandwidth(query=False)
freqs = numpy.linspace(startfreq, stopfreq, numpoints)
sweep_time = numpoints / IF_Bandwidth
print 'sending trigger to network analyzer, and wait to finish'
print 'estimated waiting time: %.2f s' % sweep_time
self.send_trigger()
qt.msleep(sweep_time)
print 'reading out network analyzer'
reply = self.read()
reply = numpy.array(reply)
qt.mend()
return (freqs, reply)
def example1(f_vec, b_vec):
'''
this example is exactly the same as 'basic_measure_script.py'
but now made into a function. The main advantage is that now
the parameters f_vec and b_vec can be provided when calling
the function: "measure_module.example1(vec1, vec2)", instead
of having to change the script.
To run the function type in the terminal:
fv=numpy.arange(0,10,0.01)
bv=numpy.arange(-5,5,0.1)
measure_module.example1(fv, bv)
'''
qt.mstart()
data = qt.Data(name='testmeasurement')
data.add_coordinate('frequency, mw src 1 [Hz]')
data.add_coordinate('Bfield, ivvi dac 3 [mV]')
data.add_value('Psw SQUID')
data.create_file()
plot2d = qt.Plot2D(data, name='measure2D')
plot3d = qt.Plot3D(data, name='measure3D', style='image')
for b in b_vec:
fake_ivvi_set_dac_3(b)
for f in f_vec:
fake_mw_src_set_freq(f)
result = fake_readout_psw()
data.add_data_point(f, b, result)
qt.msleep(0.01)
data.new_block()
data.close_file()
qt.mend()
def example3(x_vec=numpy.linspace(0, 10, 10), y_vec=numpy.linspace(0, 10, 50)):
'''
To run the function type in the terminal:
measure_module.example3()
'''
qt.mstart()
data = qt.Data(name='testmeasurement')
data.add_coordinate('x')
data.add_coordinate('y')
data.add_value('z1')
data.add_value('z2')
data.add_value('z3')
data.create_file()
plot2d_1 = qt.Plot2D(data, name='2D_1', coorddim=1, valdim=2)
plot2d_2 = qt.Plot2D(data, name='2D_2', coorddim=1, valdim=2, maxtraces=1)
plot2d_3 = qt.Plot2D(data, name='2D_3', coorddim=1, valdim=2, maxtraces=1)
plot2d_3.add_data(data, coorddim=1, valdim=3, maxtraces=1)
plot2d_3.add_data(data, coorddim=1, valdim=4, maxtraces=1)
plot2d_4 = qt.Plot2D(data, name='2D_4', coorddim=1, valdim=2, mintime=0.3)
plot2d_5 = qt.Plot2D(data,
name='2D_5',
coorddim=1,
valdim=2,
autoupdate=False)
plot3d_1 = qt.Plot3D(data, name='3D_1', style='image')
plot3d_2 = qt.Plot3D(data,
name='3D_2',
style='image',
coorddims=(1, 0),
valdim=4)
for x in x_vec:
for y in y_vec:
z1 = numpy.sin(x + y)
z2 = numpy.cos(x + y)
z3 = numpy.sin(x + 2 * y)
data.add_data_point(x, y, z1, z2, z3)
if z1 > 0:
plot2d_5.update()
qt.msleep(0.1)
data.new_block()
plot2d_1.save_png()
plot2d_1.save_gp()
plot3d_2.save_png()
plot3d_2.save_gp()
data.close_file()
qt.mend()
def example2(f_vec, b_vec):
'''
This example introduces three new features:
1) setting format and/or precision of data in the datafile.
using 'precision' will keep the default scientific notation,
'format' can be anything you like
=> add_coordinate(precision=<nr>)
=> add_coordinate(format='<format_string>')
2) specify specific filepath for the data file (in stead of
automatic filepath)
=> create_file(filepath=<filepath>)
3) turn off automatic saving of instrument-settings-file.
=> create_file(settings_file=False)
To run the function type in the terminal:
fv=numpy.arange(0,10,0.01)
bv=numpy.arange(-5,5,0.1)
measure_module.example2(fv, bv)
'''
qt.mstart()
# this shows how to change format of saved data (per column)
data = qt.Data(name='testmeasurement')
data.add_coordinate('frequency, mw src 1 [Hz]', precision=3)
data.add_coordinate('Bfield, ivvi dac 3 [mV]', format='%.12f')
data.add_value('Psw SQUID', format='%.3e')
data.create_file()
# this shows how to save to a specific path and name, and how
# to avoid a settings file to be created. The directory is first
# retreived from the previous data object
dir = data.get_dir()
maxfilepath = os.path.join(dir, 'maxvals.dat')
data_max = qt.Data(name='maxvals')
data_max.add_coordinate('Bfield, ivvi dac 3 [mV]')
data_max.add_value('resonance frequency [Hz]')
data_max.create_file(filepath=maxfilepath, settings_file=False)
plot2d = qt.Plot2D(data, name='measure2D')
plot3d = qt.Plot3D(data, name='measure3D', style='image')
plot2dmax = qt.Plot2D(data_max, name='maxvals')
for b in b_vec:
fake_ivvi_set_dac_3(b)
last_trace = []
for f in f_vec:
fake_mw_src_set_freq(f)
result = fake_readout_psw()
data.add_data_point(f, b, result)
last_trace.append(result)
qt.msleep(0.01)
data.new_block()
loc_of_max = numpy.argmax(last_trace)
freq_at_max = f_vec[loc_of_max]
data_max.add_data_point(b, freq_at_max)
data.close_file()
data_max.close_file()
qt.mend()
import numpy as np
from qtlab.source import qt
d = qt.Data()
d.add_coordinate('X')
d.add_value('Y')
d.create_file()
p = plot(d)
# Inform wrapper that a measurement has started
qt.mstart()
for x in arange(0, 40, 0.1):
y = np.sin(x) + np.random.rand() / 10
d.add_data_point(x, y)
# Sleep for 100msec and allow UI interaction
qt.msleep(0.1)
# Inform wrapper that the measurement has ended
qt.mend()
d.close_file()
import numpy as np
from qtlab.source import qt
d = qt.Data()
d.add_coordinate('X')
d.add_value('Y')
d.create_file()
p = plot(d)
# Inform wrapper that a measurement has started
qt.mstart()
for x in arange(0, 40, 0.1):
y = np.sin(x) + np.random.rand()/10
d.add_data_point(x, y)
# Sleep for 100msec and allow UI interaction
qt.msleep(0.1)
# Inform wrapper that the measurement has ended
qt.mend()
d.close_file()