THEORY2

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== Классическое восстановление энергии в моно режиме ======

Классическое восстановление энергии первичных частиц высоких энергий, регистрируемых телескопами **TAIGA-IACT** в моно режиме представляет собой линейную регрессию на основе lookup таблиц, полученных их моделирования. Весь диапазон возможных значений dist разбивается на небольшие бины, например 0.25 градуса. Для событий в каждом бине определяется линейная зависимость между log(size) и log(energy). Коэффициенты линейной зависимости записываются в lookup таблицу.

==== Структура lookup таблицы ====

^ dist ^ N_points ^ a ^ b ^
| dist - середина бина по dist | N_points - число событий в данном бине | a - коэффициент линейной связи в уравнении log(energy) = a*log(size) + b | b - коэффициент смещения в уравнении log(energy) = a*log(size) + b |

==== Функция на Python ====

<code python>
def reco_mono_energy(dist_value, size_value):
    mask = np.abs(energy_lookup['dist'] - dist_value) <= bin_size
    if not mask.any():
        raise KeyError(f"Нет подходящего интервала для dist = {dist_value}")
    row = energy_lookup.loc[mask].iloc[0] 
    return 10**(np.log10(size_value) * row['a'] + row['b'])
</code>

Энергетическое разрешение при таком подходе достигает в лучшем случае 25%, если не делать отбор событий по каким-либо параметрам.

{{.:pics:20250627-183512.png}}

Подробнее с процедурой восстановления энергии в моно режиме можно ознакомиться в jupyter notebook: \k38\taiga_pool\Hybrid_events\for_Krukov\Linear_reg_reco_en.ipynb

В нем реализована процедура создания lookup таблицы и затем её использование на тестовом наборе данных для восстановления энергии.