diff --git a/README.md b/README.md
index 60fc076..74a4612 100644
--- a/README.md
+++ b/README.md
@@ -46,8 +46,11 @@ Figure 7 can be reproduce using:
 
 Figure 8 can be reproduce using:
 
-  * TODO Chris
-  * TODO Chris
+  * `cd notebooks/real_data_figure` : move to the folder where the environment is defined
+  * `uv run sort_all_real_data.py` : Sort three datasets using four sorters
+  * `uv run generate_curation_data.py` : Use UnitRefine, Bombcell and SLAy to curate the sorting output
+  * `uv run make_drift_plots.py` : Make the drift and probe plots
+  * `figure.tpy` : Generate the plot
 
 
 Supplementary figure can be reproduce using::
diff --git a/notebooks/real_data_figure/drift_maps_and_probes/Duszkiewicz_probe.png b/notebooks/real_data_figure/drift_maps_and_probes/Duszkiewicz_probe.png
new file mode 100644
index 0000000..42d8e42
Binary files /dev/null and b/notebooks/real_data_figure/drift_maps_and_probes/Duszkiewicz_probe.png differ
diff --git a/notebooks/real_data_figure/drift_maps_and_probes/IBL_probe.png b/notebooks/real_data_figure/drift_maps_and_probes/IBL_probe.png
new file mode 100644
index 0000000..9ea22a5
Binary files /dev/null and b/notebooks/real_data_figure/drift_maps_and_probes/IBL_probe.png differ
diff --git a/notebooks/real_data_figure/drift_maps_and_probes/ucl_probe.png b/notebooks/real_data_figure/drift_maps_and_probes/ucl_probe.png
new file mode 100644
index 0000000..6358233
Binary files /dev/null and b/notebooks/real_data_figure/drift_maps_and_probes/ucl_probe.png differ
diff --git a/notebooks/real_data_figure/figure.typ b/notebooks/real_data_figure/figure.typ
new file mode 100644
index 0000000..472dc7c
--- /dev/null
+++ b/notebooks/real_data_figure/figure.typ
@@ -0,0 +1,106 @@
+#set page("us-letter")
+#set text(size: 9pt, font: "New Computer Modern")
+
+#table(
+  columns: (1fr, 0.087fr, 0.55fr),
+  // Text fits content, SVGs split remaining space
+  align: center,
+  gutter: 0pt,
+  stroke: none,
+
+  box(width: 100%)[
+    #table(
+      columns: (1.3fr, 1fr, 1fr, 1fr, 1fr, 1fr, 1.2fr, 1fr, 1fr),
+      inset: 4pt,
+      stroke: 0.5pt + gray,
+      fill: (x, y) => {
+        if x == 2 or x == 3 {
+          green.lighten(85%)
+        } else if x == 4 or x == 5 or x == 6 {
+          orange.lighten(85%)
+        } else if x == 7 or x == 8 {
+          red.lighten(85%)
+        } else if x == 0 or x == 1 {
+          gray.lighten(80%)
+        }
+      },
+      [Sorter], [Tot units], [BC good], [UR sua], [BC mua], [UR mua], [SLAy merges], [BC noise], [UR noise],
+    )
+  ],
+  align()[],
+  align()[],
+
+  // --- Row 1 ---
+
+  box(width: 100%)[
+    *Lebedeva et. al., Chronic, NP2.0, 38 mins*
+    #table(
+      columns: (1.3fr, 1fr, 1fr, 1fr, 1fr, 1fr, 1.2fr, 1fr, 1fr),
+      inset: 4pt,
+      stroke: 0.5pt + gray,
+      fill: (x, y) => {
+        if x == 2 or x == 3 {
+          green.lighten(85%)
+        } else if x == 4 or x == 5 or x == 6 {
+          orange.lighten(85%)
+        } else if x == 7 or x == 8 {
+          red.lighten(85%)
+        }
+      },
+      [KS4], [808], [246], [288], [489], [312], [5], [73], [208],
+      [Lupin], [683], [216], [259], [460], [338], [2], [7], [86],
+      [TCD2], [617], [119], [246], [491], [369], [13], [7], [2],
+      [SC2], [270], [102], [133], [166], [137], [1], [2], [0],
+    )
+  ],
+  align()[#image("drift_maps_and_probes/ucl_probe.png", height: 12.6%)],
+  image("drift_maps_and_probes/ucl_drift.svg"),
+  // --- Row 2 ---
+  box(width: 100%)[
+    *IBL, Acute, NP1.0, 67 mins*
+    #table(
+      columns: (1.3fr, 1fr, 1fr, 1fr, 1fr, 1fr, 1.2fr, 1fr, 1fr),
+      inset: 4pt,
+      stroke: 0.5pt + gray,
+      fill: (x, y) => {
+        if x == 2 or x == 3 {
+          green.lighten(85%)
+        } else if x == 4 or x == 5 or x == 6 {
+          orange.lighten(85%)
+        } else if x == 7 or x == 8 {
+          red.lighten(85%)
+        }
+      },
+      [KS4], [1050], [210], [459], [673], [354], [24], [167], [237],
+      [Lupin], [864], [209], [379], [601], [278], [6], [54], [207],
+      [TDC2], [954], [124], [417], [778], [504], [33], [52], [33],
+      [SC2], [458], [97], [170], [333], [271], [0], [28], [17],
+    )
+  ],
+  align()[#image("drift_maps_and_probes/IBL_probe.png", height: 11.9%)],
+  image("drift_maps_and_probes/IBL_drift.svg"),
+  // --- Row 3 ---
+  box(width: 100%)[
+    *Duszkiewicz et. al., Chronic, CN 156H5, 211 mins*
+    #table(
+      columns: (1.3fr, 1fr, 1fr, 1fr, 1fr, 1fr, 1.2fr, 1fr, 1fr),
+      inset: 4pt,
+      stroke: 0.5pt + gray,
+      fill: (x, y) => {
+        if x == 2 or x == 3 {
+          green.lighten(85%)
+        } else if x == 4 or x == 5 or x == 6 {
+          orange.lighten(85%)
+        } else if x == 7 or x == 8 {
+          red.lighten(85%)
+        }
+      },
+      [KS4], [174], [41], [71], [98], [68], [2], [35], [35],
+      [Lupin], [162], [56], [96], [103], [63], [4], [3], [3],
+      [TDC2], [191], [11], [60], [180], [128], [9], [0], [3],
+      [SC2], [58], [4], [9], [53], [47], [0], [1], [2],
+    )
+  ],
+  align()[#image("drift_maps_and_probes/Duszkiewicz_probe.png", height: 12.3%)],
+  image("drift_maps_and_probes/Duszkiewicz_drift.svg"),
+)
diff --git a/notebooks/real_data_figure/generate_curation_data.py b/notebooks/real_data_figure/generate_curation_data.py
new file mode 100644
index 0000000..d0998fa
--- /dev/null
+++ b/notebooks/real_data_figure/generate_curation_data.py
@@ -0,0 +1,73 @@
+"""
+Generates curation results from computed analyzers.
+
+Once you have the analyzers, run this code by `cd`ing into the `real_data_figure` 
+folder, then running
+>>> uv run generate_curation_data.py
+
+"""
+
+import spikeinterface.full as si
+import numpy as np
+import pandas as pd
+from pathlib import Path
+
+repo_folder = Path("/home/nolanlab/fromgit/sorting_components_benchmark_paper/")
+real_data_figure_folder = repo_folder / "notebooks/real_data_figure"
+analyzers_folder = real_data_figure_folder / "analyzers"
+
+dataset_protocols = {
+    'IBL': ['kilosort4_motion_correction', 'lupin_motion_correction', 'tridesclous2_motion_correction','spykingcircus2_motion_correction'],
+    'ucl': ['kilosort4_no_motion_correction', 'lupin_no_motion_correction', 'tridesclous2_no_motion_correction','spykingcircus2_no_motion_correction'],
+    'Duszkiewicz': ['kilosort4_no_motion_correction', 'lupin_no_motion_correction', 'tridesclous2_no_motion_correction','spykingcircus2_no_motion_correction'],
+}
+
+bombcell_labels = ['good', 'mua', 'noise', 'non_soma_good', 'non_soma_mua']
+unitrefine_labels = ['sua', 'mua', 'noise']
+merge_presets = ['slay']
+
+for dataset_name, protocols in dataset_protocols.items():
+    bombcell_results = []
+    unitrefine_results = []
+    all_protocols_data = []
+    for protocol in protocols:
+
+        analyzer_path = analyzers_folder / f"{dataset_name}_{protocol}_analyzer"
+        if analyzer_path.is_dir():
+            analyzer = si.load_sorting_analyzer(analyzer_path)
+        else:
+            analyzer = si.load_sorting_analyzer(str(analyzer_path) + '.zarr')
+        
+        bombcell_unit_label = si.bombcell_label_units(analyzer, split_non_somatic_good_mua=True)['bombcell_label'].values
+        bombcell_results = {label: np.sum(bombcell_unit_label == label) for label in bombcell_labels}
+
+        # You need to donwload the UnitRefine models `noise_neural_classifier_lightweight` and `sua_mua_classifier_lightweight` from
+        # https://huggingface.co/AnoushkaJain3
+        unitrefine_unit_label = si.unitrefine_label_units(analyzer, noise_neural_classifier='/home/nolanlab/Downloads/noise_neural_classifier_lightweight', sua_mua_classifier='/home/nolanlab/Downloads/sua_mua_classifier_lightweight')
+        unitrefine_results = {label: np.sum(unitrefine_unit_label['unitrefine_label'] == label) for label in unitrefine_labels}
+
+        merge_results = {merge_preset: len(si.compute_merge_unit_groups(analyzer, preset=merge_preset)) for merge_preset in merge_presets}
+
+        protocol_data = [
+            protocol,
+            analyzer.get_num_units(),
+            bombcell_results['good'] + bombcell_results['non_soma_good'],
+            unitrefine_results['sua'],
+            bombcell_results['mua'] + bombcell_results['non_soma_mua'],
+            unitrefine_results['mua'],
+            merge_results['slay'],
+            bombcell_results['noise'],
+            unitrefine_results['noise'],
+        ]
+
+        all_protocols_data.append(protocol_data)
+
+    results = pd.DataFrame(all_protocols_data, columns=["sorter", "total units", "bombcell good", "unitrefine sua", "bombcell mua", "unitrefine mua", "# slay merges", "bombcell noise", "unitrefine noise"], index=None)
+
+    results.to_csv(real_data_figure_folder / f"curation_results/{dataset_name}_results.csv", index=False)  
+
+    # render for typst rendering
+    for row in results.iterrows():
+        for cell in row[1]:
+            print(f"[{cell}], ", end="")
+        print("")
diff --git a/notebooks/real_data_figure/make_drift_plots.py b/notebooks/real_data_figure/make_drift_plots.py
new file mode 100644
index 0000000..0986fb5
--- /dev/null
+++ b/notebooks/real_data_figure/make_drift_plots.py
@@ -0,0 +1,111 @@
+import matplotlib.pyplot as plt
+import matplotlib.cm as cm
+import matplotlib.colors as mcolors
+import spikeinterface.full as si
+
+from pathlib import Path
+
+repo_folder = Path("/home/nolanlab/fromgit/sorting_components_benchmark_paper/")
+real_data_figure_folder = repo_folder / "notebooks/real_data_figure"
+analyzers_folder = real_data_figure_folder / "analyzers"
+drift_maps_folder = real_data_figure_folder / "drift_maps_and_probes"
+
+bombcell_labels = ['good', 'mua', 'noise', 'non_soma_good', 'non_soma_mua']
+
+protocol = 'no_motion_correction'
+
+FONT_SIZE = 18
+
+plotting_settings = {
+    'ucl': {
+        'protocol': 'no_motion_correction',
+        'vmin': -600,
+        'scatter_decimate': 20,
+        'cbar_ticks': [-600,-500,-400,-300,-200,-100,0],
+        'cbar_ticklabels': ['600','','400','','200','','0'],
+        'yticklabels': ['','2.9', '', '3.1', '', '3.3', '', '3.5'],
+        'xticks_s': [0,600,1200,1800],
+    },
+    'IBL': {
+        'protocol': 'motion_correction',
+        'vmin': -457.829994,
+        'scatter_decimate': 20,
+        'cbar_ticks': [-600,-500,-400,-300,-200,-100,0], 
+        'cbar_ticklabels': ['600','','400','','200','','0'],
+        'yticklabels': ['', '0', '', '1', '', '2', '', '3', ''],
+        'xticks_s': [0,900,1800,2700,3600],
+    },
+    'Duszkiewicz': {
+        'protocol': 'no_motion_correction',
+        'vmin': -380,
+        'scatter_decimate': 5,
+        'cbar_ticks': [-400,-300,-200,-100,0],
+        'cbar_ticklabels': [400,300,200,100,0],
+        'yticklabels': ['', '', '0.2', '', '0.4', '', '0.6', '', '0.8'],
+        'xticks_s': [0,3000,6000,9000,12000],
+    }
+}
+
+for dataset_name, dataset_settings in plotting_settings.items():
+   
+    protocol = dataset_settings['protocol']
+   
+    analyzer_path =  analyzers_folder / f'{dataset_name}_kilosort4_{protocol}_analyzer'
+    if analyzer_path.is_dir():
+        analyzer = si.load_sorting_analyzer(analyzer_path)
+    else:
+        analyzer = si.load_sorting_analyzer(str(analyzer_path) + '.zarr')
+
+    print(analyzer.get_total_duration())
+        
+    bombcell_unit_labels = si.bombcell_label_units(analyzer, split_non_somatic_good_mua=True)['bombcell_label'].values
+    good_units = analyzer.unit_ids[bombcell_unit_labels == 'good']
+    analyzer_good = analyzer.select_units(good_units)
+    
+    cmap_name = 'inferno'
+    
+    fig = si.plot_drift_raster_map(
+        sorting_analyzer=analyzer_good, 
+        cmap=cmap_name, 
+        alpha=0.10,
+        scatter_decimate=dataset_settings['scatter_decimate'], 
+        figsize=(8,4.5)
+    )
+    # 1. Define your parameters
+    vmin = dataset_settings['vmin']
+    vmax = 0
+    
+    # 2. Create the Normalization and Mappable objects
+    norm = mcolors.Normalize(vmin=vmin, vmax=vmax)
+    sm = cm.ScalarMappable(cmap=plt.get_cmap(cmap_name), norm=norm)
+    sm.set_array([]) # Required for the colorbar to initialize correctly
+    
+    # 3. Access your existing figure/axes and add the colorbar
+    # Assuming 'fig' is your figure object
+    ax = fig.figure.get_axes()[0] 
+    cbar = fig.figure.colorbar(sm, ax=ax)
+    
+    # 2. Find the scatter plot and rasterize it
+    # In Matplotlib, scatter plots are usually 'PathCollection' objects
+    for artist in ax.get_children():
+        if isinstance(artist, plt.matplotlib.collections.PathCollection):
+            artist.set_rasterized(True)
+    
+    # 4. (Optional) Add a label
+    cbar_ticks = dataset_settings['cbar_ticks']
+    #cbar_ticks = [-80,-60,-40,-20,0]
+    cbar.set_label('Abs peak amplitude [uV]', fontsize=FONT_SIZE)
+    cbar.set_ticklabels(dataset_settings['cbar_ticklabels']) 
+    cbar.ax.tick_params(labelsize=FONT_SIZE) # Font size for colorbar ticks
+    
+    ax.set_ylabel('Depth [mm]', fontsize=FONT_SIZE)
+    ax.set_yticklabels(dataset_settings['yticklabels'], fontsize=FONT_SIZE)
+    
+    xticks_s = dataset_settings['xticks_s']
+    ax.set_xticks(xticks_s)
+    ax.set_xticklabels([int(xtick/60) for xtick in xticks_s], fontsize=FONT_SIZE)
+    ax.set_xlabel('Time [min]', fontsize=FONT_SIZE)
+    
+    ax.set_title(label=None)
+    
+    fig.figure.savefig(drift_maps_folder / f'{dataset_name}_drift.svg',  bbox_inches='tight')
diff --git a/notebooks/real_data_figure/pyproject.toml b/notebooks/real_data_figure/pyproject.toml
new file mode 100644
index 0000000..e80d7e0
--- /dev/null
+++ b/notebooks/real_data_figure/pyproject.toml
@@ -0,0 +1,24 @@
+[project]
+name = "real-data-figure"
+version = "0.1.0"
+description = "Add your description here"
+readme = "README.md"
+requires-python = ">=3.13"
+dependencies = [
+    "h5py>=3.16.0",
+    "hdbscan>=0.8.42",
+    "ipykernel>=7.2.0",
+    "kilosort==4.1.2",
+    "matplotlib>=3.10.9",
+    "numba>=0.65.1",
+    "pandas>=3.0.2",
+    "scikit-learn==1.6",
+    "scipy>=1.17.1",
+    "skops>=0.14.0",
+    "spikeinterface>=0.104.3",
+]
+
+[dependency-groups]
+dev = [
+    "ruff>=0.15.12",
+]
diff --git a/notebooks/real_data_figure/sort_all_real_data.py b/notebooks/real_data_figure/sort_all_real_data.py
new file mode 100644
index 0000000..6a57860
--- /dev/null
+++ b/notebooks/real_data_figure/sort_all_real_data.py
@@ -0,0 +1,107 @@
+"""
+Based on code from https://github.com/MattNolanLab/nolanlab-ephys
+
+The IBL data is the file
+    sub-UCLA034_ses-3537d970-f515-4786-853f-23de525e110f_desc-raw_ecephys.nwb
+from DANDI dataset 000409 - https://dandiarchive.org/dandiset/000409
+
+The UCL data is the file
+    AL032_2020-01-07
+from https://rdr.ucl.ac.uk/articles/dataset/Chronic_recordings_from_Neuropixels_2_0_probes_in_mice/24411841
+This needs to be untarred into a folder.
+
+The Duszkiewicz data is the file
+    sub-A3702/sub-A3702_ses-191126_behavior+ecephys.nwb
+from the DANDI dataset 000939 - https://dandiarchive.org/dandiset/000939
+
+For this code to run, put all files in the "raw_data" folder and change the "repo_folder" below, to point at
+your local copy of this repo. Your file organisation should look like
+
+sorting_components_benchmark_paper/   <-- `repo_folder` points here
+   notebooks/
+       real_data_figure/
+           sort_all_real_data.py
+           raw_data/
+               sub-UCLA034_ses-3537d970-f515-4786-853f-23de525e110f_desc-raw_ecephys.nwb
+               sub-A3702_ses-191126_behavior+ecephys.nwb
+               AL032_2020-01-07/
+                   ...
+            analyzers/
+            curation_results/
+
+`cd` to the `real_data_figure` folder and run
+>>> uv run sort_all_real_data.py
+
+After this has run the `analyzers` folder should contain the following 12 analyzers:
+
+analyzers/
+    Duszkiewicz_kilosort4_no_motion_correction_analyzer
+    Duszkiewicz_lupin_no_motion_correction_analyzer
+    Duszkiewicz_spykingcircus2_no_motion_correction_analyzer
+    Duszkiewicz_tridesclous2_no_motion_correction_analyzer
+    IBL_kilosort4_motion_correction_analyzer
+    IBL_lupin_motion_correction_analyzer
+    IBL_spykingcircus2_motion_correction_analyzer
+    IBL_tridesclous2_motion_correction_analyzer
+    ucl_kilosort4_no_motion_correction_analyzer
+    ucl_lupin_no_motion_correction_analyzer
+    ucl_spykingcircus2_no_motion_correction_analyzer
+    ucl_tridesclous2_no_motion_correction_analyzer 
+
+"""
+import spikeinterface.full as si
+from pathlib import Path
+from sort_one_piece_of_data import do_sorting
+
+# edit this to point to the `sorting_components_benchmark_paper` on your own computer
+repo_folder = Path("/home/nolanlab/fromgit/sorting_components_benchmark_paper/")
+
+raw_data_folder = repo_folder / "notebooks/real_data_figure/raw_data"
+analyzers_folder = repo_folder / "notebooks/real_data_figure/raw_data"
+
+# NP1 data from IBL
+np1_protocols = [
+        'kilosort4_motion_correction',
+        'lupin_motion_correction',
+        'spykingcircus2_motion_correction'
+        'tridesclous2A_motion_correction'
+    ]
+np1_data_folder = raw_data_folder / 'sub-UCLA034_ses-3537d970-f515-4786-853f-23de525e110f_desc-raw_ecephys.nwb'
+np1_analyzer_folders = [
+    analyzers_folder / f"IBL_{protocol}_analyzer" for protocol in np1_protocols
+]
+
+recording = si.read_nwb_recording(np1_data_folder, electrical_series_path = 'acquisition/ElectricalSeriesProbe00AP')
+for protocol_name, analyzer_folder in zip(np1_protocols, np1_analyzer_folders):
+    do_sorting(recording, analyzer_folder, protocol_name)
+
+# NP2 data from UCL
+np2_protocols = [
+        'kilosort4_no_motion_correction',
+        'lupin_no_motion_correction',
+        'spykingcircus2_no_motion_correction'
+        'tridesclous2A_no_motion_correction'
+    ]
+np2_data_folder = raw_data_folder / 'AL032_2020-01-07'
+np2_analyzer_folders = [
+    analyzers_folder / f"ucl_{protocol}_analyzer" for protocol in np1_protocols
+]
+
+recording = si.read_cbin_ibl(np2_data_folder)
+for protocol_name, analyzer_folder in zip(np2_protocols, np2_analyzer_folders):
+    do_sorting(recording, analyzer_folder, protocol_name)
+
+# CN data from  Duszkiewicz
+cn_protocols = [
+        'kilosort4_no_motion_correction',
+        'lupin_no_motion_correction',
+        'spykingcircus2_no_motion_correction'
+        'tridesclous2A_no_motion_correction'
+    ]
+cn_analyzer_folders = [
+    analyzers_folder / f"Duszkiewicz_{protocol}_analyzer" for protocol in cn_protocols
+]
+cn_data_folder = raw_data_folder / 'sub-A3702_ses-191126_behavior+ecephys.nwb'
+recording = si.read_nwb_recording(cn_data_folder, electrical_series_path = 'acquisition/ElectricalSeries')
+for protocol_name, analyzer_folder in zip(cn_protocols, cn_analyzer_folders):
+    do_sorting(recording, analyzer_folder, protocol_name)
diff --git a/notebooks/real_data_figure/sort_one_piece_of_data.py b/notebooks/real_data_figure/sort_one_piece_of_data.py
new file mode 100644
index 0000000..4ea0dc8
--- /dev/null
+++ b/notebooks/real_data_figure/sort_one_piece_of_data.py
@@ -0,0 +1,149 @@
+"""
+Based on code from https://github.com/MattNolanLab/nolanlab-ephys
+"""
+import spikeinterface.full as si
+
+protocols = {
+    "kilosort4_no_motion_correction": {
+        "preprocessing": {
+        },
+        "sorting": {
+            "sorter_name": "kilosort4",
+            "do_correction": False,
+            "use_binary_file": False,
+        },
+        "preprocessing_for_analyzer": {
+            "common_reference": {},
+            "bandpass_filter": {},
+        },
+    },
+    "spykingcircus2_no_motion_correction": {
+        "preprocessing": {},
+        "sorting": {
+            "sorter_name": "spykingcircus2",
+            "apply_motion_correction": False,
+            "cache_preprocessing": {"mode": "folder", "folder": "sk2_pre"},
+        },
+        "preprocessing_for_analyzer": {
+            "bandpass_filter": {},
+            "common_reference": {},
+        },
+    },
+    "tridesclous2A_no_motion_correction": {
+        "preprocessing": {},
+        "sorting": {
+            "sorter_name": "tridesclous2",
+            "cache_preprocessing_mode": "folder",
+        },
+        "preprocessing_for_analyzer": {
+            "bandpass_filter": {},
+            "common_reference": {},
+        },
+    },
+    "lupin_no_motion_correction": {
+        "preprocessing": {},
+        "sorting": {
+            "sorter_name": "lupin",
+            "apply_motion_correction": False,
+        },
+        "preprocessing_for_analyzer": {
+            "bandpass_filter": {},
+            "common_reference": {},
+        },
+    },
+    "kilosort4_motion_correction": {
+        "preprocessing": {
+        },
+        "sorting": {
+            "sorter_name": "kilosort4",
+            "do_correction": True,
+            "use_binary_file": False,
+        },
+        "preprocessing_for_analyzer": {
+            "common_reference": {},
+            "bandpass_filter": {},
+        },
+    },
+    "spykingcircus2_motion_correction": {
+        "preprocessing": {},
+        "sorting": {
+            "sorter_name": "spykingcircus2",
+            "apply_motion_correction": True,
+            "cache_preprocessing": {"mode": "folder", "folder": "sk2_pre"},
+        },
+        "preprocessing_for_analyzer": {
+            "bandpass_filter": {},
+            "common_reference": {},
+        },
+    },
+    "tridesclous2A_motion_correction": {
+        "preprocessing": {},
+        "sorting": {
+            "sorter_name": "tridesclous2",
+            "cache_preprocessing_mode": "folder",
+            "apply_motion_correction": True,
+        },
+        "preprocessing_for_analyzer": {
+            "bandpass_filter": {},
+            "common_reference": {},
+        },
+    },
+    "lupin_motion_correction": {
+        "preprocessing": {},
+        "sorting": {
+            "sorter_name": "lupin",
+            "apply_motion_correction": True,
+        },
+        "preprocessing_for_analyzer": {
+            "bandpass_filter": {},
+            "common_reference": {},
+        },
+    },
+}
+
+postprocessing_extensions_to_compute = {
+    "unit_locations": {},
+    "random_spikes": {},
+    "noise_levels": {},
+    "waveforms": {},
+    "templates": {},
+    "spike_amplitudes": {"peak_sign": "both"},
+    "amplitude_scalings": {},
+    "isi_histograms": {},
+    "spike_locations": {"peak_sign": "both"},
+    "correlograms": {},
+    "template_similarity": {"method": "l2"},
+    "quality_metrics": {},
+    "template_metrics": {},
+}
+
+def do_sorting(recording, analyzer_path, protocol_name, n_jobs=8):
+
+    si.set_global_job_kwargs(n_jobs=n_jobs)
+
+    protocol_info = protocols[protocol_name]
+
+    pp_recording = si.apply_preprocessing_pipeline(
+        recording, protocol_info["preprocessing"]
+    )
+    sorting = si.run_sorter(
+        recording=pp_recording,
+        **protocol_info["sorting"],
+        remove_existing_folder=True,
+        verbose=True,
+    )
+
+    preprocessed_recording_for_analyzer = si.apply_preprocessing_pipeline(
+        recording, protocol_info["preprocessing_for_analyzer"]
+    )
+
+    analyzer = si.create_sorting_analyzer(
+        recording=preprocessed_recording_for_analyzer,
+        sorting=sorting,
+        folder=analyzer_path,
+        format="binary_folder",
+        peak_sign="both",
+        radius_um=70,
+    )
+
+    analyzer.compute(postprocessing_extensions_to_compute)