2015-09-07 06:08:57 +08:00
|
|
|
# Copyright (C) 2014, 2015 Robert Jordens <jordens@gmail.com>
|
|
|
|
|
2014-10-04 09:55:53 +08:00
|
|
|
import numpy as np
|
|
|
|
|
2016-01-26 05:48:58 +08:00
|
|
|
from artiq.language import *
|
2014-10-04 09:55:53 +08:00
|
|
|
|
2016-01-20 11:09:10 +08:00
|
|
|
from artiq.wavesynth.coefficients import SplineSource
|
2014-10-04 09:55:53 +08:00
|
|
|
|
2016-01-20 11:09:10 +08:00
|
|
|
transport = SplineSource(
|
|
|
|
x=np.linspace(0, 10, 101), # waveform time
|
|
|
|
y=np.random.rand(4*3*3, 101)*1e-6, # waveform data,
|
2014-12-03 18:20:30 +08:00
|
|
|
# 4 devices, 3 board each, 3 dacs each
|
|
|
|
)
|
|
|
|
|
2015-07-14 04:08:20 +08:00
|
|
|
class Transport(EnvExperiment):
|
2015-03-08 22:43:04 +08:00
|
|
|
"""Transport"""
|
2015-02-21 05:01:34 +08:00
|
|
|
|
2015-07-14 04:08:20 +08:00
|
|
|
def build(self):
|
2016-01-20 11:09:10 +08:00
|
|
|
self.core = self.get_device("core")
|
|
|
|
self.bd_sw = self.get_device("bd_sw")
|
|
|
|
self.pmt = self.get_device("pmt")
|
|
|
|
self.electrodes = self.get_device("electrodes")
|
2015-10-04 00:18:21 +08:00
|
|
|
|
2016-01-20 11:09:10 +08:00
|
|
|
self.wait_at_stop = self.get_argument("wait_at_stop",
|
|
|
|
NumberValue(100*us))
|
|
|
|
self.speed = self.get_argument("speed", NumberValue(1.5))
|
|
|
|
self.repeats = self.get_argument("repeats", NumberValue(100))
|
|
|
|
self.nbins = self.get_argument("nbins", NumberValue(100))
|
2014-10-19 13:09:10 +08:00
|
|
|
|
2015-07-14 04:08:20 +08:00
|
|
|
def calc_waveforms(self, stop):
|
2015-03-23 05:31:59 +08:00
|
|
|
self.electrodes.disarm()
|
2014-10-17 00:04:30 +08:00
|
|
|
self.tf = self.electrodes.create_frame()
|
2016-01-20 11:09:10 +08:00
|
|
|
to_stop = self.tf.create_segment("to_stop")
|
|
|
|
from_stop = self.tf.create_segment("from_stop")
|
|
|
|
transport.extend_segment(to_stop, 0, stop, scale=self.speed)
|
2014-10-04 09:55:53 +08:00
|
|
|
# append the reverse transport (from stop to 0)
|
|
|
|
# both durations are the same in this case
|
2016-01-20 11:09:10 +08:00
|
|
|
transport.extend_segment(from_stop, 0, stop, scale=self.speed)
|
2014-10-19 13:09:10 +08:00
|
|
|
# distributes frames to the sub-devices in CompoundPDQ2
|
|
|
|
# and uploads them
|
2015-03-23 05:31:59 +08:00
|
|
|
self.electrodes.arm()
|
2014-10-04 09:55:53 +08:00
|
|
|
|
|
|
|
@kernel
|
|
|
|
def cool(self):
|
2016-01-20 11:09:10 +08:00
|
|
|
self.bd_sw.pulse(1*ms)
|
2014-10-04 09:55:53 +08:00
|
|
|
|
|
|
|
@kernel
|
|
|
|
def transport(self):
|
2015-03-23 05:31:59 +08:00
|
|
|
# selects transport frame
|
2014-10-04 09:55:53 +08:00
|
|
|
# triggers pdqs to start transport frame segment
|
|
|
|
# plays the transport waveform from 0 to stop
|
|
|
|
# delay()s the core by the duration of the waveform segment
|
2014-10-17 00:04:30 +08:00
|
|
|
self.tf.to_stop.advance()
|
2014-10-04 09:55:53 +08:00
|
|
|
# leaves the ion in the dark at the transport endpoint
|
|
|
|
delay(self.wait_at_stop)
|
|
|
|
# transport back (again: trigger, delay())
|
2014-10-17 00:04:30 +08:00
|
|
|
# segments can only be advance()ed in order
|
2015-03-23 05:31:59 +08:00
|
|
|
# since this is the last segment, pdq will go back to jump table
|
2014-10-17 00:04:30 +08:00
|
|
|
self.tf.from_stop.advance()
|
2014-10-04 09:55:53 +08:00
|
|
|
|
|
|
|
@kernel
|
|
|
|
def detect(self):
|
|
|
|
with parallel:
|
2016-01-20 11:09:10 +08:00
|
|
|
self.bd_sw.pulse(100*us)
|
2014-10-04 09:55:53 +08:00
|
|
|
self.pmt.gate_rising(100*us)
|
2016-01-20 11:09:10 +08:00
|
|
|
self.bd_sw.on()
|
2014-10-04 09:55:53 +08:00
|
|
|
return self.pmt.count()
|
|
|
|
|
|
|
|
@kernel
|
|
|
|
def one(self):
|
|
|
|
self.cool()
|
|
|
|
self.transport()
|
|
|
|
return self.detect()
|
|
|
|
|
|
|
|
@kernel
|
2015-01-10 15:41:35 +08:00
|
|
|
def repeat(self):
|
2016-01-20 11:09:10 +08:00
|
|
|
self.histogram[:] = [0 for _ in range(self.nbins)]
|
2014-10-04 09:55:53 +08:00
|
|
|
|
2015-01-10 15:41:35 +08:00
|
|
|
for i in range(self.repeats):
|
2014-10-04 09:55:53 +08:00
|
|
|
n = self.one()
|
2015-01-10 15:41:35 +08:00
|
|
|
if n >= self.nbins:
|
|
|
|
n = self.nbins - 1
|
2014-12-19 15:20:19 +08:00
|
|
|
self.histogram[n] += 1
|
2014-10-04 09:55:53 +08:00
|
|
|
|
2015-01-10 15:41:35 +08:00
|
|
|
def scan(self, stops):
|
2014-10-04 09:55:53 +08:00
|
|
|
for s in stops:
|
2014-10-17 00:04:30 +08:00
|
|
|
self.histogram = []
|
2015-07-14 04:08:20 +08:00
|
|
|
# non-kernel, build frames
|
2014-10-04 09:55:53 +08:00
|
|
|
# could also be rpc'ed from repeat()
|
2015-07-14 04:08:20 +08:00
|
|
|
self.calc_waveforms(s)
|
2014-10-04 09:55:53 +08:00
|
|
|
# kernel part
|
2015-01-10 15:41:35 +08:00
|
|
|
self.repeat()
|
2014-10-05 08:25:36 +08:00
|
|
|
# live update 2d plot with current self.histogram
|
|
|
|
# broadcast(s, self.histogram)
|
2014-10-04 09:55:53 +08:00
|
|
|
|
2015-01-10 15:41:35 +08:00
|
|
|
def run(self):
|
2014-10-04 09:55:53 +08:00
|
|
|
# scan transport endpoint
|
2016-01-20 11:09:10 +08:00
|
|
|
stops = range(10, len(transport.x), 10)
|
2015-01-10 15:41:35 +08:00
|
|
|
self.scan(stops)
|