diff --git a/src/dartsort/util/analysis.py b/src/dartsort/util/analysis.py
index 8ac42841..599a6e2a 100644
--- a/src/dartsort/util/analysis.py
+++ b/src/dartsort/util/analysis.py
@@ -318,8 +318,6 @@ def unit_raw_waveforms(
             spike_length_samples=spike_length_samples,
             fill_value=np.nan,
         )
-        print(f"analysis unit_raw_waveforms {waveforms.shape=} {load_ci.shape=}")
-        print(f"analysis unit_raw_waveforms {self.shifting=}")
         if not self.shifting:
             return which, waveforms
 
@@ -434,8 +432,6 @@ def unit_shift_or_relocate_channels(
         show_chans = show_chans[show_chans < len(show_geom)]
         show_channel_index = np.broadcast_to(show_chans[None], (len(show_geom), show_chans.size))
 
-        print(f"analysis relocate {load_channel_index.shape=} {show_chans.shape=}")
-
         if not self.shifting:
             return waveforms, max_chan, show_geom, show_channel_index
 
diff --git a/src/dartsort/vis/unit.py b/src/dartsort/vis/unit.py
index 9f39ee79..26338c22 100644
--- a/src/dartsort/vis/unit.py
+++ b/src/dartsort/vis/unit.py
@@ -341,8 +341,6 @@ def draw(self, axis, sorting_analysis, unit_id):
             show_superres_templates = suptemplates.shape[0] > 1
             max_abs_amp = self.max_abs_template_scale * np.abs(suptemplates).max()
 
-        print(f"{waveforms.shape=} {ci.shape=}")
-
         ls = geomplot(
             waveforms,
             max_channels=np.full(len(waveforms), max_chan),
@@ -574,7 +572,6 @@ def make_unit_summary(
             axes = cardfig.subplots(nrows=len(card.plots), ncols=1)
             axes = np.atleast_1d(axes)
             for plot, axis in zip(card.plots, axes):
-                print(f"{plot=} {unit_id=}")
                 plot.draw(axis, sorting_analysis, unit_id)
 
     # clean up the panels, or else things get clipped