Класс TSpectrum

В 5 версии ROOT появился очень интересный класс TSpectrum. Он позволяет легко выполнять наиболее частые операции, которые требуются при анализе физических данных. Класс TSpectrum содержит следующие функции обработки спектров:
  1. Оценка фона
  2. Сглаживание
  3. Поиск пиков
  4. Деконволюция
  5. Трансформации
  6. Визуализация
Можно исследовать одномерные и двумерные спектры. Подробное описание методов класса и алгоритмов, на которых они построены, можно найти здесь. Следует отметить, что реализация некоторых методов зависит от версии ROOT.

Конструктор класса

      TSpectrum::TSpectrum(Int_t maxpositions)
            
maxpositions-максимальное число пиков, которое Вы хотите найти в спектре. Этот параметр можно не задавать.

Функции поисков пиков

      Int_t Search(const TH1* hist, Double_t sigma=2,Option_t* option, Double_t treshold = 0.05)
            

Эта функция выполняет поиск пиков в спектре источника, отсекая непрерывный фон и статистические колебания - шум. Количество найденных пиков и их позиции записываются в члены fNpeaks и fPositionX. Поиск выполняется в текущем диапазоне гистограмм. Алгоритм основан на сглаженных вторых производных, которые сравниваются со стандартными отклонениями. Поэтому сигма также передается в функцию поиска пика. Параметры:

Эта функция позволяет отделять бесполезную спектральную информацию (непрерывный фон) от пиков. Функция основана на чувствительном алгоритме нелинейного итерационного пикового отсечения, фактически представляющего собой разностный фильтр второго порядка (-1,2, -1). Алгоритм подробно описан в ...