sharp wave ripple detection

%23%20%2F%2F%2F%20script%0A%23%20requires-python%20%3D%20%22%3E%3D3.11%22%0A%23%20dependencies%20%3D%20%5B%0A%23%20%20%20%20%20%22marimo%3E%3D0.23.0%22%2C%0A%23%20%5D%0A%23%20%2F%2F%2F%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.23.6%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20html%0A%20%20%20%20import%20math%0A%20%20%20%20import%20random%0A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20html%2C%20math%2C%20mo%2C%20random%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%20Interactive%20simulated%20SWR%20detector%0A%0A%20%20%20%20This%20browser-only%20demo%20simulates%20paired%20ripple%20and%20sharp-wave%20LFP%20channels%2C%0A%20%20%20%20then%20runs%20a%20compact%20joint%20detector.%20Green%20spans%20mark%20the%20known%20simulated%20SWR%0A%20%20%20%20windows%3B%20orange%20spans%20mark%20accepted%20detections.%0A%0A%20%20%20%20This%20app%20intentionally%20avoids%20%60neuro_py%60%2C%20%60nelpy%60%2C%20SciPy%2C%20Matplotlib%2C%20and%20real%0A%20%20%20%20session%20data%20so%20the%20interactive%20example%20can%20load%20quickly%20in%20browser-WASM.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20ripple_amplitude%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.5%2C%204.0%2C%20step%3D0.1%2C%20value%3D2.0%2C%20label%3D%22Ripple%20amplitude%22%0A%20%20%20%20)%0A%20%20%20%20sharp_wave_amplitude%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.5%2C%205.0%2C%20step%3D0.1%2C%20value%3D2.8%2C%20label%3D%22Sharp-wave%20amplitude%22%0A%20%20%20%20)%0A%20%20%20%20noise_amplitude%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.2%2C%202.0%2C%20step%3D0.1%2C%20value%3D1.1%2C%20label%3D%22Noise%20amplitude%22%0A%20%20%20%20)%0A%20%20%20%20ripple_duration%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.050%2C%200.150%2C%20step%3D0.005%2C%20value%3D0.100%2C%20label%3D%22Ripple%20duration%20(s)%22%0A%20%20%20%20)%0A%20%20%20%20ripple_threshold%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.5%2C%204.0%2C%20step%3D0.1%2C%20value%3D1.8%2C%20label%3D%22Ripple%20high%20threshold%20(z)%22%0A%20%20%20%20)%0A%20%20%20%20sharp_wave_threshold%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.5%2C%204.0%2C%20step%3D0.1%2C%20value%3D1.4%2C%20label%3D%22Sharp-wave%20high%20threshold%20(z)%22%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Simulation%20and%20detector%20controls%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5Bripple_amplitude%2C%20sharp_wave_amplitude%2C%20noise_amplitude%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20justify%3D%22start%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gap%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5Bripple_duration%2C%20ripple_threshold%2C%20sharp_wave_threshold%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20justify%3D%22start%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gap%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20noise_amplitude%2C%0A%20%20%20%20%20%20%20%20ripple_amplitude%2C%0A%20%20%20%20%20%20%20%20ripple_duration%2C%0A%20%20%20%20%20%20%20%20ripple_threshold%2C%0A%20%20%20%20%20%20%20%20sharp_wave_amplitude%2C%0A%20%20%20%20%20%20%20%20sharp_wave_threshold%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell%0Adef%20_(math%2C%20random)%3A%0A%20%20%20%20FS%20%3D%201250.0%0A%20%20%20%20DURATION%20%3D%201.5%0A%20%20%20%20RIPPLE_FREQUENCY%20%3D%20165.0%0A%20%20%20%20SWR_TIMES%20%3D%20%5B0.45%2C%201.02%5D%0A%0A%20%20%20%20def%20zscore(values)%3A%0A%20%20%20%20%20%20%20%20mean%20%3D%20sum(values)%20%2F%20len(values)%0A%20%20%20%20%20%20%20%20var%20%3D%20sum((value%20-%20mean)%20**%202%20for%20value%20in%20values)%20%2F%20len(values)%0A%20%20%20%20%20%20%20%20scale%20%3D%20math.sqrt(var)%0A%20%20%20%20%20%20%20%20if%20scale%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%5B0.0%20for%20_%20in%20values%5D%0A%20%20%20%20%20%20%20%20return%20%5B(value%20-%20mean)%20%2F%20scale%20for%20value%20in%20values%5D%0A%0A%20%20%20%20def%20brown_noise(n_samples%2C%20seed)%3A%0A%20%20%20%20%20%20%20%20rng%20%3D%20random.Random(seed)%0A%20%20%20%20%20%20%20%20value%20%3D%200.0%0A%20%20%20%20%20%20%20%20values%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_%20in%20range(n_samples)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20value%20%2B%3D%20rng.gauss(0.0%2C%201.0)%0A%20%20%20%20%20%20%20%20%20%20%20%20values.append(value)%0A%20%20%20%20%20%20%20%20return%20zscore(values)%0A%0A%20%20%20%20def%20moving_average(values%2C%20window)%3A%0A%20%20%20%20%20%20%20%20window%20%3D%20max(1%2C%20int(window))%0A%20%20%20%20%20%20%20%20half%20%3D%20window%20%2F%2F%202%0A%20%20%20%20%20%20%20%20averaged%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20index%20in%20range(len(values))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20start%20%3D%20max(0%2C%20index%20-%20half)%0A%20%20%20%20%20%20%20%20%20%20%20%20stop%20%3D%20min(len(values)%2C%20index%20%2B%20half%20%2B%201)%0A%20%20%20%20%20%20%20%20%20%20%20%20averaged.append(sum(values%5Bstart%3Astop%5D)%20%2F%20(stop%20-%20start))%0A%20%20%20%20%20%20%20%20return%20averaged%0A%0A%20%20%20%20def%20threshold_bounds(values%2C%20low_threshold)%3A%0A%20%20%20%20%20%20%20%20bounds%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20start%20%3D%20None%0A%20%20%20%20%20%20%20%20for%20index%2C%20value%20in%20enumerate(values)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20value%20%3E%3D%20low_threshold%20and%20start%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20start%20%3D%20index%0A%20%20%20%20%20%20%20%20%20%20%20%20elif%20value%20%3C%20low_threshold%20and%20start%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20bounds.append((start%2C%20index%20-%201))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20start%20%3D%20None%0A%20%20%20%20%20%20%20%20if%20start%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20bounds.append((start%2C%20len(values)%20-%201))%0A%20%20%20%20%20%20%20%20return%20bounds%0A%0A%20%20%20%20def%20merge_bounds(bounds%2C%20gap_samples)%3A%0A%20%20%20%20%20%20%20%20if%20not%20bounds%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%5B%5D%0A%20%20%20%20%20%20%20%20merged%20%3D%20%5Bbounds%5B0%5D%5D%0A%20%20%20%20%20%20%20%20for%20start%2C%20stop%20in%20bounds%5B1%3A%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20prev_start%2C%20prev_stop%20%3D%20merged%5B-1%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20start%20-%20prev_stop%20%3C%3D%20gap_samples%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20merged%5B-1%5D%20%3D%20(prev_start%2C%20max(prev_stop%2C%20stop))%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20merged.append((start%2C%20stop))%0A%20%20%20%20%20%20%20%20return%20merged%0A%0A%20%20%20%20def%20bound_containing_index(bounds%2C%20index)%3A%0A%20%20%20%20%20%20%20%20for%20start%2C%20stop%20in%20bounds%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20start%20%3C%3D%20index%20%3C%3D%20stop%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20return%20start%2C%20stop%0A%20%20%20%20%20%20%20%20return%20None%0A%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20DURATION%2C%0A%20%20%20%20%20%20%20%20FS%2C%0A%20%20%20%20%20%20%20%20RIPPLE_FREQUENCY%2C%0A%20%20%20%20%20%20%20%20SWR_TIMES%2C%0A%20%20%20%20%20%20%20%20bound_containing_index%2C%0A%20%20%20%20%20%20%20%20brown_noise%2C%0A%20%20%20%20%20%20%20%20merge_bounds%2C%0A%20%20%20%20%20%20%20%20moving_average%2C%0A%20%20%20%20%20%20%20%20threshold_bounds%2C%0A%20%20%20%20%20%20%20%20zscore%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20DURATION%2C%0A%20%20%20%20FS%2C%0A%20%20%20%20RIPPLE_FREQUENCY%2C%0A%20%20%20%20SWR_TIMES%2C%0A%20%20%20%20brown_noise%2C%0A%20%20%20%20math%2C%0A%20%20%20%20noise_amplitude%2C%0A%20%20%20%20ripple_amplitude%2C%0A%20%20%20%20ripple_duration%2C%0A%20%20%20%20sharp_wave_amplitude%2C%0A)%3A%0A%20%20%20%20timestamps%20%3D%20%5Bsample%20%2F%20FS%20for%20sample%20in%20range(int(DURATION%20*%20FS))%5D%0A%20%20%20%20ripple_noise%20%3D%20brown_noise(len(timestamps)%2C%20seed%3D7)%0A%20%20%20%20sharp_noise%20%3D%20brown_noise(len(timestamps)%2C%20seed%3D8)%0A%20%20%20%20ripple_raw%20%3D%20%5Bnoise_amplitude.value%20*%20value%20for%20value%20in%20ripple_noise%5D%0A%20%20%20%20sharp_wave_raw%20%3D%20%5B%0A%20%20%20%20%20%20%20%200.75%20*%20noise_amplitude.value%20*%20value%20for%20value%20in%20sharp_noise%0A%20%20%20%20%5D%0A%0A%20%20%20%20sharp_wave_duration%20%3D%200.140%0A%20%20%20%20for%20center%20in%20SWR_TIMES%3A%0A%20%20%20%20%20%20%20%20for%20index%2C%20time%20in%20enumerate(timestamps)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20ripple_width%20%3D%20ripple_duration.value%20%2F%206.0%0A%20%20%20%20%20%20%20%20%20%20%20%20sharp_width%20%3D%20sharp_wave_duration%20%2F%206.0%0A%20%20%20%20%20%20%20%20%20%20%20%20ripple_envelope%20%3D%20math.exp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20-((time%20-%20center)%20**%202)%20%2F%20(2.0%20*%20ripple_width**2)%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20sharp_wave_envelope%20%3D%20math.exp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20-((time%20-%20center)%20**%202)%20%2F%20(2.0%20*%20sharp_width**2)%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20ripple_carrier%20%3D%20math.sin(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%202.0%20*%20math.pi%20*%20RIPPLE_FREQUENCY%20*%20time%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20ripple_raw%5Bindex%5D%20%2B%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ripple_amplitude.value%20*%20ripple_carrier%20*%20ripple_envelope%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20sharp_wave_raw%5Bindex%5D%20-%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20sharp_wave_amplitude.value%20*%20sharp_wave_envelope%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20truth%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20center%20-%20ripple_duration.value%20%2F%202.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20center%20%2B%20ripple_duration.value%20%2F%202.0%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20for%20center%20in%20SWR_TIMES%0A%20%20%20%20%5D%0A%20%20%20%20return%20ripple_raw%2C%20sharp_wave_raw%2C%20timestamps%2C%20truth%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20FS%2C%0A%20%20%20%20bound_containing_index%2C%0A%20%20%20%20merge_bounds%2C%0A%20%20%20%20moving_average%2C%0A%20%20%20%20ripple_raw%2C%0A%20%20%20%20ripple_threshold%2C%0A%20%20%20%20sharp_wave_raw%2C%0A%20%20%20%20sharp_wave_threshold%2C%0A%20%20%20%20threshold_bounds%2C%0A%20%20%20%20timestamps%2C%0A%20%20%20%20zscore%2C%0A)%3A%0A%20%20%20%20slow_ripple%20%3D%20moving_average(ripple_raw%2C%20window%3Dint(0.030%20*%20FS))%0A%20%20%20%20ripple_high%20%3D%20%5B%0A%20%20%20%20%20%20%20%20abs(raw%20-%20slow)%20for%20raw%2C%20slow%20in%20zip(ripple_raw%2C%20slow_ripple)%0A%20%20%20%20%5D%0A%20%20%20%20ripple_feature%20%3D%20zscore(%0A%20%20%20%20%20%20%20%20moving_average(ripple_high%2C%20window%3Dint(0.012%20*%20FS))%0A%20%20%20%20)%0A%0A%20%20%20%20sharp_wave_feature%20%3D%20zscore(%0A%20%20%20%20%20%20%20%20moving_average(%5B-value%20for%20value%20in%20sharp_wave_raw%5D%2C%20int(0.040%20*%20FS))%0A%20%20%20%20)%0A%0A%20%20%20%20ripple_bounds%20%3D%20merge_bounds(%0A%20%20%20%20%20%20%20%20threshold_bounds(ripple_feature%2C%20ripple_threshold.value%20*%200.45)%2C%0A%20%20%20%20%20%20%20%20gap_samples%3Dint(0.010%20*%20FS)%2C%0A%20%20%20%20)%0A%20%20%20%20sharp_wave_bounds%20%3D%20merge_bounds(%0A%20%20%20%20%20%20%20%20threshold_bounds(sharp_wave_feature%2C%20sharp_wave_threshold.value%20*%200.45)%2C%0A%20%20%20%20%20%20%20%20gap_samples%3Dint(0.010%20*%20FS)%2C%0A%20%20%20%20)%0A%0A%20%20%20%20search_radius%20%3D%20int(0.060%20*%20FS)%0A%20%20%20%20detected_events%20%3D%20%5B%5D%0A%20%20%20%20for%20ripple_start%2C%20ripple_stop%20in%20ripple_bounds%3A%0A%20%20%20%20%20%20%20%20ripple_peak%20%3D%20max(%0A%20%20%20%20%20%20%20%20%20%20%20%20range(ripple_start%2C%20ripple_stop%20%2B%201)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20key%3Dlambda%20index%3A%20ripple_feature%5Bindex%5D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20if%20ripple_feature%5Bripple_peak%5D%20%3C%20ripple_threshold.value%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%0A%20%20%20%20%20%20%20%20search_start%20%3D%20max(0%2C%20ripple_peak%20-%20search_radius)%0A%20%20%20%20%20%20%20%20search_stop%20%3D%20min(len(timestamps)%2C%20ripple_peak%20%2B%20search_radius%20%2B%201)%0A%20%20%20%20%20%20%20%20sharp_peak%20%3D%20max(%0A%20%20%20%20%20%20%20%20%20%20%20%20range(search_start%2C%20search_stop)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20key%3Dlambda%20index%3A%20sharp_wave_feature%5Bindex%5D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20if%20sharp_wave_feature%5Bsharp_peak%5D%20%3C%20sharp_wave_threshold.value%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%0A%20%20%20%20%20%20%20%20sharp_interval%20%3D%20bound_containing_index(sharp_wave_bounds%2C%20sharp_peak)%0A%20%20%20%20%20%20%20%20if%20sharp_interval%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%0A%20%20%20%20%20%20%20%20sharp_start%2C%20sharp_stop%20%3D%20sharp_interval%0A%20%20%20%20%20%20%20%20event_start%20%3D%20min(ripple_start%2C%20sharp_start)%0A%20%20%20%20%20%20%20%20event_stop%20%3D%20max(ripple_stop%2C%20sharp_stop)%0A%20%20%20%20%20%20%20%20detected_events.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22start%22%3A%20timestamps%5Bevent_start%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22stop%22%3A%20timestamps%5Bmin(event_stop%20%2B%201%2C%20len(timestamps)%20-%201)%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22peak%22%3A%20timestamps%5Bripple_peak%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22ripple_z%22%3A%20ripple_feature%5Bripple_peak%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22sharp_wave_z%22%3A%20sharp_wave_feature%5Bsharp_peak%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20return%20detected_events%2C%20ripple_feature%2C%20sharp_wave_feature%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20DURATION%2C%0A%20%20%20%20detected_events%2C%0A%20%20%20%20html%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20ripple_feature%2C%0A%20%20%20%20ripple_raw%2C%0A%20%20%20%20sharp_wave_feature%2C%0A%20%20%20%20sharp_wave_raw%2C%0A%20%20%20%20timestamps%2C%0A%20%20%20%20truth%2C%0A)%3A%0A%20%20%20%20def%20polyline(values%2C%20row_top%2C%20row_height%2C%20color%2C%20stroke_width%3D1.4)%3A%0A%20%20%20%20%20%20%20%20lo%20%3D%20min(values)%0A%20%20%20%20%20%20%20%20hi%20%3D%20max(values)%0A%20%20%20%20%20%20%20%20if%20hi%20%3D%3D%20lo%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20hi%20%3D%20lo%20%2B%201.0%0A%20%20%20%20%20%20%20%20points%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20time%2C%20value%20in%20zip(timestamps%5B%3A%3A2%5D%2C%20values%5B%3A%3A2%5D)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20x%20%3D%2070.0%20%2B%20(time%20%2F%20DURATION)%20*%20860.0%0A%20%20%20%20%20%20%20%20%20%20%20%20y%20%3D%20row_top%20%2B%20row_height%20-%20((value%20-%20lo)%20%2F%20(hi%20-%20lo))%20*%20row_height%0A%20%20%20%20%20%20%20%20%20%20%20%20points.append(f%22%7Bx%3A.1f%7D%2C%7By%3A.1f%7D%22)%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20f'%3Cpolyline%20points%3D%22%7B%22%20%22.join(points)%7D%22%20fill%3D%22none%22%20'%0A%20%20%20%20%20%20%20%20%20%20%20%20f'stroke%3D%22%7Bcolor%7D%22%20stroke-width%3D%22%7Bstroke_width%7D%22%20%2F%3E'%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20def%20span_rect(start%2C%20stop%2C%20row_top%2C%20row_height%2C%20color%2C%20opacity)%3A%0A%20%20%20%20%20%20%20%20x%20%3D%2070.0%20%2B%20(start%20%2F%20DURATION)%20*%20860.0%0A%20%20%20%20%20%20%20%20width%20%3D%20max(1.0%2C%20((stop%20-%20start)%20%2F%20DURATION)%20*%20860.0)%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20f'%3Crect%20x%3D%22%7Bx%3A.1f%7D%22%20y%3D%22%7Brow_top%3A.1f%7D%22%20width%3D%22%7Bwidth%3A.1f%7D%22%20'%0A%20%20%20%20%20%20%20%20%20%20%20%20f'height%3D%22%7Brow_height%3A.1f%7D%22%20fill%3D%22%7Bcolor%7D%22%20opacity%3D%22%7Bopacity%7D%22%20%2F%3E'%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20def%20row_label(text%2C%20row_top)%3A%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20f'%3Ctext%20x%3D%2218%22%20y%3D%22%7Brow_top%20%2B%2018%3A.1f%7D%22%20font-size%3D%2213%22%20'%0A%20%20%20%20%20%20%20%20%20%20%20%20f'font-weight%3D%22700%22%20fill%3D%22%23334155%22%3E%7Bhtml.escape(text)%7D%3C%2Ftext%3E'%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20rows%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(%22Ripple%20LFP%22%2C%20ripple_raw%2C%2028%2C%20145%2C%20%22%23111827%22)%2C%0A%20%20%20%20%20%20%20%20(%22Sharp-wave%20LFP%22%2C%20sharp_wave_raw%2C%20198%2C%20145%2C%20%22%23047857%22)%2C%0A%20%20%20%20%20%20%20%20(%22Detector%20features%22%2C%20ripple_feature%2C%20368%2C%20145%2C%20%22%23991b1b%22)%2C%0A%20%20%20%20%5D%0A%20%20%20%20overlays%20%3D%20%5B%5D%0A%20%20%20%20for%20row_top%20in%20%5B28%2C%20198%2C%20368%5D%3A%0A%20%20%20%20%20%20%20%20overlays.extend(%0A%20%20%20%20%20%20%20%20%20%20%20%20span_rect(start%2C%20stop%2C%20row_top%2C%20145%2C%20%22%2316a34a%22%2C%20%220.16%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20start%2C%20stop%20in%20truth%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20overlays.extend(%0A%20%20%20%20%20%20%20%20%20%20%20%20span_rect(event%5B%22start%22%5D%2C%20event%5B%22stop%22%5D%2C%20row_top%2C%20145%2C%20%22%23f59e0b%22%2C%20%220.28%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20event%20in%20detected_events%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20backgrounds%20%3D%20%5B%5D%0A%20%20%20%20series%20%3D%20%5B%5D%0A%20%20%20%20for%20label%2C%20values%2C%20row_top%2C%20row_height%2C%20color%20in%20rows%3A%0A%20%20%20%20%20%20%20%20backgrounds.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20f'%3Crect%20x%3D%2270%22%20y%3D%22%7Brow_top%7D%22%20width%3D%22860%22%20height%3D%22%7Brow_height%7D%22%20'%0A%20%20%20%20%20%20%20%20%20%20%20%20'rx%3D%228%22%20fill%3D%22%23f8fafc%22%20stroke%3D%22%23e2e8f0%22%20%2F%3E'%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20series.append(row_label(label%2C%20row_top))%0A%20%20%20%20%20%20%20%20series.append(polyline(values%2C%20row_top%20%2B%208%2C%20row_height%20-%2016%2C%20color))%0A%20%20%20%20series.append(%0A%20%20%20%20%20%20%20%20polyline(sharp_wave_feature%2C%20376%2C%20129%2C%20%22%23047857%22%2C%20stroke_width%3D1.1)%0A%20%20%20%20)%0A%0A%20%20%20%20ticks%20%3D%20%5B%5D%0A%20%20%20%20for%20tick%20in%20%5B0.0%2C%200.5%2C%201.0%2C%201.5%5D%3A%0A%20%20%20%20%20%20%20%20x%20%3D%2070.0%20%2B%20(tick%20%2F%20DURATION)%20*%20860.0%0A%20%20%20%20%20%20%20%20ticks.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20f'%3Cline%20x1%3D%22%7Bx%3A.1f%7D%22%20x2%3D%22%7Bx%3A.1f%7D%22%20y1%3D%22520%22%20y2%3D%22526%22%20'%0A%20%20%20%20%20%20%20%20%20%20%20%20'stroke%3D%22%2364748b%22%20%2F%3E'%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20ticks.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20f'%3Ctext%20x%3D%22%7Bx%3A.1f%7D%22%20y%3D%22546%22%20text-anchor%3D%22middle%22%20'%0A%20%20%20%20%20%20%20%20%20%20%20%20f'font-size%3D%2212%22%20fill%3D%22%23475569%22%3E%7Btick%3A.1f%7D%3C%2Ftext%3E'%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20peak_text%20%3D%20%22%2C%20%22.join(%0A%20%20%20%20%20%20%20%20f%22%7Bevent%5B'peak'%5D%3A.3f%7Ds%22%20for%20event%20in%20detected_events%0A%20%20%20%20)%0A%20%20%20%20if%20not%20peak_text%3A%0A%20%20%20%20%20%20%20%20peak_text%20%3D%20%22none%22%0A%0A%20%20%20%20svg%20%3D%20f%22%22%22%0A%20%20%20%20%3Csvg%20viewBox%3D%220%200%20980%20585%22%20width%3D%22100%25%22%20role%3D%22img%22%0A%20%20%20%20%20aria-label%3D%22Interactive%20simulated%20sharp%20wave%20ripple%20detector%20plot%22%0A%20%20%20%20%20style%3D%22font-family%3A%20ui-sans-serif%2C%20system-ui%2C%20sans-serif%3B%22%3E%0A%20%20%20%20%20%20%3Crect%20x%3D%220%22%20y%3D%220%22%20width%3D%22980%22%20height%3D%22585%22%20rx%3D%2218%22%20fill%3D%22%23ffffff%22%20%2F%3E%0A%20%20%20%20%20%20%3Ctext%20x%3D%2270%22%20y%3D%2222%22%20font-size%3D%2214%22%20font-weight%3D%22700%22%20fill%3D%22%230f172a%22%3E%0A%20%20%20%20Simulated%20joint%20SWR%20detection%0A%20%20%20%20%20%20%3C%2Ftext%3E%0A%20%20%20%20%20%20%7B%22%22.join(backgrounds)%7D%0A%20%20%20%20%20%20%7B%22%22.join(overlays)%7D%0A%20%20%20%20%20%20%7B%22%22.join(series)%7D%0A%20%20%20%20%20%20%7B%22%22.join(ticks)%7D%0A%20%20%20%20%20%20%3Ctext%20x%3D%22500%22%20y%3D%22572%22%20text-anchor%3D%22middle%22%20font-size%3D%2213%22%20fill%3D%22%23475569%22%3E%0A%20%20%20%20Time%20(s)%0A%20%20%20%20%20%20%3C%2Ftext%3E%0A%20%20%20%20%20%20%3Crect%20x%3D%22695%22%20y%3D%2212%22%20width%3D%2214%22%20height%3D%2214%22%20fill%3D%22%2316a34a%22%20opacity%3D%220.25%22%20%2F%3E%0A%20%20%20%20%20%20%3Ctext%20x%3D%22715%22%20y%3D%2224%22%20font-size%3D%2212%22%20fill%3D%22%23475569%22%3ESimulated%20SWR%3C%2Ftext%3E%0A%20%20%20%20%20%20%3Crect%20x%3D%22810%22%20y%3D%2212%22%20width%3D%2214%22%20height%3D%2214%22%20fill%3D%22%23f59e0b%22%20opacity%3D%220.38%22%20%2F%3E%0A%20%20%20%20%20%20%3Ctext%20x%3D%22830%22%20y%3D%2224%22%20font-size%3D%2212%22%20fill%3D%22%23475569%22%3EDetected%20SWR%3C%2Ftext%3E%0A%20%20%20%20%3C%2Fsvg%3E%0A%20%20%20%20%22%22%22%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**Detected%20SWRs%3A**%20%7Blen(detected_events)%7D.%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**Peak%20times%3A**%20%7Bpeak_text%7D.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.Html(svg)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A