generic line-/gridscan with profileMove

addams_bec/devices/detector.py

@@ -1,3 +1,5 @@
+import os
+
 from ophyd import (
     ADComponent as ADCpt,
     Device,
@@ -51,6 +53,14 @@ class Eiger500KSetup(CustomDetectorMixin):
             self.parent.cam.acquire.put(1, wait=False) # arm
 
         # file writer
+        # TODO: enable this when new storage path is consistent on bec server and consoles.
+        # self.parent.filepath.set(
+        #     self.parent.filewriter.compile_full_filename(f"{self.parent.name}")
+        # ).wait()
+        # file_path, file_name = os.path.split(self.parent.filepath.get())
+        # self.parent.hdf.file_path.put(file_path)
+        # self.parent.hdf.file_name.put(file_name)
+        # self.parent.hdf.file_template.put("%s%s")
         self.parent.hdf.lazy_open.put(1)
         self.parent.hdf.num_capture.put(num_points)
         self.parent.hdf.file_write_mode.put(2) # Stream
@@ -99,7 +109,8 @@ class Eiger500KSetup(CustomDetectorMixin):
 
         # publish timeseries data
         metadata = self.parent.scaninfo.scan_msg.metadata
-        metadata.update({"async_update": "append", "num_lines": self.parent.roistat.ts_current_point.get()})
+        metadata.update({"async_update": "append", "max_shape": [None, None]})
+
         msg = messages.DeviceMessage(
             signals={
                 self.parent.roistat.roi1.name_.get(): {

addams_bec/devices/profileMove.py

@@ -143,11 +143,11 @@ class ProfileMoveAxisXPS(ProfileMoveAxis):
 
 
 class ProfileMoveXPSX04SA(ProfileMoveControllerXPS):
-    alp = Cpt(ProfileMoveAxisXPS, '', motor='1')
+    ov = Cpt(ProfileMoveAxisXPS, '', motor='1')
-    delta = Cpt(ProfileMoveAxisXPS, '', motor='2')
+    xv = Cpt(ProfileMoveAxisXPS, '', motor='2')
-    gam = Cpt(ProfileMoveAxisXPS, '', motor='3')
+    oh = Cpt(ProfileMoveAxisXPS, '', motor='3')
-    ov = Cpt(ProfileMoveAxisXPS, '', motor='4')
+    phi = Cpt(ProfileMoveAxisXPS, '', motor='4')
-    xv = Cpt(ProfileMoveAxisXPS, '', motor='5')
+    nu = Cpt(ProfileMoveAxisXPS, '', motor='5')
-    phi = Cpt(ProfileMoveAxisXPS, '', motor='6')
+    alp = Cpt(ProfileMoveAxisXPS, '', motor='6')
-    oh = Cpt(ProfileMoveAxisXPS, '', motor='7')
+    gam = Cpt(ProfileMoveAxisXPS, '', motor='7')
-    nu = Cpt(ProfileMoveAxisXPS, '', motor='8')
+    delta = Cpt(ProfileMoveAxisXPS, '', motor='8')

addams_bec/scans/__init__.py

@@ -1,2 +1,4 @@
 from .hkl_scan import HklScan
-from .fly_scan import HklFlyScan
+from .fly_scan import HklFlyScan
+from .mesh_scan import HklMeshScan, HklMeshFlyScan
+from .profilemove_scan import LineProfileMove, GridProfileMove

addams_bec/scans/fly_scan.py

@@ -64,6 +64,7 @@ class HklFlyScan(AsyncFlyScanBase):
             self.positions.append(position[:-2])
 
     def scan_core(self):
+        all_motors = self.device_manager.devices[self.controller].axes.get()
         scan_motors = self.device_manager.devices[self.diffract].real_axes.get()
         hkls = numpy.linspace(self.start, self.stop, self.points).tolist()
         positions = yield from self.stubs.send_rpc_and_wait(self.diffract, 'angles_from_hkls', hkls)
@@ -80,11 +81,15 @@ class HklFlyScan(AsyncFlyScanBase):
             yield from self.stubs.send_rpc_and_wait(self.controller, 'time_mode.put', 1)
             yield from self.stubs.send_rpc_and_wait(self.controller, 'times.put', [total_time, total_time])
 
-            for index, axis_name in enumerate(scan_motors):
+            for axis_name in all_motors:
-                yield from self.stubs.send_rpc_and_wait(self.controller,
+                if axis_name not in scan_motors:
-                                                    f'{axis_name}.positions.put',
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 0)
-                                                    (positions[0, index], positions[-1, index]))
+                else:
-                yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
+                    index = scan_motors.index(axis_name)
+                    yield from self.stubs.send_rpc_and_wait(self.controller,
+                                                            f'{axis_name}.positions.put',
+                                                            (positions[0, index], positions[-1, index]))
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
 
         else:
             yield from self.stubs.send_rpc_and_wait(self.controller, 'num_points.put', num_points)
@@ -94,11 +99,15 @@ class HklFlyScan(AsyncFlyScanBase):
             yield from self.stubs.send_rpc_and_wait(self.controller, 'time_mode.put', 0)
             yield from self.stubs.send_rpc_and_wait(self.controller, 'fixed_time.put', self.exp_time)
 
-            for index, axis_name in enumerate(scan_motors):
+            for axis_name in all_motors:
-                yield from self.stubs.send_rpc_and_wait(self.controller,
+                if axis_name not in scan_motors:
-                                                    f'{axis_name}.positions.put',
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 0)
-                                                    positions[:, index])
+                else:
-                yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
+                    index = scan_motors.index(axis_name)
+                    yield from self.stubs.send_rpc_and_wait(self.controller,
+                                                            f'{axis_name}.positions.put',
+                                                            (positions[0, index], positions[-1, index]))
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
 
         yield from self.stubs.send_rpc_and_wait(self.controller, 'build_profile')
         build_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'build_status.get')
@@ -136,7 +145,7 @@ class HklFlyScan(AsyncFlyScanBase):
         logger.success(f'{self.scan_name} finished')
 
     def _publish_readbacks(self, device, readbacks):
-        metadata = {"async_update": "append", "num_lines": len(readbacks)}
+        metadata = {"async_update": "append", "max_shape": [None, None]}
         msg = messages.DeviceMessage(
             signals={device: {'value': readbacks} }, metadata=metadata
         )

addams_bec/scans/mesh_scan.py

@@ -0,0 +1,266 @@
+import itertools
+
+import numpy
+
+# from bec_lib.endpoints import MessageEndpoints
+from typing import Literal
+from bec_lib import messages
+from bec_lib.endpoints import MessageEndpoints
+from bec_lib.devicemanager import DeviceManagerBase
+from bec_lib.logger import bec_logger
+# from bec_lib import messages
+# from bec_server.scan_server.errors import ScanAbortion
+from bec_server.scan_server.scans import ScanAbortion, ScanArgType, ScanBase, AsyncFlyScanBase
+
+logger = bec_logger.logger
+
+class HklMeshScan(ScanBase):
+    scan_name = 'hklmesh_scan'
+    arg_input = {
+        'index': ScanArgType.DEVICE,
+        'start': ScanArgType.FLOAT,
+        'stop': ScanArgType.FLOAT,
+        'points': ScanArgType.INT
+    }
+    arg_bundle_size = {"bundle": len(arg_input), "min": 2, "max": 2}
+    required_kwargs = ['diffract', 'exp_time']
+
+    def __init__(self, *args,  diffract: str, **kwargs):
+        self.diffract = diffract
+        super().__init__(**kwargs)
+
+        if any(m not in ['h', 'k', 'l'] for m in self.caller_args):
+            raise ValueError("Invalid name. Must be 'h', 'k', or 'l'.")
+        if len(self.caller_args) != 2:
+            raise ValueError("Only 2 names can be given.")
+
+        self.scan_report_devices = ['h', 'k', 'l'] + self.scan_motors + self.readout_priority['monitored']
+
+    def update_scan_motors(self):
+        self.scan_motors = self.device_manager.devices[self.diffract].real_axes.get()
+
+    def prepare_positions(self):
+        """
+        Override base method to yield from _calculate_position method
+        """
+        yield from self._calculate_positions()
+        self._optimize_trajectory()
+        self.num_pos = len(self.positions) * self.burst_at_each_point
+        yield from self._set_position_offset()
+        self._check_limits()
+
+    def _calculate_positions(self):
+        inner_name, outer_name =  [x for x in self.caller_args if x in ['h','k','l']]
+        fixed_name = [x for x in ['h', 'k', 'l'] if x not in self.caller_args][0]
+        fixed_value = yield from self.stubs.send_rpc_and_wait(self.diffract, f'{fixed_name}.position')
+
+        print(f"{inner_name=}, {outer_name=}, {fixed_name=}")
+
+        inner_ind, outer_ind, fixed_ind = (['h', 'k', 'l'].index(x) for x in [inner_name, outer_name, fixed_name])
+
+        hkls = []
+        outer_start, outer_stop, outer_points = self.caller_args[outer_name]
+        inner_start, inner_stop, inner_points = self.caller_args[inner_name]
+
+        outer_vect = numpy.linspace(outer_start, outer_stop, outer_points, dtype=float)
+        for i, outer in enumerate(outer_vect):
+            if i % 2 == 0:
+                inner_vect = numpy.linspace(inner_start, inner_stop, inner_points, dtype=float)
+            else:
+                inner_vect = numpy.linspace(inner_stop, inner_start, inner_points, dtype=float)
+
+            for inner in inner_vect:
+                hkl = [0, 0, 0]
+                hkl[inner_ind] = inner
+                hkl[outer_ind] = outer
+                hkl[fixed_ind] = fixed_value
+                hkls.append(hkl)
+
+        positions = yield from self.stubs.send_rpc_and_wait(self.diffract, 'angles_from_hkls', hkls)
+        # the last two positions 'betaIn' and 'betaOut' are not real motors
+        self.positions = []
+        for position in positions:
+            self.positions.append(position[:-2])
+
+
+class HklMeshFlyScan(AsyncFlyScanBase):
+    scan_name = 'hklmesh_flyscan'
+    scan_type = 'fly'
+    arg_input = {
+        'index': ScanArgType.DEVICE,
+        'start': ScanArgType.FLOAT,
+        'stop': ScanArgType.FLOAT,
+        'points': ScanArgType.INT
+    }
+    arg_bundle_size = {"bundle": len(arg_input), "min": 2, "max": 2}
+    required_kwargs = ['diffract', 'controller']
+    use_scan_progress_report = False
+
+    def __init__(self, *args, diffract: str, controller: str, **kwargs):
+        self.diffract = diffract
+        self.controller = controller
+        super().__init__(**kwargs)
+
+        if any(m not in ['h', 'k', 'l'] for m in self.caller_args):
+            raise ValueError("Invalid name. Must be 'h', 'k', or 'l'.")
+        if len(self.caller_args) != 2:
+            raise ValueError("Only 2 names can be given.")
+
+        self.scan_report_devices = ['h', 'k', 'l'] + self.scan_motors + self.readout_priority['monitored']
+
+    @property
+    def monitor_sync(self):
+        return self.diffract
+
+    def update_scan_motors(self):
+        self.scan_motors = self.device_manager.devices[self.diffract].real_axes.get()
+
+    def update_readout_priority(self):
+        self.readout_priority['async'].extend(['h', 'k', 'l'])
+        self.readout_priority['async'].extend(self.scan_motors)
+
+    def prepare_positions(self):
+        """
+        Override base method to yield from _calculate_position method
+        """
+        yield from self._calculate_positions()
+        self._optimize_trajectory()
+        self.num_pos = len(self.positions) * self.burst_at_each_point
+        yield from self._set_position_offset()
+        self._check_limits()
+
+    def _calculate_positions(self):
+        inner_name, outer_name =  [x for x in self.caller_args if x in ['h','k','l']]
+        fixed_name = [x for x in ['h', 'k', 'l'] if x not in self.caller_args][0]
+        fixed_value = yield from self.stubs.send_rpc_and_wait(self.diffract, f'{fixed_name}.position')
+
+        inner_ind, outer_ind, fixed_ind = (['h', 'k', 'l'].index(x) for x in [inner_name, outer_name, fixed_name])
+
+        hkls = []
+        outer_start, outer_stop, outer_points = self.caller_args[outer_name]
+        inner_start, inner_stop, inner_points = self.caller_args[inner_name]
+
+        outer_vect = numpy.linspace(outer_start, outer_stop, outer_points, dtype=float)
+        for i, outer in enumerate(outer_vect):
+            if i % 2 == 0:
+                inner_vect = numpy.linspace(inner_start, inner_stop, inner_points, dtype=float)
+            else:
+                inner_vect = numpy.linspace(inner_stop, inner_start, inner_points, dtype=float)
+
+            for inner in inner_vect:
+                hkl = [0, 0, 0]
+                hkl[inner_ind] = inner
+                hkl[outer_ind] = outer
+                hkl[fixed_ind] = fixed_value
+                hkls.append(hkl)
+
+        positions = yield from self.stubs.send_rpc_and_wait(self.diffract, 'angles_from_hkls', hkls)
+        # the last two positions 'betaIn' and 'betaOut' are not real motors
+        self.positions = []
+        for position in positions:
+            self.positions.append(position[:-2])
+
+    def scan_core(self):
+        all_motors = self.device_manager.devices[self.controller].axes.get()
+
+        inner_name, outer_name =  [x for x in self.caller_args if x in ['h','k','l']]
+        fixed_name = [x for x in ['h', 'k', 'l'] if x not in self.caller_args][0]
+        fixed_value = yield from self.stubs.send_rpc_and_wait(self.diffract, f'{fixed_name}.position')
+
+        inner_ind, outer_ind, fixed_ind = (['h', 'k', 'l'].index(x) for x in [inner_name, outer_name, fixed_name])
+
+        num_pos = 0
+        hkls = []
+        outer_start, outer_stop, outer_points = self.caller_args[outer_name]
+        inner_start, inner_stop, inner_points = self.caller_args[inner_name]
+
+        outer_vect = numpy.linspace(outer_start, outer_stop, outer_points, dtype=float)
+        for i, outer in enumerate(outer_vect):
+            if i % 2 == 0:
+                start, stop = inner_start, inner_stop
+            else:
+                start, stop = inner_stop, inner_start
+
+            hkls = [[0, 0, 0], [0, 0, 0]]
+            hkls[0][inner_ind] = start
+            hkls[0][outer_ind] = outer
+            hkls[0][fixed_ind] = fixed_value
+
+            hkls[1][inner_ind] = stop
+            hkls[1][outer_ind] = outer
+            hkls[1][fixed_ind] = fixed_value
+
+            positions = yield from self.stubs.send_rpc_and_wait(self.diffract, 'angles_from_hkls', hkls)
+
+            # Move motors to start position, the last two positions 'betaIn' and 'betaOut' are not real motors
+            yield from self.stubs.set(device=self.scan_motors, value=positions[0][:-2])
+
+            inner_time = inner_points * self.exp_time
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'num_points.put', 2)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'num_pulses.put', inner_points)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'start_pulses.put', 1)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'end_pulses.put', 2)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'time_mode.put', 1)
+            print(f"{inner_time=} {inner_points=} {self.exp_time=}")
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'times.put', [inner_time, inner_time])
+
+            for axis_name in all_motors:
+                if axis_name not in self.scan_motors:
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 0)
+                else:
+                    index = self.scan_motors.index(axis_name)
+                    yield from self.stubs.send_rpc_and_wait(self.controller,
+                                                            f'{axis_name}.positions.put',
+                                                            (positions[0][index], positions[-1][index]))
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
+
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'build_profile')
+            build_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'build_status.get')
+            if build_status != 1:
+                raise ScanAbortion('Profile build failed')
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'execute_profile')
+            execute_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'execute_status.get')
+            if  execute_status != 1:
+                raise ScanAbortion('Profile execute failed')
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'readback_profile')
+            readback_status = yield from self.stubs.send_rpc_and_wait(self.controller,'readback_status.get')
+            if readback_status != 1:
+                raise ScanAbortion('Profile readback failed')
+
+            angle_readbacks = []
+            for index, axis_name in enumerate(self.scan_motors):
+                readbacks = yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.readbacks.get')
+                self._publish_readbacks(axis_name, readbacks)
+                angle_readbacks.append(readbacks)
+            # motor readbacks are aranged column-wise
+            angle_readbacks = [list(x) for x in zip(*angle_readbacks)]
+
+            hkls = yield from self.stubs.send_rpc_and_wait(self.diffract, 'hkls_from_angles', angle_readbacks)
+            hkls = numpy.array(hkls)
+
+            for index, name in enumerate(['h', 'k', 'l', 'betaIn', 'betaOut']):
+                self._publish_readbacks(name, hkls[:, index])
+
+            # motor readbacks have more points than generated pulses
+            num_pos += len(angle_readbacks)
+
+        self.num_pos = num_pos
+        logger.success(f'{self.scan_name} finished')
+
+    def _publish_readbacks(self, device, readbacks):
+        metadata = {"async_update": "append", "max_shape": [None, None]}
+        msg = messages.DeviceMessage(
+            signals={device: {'value': readbacks} }, metadata=metadata
+        )
+
+        self.stubs.connector.xadd(
+            topic=MessageEndpoints.device_async_readback(
+                scan_id=self.scan_id, device=device
+            ),
+            msg_dict={"data": msg},
+            expire=1800,
+        )

addams_bec/scans/profilemove_scan.py

@@ -0,0 +1,233 @@
+import numpy
+
+from bec_lib import messages
+from bec_lib.endpoints import MessageEndpoints
+from bec_lib.devicemanager import DeviceManagerBase
+from bec_lib.logger import bec_logger
+from bec_server.scan_server.scans import ScanAbortion, ScanArgType, ScanBase, AsyncFlyScanBase
+
+logger = bec_logger.logger
+
+class LineProfileMove(AsyncFlyScanBase):
+    scan_name = 'line_profilemove'
+    scan_type = 'fly'
+    arg_input = {
+        'motor': ScanArgType.DEVICE,
+        'start': ScanArgType.FLOAT,
+        'stop': ScanArgType.FLOAT,
+    }
+    arg_bundle_size = {"bundle": len(arg_input), "min": 1, "max": None}
+    required_kwargs = ['controller', 'points']
+    use_scan_progress_report = False
+
+    def __init__(self, *args,  controller: str, points: int, **kwargs):
+        self.controller = controller
+        self.points = self.num_pos = points
+        super().__init__(**kwargs)
+
+        self.all_motors = self.device_manager.devices[self.controller].axes.get()
+        if not all(m in self.all_motors for m in self.scan_motors):
+            raise ValueError(f"Scan motors {self.scan_motors} not in {self.all_motors}")
+
+    def update_readout_priority(self):
+        self.readout_priority['async'].extend(self.scan_motors)
+
+    def update_scan_motors(self):
+        ScanBase.update_scan_motors(self)
+
+    def prepare_positions(self):
+        """prepare the positions for the scan"""
+        self._calculate_positions()
+        self._optimize_trajectory()
+        self.num_pos = len(self.positions) * self.burst_at_each_point
+        yield from self._set_position_offset()
+        self._check_limits()
+
+    def _calculate_positions(self):
+        axis = []
+        for _, val in self.caller_args.items():
+            ax_pos = numpy.linspace(val[0], val[1], self.points, dtype=float)
+            axis.append(ax_pos)
+        self.positions = numpy.array(list(zip(*axis)), dtype=float)
+
+    def scan_core(self):
+
+        total_time = self.exp_time * self.points
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'num_points.put', 2)
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'num_pulses.put', self.points)
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'start_pulses.put', 1)
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'end_pulses.put', 2)
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'time_mode.put', 1)
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'times.put', [total_time, total_time])
+
+        for axis_name in self.all_motors:
+            if axis_name not in self.scan_motors:
+                yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 0)
+            else:
+                index = self.scan_motors.index(axis_name)
+                yield from self.stubs.send_rpc_and_wait(self.controller,
+                                                        f'{axis_name}.positions.put',
+                                                        (self.caller_args[axis_name][0], self.caller_args[axis_name][1]))
+                yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
+
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'build_profile')
+        build_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'build_status.get')
+        if build_status != 1:
+            raise ScanAbortion('Profile build failed')
+
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'execute_profile')
+        execute_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'execute_status.get')
+        if  execute_status != 1:
+            raise ScanAbortion('Profile execute failed')
+
+        yield from self.stubs.send_rpc_and_wait(self.controller, 'readback_profile')
+        readback_status = yield from self.stubs.send_rpc_and_wait(self.controller,'readback_status.get')
+        if readback_status != 1:
+            raise ScanAbortion('Profile readback failed')
+
+        for index, axis_name in enumerate(self.scan_motors):
+            readbacks = yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.readbacks.get')
+            self._publish_readbacks(axis_name, readbacks)
+
+        logger.success(f'{self.scan_name} finished')
+
+    def _publish_readbacks(self, device, readbacks):
+        metadata = {"async_update": "append", "max_shape": [None, None]}
+        msg = messages.DeviceMessage(
+            signals={device: {'value': readbacks} }, metadata=metadata
+        )
+
+        self.stubs.connector.xadd(
+            topic=MessageEndpoints.device_async_readback(
+                scan_id=self.scan_id, device=device
+            ),
+            msg_dict={"data": msg},
+            expire=1800,
+        )
+
+
+class GridProfileMove(AsyncFlyScanBase):
+    scan_name = 'grid_profilemove'
+    scan_type = 'fly'
+    arg_input = {
+        'motor': ScanArgType.DEVICE,
+        'start': ScanArgType.FLOAT,
+        'stop': ScanArgType.FLOAT,
+        'points': ScanArgType.INT,
+    }
+    arg_bundle_size = {"bundle": len(arg_input), "min": 2, "max": 2}
+    required_kwargs = ['controller']
+    use_scan_progress_report = False
+
+    def __init__(self, *args,  controller: str, **kwargs):
+        self.controller = controller
+        super().__init__(**kwargs)
+
+        self.all_motors = self.device_manager.devices[self.controller].axes.get()
+        if not all(m in self.all_motors for m in self.scan_motors):
+            raise ValueError(f"Scan motors {self.scan_motors} not in {self.all_motors}")
+
+    def update_readout_priority(self):
+        self.readout_priority['async'].extend(self.scan_motors)
+
+    def update_scan_motors(self):
+        ScanBase.update_scan_motors(self)
+
+    def prepare_positions(self):
+        """prepare the positions for the scan"""
+        self._calculate_positions()
+        self._optimize_trajectory()
+        self.num_pos = len(self.positions) * self.burst_at_each_point
+        yield from self._set_position_offset()
+        self._check_limits()
+
+    def _calculate_positions(self):
+        inner_motor, outer_motor =  self.scan_motors
+        outer_start, outer_stop, outer_points = self.caller_args[outer_motor]
+        inner_start, inner_stop, inner_points = self.caller_args[inner_motor]
+
+        self.positions = []
+        outer_vect = numpy.linspace(outer_start, outer_stop, outer_points, dtype=float)
+        for i, outer in enumerate(outer_vect):
+            if i % 2 == 0:
+                inner_vect = numpy.linspace(inner_start, inner_stop, inner_points, dtype=float)
+            else:
+                inner_vect = numpy.linspace(inner_stop, inner_start, inner_points, dtype=float)
+
+            for inner in inner_vect:
+                self.positions.append((inner, outer))
+
+    def scan_core(self):
+        inner_motor, outer_motor =  self.scan_motors
+        outer_start, outer_stop, outer_points = self.caller_args[outer_motor]
+        inner_start, inner_stop, inner_points = self.caller_args[inner_motor]
+        
+        num_pos = 0
+        outer_vect = numpy.linspace(outer_start, outer_stop, outer_points, dtype=float)
+        for i, outer in enumerate(outer_vect):
+            if i % 2 == 0:
+                start, stop = inner_start, inner_stop
+            else:
+                start, stop = inner_stop, inner_start
+
+            # Move outer motor to start position
+            yield from self.stubs.set(device=outer_motor, value=outer)
+            outer_readback = yield from self.stubs.send_rpc_and_wait(outer_motor, "position")
+
+            inner_time = self.exp_time * inner_points
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'num_points.put', 2)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'num_pulses.put', inner_points)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'start_pulses.put', 1)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'end_pulses.put', 2)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'time_mode.put', 1)
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'times.put', [inner_time, inner_time])
+
+            for axis_name in self.all_motors:
+                if axis_name == inner_motor:
+                    yield from self.stubs.send_rpc_and_wait(self.controller,
+                                                            f'{axis_name}.positions.put',
+                                                            (start, stop))
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 1)
+                else:
+                    yield from self.stubs.send_rpc_and_wait(self.controller, f'{axis_name}.use_axis.put', 0)
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'build_profile')
+            build_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'build_status.get')
+            if build_status != 1:
+                raise ScanAbortion('Profile build failed')
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'execute_profile')
+            execute_status = yield from self.stubs.send_rpc_and_wait(self.controller, 'execute_status.get')
+            if  execute_status != 1:
+                raise ScanAbortion('Profile execute failed')
+
+            yield from self.stubs.send_rpc_and_wait(self.controller, 'readback_profile')
+            readback_status = yield from self.stubs.send_rpc_and_wait(self.controller,'readback_status.get')
+            if readback_status != 1:
+                raise ScanAbortion('Profile readback failed')
+
+            readbacks = yield from self.stubs.send_rpc_and_wait(self.controller, f'{inner_motor}.readbacks.get')
+            self._publish_readbacks(inner_motor, readbacks)
+
+            readbacks = [outer_readback] * len(readbacks)
+            self._publish_readbacks(outer_motor, readbacks)
+
+            # motor readbacks have more points than generated pulses
+            num_pos += len(readbacks)
+
+        self.num_pos = num_pos
+        logger.success(f'{self.scan_name} finished')
+
+    def _publish_readbacks(self, device, readbacks):
+        metadata = {"async_update": "append", "max_shape": [None]}
+        msg = messages.DeviceMessage(
+            signals={device: {'value': readbacks} }, metadata=metadata
+        )
+
+        self.stubs.connector.xadd(
+            topic=MessageEndpoints.device_async_readback(
+                scan_id=self.scan_id, device=device
+            ),
+            msg_dict={"data": msg},
+            expire=1800,
+        )