This page was generated from
docs/examples/DataSet/Paramtypes explained.ipynb.
Interactive online version:
.
Paramtypes explained¶
Internally in the SQLite database on disk, data are registered as being of one of (currently) four allowed types:
numericarraytextcomplex
This notebook seeks to exemplify when each type should be used, and how differently the Measurement object treats data of each type.
We start with necessary imports, and then initialising our database and creating an experiment.
[1]:
import os
import time
import numpy as np
from qcodes.dataset import (
Measurement,
initialise_or_create_database_at,
load_or_create_experiment,
)
from qcodes.instrument_drivers.mock_instruments import DummyInstrument
from qcodes.parameters import ArrayParameter
Logging hadn't been started.
Activating auto-logging. Current session state plus future input saved.
Filename : /home/runner/.qcodes/logs/command_history.log
Mode : append
Output logging : True
Raw input log : False
Timestamping : True
State : active
Qcodes Logfile : /home/runner/.qcodes/logs/241218-15940-qcodes.log
[2]:
initialise_or_create_database_at(os.path.join(os.getcwd(), "paramtypes_explained.db"))
exp = load_or_create_experiment("paramtypes", sample_name="not_available")
Let us, now, create two dummy instruments to be used in our experiment.
[3]:
dac = DummyInstrument("dac", gates=["ch1", "ch2"])
SA = DummyInstrument("SA")
[4]:
# some array-like data types
class Spectrum(ArrayParameter):
def __init__(self, name, instrument):
self.N = 7
setpoints = (np.linspace(0, 1, self.N),)
super().__init__(
name=name,
instrument=instrument,
setpoints=setpoints,
shape=(20,),
label="Noisy spectrum",
unit="V/sqrt(Hz)",
setpoint_names=("Frequency",),
setpoint_units=("Hz",),
)
def get_raw(self):
return np.random.randn(self.N)
class MultiDimSpectrum(ArrayParameter):
def __init__(self, name, instrument):
self.start = 0
self.stop = 1
self.npts = (2, 5, 3)
sp1 = np.linspace(self.start, self.stop, self.npts[0])
sp2 = np.linspace(self.start, self.stop, self.npts[1])
sp3 = np.linspace(self.start, self.stop, self.npts[2])
setpoints = (
sp1,
np.tile(sp2, (len(sp1), 1)),
np.tile(sp3, (len(sp1), len(sp2), 1)),
)
super().__init__(
name=name,
instrument=instrument,
setpoints=setpoints,
shape=(100, 50, 20),
label="Flower Power Spectrum in 3D",
unit="V/sqrt(Hz)",
setpoint_names=("Frequency0", "Frequency1", "Frequency2"),
setpoint_units=("Hz", "Other Hz", "Third Hz"),
)
def get_raw(self):
a = self.npts[0]
b = self.npts[1]
c = self.npts[2]
return np.reshape(np.arange(a * b * c), (a, b, c))
# a string-valued parameter
def dac1_too_high():
return "Too high" if dac.ch1() > 5 else "OK"
Finally, we add our parameters to the dummy instruments:
[5]:
dac.add_parameter("control", get_cmd=dac1_too_high)
SA.add_parameter("spectrum", parameter_class=Spectrum)
SA.add_parameter("spectrum3D", parameter_class=MultiDimSpectrum)
2024-12-18 07:06:32,482 ¦ py.warnings ¦ WARNING ¦ warnings ¦ _showwarnmsg ¦ 110 ¦ /opt/hostedtoolcache/Python/3.11.11/x64/lib/python3.11/site-packages/qcodes/instrument/instrument_base.py:173: QCoDeSDeprecationWarning: Parameter spectrum on instrument SA does not correctly pass kwargs to its baseclass. A Parameter class must take `**kwargs` and forward them to its baseclass.
warnings.warn(
2024-12-18 07:06:32,483 ¦ py.warnings ¦ WARNING ¦ warnings ¦ _showwarnmsg ¦ 110 ¦ /opt/hostedtoolcache/Python/3.11.11/x64/lib/python3.11/site-packages/qcodes/instrument/instrument_base.py:173: QCoDeSDeprecationWarning: Parameter spectrum3D on instrument SA does not correctly pass kwargs to its baseclass. A Parameter class must take `**kwargs` and forward them to its baseclass.
warnings.warn(
[5]:
<__main__.MultiDimSpectrum: spectrum3D at 139863537846160>
Numeric¶
The numeric datatype is simply a number. Data registered with this type are saved as individual numbers. This is the default datatype when registering parameters.
Numeric example 1¶
In this example, all parameters get registered as numeric type. This entails that the array in unraveled and inserted point-by-point.
[6]:
meas = Measurement(exp=exp)
meas.register_parameter(dac.ch1)
meas.register_parameter(SA.spectrum, setpoints=(dac.ch1,))
t0 = time.perf_counter()
with meas.run() as datasaver:
for dac_v in np.linspace(0, 2, 5):
dac.ch1(dac_v)
datasaver.add_result((dac.ch1, dac_v), (SA.spectrum, SA.spectrum()))
t1 = time.perf_counter()
print(f"Finished run in {(t1-t0):.3f} s")
dataset1 = datasaver.dataset
Starting experimental run with id: 1.
Finished run in 0.015 s
The data may be retrieved using the get_parameter_data method. This function will bring back the data in a way that reflects the datastructure as it is stored.
[7]:
dataset1.get_parameter_data()
[7]:
{'SA_spectrum': {'SA_spectrum': array([[-0.6814824 , -0.74946626, -0.38597474, -1.02153739, -2.04081094,
-0.09279146, -0.6864718 ],
[-0.4216081 , -1.01521131, 1.21051307, 0.6338473 , -0.55575643,
0.61161355, -0.50125294],
[ 1.57947322, 0.2580065 , -0.98999669, -0.40029859, -0.61697633,
1.10295331, -0.89180732],
[-1.72429486, 1.40686508, 1.38774825, -0.13346496, -0.61188738,
-1.13295855, -0.05560206],
[-2.05669888, 0.96469791, 0.86816617, -0.20606045, 0.51835868,
-0.18693992, 0.46001545]]),
'dac_ch1': array([[0. , 0. , 0. , 0. , 0. , 0. , 0. ],
[0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5],
[1. , 1. , 1. , 1. , 1. , 1. , 1. ],
[1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5],
[2. , 2. , 2. , 2. , 2. , 2. , 2. ]]),
'SA_Frequency': array([[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ]])}}
Array¶
The array paramtype stores data as binary blobs in the database. Insertion is faster (asymptotically much faster) this way, but the data are “dead” to SQL queries inside the database. Be informed that a BLOB in sqlite has a default max length limit set at 1 billion (1,000,000,000) bytes (for more information, refer to Sqlite docs).
Array example 1¶
Let us repeat the above measurement, but this time using array paramtypes.
[8]:
meas = Measurement(exp=exp)
meas.register_parameter(dac.ch1)
meas.register_parameter(SA.spectrum, setpoints=(dac.ch1,), paramtype="array")
t0 = time.perf_counter()
with meas.run() as datasaver:
for dac_v in np.linspace(0, 2, 5):
dac.ch1(dac_v)
datasaver.add_result((dac.ch1, dac_v), (SA.spectrum, SA.spectrum()))
t1 = time.perf_counter()
print(f"Finished run in {(t1-t0):.3f} s")
dataset2 = datasaver.dataset
Starting experimental run with id: 2.
Finished run in 0.011 s
[9]:
dataset2.get_parameter_data()
[9]:
{'SA_spectrum': {'SA_spectrum': array([[-7.60597039e-01, -3.09921226e-02, -1.70332171e-01,
-8.96008831e-01, 5.19751412e-03, 1.14986085e+00,
-7.50423739e-01],
[ 4.43093487e-01, -1.01785899e+00, 1.54686958e+00,
8.80969746e-01, -4.60650193e-01, 2.17418621e+00,
8.10620681e-02],
[-3.02749971e-01, -8.34551275e-01, 1.03108295e-03,
-1.20783791e+00, 4.14869095e-01, -2.16585717e+00,
-1.04032049e+00],
[ 1.46982837e+00, -1.03979531e+00, 5.10961311e-01,
1.72782342e+00, 1.56380323e+00, 2.55651901e+00,
1.14032940e+00],
[-6.72411169e-01, 6.46030189e-01, 6.60892574e-01,
-1.28347560e+00, 2.24078612e-01, -2.75634847e-01,
1.16169994e+00]]),
'dac_ch1': array([[0. , 0. , 0. , 0. , 0. , 0. , 0. ],
[0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5],
[1. , 1. , 1. , 1. , 1. , 1. , 1. ],
[1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5],
[2. , 2. , 2. , 2. , 2. , 2. , 2. ]]),
'SA_Frequency': array([[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ],
[0. , 0.16666667, 0.33333333, 0.5 , 0.66666667,
0.83333333, 1. ]])}}
Array example 2¶
When storing multidimensional array data (think: Alazar cards), both numeric and array types can be used.
[10]:
meas = Measurement(exp=exp)
meas.register_parameter(SA.spectrum3D, paramtype="array")
with meas.run() as datasaver:
datasaver.add_result((SA.spectrum3D, SA.spectrum3D()))
dataset3 = datasaver.dataset
Starting experimental run with id: 3.
The data come out the way we expect them to.
[11]:
dataset3.get_parameter_data()
[11]:
{'SA_spectrum3D': {'SA_spectrum3D': array([[[[ 0, 1, 2],
[ 3, 4, 5],
[ 6, 7, 8],
[ 9, 10, 11],
[12, 13, 14]],
[[15, 16, 17],
[18, 19, 20],
[21, 22, 23],
[24, 25, 26],
[27, 28, 29]]]]),
'SA_Frequency0': array([[[[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.]],
[[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.]]]]),
'SA_Frequency1': array([[[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]],
[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]]]]),
'SA_Frequency2': array([[[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]],
[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]]]])}}
Array example 3¶
For completeness, here, we provide an example where the multidimensional array has an auxiliary setpoint.
[12]:
meas = Measurement(exp=exp)
meas.register_parameter(dac.ch1)
meas.register_parameter(SA.spectrum3D, paramtype="array", setpoints=(dac.ch1,))
with meas.run() as datasaver:
for dac_v in [3, 4, 5]:
dac.ch1(dac_v)
datasaver.add_result((dac.ch1, dac_v), (SA.spectrum3D, SA.spectrum3D()))
dataset4 = datasaver.dataset
Starting experimental run with id: 4.
[13]:
dataset4.get_parameter_data()
[13]:
{'SA_spectrum3D': {'SA_spectrum3D': array([[[[ 0, 1, 2],
[ 3, 4, 5],
[ 6, 7, 8],
[ 9, 10, 11],
[12, 13, 14]],
[[15, 16, 17],
[18, 19, 20],
[21, 22, 23],
[24, 25, 26],
[27, 28, 29]]],
[[[ 0, 1, 2],
[ 3, 4, 5],
[ 6, 7, 8],
[ 9, 10, 11],
[12, 13, 14]],
[[15, 16, 17],
[18, 19, 20],
[21, 22, 23],
[24, 25, 26],
[27, 28, 29]]],
[[[ 0, 1, 2],
[ 3, 4, 5],
[ 6, 7, 8],
[ 9, 10, 11],
[12, 13, 14]],
[[15, 16, 17],
[18, 19, 20],
[21, 22, 23],
[24, 25, 26],
[27, 28, 29]]]]),
'dac_ch1': array([[[[3., 3., 3.],
[3., 3., 3.],
[3., 3., 3.],
[3., 3., 3.],
[3., 3., 3.]],
[[3., 3., 3.],
[3., 3., 3.],
[3., 3., 3.],
[3., 3., 3.],
[3., 3., 3.]]],
[[[4., 4., 4.],
[4., 4., 4.],
[4., 4., 4.],
[4., 4., 4.],
[4., 4., 4.]],
[[4., 4., 4.],
[4., 4., 4.],
[4., 4., 4.],
[4., 4., 4.],
[4., 4., 4.]]],
[[[5., 5., 5.],
[5., 5., 5.],
[5., 5., 5.],
[5., 5., 5.],
[5., 5., 5.]],
[[5., 5., 5.],
[5., 5., 5.],
[5., 5., 5.],
[5., 5., 5.],
[5., 5., 5.]]]]),
'SA_Frequency0': array([[[[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.]],
[[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.]]],
[[[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.]],
[[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.]]],
[[[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.]],
[[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.],
[1., 1., 1.]]]]),
'SA_Frequency1': array([[[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]],
[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]]],
[[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]],
[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]]],
[[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]],
[[0. , 0. , 0. ],
[0.25, 0.25, 0.25],
[0.5 , 0.5 , 0.5 ],
[0.75, 0.75, 0.75],
[1. , 1. , 1. ]]]]),
'SA_Frequency2': array([[[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]],
[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]]],
[[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]],
[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]]],
[[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]],
[[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ],
[0. , 0.5, 1. ]]]])}}
Text¶
Text is strings. Sometimes it may be useful to capture categorial data that is represented as string values, or a log message, or else.
Note that the paramtype setting is important. The datasaver will not allow to save numeric data for a parameter that was registered as text. The opposite it also true: the datasaver will not allow to save strings for a parameter what was registered as non-text (numeric or array).
[14]:
meas = Measurement(exp=exp)
meas.register_parameter(dac.ch1)
meas.register_parameter(dac.control, setpoints=(dac.ch1,), paramtype="text")
with meas.run() as datasaver:
for dac_v in np.linspace(4, 6, 10):
dac.ch1(dac_v)
datasaver.add_result((dac.ch1, dac_v), (dac.control, dac.control()))
dataset5 = datasaver.dataset
Starting experimental run with id: 5.
[15]:
dataset5.get_parameter_data()
[15]:
{'dac_control': {'dac_control': array(['OK', 'OK', 'OK', 'OK', 'OK', 'Too high', 'Too high', 'Too high',
'Too high', 'Too high'], dtype='<U8'),
'dac_ch1': array([4. , 4.22222222, 4.44444444, 4.66666667, 4.88888889,
5.11111111, 5.33333333, 5.55555556, 5.77777778, 6. ])}}
[ ]: