Métricas para comparar inyecciones de distintas especies (Fe, Ni, Mo, W) sin
depender de la escala absoluta de cada descarga. Se usa en
notebooks/10W1-10W2/10W2.ipynb.
from aug.analysis.compare import comparison_table, residual_profiles, collapsed_curves, multiline_timing
res = group.analyze_injections()
cmp = comparison_table(res, group) # res.table + métricas nuevas
cmp[cmp.clean].groupby(["species", "regime"])[["curvature", "rise_over_tau", "tau_over_tauE"]].median()
Columnas que agrega comparison_table
| Columna | Definición | Lectura |
|---|---|---|
Z |
número atómico (SPECIES_Z) |
|
rise_over_tau |
t_onset_to_peak_ms / tau_imp_ms |
subida relativa al decaimiento; igual en todos los grupos si solo cambia la escala de tiempo |
tau_A_ms, tau_B_ms |
τ ajustado en [0, τ] y [τ, 2.5τ] desde el pico, con el offset b fijo |
|
curvature (κ) |
tau_B / tau_A |
1 = exponencial pura; < 1 = hombro (cae más rápido después); > 1 = cola lenta |
tauE_ms |
W_mhd / (P_ECRH + P_NBI + P_OH) |
orden de magnitud (sin radiación ni dW/dt) |
tau_over_tauE |
tau_imp / tau_E |
τ_imp normalizado |
dTe_pct, dW_pct |
cambio % de Te central (ECE) y W_mhd entre −100, 0 ms y 30, 120 ms del onset | enfriamiento por la inyección |
Otras funciones
decay_shape(detail): κ de una inyección (res.details[(shot, inj)]).segment_tau(tr, z, t0, t1): τ deA exp(-(t-t0)/τ)en un tramo (con el
offset ya restado).residual_profiles(res, x_edges): residuo medio(y − modelo)/adel ajuste
simple por intervalo dex = t/τ. Promediado por grupo, muestra
desviaciones sistemáticas de la forma exponencial.collapsed_curves(res, x_grid):(y − b)/avs.x = t/τ. Si todo es
exponencial, las curvas colapsan enexp(−x).multiline_timing(shot, t_ref, t_onset_core): onset y pico de cada línea de
la especie (SPECIES_LINES, de carga baja a alta) relativos al onset de la
línea de núcleo. Mide la penetración desde el borde (carga baja) hasta el
núcleo (carga alta). Para W solo existeW44_133.tau_E_estimate(shot, t_onset),plasma_response(shot, t_onset).
Notas
P_OHviene en MW en casi todos los-rest.txt, y en W en I-mode Mo
#38858/9. Los valores > 1e3 se toman como W.- Las métricas de forma solo tienen sentido en inyecciones limpias (
clean) y
con ventana ≥ 2.5 τ, para quetau_Besté definido.