model comparison 2017 vs 2025

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.15%22%0Aapp%20%3D%20marimo.App(width%3D%22full%22)%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%20if%20%22pyodide%22%20in%20sys.modules%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%20both%20plausible%20depths%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..%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%20str(_mo.notebook_location())%20%2B%20%22%2F%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_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%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%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%20%23%20TEMPORARY%20pre-publication%20marker%3A%20stamp%20every%20figure%20this%20notebook%0A%20%20%20%20%23%20shows%20with%20a%20diagonal%20PRELIMINARY%20watermark.%20Delete%20at%20the%20GSF%202026%0A%20%20%20%20%23%20release%20(guarded%20so%20a%20cell%20re-run%20does%20not%20stack%20wrappers).%0A%20%20%20%20if%20not%20getattr(plt.show%2C%20%22_gsf_preliminary%22%2C%20False)%3A%0A%20%20%20%20%20%20%20%20_orig_show%20%3D%20plt.show%0A%0A%20%20%20%20%20%20%20%20def%20_show_with_watermark(*args%2C%20**kwargs)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20plt.gcf().text(%0A%20%20%20%20%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%20%20%20%200.5%2C%0A%20%20%20%20%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%20%20%20%20%20ha%3D%22center%22%2C%0A%20%20%20%20%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%20%20%20%20%20rotation%3D30%2C%0A%20%20%20%20%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%20%20%20%20%20fontweight%3D%22bold%22%2C%0A%20%20%20%20%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%20%20%20%20%20alpha%3D0.18%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D1000%2C%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%20return%20_orig_show(*args%2C%20**kwargs)%0A%0A%20%20%20%20%20%20%20%20_show_with_watermark._gsf_preliminary%20%3D%20True%0A%20%20%20%20%20%20%20%20plt.show%20%3D%20_show_with_watermark%0A%20%20%20%20return%20GSFEnergy%2C%20GSFEnergyPerNucleon%2C%20np%2C%20plt%0A%0A%0A%40app.cell%0Adef%20_(GSFEnergy%2C%20GSFEnergyPerNucleon%2C%20np)%3A%0A%20%20%20%20energy%20%3D%20np.logspace(np.log10(0.5)%2C%20np.log10(100000000000.0)%2C%20500)%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%2C%20%222026%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%22Proton%22%2C%20%22Helium%22%2C%20%22Oxygen*%22%2C%20%22Iron*%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%222025%22%3A%20%22--%22%2C%20%222017%22%3A%20%22%3A%22%2C%20%222026%22%3A%20%22-%22%7D%0A%20%20%20%20energy_exponent%20%3D%202.7%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_exponent%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%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)%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%20_ax.set_ylim(8.0%2C%2090000.0)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%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)%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%20energy_errors%5B_v%5D%5B_group%5D%20%2F%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_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%20energy_errors%5B_v%5D%5B%22Total%22%5D%20%2F%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_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.001%2C%200.98)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%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)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots()%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%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%20_ax.set_ylim(8.0%2C%2090000.0)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%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)%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%20energy_errors%5B_v%5D%5B_group%5D%20%2F%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_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%20energy_errors%5B_v%5D%5B%22Total%22%5D%20%2F%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_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.001%2C%200.98)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%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_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)%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%20_ax.set_ylim(8.0%2C%2090000.0)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%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%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)%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%20_ax.set_ylim(80.0%2C%2090000.0)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%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)%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%20nucleon_errors%5B_v%5D%5B_group%5D%20%2F%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_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%20nucleon_errors%5B_v%5D%5B%22Total%22%5D%20%2F%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_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%201)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%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)%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%20_ax.set_ylim(8.0%2C%2090000.0)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%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%20nucleon_errors%2C%20nucleon_fluxes%2C%20plt)%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%20flux_2017%20%3D%20nucleon_fluxes%5B%222017%22%5D%5B%22Total%22%5D%0A%20%20%20%20flux_2025%20%3D%20nucleon_fluxes%5B%222025%22%5D%5B%22Total%22%5D%0A%20%20%20%20err_2017%20%3D%20nucleon_errors%5B%222017%22%5D%5B%22Total%22%5D%0A%20%20%20%20err_2025%20%3D%20nucleon_errors%5B%222025%22%5D%5B%22Total%22%5D%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%20plt.show()%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