Stimulator: scan ALL underlying epochprobemaps for channels

MATLAB's getchanneldevinfo iterates ALL epochprobemaps in the
underlying epoch and extracts channels from every matching one.
The Python version only looked at the single matching epm stored
in the probe's epoch table entry.

Also print all underlying epochprobemaps and their devicestrings
in the test diagnostic to understand what channels are available.

https://claude.ai/code/session_01A7rAxYf5pSvs19iVJe3ncL

Author: claude

src/ndi/probe/timeseries_stimulator.py

@@ -431,13 +431,27 @@ def getchanneldevinfo(
         if dev is None:
             return None
 
-        # Get channel info from the epochprobemap's devicestring
-        epms = base_info.get("epochprobemap", [])
+        # Get channel info from ALL epochprobemaps in the underlying epoch
+        # that match this probe. MATLAB iterates all maps in the underlying
+        # epoch and extracts channels from every matching one.
+        et, _ = self.epochtable()
+        entry = et[epoch - 1] if isinstance(epoch, int) and epoch <= len(et) else None
+        underlying = entry.get("underlying_epochs", {}) if entry else {}
+        all_epms = underlying.get("epochprobemap", base_info.get("epochprobemap", []))
+        if not isinstance(all_epms, list):
+            all_epms = [all_epms]
+
         channeltype = []
         channel = []
 
-        for epm in epms:
+        for epm in all_epms:
+            if not self.epochprobemapmatch(epm):
+                continue
             if hasattr(epm, "devicestring") and epm.devicestring:
+                logger.debug(
+                    "stimulator: matched epm devicestring='%s'",
+                    epm.devicestring,
+                )
                 try:
                     from ..daq.daqsystemstring import ndi_daq_daqsystemstring
 

tests/test_cloud_read_ingested.py

@@ -283,13 +283,20 @@ def test_stimulator_probe_timeseries(self, session):
         ch = devinfo.get("channel", [])
         print(f"  dev={type(dev).__name__}, devepoch={devepoch}")
         print(f"  channeltype={ct}, channel={ch}")
-        # Print the epochprobemap devicestring
-        et, _ = stim.epochtable()
-        if et:
-            entry_epm = et[0].get("epochprobemap", [])
-            for m in (entry_epm if isinstance(entry_epm, list) else [entry_epm]):
-                if hasattr(m, "devicestring"):
-                    print(f"  devicestring: {m.devicestring}")
+        # Print ALL epochprobemaps from underlying epoch
+        et_stim, _ = stim.epochtable()
+        if et_stim:
+            underlying = et_stim[0].get("underlying_epochs", {})
+            all_epms = underlying.get("epochprobemap", [])
+            if not isinstance(all_epms, list):
+                all_epms = [all_epms]
+            print(f"  underlying epochprobemaps count: {len(all_epms)}")
+            for i, m in enumerate(all_epms):
+                ds = getattr(m, "devicestring", "?")
+                nm = getattr(m, "name", "?")
+                print(f"  epm[{i}]: name={nm} devicestring={ds}")
+                match = stim.epochprobemapmatch(m) if hasattr(stim, "epochprobemapmatch") else "?"
+                print(f"    matches this probe: {match}")
 
         # Try readevents directly to see the error
         if dev is not None and ct: