This is an old revision of the document!
Классическое восстановление энергии первичных частиц высоких энергий, регистрируемых телескопами TAIGA-IACT в моно режиме представляет собой линейную регрессию на основе lookup таблиц, полученных их моделирования. Весь диапазон звозможных значений dist разбивается на небольшие бины, например 0.25 градуса. Для событий в каждом бине определяется линейная зависимость между log(size) и log(energy). Коэффициенты линейной зависимости записываются в lookup таблицу. Заголовок lookup таблицы следующий: dist,N_points,a,b dist - середина бина по dist N_points - число событий в данном бине a и b - коэффициенты линейной связи log(energy) = a*log(size) + b Функция на python для восстановления энергии на основе lookup таблицы может выглядеть следующим образом: def reco_mono_energy(dist_value, size_value):
# маска строк, где |lookup.dist - val| <= tol
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'])
Энергетическое разрешение при таком подходе достигает в лучшем случае 25%, если не делать отбор событий по каким-либо параметрам.
Подробнее с процедурой восстановления энергии в моно режиме можно ознакомиться в jupyter notebook: L:\k38\taiga_pool\Hybrid_events\for_Krukov\Linear_reg_reco_en.ipynb
В нем реализована процедура создания lookup таблицы и затем её использование на тестовом наборе данных для восстановления энергии.