forked from M-Labs/artiq
1
0
Fork 0
artiq/examples/master/repository/photon_histogram.py

58 lines
1.3 KiB
Python
Raw Normal View History

2014-09-30 17:38:02 +08:00
from artiq import *
2014-09-13 19:38:15 +08:00
class PhotonHistogram(Experiment, AutoDB):
"""Photon histogram"""
class DBKeys:
2015-03-08 18:37:53 +08:00
core = Device()
bd = Device()
bdd = Device()
pmt = Device()
2014-09-13 19:38:15 +08:00
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()
2014-09-13 19:38:15 +08:00
@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)
2014-09-13 19:38:15 +08:00
with parallel:
self.bd.pulse(self.detect_f, self.detect_t)
self.pmt.gate_rising(self.detect_t)
2014-09-13 19:38:15 +08:00
self.bd.on(200*MHz)
self.bdd.on(300*MHz)
return self.pmt.count()
2014-09-13 19:38:15 +08:00
@kernel
def run(self):
hist = [0 for _ in range(self.nbins)]
total = 0
2014-09-13 19:38:15 +08:00
for i in range(self.repeats):
2014-09-13 19:38:15 +08:00
n = self.cool_detect()
if n >= self.nbins:
n = self.nbins - 1
2014-09-13 19:38:15 +08:00
hist[n] += 1
total += n
self.hist = hist
self.total = total
self.ion_present = total > 5*self.repeats
2014-09-13 19:38:15 +08:00
if __name__ == "__main__":
from artiq.frontend.artiq_run import run
run()