[#2026_RSI] Created

.gitattributes

@@ -0,0 +1,12 @@
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/d2cross_Ee0100eV-10MeV-80log_Eph100eV-10MeV-81log_ntheta61_EH.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/d2cross_Ee01eV-100MeV-80log_Eph1eV-100MeV-81log_ntheta61_EH.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/d2cross_phi_Ee01eV-100MeV-80log_Eph1eV-100MeV-81log_nthetaph61_nthetae060_EH.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/responsivities.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/d2cross_phi_Ee01eV-100MeV-240log_Eph1eV-100MeV-241log_nthetaph61_nthetae060_EH.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260616_SCRAM86_Ar_data.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260616_SCRAM86_Fe_data.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260616_SCRAM86_Kr_data.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260616_SCRAM86_O_data.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260616_SCRAM86_W_data.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260625_FLYCHK_H_data.npz filter=lfs diff=lfs merge=lfs -text
+articles/RSI_2026_RunawayBremsstrahlungDetection/inputs/260625_FLYCHK_He_data.npz filter=lfs diff=lfs merge=lfs -text

articles/RSI_2026_RunawayBremsstrahlungDetection/.gitattributes

@@ -0,0 +1,5 @@
+figures/fig01_cross.png filter=lfs diff=lfs merge=lfs -text
+figures/fig02_dist.png filter=lfs diff=lfs merge=lfs -text
+figures/fig04_bremsstrahlung.png filter=lfs diff=lfs merge=lfs -text
+figures/fig05_tokamak.png filter=lfs diff=lfs merge=lfs -text
+figures/fig06_responsivities.png filter=lfs diff=lfs merge=lfs -text

articles/RSI_2026_RunawayBremsstrahlungDetection/__init__.py

@@ -0,0 +1 @@
+from .code import *

articles/RSI_2026_RunawayBremsstrahlungDetection/code/__init__.py

@@ -0,0 +1,9 @@
+from ._main import main
+from ._load_spect import main as load_spect
+from ._fig01_cross import main as fig01_cross
+from ._fig02_dist import main as fig02_dist
+from ._fig03_validate_ff import main as fig03_validate_ff
+from ._fig04_bremsstrahlung import main as fig04_bremsstrahlung
+from ._fig05_emiss import main as fig05_emiss
+from ._fig06_tokamak import main as fig06_tokamak
+from ._fig07_responsivities import main as fig07_responsivities

articles/RSI_2026_RunawayBremsstrahlungDetection/code/_fig01_cross.py

@@ -0,0 +1,295 @@
+
+
+import os
+
+
+import numpy as np
+import matplotlib.pyplot as plt
+import matplotlib.gridspec as gridspec
+import tofu as tf
+
+
+from ._savefig import main as savefig
+
+
+tfphysemis = tf.physics_tools.electrons.emission
+
+
+# #####################################################
+# #####################################################
+#               DEFAULTS
+# #####################################################
+
+
+_PATH_HERE = os.path.dirname(__file__)
+_PATH_PAPER = os.path.dirname(_PATH_HERE)
+
+
+_DPFE_DCROSS = {
+    'EH0': os.path.join(
+        _PATH_PAPER,
+        'inputs',
+        'd2cross_Ee01eV-100MeV-80log_Eph1eV-100MeV-81log_ntheta61_EH.npz',
+    ),
+    'EH1': os.path.join(
+        _PATH_PAPER,
+        'inputs',
+        'd2cross_Ee0100eV-10MeV-80log_Eph100eV-10MeV-81log_ntheta61_EH.npz',
+    ),
+    'BHE': os.path.join(
+        _PATH_PAPER,
+        'inputs',
+        'd2cross_Ee01eV-100MeV-400log_Eph1eV-100MeV-401log_ntheta181_BHE.npz',
+    ),
+}
+
+
+# #####################################################
+# #####################################################
+#       Fig 1 - cross-section
+# #####################################################
+
+
+def main(
+    figsize=(15, 7),
+    pfe_cross='EH0',
+    version='EH',
+    Eph_eV=np.r_[1e3, 10e3, 500e3],
+    Ee0_eV=np.r_[20e3, 1e6],
+    fontsize=14,
+    # save
+    path_save=None,
+    pfe_save=None,
+    # unused
+    **kwdargs,
+):
+    """ Plot EH cross-section over a wide range of Eph and Ee0
+
+    Add 4 hand-picked cases vs angle of emission theta_ph
+
+    """
+
+    # --------------
+    # load
+    # --------------
+
+    pfe = _DPFE_DCROSS[pfe_cross]
+
+    dout = {
+        kk: vv.tolist()
+        for kk, vv in dict(np.load(pfe, allow_pickle=True)).items()
+    }
+    units = dout['cross'][version]['units']
+    Z = dout.get('Z', {'data': 1})['data']
+
+    # --------------
+    # prepare axes
+    # --------------
+
+    dmargin = {
+        'left': 0.06, 'right': 0.99,
+        'bottom': 0.08, 'top': 0.93,
+        'wspace': 0.25, 'hspace': 0.10,
+    }
+
+    fig = plt.figure(figsize=figsize)
+
+    gs = gridspec.GridSpec(ncols=4, nrows=2, **dmargin)
+    dax = {}
+
+    # --------------
+    # prepare axes
+    # --------------
+
+    # --------------
+    # ax - isolines
+
+    ax = fig.add_subplot(
+        gs[:, -2:],
+        xscale='log',
+        yscale='log',
+        aspect='equal',
+    )
+    ax.set_xlabel(
+        r"$E_{e,0}$ (keV)",
+        size=fontsize,
+        fontweight='bold',
+    )
+    ax.set_ylabel(
+        r"$E_{ph}$ (keV)",
+        size=fontsize,
+        fontweight='bold',
+    )
+    ax.set_title(
+        r"$d^2\sigma(E_{e0}, E_{ph}, \theta_{ph}, Z)$"
+        + f"\n Z = {Z} - version = {version}",
+        size=fontsize,
+        fontweight='bold',
+    )
+
+    # store
+    dax['map'] = {'handle': ax, 'type': 'isolines'}
+
+    # --------------
+    # ax - theta_norm
+
+    # theta_norm0
+    ax = fig.add_subplot(
+        gs[0, 0],
+        xscale='linear',
+    )
+    ax.set_ylabel(
+        "normalized cross-section (adim.)",
+        size=fontsize,
+        fontweight='bold',
+    )
+    ax.set_title(
+        r"$E_{e0}$" + f" = {Ee0_eV[0]*1e-3:2.0f} keV",
+        size=fontsize,
+        fontweight='bold',
+    )
+
+    # store
+    dax['theta_norm0'] = {'handle': ax, 'type': 'isolines'}
+
+    # theta_norm1
+    ax = fig.add_subplot(
+        gs[0, 1],
+        sharex=dax['theta_norm0']['handle'],
+        sharey=dax['theta_norm0']['handle'],
+    )
+    ax.set_title(
+        r"$E_{e0}$" + f" = {Ee0_eV[1]*1e-6:2.0f} MeV",
+        size=fontsize,
+        fontweight='bold',
+    )
+
+    # store
+    dax['theta_norm1'] = {'handle': ax, 'type': 'isolines'}
+
+    # --------------
+    # ax - theta_abs
+
+    # theta_abs0
+    ax = fig.add_subplot(
+        gs[1, 0],
+        sharex=dax['theta_norm0']['handle'],
+    )
+    ax.set_xlabel(
+        r"$\theta_{ph}$ (deg)",
+        size=fontsize,
+        fontweight='bold',
+    )
+    ax.set_ylabel(
+        r"$d\sigma$" + f"({units})",
+        size=fontsize,
+        fontweight='bold',
+    )
+
+    # store
+    dax['theta_abs0'] = {'handle': ax, 'type': 'isolines'}
+
+    # theta_abs1
+    ax = fig.add_subplot(
+        gs[1, 1],
+        sharex=dax['theta_norm0']['handle'],
+        sharey=dax['theta_abs0']['handle'],
+    )
+    ax.set_xlabel(
+        r"$\theta_{ph}$ (deg)",
+        size=fontsize,
+        fontweight='bold',
+    )
+
+    # store
+    dax['theta_abs1'] = {'handle': ax, 'type': 'isolines'}
+
+    # ------------------
+    # call built-in
+    # ------------------
+
+    # cases 100 keV
+    lc = ['r', 'g', 'b']
+    for i0, e0 in enumerate(Ee0_eV):
+        iphok = Eph_eV < e0
+        dcases = {
+            i1: {
+                'E_e0_eV': e0,
+                'E_ph_eV': eph,
+                'color': lc[i1],
+                'label': f"Eph = {eph*1e-3:3.0f} keV",
+            }
+            for i1, eph in enumerate(Eph_eV[iphok])
+        }
+        _dax, _d2cross = tfphysemis.plot_xray_thin_d2cross_ei_anisotropy(
+            d2cross=pfe,
+            dcases=dcases,
+            dax={
+                'map': dax['map']['handle'],
+                'theta_norm': dax[f'theta_norm{i0}']['handle'],
+                'theta_abs': dax[f'theta_abs{i0}']['handle'],
+            },
+        )
+
+        # remove contour plots
+        if i0 == 0:
+            for cc in dax['map']['handle'].get_children():
+                if cc.__class__.__name__ == 'QuadContourSet':
+                    cc.remove()
+
+            # remove legend
+            dax['theta_norm0']['handle'].get_legend().remove()
+
+    # ---------------------
+    # Adjust map x/y scales
+    # ---------------------
+
+    dax['map']['handle'].set_xlim(0.1, 10e3)
+    dax['map']['handle'].set_ylim(0.1, 10e3)
+    dax['theta_abs1']['handle'].set_ylabel('')
+    dax['theta_norm1']['handle'].legend(loc='lower right')
+
+    # --------------
+    # add a, b, c, d, e
+    # --------------
+
+    dabc = {
+        'theta_norm0': '(a)',
+        'theta_norm1': '(c)',
+        'theta_abs0': '(b)',
+        'theta_abs1': '(d)',
+    }
+    for kax, abc in dabc.items():
+        dax[kax]['handle'].grid(visible=True, which='major', axis='both')
+        dax[kax]['handle'].text(
+            0.95, 0.95,
+            abc,
+            fontsize=fontsize,
+            fontweight='bold',
+            horizontalalignment='right',
+            verticalalignment='top',
+            transform=dax[kax]['handle'].transAxes,
+        )
+
+    dax['map']['handle'].text(
+        0., 1.02,
+        "(e)",
+        fontsize=fontsize,
+        fontweight='bold',
+        horizontalalignment='left',
+        verticalalignment='bottom',
+        transform=dax['map']['handle'].transAxes,
+    )
+
+    # --------------
+    # save
+    # --------------
+
+    savefig(
+        fig=fig,
+        pfe_save=pfe_save,
+        path_save=path_save,
+        file=__file__,
+    )
+
+    return dax