CC BY
19
84
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2013. № 1
КРАТКИЕ СООБЩЕНИЯ ФИЗИКА АТОМНОГО ЯДРА И ЭЛЕМЕНТАРНЫХ ЧАСТИЦ
Автоматизация измерений и обработки результатов в экспериментах
по методике наведенной активности
С. С. Белышев1,0, К. А. Стопани2,6
Московский государственный университет имени М. В. Ломоносова, 1 физический факультет, кафедра общей ядерной физики;
2 Научно-исследовательский институт ядерной физики имени Д. В. Скобельцына.
Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
E-mail: а belyshev@deprii.sinp.msu.ru, ь hatta@depni.sinp.msu.ru Статья поступила 19.07.2012, подписана в печать 26.07.2012.
Описана система сбора и анализа гамма-спектров, используемая при исследовании фотоядерных реакций при помощи метода наведенной активности. Она состоит из базы данных, которая осуществляет автоматический набор спектров и предоставляет сетевой доступ к полученным данным, и программы для анализа спектров. Использование данной системы существенно упрощает проведение экспериментов и обработку их результатов.
Ключевые слова: ядерные реакции, активационный анализ, гамма-спектрометрия.
УДК: 539.1.074.55. PACS: 25.85,+с.
Введение
Метод наведенной активности часто применяется при изучении ядерных реакций, таких, как реакции под воздействием 7-квантов [1], и для проведения активационного анализа состава вещества [2]. Данная методика основана на том, что многие изотопы, образующиеся в результате ядерных реакций, являются нестабильными и распадаются в течение некоторого времени после образования. Поскольку наиболее типичным видом распада для большинства ядер является бета-распад, в котором образовавшееся в возбужденном состоянии дочернее ядро может испускать гамма-кванты, то по спектру гамма-квантов возможны точная идентификация и определение количества распавшихся ядер.
Эксперимент по методике наведенной активности состоит из двух этапов. Вначале проводится облучение образца потоком налетающих частиц (например, фотонами в гамма-активационном анализе), во время которого в образце протекают ядерные реакции с образованием нестабильных продуктов. На втором этапе эксперимента с помощью детектора гамма-квантов измеряются спектры фотонов, испускаемых облученным образцом. Для нахождения активности образовавшихся нестабильных изотопов в спектрах проводятся поиск и определение площадей пиков. Затем полученные данные могут использоваться для определения выходов ядерных реакций, количества исходных ядер в составе образца и т. д.
При практическом использовании метода наведенной активности для повышения точности результатов вместо измерения одного спектра гамма-квантов обычно делается серия последовательных измерений, что позволяет получить зависимость интенсивности распадов от времени, т. е. кривую распада. Для макси-
мальной точности результатов обработки необходимо, чтобы время измерения каждого спектра было значительно меньше периода полураспада, а суммарная продолжительность измерений должна быть велика по сравнению с ним [3]. Таким образом, для получения наиболее качественных результатов требуется провести большое количество последовательных измерений спектров, каждое из которых имеет небольшую длительность. Отметим, что практически в методе наведенной активности доступны для наблюдения распады с периодами от десятков секунд до нескольких лет, поэтому число измеряемых и обрабатываемых спектров может достигать нескольких тысяч [4].
В настоящей работе описывается система, используемая для автоматизации процесса сбора и анализа большого числа спектров при проведении экспериментов по методике наведенной активности в НИИЯФ МГУ. Схема взаимодействия компонентов системы показана на рис. 1. В ее состав входит база данных, которая осуществляет автоматический набор спектров и предоставляет сетевой доступ к полученным данным, и программа для анализа спектров и определения параметров пиков.
Сбор и хранение спектров в базе данных
Спектры гамма-квантов измеряются с помощью германиевого детектора Canberra GC 3019, подключенного к блоку многоканального анализатора InSpector 1250 (МКА). Обмен данными с МКА осуществляется по закрытому протоколу через интерфейс RS232. Для управления детектором, получения с него данных и для анализа спектров фирмой-производителем поставляется программа Genie 2000 [5], которая предназначена для работы в интерактивном режиме, удобном для повседневных измерений, но плохо подходящем для про-
Клиент
Рис. 1. Система накопления и анализа гамма-спектров
ведения большого количества однотипных измерений. Вместо нее для последовательных измерений большого количества спектров с минимальными промежутками между измерениями в системе сбора данных применяется программа-сценарий, использующая отдельные компоненты Genie 2000 для управления МКА, запуска и остановки набора спектров и для получения измеренных спектров из МКА. Результатом ее работы являются файлы со спектрами во внутреннем формате Genie 2000, сохраняемые на компьютере управления детектором.
Для хранения спектров используется база данных PostgreSQL [6], размещенная на отдельном сетевом сервере. Полученные с детектора файлы спектров преобразуются в массив целых чисел и сохраняются в таблице БД. Также сохраняется информация о времени измерения спектра, идентификаторы исследуемого образца и текущего облучения, мертвое время детектора и его идентификатор. Преимущества использования СУБД для хранения спектров заключаются в возможности гибкой и быстрой работы с большим количеством данных (в данный момент в БД содержится около 30000 спектров), возможности сортировки и поиска информации.
Доступ к базе данных возможен через интернет по протоколу HTTP. Запрос специального формата позволяет получить спектр в текстовом или графическом виде, а также провести сложение спектров в некотором диапазоне времени измерения. Также для этих операций был написан пользовательский web-интерфейс, при помощи которого доступ к спектрам может осу-
ществляться через \уеЬ-браузер. В качестве \уеЬ-сервера используется Ц^Ы^рс! [7], формирование интерфейса, обработка запросов, и генерация данных выполняются модулем РаБ1:С01 [8]. Для обеспечения достаточной скорости сложения и преобразования в текстовый формат большого количества спектров данный модуль написан на языке С++.
Автоматический анализ спектров
Программа анализа спектров позволяет в измеренных спектрах проводить поиск отдельных пиков, определять их характеристики и строить кривые распада. Анализ выполняется в интерфейсе \уеЬ-браузера, и переход к нему возможен непосредственно из просмотра спектров в базе данных после выбора необходимых спектров. Страница анализа с результатами поиска пиков в спектрах 60 Со, измеренных в течение 10 мин на разных расстояниях от детектора, схематически изображена на рис. 2.
Для каждого выбранного спектра приводится таблица, содержащая параметры найденных пиков: энергия, ширина, площадь, а также параметры, определяющие форму пика и характеристики фоновой подложки. Кроме того, для контроля правильности поиска результаты анализа каждого спектра показываются на графике. Под результатами анализа спектров выводится сводная таблица, в которой приводятся площади всех найденных пиков, упорядоченные по номеру спектра. При нажатии на кнопку «Построить кривую распада» информация из таблицы выводится в отдельном окне
86
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2013. № 1
Спектр 1, 2012-05-31, 600.0 сек, живое время 249.39 сек:
Найденные пики: Е, кэВ Г, кэВ 1173.04 1.92 1332.28 2.01 2505.12 2.61
Площадь пика 457176.54 407059.14 22281.96
Погрешность 689.579 613.489 154.100
1400
Спектр 2, 2012-06-01, 600.0 сек, живое время 579.55 сек:
1100 Найденные пики:
Е, кэВ Г, кэВ 1172.90 1.71 1332.13 1.82
1150
1200
Сводная таблица: Е, кэВ Спектр 1 Спектр 2 1172.94 457176.54 42314.45 1332.17 407059.14 38515.31 2505.12 22281.96
1250
1300
Площадь пика
42314.45
38515.31
Погрешность
208.315
196.434
1350
1400
Постр. кр. распада Постр. кр. распада
Рис. 2. Результат анализа двух спектров Со (схематическое изображение работы программы РЕАКОЕСО)
в виде графика скорости счета в пике в зависимости от времени вместе с результатами аппроксимации этих данных экспоненциальной функцией (для экономии места на рисунке показана обработка двух спектров, но возможность построения кривой распада предназначена для одновременной обработки большого количества спектров).
Программа анализа спектров выполняется на отдельном web-cepвepe в виде СХН-сценария сШМ, который формирует интерфейс и запускает написанную на С++ программу РЕАКОЕСО, содержащую математическую часть алгоритма обработки данных. Нахождение параметров пиков разбито на три этапа: выделение фона, поиск максимумов характерной формы, выделяющихся над фоном, и аппроксимация найденных максимумов модельной функцией.
Для нахождения приближенного значения фона в каждой точке спектра используется следующий метод. Выбирается ширина окна N. которая заведомо больше ширины пика. В диапазон точек шириной N в таком случае попадет какая-то часть фона. В данном диапазоне находится точка с минимальным значением числа отсчетов Ьт-т. Если бы в данном диапазоне отсутствовали пики, а среднее значение фона было бы одинаковым и равным Ь, то вероятность получить данное Ьт-т определялась бы следующим выражением:
Р(Ьт[п) =Мр(Ьт1п)(р^(Ьт1п))ы-\ (1)
где р(Ьгшп) — вероятность получить Ьт-т в одном канале спектра, равная в соответствии с законом Пуассона
е~ьЬь™п/(Ьт-т\). При достаточно больших значениях Ьтт ^ Ю от дискретного распределения Пуассона можно перейти к непрерывному нормальному распределению со средним значением М = Ьт-т и стандартным отклонением а = \/Ьт-т. Тогда среднее значение Ьт-т можно рассчитать по формуле
Ьгп\гл -
ЬттР{Ьтт)
(2)
а затем, приравняв Ьт-т = Ьт-т, получить уравнение, которое численно решается относительно Ь, являющегося искомым приближенным значением фона.
Вторым этапом анализа спектра является поиск в за-шумленном спектре максимумов, соответствующих по форме искомым пикам и обладающих шириной, близкой к разрешению детектора. Для этого используется алгоритм, аналогичный использованному в программе БАМР080 [9]. В каждой точке спектра вычисляется
1+Ь
величина р(£/) = X) А/&5(£ь £/), где
5(£,£0) =
(Е—Е0)2/2о2
сг*
(3)
— вторая производная функции Гаусса, используемая как критерий соответствия спектра в данной точке пику с центром в точке Ео и дисперсией а2. Пикам в исходном спектре соответствуют максимумы функции р(£/).
На третьем этапе параметры найденных максимумов используются в качестве начальных значений для процедуры аппроксимации участков спектра функциями следующего вида, содержащими слагаемые, которые характеризуют общий для всех пиков линейный фон, сами пики, имеющие форму функции Гаусса, и фон под каждым пиком по отдельности:
f(E) = aE + b +
где erf — функция ошибок; М — число пиков на данном участке спектра; а, Ь, /3; — параметры аппроксимации фона; Аи E-L, a-L — площадь, средняя энергия и стандартное отклонение каждого пика соответственно. Для нахождения значений параметров используется метод максимального правдоподобия. Аппроксимация каждого участка спектра выполняется несколько раз с различным набором пиков в аппроксимирующей функции, после чего выбирается наилучший вариант, что позволяет разрешать не только одиночные пики, но и мультиплеты, в которых несколько пиков накладываются друг на друга. Дополнительно для разделения мультиплетов может быть задействован метод деконволюции, описанный в [10]. Для численной максимизации функции правдоподобия используется пакет Minuit2 [11]. Кроме автоматического поиска пиков возможен ручной ввод начальных параметров аппроксимации.
Заключение
В статье описана компьютерная система сбора и анализа данных активационных экспериментов, применяемая в НИИЯФ МГУ. Архитектурно данная система состоит из набора модулей, каждый из которых вы-
полняет ограниченный набор функций и может выполняться независимо на отдельном компьютере в сети. Достоинствами данной схемы являются гибкость и расширяемость — за счет увеличения числа узлов она может быть сделана произвольно сложной. Так, в данную систему тривиально включается еще один детектор или компьютер обработки.
Построение системы сбора и анализа данных в виде веб-приложения позволяет проводить мониторинг и обработку экспериментальных данных удаленно, что упрощает и ускоряет этот процесс, а централизация хранения измеренных спектров в единой базе данных упрощает работу с большим количеством измеренных спектров.
Список литературы
1. Варламов В.В., Ишханов B.C., Четверткова В.А. // Изв. РАН. Сер. физ. 2010. 74, № 6. С. 875.
2. Bujdosö Е. // J. Radioanal. Nucí. Chem. 1997. 223, N 1-2. P. 251.
3. Велышев С.С., Стопани К.А., Трощиев С.Ю. и др. // Вестн. Моск. ун-та. Физ. Астрон. 2011. № 4. С. 42.
4. Стопани К.А., Велышев С.С., Кузнецов A.A., Кури-лик A.C. // 62-е Междунар. совещ. по ядер, спектроскопии и структуре ядра «Ядро-2012». Воронеж, 2012.
5. Genie 2000 Applications Software. Canberra Industries Inc. http://www.canberra.com/products/831.asp.
6. PostgreSQL: The world's most advanced open source database, http://www.postgresql.org/.
7. Lighttpd. http://www.lighttpd.net/.
8. FastCGI Specification // Mark R. Brown Open Market, Inc. http://fastcgi.com/devkit/doc/fcgi-spec.html.
9. Aarnio P., Nikkinen M., Routti J. SAMPO Advanced Gamma Spectrum Analysis Software. Version 3.62, 18.12.1998. Logion Oy, Helsinki, Finland.
10. Morhac M. et al. // Nucí. Instr. and Meth. A. 1997. 401. P. 385.
11. James F., Winkler M. MINUIT User's Guide. Geneva: CERN, 2004.
Automatic data acquisition and analysis in activation experiments S.S. Belyshev10, K. A. Stopani2 ft
1 Department of General Nuclear Physics, Faculty of Physics;
2 D. V. Skobeltsyn Institute of Nuclear Physics,
M. V. Lomonosov Moscow State University, Moscow 119991, Russia. E-mail: a belyshev@depni.sinp.msu.ru, b hatta@depni.sinp.msu.ru.
The system for automatic acquisition and analysis of gamma-ray spectra that is used in the studies of photonuclear reactions with activation technique is described. The system is comprised of the database which is used to access the data over the network and of the spectrum analysis program. Experimental procedure and data analysis are greatly simplified with the use of the described system.
Keywords: nuclear reactions, activation analysis, gamma-ray spectrometry. PACS: 25.85,+c. Received 19 July 2012.
English version: Moscow University Physics Bulletin 1(2013).
Сведения об авторах
1. Белышев Сергей Сергеевич — физик; тел.: (495) 939-25-58, e-mail: belyshev@depni.sinp.msu.ru.
2. Стопани Константин Александрович — мл. науч. сотрудник, тел. (495) 939-25-58, e-mail: hatta@depni.sinp.msu.ru.
