GSF — Model Comparison

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.15%22%0Aapp%20%3D%20marimo.App(%0A%20%20%20%20width%3D%22full%22%2C%0A%20%20%20%20app_title%3D%22GSF%20%E2%80%94%20Model%20Comparison%22%2C%0A%20%20%20%20css_file%3D%22gallery.css%22%2C%0A%20%20%20%20html_head_file%3D%22gallery_head.html%22%2C%0A)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20(mo%2C)%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%20Model%20Comparison%3A%20GSF%202017%20vs%202025%0A%0A%20%20%20%20The%202025%20update%20extends%20the%20input%20dataset%20and%20refits%20the%20splines%2C%20changing%20both%20central%20fluxes%20and%20uncertainties.%20We%20compare%20the%20two%20versions%20for%20%60GSFEnergy%60%20and%20%60GSFEnergyPerNucleon%60%2C%20scaled%20by%20%24E%5E%7B2.7%7D%24.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Aasync%20def%20_()%3A%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20NOTEBOOK%20%3D%20(%0A%20%20%20%20%20%20%20%20%22model_comparison_2017_vs_2025%22%20%20%23%20used%20by%20the%20download%20buttons%20at%20the%20end%0A%20%20%20%20)%0A%20%20%20%20IN_WASM%20%3D%20%22pyodide%22%20in%20sys.modules%0A%20%20%20%20SITE%20%3D%20%22%22%0A%0A%20%20%20%20if%20IN_WASM%3A%0A%20%20%20%20%20%20%20%20%23%20WASM%20(docs%20gallery)%3A%20TEMPORARY%20until%20globalsplinefit%20is%20on%20PyPI%20%E2%80%94%0A%20%20%20%20%20%20%20%20%23%20install%20the%20wheel%20published%20under%20%2Fwheels%2F%20on%20the%20same%20site%20the%0A%20%20%20%20%20%20%20%20%23%20notebook%20is%20served%20from%20(also%20works%20on%20staging%20or%20local%20hosts).%0A%20%20%20%20%20%20%20%20from%20urllib.parse%20import%20urljoin%0A%0A%20%20%20%20%20%20%20%20import%20marimo%20as%20_mo%0A%20%20%20%20%20%20%20%20import%20micropip%0A%0A%20%20%20%20%20%20%20%20%23%20notebook_location()%20resolution%20depth%20differs%20between%20marimo%20export%0A%20%20%20%20%20%20%20%20%23%20layouts%20%E2%80%94%20try%20site-root%20%2Fwheels%2F%20from%20every%20plausible%20depth%2C%20and%0A%20%20%20%20%20%20%20%20%23%20keep%20the%20depth%20that%20worked%20as%20the%20site%20root%20for%20other%20assets.%0A%20%20%20%20%20%20%20%20_base%20%3D%20str(_mo.notebook_location())%20%2B%20%22%2F%22%0A%20%20%20%20%20%20%20%20_err%20%3D%20None%0A%20%20%20%20%20%20%20%20for%20_up%20in%20(%22%22%2C%20%22..%2F%22%2C%20%22..%2F..%2F%22%2C%20%22..%2F..%2F..%2F%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20await%20micropip.install(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20urljoin(%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_base%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_up%20%2B%20%22wheels%2Fglobalsplinefit-2.0.0a1-py3-none-any.whl%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20SITE%20%3D%20urljoin(_base%2C%20_up)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20None%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%0A%20%20%20%20%20%20%20%20%20%20%20%20except%20Exception%20as%20_e%3A%20%20%23%20noqa%3A%20BLE001%20%E2%80%94%20404%20lands%20as%20generic%20error%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20_e%0A%20%20%20%20%20%20%20%20if%20_err%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20_err%0A%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%0A%20%20%20%20from%20globalsplinefit%20import%20MODEL_VERSIONS%2C%20GSFEnergy%2C%20GSFEnergyPerNucleon%0A%0A%20%20%20%20plt.rcParams%5B%22figure.figsize%22%5D%20%3D%20(8%2C%205)%0A%20%20%20%20plt.rcParams%5B%22font.size%22%5D%20%3D%2011%0A%20%20%20%20plt.rcParams%5B%22lines.linewidth%22%5D%20%3D%201.2%0A%20%20%20%20plt.rcParams%5B%22grid.alpha%22%5D%20%3D%200.3%0A%0A%20%20%20%20def%20show(fig%3DNone)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Stamp%20and%20return%20a%20figure%20as%20the%20cell%20output.%0A%0A%20%20%20%20%20%20%20%20marimo%20renders%20a%20cell's%20last%20expression%3B%20%60plt.show()%60%20produces%20no%0A%20%20%20%20%20%20%20%20output%20at%20all%20in%20app%20view%2C%20so%20every%20plot%20cell%20ends%20with%20%60show(fig)%60.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20fig%20%3D%20plt.gcf()%20if%20fig%20is%20None%20else%20fig%0A%20%20%20%20%20%20%20%20%23%20TEMPORARY%20pre-publication%20marker%20(_gsf_preliminary)%3A%20diagonal%0A%20%20%20%20%20%20%20%20%23%20PRELIMINARY%20watermark%20on%20every%20figure.%20Delete%20this%20block%20at%20the%0A%20%20%20%20%20%20%20%20%23%20GSF%202026%20release%20%E2%80%94%20keep%20the%20%60return%20fig%60.%0A%20%20%20%20%20%20%20%20fig.text(%0A%20%20%20%20%20%20%20%20%20%20%20%200.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%200.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22PRELIMINARY%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22center%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20va%3D%22center%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20rotation%3D30%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D42%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontweight%3D%22bold%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22gray%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.18%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D1000%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20return%20fig%0A%0A%20%20%20%20def%20autoscale(ax%2C%20*series%2C%20pad%3D1.6%2C%20decades%3DNone)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Fit%20the%20y-axis%20to%20the%20positive%2C%20finite%20values%20in%20%60series%60.%0A%0A%20%20%20%20%20%20%20%20With%20%60decades%60%2C%20the%20lower%20limit%20is%20that%20many%20decades%20below%20the%20top%0A%20%20%20%20%20%20%20%20instead%20of%20the%20smallest%20value%3A%20every%20group%20here%20ends%20in%20a%20cutoff%0A%20%20%20%20%20%20%20%20tail%20that%20falls%20another%20eight%20decades%2C%20and%20letting%20those%20set%20the%0A%20%20%20%20%20%20%20%20bottom%20of%20the%20axis%20flattens%20the%20entire%20comparison%20into%20a%20strip.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20_v%20%3D%20np.concatenate(%5Bnp.asarray(s%2C%20float).ravel()%20for%20s%20in%20series%5D)%0A%20%20%20%20%20%20%20%20_v%20%3D%20_v%5Bnp.isfinite(_v)%20%26%20(_v%20%3E%200)%5D%0A%20%20%20%20%20%20%20%20if%20_v.size%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_top%20%3D%20_v.max()%20*%20pad%0A%20%20%20%20%20%20%20%20%20%20%20%20ax.set_ylim(_top%20%2F%2010**decades%20if%20decades%20else%20_v.min()%20%2F%20pad%2C%20_top)%0A%0A%20%20%20%20def%20relative(flux%2C%20error%2C%20cap%3D3.0)%3A%0A%20%20%20%20%20%20%20%20%22%22%22sigma%2Fflux%20up%20to%20%60cap%60%2C%20blank%20above%20it.%0A%0A%20%20%20%20%20%20%20%20Past%20the%20end%20of%20a%20group's%20fitted%20range%20the%20flux%20collapses%20(and%0A%20%20%20%20%20%20%20%20eventually%20underflows%20to%20exactly%20zero)%20while%20its%20error%20stays%20finite%2C%0A%20%20%20%20%20%20%20%20so%20the%20ratio%20runs%20away.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20with%20np.errstate(divide%3D%22ignore%22%2C%20invalid%3D%22ignore%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_r%20%3D%20np.where(flux%20%3E%200%2C%20error%20%2F%20flux%2C%20np.inf)%0A%20%20%20%20%20%20%20%20return%20np.where(np.isfinite(_r)%20%26%20(_r%20%3C%3D%20cap)%2C%20_r%2C%20np.nan)%0A%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20GSFEnergy%2C%0A%20%20%20%20%20%20%20%20GSFEnergyPerNucleon%2C%0A%20%20%20%20%20%20%20%20NOTEBOOK%2C%0A%20%20%20%20%20%20%20%20SITE%2C%0A%20%20%20%20%20%20%20%20autoscale%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%20relative%2C%0A%20%20%20%20%20%20%20%20show%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20exponent%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20start%3D0.0%2C%0A%20%20%20%20%20%20%20%20stop%3D3.0%2C%0A%20%20%20%20%20%20%20%20step%3D0.1%2C%0A%20%20%20%20%20%20%20%20value%3D2.7%2C%0A%20%20%20%20%20%20%20%20label%3D%22flux%20scaling%20%20E%5Ek%22%2C%0A%20%20%20%20%20%20%20%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20plot_size%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20start%3D0.6%2C%0A%20%20%20%20%20%20%20%20stop%3D1.6%2C%0A%20%20%20%20%20%20%20%20step%3D0.1%2C%0A%20%20%20%20%20%20%20%20value%3D1.0%2C%0A%20%20%20%20%20%20%20%20label%3D%22plot%20size%22%2C%0A%20%20%20%20%20%20%20%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20mo.hstack(%5Bexponent%2C%20plot_size%5D%2C%20justify%3D%22start%22%2C%20gap%3D2)%0A%20%20%20%20return%20exponent%2C%20plot_size%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20plot_size)%3A%0A%20%20%20%20%23%20Drives%20--gsf-plot-scale%2C%20which%20gallery.css%20multiplies%20into%20every%20figure's%0A%20%20%20%20%23%20width.%20Scaling%20the%20rendered%20raster%20(rather%20than%20the%20figsize)%20is%20what%20keeps%0A%20%20%20%20%23%20it%20proportional%3A%20fonts%2C%20line%20widths%20and%20the%20aspect%20ratio%20all%20ride%20along.%0A%20%20%20%20mo.Html(f%22%3Cstyle%3E%3Aroot%7B%7B--gsf-plot-scale%3A%7Bplot_size.value%7D%7D%7D%3C%2Fstyle%3E%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(exponent)%3A%0A%20%20%20%20energy_exponent%20%3D%20exponent.value%20%20%23%20display%20scaling%2C%20driven%20by%20the%20slider%0A%20%20%20%20return%20(energy_exponent%2C)%0A%0A%0A%40app.cell%0Adef%20_(GSFEnergy%2C%20GSFEnergyPerNucleon%2C%20np)%3A%0A%20%20%20%20%23%20200%20points%20is%20plenty%20on%20a%20log%20grid%2C%20and%20the%20browser%20runtime%20pays%20for%0A%20%20%20%20%23%20every%20one%20of%20them%20twice%20over%20(two%20versions%20x%20flux%20%2B%20covariance%20error)%0A%20%20%20%20energy%20%3D%20np.logspace(np.log10(0.5)%2C%20np.log10(100000000000.0)%2C%20200)%0A%20%20%20%20energy_per_nucleon%20%3D%20energy%0A%20%20%20%20compared_versions%20%3D%20%5B%222017%22%2C%20%222025%22%5D%0A%20%20%20%20models_e%20%3D%20%7B_v%3A%20GSFEnergy(version%3D_v)%20for%20_v%20in%20compared_versions%7D%0A%20%20%20%20models_n%20%3D%20%7B_v%3A%20GSFEnergyPerNucleon(version%3D_v)%20for%20_v%20in%20compared_versions%7D%0A%20%20%20%20groups%20%3D%20%5B%22p%22%2C%20%22He%22%2C%20%22O*%22%2C%20%22Fe*%22%5D%0A%20%20%20%20group_labels%20%3D%20%5B%22H*%22%2C%20%22He*%22%2C%20%22O*%22%2C%20%22Fe*%22%5D%0A%20%20%20%20group_colors%20%3D%20%5B%22red%22%2C%20%22orange%22%2C%20%22green%22%2C%20%22blue%22%5D%0A%20%20%20%20line_styles%20%3D%20%7B%222017%22%3A%20%22--%22%2C%20%222025%22%3A%20%22-%22%7D%0A%0A%20%20%20%20def%20collect(models%2C%20grid)%3A%0A%20%20%20%20%20%20%20%20_flux%20%3D%20%7B_v%3A%20%7Bg%3A%20m.flux(grid%2C%20g)%20for%20g%20in%20groups%7D%20for%20_v%2C%20m%20in%20models.items()%7D%0A%20%20%20%20%20%20%20%20err%20%3D%20%7B_v%3A%20%7Bg%3A%20m.error(grid%2C%20g)%20for%20g%20in%20groups%7D%20for%20_v%2C%20m%20in%20models.items()%7D%0A%20%20%20%20%20%20%20%20for%20_v%2C%20m%20in%20models.items()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%5B_v%5D%5B%22Total%22%5D%20%3D%20m.total_flux(grid)%0A%20%20%20%20%20%20%20%20%20%20%20%20err%5B_v%5D%5B%22Total%22%5D%20%3D%20m.total_error(grid)%0A%20%20%20%20%20%20%20%20return%20(_flux%2C%20err)%0A%0A%20%20%20%20energy_fluxes%2C%20energy_errors%20%3D%20collect(models_e%2C%20energy)%0A%20%20%20%20nucleon_fluxes%2C%20nucleon_errors%20%3D%20collect(models_n%2C%20energy_per_nucleon)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20energy_errors%2C%0A%20%20%20%20%20%20%20%20energy_fluxes%2C%0A%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20group_colors%2C%0A%20%20%20%20%20%20%20%20group_labels%2C%0A%20%20%20%20%20%20%20%20groups%2C%0A%20%20%20%20%20%20%20%20line_styles%2C%0A%20%20%20%20%20%20%20%20nucleon_errors%2C%0A%20%20%20%20%20%20%20%20nucleon_fluxes%2C%0A%20%20%20%20)%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%23%20Total-energy%20flux%20comparison%0A%0A%20%20%20%20Solid%3A%202025.%20Dashed%3A%202017.%0A%0A%20%20%20%201.%20Per-group%20%2B%20total%20flux.%0A%20%20%20%202.%20Relative%20uncertainty%20%24%5Csigma%2F%5CPhi%24.%0A%20%20%20%203.%20Flux%20with%201-sigma%20error%20bands.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20autoscale%2C%0A%20%20%20%20energy%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_fluxes%2C%0A%20%20%20%20group_colors%2C%0A%20%20%20%20group_labels%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%0A%20%20%20%20for%20_group%2C%20_label%2C%20_color%20in%20zip(groups%2C%20group_labels%2C%20group_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%20%3D%20energy_fluxes%5B_v%5D%5B_group%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20energy**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22%7B_label%7D%20(%7B_v%7D)%22%20if%20_v%20%3D%3D%20%222025%22%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20_flux%20%3D%20energy_fluxes%5B_v%5D%5B%22Total%22%5D%0A%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20energy**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.9%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22Total%20(%7B_v%7D)%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20%5BGeV%5D%22)%0A%20%20%20%20_ax.set_ylabel(f%22Flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20center%22%2C%20ncol%3D2)%0A%20%20%20%20_ax.set_xlim(energy%5B0%5D%2C%20energy%5B-1%5D)%0A%20%20%20%20autoscale(_ax%2C%20*%5B_l.get_ydata()%20for%20_l%20in%20_ax.get_lines()%5D%2C%20decades%3D3)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20energy%2C%0A%20%20%20%20energy_errors%2C%0A%20%20%20%20energy_fluxes%2C%0A%20%20%20%20group_colors%2C%0A%20%20%20%20group_labels%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20relative%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%0A%20%20%20%20for%20_group%2C%20_label%2C%20_color%20in%20zip(groups%2C%20group_labels%2C%20group_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_rel%20%3D%20relative(energy_fluxes%5B_v%5D%5B_group%5D%2C%20energy_errors%5B_v%5D%5B_group%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_rel%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22%7B_label%7D%20(%7B_v%7D)%22%20if%20_v%20%3D%3D%20%222025%22%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20_rel%20%3D%20relative(energy_fluxes%5B_v%5D%5B%22Total%22%5D%2C%20energy_errors%5B_v%5D%5B%22Total%22%5D)%0A%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_rel%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.9%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22Total%20(%7B_v%7D)%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20%5BGeV%5D%22)%0A%20%20%20%20_ax.set_ylabel(%22Relative%20error%20%CF%83%2F%CE%A6%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20right%22%2C%20ncol%3D2)%0A%20%20%20%20_ax.set_xlim(energy%5B0%5D%2C%20energy%5B-1%5D)%0A%20%20%20%20_ax.set_ylim(0.003%2C%203)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20autoscale%2C%0A%20%20%20%20energy%2C%0A%20%20%20%20energy_errors%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_fluxes%2C%0A%20%20%20%20group_colors%2C%0A%20%20%20%20group_labels%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%0A%0A%20%20%20%20def%20plot_with_band(ax%2C%20x%2C%20flux%2C%20err%2C%20color%2C%20ls%2C%20lw%3D1%2C%20label%3DNone)%3A%0A%20%20%20%20%20%20%20%20scaled%20%3D%20flux%20*%20x**energy_exponent%0A%20%20%20%20%20%20%20%20scaled_err%20%3D%20err%20*%20x**energy_exponent%0A%20%20%20%20%20%20%20%20ax.loglog(x%2C%20scaled%2C%20color%3Dcolor%2C%20ls%3Dls%2C%20lw%3Dlw%2C%20alpha%3D0.8%2C%20label%3Dlabel)%0A%20%20%20%20%20%20%20%20ax.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%2C%20scaled%20-%20scaled_err%2C%20scaled%20%2B%20scaled_err%2C%20color%3Dcolor%2C%20alpha%3D0.2%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20for%20_group%2C%20_label%2C%20_color%20in%20zip(groups%2C%20group_labels%2C%20group_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20plot_with_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_fluxes%5B%222017%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_errors%5B%222017%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22--%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20plot_with_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_fluxes%5B%222025%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_errors%5B%222025%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22-%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3D_label%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20for%20_v%2C%20_ls%20in%20line_styles.items()%3A%0A%20%20%20%20%20%20%20%20plot_with_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_fluxes%5B_v%5D%5B%22Total%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_errors%5B_v%5D%5B%22Total%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22black%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_ls%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22Total%20%7B_v%7D%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20%5BGeV%5D%22)%0A%20%20%20%20_ax.set_ylabel(f%22Flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20center%22%2C%20ncol%3D2)%0A%20%20%20%20_ax.set_xlim(energy%5B0%5D%2C%20energy%5B-1%5D)%0A%20%20%20%20autoscale(_ax%2C%20*%5B_l.get_ydata()%20for%20_l%20in%20_ax.get_lines()%5D%2C%20decades%3D3)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%20(plot_with_band%2C)%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%23%20Energy-per-nucleon%20flux%20comparison%0A%0A%20%20%20%20Same%20three%20plots%20in%20the%20%60GSFEnergyPerNucleon%60%20view.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20autoscale%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20group_colors%2C%0A%20%20%20%20group_labels%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20nucleon_fluxes%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%0A%20%20%20%20for%20_group%2C%20_label%2C%20_color%20in%20zip(groups%2C%20group_labels%2C%20group_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%20%3D%20nucleon_fluxes%5B_v%5D%5B_group%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20energy_per_nucleon**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22%7B_label%7D%20(%7B_v%7D)%22%20if%20_v%20%3D%3D%20%222025%22%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20_flux%20%3D%20nucleon_fluxes%5B_v%5D%5B%22Total%22%5D%0A%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20energy_per_nucleon**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.9%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22Total%20(%7B_v%7D)%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20per%20nucleon%20%5BGeV%2Fnucleon%5D%22)%0A%20%20%20%20_ax.set_ylabel(f%22Nucleon%20flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22upper%20right%22%2C%20ncol%3D2)%0A%20%20%20%20_ax.set_xlim(energy_per_nucleon%5B0%5D%2C%20energy_per_nucleon%5B-1%5D)%0A%20%20%20%20autoscale(_ax%2C%20*%5B_l.get_ydata()%20for%20_l%20in%20_ax.get_lines()%5D%2C%20decades%3D3)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20group_colors%2C%0A%20%20%20%20group_labels%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20nucleon_errors%2C%0A%20%20%20%20nucleon_fluxes%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20relative%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%0A%20%20%20%20for%20_group%2C%20_label%2C%20_color%20in%20zip(groups%2C%20group_labels%2C%20group_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_rel%20%3D%20relative(nucleon_fluxes%5B_v%5D%5B_group%5D%2C%20nucleon_errors%5B_v%5D%5B_group%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_rel%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22%7B_label%7D%20(%7B_v%7D)%22%20if%20_v%20%3D%3D%20%222025%22%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20for%20_v%20in%20(%222017%22%2C%20%222025%22)%3A%0A%20%20%20%20%20%20%20%20_rel%20%3D%20relative(nucleon_fluxes%5B_v%5D%5B%22Total%22%5D%2C%20nucleon_errors%5B_v%5D%5B%22Total%22%5D)%0A%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_rel%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dline_styles%5B_v%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22black%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.9%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22Total%20(%7B_v%7D)%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20per%20nucleon%20%5BGeV%2Fnucleon%5D%22)%0A%20%20%20%20_ax.set_ylabel(%22Relative%20error%20%CF%83%2F%CE%A6%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20right%22%2C%20ncol%3D2)%0A%20%20%20%20_ax.set_xlim(energy_per_nucleon%5B0%5D%2C%20energy_per_nucleon%5B-1%5D)%0A%20%20%20%20_ax.set_ylim(0.003%2C%203)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20autoscale%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20group_colors%2C%0A%20%20%20%20group_labels%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20nucleon_errors%2C%0A%20%20%20%20nucleon_fluxes%2C%0A%20%20%20%20plot_with_band%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%0A%20%20%20%20for%20_group%2C%20_label%2C%20_color%20in%20zip(groups%2C%20group_labels%2C%20group_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20plot_with_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nucleon_fluxes%5B%222017%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nucleon_errors%5B%222017%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22--%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20plot_with_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nucleon_fluxes%5B%222025%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nucleon_errors%5B%222025%22%5D%5B_group%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22-%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3D_label%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20for%20_v%2C%20_ls%20in%20line_styles.items()%3A%0A%20%20%20%20%20%20%20%20plot_with_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nucleon_fluxes%5B_v%5D%5B%22Total%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20nucleon_errors%5B_v%5D%5B%22Total%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22black%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_ls%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22Total%20%7B_v%7D%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20per%20nucleon%20%5BGeV%2Fnucleon%5D%22)%0A%20%20%20%20_ax.set_ylabel(f%22Nucleon%20flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20center%22%2C%20ncol%3D2)%0A%20%20%20%20_ax.set_xlim(energy_per_nucleon%5B0%5D%2C%20energy_per_nucleon%5B-1%5D)%0A%20%20%20%20autoscale(_ax%2C%20*%5B_l.get_ydata()%20for%20_l%20in%20_ax.get_lines()%5D%2C%20decades%3D3)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%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%23%20Ratio%20of%202017%20to%202025%0A%0A%20%20%20%20Top%3A%20total%20nucleon%20flux%20ratio%20%24%5CPhi_%7B2017%7D%2F%5CPhi_%7B2025%7D%24%20with%20the%20relative%20uncertainty%20of%20each%20version%20shaded.%20Bottom%3A%20ratio%20of%20absolute%20errors.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(energy_per_nucleon%2C%20np%2C%20nucleon_errors%2C%20nucleon_fluxes%2C%20plt%2C%20show)%3A%0A%20%20%20%20_fig%2C%20(ax1%2C%20ax2)%20%3D%20plt.subplots(2%2C%201%2C%20sharex%3DTrue)%0A%20%20%20%20%23%20NaN%20out%20the%20underflowed%20tail%20(both%20totals%20reach%20zero%20past%20the%20fitted%0A%20%20%20%20%23%20range)%20to%20keep%20the%20ratios%20finite%3B%20the%20axis%20stops%20well%20below%20it%20anyway%0A%20%20%20%20_pos%20%3D%20(nucleon_fluxes%5B%222017%22%5D%5B%22Total%22%5D%20%3E%200)%20%26%20(nucleon_fluxes%5B%222025%22%5D%5B%22Total%22%5D%20%3E%200)%0A%20%20%20%20flux_2017%20%3D%20np.where(_pos%2C%20nucleon_fluxes%5B%222017%22%5D%5B%22Total%22%5D%2C%20np.nan)%0A%20%20%20%20flux_2025%20%3D%20np.where(_pos%2C%20nucleon_fluxes%5B%222025%22%5D%5B%22Total%22%5D%2C%20np.nan)%0A%20%20%20%20err_2017%20%3D%20np.where(_pos%2C%20nucleon_errors%5B%222017%22%5D%5B%22Total%22%5D%2C%20np.nan)%0A%20%20%20%20err_2025%20%3D%20np.where(_pos%2C%20nucleon_errors%5B%222025%22%5D%5B%22Total%22%5D%2C%20np.nan)%0A%20%20%20%20flux_ratio%20%3D%20flux_2017%20%2F%20flux_2025%0A%20%20%20%20err_ratio%20%3D%20err_2017%20%2F%20err_2025%0A%20%20%20%20ax1.semilogx(energy_per_nucleon%2C%20flux_ratio%2C%20%22b-%22%2C%20label%3D%22GSF2017%22)%0A%20%20%20%20ax1.fill_between(%0A%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20flux_ratio%20-%20err_2017%20%2F%20flux_2017%2C%0A%20%20%20%20%20%20%20%20flux_ratio%20%2B%20err_2017%20%2F%20flux_2017%2C%0A%20%20%20%20%20%20%20%20color%3D%22b%22%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.2%2C%0A%20%20%20%20)%0A%20%20%20%20ax1.fill_between(%0A%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%201%20-%20err_2025%20%2F%20flux_2025%2C%0A%20%20%20%20%20%20%20%201%20%2B%20err_2025%20%2F%20flux_2025%2C%0A%20%20%20%20%20%20%20%20color%3D%22gray%22%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.2%2C%0A%20%20%20%20)%0A%20%20%20%20ax1.axhline(1%2C%20color%3D%22k%22%2C%20ls%3D%22--%22%2C%20alpha%3D0.5%2C%20label%3D%22GSF2025%22)%0A%20%20%20%20ax1.set_ylabel(%22Flux%20ratio%20(2017%2F2025)%22)%0A%20%20%20%20ax1.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20ax1.legend()%0A%20%20%20%20ax1.set_xlim(1.0%2C%203000000000.0)%0A%20%20%20%20ax1.set_ylim(0.75%2C%201.25)%0A%20%20%20%20ax2.semilogx(energy_per_nucleon%2C%20err_ratio%2C%20%22k-%22)%0A%20%20%20%20ax2.axhline(1%2C%20color%3D%22k%22%2C%20ls%3D%22--%22%2C%20alpha%3D0.5)%0A%20%20%20%20ax2.set_xlabel(%22Energy%20per%20nucleon%20%5BGeV%2Fnucleon%5D%22)%0A%20%20%20%20ax2.set_ylabel(%22Error%202017%20%2F%20Error%202025%22)%0A%20%20%20%20ax2.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20ax2.set_ylim(0.0%2C%205.5)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(NOTEBOOK%2C%20SITE%2C%20mo)%3A%0A%20%20%20%20if%20SITE%3A%20%20%23%20running%20in%20the%20browser%3A%20link%20the%20sources%20published%20next%20to%20it%0A%20%20%20%20%20%20%20%20_py%20%3D%20f%22%7BSITE%7Dgallery%2F%7BNOTEBOOK%7D%2F%7BNOTEBOOK%7D.py%22%0A%20%20%20%20%20%20%20%20_ipynb%20%3D%20f%22%7BSITE%7Dgallery%2F%7BNOTEBOOK%7D%2F%7BNOTEBOOK%7D.ipynb%22%0A%20%20%20%20else%3A%20%20%23%20running%20locally%3A%20hand%20over%20the%20file%20on%20disk%0A%20%20%20%20%20%20%20%20_py%20%3D%20(mo.notebook_dir()%20%2F%20f%22%7BNOTEBOOK%7D.py%22).read_bytes()%0A%20%20%20%20%20%20%20%20_ipynb%20%3D%20None%0A%0A%20%20%20%20_buttons%20%3D%20%5B%0A%20%20%20%20%20%20%20%20mo.download(%0A%20%20%20%20%20%20%20%20%20%20%20%20_py%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20filename%3Df%22%7BNOTEBOOK%7D.py%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22marimo%20notebook%20(.py)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mimetype%3D%22text%2Fx-python%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%5D%0A%20%20%20%20if%20_ipynb%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20_buttons.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.download(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ipynb%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20filename%3Df%22%7BNOTEBOOK%7D.ipynb%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22Jupyter%20notebook%20(.ipynb)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mimetype%3D%22application%2Fx-ipynb%2Bjson%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20)%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(r%22%22%22%23%23%23%20Take%20this%20notebook%20with%20you%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(_buttons%2C%20justify%3D%22start%22)%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