artiq/examples/master/repository/photon_histogram.py

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()