forked from M-Labs/artiq
58 lines
1.3 KiB
Python
58 lines
1.3 KiB
Python
from artiq import *
|
|
|
|
|
|
class PhotonHistogram(Experiment, AutoDB):
|
|
"""Photon histogram"""
|
|
|
|
class DBKeys:
|
|
core = Device()
|
|
bd = Device()
|
|
bdd = Device()
|
|
pmt = Device()
|
|
|
|
nbins = Argument(100)
|
|
repeats = Argument(100)
|
|
|
|
cool_f = Parameter(230*MHz)
|
|
detect_f = Parameter(220*MHz)
|
|
detect_t = Parameter(100*us)
|
|
|
|
ion_present = Parameter(True)
|
|
|
|
hist = Result()
|
|
total = Result()
|
|
|
|
@kernel
|
|
def cool_detect(self):
|
|
with parallel:
|
|
self.bd.pulse(200*MHz, 1*ms)
|
|
self.bdd.pulse(300*MHz, 1*ms)
|
|
self.bd.pulse(self.cool_f, 100*us)
|
|
with parallel:
|
|
self.bd.pulse(self.detect_f, self.detect_t)
|
|
self.pmt.gate_rising(self.detect_t)
|
|
self.bd.on(200*MHz)
|
|
self.bdd.on(300*MHz)
|
|
return self.pmt.count()
|
|
|
|
@kernel
|
|
def run(self):
|
|
hist = [0 for _ in range(self.nbins)]
|
|
total = 0
|
|
|
|
for i in range(self.repeats):
|
|
n = self.cool_detect()
|
|
if n >= self.nbins:
|
|
n = self.nbins - 1
|
|
hist[n] += 1
|
|
total += n
|
|
|
|
self.hist = hist
|
|
self.total = total
|
|
self.ion_present = total > 5*self.repeats
|
|
|
|
|
|
if __name__ == "__main__":
|
|
from artiq.frontend.artiq_run import run
|
|
run()
|