dt5810 2.0.1__tar.gz

dt5810-2.0.1/LICENSE

@@ -0,0 +1,24 @@
+MIT License
+
+Copyright (c) 2026 Nuclear Instruments SRL
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.
+
+The CyAPI library (sdk/cpp/third_party/cyapi) is Copyright Cypress Semiconductor
+Corporation (Infineon Technologies) and is distributed under its own license terms.

dt5810-2.0.1/MANIFEST.in

@@ -0,0 +1,6 @@
+include README.md
+include build_wheel.py
+graft tests
+graft examples
+prune src/dt5810/_native
+global-exclude *.py[cod]

dt5810-2.0.1/PKG-INFO

@@ -0,0 +1,88 @@
+Metadata-Version: 2.4
+Name: dt5810
+Version: 2.0.1
+Summary: Python binding of the DT5810 SDK: control of the CAEN DT5810 / DT5810B Digital Detector Emulator
+Author-email: Nuclear Instruments <support@nuclearinstruments.eu>
+License: MIT
+Project-URL: Homepage, https://gitlab.nuclearinstruments.eu/public-repo/dt5810/dt5810-sdk
+Project-URL: Documentation, https://public-repo.pages.nuclearinstruments.eu/dt5810/dt5810-sdk/
+Project-URL: Source, https://gitlab.nuclearinstruments.eu/public-repo/dt5810/dt5810-sdk
+Project-URL: Issues, https://gitlab.nuclearinstruments.eu/public-repo/dt5810/dt5810-sdk/-/issues
+Keywords: DT5810,DT5810B,CAEN,Nuclear Instruments,detector emulator,pulser,nuclear physics
+Classifier: Development Status :: 5 - Production/Stable
+Classifier: Intended Audience :: Science/Research
+Classifier: Intended Audience :: Developers
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: Microsoft :: Windows
+Classifier: Operating System :: POSIX :: Linux
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3 :: Only
+Classifier: Programming Language :: Python :: 3.8
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
+Classifier: Topic :: Scientific/Engineering :: Physics
+Classifier: Topic :: System :: Hardware :: Hardware Drivers
+Classifier: Typing :: Typed
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Provides-Extra: numpy
+Requires-Dist: numpy; extra == "numpy"
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Dynamic: license-file
+
+# dt5810 — Python binding of the DT5810 SDK
+
+Control the **CAEN DT5810 / DT5810B Digital Detector Emulator** (Nuclear Instruments) from
+Python: energy (fixed or spectrum), pulse shape, time base, noise, baseline drift,
+correlation, digital I/O, run control, monitoring, AWG / sequences and the analog input
+MCA. The package calls the native DT5810 library (`dt5810.dll` / `libdt5810.so`) through
+`ctypes`; a built-in simulator (`sim://`) lets you develop and test without hardware.
+
+## Installation
+
+```
+pip install dt5810
+```
+
+The wheels for Windows x64 and Linux x86_64 (glibc 2.34+) include the native library.
+Python 3.8 or later; no mandatory dependency (`pip install dt5810[numpy]` to get numpy
+arrays from the read functions with `as_numpy=True`).
+
+With the source distribution, or to use another build of the library, point
+`DT5810_LIBRARY` at it (full path of `dt5810.dll` / `libdt5810.so`); otherwise the library
+is looked up in the package, then by the system loader.
+
+## Quick start
+
+```python
+import dt5810
+
+print(dt5810.version())
+for d in dt5810.enumerate(dt5810.Bus.USB):
+    print(d.serial, d.address)
+
+with dt5810.open("sim://") as dev:          # or "usb://", "usb://NI050123"
+    print(dev.info.model_name, dev.info.serial)
+    ch = dev.channel(0)
+    ch.set_fixed_energy(20000)              # LSB, 0..32767
+    ch.set_drc(rise_ns=20, decay_ns=50000)
+    ch.set_timebase(dt5810.TimebaseMode.POISSON, rate_hz=10000)
+    ch.enable()
+    print(ch.stats())
+```
+
+Errors raise `dt5810.Dt5810Error` subclasses (`OutOfRangeError`, `InvalidParameterError`,
+`NoDataError`, `NotFoundError`, `DeviceTimeoutError`, ...) carrying `.status` and
+`.message`. The `examples/` folder has runnable scripts for every feature
+(`python examples/01_info.py usb://`).
+
+## Documentation
+
+User guide and API reference: <https://public-repo.pages.nuclearinstruments.eu/dt5810/dt5810-sdk/>
+
+Source: <https://gitlab.nuclearinstruments.eu/public-repo/dt5810/dt5810-sdk> — MIT License.

dt5810-2.0.1/README.md

@@ -0,0 +1,51 @@
+# dt5810 — Python binding of the DT5810 SDK
+
+Control the **CAEN DT5810 / DT5810B Digital Detector Emulator** (Nuclear Instruments) from
+Python: energy (fixed or spectrum), pulse shape, time base, noise, baseline drift,
+correlation, digital I/O, run control, monitoring, AWG / sequences and the analog input
+MCA. The package calls the native DT5810 library (`dt5810.dll` / `libdt5810.so`) through
+`ctypes`; a built-in simulator (`sim://`) lets you develop and test without hardware.
+
+## Installation
+
+```
+pip install dt5810
+```
+
+The wheels for Windows x64 and Linux x86_64 (glibc 2.34+) include the native library.
+Python 3.8 or later; no mandatory dependency (`pip install dt5810[numpy]` to get numpy
+arrays from the read functions with `as_numpy=True`).
+
+With the source distribution, or to use another build of the library, point
+`DT5810_LIBRARY` at it (full path of `dt5810.dll` / `libdt5810.so`); otherwise the library
+is looked up in the package, then by the system loader.
+
+## Quick start
+
+```python
+import dt5810
+
+print(dt5810.version())
+for d in dt5810.enumerate(dt5810.Bus.USB):
+    print(d.serial, d.address)
+
+with dt5810.open("sim://") as dev:          # or "usb://", "usb://NI050123"
+    print(dev.info.model_name, dev.info.serial)
+    ch = dev.channel(0)
+    ch.set_fixed_energy(20000)              # LSB, 0..32767
+    ch.set_drc(rise_ns=20, decay_ns=50000)
+    ch.set_timebase(dt5810.TimebaseMode.POISSON, rate_hz=10000)
+    ch.enable()
+    print(ch.stats())
+```
+
+Errors raise `dt5810.Dt5810Error` subclasses (`OutOfRangeError`, `InvalidParameterError`,
+`NoDataError`, `NotFoundError`, `DeviceTimeoutError`, ...) carrying `.status` and
+`.message`. The `examples/` folder has runnable scripts for every feature
+(`python examples/01_info.py usb://`).
+
+## Documentation
+
+User guide and API reference: <https://public-repo.pages.nuclearinstruments.eu/dt5810/dt5810-sdk/>
+
+Source: <https://gitlab.nuclearinstruments.eu/public-repo/dt5810/dt5810-sdk> — MIT License.

dt5810-2.0.1/build_wheel.py

@@ -0,0 +1,195 @@
+#!/usr/bin/env python3
+"""Builds a platform wheel of the dt5810 package with the native library inside.
+
+    python build_wheel.py --lib <dt5810.dll | libdt5810.so> --plat-name <tag> [--outdir dist]
+
+Examples of `--plat-name`: ``win_amd64``, ``manylinux_2_34_x86_64``.
+
+Steps: copy the library into ``src/dt5810/_native/``, build the (pure) wheel with
+``python -m build --wheel`` (or ``pip wheel --no-deps`` when ``build`` is not installed),
+retag it ``py3-none-<plat>`` (WHEEL file, RECORD and file name rewritten here: no extra
+tool needed), check that it carries the library, and empty ``_native/`` again.
+
+The build does not use an isolated environment when setuptools is already installed
+(no network needed); pass ``--isolation`` to force one.
+"""
+
+from __future__ import annotations
+
+import argparse
+import base64
+import glob
+import hashlib
+import importlib.util
+import os
+import shutil
+import subprocess
+import sys
+import tempfile
+import zipfile
+
+HERE = os.path.dirname(os.path.abspath(__file__))
+NATIVE = os.path.join(HERE, "src", "dt5810", "_native")
+NATIVE_NAMES = ("dt5810.dll", "libdt5810.so", "libdt5810.dylib")
+
+
+def log(msg: str) -> None:
+    print("[build_wheel] " + msg, flush=True)
+
+
+def native_name(lib: str, plat: str) -> str:
+    if plat.startswith("win"):
+        return "dt5810.dll"
+    if plat.startswith("macosx"):
+        return "libdt5810.dylib"
+    if lib.lower().endswith(".dll"):
+        raise SystemExit("a .dll cannot go in a %s wheel" % plat)
+    return "libdt5810.so"
+
+
+def clean_native() -> None:
+    for name in NATIVE_NAMES:
+        p = os.path.join(NATIVE, name)
+        if os.path.exists(p):
+            os.remove(p)
+
+
+def have(module: str) -> bool:
+    return importlib.util.find_spec(module) is not None
+
+
+def setuptools_can_build_wheels() -> bool:
+    if not have("setuptools"):
+        return False
+    if have("wheel"):
+        return True
+    try:  # setuptools >= 70.1 has its own bdist_wheel
+        import setuptools
+
+        major, minor = (int(x) for x in setuptools.__version__.split(".")[:2])
+        return (major, minor) >= (70, 1)
+    except Exception:
+        return False
+
+
+def build_pure_wheel(tmpdir: str, isolation: bool) -> str:
+    for d in ("build", os.path.join("src", "dt5810.egg-info")):
+        shutil.rmtree(os.path.join(HERE, d), ignore_errors=True)
+    no_isolation = not isolation and setuptools_can_build_wheels()
+    if have("build"):
+        cmd = [sys.executable, "-m", "build", "--wheel", "--outdir", tmpdir]
+        if no_isolation:
+            cmd.append("--no-isolation")
+    else:
+        cmd = [sys.executable, "-m", "pip", "wheel", "--no-deps", "-w", tmpdir]
+        if no_isolation:
+            cmd.append("--no-build-isolation")
+        cmd.append(HERE)
+    log(" ".join(cmd))
+    subprocess.check_call(cmd, cwd=HERE)
+    wheels = glob.glob(os.path.join(tmpdir, "dt5810-*.whl"))
+    if len(wheels) != 1:
+        raise SystemExit("expected one wheel in %s, found %r" % (tmpdir, wheels))
+    for d in ("build", os.path.join("src", "dt5810.egg-info")):
+        shutil.rmtree(os.path.join(HERE, d), ignore_errors=True)
+    return wheels[0]
+
+
+def record_hash(data: bytes) -> str:
+    digest = hashlib.sha256(data).digest()
+    return "sha256=" + base64.urlsafe_b64encode(digest).rstrip(b"=").decode("ascii")
+
+
+def retag(wheel: str, plat: str, outdir: str) -> str:
+    """Rewrites `wheel` as py3-none-<plat> into `outdir`; returns the new path."""
+    base = os.path.basename(wheel)
+    parts = base[:-4].split("-")
+    if len(parts) < 5:
+        raise SystemExit("unexpected wheel name " + base)
+    name, ver = parts[0], parts[1]
+    tag = "py3-none-" + plat
+    target = os.path.join(outdir, "%s-%s-%s.whl" % (name, ver, tag))
+    dist_info = "%s-%s.dist-info" % (name, ver)
+    with zipfile.ZipFile(wheel) as zin:
+        entries = [(i, zin.read(i.filename)) for i in zin.infolist()]
+    files = {}
+    order = []
+    for info, data in entries:
+        order.append(info)
+        files[info.filename] = data
+    wheel_meta = files[dist_info + "/WHEEL"].decode("utf-8").splitlines()
+    out_lines = []
+    for line in wheel_meta:
+        if line.startswith("Tag:"):
+            continue
+        if line.startswith("Root-Is-Purelib:"):
+            line = "Root-Is-Purelib: false"
+        out_lines.append(line)
+    while out_lines and not out_lines[-1].strip():
+        out_lines.pop()
+    out_lines.append("Tag: " + tag)
+    files[dist_info + "/WHEEL"] = ("\n".join(out_lines) + "\n").encode("utf-8")
+    record_name = dist_info + "/RECORD"
+    record = []
+    for info in order:
+        if info.filename == record_name or info.filename.endswith("/"):
+            continue
+        data = files[info.filename]
+        record.append("%s,%s,%d" % (info.filename, record_hash(data), len(data)))
+    record.append(record_name + ",,")
+    files[record_name] = ("\n".join(record) + "\n").encode("utf-8")
+    if os.path.exists(target):
+        os.remove(target)
+    with zipfile.ZipFile(target, "w", zipfile.ZIP_DEFLATED) as zout:
+        for info in order:
+            zi = zipfile.ZipInfo(info.filename, date_time=info.date_time)
+            zi.external_attr = info.external_attr
+            zi.compress_type = zipfile.ZIP_DEFLATED
+            zout.writestr(zi, files[info.filename])
+    return target
+
+
+def main() -> int:
+    ap = argparse.ArgumentParser(description=__doc__.split("\n\n")[0])
+    ap.add_argument("--lib", required=True, help="path of dt5810.dll or libdt5810.so")
+    ap.add_argument("--plat-name", required=True,
+                    help="platform tag, e.g. win_amd64 or manylinux_2_34_x86_64")
+    ap.add_argument("--outdir", default=os.path.join(HERE, "dist"))
+    ap.add_argument("--isolation", action="store_true",
+                    help="build in an isolated environment (downloads setuptools)")
+    ap.add_argument("--keep-native", action="store_true",
+                    help="leave the library in src/dt5810/_native/ afterwards")
+    args = ap.parse_args()
+
+    lib = os.path.abspath(args.lib)
+    if not os.path.isfile(lib):
+        raise SystemExit("library not found: " + lib)
+    plat = args.plat_name.replace("-", "_").replace(".", "_")
+    outdir = os.path.abspath(args.outdir)
+    os.makedirs(outdir, exist_ok=True)
+
+    os.makedirs(NATIVE, exist_ok=True)
+    clean_native()
+    dest = os.path.join(NATIVE, native_name(lib, plat))
+    shutil.copy2(lib, dest)
+    log("copied %s -> %s" % (lib, dest))
+    try:
+        with tempfile.TemporaryDirectory() as tmp:
+            pure = build_pure_wheel(tmp, args.isolation)
+            target = retag(pure, plat, outdir)
+    finally:
+        if not args.keep_native:
+            clean_native()
+
+    with zipfile.ZipFile(target) as z:
+        names = z.namelist()
+    inside = "dt5810/_native/" + os.path.basename(dest)
+    if inside not in names:
+        raise SystemExit("the wheel does not contain " + inside)
+    log("wheel: " + target)
+    print(target)
+    return 0
+
+
+if __name__ == "__main__":
+    sys.exit(main())

dt5810-2.0.1/examples/01_info.py

@@ -0,0 +1,38 @@
+"""Lists the instruments and prints the identity of one of them.
+
+    python 01_info.py [address]        (default "sim://", e.g. "usb://" or "usb://NI050123")
+"""
+
+import sys
+
+import dt5810
+
+
+def main(address: str) -> None:
+    print("dt5810 package %s, native library %s" % (dt5810.__version__, dt5810.version()))
+    print("buses compiled in: %s" % ", ".join(b.name for b in dt5810.supported_buses()))
+
+    buses = dt5810.Bus.SIMULATOR if address.startswith("sim://") else dt5810.Bus.USB
+    for d in dt5810.enumerate(buses):
+        print("found %-10s %-20s %s" % (d.serial or "?", d.address, d.description))
+
+    with dt5810.open(address) as dev:
+        info = dev.info
+        print("model          %s (code 0x%08X)" % (info.model_name, info.model_code))
+        print("serial         %s" % info.serial)
+        print("firmware       0x%08X" % info.firmware_release)
+        print("outputs/inputs %d / %d" % (info.output_channels, info.input_channels))
+        print("clock          %.0f MHz, DAC %.0f MS/s"
+              % (info.clock_hz / 1e6, info.sample_rate_hz / 1e6))
+        print("running        %s" % dt5810.RunningMode(dev.running_mode()).name)
+        act = dev.get_activation()
+        print("activation     active=%s trial_counter=%d expired=%s" % tuple(act))
+        clk = dev.get_clock_config()
+        print("clock config   external_in=%s clock_out=%s" % tuple(clk))
+        for i in range(4):
+            cal = dev.read_calibration(i)
+            print("calibration %d  offset=%g gain=%g lsb_to_volt=%g" % ((i,) + tuple(cal)))
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://")

dt5810-2.0.1/examples/02_fixed_energy_pulser.py

@@ -0,0 +1,31 @@
+"""A periodic pulser: fixed amplitude, exponential shape, constant rate.
+
+    python 02_fixed_energy_pulser.py [address]      (default "sim://")
+"""
+
+import sys
+import time
+
+import dt5810
+
+
+def main(address: str) -> None:
+    with dt5810.open(address) as dev:
+        ch = dev.channel(0)
+        ch.set_output(gain=1.0, offset=0, invert=False, filter=True,
+                      range=dt5810.OutputRange.HIGH_PERFORMANCE)
+        ch.set_fixed_energy(16000)                         # LSB (0..32767)
+        ch.set_drc(rise_ns=20, decay_ns=50000)             # exponential pulse
+        ch.set_timebase(dt5810.TimebaseMode.CONSTANT, rate_hz=1000)
+        ch.set_noise()                                     # no noise
+        ch.enable()
+        print("amplitude: 16000 LSB = %.3f V" % ch.lsb_to_volt(16000))
+
+        time.sleep(1.0)
+        st = ch.stats()
+        print("rate %d cps, %d events" % (st.count_rate, st.event_counter))
+        ch.disable()
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://")

dt5810-2.0.1/examples/03_spectrum_emulation.py

@@ -0,0 +1,51 @@
+"""Emulates a gamma spectrum: three Gaussian peaks on an exponential background,
+Poisson arrivals and some Gaussian noise.
+
+    python 03_spectrum_emulation.py [address]       (default "sim://")
+"""
+
+import math
+import sys
+import time
+
+import dt5810
+
+PEAKS = [  # (bin, sigma in bins, relative intensity)
+    (2000, 25, 1.0),
+    (6000, 40, 0.6),
+    (11000, 60, 0.3),
+]
+
+
+def make_spectrum():
+    s = []
+    for i in range(dt5810.SPECTRUM_BINS):
+        v = 0.02 * math.exp(-i / 3000.0)                    # background
+        for centre, sigma, height in PEAKS:
+            v += height * math.exp(-0.5 * ((i - centre) / sigma) ** 2)
+        s.append(v)
+    return s
+
+
+def main(address: str) -> None:
+    with dt5810.open(address) as dev:
+        ch = dev.channel(0)
+        ch.set_spectrum(make_spectrum())                    # program + select the spectrum
+        ch.set_drc(rise_ns=30, decay_ns=40000)
+        ch.set_timebase(dt5810.TimebaseMode.POISSON, rate_hz=20000, dead_time_ns=200)
+        ch.set_noise(gauss=100)
+        ch.reset_spectrum_loopback()
+        ch.enable()
+
+        time.sleep(2.0)
+        spectrum = ch.spectrum_loopback()                   # 4096 bins
+        total = sum(spectrum)
+        print("emitted spectrum: %d counts in %d bins" % (total, len(spectrum)))
+        if total:
+            top = max(range(len(spectrum)), key=spectrum.__getitem__)
+            print("highest bin: %d (%d counts)" % (top, spectrum[top]))
+        ch.disable()
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://")

dt5810-2.0.1/examples/04_custom_shape.py

@@ -0,0 +1,36 @@
+"""Custom pulse shapes: one shape from samples, then two shapes with probabilities.
+
+    python 04_custom_shape.py [address]      (default "sim://")
+"""
+
+import math
+import sys
+
+import dt5810
+
+
+def main(address: str) -> None:
+    with dt5810.open(address) as dev:
+        ch = dev.channel(0)
+        ts = 1.0 / dev.info.sample_rate_hz                  # sample period in s
+
+        # the "fast" exponential of the Control Software: 50 ns rise, 5 us decay
+        fast = dt5810.make_fast_exponential(rise_us=0.05, decay_us=5.0, sample_period_s=ts,
+                                            length=dt5810.CUSTOM_SHAPE_MAX_POINTS)
+        ch.set_custom_shape(fast, pileup_limit=16)
+        print("custom shape: %d samples, peak %.0f" % (len(fast), max(fast)))
+
+        # a slower, Gaussian-like shape mixed with the first one: 70 % / 30 %
+        slow = [32575 * math.exp(-0.5 * ((i - 300) / 80.0) ** 2) for i in range(1200)]
+        ch.set_multi_shapes([fast, slow], [70, 30], pileup_limit=8)
+        print("multishape: 70 % fast exponential, 30 % gaussian")
+
+        ch.set_fixed_energy(20000)
+        ch.set_timebase(dt5810.TimebaseMode.POISSON, rate_hz=5000)
+        ch.enable()
+        ch.disable()
+        ch.disable_custom_shape()                           # back to the digital RC shaper
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://")

dt5810-2.0.1/examples/05_monitor.py

@@ -0,0 +1,42 @@
+"""Monitoring: statistics, run-control live data and the signal / spectrum loopback.
+
+    python 05_monitor.py [address] [seconds]      (default "sim://", 3 s)
+"""
+
+import sys
+import time
+
+import dt5810
+
+
+def main(address: str, seconds: float) -> None:
+    with dt5810.open(address) as dev:
+        ch = dev.channel(0)
+        ch.set_fixed_energy(12000)
+        ch.set_drc(rise_ns=20, decay_ns=20000)
+        ch.set_timebase(dt5810.TimebaseMode.POISSON, rate_hz=50000)
+        ch.set_run_control(dt5810.RunMode.FREE_RUNNING)
+        ch.run_control_reset_stats()
+        ch.reset_counter()
+        ch.enable()
+
+        end = time.time() + seconds
+        while time.time() < end:
+            time.sleep(0.5)
+            st = ch.stats()
+            ld = ch.live_data()
+            print("rate %6d cps  events %9d  real %.2f s  run %.2f s  busy=%s sat=%s"
+                  % (st.count_rate, st.event_counter, ld.real_time_s, ld.run_time_s,
+                     ld.busy, ld.saturated))
+
+        signal = ch.signal_loopback()          # 8192 samples
+        spectrum = ch.spectrum_loopback()      # 4096 bins
+        print("signal loopback: %d samples, min %d, max %d"
+              % (len(signal), min(signal), max(signal)))
+        print("spectrum loopback: %d bins, %d counts" % (len(spectrum), sum(spectrum)))
+        ch.disable()
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://",
+         float(sys.argv[2]) if len(sys.argv) > 2 else 3.0)

dt5810-2.0.1/examples/06_mca_input.py

@@ -0,0 +1,51 @@
+"""The analog input MCA: trapezoidal filter, spectrum acquisition and preview.
+
+    python 06_mca_input.py [address] [seconds]      (default "sim://", 2 s)
+"""
+
+import sys
+import time
+
+import dt5810
+
+
+def main(address: str, seconds: float) -> None:
+    with dt5810.open(address) as dev:
+        dev.mca_set_input(impedance=0, divider=0, mode=0, scale=0, gain=0, offset=0)
+
+        cfg = dt5810.mca_default_config()
+        cfg.peaking_us = 2.0
+        cfg.flat_top_us = 0.5
+        cfg.decay_tau_us = 50.0
+        dev.mca_configure(cfg)
+
+        dev.mca_reset()
+        dev.mca_clear_spectrum()
+        dev.mca_run(True)
+        time.sleep(seconds)
+        spectrum = dev.mca_get_spectrum()                   # 16384 bins
+        dev.mca_stop()
+        total = sum(spectrum)
+        print("MCA spectrum: %d counts in %d bins" % (total, len(spectrum)))
+        if total:
+            peak = max(range(len(spectrum)), key=spectrum.__getitem__)
+            print("peak at bin %d" % peak)
+
+        dev.mca_preview_configure(analog1=0, analog2=1, digital1=0, digital2=1,
+                                  trigger_source=0, trigger_level=500, post_trigger=1024)
+        for _ in range(20):
+            try:
+                p = dev.mca_preview_read()
+                break
+            except dt5810.NoDataError:                      # not triggered yet
+                time.sleep(0.1)
+        else:
+            print("the preview did not trigger")
+            return
+        print("preview: %d samples, analog1 max %d, analog2 max %d"
+              % (len(p.analog1), max(p.analog1), max(p.analog2)))
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://",
+         float(sys.argv[2]) if len(sys.argv) > 2 else 2.0)

dt5810-2.0.1/examples/07_awg.py

@@ -0,0 +1,37 @@
+"""Arbitrary waveform generator: a sine wave on channel 0, a pulse sequence on channel 1.
+
+    python 07_awg.py [address]      (default "sim://")
+"""
+
+import math
+import sys
+
+import dt5810
+
+
+def main(address: str) -> None:
+    with dt5810.open(address) as dev:
+        awg = dev.channel(0)
+        samples = [int(15000 * math.sin(2 * math.pi * i / 256)) for i in range(4096)]
+        awg.set_awg_waveform(samples)                 # padded to a multiple of 16
+        awg.set_awg_scaler(0)                         # one sample per clock period
+        awg.set_mode(dt5810.ChannelMode.AWG)
+        awg.enable()
+        print("channel 0: AWG, %d samples" % len(samples))
+
+        seq = dev.channel(1)
+        energies = [4000, 8000, 12000, 16000]
+        times = [2500, 2500, 2500, 25000]             # clock periods after the previous pulse
+        seq.set_sequence(energy=energies, time=times)
+        seq.set_drc(rise_ns=20, decay_ns=10000)
+        seq.enable()
+        print("channel 1: sequence of %d pulses" % len(energies))
+
+        seq.disable()
+        seq.clear_sequence()
+        awg.disable()
+        awg.set_mode(dt5810.ChannelMode.PULSER)
+
+
+if __name__ == "__main__":
+    main(sys.argv[1] if len(sys.argv) > 1 else "sim://")