GSF — Paper Figures

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%20Paper%20Figures%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%20GSF%202026%20%E2%80%94%20the%20paper%20figures%0A%0A%20%20%20%20The%20figures%20of%20the%20GSF%202026%20paper%20that%20the%20model%20draws%20on%20its%20own%2C%20in%0A%20%20%20%20paper%20order%2C%20at%20the%20paper's%20figure%20sizes%2C%20weightings%2C%20axis%20ranges%2C%20tick%0A%20%20%20%20locators%20and%20colours%20%E2%80%94%20ported%20from%20the%20code%20that%20produces%20them%20(the%20deck%0A%20%20%20%20in%20%60gsf_apps%2Fdeck%2F%60%2C%20%60gsf_apps%2Fpaper_figures%2F%60%20and%20the%20per-figure%20run%0A%20%20%20%20scripts%20in%20the%20fitter%20repo).%0A%0A%20%20%20%20**The%20measurements%20are%20left%20out.**%20They%20belong%20to%20the%20experiments%20that%0A%20%20%20%20published%20them%3B%20the%20paper%20is%20where%20they%20are%20shown%20together%20with%20the%20fit.%0A%20%20%20%20Everything%20else%20is%20the%20paper's%20recipe%3A%20the%20model%20is%20drawn%20as%20the%20fitted%0A%20%20%20%20**local%20interstellar%20spectrum**%20(%24J_%5Cmathrm%7BLIS%7D%24%2C%20no%20solar%20modulation)%2C%0A%20%20%20%20against%20**kinetic%20energy**%20or%20rigidity%20as%20each%20figure%20requires.%0A%0A%20%20%20%20Two%20typographic%20families%20are%20in%20play%2C%20exactly%20as%20in%20the%20paper%3A%20the%20deck%0A%20%20%20%20figures%20are%20set%20in%20a%20Times-like%20serif%2C%20the%20stand-alone%20figures%20in%20the%0A%20%20%20%20stock%20sans-serif.%20Each%20cell%20sets%20its%20own%2C%20so%20the%20order%20the%20notebook%20runs%0A%20%20%20%20in%20cannot%20change%20the%20result.%0A%0A%20%20%20%20Pick%20a%20parameter%20set%20below%20%E2%80%94%20every%20figure%20redraws.%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%22paper_figures%22%20%20%23%20used%20by%20the%20download%20buttons%20at%20the%20end%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%20%20%20%20%20%23%20crflux%20(H3a%2C%20GST%2C%20poly-gonato%2C%20...)%20for%20the%20comparison%20figures%20%E2%80%94%20a%0A%20%20%20%20%20%20%20%20%23%20pure%20python%20wheel%2C%20so%20micropip%20takes%20it%20straight%20from%20PyPI%0A%20%20%20%20%20%20%20%20await%20micropip.install(%22crflux%22)%0A%0A%20%20%20%20import%20crflux.models%20as%20cm%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%20(%0A%20%20%20%20%20%20%20%20MODEL_VERSIONS%2C%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%20GSFKineticEnergy%2C%0A%20%20%20%20%20%20%20%20GSFKineticEnergyPerNucleon%2C%0A%20%20%20%20%20%20%20%20GSFRigidity%2C%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20------------------------------------------------------------%20raster%20dpi%20--%0A%20%20%20%20%23%20These%20figsizes%20are%20paper%20column%20widths%20(3.4-3.8%20in)%2C%20so%20marimo's%0A%20%20%20%20%23%20%60figsize%20x%20100%60%20display%20rule%20makes%20them%20the%20smallest%20plots%20in%20the%0A%20%20%20%20%23%20gallery%3B%20gallery.css%20scales%20every%20figure%20up%20to%20a%20common%20width.%20dpi%20does%0A%20%20%20%20%23%20not%20affect%20the%20display%20size%20(marimo%20normalises%20it%20away)%20%E2%80%94%20only%20the%20pixel%0A%20%20%20%20%23%20count%20%E2%80%94%20so%20raise%20it%20here%20to%20keep%20the%20upscaled%20figure%20sharp%3A%20at%20150%20a%0A%20%20%20%20%23%203.4%20in%20figure%20rasterises%201020%20px%20wide%20for%20a%20900%20px%20box.%20Do%20NOT%20express%0A%20%20%20%20%23%20the%20size%20change%20as%20a%20bigger%20figsize%3A%20that%20would%20rescale%20the%20fonts%0A%20%20%20%20%23%20against%20the%20axes%20and%20break%20the%20match%20to%20the%20published%20figures.%0A%20%20%20%20plt.rcParams%5B%22figure.dpi%22%5D%20%3D%20150%0A%0A%20%20%20%20%23%20----------------------------------------------------------------%20fonts%20--%0A%20%20%20%20%23%20The%20paper's%20figures%20come%20from%20two%20families%20of%20scripts%20and%20they%20do%20NOT%0A%20%20%20%20%23%20share%20a%20font.%20deck%2Fconfig.py%20asks%20for%20Nimbus%20Roman%20(else%20a%20Times-like%0A%20%20%20%20%23%20serif)%3B%20the%20stand-alone%20scripts%20(paper_figures%2F*%2C%20the%20SM%20run%20scripts)%0A%20%20%20%20%23%20never%20touch%20the%20font%20and%20get%20matplotlib's%20stock%20DejaVu%20Sans.%20Neither%0A%20%20%20%20%23%20Nimbus%20Roman%20nor%20Times%20exists%20in%20pyodide%2C%20but%20matplotlib%20bundles%0A%20%20%20%20%23%20STIXGeneral%20%E2%80%94%20a%20Times%20metric%20clone%20%E2%80%94%20so%20the%20browser%20gets%20the%20paper's%0A%20%20%20%20%23%20letterforms.%20Every%20plot%20cell%20opens%20%60plt.rc_context(...)%60%20with%20the%20family%0A%20%20%20%20%23%20its%20figure%20was%20drawn%20in%2C%20so%20cell%20execution%20order%20cannot%20leak%20a%20font.%0A%20%20%20%20RC_SERIF%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22font.family%22%3A%20%22serif%22%2C%0A%20%20%20%20%20%20%20%20%22font.serif%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Nimbus%20Roman%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Times%20New%20Roman%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22STIXGeneral%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22DejaVu%20Serif%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22serif%22%2C%0A%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%22mathtext.fontset%22%3A%20%22stix%22%2C%0A%20%20%20%20%7D%0A%20%20%20%20RC_SANS%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22font.family%22%3A%20%22sans-serif%22%2C%0A%20%20%20%20%20%20%20%20%22font.sans-serif%22%3A%20%5B%22DejaVu%20Sans%22%2C%20%22sans-serif%22%5D%2C%0A%20%20%20%20%20%20%20%20%22mathtext.fontset%22%3A%20%22dejavusans%22%2C%0A%20%20%20%20%7D%0A%20%20%20%20%23%20the%20size%20block%20shared%20by%20paper_figures%2Fmassgroups.py%2C%20comparison.py%20and%0A%20%20%20%20%23%20the%20SM%20fidelity%2Fsampling%20scripts%0A%20%20%20%20RC_SMALL%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22font.size%22%3A%208.5%2C%0A%20%20%20%20%20%20%20%20%22axes.labelsize%22%3A%209.5%2C%0A%20%20%20%20%20%20%20%20%22axes.titlesize%22%3A%2010%2C%0A%20%20%20%20%20%20%20%20%22xtick.labelsize%22%3A%207.5%2C%0A%20%20%20%20%20%20%20%20%22ytick.labelsize%22%3A%207.5%2C%0A%20%20%20%20%20%20%20%20%22legend.fontsize%22%3A%207.5%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20%23%20-------------------------------------------%20the%20deck's%20drawing%20vocabulary%0A%20%20%20%20COLORS%20%3D%20(%22r%22%2C%20%22y%22%2C%20%22g%22%2C%20%22b%22)%20%20%23%20style.colors%0A%20%20%20%20LINES%20%3D%20(%22-%22%2C%20%22--%22%2C%20%22-.%22%2C%20%22%3A%22)%20%20%23%20style.lines%0A%20%20%20%20LABELS%20%3D%20(%22H*%22%2C%20%22He*%22%2C%20%22O*%22%2C%20%22Fe*%22)%20%20%23%20style.labels%0A%20%20%20%20ELEM_LABELS%20%3D%20(%22%24p%24%22%2C%20%22He%22%2C%20%22O%22%2C%20%22Fe%22)%20%20%23%20config.prim_label%20of%20the%20leaders%0A%20%20%20%20GROUPS%20%3D%20(%22H*%22%2C%20%22He%22%2C%20%22O*%22%2C%20%22Fe*%22)%20%20%23%20the%20model's%20mass-GROUP%20target%20keys%0A%20%20%20%20LEADER_Z%20%3D%20(1%2C%202%2C%208%2C%2026)%0A%20%20%20%20FIGSIZE_SMALL%20%3D%20(3.8%2C%203.2)%20%20%23%20config.figsize_small%0A%20%20%20%20FIGSIZE_WIDE%20%3D%20(6%2C%203.6)%20%20%23%20config.figsize_wide%0A%20%20%20%20LIGHT_COLOR%20%3D%20(1.0%2C%200.4%2C%200.0)%20%20%23%20style.light_color%0A%20%20%20%20HEAVY_COLOR%20%3D%20(0.0%2C%200.7%2C%200.7)%20%20%23%20style.heavy_color%0A%0A%20%20%20%20SYMBOLS%20%3D%20(%0A%20%20%20%20%20%20%20%20%22H%22%2C%0A%20%20%20%20%20%20%20%20%22He%22%2C%0A%20%20%20%20%20%20%20%20%22Li%22%2C%0A%20%20%20%20%20%20%20%20%22Be%22%2C%0A%20%20%20%20%20%20%20%20%22B%22%2C%0A%20%20%20%20%20%20%20%20%22C%22%2C%0A%20%20%20%20%20%20%20%20%22N%22%2C%0A%20%20%20%20%20%20%20%20%22O%22%2C%0A%20%20%20%20%20%20%20%20%22F%22%2C%0A%20%20%20%20%20%20%20%20%22Ne%22%2C%0A%20%20%20%20%20%20%20%20%22Na%22%2C%0A%20%20%20%20%20%20%20%20%22Mg%22%2C%0A%20%20%20%20%20%20%20%20%22Al%22%2C%0A%20%20%20%20%20%20%20%20%22Si%22%2C%0A%20%20%20%20%20%20%20%20%22P%22%2C%0A%20%20%20%20%20%20%20%20%22S%22%2C%0A%20%20%20%20%20%20%20%20%22Cl%22%2C%0A%20%20%20%20%20%20%20%20%22Ar%22%2C%0A%20%20%20%20%20%20%20%20%22K%22%2C%0A%20%20%20%20%20%20%20%20%22Ca%22%2C%0A%20%20%20%20%20%20%20%20%22Sc%22%2C%0A%20%20%20%20%20%20%20%20%22Ti%22%2C%0A%20%20%20%20%20%20%20%20%22V%22%2C%0A%20%20%20%20%20%20%20%20%22Cr%22%2C%0A%20%20%20%20%20%20%20%20%22Mn%22%2C%0A%20%20%20%20%20%20%20%20%22Fe%22%2C%0A%20%20%20%20%20%20%20%20%22Co%22%2C%0A%20%20%20%20%20%20%20%20%22Ni%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20def%20species_name(sid)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Return%20the%20deck's%20species%20name.%0A%0A%20%20%20%20%20%20%20%20The%20element%20symbol%2C%20but%20'p'%20for%20H%20and%20'D'%20for%20the%20deuteron%2C%20which%0A%20%20%20%20%20%20%20%20carries%20its%20own%20spline%20under%20the%20proton%20group.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20z%2C%20a%20%3D%20int(sid%5B0%5D)%2C%20round(float(sid%5B1%5D))%0A%20%20%20%20%20%20%20%20if%20z%20%3D%3D%201%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%22D%22%20if%20a%20%3D%3D%202%20else%20%22p%22%0A%20%20%20%20%20%20%20%20return%20SYMBOLS%5Bz%20-%201%5D%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%3D28%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%20draw_band(ax%2C%20x%2C%20val%2C%20err%2C%20color%2C%20line%3DNone%2C%20lw%3DNone%2C%20alpha%3D0.3%2C%20**kw)%3A%0A%20%20%20%20%20%20%20%20%22%22%22deck%2Fdraw.py%3A%20model%20central%20line%20(optional)%20%2B%20translucent%20%C2%B11%CF%83%20band.%22%22%22%0A%20%20%20%20%20%20%20%20if%20line%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20ax.plot(x%2C%20val%2C%20line%2C%20color%3Dcolor%2C%20lw%3Dlw)%0A%20%20%20%20%20%20%20%20return%20ax.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%2C%20val%20-%20err%2C%20val%20%2B%20err%2C%20color%3Dcolor%2C%20alpha%3Dalpha%2C%20linewidth%3D0%2C%20**kw%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20def%20band_line_handle(color%2C%20ls%3D%22-%22%2C%20lw%3D1.2%2C%20alpha%3D0.3)%3A%0A%20%20%20%20%20%20%20%20%22%22%22deck%2Fdraw.py%3A%20ONE%20legend%20key%20showing%20the%20band%20with%20its%20line%20on%20top.%22%22%22%0A%20%20%20%20%20%20%20%20from%20matplotlib.lines%20import%20Line2D%0A%20%20%20%20%20%20%20%20from%20matplotlib.patches%20import%20Patch%0A%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20Patch(facecolor%3Dcolor%2C%20alpha%3Dalpha%2C%20edgecolor%3D%22none%22%2C%20zorder%3D1)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20Line2D(%5B%5D%2C%20%5B%5D%2C%20color%3Dcolor%2C%20ls%3Dls%2C%20lw%3Dlw%2C%20zorder%3D2)%2C%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20def%20logx_ticks(ax%2C%20numticks%3D20)%3A%0A%20%20%20%20%20%20%20%20%22%22%22deck%2Fdraw.py%3A%20the%20dense%20log-axis%20tick%20locators%20used%20on%20every%20panel.%22%22%22%0A%20%20%20%20%20%20%20%20from%20matplotlib.ticker%20import%20LogLocator%0A%0A%20%20%20%20%20%20%20%20ax.xaxis.set_major_locator(LogLocator(numticks%3Dnumticks))%0A%20%20%20%20%20%20%20%20ax.xaxis.set_minor_locator(LogLocator(subs%3D%22auto%22%2C%20numticks%3Dnumticks))%0A%0A%20%20%20%20def%20cornertext(ax%2C%20text%2C%20loc%3D%22upper%20left%22%2C%20fontsize%3DNone)%3A%0A%20%20%20%20%20%20%20%20%22%22%22deck%2Fdraw.py%3A%20anchored%20corner%20label%20(pyik.mplext.cornertext).%22%22%22%0A%20%20%20%20%20%20%20%20from%20matplotlib.offsetbox%20import%20AnchoredText%0A%0A%20%20%20%20%20%20%20%20codes%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22upper%20left%22%3A%202%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22upper%20right%22%3A%201%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22lower%20left%22%3A%203%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22lower%20right%22%3A%204%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20at%20%3D%20AnchoredText(%0A%20%20%20%20%20%20%20%20%20%20%20%20text%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3Dcodes.get(loc%2C%202)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20prop%3D%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22size%22%3A%20fontsize%20if%20fontsize%20is%20not%20None%20else%20plt.rcParams%5B%22font.size%22%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20ax.add_artist(at)%0A%20%20%20%20%20%20%20%20return%20at%0A%0A%20%20%20%20def%20pair_error(model%2C%20x%2C%20a%2C%20b)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Return%20the%20band%20for%20the%20sum%20of%20two%20groups.%0A%0A%20%20%20%20%20%20%20%20sqrt(sum_ij%20Cov(g_i%2C%20g_j))%20over%20the%20light%20(p%2BHe)%20or%20heavy%20(O*%2BFe*)%0A%20%20%20%20%20%20%20%20pair%20%E2%80%94%20gsf_backend.py%3A%20pair_flux_error.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20return%20np.sqrt(%0A%20%20%20%20%20%20%20%20%20%20%20%20np.clip(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20sum(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.diag(np.asarray(model.covariance(_a%2C%20_b%2C%20x)%2C%20float))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20_a%20in%20(a%2C%20b)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20_b%20in%20(a%2C%20b)%0A%20%20%20%20%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%20%20%20%20%200.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20def%20gsf_label(version)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Return%20the%20figure%20label%20for%20a%20parameter-set%20version.%0A%0A%20%20%20%20%20%20%20%20The%20paper%20writes%20the%20current%20fit%20as%20'GSF2026'%3B%20keep%20that%20verbatim%20for%0A%20%20%20%20%20%20%20%20the%20released%20set%20and%20fall%20back%20to%20the%20version%20string%20otherwise.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20return%20%22GSF2026%22%20if%20version%20%3D%3D%20%222026.0%22%20else%20f%22GSF%20%7Bversion%7D%22%0A%0A%20%20%20%20%23%20-------------------------------------------------------%20data-free%20tables%20--%0A%20%20%20%20%23%20Fig.%203%20draws%20each%20species%20over%20the%20rigidity%20interval%20its%20DIRECT%20data%0A%20%20%20%20%23%20covers%2C%20not%20over%20its%20knot%20range%20%E2%80%94%20so%20without%20the%20data%20the%20extent%20has%20to%0A%20%20%20%20%23%20come%20from%20somewhere.%20These%20are%20the%20(min%2C%20max)%20measured%20rigidities%20per%0A%20%20%20%20%23%20species%20in%20the%20GSF2026%20fit%2C%20i.e.%20the%20coverage%20of%20the%20published%20spectra%2C%0A%20%20%20%20%23%20not%20the%20measurements%20themselves.%20Read%20once%20off%20the%20fitter's%20data%20zoo%0A%20%20%20%20%23%20(gsf_apps%2Fdeck%2Fflux.py%20builds%20the%20same%20dict%20at%20draw%20time).%0A%20%20%20%20DIRECT_RANGE_GV%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22p%22%3A%20(1.64192%2C%20529880)%2C%0A%20%20%20%20%20%20%20%20%22D%22%3A%20(2.78832%2C%2020.8802)%2C%0A%20%20%20%20%20%20%20%20%22He%22%3A%20(1.86055%2C%20424948)%2C%0A%20%20%20%20%20%20%20%20%22Li%22%3A%20(2.88908%2C%202071.87)%2C%0A%20%20%20%20%20%20%20%20%22Be%22%3A%20(2.87991%2C%202071.87)%2C%0A%20%20%20%20%20%20%20%20%22B%22%3A%20(1.51943%2C%2011997)%2C%0A%20%20%20%20%20%20%20%20%22C%22%3A%20(1.49172%2C%2084028.8)%2C%0A%20%20%20%20%20%20%20%20%22N%22%3A%20(1.49135%2C%202071.87)%2C%0A%20%20%20%20%20%20%20%20%22O%22%3A%20(1.53236%2C%2084028.8)%2C%0A%20%20%20%20%20%20%20%20%22F%22%3A%20(1.5842%2C%201897.96)%2C%0A%20%20%20%20%20%20%20%20%22Ne%22%3A%20(1.58835%2C%201897.96)%2C%0A%20%20%20%20%20%20%20%20%22Na%22%3A%20(1.62407%2C%201897.98)%2C%0A%20%20%20%20%20%20%20%20%22Mg%22%3A%20(1.59205%2C%201897.96)%2C%0A%20%20%20%20%20%20%20%20%22Al%22%3A%20(1.6169%2C%201897.98)%2C%0A%20%20%20%20%20%20%20%20%22Si%22%3A%20(1.64413%2C%201897.96)%2C%0A%20%20%20%20%20%20%20%20%22P%22%3A%20(1.67553%2C%2071.917)%2C%0A%20%20%20%20%20%20%20%20%22S%22%3A%20(1.64304%2C%201897.96)%2C%0A%20%20%20%20%20%20%20%20%22Cl%22%3A%20(1.68648%2C%2072.6246)%2C%0A%20%20%20%20%20%20%20%20%22Ar%22%3A%20(1.75753%2C%2077.2329)%2C%0A%20%20%20%20%20%20%20%20%22K%22%3A%20(1.74014%2C%2071.6722)%2C%0A%20%20%20%20%20%20%20%20%22Ca%22%3A%20(1.70923%2C%2069.8164)%2C%0A%20%20%20%20%20%20%20%20%22Sc%22%3A%20(1.78756%2C%2074.5279)%2C%0A%20%20%20%20%20%20%20%20%22Ti%22%3A%20(1.80754%2C%20445.819)%2C%0A%20%20%20%20%20%20%20%20%22V%22%3A%20(1.82979%2C%2077.0785)%2C%0A%20%20%20%20%20%20%20%20%22Cr%22%3A%20(1.80225%2C%20443.922)%2C%0A%20%20%20%20%20%20%20%20%22Mn%22%3A%20(1.81993%2C%2076.4821)%2C%0A%20%20%20%20%20%20%20%20%22Fe%22%3A%20(3.53461%2C%2050888.9)%2C%0A%20%20%20%20%20%20%20%20%22Co%22%3A%20(1.90252%2C%2075.9722)%2C%0A%20%20%20%20%20%20%20%20%22Ni%22%3A%20(1.84893%2C%20394.116)%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20%23%20Solar-system%20abundances%20of%20Lodders%2C%20ApJ%20591%2C%201220%20(2003)%2C%20summed%20over%20the%0A%20%20%20%20%23%20isotopes%20of%20each%20element%20(Si%20%3D%2010%5E6%20in%20the%20source%20normalisation).%20SM%20S2%0A%20%20%20%20%23%20overlays%20these%3B%20they%20are%20a%20published%20reference%20table%2C%20so%20they%20travel%20with%0A%20%20%20%20%23%20the%20notebook%20rather%20than%20with%20the%20fitter.%0A%20%20%20%20LODDERS_2003%20%3D%20(%0A%20%20%20%20%20%20%20%202.43105e10%2C%0A%20%20%20%20%20%20%20%202.34339e9%2C%0A%20%20%20%20%20%20%20%2055.47%2C%0A%20%20%20%20%20%20%20%200.7374%2C%0A%20%20%20%20%20%20%20%2017.32%2C%0A%20%20%20%20%20%20%20%207.07942e6%2C%0A%20%20%20%20%20%20%20%201.95014e6%2C%0A%20%20%20%20%20%20%20%201.41335e7%2C%0A%20%20%20%20%20%20%20%20841.1%2C%0A%20%20%20%20%20%20%20%202.14759e6%2C%0A%20%20%20%20%20%20%20%205.7510e4%2C%0A%20%20%20%20%20%20%20%201.02000e6%2C%0A%20%20%20%20%20%20%20%208.4100e4%2C%0A%20%20%20%20%20%20%20%201.0e6%2C%0A%20%20%20%20%20%20%20%208373.0%2C%0A%20%20%20%20%20%20%20%204.44866e5%2C%0A%20%20%20%20%20%20%20%205237.0%2C%0A%20%20%20%20%20%20%20%201.02500e5%2C%0A%20%20%20%20%20%20%20%203691.93%2C%0A%20%20%20%20%20%20%20%206.28704e4%2C%0A%20%20%20%20%20%20%20%2034.2%2C%0A%20%20%20%20%20%20%20%202422.0%2C%0A%20%20%20%20%20%20%20%20288.4%2C%0A%20%20%20%20%20%20%20%201.2860e4%2C%0A%20%20%20%20%20%20%20%209168.0%2C%0A%20%20%20%20%20%20%20%208.38000e5%2C%0A%20%20%20%20%20%20%20%202323.0%2C%0A%20%20%20%20%20%20%20%204.7797e4%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20COLORS%2C%0A%20%20%20%20%20%20%20%20DIRECT_RANGE_GV%2C%0A%20%20%20%20%20%20%20%20ELEM_LABELS%2C%0A%20%20%20%20%20%20%20%20FIGSIZE_SMALL%2C%0A%20%20%20%20%20%20%20%20FIGSIZE_WIDE%2C%0A%20%20%20%20%20%20%20%20GROUPS%2C%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%20GSFKineticEnergy%2C%0A%20%20%20%20%20%20%20%20GSFKineticEnergyPerNucleon%2C%0A%20%20%20%20%20%20%20%20GSFRigidity%2C%0A%20%20%20%20%20%20%20%20HEAVY_COLOR%2C%0A%20%20%20%20%20%20%20%20LABELS%2C%0A%20%20%20%20%20%20%20%20LEADER_Z%2C%0A%20%20%20%20%20%20%20%20LIGHT_COLOR%2C%0A%20%20%20%20%20%20%20%20LINES%2C%0A%20%20%20%20%20%20%20%20LODDERS_2003%2C%0A%20%20%20%20%20%20%20%20MODEL_VERSIONS%2C%0A%20%20%20%20%20%20%20%20NOTEBOOK%2C%0A%20%20%20%20%20%20%20%20RC_SANS%2C%0A%20%20%20%20%20%20%20%20RC_SERIF%2C%0A%20%20%20%20%20%20%20%20RC_SMALL%2C%0A%20%20%20%20%20%20%20%20SITE%2C%0A%20%20%20%20%20%20%20%20SYMBOLS%2C%0A%20%20%20%20%20%20%20%20band_line_handle%2C%0A%20%20%20%20%20%20%20%20cm%2C%0A%20%20%20%20%20%20%20%20cornertext%2C%0A%20%20%20%20%20%20%20%20draw_band%2C%0A%20%20%20%20%20%20%20%20gsf_label%2C%0A%20%20%20%20%20%20%20%20logx_ticks%2C%0A%20%20%20%20%20%20%20%20np%2C%0A%20%20%20%20%20%20%20%20pair_error%2C%0A%20%20%20%20%20%20%20%20plt%2C%0A%20%20%20%20%20%20%20%20show%2C%0A%20%20%20%20%20%20%20%20species_name%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(MODEL_VERSIONS%2C%20mo)%3A%0A%20%20%20%20version%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dlist(MODEL_VERSIONS)%2C%20value%3D%222026.0%22%2C%20label%3D%22parameter%20set%22%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(%5Bversion%2C%20plot_size%5D%2C%20justify%3D%22start%22%2C%20gap%3D2)%0A%20%20%20%20return%20plot_size%2C%20version%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_(%0A%20%20%20%20GSFEnergy%2C%0A%20%20%20%20GSFEnergyPerNucleon%2C%0A%20%20%20%20GSFKineticEnergy%2C%0A%20%20%20%20GSFKineticEnergyPerNucleon%2C%0A%20%20%20%20GSFRigidity%2C%0A%20%20%20%20gsf_label%2C%0A%20%20%20%20version%2C%0A)%3A%0A%20%20%20%20%23%20The%20deck%20draws%20EVERY%20model%20curve%20as%20the%20fitted%20LIS%20(gsf_backend.py%3A%0A%20%20%20%20%23%20ModelAdapter%20pins%20default_time_interval%3D%22LIS%22%20on%20all%20three%20model%0A%20%20%20%20%23%20classes)%2C%20against%20kinetic%20energy%20per%20nucleus%2C%20rigidity%2C%20or%20kinetic%0A%20%20%20%20%23%20energy%20per%20nucleon%20as%20the%20figure%20requires.%0A%20%20%20%20gsf_k%20%3D%20GSFKineticEnergy(version%3Dversion.value%2C%20default_time_interval%3D%22LIS%22)%0A%20%20%20%20gsf_r%20%3D%20GSFRigidity(version%3Dversion.value%2C%20default_time_interval%3D%22LIS%22)%0A%20%20%20%20gsf_kn%20%3D%20GSFKineticEnergyPerNucleon(%0A%20%20%20%20%20%20%20%20version%3Dversion.value%2C%20default_time_interval%3D%22LIS%22%0A%20%20%20%20)%0A%20%20%20%20%23%20The%20comparison%20figures%20(9%2C%2010%2C%2012)%20and%20the%20SM%20covariance%20figures%20work%20in%0A%20%20%20%20%23%20TOTAL%20energy%20and%20construct%20their%20models%20without%20a%20time%20interval%2C%20which%20is%0A%20%20%20%20%23%20the%20frame%20the%20external%20parametrizations%20are%20defined%20in%20%E2%80%94%20kept%20verbatim.%0A%20%20%20%20gsf_e%20%3D%20GSFEnergy(version%3Dversion.value)%0A%20%20%20%20gsf_en%20%3D%20GSFEnergyPerNucleon(version%3Dversion.value)%0A%20%20%20%20VLAB%20%3D%20gsf_label(version.value)%0A%20%20%20%20return%20VLAB%2C%20gsf_e%2C%20gsf_en%2C%20gsf_k%2C%20gsf_kn%2C%20gsf_r%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%20Fig.%201%20%E2%80%94%20Elements%20and%20mass%20groups%0A%0A%20%20%20%20A%20group%20is%20named%20after%20its%20leading%20element%3B%20groups%20divide%20the%20%24%5Cln%20A%24%20axis%0A%20%20%20%20into%20four%20roughly%20equal%20parts.%20The%20area%20of%20each%20circle%20is%20proportional%20to%0A%20%20%20%20the%20flux%20ratio%20of%20the%20element%20to%20its%20group%20leader%20(values%20tabulated%20in%20the%0A%20%20%20%20SM).%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(FIGSIZE_SMALL%2C%20LEADER_Z%2C%20RC_SERIF%2C%20gsf_r%2C%20np%2C%20plt%2C%20show%2C%20species_name)%3A%0A%20%20%20%20%23%20deck%2Fcomposition.py%3A%20mass_groups_lnA_vs_Z%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_SMALL)%0A%20%20%20%20%20%20%20%20for%20_sid%20in%20sorted(gsf_r.species%2C%20key%3Dlambda%20s%3A%20(s%5B0%5D%2C%20s%5B1%5D))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_z%2C%20_a%20%3D%20_sid%0A%20%20%20%20%20%20%20%20%20%20%20%20_lead%20%3D%20gsf_r.flux_ratio%5B_sid%5D%5B0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_sid%20!%3D%20_lead%3A%20%20%23%20ratio%20at%20the%20species'%20last%20knot%2C%20as%20in%20the%20deck%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_rmax%20%3D%20float(np.exp(np.asarray(gsf_r.kx%5B_sid%5D%2C%20float)%5B-1%5D))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_num%20%3D%20float(np.asarray(gsf_r.flux(_rmax%2C%20_sid)).ravel()%5B0%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_den%20%3D%20float(np.asarray(gsf_r.flux(_rmax%2C%20_lead)).ravel()%5B0%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ratio%20%3D%20_num%20%2F%20_den%20if%20_den%20%3E%200%20else%200.0%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%20_ratio%20%3D%201.0%0A%20%20%20%20%20%20%20%20%20%20%20%20_col%20%3D%20%22rygb%22%5BLEADER_Z.index(_lead%5B0%5D)%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_name%20%3D%20species_name(_sid)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%5Bnp.log(_a)%5D%2C%20%5B_z%5D%2C%20%22o%22%20%2B%20_col%2C%20ms%3D0.3%20%2B%206%20*%20_ratio**0.5)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_name%20%3D%3D%20%22p%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dx%2C%20_dy%20%3D%20-0.02%2C%201.4%0A%20%20%20%20%20%20%20%20%20%20%20%20elif%20_name%20%3D%3D%20%22He%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dx%2C%20_dy%20%3D%20-0.24%2C%201.0%0A%20%20%20%20%20%20%20%20%20%20%20%20elif%20_name%20in%20(%22Ne%22%2C%20%22Mg%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dx%2C%20_dy%20%3D%20-0.31%2C%200.5%0A%20%20%20%20%20%20%20%20%20%20%20%20elif%20_name%20in%20(%22C%22%2C%20%22O%22%2C%20%22Si%22%2C%20%22Fe%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dx%2C%20_dy%20%3D%20-0.23%2C%200.5%0A%20%20%20%20%20%20%20%20%20%20%20%20elif%20_name%20in%20(%22Be%22%2C%20%22S%22%2C%20%22K%22%2C%20%22Ca%22%2C%20%22Ti%22%2C%20%22Cr%22%2C%20%22Ni%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dx%2C%20_dy%20%3D%20-0.20%2C%200.5%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%20_dx%2C%20_dy%20%3D%200.05%2C%20-0.6%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.text(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.log(_a)%20%2B%20_dx%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_z%20%2B%20_dy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_name%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%20fontsize%3D%22x-small%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.text(0.44%2C%2030.55%2C%20%22H*%22%2C%20color%3D%22r%22)%0A%20%20%20%20%20%20%20%20_ax.text(1.38%2C%2030.5%2C%20%22He*%22%2C%20color%3D%22y%22)%0A%20%20%20%20%20%20%20%20_ax.text(2.37%2C%2030.5%2C%20%22O*%22%2C%20color%3D%22g%22)%0A%20%20%20%20%20%20%20%20_ax.text(3.45%2C%2030.5%2C%20%22Fe*%22%2C%20color%3D%22b%22)%0A%20%20%20%20%20%20%20%20_ax.axvspan(0%2C%201%2C%20alpha%3D0.1%2C%20color%3D%22r%22)%0A%20%20%20%20%20%20%20%20_ax.axvspan(1%2C%201.93%2C%20alpha%3D0.1%2C%20color%3D%22y%22)%0A%20%20%20%20%20%20%20%20_ax.axvspan(1.93%2C%203%2C%20alpha%3D0.1%2C%20color%3D%22g%22)%0A%20%20%20%20%20%20%20%20_ax.axvspan(3%2C%204.1%2C%20alpha%3D0.1%2C%20color%3D%22b%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(r%22Z%22)%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24%5Cln%20A%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(-0.1%2C%204.5)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(0%2C%2030)%0A%20%20%20%20show(_fig)%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%20Fig.%202%20%E2%80%94%20Groups%20against%20their%20leading%20elements%0A%0A%20%20%20%20Flux%20of%20the%20oxygen%20(top)%20and%20iron%20(bottom)%20groups%20(thick%20lines)%20compared%20to%0A%20%20%20%20the%20flux%20of%20the%20leading%20elements%20alone%20(dashed)%20and%20of%20the%20sub-leading%0A%20%20%20%20elements%20(thin%20lines).%20The%20offset%20between%20group%20and%20leading%20element%0A%20%20%20%20underlines%20the%20role%20of%20the%20sub-leading%20abundances.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20FIGSIZE_SMALL%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20cornertext%2C%0A%20%20%20%20gsf_k%2C%0A%20%20%20%20gsf_r%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Fspectra.py%3A%20flux_from_group_components%20(one%20figure%20per%20group)%0A%20%20%20%20def%20group_panel(leader_z%2C%20group_key%2C%20tag%2C%20color)%3A%0A%20%20%20%20%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_SMALL)%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.14%2C%20right%3D0.95%2C%20top%3D0.97%2C%20bottom%3D0.12)%0A%20%20%20%20%20%20%20%20%20%20%20%20_lsid%20%3D%20gsf_r.z_to_sids%5Bleader_z%5D%5B0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_kn%20%3D%20np.exp(np.asarray(gsf_r.kx%5B_lsid%5D%2C%20float))%0A%20%20%20%20%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(_kn%5B0%5D%2C%20_kn%5B-1%5D%2C%201000)%0A%20%20%20%20%20%20%20%20%20%20%20%20_sc%20%3D%20_e**3.0%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_k.flux(_e%2C%20group_key)%2C%20float)%20*%20_sc%2C%0A%20%20%20%20%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%20%20%20%20%20color%3Dcolor%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D3%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%20for%20_sid%20in%20gsf_r.species%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20gsf_r.flux_ratio%5B_sid%5D%5B0%5D%20!%3D%20_lsid%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_is_leader%20%3D%20_sid%20%3D%3D%20_lsid%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_k.flux(_e%2C%20_sid)%2C%20float)%20*%20_sc%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22--%22%20if%20_is_leader%20else%20%22-%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3Dcolor%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%20if%20_is_leader%20else%201%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D1.0%20if%20_is_leader%20else%200.2%2C%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_ax.loglog()%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24E_%5Cmathrm%7Bkin%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%22%24J_%5Cmathrm%7BLIS%7D%20%2F%20(%5Cmathrm%7BGeV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(E_%5Cmathrm%7Bkin%7D%2F%5Cmathrm%7BGeV%7D)%5E%7B3%7D%24%22%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%20cornertext(_ax%2C%20tag%2C%20loc%3D%22upper%20right%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20logx_ticks(_ax)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xlim(5%2C%2010**11.5)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_ylim(1e2%2C%201e7)%0A%20%20%20%20%20%20%20%20return%20_fig%0A%0A%20%20%20%20show(group_panel(8%2C%20%22O*%22%2C%20%22O*%22%2C%20%22g%22))%0A%20%20%20%20return%20(group_panel%2C)%0A%0A%0A%40app.cell%0Adef%20_(group_panel%2C%20show)%3A%0A%20%20%20%20show(group_panel(26%2C%20%22Fe*%22%2C%20%22Fe*%22%2C%20%22b%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%23%20Fig.%203%20%E2%80%94%20The%20leading%20elements%20versus%20rigidity%0A%0A%20%20%20%20The%20frame%20of%20the%20direct%20measurements%3A%20fitted%20LIS%20rigidity%20spectra%20of%20the%0A%20%20%20%20leading%20elements%20(%24%5Cpm1%5Csigma%24%20band)%20and%20of%20the%20sub-leading%20members%20of%20the%0A%20%20%20%20oxygen%20and%20iron%20groups%20(pale)%3B%20the%20deuteron%2C%20carried%20as%20a%20sub-leading%0A%20%20%20%20species%20of%20the%20proton%20group%2C%20is%20dimmed.%20Each%20curve%20spans%20the%20rigidity%0A%20%20%20%20range%20its%20direct%20measurements%20cover%2C%20as%20in%20the%20paper.%20In%20the%20paper%20the%0A%20%20%20%20data%20are%20drawn%20on%20top%2C%20each%20on%20its%20fitted%20energy%20scale.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20COLORS%2C%0A%20%20%20%20DIRECT_RANGE_GV%2C%0A%20%20%20%20FIGSIZE_SMALL%2C%0A%20%20%20%20LEADER_Z%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20draw_band%2C%0A%20%20%20%20gsf_r%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A%20%20%20%20species_name%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Fflux.py%3A%20flux_direct%20(show_data%3DFalse)%0A%20%20%20%20def%20_direct_color(z%2C%20leader_z)%3A%0A%20%20%20%20%20%20%20%20_i%20%3D%20LEADER_Z.index(leader_z)%0A%20%20%20%20%20%20%20%20if%20leader_z%20in%20(1%2C%202)%20or%20z%20%3D%3D%20leader_z%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20COLORS%5B_i%5D%0A%20%20%20%20%20%20%20%20return%20(0.7%2C%200.9%2C%200.7)%20if%20leader_z%20%3D%3D%208%20else%20(0.7%2C%200.7%2C%201.0)%0A%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_SMALL)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.17%2C%20bottom%3D0.12%2C%20right%3D0.95%2C%20top%3D0.85)%0A%20%20%20%20%20%20%20%20for%20_sid%20in%20sorted(gsf_r.species%2C%20key%3Dlambda%20s%3A%20(s%5B0%5D%2C%20s%5B1%5D))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_lsid%20%3D%20gsf_r.flux_ratio%5B_sid%5D%5B0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_name%20%3D%20species_name(_sid)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_name%20not%20in%20DIRECT_RANGE_GV%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20_a%2C%20_b%20%3D%20DIRECT_RANGE_GV%5B_name%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_col%20%3D%20_direct_color(_sid%5B0%5D%2C%20_lsid%5B0%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20_r%20%3D%20np.geomspace(_a%2C%20_b%2C%201000)%0A%20%20%20%20%20%20%20%20%20%20%20%20_sc%20%3D%20_r**2.6%0A%20%20%20%20%20%20%20%20%20%20%20%20_val%20%3D%20np.asarray(gsf_r.flux(_r%2C%20_sid)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_r%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_val%2C%0A%20%20%20%20%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%20%20%20%20%20color%3D_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D11%20if%20_sid%20%3D%3D%20_lsid%20else%206%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.4%20if%20_name%20%3D%3D%20%22D%22%20else%201.0%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%20if%20_sid%20%3D%3D%20_lsid%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20np.asarray(gsf_r.error(_r%2C%20_sid)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_r%2C%20_val%2C%20_err%2C%20_col%2C%20alpha%3D0.1%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20_ax.text(1200%2C%209000%2C%20%22%24p%24%22%2C%20color%3DCOLORS%5B0%5D%2C%20zorder%3D10)%0A%20%20%20%20%20%20%20%20_ax.text(1200%2C%204e2%2C%20%22He%22%2C%20color%3DCOLORS%5B1%5D%2C%20zorder%3D10)%0A%20%20%20%20%20%20%20%20_ax.text(8e2%2C%2061%2C%20%22O%22%2C%20color%3DCOLORS%5B2%5D%2C%20zorder%3D10)%0A%20%20%20%20%20%20%20%20_ax.text(4e2%2C%2010%2C%20%22Fe%22%2C%20color%3DCOLORS%5B3%5D%2C%20zorder%3D10)%0A%20%20%20%20%20%20%20%20_ax.loglog()%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24R%20%2F%20%5Cmathrm%7BGV%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%24J_%5Cmathrm%7BLIS%7D%20%2F%20(%5Cmathrm%7BGV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(R%2F%5Cmathrm%7BGV%7D)%5E%7B2.6%7D%24%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20logx_ticks(_ax)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(1%2C%201e6)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(1e-2%2C%202e4)%0A%20%20%20%20show(_fig)%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%20Fig.%205%20%E2%80%94%20Mass-group%20fractions%0A%0A%20%20%20%20Fractions%20of%20the%20four%20mass%20groups%20as%20a%20function%20of%20energy%2C%20%24%5Cpm1%5Csigma%24.%0A%20%20%20%20In%20the%20paper%20each%20panel%20carries%20the%20composition-sensitive%20air-shower%0A%20%20%20%20measurements%2C%20adjusted%20to%20the%20common%20energy%20scale%2C%20and%20the%20two%20lower%0A%20%20%20%20panels%20a%20note%20about%20summing%20elemental%20flux%20points%20%E2%80%94%20both%20are%20data%0A%20%20%20%20annotations%20and%20are%20left%20out%20here.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20COLORS%2C%0A%20%20%20%20FIGSIZE_WIDE%2C%0A%20%20%20%20GROUPS%2C%0A%20%20%20%20LABELS%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20draw_band%2C%0A%20%20%20%20gsf_k%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Ffractions.py%3A%20fractions%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(4%2C%201%2C%20figsize%3DFIGSIZE_WIDE)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(hspace%3D0%2C%20right%3D0.97%2C%20left%3D0.1%2C%20top%3D0.91%2C%20bottom%3D0.12)%0A%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(1%2C%205e11%2C%201000)%0A%20%20%20%20%20%20%20%20_all%20%3D%20%5Bnp.asarray(gsf_k.flux(_e%2C%20_g)%2C%20float)%20for%20_g%20in%20GROUPS%5D%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_g%20in%20enumerate(GROUPS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%20%3D%20_axes%5B_i%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_val%20%3D%20_all%5B_i%5D%20%2F%20np.sum(_all%2C%20axis%3D0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20np.asarray(gsf_k.fraction_error(_e%2C%20_g)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_e%2C%20_val%2C%20_err%2C%20COLORS%5B_i%5D%2C%20line%3D%22-%22%2C%20lw%3D1%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_ylabel(LABELS%5B_i%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_i%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_yticks((0%2C%200.25%2C%200.5%2C%200.75%2C%201.0))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_yticklabels((%22%240%24%22%2C%20%22%241%2F4%24%22%2C%20%22%241%2F2%24%22%2C%20%22%243%2F4%24%22%2C%20%22%241%24%22))%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%20_ax.set_yticks((0%2C%200.25%2C%200.5%2C%200.75))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_yticklabels((%22%240%24%22%2C%20%22%241%2F4%24%22%2C%20%22%241%2F2%24%22%2C%20%22%243%2F4%24%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_ylim(0%2C%201)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.semilogx()%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_i%20!%3D%203%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ax.tick_params(labelbottom%3DFalse%2C%20bottom%3DFalse)%0A%20%20%20%20%20%20%20%20%20%20%20%20logx_ticks(_ax)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xlim(1e0%2C%201e11)%0A%20%20%20%20%20%20%20%20_axes%5B3%5D.set_xlabel(r%22%24E_%5Cmathrm%7Bkin%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20show(_fig)%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%20Fig.%206%20(top)%20%E2%80%94%20The%20cosmic-ray%20flux%2C%20~1%20GeV%20to%20%2410%5E%7B11%7D%24%20GeV%0A%0A%20%20%20%20All-particle%20and%20mass-group%20fluxes%20with%20their%20%241%5Csigma%24%20uncertainties%2C%0A%20%20%20%20weighted%20by%20%24E_%5Cmathrm%7Bkin%7D%5E%7B2.6%7D%24%2C%20with%20the%20elemental%20oxygen%20and%20iron%0A%20%20%20%20spectra%20over%20the%20range%20their%20direct%20data%20cover%20and%20the%20light%20(%24p%2B%24He)%0A%20%20%20%20component%20drawn%20over%20the%20range%20the%20paper%20centres%20it%20on.%20In%20the%20paper%20this%0A%20%20%20%20panel%20carries%20the%20direct%20measurements%20of%20the%20leading%20elements%20and%20the%0A%20%20%20%20air-shower%20measurements%20of%20the%20all-particle%20flux%20and%20of%20the%20mass%20groups%2C%0A%20%20%20%20all%20adjusted%20to%20the%20common%20energy%20scale%20found%20by%20the%20fit%20and%20demodulated%0A%20%20%20%20to%20the%20LIS%20%E2%80%94%20the%20frame%20in%20which%20the%20model%20is%20drawn.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20COLORS%2C%0A%20%20%20%20ELEM_LABELS%2C%0A%20%20%20%20FIGSIZE_WIDE%2C%0A%20%20%20%20GROUPS%2C%0A%20%20%20%20LABELS%2C%0A%20%20%20%20LIGHT_COLOR%2C%0A%20%20%20%20LINES%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20band_line_handle%2C%0A%20%20%20%20draw_band%2C%0A%20%20%20%20gsf_k%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20pair_error%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Fflux.py%3A%20flux_log%20(energy_exponent%3D2.6%2C%20show_composition%3D4%2C%0A%20%20%20%20%23%20show_heavy%3DFalse)%2C%20as%20the%20deck%20registry%20calls%20it%20for%20cr_flux_traditional%0A%20%20%20%20from%20matplotlib.legend_handler%20import%20HandlerTuple%20as%20_HandlerTuple1%0A%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_WIDE%2C%20facecolor%3D%22white%22)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.1%2C%20bottom%3D0.15%2C%20right%3D0.97%2C%20top%3D0.97)%0A%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(1%2C%205e11%2C%201000)%0A%20%20%20%20%20%20%20%20_sc%20%3D%20_e**2.6%0A%20%20%20%20%20%20%20%20_keys%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_g%20in%20enumerate(GROUPS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_val%20%3D%20np.asarray(gsf_k.flux(_e%2C%20_g)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20np.asarray(gsf_k.error(_e%2C%20_g)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_e%2C%20_val%2C%20_err%2C%20COLORS%5B_i%5D%2C%20line%3DLINES%5B_i%5D%2C%20lw%3D0.75)%0A%20%20%20%20%20%20%20%20%20%20%20%20_keys.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(band_line_handle(COLORS%5B_i%5D%2C%20ls%3DLINES%5B_i%5D%2C%20lw%3D0.75)%2C%20LABELS%5B_i%5D)%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%20%23%20the%20single%20ELEMENT%20oxygen%20%2F%20iron%20spectra%20over%20their%20data%20range%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_i%20in%20(2%2C%203)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e1%20%3D%20np.geomspace(1%2C%2010**5.5%20if%20_i%20%3D%3D%202%20else%2010**5.8%2C%20200)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_k.flux(_e1%2C%20ELEM_LABELS%5B_i%5D)%2C%20float)%20*%20_e1**2.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20LINES%5B_i%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3DCOLORS%5B_i%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0.75%2C%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%23%20light%20component%20(p%20%2B%20He)%2C%20drawn%20over%20the%20range%20the%20paper%20centres%20it%20on%0A%20%20%20%20%20%20%20%20_mask%20%3D%20(_e%20%3E%202e3)%20%26%20(_e%20%3C%202e5)%0A%20%20%20%20%20%20%20%20_light%20%3D%20sum(np.asarray(gsf_k.flux(_e%2C%20_g)%2C%20float)%20for%20_g%20in%20GROUPS%5B%3A2%5D)%20*%20_sc%0A%20%20%20%20%20%20%20%20_lerr%20%3D%20pair_error(gsf_k%2C%20_e%2C%20%22H*%22%2C%20%22He%22)%20*%20_sc%0A%20%20%20%20%20%20%20%20_ax.plot(_e%2C%20np.where(_mask%2C%20_light%2C%20np.nan)%2C%20%22-%22%2C%20color%3DLIGHT_COLOR%2C%20lw%3D0.75)%0A%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_e%5B_mask%5D%2C%20_light%5B_mask%5D%2C%20_lerr%5B_mask%5D%2C%20LIGHT_COLOR%2C%20zorder%3D0)%0A%0A%20%20%20%20%20%20%20%20_tot%20%3D%20np.asarray(gsf_k.total_flux(_e)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20_terr%20%3D%20np.asarray(gsf_k.total_error(_e)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20_ax.plot(_e%2C%20_tot%2C%20%22k-%22%2C%20lw%3D1.2)%0A%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_e%2C%20_tot%2C%20_terr%2C%20%22k%22%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20_keys.append((band_line_handle(%22k%22%2C%20lw%3D1.2)%2C%20%22total%22))%0A%0A%20%20%20%20%20%20%20%20_ax.loglog()%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24E_%5Cmathrm%7Bkin%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%24J_%5Cmathrm%7BLIS%7D%20%2F%20(%5Cmathrm%7BGeV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(E_%5Cmathrm%7Bkin%7D%2F%5Cmathrm%7BGeV%7D)%5E%7B2.6%7D%24%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20handles%3D%5B_h%20for%20_h%2C%20_%20in%20_keys%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20labels%3D%5B_l%20for%20_%2C%20_l%20in%20_keys%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handler_map%3D%7Btuple%3A%20_HandlerTuple1(ndivide%3D1)%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22lower%20left%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ncol%3Dlen(_keys)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20columnspacing%3D0.7%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D%22small%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handletextpad%3D0.3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20bbox_to_anchor%3D(0.15%2C%200.08%2C%201.0%2C%200.8)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20logx_ticks(_ax)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(1e0%2C%205e11)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(5e0%2C%203e5)%0A%20%20%20%20show(_fig)%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%20Fig.%206%20(bottom)%20%E2%80%94%20The%20knee-to-suppression%20region%0A%0A%20%20%20%20The%20same%20fluxes%2C%20scaled%20by%20%24E_%5Cmathrm%7Bkin%7D%5E%7B3%7D%24%20and%20drawn%20with%20a%20**linear**%0A%20%20%20%20flux%20axis%2C%20which%20emphasizes%20the%20structures%20between%20the%20knee%20and%20the%0A%20%20%20%20suppression.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20COLORS%2C%0A%20%20%20%20FIGSIZE_WIDE%2C%0A%20%20%20%20GROUPS%2C%0A%20%20%20%20LINES%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20draw_band%2C%0A%20%20%20%20gsf_k%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Fflux.py%3A%20flux_linear%20(energy_exponent%3D3%2C%20linear%20y%20axis)%0A%20%20%20%20from%20matplotlib.ticker%20import%20EngFormatter%20as%20_EngFormatter1%0A%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_WIDE)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.1%2C%20bottom%3D0.15%2C%20right%3D0.97%2C%20top%3D0.97)%0A%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(1e1%2C%205e11%2C%201000)%0A%20%20%20%20%20%20%20%20_sc%20%3D%20_e**3.0%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_g%20in%20enumerate(GROUPS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_val%20%3D%20np.asarray(gsf_k.flux(_e%2C%20_g)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20np.asarray(gsf_k.error(_e%2C%20_g)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_e%2C%20_val%2C%20_err%2C%20COLORS%5B_i%5D%2C%20line%3DLINES%5B_i%5D%2C%20lw%3D0.75)%0A%20%20%20%20%20%20%20%20_tot%20%3D%20np.asarray(gsf_k.total_flux(_e)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20_terr%20%3D%20np.asarray(gsf_k.total_error(_e)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20draw_band(_ax%2C%20_e%2C%20_tot%2C%20_terr%2C%20%22k%22%2C%20line%3D%22-%22%2C%20lw%3D1.2%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20_ax.semilogx()%0A%20%20%20%20%20%20%20%20_ax.yaxis.set_major_formatter(_EngFormatter1())%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24E_%5Cmathrm%7Bkin%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%24J_%5Cmathrm%7BLIS%7D%20%2F%20(%5Cmathrm%7BGeV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(E_%5Cmathrm%7Bkin%7D%2F%5Cmathrm%7BGeV%7D)%5E%7B3%7D%24%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20logx_ticks(_ax)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(1e3%2C%2010**11.5)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(0%2C%206e6)%0A%20%20%20%20show(_fig)%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%20Fig.%208%20%E2%80%94%20Mean%20logarithmic%20mass%20and%20its%20variance%0A%0A%20%20%20%20%24%5Clangle%20%5Cln%20A%20%5Crangle%24%20(top)%20and%20%24%5Csigma%5E2(%5Cln%20A)%24%20(bottom)%20computed%20from%0A%20%20%20%20the%20model.%20The%20dashed%20lines%20mark%20the%20four%20group%20leaders%20(top)%20and%20the%0A%20%20%20%20variance%20of%20an%20equal%20mix%20of%20protons%20with%20each%20of%20them%20(bottom).%20In%20the%0A%20%20%20%20paper%20the%20%24%5Clangle%20%5Cln%20A%20%5Crangle%24%20data%20used%20in%20the%20fit%20and%20the%20Auger%0A%20%20%20%20%24X_%5Cmathrm%7Bmax%7D%24%20moments%20are%20overlaid.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(COLORS%2C%20RC_SERIF%2C%20draw_band%2C%20gsf_k%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20deck%2Flna.py%3A%20lna_moments%20(the%20stacked%20two-panel%20figure%20%3D%20the%20paper%20figure)%0A%20%20%20%20from%20matplotlib.ticker%20import%20LogLocator%20as%20_LogLocator1%0A%20%20%20%20from%20matplotlib.ticker%20import%20NullFormatter%20as%20_NullFormatter1%0A%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(2%2C%201%2C%20figsize%3D(3.4%2C%204.0)%2C%20sharex%3DTrue)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(hspace%3D0.1%2C%20left%3D0.15%2C%20right%3D0.85)%0A%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(1%2C%201e12%2C%20200)%0A%20%20%20%20%20%20%20%20_lna_main%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20float(np.log(gsf_k.mass_number%5Bgsf_k.z_to_sids%5B_z%5D%5B0%5D%5D))%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_z%20in%20(1%2C%202%2C%208%2C%2026)%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20_m%20%3D%20np.asarray(gsf_k.mean_lnA(_e)%2C%20float)%0A%20%20%20%20%20%20%20%20_me%20%3D%20np.asarray(gsf_k.mean_lnA_error(_e)%2C%20float)%0A%20%20%20%20%20%20%20%20_v%20%3D%20np.asarray(gsf_k.var_lnA(_e)%2C%20float)%0A%20%20%20%20%20%20%20%20_ve%20%3D%20np.asarray(gsf_k.var_lnA_error(_e)%2C%20float)%0A%0A%20%20%20%20%20%20%20%20draw_band(_axes%5B0%5D%2C%20_e%2C%20_m%2C%20_me%2C%20%22k%22%2C%20line%3D%22-%22%2C%20lw%3D1.0)%0A%20%20%20%20%20%20%20%20draw_band(_axes%5B1%5D%2C%20_e%2C%20_v%2C%20_ve%2C%20%22k%22%2C%20line%3D%22-%22%2C%20lw%3D1.0)%0A%20%20%20%20%20%20%20%20for%20_j%2C%20_lna%20in%20enumerate(_lna_main)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B0%5D.axhline(_lna%2C%20ls%3D%22--%22%2C%20color%3DCOLORS%5B_j%5D%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B0%5D.text(3e11%2C%20_lna%2C%20(%22%24p%24%22%2C%20%22He%22%2C%20%22O%22%2C%20%22Fe%22)%5B_j%5D%2C%20va%3D%22center%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_j%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20_var%20%3D%200.5%20*%20_lna**2%20-%20(0.5%20*%20_lna)%20**%202%0A%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B1%5D.axhline(_var%2C%20ls%3D%22--%22%2C%20color%3DCOLORS%5B_j%5D%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B1%5D.text(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%203e11%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_var%2C%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%20%20%20%20%22%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%22%24%5Cfrac%7B1%7D%7B2%7D%24%24p%24%2C%24%5Cfrac%7B1%7D%7B2%7D%24He%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%22%24%5Cfrac%7B1%7D%7B2%7D%24%24p%24%2C%24%5Cfrac%7B1%7D%7B2%7D%24O%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%22%24%5Cfrac%7B1%7D%7B2%7D%24%24p%24%2C%24%5Cfrac%7B1%7D%7B2%7D%24Fe%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%5B_j%5D%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)%0A%20%20%20%20%20%20%20%20_axes%5B0%5D.set_ylabel(r%22%24%5Clangle%20%5Cln%20A%20%5Crangle%24%22)%0A%20%20%20%20%20%20%20%20_axes%5B0%5D.set_ylim(-0.02%2C%20_lna_main%5B-1%5D%20%2B%200.05)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_ylabel(r%22%24%5Csigma%5E2(%5Cln%20A)%24%22)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_ylim(0%2C%204.2)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_xlabel(r%22%24E_%5Cmathrm%7Bkin%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20%20%20%20%20for%20_ax%20in%20_axes%3A%20%20%23%20deck%2Flna.py%20_energy_ticks%3A%20label%20every%20second%20decade%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xscale(%22log%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_major_locator(_LogLocator1(base%3D100))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_minor_locator(_LogLocator1(base%3D100%2C%20subs%3D(10.0%2C)))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_minor_formatter(_NullFormatter1())%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xlim(5%2C%202e11)%0A%20%20%20%20show(_fig)%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%20Fig.%209%20%E2%80%94%20Mass%20groups%20against%20other%20parametrizations%0A%0A%20%20%20%20The%20GSF%202026%20fit%20(solid%2C%20filled%20%24%5Cpm1%5Csigma%24)%20against%20the%20earlier%20GSF%202019%0A%20%20%20%20fit%20(dashed%2C%20hatched%20band)%20and%20the%20parametrizations%20of%20Gaisser%20(H3a)%2C%0A%20%20%20%20Gaisser%E2%80%93Stanev%E2%80%93Tilav%20(GST-4)%20and%20Lv%20et%20al.%20(2024)%20(drawn%20only%20above%20its%0A%20%20%20%20100%20GeV%20all-particle%20validity%20floor).%20Rows%3A%20all-particle%20flux%2C%20H%5C*%2C%20He%5C*%2C%0A%20%20%20%20O%5C*%20and%20Fe%5C*%2C%20and%20%24%5Clangle%5Cln%20A%5Crangle%24.%20Left%3A%20spectra%20weighted%20by%0A%20%20%20%20%24E%5E%7B2.7%7D%24%3B%20right%3A%20ratio%20to%20the%20GSF%202026.%20Each%20external%20composition%20is%0A%20%20%20%20summed%20into%20the%20four%20GSF%20leading%20groups%20using%20the%20fit's%20own%20charge%E2%86%92leader%0A%20%20%20%20map%2C%20and%20its%20%24%5Clangle%5Cln%20A%5Crangle%24%20uses%20the%20true%20mass%20of%20each%20component.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(RC_SERIF%2C%20RC_SMALL%2C%20VLAB%2C%20cm%2C%20gsf_e%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20gsf_apps%2Fpaper_figures%2Fmassgroups.py%2C%20verbatim%0A%20%20%20%20from%20matplotlib.legend_handler%20import%20HandlerTuple%20as%20_HandlerTuple2%0A%20%20%20%20from%20matplotlib.lines%20import%20Line2D%20as%20_Line2D2%0A%20%20%20%20from%20matplotlib.patches%20import%20Patch%20as%20_Patch2%0A%0A%20%20%20%20from%20globalsplinefit%20import%20GSFEnergy%20as%20_GSFEnergy2%0A%0A%20%20%20%20_DASHDOTDOT%20%3D%20(0%2C%20(3%2C%201.2%2C%201%2C%201.2%2C%201%2C%201.2))%0A%20%20%20%20%23%20the%20palette%20this%20figure%20was%20drawn%20with%20(globalsplinefit.plotting.COLORS)%0A%20%20%20%20_MG%20%3D%20%7B%22H*%22%3A%20%22%23d62728%22%2C%20%22He*%22%3A%20%22%23ff7f0e%22%2C%20%22O*%22%3A%20%22%232ca02c%22%2C%20%22Fe*%22%3A%20%22%231f77b4%22%7D%0A%20%20%20%20_RATIO_YLIM%20%3D%20(0.5%2C%201.7)%0A%20%20%20%20_RPAD%20%3D%200.5%20*%20(_RATIO_YLIM%5B1%5D%20-%20_RATIO_YLIM%5B0%5D)%0A%0A%20%20%20%20%23%20Lv%20et%20al.%20(2024)%2C%20ApJ%20979%2C%20225%20%E2%80%94%20a%20five-component%20power-law%20model%20with%0A%20%20%20%20%23%20rigidity-dependent%20cutoffs.%20It%20ships%20with%20the%20fitter%2C%20not%20with%20crflux%2C%20so%0A%20%20%20%20%23%20the%20parametrization%20is%20carried%20here%20as%20a%20crflux%20PrimaryFlux%20subclass%0A%20%20%20%20%23%20(gsf_apps%2Fpaper_figures%2Flv2024.py%2C%20transcribed%20from%20the%20arXiv%20source).%0A%20%20%20%20_P%2C%20_HE%2C%20_C%2C%20_N%2C%20_O%2C%20_NE%2C%20_MG12%2C%20_SI%2C%20_FE%2C%20_Z53%20%3D%20(%0A%20%20%20%20%20%20%20%2014%2C%0A%20%20%20%20%20%20%20%20402%2C%0A%20%20%20%20%20%20%20%201206%2C%0A%20%20%20%20%20%20%20%201407%2C%0A%20%20%20%20%20%20%20%201608%2C%0A%20%20%20%20%20%20%20%202010%2C%0A%20%20%20%20%20%20%20%202412%2C%0A%20%20%20%20%20%20%20%202814%2C%0A%20%20%20%20%20%20%20%205626%2C%0A%20%20%20%20%20%20%20%2012753%2C%0A%20%20%20%20)%0A%20%20%20%20_POP1%20%3D%20%7B%0A%20%20%20%20%20%20%20%20_P%3A%20(4.5e-5%2C%202.90)%2C%0A%20%20%20%20%20%20%20%20_HE%3A%20(4.4e-5%2C%202.81)%2C%0A%20%20%20%20%20%20%20%20_C%3A%20(1.1e-5%2C%202.77)%2C%0A%20%20%20%20%20%20%20%20_O%3A%20(2.1e-5%2C%202.77)%2C%0A%20%20%20%20%20%20%20%20_NE%3A%20(5.0e-6%2C%202.75)%2C%0A%20%20%20%20%20%20%20%20_MG12%3A%20(8.3e-6%2C%202.75)%2C%0A%20%20%20%20%20%20%20%20_SI%3A%20(9.4e-6%2C%202.75)%2C%0A%20%20%20%20%20%20%20%20_FE%3A%20(2.9e-5%2C%202.75)%2C%0A%20%20%20%20%7D%0A%20%20%20%20_LS%20%3D%20%7B%0A%20%20%20%20%20%20%20%20_P%3A%20(2.3e-4%2C%202.77)%2C%0A%20%20%20%20%20%20%20%20_HE%3A%20(1.5e-4%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_C%3A%20(3.2e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_O%3A%20(4.5e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_NE%3A%20(1.5e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_MG12%3A%20(1.7e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_SI%3A%20(2.0e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_FE%3A%20(1.8e-5%2C%202.45)%2C%0A%20%20%20%20%7D%0A%20%20%20%20_POP2%20%3D%20%7B%0A%20%20%20%20%20%20%20%20_P%3A%20(2.8e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_HE%3A%20(2.1e-5%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_C%3A%20(2.4e-6%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_O%3A%20(3.1e-6%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_NE%3A%20(6.0e-7%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_MG12%3A%20(7.8e-7%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_SI%3A%20(8.3e-7%2C%202.45)%2C%0A%20%20%20%20%20%20%20%20_FE%3A%20(1.6e-6%2C%202.38)%2C%0A%20%20%20%20%20%20%20%20_Z53%3A%20(1.5e-8%2C%202.10)%2C%0A%20%20%20%20%7D%0A%20%20%20%20_XGLE%20%3D%20%7B_P%3A%20(8.0e-19%2C%202.67)%2C%20_HE%3A%20(8.5e-19%2C%202.67)%7D%0A%20%20%20%20_XGHE%20%3D%20%7B_HE%3A%20(3.1e-20%2C%202.0)%2C%20_N%3A%20(1.7e-20%2C%202.0)%7D%0A%20%20%20%20_TEV%2C%20_100PEV%20%3D%201.0e3%2C%201.0e8%0A%20%20%20%20_EC_POP1%2C%20_EC_LS%2C%20_EC_POP2%20%3D%203.0e4%2C%203.6e4%2C%204.8e6%0A%20%20%20%20_EC_XGLE%2C%20_EL_XGLE%20%3D%201.0e9%2C%201.9e7%0A%20%20%20%20_EC_XGHE%2C%20_EL_XGHE%20%3D%204.0e9%2C%201.5e9%0A%20%20%20%20_PC%2C%20_YR%20%3D%203.0857e18%2C%203.1557e7%0A%20%20%20%20_KSUP%20%3D%20(500.0%20*%20_PC)%20**%202%20%2F%20(4.0%20*%205.0e27%20*%202.0e5%20*%20_YR)%0A%20%20%20%20_DELTA%20%3D%200.33%0A%0A%20%20%20%20class%20LvEtAl2024(cm.PrimaryFlux)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Lv%20et%20al.%20(2024)%20LHAASO-anchored%20multi-component%20flux%20model.%22%22%22%0A%0A%20%20%20%20%20%20%20%20def%20__init__(self%2C%20model_config%3D%22default%22%2C%20geomagnetic_cutoff%3DNone)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20super().__init__(geomagnetic_cutoff%3Dgeomagnetic_cutoff)%0A%20%20%20%20%20%20%20%20%20%20%20%20self.name%20%3D%20%22Lv%20et%20al.%20(2024)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20self.sname%20%3D%20%22Lv24%22%0A%20%20%20%20%20%20%20%20%20%20%20%20self.model_config%20%3D%20model_config%0A%20%20%20%20%20%20%20%20%20%20%20%20_ids%20%3D%20set()%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_tab%20in%20(_POP1%2C%20_LS%2C%20_POP2%2C%20_XGLE%2C%20_XGHE)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ids.update(_tab)%0A%20%20%20%20%20%20%20%20%20%20%20%20self.nucleus_ids%20%3D%20sorted(_ids)%0A%0A%20%20%20%20%20%20%20%20%40staticmethod%0A%20%20%20%20%20%20%20%20def%20_plaw(E%2C%20phi0%2C%20gamma%2C%20pivot%2C%20Z%2C%20ec)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20phi0%20*%20(E%20%2F%20pivot)%20**%20(-gamma)%20*%20np.exp(-E%20%2F%20(Z%20*%20ec))%0A%0A%20%20%20%20%20%20%20%20%40staticmethod%0A%20%20%20%20%20%20%20%20def%20_sech(x)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%201.0%20%2F%20np.cosh(np.minimum(x%2C%20700.0))%0A%0A%20%20%20%20%20%20%20%20def%20_nucleus_flux(self%2C%20corsika_id%2C%20E)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20E%20%3D%20np.atleast_1d(np.asarray(E%2C%20dtype%3Dfloat))%0A%20%20%20%20%20%20%20%20%20%20%20%20Z%2C%20_%20%3D%20self.Z_A(corsika_id)%0A%20%20%20%20%20%20%20%20%20%20%20%20flux%20%3D%20np.zeros_like(E)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20corsika_id%20in%20_POP1%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_phi0%2C%20_g%20%3D%20_POP1%5Bcorsika_id%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20flux%20%2B%3D%20self._plaw(E%2C%20_phi0%2C%20_g%2C%20_TEV%2C%20Z%2C%20_EC_POP1)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20corsika_id%20in%20_LS%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_phi0%2C%20_g%20%3D%20_LS%5Bcorsika_id%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_sup%20%3D%20np.exp(-_KSUP%20*%20((E%20%2F%20Z)%20%2F%204.0)%20**%20(-_DELTA))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20flux%20%2B%3D%20self._plaw(E%2C%20_phi0%2C%20_g%2C%20_TEV%2C%20Z%2C%20_EC_LS)%20*%20_sup%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20corsika_id%20in%20_POP2%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_phi0%2C%20_g%20%3D%20_POP2%5Bcorsika_id%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20flux%20%2B%3D%20self._plaw(E%2C%20_phi0%2C%20_g%2C%20_TEV%2C%20Z%2C%20_EC_POP2)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20corsika_id%20in%20_XGLE%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_phi0%2C%20_g%20%3D%20_XGLE%5Bcorsika_id%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20flux%20%2B%3D%20self._plaw(E%2C%20_phi0%2C%20_g%2C%20_100PEV%2C%20Z%2C%20_EC_XGLE)%20*%20self._sech(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Z%20*%20_EL_XGLE%20%2F%20E%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%20if%20corsika_id%20in%20_XGHE%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_phi0%2C%20_g%20%3D%20_XGHE%5Bcorsika_id%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20flux%20%2B%3D%20self._plaw(E%2C%20_phi0%2C%20_g%2C%20_100PEV%2C%20Z%2C%20_EC_XGHE)%20*%20self._sech(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Z%20*%20_EL_XGHE%20%2F%20E%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%20return%20flux%0A%0A%20%20%20%20_E%20%3D%20np.logspace(np.log10(5.0)%2C%20np.log10(5e11)%2C%20500)%0A%20%20%20%20_w%20%3D%20_E**2.7%0A%20%20%20%20_rows%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(%22all-particle%22%2C%20None%2C%20None%2C%20%22k%22%2C%20%220.55%22%2C%20%22flux%22)%2C%0A%20%20%20%20%20%20%20%20(%22H*%22%2C%20%22H*%22%2C%201%2C%20_MG%5B%22H*%22%5D%2C%20_MG%5B%22H*%22%5D%2C%20%22flux%22)%2C%0A%20%20%20%20%20%20%20%20(%22He*%22%2C%20%22He%22%2C%202%2C%20_MG%5B%22He*%22%5D%2C%20_MG%5B%22He*%22%5D%2C%20%22flux%22)%2C%0A%20%20%20%20%20%20%20%20(%22O*%22%2C%20%22O*%22%2C%208%2C%20_MG%5B%22O*%22%5D%2C%20_MG%5B%22O*%22%5D%2C%20%22flux%22)%2C%0A%20%20%20%20%20%20%20%20(%22Fe*%22%2C%20%22Fe*%22%2C%2026%2C%20_MG%5B%22Fe*%22%5D%2C%20_MG%5B%22Fe*%22%5D%2C%20%22flux%22)%2C%0A%20%20%20%20%20%20%20%20(r%22%24%5Clangle%5Cln%20A%5Crangle%24%22%2C%20None%2C%20None%2C%20%22k%22%2C%20%220.55%22%2C%20%22lna%22)%2C%0A%20%20%20%20%5D%0A%20%20%20%20_model_ctors%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.HillasGaisser2012(%22H3a%22)%2C%20%22H3a%22%2C%20%22%3A%22%2C%20None)%2C%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.GaisserStanevTilav(%224-gen%22)%2C%20%22GST-4%22%2C%20%22-.%22%2C%20None)%2C%0A%20%20%20%20%20%20%20%20(LvEtAl2024%2C%20%22Lv%20et%20al.%202024%22%2C%20_DASHDOTDOT%2C%20100.0)%2C%0A%20%20%20%20%5D%0A%20%20%20%20_hist%20%3D%20_GSFEnergy2(version%3D%222019%22)%0A%20%20%20%20_zungroup%20%3D%20gsf_e.z_ungroup%0A%0A%20%20%20%20def%20_safe_ratio(num%2C%20den)%3A%0A%20%20%20%20%20%20%20%20_r%20%3D%20np.full_like(den%2C%20np.nan)%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_ok%20%3D%20den%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20%20%20_r%5B_ok%5D%20%3D%20num%5B_ok%5D%20%2F%20den%5B_ok%5D%0A%20%20%20%20%20%20%20%20return%20_r%0A%0A%20%20%20%20def%20_clip_ratio(y)%3A%0A%20%20%20%20%20%20%20%20return%20np.clip(y%2C%20_RATIO_YLIM%5B0%5D%20-%20_RPAD%2C%20_RATIO_YLIM%5B1%5D%20%2B%20_RPAD)%0A%0A%20%20%20%20def%20_model_group_flux(pf%2C%20leader)%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20np.zeros_like(_E)%0A%20%20%20%20%20%20%20%20for%20_cid%20in%20pf.nucleus_ids%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_zungroup.get(pf.Z_A(_cid)%5B0%5D%2C%2026)%20%3D%3D%20leader%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_out%20%2B%3D%20pf.nucleus_flux(_cid%2C%20_E)%0A%20%20%20%20%20%20%20%20return%20_out%0A%0A%20%20%20%20def%20_model_mean_lnA(pf)%3A%0A%20%20%20%20%20%20%20%20_num%2C%20_den%20%3D%20np.zeros_like(_E)%2C%20np.zeros_like(_E)%0A%20%20%20%20%20%20%20%20for%20_cid%20in%20pf.nucleus_ids%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_f%20%3D%20pf.nucleus_flux(_cid%2C%20_E)%0A%20%20%20%20%20%20%20%20%20%20%20%20_num%20%2B%3D%20_f%20*%20np.log(pf.Z_A(_cid)%5B1%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20_den%20%2B%3D%20_f%0A%20%20%20%20%20%20%20%20return%20_safe_ratio(_num%2C%20_den)%0A%0A%20%20%20%20def%20_gsf_quantity(model%2C%20kind%2C%20target)%3A%0A%20%20%20%20%20%20%20%20if%20kind%20%3D%3D%20%22lna%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(model.mean_lnA(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(model.mean_lnA_error(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20if%20target%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(model.total_flux(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(model.total_error(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(model.flux(_E%2C%20target)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(model.error(_E%2C%20target)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20def%20_crflux_quantity(pf%2C%20kind%2C%20leader%2C%20target)%3A%0A%20%20%20%20%20%20%20%20if%20kind%20%3D%3D%20%22lna%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20_model_mean_lnA(pf)%0A%20%20%20%20%20%20%20%20if%20target%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20np.asarray(pf.total_flux(_E)%2C%20float)%0A%20%20%20%20%20%20%20%20return%20_model_group_flux(pf%2C%20leader)%0A%0A%20%20%20%20def%20_hist_band(ax%2C%20x%2C%20lo%2C%20hi%2C%20col)%3A%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%0A%20%20%20%20%20%20%20%20%20%20%20%20lo%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20hi%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20facecolor%3D%22none%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20edgecolor%3Dcol%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20hatch%3D%22%2F%2F%2F%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20linewidth%3D0.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.55%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D2.4%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20ax.plot(x%2C%20lo%2C%20color%3Dcol%2C%20ls%3D%22-%22%2C%20lw%3D0.5%2C%20alpha%3D0.6%2C%20zorder%3D2.5)%0A%20%20%20%20%20%20%20%20ax.plot(x%2C%20hi%2C%20color%3Dcol%2C%20ls%3D%22-%22%2C%20lw%3D0.5%2C%20alpha%3D0.6%2C%20zorder%3D2.5)%0A%0A%20%20%20%20_models%20%3D%20%5B(ctor()%2C%20tag%2C%20ls%2C%20floor)%20for%20ctor%2C%20tag%2C%20ls%2C%20floor%20in%20_model_ctors%5D%0A%0A%20%20%20%20with%20plt.rc_context(%7B**RC_SERIF%2C%20**RC_SMALL%7D)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(%0A%20%20%20%20%20%20%20%20%20%20%20%206%2C%0A%20%20%20%20%20%20%20%20%20%20%20%202%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20figsize%3D(6.9%2C%207.2)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sharex%3DTrue%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20gridspec_kw%3D%7B%22hspace%22%3A%200.09%2C%20%22wspace%22%3A%200.06%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20for%20_i%2C%20(_label%2C%20_target%2C%20_leader%2C%20_col%2C%20_bandcol%2C%20_kind)%20in%20enumerate(_rows)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_axL%2C%20_axR%20%3D%20_axes%5B_i%2C%200%5D%2C%20_axes%5B_i%2C%201%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_weight%20%3D%201.0%20if%20_kind%20%3D%3D%20%22lna%22%20else%20_w%0A%20%20%20%20%20%20%20%20%20%20%20%20_c26%2C%20_e26%20%3D%20_gsf_quantity(gsf_e%2C%20_kind%2C%20_target)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axL.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_E%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_weight%20*%20(_c26%20-%20_e26)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_weight%20*%20(_c26%20%2B%20_e26)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_bandcol%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.30%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D1%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%20_rel26%20%3D%20_safe_ratio(_e26%2C%20_c26)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axR.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_E%2C%201%20-%20_rel26%2C%201%20%2B%20_rel26%2C%20color%3D_bandcol%2C%20alpha%3D0.30%2C%20lw%3D0%2C%20zorder%3D1%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%20_axR.axhline(1%2C%20color%3D%22k%22%2C%20lw%3D0.8%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_cH%2C%20_eH%20%3D%20_gsf_quantity(_hist%2C%20_kind%2C%20_target)%0A%20%20%20%20%20%20%20%20%20%20%20%20_hist_band(_axL%2C%20_E%2C%20_weight%20*%20(_cH%20-%20_eH)%2C%20_weight%20*%20(_cH%20%2B%20_eH)%2C%20_col)%0A%20%20%20%20%20%20%20%20%20%20%20%20_hist_band(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axR%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_E%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_clip_ratio(_safe_ratio(_cH%20-%20_eH%2C%20_c26))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_clip_ratio(_safe_ratio(_cH%20%2B%20_eH%2C%20_c26))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_col%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%20_axL.plot(_E%2C%20_weight%20*%20_c26%2C%20color%3D_col%2C%20ls%3D%22-%22%2C%20lw%3D1.7%2C%20zorder%3D5)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axL.plot(_E%2C%20_weight%20*%20_cH%2C%20color%3D_col%2C%20ls%3D%22--%22%2C%20lw%3D1.2%2C%20zorder%3D4)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axR.plot(_E%2C%20_safe_ratio(_cH%2C%20_c26)%2C%20color%3D_col%2C%20ls%3D%22--%22%2C%20lw%3D1.2%2C%20zorder%3D4)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_pf%2C%20_tag%2C%20_ls%2C%20_floor%20in%20_models%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_y%20%3D%20_crflux_quantity(_pf%2C%20_kind%2C%20_leader%2C%20_target)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_m%20%3D%20np.ones_like(_E%2C%20dtype%3Dbool)%20if%20_floor%20is%20None%20else%20(_floor%20%3C%3D%20_E)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axL.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_E%5B_m%5D%2C%20(_weight%20*%20_y)%5B_m%5D%2C%20color%3D_col%2C%20ls%3D_ls%2C%20lw%3D1.1%2C%20zorder%3D3%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_axR.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_E%5B_m%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_safe_ratio(_y%2C%20_c26)%5B_m%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3D_ls%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D3%2C%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_axL.set_xscale(%22log%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_kind%20%3D%3D%20%22flux%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axL.set_yscale(%22log%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axL.text(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.035%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.94%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_label%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20transform%3D_axL.transAxes%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D10%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_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20va%3D%22top%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22left%22%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%20_axR.text(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.035%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.94%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_label%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20transform%3D_axR.transAxes%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D9%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_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20va%3D%22top%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22left%22%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%20_axR.yaxis.tick_right()%0A%20%20%20%20%20%20%20%20%20%20%20%20_axR.yaxis.set_label_position(%22right%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axR.tick_params(labelright%3DTrue)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axR.set_ylim(*_RATIO_YLIM)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_kind%20%3D%3D%20%22lna%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axL.set_ylim(0.0%2C%204.2)%0A%20%20%20%20%20%20%20%20%20%20%20%20elif%20_target%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axL.set_ylim(2e2%2C%202e5)%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%20_axL.set_ylim(1.5e2%2C%205e4)%0A%20%20%20%20%20%20%20%20_axes%5B0%2C%200%5D.set_xlim(5.0%2C%205e11)%0A%20%20%20%20%20%20%20%20_axes%5B2%2C%200%5D.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%24E%5E%7B2.7%7D%20J(E)%24%20%20%5BGeV%24%5E%7B1.7%7D%24%20m%24%5E%7B-2%7D%24%20s%24%5E%7B-1%7D%24%20sr%24%5E%7B-1%7D%24%5D%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_axes%5B5%2C%200%5D.set_ylabel(r%22%24%5Clangle%5Cln%20A%5Crangle%24%22)%0A%20%20%20%20%20%20%20%20_axes%5B2%2C%201%5D.set_ylabel(f%22model%20%2F%20%7BVLAB%7D%22)%0A%20%20%20%20%20%20%20%20_axes%5B5%2C%200%5D.set_xlabel(r%22%24E%24%20%20%5BGeV%5D%22)%0A%20%20%20%20%20%20%20%20_axes%5B5%2C%201%5D.set_xlabel(r%22%24E%24%20%20%5BGeV%5D%22)%0A%20%20%20%20%20%20%20%20_handles%20%3D%20%5B%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%20%20%20%20%20_Patch2(facecolor%3D%220.55%22%2C%20alpha%3D0.5)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_Line2D2(%5B0%5D%2C%20%5B0%5D%2C%20color%3D%22k%22%2C%20ls%3D%22-%22%2C%20lw%3D1.6)%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%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_Patch2(facecolor%3D%22none%22%2C%20edgecolor%3D%22k%22%2C%20hatch%3D%22%2F%2F%2F%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_Line2D2(%5B0%5D%2C%20%5B0%5D%2C%20color%3D%22k%22%2C%20ls%3D%22--%22%2C%20lw%3D1.2)%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%20_Line2D2(%5B0%5D%2C%20%5B0%5D%2C%20color%3D%22k%22%2C%20ls%3D%22%3A%22%2C%20lw%3D1.1)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_Line2D2(%5B0%5D%2C%20%5B0%5D%2C%20color%3D%22k%22%2C%20ls%3D%22-.%22%2C%20lw%3D1.1)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_Line2D2(%5B0%5D%2C%20%5B0%5D%2C%20color%3D%22k%22%2C%20ls%3D_DASHDOTDOT%2C%20lw%3D1.1)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20_fig.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20_handles%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5BVLAB%2C%20%22GSF2019%22%2C%20%22H3a%22%2C%20%22GST-4%22%2C%20%22Lv%20et%20al.%202024%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22upper%20center%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20bbox_to_anchor%3D(0.5%2C%200.985)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ncol%3D5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D2.2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20columnspacing%3D1.3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handletextpad%3D0.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handler_map%3D%7Btuple%3A%20_HandlerTuple2(ndivide%3D1)%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(top%3D0.955)%0A%20%20%20%20show(_fig)%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%20Fig.%2010%20%E2%80%94%20All-particle%20flux%20against%20common%20parametrizations%0A%0A%20%20%20%20The%20GSF%20all-particle%20flux%20(black%2C%20%24%5Cpm1%5Csigma%24)%20compared%20with%20the%0A%20%20%20%20broken-power-law%20parametrizations%20in%20common%20use%2C%20and%20with%20the%20earlier%20GSF%0A%20%20%20%202019%20fit%20(hatched).%20Bottom%3A%20ratio%20to%20the%20GSF.%20Zatsepin%E2%80%93Sokolskaya%20is%0A%20%20%20%20truncated%20at%20%245%5Ctimes10%5E%7B7%7D%24%20GeV%20and%20Gaisser%E2%80%93Honda%20at%20%2410%5E%7B7%7D%24%20GeV%2C%20the%0A%20%20%20%20limits%20of%20their%20validity.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(GSFEnergy%2C%20RC_SERIF%2C%20RC_SMALL%2C%20VLAB%2C%20cm%2C%20gsf_e%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20gsf_apps%2Fpaper_figures%2Fcomparison.py%3A%20_two_panel%2C%20all-particle%0A%20%20%20%20from%20matplotlib.legend_handler%20import%20HandlerTuple%0A%0A%20%20%20%20MODEL_CTORS%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.HillasGaisser2012(%22H3a%22)%2C%20%22H3a%22)%2C%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.HillasGaisser2012(%22H4a%22)%2C%20%22H4a%22)%2C%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.GaisserStanevTilav(%224-gen%22)%2C%20%22GST%20(4-gen)%22)%2C%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.PolyGonato()%2C%20%22poly-gonato%22)%2C%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.GaisserHonda()%2C%20%22Gaisser-Honda%22)%2C%0A%20%20%20%20%20%20%20%20(lambda%3A%20cm.ZatsepinSokolskaya(%22default%22)%2C%20%22Zatsepin-Sokolskaya%22)%2C%0A%20%20%20%20%5D%0A%20%20%20%20TRUNC%20%3D%20%7B%22Zatsepin-Sokolskaya%22%3A%205e7%2C%20%22Gaisser-Honda%22%3A%201e7%7D%0A%20%20%20%20OKABE_ITO%20%3D%20%5B%22%23E69F00%22%2C%20%22%2356B4E9%22%2C%20%22%23009E73%22%2C%20%22%230072B2%22%2C%20%22%23D55E00%22%2C%20%22%23CC79A7%22%5D%0A%20%20%20%20HIST_COL%20%3D%20%220.40%22%0A%0A%20%20%20%20def%20two_panel(%0A%20%20%20%20%20%20%20%20x%2C%20gsf_flux%2C%20gsf_err%2C%20model_curves%2C%20xlabel%2C%20ylabel%2C%20ylim1%2C%20ylim2%2C%20hist%2C%20wexp%3D2.7%0A%20%20%20%20)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20(_ax1%2C%20_ax2)%20%3D%20plt.subplots(%0A%20%20%20%20%20%20%20%20%20%20%20%202%2C%0A%20%20%20%20%20%20%20%20%20%20%20%201%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20figsize%3D(3.4%2C%204.2)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sharex%3DTrue%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20gridspec_kw%3D%7B%22height_ratios%22%3A%20%5B2.0%2C%201.0%5D%2C%20%22hspace%22%3A%200.07%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_w%20%3D%20x**wexp%0A%20%20%20%20%20%20%20%20_band%20%3D%20_ax1.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_w%20*%20(gsf_flux%20-%20gsf_err)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_w%20*%20(gsf_flux%20%2B%20gsf_err)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%220.55%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20(_line%2C)%20%3D%20_ax1.loglog(x%2C%20_w%20*%20gsf_flux%2C%20color%3D%22k%22%2C%20lw%3D1.6)%0A%20%20%20%20%20%20%20%20_handles%2C%20_labels%20%3D%20%5B(_band%2C%20_line)%5D%2C%20%5BVLAB%5D%0A%20%20%20%20%20%20%20%20_yh%2C%20_eh%20%3D%20hist%0A%20%20%20%20%20%20%20%20_b19%20%3D%20_ax1.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_w%20*%20(_yh%20-%20_eh)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_w%20*%20(_yh%20%2B%20_eh)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20facecolor%3D%22none%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20hatch%3D%22%2F%2F%2F%2F%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20edgecolor%3DHIST_COL%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0%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)%0A%20%20%20%20%20%20%20%20(_l19%2C)%20%3D%20_ax1.loglog(x%2C%20_w%20*%20_yh%2C%20color%3DHIST_COL%2C%20lw%3D1.3%2C%20ls%3D%22--%22)%0A%20%20%20%20%20%20%20%20_ax2.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(_yh%20-%20_eh)%20%2F%20gsf_flux%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(_yh%20%2B%20_eh)%20%2F%20gsf_flux%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20facecolor%3D%22none%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20hatch%3D%22%2F%2F%2F%2F%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20edgecolor%3DHIST_COL%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0%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)%0A%20%20%20%20%20%20%20%20_ax2.semilogx(x%2C%20_yh%20%2F%20gsf_flux%2C%20color%3DHIST_COL%2C%20lw%3D1.3%2C%20ls%3D%22--%22)%0A%20%20%20%20%20%20%20%20_handles.append((_b19%2C%20_l19))%0A%20%20%20%20%20%20%20%20_labels.append(%22GSF2019%22)%0A%20%20%20%20%20%20%20%20for%20(_y%2C%20_tag)%2C%20_c%20in%20zip(model_curves%2C%20OKABE_ITO%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20(_ln%2C)%20%3D%20_ax1.loglog(x%2C%20_w%20*%20_y%2C%20lw%3D1.2%2C%20color%3D_c)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax2.semilogx(x%2C%20_y%20%2F%20gsf_flux%2C%20lw%3D1.2%2C%20color%3D_c)%0A%20%20%20%20%20%20%20%20%20%20%20%20_handles.append(_ln)%0A%20%20%20%20%20%20%20%20%20%20%20%20_labels.append(_tag)%0A%20%20%20%20%20%20%20%20_rel%20%3D%20gsf_err%20%2F%20gsf_flux%0A%20%20%20%20%20%20%20%20_ax2.fill_between(x%2C%201%20-%20_rel%2C%201%20%2B%20_rel%2C%20color%3D%220.55%22%2C%20alpha%3D0.5%2C%20lw%3D0)%0A%20%20%20%20%20%20%20%20_ax2.axhline(1%2C%20color%3D%22k%22%2C%20lw%3D0.9)%0A%20%20%20%20%20%20%20%20_ax1.set_ylabel(ylabel)%0A%20%20%20%20%20%20%20%20_ax1.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20_handles%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_labels%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D6.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ncol%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22lower%20left%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D1.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20columnspacing%3D1.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handletextpad%3D0.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20labelspacing%3D0.3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20borderpad%3D0.4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handler_map%3D%7Btuple%3A%20HandlerTuple(ndivide%3D1)%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax2.set_ylabel(f%22model%20%2F%20%7BVLAB%7D%22)%0A%20%20%20%20%20%20%20%20_ax2.set_xlabel(xlabel)%0A%20%20%20%20%20%20%20%20_ax2.set_ylim(*ylim2)%0A%20%20%20%20%20%20%20%20_ax1.set_ylim(*ylim1)%0A%20%20%20%20%20%20%20%20return%20_fig%0A%0A%20%20%20%20with%20plt.rc_context(%7B**RC_SERIF%2C%20**RC_SMALL%7D)%3A%0A%20%20%20%20%20%20%20%20_E%20%3D%20np.logspace(1%2C%2011%2C%20400)%0A%20%20%20%20%20%20%20%20_curves%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_ctor%2C%20_tag%20in%20MODEL_CTORS%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_y%20%3D%20np.asarray(_ctor().total_flux(_E)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_tag%20in%20TRUNC%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_y%20%3D%20np.where(TRUNC%5B_tag%5D%20%3E%3D%20_E%2C%20_y%2C%20np.nan)%0A%20%20%20%20%20%20%20%20%20%20%20%20_curves.append((_y%2C%20_tag))%0A%20%20%20%20%20%20%20%20_g19%20%3D%20GSFEnergy(version%3D%222019%22)%0A%20%20%20%20%20%20%20%20_fig10%20%3D%20two_panel(%0A%20%20%20%20%20%20%20%20%20%20%20%20_E%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_e.total_flux(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_e.total_error(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_curves%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20xlabel%3D%22E%20(GeV)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ylabel%3Dr%22%24E%5E%7B2.7%7D%20J(E)%24%20%5BGeV%24%5E%7B1.7%7D%24%20m%24%5E%7B-2%7D%24%20s%24%5E%7B-1%7D%24%20sr%24%5E%7B-1%7D%24%5D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ylim1%3D(8e2%2C%208e4)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ylim2%3D(0.51%2C%201.49)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20hist%3D(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(_g19.total_flux(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(_g19.total_error(_E)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20show(_fig10)%0A%20%20%20%20return%20MODEL_CTORS%2C%20TRUNC%2C%20two_panel%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%20Fig.%2011%20%E2%80%94%20Nucleon%20flux%0A%0A%20%20%20%20Nucleon%20flux%20versus%20energy%20per%20nucleon%20from%20the%20GSF%20(thick%20line%2C%0A%20%20%20%20%24%5Cpm1%5Csigma%24%20band)%20with%20the%20contributions%20of%20the%20four%20mass%20groups.%20The%0A%20%20%20%20nucleon%20flux%20is%20not%20directly%20observed%3A%20it%20is%20an%20elementary%20derived%0A%20%20%20%20quantity%20of%20the%20model%2C%20with%20uncertainties%20propagated%20through%20the%20full%0A%20%20%20%20covariance.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(COLORS%2C%20GROUPS%2C%20LABELS%2C%20LINES%2C%20RC_SERIF%2C%20gsf_kn%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20deck%2Fspectra.py%3A%20nucleon_flux%20(energy_exponent%3D2.7%2C%20show_lines%3DTrue)%0A%20%20%20%20from%20matplotlib.ticker%20import%20FixedLocator%20as%20_FixedLocator1%0A%20%20%20%20from%20matplotlib.ticker%20import%20NullFormatter%20as%20_NullFormatter2%0A%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(3.4%2C%202.9))%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.13%2C%20right%3D0.95%2C%20top%3D0.95)%0A%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(1e-1%2C%201e12%2C%201000)%0A%20%20%20%20%20%20%20%20_sc%20%3D%20_e**2.7%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_g%20in%20enumerate(GROUPS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_fl%20%3D%20np.asarray(gsf_kn.flux(_e%2C%20_g)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%20_fl%20*%20_sc%2C%20LINES%5B_i%5D%2C%20color%3DCOLORS%5B_i%5D%2C%20lw%3D0.75%2C%20label%3DLABELS%5B_i%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_tot%20%3D%20np.asarray(gsf_kn.total_flux(_e)%2C%20float)%0A%20%20%20%20%20%20%20%20_terr%20%3D%20np.asarray(gsf_kn.total_error(_e)%2C%20float)%0A%20%20%20%20%20%20%20%20_ax.plot(_e%2C%20_tot%20*%20_sc%2C%20%22-k%22%2C%20lw%3D1.2%2C%20label%3D%22Total%22)%0A%20%20%20%20%20%20%20%20_ax.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%20(_tot%20-%20_terr)%20*%20_sc%2C%20(_tot%20%2B%20_terr)%20*%20_sc%2C%20color%3D%22k%22%2C%20alpha%3D0.3%2C%20lw%3D0%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.loglog()%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24E_%7B%5Cmathrm%7Bkin%2Cn%7D%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%24J_%7B%5Cmathrm%7BLIS%2Cn%7D%7D%20%2F%20(%5Cmathrm%7BGeV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(E_%7B%5Cmathrm%7Bkin%2Cn%7D%7D%2F%5Cmathrm%7BGeV%7D)%5E%7B2.7%7D%24%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22upper%20right%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ncol%3D5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20columnspacing%3D1.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D1.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handletextpad%3D0.4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D%22small%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%23%20major%20ticks%20every%20second%20power%20of%20ten%20(labelled)%2C%20minor%20at%20every%20power%0A%20%20%20%20%20%20%20%20_ax.xaxis.set_major_locator(%0A%20%20%20%20%20%20%20%20%20%20%20%20_FixedLocator1(%5B10.0**_k%20for%20_k%20in%20range(0%2C%2012%2C%202)%5D)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.xaxis.set_minor_locator(_FixedLocator1(%5B10.0**_k%20for%20_k%20in%20range(0%2C%2012)%5D))%0A%20%20%20%20%20%20%20%20_ax.xaxis.set_minor_formatter(_NullFormatter2())%0A%20%20%20%20%20%20%20%20_ax.set_xlim(1%2C%2010**11.5)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(1e1%2C%201.001e5)%0A%20%20%20%20show(_fig)%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%20Fig.%2012%20%E2%80%94%20Nucleon%20flux%20against%20the%20same%20parametrizations%0A%0A%20%20%20%20Bottom%3A%20ratio%20to%20the%20GSF%3B%20the%20grey%20band%20is%20the%20GSF%20%24%5Cpm1%5Csigma%24.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20GSFEnergyPerNucleon%2C%0A%20%20%20%20MODEL_CTORS%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20RC_SMALL%2C%0A%20%20%20%20TRUNC%2C%0A%20%20%20%20gsf_en%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A%20%20%20%20two_panel%2C%0A)%3A%0A%20%20%20%20%23%20gsf_apps%2Fpaper_figures%2Fcomparison.py%3A%20_two_panel%2C%20nucleon%20flux%0A%20%20%20%20with%20plt.rc_context(%7B**RC_SERIF%2C%20**RC_SMALL%7D)%3A%0A%20%20%20%20%20%20%20%20_EN%20%3D%20np.logspace(1%2C%209.5%2C%20400)%0A%20%20%20%20%20%20%20%20_ncurves%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_ctor%2C%20_tag%20in%20MODEL_CTORS%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_%2C%20_p%2C%20_n%20%3D%20_ctor().p_and_n_flux(_EN)%0A%20%20%20%20%20%20%20%20%20%20%20%20_y%20%3D%20np.asarray(_p%20%2B%20_n%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_tag%20in%20TRUNC%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_y%20%3D%20np.where(TRUNC%5B_tag%5D%20%3E%3D%20_EN%2C%20_y%2C%20np.nan)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ncurves.append((_y%2C%20_tag))%0A%20%20%20%20%20%20%20%20_g19n%20%3D%20GSFEnergyPerNucleon(version%3D%222019%22)%0A%20%20%20%20%20%20%20%20_fig12%20%3D%20two_panel(%0A%20%20%20%20%20%20%20%20%20%20%20%20_EN%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_en.total_flux(_EN)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_en.total_error(_EN)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_ncurves%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20xlabel%3Dr%22%24E_N%24%20(GeV)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ylabel%3Dr%22%24E_N%5E%7B2.7%7D%20J_N%24%20%5BGeV%24%5E%7B1.7%7D%24%20m%24%5E%7B-2%7D%24%20s%24%5E%7B-1%7D%24%20sr%24%5E%7B-1%7D%24%5D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ylim1%3D(2e2%2C%205e4)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ylim2%3D(0.3%2C%201.9)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20hist%3D(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(_g19n.total_flux(_EN)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(_g19n.total_error(_EN)%2C%20float)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20show(_fig12)%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%20SM%20S1%20%E2%80%94%20The%20proton-leader%20spline%0A%0A%20%20%20%20The%20fitted%20flux%20(thick%20line)%2C%20the%20knot%20positions%20(dots)%2C%20and%20each%0A%20%20%20%20single-knot%20basis%20contribution%20(thin%20dashed%2Fdotted)%2C%20which%20sum%20to%20the%20flux.%0A%0A%20%20%20%20%3E%20The%20published%20version%20of%20this%20panel%20shows%20only%20the%20flux%20and%20the%20knots%3A%0A%20%20%20%20%3E%20the%20deck%20draws%20each%20basis%20contribution%20by%20zeroing%20every%20other%20amplitude%0A%20%20%20%20%3E%20and%20re-evaluating%2C%20but%20the%20frozen%20model%20it%20renders%20from%20ignores%20the%0A%20%20%20%20%3E%20parameter%20vector%2C%20so%20all%2025%20%22basis%22%20curves%20come%20out%20equal%20to%20the%20full%0A%20%20%20%20%3E%20spline%20and%20hide%20under%20it.%20Here%20the%20contributions%20are%20taken%20from%20the%0A%20%20%20%20%3E%20spline's%20own%20Jacobian%2C%20so%20the%20figure%20shows%20what%20its%20caption%20describes.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20COLORS%2C%0A%20%20%20%20ELEM_LABELS%2C%0A%20%20%20%20FIGSIZE_SMALL%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20cornertext%2C%0A%20%20%20%20gsf_r%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Fspectra.py%3A%20spline_plots%20(proton%20panel%2C%20rigidity_exponent%3D3)%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_SMALL)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.15%2C%20bottom%3D0.12%2C%20right%3D0.95%2C%20top%3D0.98)%0A%20%20%20%20%20%20%20%20_sid%20%3D%20gsf_r.z_to_sids%5B1%5D%5B0%5D%0A%20%20%20%20%20%20%20%20_kn%20%3D%20np.exp(np.asarray(gsf_r.kx%5B_sid%5D%2C%20float))%0A%20%20%20%20%20%20%20%20_r%20%3D%20np.geomspace(_kn%5B0%5D%2C%20_kn%5B-1%5D%2C%201000)%0A%20%20%20%20%20%20%20%20_sc%20%3D%20_r**3.0%0A%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20_r%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_r.flux(_r%2C%20_sid)%2C%20float)%20*%20_sc%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%20color%3DCOLORS%5B0%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.5%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_knots%20%3D%20np.unique(_kn)%0A%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20_knots%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20np.asarray(gsf_r.flux(_knots%2C%20_sid)%2C%20float)%20*%20_knots**3.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22o%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3DCOLORS%5B0%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20markersize%3D3%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%23%20each%20basis%20contribution%20%3D%20one%20amplitude%20x%20its%20basis%20function%3B%20the%0A%20%20%20%20%20%20%20%20%23%20jacobian%20columns%20ARE%20the%20basis%20functions%2C%20so%20this%20is%20the%20deck's%0A%20%20%20%20%20%20%20%20%23%20%22all%20other%20coefficients%20zeroed%22%20curve%20without%20touching%20the%20parameters%0A%20%20%20%20%20%20%20%20_jac%20%3D%20np.asarray(gsf_r.jacobian(_r%2C%20_sid)%2C%20float)%0A%20%20%20%20%20%20%20%20_amp%20%3D%20np.asarray(gsf_r.pars%5B_sid%5D%2C%20float)%5B%3A%20_jac.shape%5B1%5D%5D%0A%20%20%20%20%20%20%20%20for%20_i%20in%20range(_jac.shape%5B1%5D)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_r%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_jac%5B%3A%2C%20_i%5D%20*%20_amp%5B_i%5D%20*%20_sc%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22--%22%20if%20_i%20%25%202%20%3D%3D%200%20else%20%22%3A%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3DCOLORS%5B0%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0.75%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.loglog()%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24R%20%2F%20%5Cmathrm%7BGV%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%24J_%5Cmathrm%7BLIS%7D%20%2F%20(%5Cmathrm%7BGV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(R%2F%5Cmathrm%7BGV%7D)%5E%7B3%7D%24%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20cornertext(_ax%2C%20ELEM_LABELS%5B0%5D%2C%20loc%3D%22upper%20left%22)%0A%20%20%20%20%20%20%20%20logx_ticks(_ax%2C%20numticks%3D30)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(_r%5B0%5D%2C%20_r%5B-1%5D)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(1e1%2C%205e6)%0A%20%20%20%20show(_fig)%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%20SM%20S2%20%E2%80%94%20Elemental%20abundances%0A%0A%20%20%20%20Abundances%20at%20fixed%20kinetic%20energy%20per%20nucleon%20from%20the%20GSF%20(coloured%0A%20%20%20%20lines%2C%20%241%24%E2%80%93%2410%5E%7B9%7D%24%20GeV%2Fn)%2C%20normalized%20to%20%24%5Cmathrm%7BSi%7D%5Cequiv100%24%2C%20compared%0A%20%20%20%20with%20the%20solar-system%20abundances%20of%20Lodders%20(2003).%20With%20the%20constant-ratio%0A%20%20%20%20extrapolation%20of%20previous%20releases%20all%20curves%20above%20~1%20TeV%2Fn%20would%0A%20%20%20%20coincide%3B%20the%20power-law%20tilts%20drive%20the%20secondary%20valleys%20toward%20the%0A%20%20%20%20(source-like)%20solar%20pattern%20as%20the%20energy%20increases%2C%20until%20the%20pattern%0A%20%20%20%20freezes%20at%20the%20saturation%20rigidity%20%24R_%5Cmathrm%7Bsat%7D%3D5%24%20PV.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(LODDERS_2003%2C%20RC_SANS%2C%20SYMBOLS%2C%20gsf_kn%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20runs%2F2026-07-26_boron-toa-refit%2Finputs%2Fplot_abundance_vs_z.py%0A%20%20%20%20with%20plt.rc_context(RC_SANS)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(5.6%2C%203.4))%0A%20%20%20%20%20%20%20%20_zs%20%3D%20np.arange(1%2C%2029)%0A%20%20%20%20%20%20%20%20_ek%20%3D%20%5B10.0**_k%20for%20_k%20in%20range(10)%5D%0A%20%20%20%20%20%20%20%20_colors%20%3D%20plt.cm.viridis(np.linspace(0.88%2C%200.02%2C%20len(_ek)))%0A%0A%20%20%20%20%20%20%20%20def%20_elabel(e)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20e%20%3C%201e4%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20return%20f%22%7Be%3Ag%7D%20GeV%2Fn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20rf%22%2410%5E%7B%7B%7Bint(np.log10(e))%7D%7D%7D%24%20GeV%2Fn%22%0A%0A%20%20%20%20%20%20%20%20for%20_e%2C%20_col%20in%20zip(_ek%2C%20_colors%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_f%20%3D%20np.array(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20float(np.asarray(gsf_kn.flux(np.array(%5B_e%5D)%2C%20int(_z))).ravel()%5B0%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20_z%20in%20_zs%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%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%20_si%20%3D%20_f%5B13%5D%20%20%23%20Z%20%3D%2014%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_si%20%3C%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20_rel%20%3D%20100.0%20*%20_f%20%2F%20_si%0A%20%20%20%20%20%20%20%20%20%20%20%20_ok%20%3D%20_rel%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_zs%5B_ok%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_rel%5B_ok%5D%2C%0A%20%20%20%20%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%20%20%20%20%20color%3D_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D_elabel(_e)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20_solar%20%3D%20np.asarray(LODDERS_2003%2C%20float)%0A%20%20%20%20%20%20%20%20_solar%20%3D%20_solar%20*%20(100.0%20%2F%20_solar%5B13%5D)%0A%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20_zs%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_solar%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22--o%22%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%3D1.1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ms%3D2.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mfc%3D%22white%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mew%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22Solar%20system%20(Lodders%202003)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D2%2C%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20_ax.set_yscale(%22log%22)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(0.5%2C%2028.8)%0A%20%20%20%20%20%20%20%20_ax.set_xticks(_zs)%0A%20%20%20%20%20%20%20%20_ax.set_xticklabels(%5Bstr(_z)%20for%20_z%20in%20_zs%5D%2C%20fontsize%3D6.5)%0A%20%20%20%20%20%20%20%20_sec%20%3D%20_ax.secondary_xaxis(%22top%22)%0A%20%20%20%20%20%20%20%20_sec.set_xticks(_zs)%0A%20%20%20%20%20%20%20%20_sec.set_xticklabels(SYMBOLS%2C%20fontsize%3D6.5)%0A%20%20%20%20%20%20%20%20_sec.tick_params(length%3D0%2C%20pad%3D2)%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(%22nuclear%20charge%20%24Z%24%22%2C%20fontsize%3D9)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(%22relative%20abundance%20(Si%20%24%5C%5Cequiv%24%20100)%22%2C%20fontsize%3D9)%0A%20%20%20%20%20%20%20%20_ax.tick_params(axis%3D%22y%22%2C%20labelsize%3D7)%0A%20%20%20%20%20%20%20%20_ax.grid(True%2C%20which%3D%22major%22%2C%20color%3D%22%23e3e5e8%22%2C%20lw%3D0.5)%0A%20%20%20%20%20%20%20%20_ax.set_axisbelow(True)%0A%20%20%20%20%20%20%20%20_ax.spines%5B%22top%22%5D.set_visible(False)%0A%20%20%20%20%20%20%20%20_ax.spines%5B%22right%22%5D.set_visible(True)%0A%20%20%20%20%20%20%20%20_ax.tick_params(axis%3D%22y%22%2C%20which%3D%22both%22%2C%20right%3DTrue%2C%20direction%3D%22in%22)%0A%20%20%20%20%20%20%20%20_ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D7%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ncol%3D3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22upper%20right%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20title%3D%22kinetic%20energy%20per%20nucleon%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20title_fontsize%3D7%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20columnspacing%3D1.2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D1.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20borderaxespad%3D0.3%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_fig.tight_layout(pad%3D0.4)%0A%20%20%20%20show(_fig)%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%20SM%20S3%20%E2%80%94%20Sub-leading%20to%20leader%20flux%20ratios%0A%0A%20%20%20%20One%20panel%20per%20sub-leading%20species%3A%20the%20fitted%20spline%20ratio%20to%20its%20group%0A%20%20%20%20leader%20(black)%20over%20the%20measured%20range%2C%20the%20slope-fit%20window%20(shaded)%2C%20the%0A%20%20%20%20adopted%20power-law%20extrapolation%20above%20the%20matching%20point%20with%20the%20fitted%0A%20%20%20%20index%20%24s_j%24%20quoted%20per%20panel%20(blue%2C%20saturating%20to%20a%20constant%20ratio%20at%0A%20%20%20%20%24R_%5Cmathrm%7Bsat%7D%3D5%24%20PV)%2C%20and%20the%20constant-ratio%20extrapolation%20of%20previous%0A%20%20%20%20GSF%20releases%20(grey%20dashed).%20The%20paper%20adds%20the%20offset-%20and%0A%20%20%20%20demodulation-corrected%20data%20ratios%20that%20enter%20the%20likelihood.%0A%0A%20%20%20%20The%20window%20is%20the%20last%20decade%20of%20the%20species'%20direct%20data%20below%20the%0A%20%20%20%20matching%20point%3B%20here%20it%20is%20taken%20as%20the%20decade%20below%20the%20top%20knot%2C%20which%0A%20%20%20%20the%20fit%20ties%20to%20the%20last%20datum%20(within%20a%20few%20percent%20of%20the%20fitted%20window%0A%20%20%20%20for%20every%20species).%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(RC_SANS%2C%20gsf_r%2C%20np%2C%20plt%2C%20show%2C%20species_name)%3A%0A%20%20%20%20%23%20runs%2F2026-07-26_boron-toa-refit%2Finputs%2Fcapture_subleading.py%3A%20render()%0A%20%20%20%20import%20matplotlib.ticker%20as%20_mticker%0A%20%20%20%20from%20matplotlib.lines%20import%20Line2D%20as%20_Line2D3%0A%20%20%20%20from%20matplotlib.patches%20import%20Patch%20as%20_Patch3%0A%0A%20%20%20%20_C_SPLINE%2C%20_C_NEW%2C%20_C_OLD%20%3D%20%22%2330363d%22%2C%20%22%232f6fde%22%2C%20%22%239aa0a6%22%0A%20%20%20%20_X_END%2C%20_X_SAT%20%3D%20np.log(1e8)%2C%20np.log(5.0e6)%0A%0A%20%20%20%20_subs%20%3D%20%5B%5D%0A%20%20%20%20for%20_sid%20in%20gsf_r.species%3A%0A%20%20%20%20%20%20%20%20_lsid%2C%20_norm%20%3D%20gsf_r.flux_ratio%5B_sid%5D%0A%20%20%20%20%20%20%20%20if%20_sid%20!%3D%20_lsid%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_subs.append((_sid%2C%20_lsid%2C%20_norm%2C%20gsf_r.flux_slope%5B_sid%5D))%0A%20%20%20%20_subs.sort(key%3Dlambda%20t%3A%20(t%5B1%5D%5B0%5D%2C%20t%5B0%5D%5B0%5D%2C%20t%5B0%5D%5B1%5D))%0A%0A%20%20%20%20with%20plt.rc_context(RC_SANS)%3A%0A%20%20%20%20%20%20%20%20_ncol%2C%20_nrow%20%3D%204%2C%207%0A%20%20%20%20%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(_nrow%2C%20_ncol%2C%20figsize%3D(7.08%2C%208.15)%2C%20sharex%3DTrue)%0A%20%20%20%20%20%20%20%20for%20_k%2C%20(_sid%2C%20_lsid%2C%20_norm%2C%20_slope)%20in%20enumerate(_subs)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%20%3D%20_axes%5B_k%20%2F%2F%20_ncol%2C%20_k%20%25%20_ncol%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_kx%20%3D%20np.asarray(gsf_r.kx%5B_sid%5D%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20_xmin%2C%20_xmax%20%3D%20_kx%5B0%5D%2C%20_kx%5B-1%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_x0%20%3D%20max(_xmin%2C%20float(np.asarray(gsf_r.kx%5B_lsid%5D%2C%20float)%5B0%5D))%20%2B%200.02%0A%20%20%20%20%20%20%20%20%20%20%20%20_x_lo%20%3D%20np.linspace(_x0%2C%20_xmax%2C%20200)%0A%20%20%20%20%20%20%20%20%20%20%20%20_R_lo%20%3D%20np.exp(_x_lo)%0A%20%20%20%20%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%20%20%20%20_r_lo%20%3D%20np.asarray(gsf_r.flux(_R_lo%2C%20_sid)%2C%20float)%20%2F%20np.asarray(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20gsf_r.flux(_R_lo%2C%20_lsid)%2C%20float%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_r_lo%5B~np.isfinite(_r_lo)%5D%20%3D%200.0%0A%20%20%20%20%20%20%20%20%20%20%20%20_x_hi%20%3D%20np.linspace(_xmax%2C%20_X_END%2C%20160)%0A%20%20%20%20%20%20%20%20%20%20%20%20_r_new%20%3D%20_norm%20*%20np.exp(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_slope%20*%20np.clip(_x_hi%20-%20_xmax%2C%200.0%2C%20max(_X_SAT%20-%20_xmax%2C%200.0))%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%20_ax.axvspan(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.exp(_xmax)%20%2F%2010.0%2C%20np.exp(_xmax)%2C%20color%3D_C_NEW%2C%20alpha%3D0.10%2C%20lw%3D0%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%20_ok%20%3D%20_r_lo%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(_R_lo%5B_ok%5D%2C%20_r_lo%5B_ok%5D%2C%20color%3D_C_SPLINE%2C%20lw%3D1.1%2C%20zorder%3D3)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.exp(_x_hi)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.full_like(_x_hi%2C%20_norm)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_C_OLD%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0.9%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3D%22--%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D2%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%20(_ln_new%2C)%20%3D%20_ax.plot(np.exp(_x_hi)%2C%20_r_new%2C%20color%3D_C_NEW%2C%20lw%3D1.4%2C%20zorder%3D4)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xscale(%22log%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_yscale(%22log%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ymax%20%3D%20max(np.max(_r_lo%5B_ok%5D)%20if%20np.any(_ok)%20else%20_norm%2C%20_norm)%20*%203%0A%20%20%20%20%20%20%20%20%20%20%20%20_ymin%20%3D%20min(np.min(_r_new)%2C%20_norm)%20%2F%205%0A%20%20%20%20%20%20%20%20%20%20%20%20_ymin%20%3D%20max(_ymin%2C%20_ymax%20*%201e-6)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_ymax%20%2F%20_ymin%20%3C%2010**1.5%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_c%20%3D%20np.sqrt(_ymax%20*%20_ymin)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ymin%2C%20_ymax%20%3D%20_c%20%2F%2010**0.75%2C%20_c%20*%2010**0.75%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_ylim(_ymin%2C%20_ymax)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.yaxis.set_major_locator(_mticker.LogLocator(base%3D10.0%2C%20numticks%3D4))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.yaxis.set_major_formatter(_mticker.LogFormatterMathtext(base%3D10.0))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.yaxis.set_minor_locator(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_mticker.LogLocator(base%3D10.0%2C%20subs%3Dnp.arange(2%2C%2010)%20*%200.1%2C%20numticks%3D99)%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%20_ax.yaxis.set_minor_formatter(_mticker.NullFormatter())%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B_ln_new%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5Bf%22%7Bspecies_name(_sid)%7D%2F%7Bspecies_name(_lsid)%7D%20%20s%3D%7B_slope%3A%2B.3f%7D%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22lower%20left%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D6.3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D1.2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20borderaxespad%3D0.15%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%20_ax.tick_params(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20which%3D%22both%22%2C%20direction%3D%22in%22%2C%20top%3DTrue%2C%20right%3DTrue%2C%20labelsize%3D6%2C%20pad%3D2%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%20_ax.grid(True%2C%20axis%3D%22x%22%2C%20which%3D%22major%22%2C%20color%3D%22%23e3e5e8%22%2C%20lw%3D0.5)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_axisbelow(True)%0A%20%20%20%20%20%20%20%20_axes%5B0%2C%200%5D.set_xlim(0.8%2C%201.3e8)%0A%20%20%20%20%20%20%20%20for%20_ax%20in%20_axes.ravel()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_major_locator(_mticker.LogLocator(base%3D100.0))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_major_formatter(_mticker.LogFormatterMathtext(base%3D10.0))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_minor_locator(_mticker.LogLocator(base%3D10.0%2C%20numticks%3D99))%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.xaxis.set_minor_formatter(_mticker.NullFormatter())%0A%20%20%20%20%20%20%20%20_leg_ax%20%3D%20_axes%5B6%2C%201%5D%0A%20%20%20%20%20%20%20%20_leg_ax.axis(%22off%22)%0A%20%20%20%20%20%20%20%20_leg_ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20handles%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_Line2D3(%5B%5D%2C%20%5B%5D%2C%20color%3D_C_SPLINE%2C%20lw%3D1.1%2C%20label%3D%22fitted%20spline%20ratio%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_Patch3(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20facecolor%3D_C_NEW%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.10%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22slope-fit%20window%20(last%20data%20decade)%22%2C%0A%20%20%20%20%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%20%20%20%20%20_Line2D3(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_C_NEW%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Dr%22extrapolation%20norm%24%5Ccdot(%5Cmin(R%2CR_%7B%5Crm%20sat%7D)%2F%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20r%22R_%7B%5Crm%20max%7D)%5E%7Bs%7D%24%2C%20%24R_%7B%5Crm%20sat%7D%3D5%24%20PV%22%2C%0A%20%20%20%20%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%20%20%20%20%20_Line2D3(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%5D%2C%20%5B%5D%2C%20color%3D_C_OLD%2C%20lw%3D0.9%2C%20ls%3D%22--%22%2C%20label%3D%22constant%20ratio%20(old)%22%0A%20%20%20%20%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%20%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22center%20left%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20bbox_to_anchor%3D(0.02%2C%200.5)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D6.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D1.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20borderaxespad%3D0.0%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20for%20_k%20in%20range(len(_subs)%2C%20_nrow%20*%20_ncol)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_k%20!%3D%2025%3A%20%20%23%2025%20%3D%20the%20legend%20panel%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B_k%20%2F%2F%20_ncol%2C%20_k%20%25%20_ncol%5D.axis(%22off%22)%0A%20%20%20%20%20%20%20%20for%20_j%20in%20range(_ncol)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax%20%3D%20_axes%5B6%2C%20_j%5D%20if%20_axes%5B6%2C%20_j%5D.axison%20else%20_axes%5B5%2C%20_j%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.tick_params(labelbottom%3DTrue)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xlabel(r%22%24R%24%20%5BGV%5D%22%2C%20fontsize%3D7%2C%20labelpad%3D1.5)%0A%20%20%20%20%20%20%20%20for%20_i%20in%20range(_nrow)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_axes%5B_i%2C%200%5D.axison%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B_i%2C%200%5D.set_ylabel(%22flux%20ratio%22%2C%20fontsize%3D7%2C%20labelpad%3D1.5)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(%0A%20%20%20%20%20%20%20%20%20%20%20%20left%3D0.075%2C%20right%3D0.982%2C%20top%3D0.995%2C%20bottom%3D0.05%2C%20wspace%3D0.28%2C%20hspace%3D0.10%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20show(_fig)%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%20SM%20S4%20%E2%80%94%20Impact%20of%20the%20solar-modulation%20potential%0A%0A%20%20%20%20Fitted%20LIS%20of%20the%20four%20mass%20groups%20relative%20to%20the%20default%0A%20%20%20%20Ghelfi%E2%80%93Maurin%E2%80%93Derome%20reconstruction%20(parameter%20set%20%602026.0%60)%3B%20the%20grey%20band%0A%20%20%20%20is%20the%20%241%5Csigma%24%20uncertainty%20of%20the%20default.%20Blue%3A%20the%20Usoskin%202017%0A%20%20%20%20potential%20with%20the%20same%20window%20average%20and%20%24%5Cdelta%5Cphi%24%20nuisances%0A%20%20%20%20(parameter%20set%20%602026.0-USO%60).%0A%0A%20%20%20%20The%20paper's%20third%20curve%20%E2%80%94%20the%20single-epoch%20demodulation%20used%20by%20previous%0A%20%20%20%20GSF%20releases%20%E2%80%94%20comes%20from%20a%20third%20re-fit%20that%20is%20not%20part%20of%20the%20released%0A%20%20%20%20model%2C%20so%20it%20cannot%20be%20drawn%20here.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(GROUPS%2C%20GSFRigidity%2C%20RC_SERIF%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20runs%2F2026-07-26_boron-toa-refit%2Finputs%2Ffig_phi_source.py%0A%20%20%20%20import%20matplotlib.ticker%20as%20_mticker2%0A%0A%20%20%20%20_C_USO%20%3D%20%22%231f77b4%22%0A%0A%20%20%20%20with%20plt.rc_context(%7B**RC_SERIF%2C%20%22font.size%22%3A%208%7D)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(4%2C%201%2C%20figsize%3D(3.4%2C%205.7)%2C%20sharex%3DTrue%2C%20sharey%3DTrue)%0A%20%20%20%20%20%20%20%20_R%20%3D%20np.logspace(np.log10(1.3)%2C%20np.log10(50)%2C%20160)%0A%20%20%20%20%20%20%20%20_gmd%20%3D%20GSFRigidity(version%3D%222026.0%22)%0A%20%20%20%20%20%20%20%20_uso%20%3D%20GSFRigidity(version%3D%222026.0-USO%22)%0A%20%20%20%20%20%20%20%20_panel_labels%20%3D%20(%22H*%22%2C%20%22He*%22%2C%20%22O*%22%2C%20%22Fe*%22)%0A%20%20%20%20%20%20%20%20for%20_i%2C%20(_ax%2C%20_g)%20in%20enumerate(zip(_axes%2C%20GROUPS%2C%20strict%3DFalse))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_jg%20%3D%20np.asarray(_gmd.flux(_R%2C%20_g%2C%20time_interval%3D%22LIS%22)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20_eg%20%3D%20np.asarray(_gmd.error(_R%2C%20_g%2C%20time_interval%3D%22LIS%22)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ref%20%3D%20np.where(_jg%20%3E%200%2C%20_jg%2C%201.0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.fill_between(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_R%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20-100%20*%20_eg%20%2F%20_ref%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20100%20*%20_eg%20%2F%20_ref%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D%220.55%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.55%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Dr%22GMD%20(default)%20%24%5Cpm1%5Csigma%24%22%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%20_ju%20%3D%20np.asarray(_uso.flux(_R%2C%20_g%2C%20time_interval%3D%22LIS%22)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_R%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20100%20*%20(_ju%20%2F%20_ref%20-%201.0)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_C_USO%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22Usoskin%20(USO)%22%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%20_ax.axhline(0%2C%20color%3D%22k%22%2C%20lw%3D0.5)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.set_xscale(%22log%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.text(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.035%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.88%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_panel_labels%5B_i%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20transform%3D_ax.transAxes%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D9%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%20va%3D%22top%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22left%22%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%20_ax.grid(alpha%3D0.3)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.tick_params(labelsize%3D7.5)%0A%20%20%20%20%20%20%20%20_axes%5B-1%5D.set_xticks(%5B2%2C%205%2C%2010%2C%2020%2C%2050%5D)%0A%20%20%20%20%20%20%20%20_axes%5B-1%5D.xaxis.set_major_formatter(_mticker2.ScalarFormatter())%0A%20%20%20%20%20%20%20%20_axes%5B-1%5D.xaxis.set_minor_formatter(_mticker2.NullFormatter())%0A%20%20%20%20%20%20%20%20_axes%5B-1%5D.set_xlabel(r%22rigidity%20%24R%24%20%5BGV%5D%22)%0A%20%20%20%20%20%20%20%20_axes%5B0%5D.set_ylim(-22%2C%2034)%0A%20%20%20%20%20%20%20%20_axes%5B0%5D.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20frameon%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D6.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22upper%20right%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20handlelength%3D1.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20borderaxespad%3D0.4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20labelspacing%3D0.3%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_fig.supylabel(%22deviation%20from%20GMD%20default%20%20%5B%25%5D%22%2C%20fontsize%3D8%2C%20x%3D0.02)%0A%20%20%20%20%20%20%20%20_fig.tight_layout(pad%3D0.3%2C%20h_pad%3D0.4)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.155)%0A%20%20%20%20show(_fig)%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%20SM%20S7%20%E2%80%94%20Correlation%20matrix%20of%20the%20fitted%20spline%20amplitudes%0A%0A%20%20%20%20Correlation%20of%20the%20leading%20spline%20amplitudes%2C%20grouped%20into%20the%20%24p%24%2C%20He%2C%0A%20%20%20%20O%5C*%20and%20Fe%5C*%20blocks.%20The%20energy-scale%20offsets%20are%20not%20part%20of%20the%20published%0A%20%20%20%20covariance%20and%20are%20not%20shown.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(ELEM_LABELS%2C%20FIGSIZE_SMALL%2C%20RC_SERIF%2C%20gsf_e%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20deck%2Fmisc.py%3A%20covariance_plot%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_leaders%20%3D%20%5Bgsf_e.z_to_sids%5B_z%5D%5B0%5D%20for%20_z%20in%20(1%2C%202%2C%208%2C%2026)%5D%0A%20%20%20%20%20%20%20%20_sizes%20%3D%20%5Bnp.shape(gsf_e.cov%5B(_l%2C%20_l)%5D)%5B0%5D%20for%20_l%20in%20_leaders%5D%0A%20%20%20%20%20%20%20%20_n1%20%3D%20sum(_sizes)%0A%20%20%20%20%20%20%20%20_cov%20%3D%20np.zeros((_n1%2C%20_n1))%0A%20%20%20%20%20%20%20%20_off%20%3D%20np.cumsum(%5B0%2C%20*_sizes%5D)%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_li%20in%20enumerate(_leaders)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_j%2C%20_lj%20in%20enumerate(_leaders)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_block%20%3D%20gsf_e.cov.get((_li%2C%20_lj))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20_block%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cov%5B_off%5B_i%5D%20%3A%20_off%5B_i%20%2B%201%5D%2C%20_off%5B_j%5D%20%3A%20_off%5B_j%20%2B%201%5D%5D%20%3D%20np.asarray(%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_block%2C%20float%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_m%20%3D%20np.zeros((_n1%2C%20_n1))%0A%20%20%20%20%20%20%20%20for%20_i%20in%20range(_n1)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_j%20in%20range(_i%20%2B%201%2C%20_n1)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_m%5B_i%2C%20_j%5D%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cov%5B_i%2C%20_j%5D%20%2F%20abs(_cov%5B_i%2C%20_i%5D%20*%20_cov%5B_j%2C%20_j%5D%20%2B%201e-300)%20**%200.5%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_m%5B_j%2C%20_i%5D%20%3D%20_m%5B_i%2C%20_j%5D%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3DFIGSIZE_SMALL)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(left%3D0.12%2C%20right%3D0.95%2C%20top%3D0.95%2C%20bottom%3D0.12)%0A%20%20%20%20%20%20%20%20_pc%20%3D%20_ax.pcolormesh(_m.T%2C%20cmap%3D%22bwr%22%2C%20vmin%3D-1%2C%20vmax%3D1)%0A%20%20%20%20%20%20%20%20_fig.colorbar(_pc%2C%20ax%3D_ax)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(0%2C%20_n1)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(0%2C%20_n1)%0A%20%20%20%20%20%20%20%20for%20_i%20in%20range(len(_leaders))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ia%2C%20_ib%20%3D%20_off%5B_i%5D%2C%20_off%5B_i%20%2B%201%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.plot(%5B_ia%2C%20_ib%2C%20_ib%2C%20_ia%2C%20_ia%5D%2C%20%5B_ia%2C%20_ia%2C%20_ib%2C%20_ib%2C%20_ia%5D%2C%20%22-k%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.text(0.5%20*%20(_ia%20%2B%20_ib)%2C%20-2.5%2C%20ELEM_LABELS%5B_i%5D%2C%20va%3D%22top%22%2C%20ha%3D%22center%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.text(-2.5%2C%200.5%20*%20(_ia%20%2B%20_ib)%2C%20ELEM_LABELS%5B_i%5D%2C%20va%3D%22center%22%2C%20ha%3D%22right%22)%0A%20%20%20%20%20%20%20%20_ax.axis(%22off%22)%0A%20%20%20%20show(_fig)%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%20SM%20S8%20%E2%80%94%20Pivot%20components%20against%20the%20full%20covariance%0A%0A%20%20%20%20Relative%20%241%5Csigma%24%20uncertainty%20of%20the%20total%20proton%20and%20neutron%20fluxes%0A%20%20%20%20from%20the%20full%20amplitude%20covariance%20(solid)%20and%20from%20the%2024%20pivot%0A%20%20%20%20components%20(dashed)%3B%20dots%20mark%20the%20pivot%20energies.%20Bottom%3A%20ratio%20of%20the%0A%20%20%20%20two%20%E2%80%94%20within%20the%20shaded%20band%20the%20reduced%20uncertainty%20stays%20inside%20the%0A%20%20%20%20worst-case%20coverage%20factor%2C%20and%20it%20is%20exact%20at%20the%20pivots.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(RC_SANS%2C%20RC_SMALL%2C%20gsf_en%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20runs%2F2026-08-05_reduced-fidelity-figure%2Finputs%2Fmake_reduced_fidelity_figure.py%0A%20%20%20%20from%20globalsplinefit.reduced%20import%20ReducedGSF%2C%20_relative_species_system%0A%0A%20%20%20%20_COL%20%3D%20%7B%22p%22%3A%20%22%230072B2%22%2C%20%22n%22%3A%20%22%23D55E00%22%7D%20%20%23%20Okabe-Ito%20blue%20%2F%20vermillion%0A%20%20%20%20_red%20%3D%20ReducedGSF(gsf_en)%0A%20%20%20%20_lo%2C%20_hi%20%3D%20_red.pivot_energies%5B0%5D%2C%20_red.pivot_energies%5B-1%5D%0A%20%20%20%20_e%20%3D%20np.logspace(np.log10(_lo)%2C%20np.log10(_hi)%2C%20601)%0A%20%20%20%20%23%20exact%20reference%20from%20the%20construction's%20own%20trimmed%20stacked%20system%2C%20so%0A%20%20%20%20%23%20the%20ratio%20panel%20reproduces%20the%20coverage-factor%20definition%20verbatim%0A%20%20%20%20_jac_rel%2C%20_cov_par%2C%20_central%2C%20_species%20%3D%20_relative_species_system(%0A%20%20%20%20%20%20%20%20gsf_en%2C%20_e%2C%20per_group%3DFalse%0A%20%20%20%20)%0A%20%20%20%20_exact%20%3D%20np.sqrt(np.einsum(%22ij%2Cjk%2Cik-%3Ei%22%2C%20_jac_rel%2C%20_cov_par%2C%20_jac_rel)).reshape(%0A%20%20%20%20%20%20%20%20len(_species)%2C%20len(_e)%0A%20%20%20%20)%0A%20%20%20%20_reduced%20%3D%20_red.error(_e)%20%2F%20_red.flux(_e)%0A%20%20%20%20_ratio%20%3D%20_reduced%20%2F%20_exact%0A%20%20%20%20_worst%20%3D%20float(np.exp(np.max(np.abs(np.log(_ratio)))))%0A%20%20%20%20_piv%20%3D%20_red.error(_red.pivot_energies)%20%2F%20_red.flux(_red.pivot_energies)%0A%0A%20%20%20%20with%20plt.rc_context(%7B**RC_SANS%2C%20**RC_SMALL%7D)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20(_ax1%2C%20_ax2)%20%3D%20plt.subplots(%0A%20%20%20%20%20%20%20%20%20%20%20%202%2C%0A%20%20%20%20%20%20%20%20%20%20%20%201%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20figsize%3D(3.4%2C%203.6)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20sharex%3DTrue%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20gridspec_kw%3D%7B%22height_ratios%22%3A%20%5B2.0%2C%201.0%5D%2C%20%22hspace%22%3A%200.06%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20for%20_s%2C%20_name%20in%20enumerate(_species)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax1.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%20100%20*%20_exact%5B_s%5D%2C%20color%3D_COL%5B_name%5D%2C%20lw%3D1.4%2C%20label%3Df%22%7B_name%7D%2C%20full%22%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%20_ax1.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20100%20*%20_reduced%5B_s%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_COL%5B_name%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3D%22--%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22%7B_name%7D%2C%20reduced%22%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%20_ax1.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_red.pivot_energies%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20100%20*%20_piv%5B_s%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22o%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_COL%5B_name%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ms%3D3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mec%3D%22none%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax1.set_ylabel(r%22relative%20%241%5Csigma%24%20uncertainty%20(%25)%22)%0A%20%20%20%20%20%20%20%20_ax1.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22upper%20left%22%2C%20frameon%3DFalse%2C%20ncol%3D2%2C%20columnspacing%3D1.2%2C%20handlelength%3D1.8%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax1.tick_params(which%3D%22both%22%2C%20direction%3D%22in%22%2C%20top%3DTrue%2C%20right%3DTrue)%0A%20%20%20%20%20%20%20%20_ax2.axhspan(1%20%2F%20_worst%2C%20_worst%2C%20color%3D%220.90%22%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20_ax2.axhline(1.0%2C%20color%3D%22k%22%2C%20lw%3D0.7)%0A%20%20%20%20%20%20%20%20for%20_s%2C%20_name%20in%20enumerate(_species)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax2.semilogx(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_e%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ratio%5B_s%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_COL%5B_name%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3D%22-%22%20if%20_name%20%3D%3D%20%22p%22%20else%20%22--%22%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%20_ax2.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_red.pivot_energies%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20np.ones_like(_red.pivot_energies)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22o%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_COL%5B_name%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ms%3D3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mec%3D%22none%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax2.set_xlabel(%22kinetic%20energy%20per%20nucleon%20(GeV)%22)%0A%20%20%20%20%20%20%20%20_ax2.set_ylabel(r%22%24%5Csigma_%7B%5Crm%20red%7D%2F%5Csigma_%7B%5Crm%20full%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax2.set_ylim(0.65%2C%201.45)%0A%20%20%20%20%20%20%20%20_ax2.set_xlim(_lo%2C%20_hi)%0A%20%20%20%20%20%20%20%20_ax2.tick_params(which%3D%22both%22%2C%20direction%3D%22in%22%2C%20top%3DTrue%2C%20right%3DTrue)%0A%20%20%20%20show(_fig)%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%20SM%20S9%20%E2%80%94%20Sampling%20check%20for%20pseudo-experiments%0A%0A%20%20%20%20Empirical%20standard%20deviation%20of%202000%20draws%20from%20the%20native%20amplitude%0A%20%20%20%20covariance%2C%20divided%20by%20the%20exact%20Jacobian-propagated%20uncertainty%2C%20for%20the%0A%20%20%20%20four%20mass%20groups%20and%20the%20all-particle%20flux%2C%20in%20total%20energy%20per%20nucleus.%0A%20%20%20%20The%20residual%20scatter%20is%20the%20statistical%20fluctuation%20of%20the%20ensemble.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(GROUPS%2C%20RC_SANS%2C%20RC_SMALL%2C%20gsf_e%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20%23%20runs%2F2026-08-06_reduced-sampling-nucleus%2Finputs%2Fmake_native_sampling_figure.py%0A%20%20%20%20_COLG%20%3D%20%7B%22H*%22%3A%20%22%23D55E00%22%2C%20%22He%22%3A%20%22%23E69F00%22%2C%20%22O*%22%3A%20%22%23009E73%22%2C%20%22Fe*%22%3A%20%22%230072B2%22%7D%0A%20%20%20%20_E_LO%2C%20_E_HI%2C%20_N_GRID%2C%20_N_DRAWS%20%3D%201e2%2C%201e11%2C%20601%2C%202000%0A%0A%20%20%20%20_energy%20%3D%20np.logspace(np.log10(_E_LO)%2C%20np.log10(_E_HI)%2C%20_N_GRID)%0A%20%20%20%20_central%20%3D%20np.array(%5Bnp.asarray(gsf_e.flux(_energy%2C%20_g)%2C%20float)%20for%20_g%20in%20GROUPS%5D)%0A%20%20%20%20_jacs%20%3D%20%5Bnp.asarray(gsf_e.jacobian(_energy%2C%20_g)%2C%20float)%20for%20_g%20in%20GROUPS%5D%0A%20%20%20%20_widths%20%3D%20%5B_j.shape%5B1%5D%20for%20_j%20in%20_jacs%5D%0A%20%20%20%20_sids%20%3D%20%5Bgsf_e.z_to_sids%5B_z%5D%5B0%5D%20for%20_z%20in%20(1%2C%202%2C%208%2C%2026)%5D%0A%0A%20%20%20%20_P%20%3D%20sum(_widths)%0A%20%20%20%20_C%20%3D%20np.zeros((_P%2C%20_P))%0A%20%20%20%20_off%20%3D%20np.cumsum(%5B0%2C%20*_widths%5D)%0A%20%20%20%20for%20_i%2C%20_si%20in%20enumerate(_sids)%3A%0A%20%20%20%20%20%20%20%20for%20_j%2C%20_sj%20in%20enumerate(_sids)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_block%20%3D%20gsf_e.cov.get((_si%2C%20_sj))%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_block%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_C%5B_off%5B_i%5D%20%3A%20_off%5B_i%20%2B%201%5D%2C%20_off%5B_j%5D%20%3A%20_off%5B_j%20%2B%201%5D%5D%20%3D%20_block%0A%20%20%20%20_C%20%3D%200.5%20*%20(_C%20%2B%20_C.T)%0A%20%20%20%20_w%2C%20_V%20%3D%20np.linalg.eigh(_C)%0A%20%20%20%20_factor%20%3D%20_V%20*%20np.sqrt(np.clip(_w%2C%200.0%2C%20None))%0A%0A%20%20%20%20_sig_g%20%3D%20np.array(%5Bnp.asarray(gsf_e.error(_energy%2C%20_g)%2C%20float)%20for%20_g%20in%20GROUPS%5D)%0A%20%20%20%20_var_all%20%3D%20np.zeros(_N_GRID)%0A%20%20%20%20for%20_gi%20in%20GROUPS%3A%0A%20%20%20%20%20%20%20%20for%20_gj%20in%20GROUPS%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_var_all%20%2B%3D%20np.diag(np.asarray(gsf_e.covariance(_gi%2C%20_gj%2C%20_energy)%2C%20float))%0A%20%20%20%20_sig_all%20%3D%20np.sqrt(_var_all)%0A%0A%20%20%20%20_rng%20%3D%20np.random.default_rng(0)%0A%20%20%20%20_dc%20%3D%20_rng.standard_normal((_N_DRAWS%2C%20_P))%20%40%20_factor.T%0A%20%20%20%20_draws%20%3D%20np.empty((_N_DRAWS%2C%20len(GROUPS)%2C%20_N_GRID))%0A%20%20%20%20for%20_i%20in%20range(len(GROUPS))%3A%0A%20%20%20%20%20%20%20%20_draws%5B%3A%2C%20_i%5D%20%3D%20_central%5B_i%5D%20%2B%20_dc%5B%3A%2C%20_off%5B_i%5D%20%3A%20_off%5B_i%20%2B%201%5D%5D%20%40%20_jacs%5B_i%5D.T%0A%20%20%20%20_r_g%20%3D%20_draws.std(axis%3D0)%20%2F%20_sig_g%0A%20%20%20%20_r_all%20%3D%20_draws.sum(axis%3D1).std(axis%3D0)%20%2F%20_sig_all%0A%0A%20%20%20%20with%20plt.rc_context(%7B**RC_SANS%2C%20**RC_SMALL%7D)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(3.4%2C%202.3))%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_g%20in%20enumerate(GROUPS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.semilogx(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_energy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_r_g%5B_i%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_COLG%5B_g%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22H%22%20if%20_g%20%3D%3D%20%22H*%22%20else%20_g%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.semilogx(_energy%2C%20_r_all%2C%20color%3D%22k%22%2C%20lw%3D1.4%2C%20label%3D%22all-particle%22)%0A%20%20%20%20%20%20%20%20_ax.axhline(1%2C%20color%3D%22k%22%2C%20lw%3D0.7)%0A%20%20%20%20%20%20%20%20_ax.set_xlabel(%22total%20energy%20per%20nucleus%20(GeV)%22)%0A%20%20%20%20%20%20%20%20_ax.set_ylabel(r%22%24%5Csigma_%7B%5Crm%20sample%7D%2F%5Csigma_%7B%5Crm%20exact%7D%24%22)%0A%20%20%20%20%20%20%20%20_ax.set_xlim(_E_LO%2C%20_E_HI)%0A%20%20%20%20%20%20%20%20_ax.set_ylim(0.85%2C%201.15)%0A%20%20%20%20%20%20%20%20_ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22upper%20left%22%2C%20frameon%3DFalse%2C%20ncol%3D3%2C%20columnspacing%3D1.0%2C%20handlelength%3D1.4%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.tick_params(which%3D%22both%22%2C%20direction%3D%22in%22%2C%20top%3DTrue%2C%20right%3DTrue)%0A%20%20%20%20show(_fig)%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%20SM%20S14%20%E2%80%94%20Light%20and%20heavy%20fluxes%0A%0A%20%20%20%20Light%20(%24p%2B%24He)%20and%20heavy%20(O%5C*%24%2B%24Fe%5C*)%20fluxes%2C%20weighted%20by%0A%20%20%20%20%24E_%5Cmathrm%7Bkin%7D%5E%7B3%7D%24%2C%20against%20the%20all-particle%20flux%20(black).%20In%20the%20paper%0A%20%20%20%20both%20panels%20carry%20the%20experiments%20that%20report%20two-component%20composition.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20FIGSIZE_WIDE%2C%0A%20%20%20%20HEAVY_COLOR%2C%0A%20%20%20%20LIGHT_COLOR%2C%0A%20%20%20%20RC_SERIF%2C%0A%20%20%20%20draw_band%2C%0A%20%20%20%20gsf_k%2C%0A%20%20%20%20logx_ticks%2C%0A%20%20%20%20np%2C%0A%20%20%20%20pair_error%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20%23%20deck%2Ffractions.py%3A%20light_heavy%0A%20%20%20%20from%20matplotlib.ticker%20import%20EngFormatter%20as%20_EngFormatter2%0A%0A%20%20%20%20with%20plt.rc_context(RC_SERIF)%3A%0A%20%20%20%20%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(2%2C%201%2C%20figsize%3DFIGSIZE_WIDE%2C%20sharex%3DTrue%2C%20sharey%3DTrue)%0A%20%20%20%20%20%20%20%20_fig.subplots_adjust(hspace%3D0%2C%20right%3D0.95%2C%20left%3D0.13%2C%20top%3D0.95%2C%20bottom%3D0.12)%0A%20%20%20%20%20%20%20%20_e%20%3D%20np.geomspace(1e2%2C%201e11%2C%201000)%0A%20%20%20%20%20%20%20%20_sc%20%3D%20_e**3.0%0A%20%20%20%20%20%20%20%20_tot%20%3D%20np.asarray(gsf_k.total_flux(_e)%2C%20float)%20*%20_sc%0A%20%20%20%20%20%20%20%20_pairs%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20((%22H*%22%2C%20%22He%22)%2C%20LIGHT_COLOR%2C%20r%22%24p%20%2B%20%5Cmathrm%7BHe%7D%24%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20((%22O*%22%2C%20%22Fe*%22)%2C%20HEAVY_COLOR%2C%20r%22%24%5Cmathrm%7BO*%7D%20%2B%20%5Cmathrm%7BFe*%7D%24%22)%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20for%20_i%2C%20(_pair%2C%20_col%2C%20_label)%20in%20enumerate(_pairs)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_val%20%3D%20sum(np.asarray(gsf_k.flux(_e%2C%20_g)%2C%20float)%20for%20_g%20in%20_pair)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20pair_error(gsf_k%2C%20_e%2C%20*_pair)%20*%20_sc%0A%20%20%20%20%20%20%20%20%20%20%20%20draw_band(_axes%5B_i%5D%2C%20_e%2C%20_val%2C%20_err%2C%20_col%2C%20line%3D%22-%22%2C%20lw%3D0.9%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B_i%5D.plot(_e%2C%20_tot%2C%20%22-%22%2C%20color%3D%22k%22%2C%20lw%3D1.5)%0A%20%20%20%20%20%20%20%20%20%20%20%20_axes%5B_i%5D.set_ylabel(_label)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.yaxis.set_major_formatter(_EngFormatter2())%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_xscale(%22log%22)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_xlabel(r%22%24E_%5Cmathrm%7Bkin%7D%20%2F%20%5Cmathrm%7BGeV%7D%24%22)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_yticks(np.arange(0%2C%204.5e6%2C%201e6))%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_ylim(0%2C%204.8e6)%0A%20%20%20%20%20%20%20%20_axes%5B1%5D.set_xlim(1e3%2C%201e11)%0A%20%20%20%20%20%20%20%20logx_ticks(_axes%5B1%5D)%0A%20%20%20%20%20%20%20%20_fig.text(%0A%20%20%20%20%20%20%20%20%20%20%20%200.010%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%20r%22%24J_%5Cmathrm%7BLIS%7D%20%2F%20(%5Cmathrm%7BGeV%7D%5C%2C%5Cmathrm%7Bm%7D%5E2%5C%2C%5Cmathrm%7Bs%7D%5C%2C%5Cmathrm%7Bsr%7D)%5E%7B-1%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20r%22%5Ctimes%20(E_%5Cmathrm%7Bkin%7D%2F%5Cmathrm%7BGeV%7D)%5E%7B3%7D%24%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20rotation%3D90%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)%0A%20%20%20%20show(_fig)%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