gw170817 strategies

import%20marimo%0A%0A__generated_with%20%3D%20%220.24.0%22%0Aapp%20%3D%20marimo.App(width%3D%22full%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Aasync%20def%20_()%3A%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20%23%20WASM%20has%20no%20site-packages%2C%20so%20divtel%20has%20to%20come%20from%20a%20wheel.%0A%20%20%20%20if%20sys.platform%20%3D%3D%20%22emscripten%22%3A%0A%20%20%20%20%20%20%20%20import%20json%0A%0A%20%20%20%20%20%20%20%20import%20micropip%0A%20%20%20%20%20%20%20%20from%20pyodide.http%20import%20pyfetch%0A%0A%20%20%20%20%20%20%20%20%23%20Pyodide%20runs%20from%20a%20worker%20under%20assets%2F%2C%20so%20step%20up%20one%20level%0A%20%20%20%20%20%20%20%20%23%20to%20reach%20the%20site%20root%20where%20the%20wheel%20is%20published.%0A%20%20%20%20%20%20%20%20base%20%3D%20%22..%2Fpypi%2F%22%0A%20%20%20%20%20%20%20%20%23%20The%20wheel's%20filename%20carries%20a%20version%2C%20so%20read%20it%20from%20a%0A%20%20%20%20%20%20%20%20%23%20manifest%20instead%20of%20hardcoding%20it.%0A%20%20%20%20%20%20%20%20response%20%3D%20await%20pyfetch(base%20%2B%20%22manifest.json%22)%0A%20%20%20%20%20%20%20%20if%20response.status%20!%3D%20200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20RuntimeError(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22could%20not%20fetch%20%7Bbase%7Dmanifest.json%20(HTTP%20%7Bresponse.status%7D).%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22The%20divtel%20wheel%20is%20missing%20from%20the%20published%20site.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20manifest%20%3D%20json.loads(await%20response.string())%0A%20%20%20%20%20%20%20%20await%20micropip.install(base%20%2B%20manifest%5B%22divtel%22%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20from%20importlib.resources%20import%20files%0A%0A%20%20%20%20import%20astropy.units%20as%20u%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%0A%20%20%20%20from%20divtel%20import%20strategy%2C%20visualization%0A%20%20%20%20from%20divtel.layout%20import%20load_array%0A%20%20%20%20from%20divtel.pointing%20import%20best_pointing%0A%20%20%20%20from%20divtel.region%20import%20SkyRegion%0A%0A%20%20%20%20%23%20SVG%20scales%20to%20its%20container%3B%20marimo's%20PNG%20path%20stamps%20a%20fixed%20pixel%0A%20%20%20%20%23%20width%20that%20overflows%20a%20frame%20narrower%20than%20the%20figure.%0A%20%20%20%20plt.rcParams%5B%22savefig.format%22%5D%20%3D%20%22svg%22%0A%20%20%20%20plt.rcParams%5B%22font.size%22%5D%20%3D%209%0A%0A%20%20%20%20DATA%20%3D%20files(%22divtel%22)%20%2F%20%22data%22%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20DATA%2C%0A%20%20%20%20%20%20%20%20SkyRegion%2C%0A%20%20%20%20%20%20%20%20best_pointing%2C%0A%20%20%20%20%20%20%20%20load_array%2C%0A%20%20%20%20%20%20%20%20mo%2C%0A%20%20%20%20%20%20%20%20np%2C%0A%20%20%20%20%20%20%20%20plt%2C%0A%20%20%20%20%20%20%20%20strategy%2C%0A%20%20%20%20%20%20%20%20u%2C%0A%20%20%20%20%20%20%20%20visualization%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.Html(%22%22%22%3Cstyle%3E%20img%2C%20svg%20%7B%20max-width%3A%20100%25%3B%20height%3A%20auto%3B%20%7D%20%3C%2Fstyle%3E%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Covering%20GW170817%0A%0A%20%20%20%20A%20gravitational-wave%20localization%20is%20a%20probability%20distribution%2C%20not%20a%0A%20%20%20%20point%2C%20and%20the%20early%20ones%20are%20**arcs**%3A%20long%2C%20thin%2C%20and%20lumpy.%20There%20are%0A%20%20%20%20four%20ways%20to%20put%20an%20array%20on%20one%2C%20and%20they%20are%20not%20equally%20good.%0A%0A%20%20%20%20Pick%20a%20map%2C%20a%20site%20and%20a%20strategy%2C%20and%20watch%20the%20array's%20*multiplicity*%20%E2%80%94%0A%20%20%20%20how%20many%20cameras%20have%20each%20direction%20in%20frame%20%E2%80%94%20change%20shape.%20Two%20or%20more%0A%20%20%20%20telescopes%20can%20reconstruct%20a%20shower%3B%20one%20cannot%2C%20so%20grey%20is%20blind%20and%20the%0A%20%20%20%20palest%20blue%20is%20nearly%20so.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20MAPS%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22BAYESTAR%20H%2BL%20%E2%80%94%20the%20first%20alert%2C%20%2B0.7%20h%22%3A%20%22bayestar_hl%22%2C%0A%20%20%20%20%20%20%20%20%22BAYESTAR%20H%2BL%2BV%20%E2%80%94%20with%20Virgo%2C%20%2B5.2%20h%22%3A%20%22bayestar_hlv%22%2C%0A%20%20%20%20%20%20%20%20%22LALInference%20(preliminary)%2C%20%2B16.3%20h%22%3A%20%22lalinference_prelim%22%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20which_map%20%3D%20mo.ui.dropdown(MAPS%2C%20value%3Dlist(MAPS)%5B0%5D%2C%20label%3D%22localization%22)%0A%20%20%20%20level%20%3D%20mo.ui.dropdown(%7B%2250%20%25%22%3A%2050%2C%20%2290%20%25%22%3A%2090%2C%20%2295%20%25%22%3A%2095%7D%2C%20value%3D%2290%20%25%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22credible%20level%22)%0A%20%20%20%20site%20%3D%20mo.ui.dropdown(%7B%22CTAO-South%20(14%20MST%20%2B%2037%20SST)%22%3A%20%22south%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22CTAO-North%20(4%20LST%20%2B%209%20MST)%22%3A%20%22north%22%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20value%3D%22CTAO-South%20(14%20MST%20%2B%2037%20SST)%22%2C%20label%3D%22array%22)%0A%20%20%20%20altitude%20%3D%20mo.ui.slider(20%2C%2085%2C%20value%3D60%2C%20step%3D5%2C%20label%3D%22placed%20at%20altitude%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20show_value%3DTrue)%0A%0A%20%20%20%20mo.hstack(%5Bwhich_map%2C%20level%2C%20site%2C%20altitude%5D%2C%20justify%3D%22start%22%2C%20gap%3D1.5)%0A%20%20%20%20return%20altitude%2C%20level%2C%20site%2C%20which_map%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(DATA%2C%20SkyRegion%2C%20altitude%2C%20level%2C%20load_array%2C%20site%2C%20u%2C%20which_map)%3A%0A%20%20%20%20LAYOUTS%20%3D%20%7B%22north%22%3A%20%22cta-north-lapalma-alpha-prod6.ecsv%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22south%22%3A%20%22cta-south-paranal-alpha-prod6.ecsv%22%7D%0A%0A%20%20%20%20array%20%3D%20load_array(DATA%20%2F%20LAYOUTS%5Bsite.value%5D)%0A%0A%20%20%20%20%23%20The%20region%20is%20rotated%20rigidly%20to%20the%20chosen%20altitude%2C%20then%20cut%20at%20the%0A%20%20%20%20%23%20horizon%3A%20an%20arc%20a%20hundred%20degrees%20long%20can%20run%20underground%2C%20and%20an%0A%20%20%20%20%23%20optimiser%20handed%20that%20will%20aim%20telescopes%20there%20and%20report%20the%20coverage.%0A%20%20%20%20_whole%20%3D%20SkyRegion.from_table(%0A%20%20%20%20%20%20%20%20DATA%20%2F%20%22gw170817%22%20%2F%20f%22%7Bwhich_map.value%7D_%7Blevel.value%7D.ecsv.gz%22)%0A%20%20%20%20region%2C%20underground%20%3D%20_whole.place(altitude.value%20*%20u.deg%2C%20180%20*%20u.deg).visible_part()%0A%20%20%20%20return%20array%2C%20region%2C%20underground%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20STRATEGIES%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22parallel%20%E2%80%94%20everything%20on%20one%20spot%22%3A%20%22parallel%22%2C%0A%20%20%20%20%20%20%20%20%22divergent%20%E2%80%94%20one%20number%2C%20fanned%20out%22%3A%20%22divergent%22%2C%0A%20%20%20%20%20%20%20%20%22tiling%20%E2%80%94%20one%20array%2C%20moved%20in%20steps%22%3A%20%22tiling%22%2C%0A%20%20%20%20%20%20%20%20%22sub-arrays%20%E2%80%94%20split%20it%2C%20point%20the%20pieces%22%3A%20%22subarrays%22%2C%0A%20%20%20%20%20%20%20%20%22sub-arrays%2C%20sized%20by%20probability%22%3A%20%22weighted%22%2C%0A%20%20%20%20%20%20%20%20%22shaped%20%E2%80%94%20every%20telescope%20on%20its%20own%22%3A%20%22shaped%22%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20how%20%3D%20mo.ui.dropdown(STRATEGIES%2C%20value%3Dlist(STRATEGIES)%5B3%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22pointing%20strategy%22)%0A%20%20%20%20div%20%3D%20mo.ui.slider(0.0%2C%200.30%2C%20value%3D0.03%2C%20step%3D0.005%2C%20label%3D%22divergence%20d%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20show_value%3DTrue)%0A%20%20%20%20groups%20%3D%20mo.ui.slider(2%2C%2017%2C%20value%3D8%2C%20step%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22tiles%2C%20or%20sub-arrays%22%2C%20show_value%3DTrue)%0A%20%20%20%20beta%20%3D%20mo.ui.slider(0.0%2C%202.0%2C%20value%3D1.0%2C%20step%3D0.25%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22%CE%B2%20%E2%80%94%20how%20hard%20to%20track%20the%20probability%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20show_value%3DTrue)%0A%0A%20%20%20%20mo.hstack(%5Bhow%2C%20div%2C%20groups%2C%20beta%5D%2C%20justify%3D%22start%22%2C%20gap%3D1.5)%0A%20%20%20%20return%20beta%2C%20div%2C%20groups%2C%20how%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(how%2C%20mo)%3A%0A%20%20%20%20%23%20Shaped%20pointing%20is%20a%20coordinate%20ascent%20over%20every%20telescope%20at%20once.%20In%20a%0A%20%20%20%20%23%20browser%20that%20is%20seconds%20rather%20than%20milliseconds%2C%20so%20it%20does%20not%20run%20on%20a%0A%20%20%20%20%23%20slider%20drag%20--%20and%20the%20settings%20below%20are%20cut%20down%20from%20the%20ones%20the%0A%20%20%20%20%23%20study%20page's%20figures%20were%20computed%20with.%0A%20%20%20%20solve%20%3D%20mo.ui.run_button(label%3D%22solve%20the%20shaped%20pointing%22)%0A%20%20%20%20solve%20if%20how.value%20%3D%3D%20%22shaped%22%20else%20mo.md(%22%22)%0A%20%20%20%20return%20(solve%2C)%0A%0A%0A%40app.cell%0Adef%20_(array%2C%20best_pointing%2C%20beta%2C%20div%2C%20groups%2C%20how%2C%20region%2C%20solve%2C%20strategy)%3A%0A%20%20%20%20exposures%2C%20note%20%3D%201%2C%20%22%22%0A%0A%20%20%20%20if%20how.value%20%3D%3D%20%22parallel%22%3A%0A%20%20%20%20%20%20%20%20aim%20%3D%20best_pointing(array%2C%20region%2C%200.0)%0A%20%20%20%20%20%20%20%20note%20%3D%20f%22aimed%20at%20alt%20%7Baim%5B'alt'%5D%3A.1f%7D%2C%20az%20%7Baim%5B'az'%5D%3A.1f%7D%22%0A%0A%20%20%20%20elif%20how.value%20%3D%3D%20%22divergent%22%3A%0A%20%20%20%20%20%20%20%20best_pointing(array%2C%20region%2C%20div.value)%0A%0A%20%20%20%20elif%20how.value%20%3D%3D%20%22tiling%22%3A%0A%20%20%20%20%20%20%20%20tiles%20%3D%20strategy.tile_region(array%2C%20region%2C%20pointings%3Dgroups.value)%0A%20%20%20%20%20%20%20%20%23%20The%20picture%20shows%20one%20tile%3B%20the%20coverage%20below%20is%20what%20the%20whole%0A%20%20%20%20%20%20%20%20%23%20sequence%20reaches%2C%20which%20costs%20one%20exposure%20per%20tile.%0A%20%20%20%20%20%20%20%20array.divergent_pointing(0.0%2C%20tiles%5B0%5D%5B%22alt%22%5D%2C%20tiles%5B0%5D%5B%22az%22%5D)%0A%20%20%20%20%20%20%20%20exposures%20%3D%20len(tiles)%0A%20%20%20%20%20%20%20%20note%20%3D%20(f%22%7B100%20*%20tiles%5B-1%5D%5B'cumulative'%5D%3A.0f%7D%20%25%20after%20all%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Blen(tiles)%7D%20tiles%3B%20the%20panel%20shows%20the%20first%22)%0A%0A%20%20%20%20elif%20how.value%20%3D%3D%20%22subarrays%22%3A%0A%20%20%20%20%20%20%20%20summary%20%3D%20strategy.point_subarrays(array%2C%20region%2C%20groups.value)%0A%20%20%20%20%20%20%20%20note%20%3D%20f%22group%20sizes%20%7Bsummary%5B'sizes'%5D%7D%22%0A%0A%20%20%20%20elif%20how.value%20%3D%3D%20%22weighted%22%3A%0A%20%20%20%20%20%20%20%20summary%20%3D%20strategy.weighted_split(array%2C%20region%2C%20groups.value)%0A%20%20%20%20%20%20%20%20note%20%3D%20f%22group%20sizes%20%7Bsummary%5B'sizes'%5D%7D%22%0A%0A%20%20%20%20elif%20how.value%20%3D%3D%20%22shaped%22%3A%0A%20%20%20%20%20%20%20%20if%20solve.value%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Cut%20down%20for%20the%20browser%3A%20fewer%20candidate%20pointings%2C%20fewer%20scored%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20directions%2C%20fewer%20sweeps%20and%20one%20warm%20start%20instead%20of%20seven.%0A%20%20%20%20%20%20%20%20%20%20%20%20strategy.shaped_pointing(array%2C%20region%2C%20beta%3Dbeta.value%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20candidates%3D400%2C%20search_pixels%3D600%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20sweeps%3D4%2C%20rounds%3D1%2C%20warm_starts%3D(6%2C))%0A%20%20%20%20%20%20%20%20%20%20%20%20note%20%3D%20%22cut-down%20solve%3B%20the%20study%20page%20uses%20the%20full%20budget%22%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20best_pointing(array%2C%20region%2C%200.0)%0A%20%20%20%20%20%20%20%20%20%20%20%20note%20%3D%20%22press%20the%20button%20to%20solve%22%0A%0A%20%20%20%20scored%20%3D%20strategy.describe(array%2C%20region%2C%20m_cut%3D2)%0A%20%20%20%20profile%20%3D%20strategy.multiplicity_by_probability(array%2C%20region)%0A%0A%20%20%20%20%23%20Every%20strategy%20above%20re-points%20%60array%60%20in%20place%2C%20and%20marimo%20re-runs%20a%20cell%0A%20%20%20%20%23%20when%20a%20variable%20it%20reads%20is%20*redefined*%20--%20mutation%20is%20invisible%20to%20it.%20So%0A%20%20%20%20%23%20hand%20the%20pointed%20array%20on%20under%20a%20new%20name%2C%20or%20the%20figures%20below%20go%20stale%0A%20%20%20%20%23%20while%20the%20numbers%20beside%20them%20update.%0A%20%20%20%20pointed%20%3D%20array%0A%20%20%20%20return%20exposures%2C%20note%2C%20pointed%2C%20profile%2C%20scored%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(exposures%2C%20mo%2C%20np%2C%20pointed%2C%20region%2C%20scored%2C%20u%2C%20underground)%3A%0A%20%20%20%20%23%20How%20many%20telescopes%20have%20any%20part%20of%20the%20region%20in%20frame.%20It%20is%20the%0A%20%20%20%20%23%20number%20that%20separates%20the%20strategies%20most%20sharply%3A%20divergence%20spreads%0A%20%20%20%20%23%20the%20array%20so%20wide%20that%20a%20third%20of%20it%20ends%20up%20staring%20at%20empty%20sky.%0A%20%20%20%20_cos_radii%20%3D%20np.cos(%5Bt.fov_radius.to_value(u.rad)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20t%20in%20pointed.telescopes%5D)%0A%20%20%20%20_on_target%20%3D%20int(((region.directions%20%40%20pointed.pointing_vectors.T)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%3E%3D%20_cos_radii).any(axis%3D0).sum())%0A%0A%20%20%20%20mo.hstack(%5B%0A%20%20%20%20%20%20%20%20mo.stat(label%3D%22covered%20in%20stereo%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20value%3Df%22%7B100%20*%20scored%5B'covered_stereo'%5D%3A.0f%7D%20%25%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20caption%3Df%22of%20the%20region%20above%20the%20horizon%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22(%7B100%20*%20underground%3A.0f%7D%20%25%20of%20the%20map%20is%20not)%22)%2C%0A%20%20%20%20%20%20%20%20mo.stat(label%3D%22telescopes%20per%20shower%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20value%3Df%22%7Bscored%5B'mean_covered'%5D%3A.1f%7D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20caption%3D%22averaged%20over%20the%20part%20covered%22)%2C%0A%20%20%20%20%20%20%20%20mo.stat(label%3D%22telescopes%20on%20target%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20value%3Df%22%7B_on_target%7D%20of%20%7Blen(pointed.telescopes)%7D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20caption%3D%22the%20rest%20are%20pointed%20at%20empty%20sky%22)%2C%0A%20%20%20%20%20%20%20%20mo.stat(label%3D%22exposures%22%2C%20value%3Dstr(exposures)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20caption%3D%22the%20window%20is%20divided%20between%20them%22)%2C%0A%20%20%20%20%5D%2C%20justify%3D%22start%22%2C%20gap%3D2)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20note)%3A%0A%20%20%20%20mo.md(f%22*%7Bnote%7D*%22)%20if%20note%20else%20mo.md(%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(plt%2C%20pointed%2C%20region%2C%20visualization)%3A%0A%20%20%20%20%23%20Where%20the%20cameras%20actually%20are%2C%20over%20the%20probability%20they%20are%20covering.%0A%20%20%20%20%23%20The%20panel%20below%20shades%20the%20result%3B%20this%20one%20shows%20the%20cause.%0A%20%20%20%20_frame%20%3D%20visualization.projection_frame(region)%0A%20%20%20%20_figure%2C%20_ax%20%3D%20plt.subplots(figsize%3D(11.0%2C%204.0))%0A%0A%20%20%20%20visualization.probability_over_region(region%2C%20ax%3D_ax%2C%20frame%3D_frame%2C%20size%3D3.0)%0A%20%20%20%20_rims%20%3D%20visualization.camera_rims(pointed%2C%20ax%3D_ax%2C%20frame%3D_frame%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20colors%3D%22type%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20type_names%3D(%22SST%22%2C%20%22MST%22))%0A%20%20%20%20visualization.frame_on(%0A%20%20%20%20%20%20%20%20_ax%2C%20_rims%2C%20visualization.region_extent(region%2C%20frame%3D_frame))%0A%0A%20%20%20%20_legend%20%3D%20_ax.legend(frameon%3DFalse%2C%20fontsize%3D8%2C%20loc%3D%22upper%20left%22%2C%20ncol%3D2)%0A%20%20%20%20for%20_handle%20in%20_legend.legend_handles%3A%0A%20%20%20%20%20%20%20%20if%20hasattr(_handle%2C%20%22set_sizes%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_handle.set_sizes(%5B26%5D)%0A%20%20%20%20_ax.set_title(%22the%20cameras%2C%20over%20the%20probability%20they%20cover%22%2C%20loc%3D%22left%22)%0A%20%20%20%20_figure.tight_layout()%0A%20%20%20%20_figure%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(plt%2C%20pointed%2C%20region%2C%20visualization)%3A%0A%20%20%20%20%23%20Its%20own%20figure%2C%20so%20the%20equal-aspect%20panel%20can%20take%20the%20full%20width%20rather%0A%20%20%20%20%23%20than%20sharing%20it%20with%20a%20bar%20chart%20that%20wants%20a%20different%20shape.%0A%20%20%20%20_figure%2C%20_ax%20%3D%20plt.subplots(%0A%20%20%20%20%20%20%20%20figsize%3D(11.0%2C%201.6%20%2B%2011.0%20*%20visualization.region_span(region)))%0A%20%20%20%20visualization.multiplicity_over_region(pointed%2C%20region%2C%20ax%3D_ax%2C%20size%3D6.0)%0A%20%20%20%20_ax.set_title(%22what%20that%20delivers%3A%20telescopes%20per%20direction%22%2C%20loc%3D%22left%22)%0A%20%20%20%20_figure.tight_layout()%0A%20%20%20%20_figure%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(plt%2C%20profile%2C%20visualization)%3A%0A%20%20%20%20_figure%2C%20_ax%20%3D%20plt.subplots(figsize%3D(11.0%2C%203.0))%0A%20%20%20%20visualization.multiplicity_by_probability(%7B%22this%20configuration%22%3A%20profile%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ax%3D_ax)%0A%20%20%20%20_figure.tight_layout()%0A%20%20%20%20_figure%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20profile%2C%20scored)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20**Reading%20the%20three%20panels.**%20The%20first%20draws%20every%20telescope's%20field%20of%0A%20%20%20%20view%20over%20the%20localization%2C%20shaded%20darkest%20where%20the%20source%20most%20likely%20is.%0A%20%20%20%20Rims%20are%20coloured%20by%20camera%20size%2C%20so%20you%20can%20watch%20free%20pointing%20send%20the%0A%20%20%20%20narrow%20cameras%20where%20the%20map%20is%20sharp%20%E2%80%94%20something%20neither%20divergence%20nor%20an%0A%20%20%20%20even%20split%20can%20do.%20All%20three%20panels%20are%20equal-area%2C%20so%20a%20wide%20patch%20of%0A%20%20%20%20shallow%20coverage%20looks%20as%20big%20as%20it%20is.%0A%0A%20%20%20%20The%20second%20shades%20the%20map%20by%20how%20many%20cameras%20land%20on%20each%20direction%3B%20grey%0A%20%20%20%20is%20sky%20no%20camera%20reaches.%20The%20third%20cuts%20the%20map%20into%20fifths%20of%0A%20%20%20%20*probability*%20%E2%80%94%20the%20last%20fifth%20is%20the%20small%20bright%20core%2C%20the%20first%20is%20the%0A%20%20%20%20wide%20faint%20skirt%20%E2%80%94%20and%20asks%20how%20deeply%20each%20is%20watched.%0A%0A%20%20%20%20A%20flat%20profile%20is%20an%20array%20spread%20evenly%20over%20the%20region%2C%20which%20is%20only%20the%0A%20%20%20%20right%20answer%20if%20the%20map%20is%20flat.%20This%20one%20runs%0A%20%20%20%20**%7Bprofile%5B0%5D%3A.1f%7D%20%E2%86%92%20%7Bprofile%5B-1%5D%3A.1f%7D**%2C%20and%20its%20mismatch%20against%20the%0A%20%20%20%20probability%20is%20**%7Bscored%5B'mismatch'%5D%3A.3f%7D**%20(zero%20would%20be%20depth%20exactly%0A%20%20%20%20proportional%20to%20probability%20over%20the%20part%20that%20clears%20the%20stereo%20cut).%0A%0A%20%20%20%20Read%20that%20number%20beside%20the%20coverage%20and%20never%20alone%3A%20an%20array%20that%0A%20%20%20%20abandons%20everything%20but%20the%20core%20matches%20the%20shape%20of%20what%20it%20kept%2C%0A%20%20%20%20perfectly%2C%20and%20is%20useless.%0A%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20Everything%20here%20is%20geometry%3A%20where%20cameras%20point%20and%20what%20they%20cover.%20There%0A%20%20%20%20is%20no%20effective%20area%2C%20no%20energy%20threshold%20and%20no%20background%2C%20so%20none%20of%0A%20%20%20%20these%20percentages%20is%20a%20sensitivity.%20Coverage%20is%20necessary%20for%20a%20detection%0A%20%20%20%20and%20nowhere%20near%20sufficient.%0A%0A%20%20%20%20The%20full%20argument%2C%20with%20the%20numbers%2C%20is%20in%0A%20%20%20%20%5Bthe%20study%20page%5D(..%2F..%2Fgw170817.html).%0A%20%20%20%20%22%22%22)%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