Developing #184
Analyzing #183: Express analysis of experimental data
Калибрация тонково кремневого детектора
Status: | Открыта | Start date: | 04/13/2018 | ||
---|---|---|---|---|---|
Priority: | Нормальный | Due date: | 04/16/2018 | ||
Assignee: | - | % Done: | 30% | ||
Category: | Software | ||||
Target version: | - |
Description
Получить калибровочные параметры для 20 микронного детектора. Калибровочные измерения проводились с источником альфа частиц 226Ra.
Измерение 1:¶
Источник ориентирован к поверхности детектора под углом ~90° на расстоянии 40 cm. Для некоторых пиков ожидается прострел тонково детектора и выделение энергии в милиметровом детекторе. Записывались два тригера:
- от тонково слоя - 1
- от милиметрового детектора - 3.
файл:
159.93.80.161:/home/LynxOS/mbsusr/mbsdaq/mbsrun/exp201804/data/calib/si_20_03_*.lmd
Измерение 2:¶
Источник на расстоянии 40 cm от тонкого детектора ориентирован под углом примерно 58°, чтобы обеспечить полную остановку трех пиков в первом слою. Измерение записано в файле
159.93.80.161:/home/LynxOS/mbsusr/mbsdaq/mbsrun/exp201804/data/calib/si_20_04_*.lmd
History
#1 Updated by Vratislav Chudoba over 6 years ago
- Description updated (diff)
#2 Updated by Vratislav Chudoba over 6 years ago
- File sqx_l_ec.clb added
- Description updated (diff)
Измерение 1:¶
Калибровать будет сложно.
На картинке видим соотношение тригеров (1 и 3), множественность сработавших стрипов в 1 мм детекторе (при условии офлайн порога 0.2 МэВ для каждого стрипа) и суммарную энергию по всем стрипам (калибрация проведена с помощью файла, который использовался во время эксперимента для онлайн просмотра, см. прикрепленный файл sqx_l_ec.clb ).
После отсечения офлайн порогов в стрипах тонково детектора (для каждого 120 каналов АЦП) видим мизерную статистику для проведения калибровки (около 1800 физических событий в стрипе) и очень плохо разделенные пики.
текст
#3 Updated by Vratislav Chudoba over 6 years ago
- Description updated (diff)
#4 Updated by Ivan Muzalevsky over 6 years ago
- File parforcal.txt added
- File calPar.txt added
- File instruction.pdf added
#5 Updated by Ivan Muzalevsky over 6 years ago
- % Done changed from 0 to 30
Для получения мёртвого слоя делалось следующее:
1) находились калибровочные параметры для трёх пиков альфа частиц в данных, собранных при облучении 20 мк детектора под углом 58 градусов.
Находить положения трёх пиков пришлось вручную, не используя программу автоматической калибровки AculUtils, из-за малой статистики и сложности выделения пиков из остального забора. При определении позиции пиков учитывалась, что все три пика во всех каналах лежат примерно одних и тех же диапазонах.
Из картинок видно, что три пика можно наблюдать только в первых 14 стрипах, вероятно, что последние два стирпа простреливались альфой с энергией 6.0024 [МэВ], так как для них угол, относительно которого было произведено облучение был самым маленьким, по отношению к остальным в силу законов геометрии. соответственно, и длина материала, через который они пробегали был меньше.
С помощью AculUtils было посчитано, что длина пробега альфа частицы с энергией 6.0024 МэВ равна 33.1811 микрон. ( точность рассчёта приемлимая. В хорошем согласии с LISE)
Учитывая, что ширина изучаемого детектора 6 см, было рассчитано, что минимальный угол, под которым альфа частицы влетают в материал детектора ~55 градусов. При такой траектории альфа частица проходит через 35 микрон кремния. Величина, близка к 33, к тому же, не учитывалось пространственное размытие источника альфа частиц, потери в защитной коробочке источника...
2) Для калибровки брались следующие энергии альфа частиц : 4.7843, 5.48948, 6.0024 [МэВ] (http://www.khlopin.ru/docs/products/source/E_OSAI.pdf).
Для них были получены предварительные калибровочные параметры (preCalPar.txt).
На картинке представлены линейные фиты по трём точкам ( мы не учитывали пьедестал, не учитывали точку (0,0), делали калибровку только по трём пикам)
3) С помощью этих параметров и положений пиков, найденных в пункте 1, были пересчитаны энергии трёх альфа частиц для тех же самых данных для каждого стрипа. (эти мелкие поправки в конце рассчётов дают достаточно существенный вклад)
4) Для данных, собранных под нулём градусов, было обнаружено, что единственный пик, соответствует единственной непростреленной альфе с минимальной энергией. Было обнаружено, что события, заполняющие этот пик не дают сигнала в 1-мм детекторе, что не справедливо для остальных событий. Такжеб ыло рассчитано, что длина пробега альфа частицы с энергией 4.7843 МэВ равняется 23.9217 микрон (также не учитываются некоторые потери).
Находилось положение пика в каналах. Было обнаружено, что пик наблюдается в стрипах с 2 по 11. Положение пика также находилось вручную.
5) С помощью калибровочных параметров (preCalPar.txt), получалось значение энергии альфы, соответствующей этому пику.
6) Из средней разницы энергий для альфы, в измерениях под нулём и под 58 градусов высчитывалось разница длин пробега альфа частиц и, соответственно значение толщины мёртвого слоя. Которое оказалось равно
0.519412 мкм
7) С помощью рассчитанного мёртвого слоя, высчитывались поправки к энергиям трёх альфа частиц. Полученные энергии: 4.71199, 5.42139 и 5.93722 МэВ. С новыим энергиями были рассчитаны окончательные, на данный момент, калибровочные параметры (SQ20_58.cal). Линейный фиты выглядели примерно также как и при получении предварительных параметров.
все рассчёты представлены здесь. лист в таблице SQ20 : https://docs.google.com/spreadsheets/d/1RjBBLo_H3SWNhnWMjvQQSxE3np8uQBSFwvwXzX3idto/edit?usp=sharing
#6 Updated by Ivan Muzalevsky over 6 years ago
- File deleted (
parforcal.txt)
#7 Updated by Ivan Muzalevsky over 6 years ago
- File deleted (
calPar.txt)
#8 Updated by Ivan Muzalevsky over 6 years ago
- File deleted (
instruction.pdf)
#9 Updated by Vratislav Chudoba over 6 years ago
После оценки переднего мертвого слоя в тонком детекторе можно предположить, что задний мертвый слой будет такой же.
С энергиями поправленными на толщину мертвого слоя (под углом будет эффективная толщина отличаться) получить окончательные калибровочные параметры тонкого детектора на основе измерений под углом (там, где видны три пика).
Делаем следующие упрощения:
- толщину мертвого слоя берем однорудную по всей площади детектора
- поправку энергий делаем одинаковую (в зависимости от среднего угла) для всех стрипов
Проверку калибровки можно сделать следующим образом: выводим картинку сыммы депозита в тонком и депозита в толстом детекторе и должны увидеть 4 пика соотвествующие 226Ra на своих местах.
#10 Updated by Ivan Muzalevsky over 6 years ago
- File preCalPar.txt added
#11 Updated by Ivan Muzalevsky over 6 years ago
- File SQ20_58.cal added
#12 Updated by Ivan Muzalevsky over 6 years ago
- File fill1mm.C added
Vratislav Chudoba wrote:
После оценки переднего мертвого слоя в тонком детекторе можно предположить, что задний мертвый слой будет такой же.С энергиями поправленными на толщину мертвого слоя (под углом будет эффективная толщина отличаться) получить окончательные калибровочные параметры тонкого детектора на основе измерений под углом (там, где видны три пика).
Делаем следующие упрощения:
- толщину мертвого слоя берем однорудную по всей площади детектора
- поправку энергий делаем одинаковую (в зависимости от среднего угла) для всех стрипов
Проверку калибровки можно сделать следующим образом: выводим картинку сыммы депозита в тонком и депозита в толстом детекторе и должны увидеть 4 пика соотвествующие 226Ra на своих местах.
Попытался проверить так, как вы описали (метод заполнения деревьев можно увидеть в макросе fill1mm.C), получил следующую картинку:
#13 Updated by Ivan Muzalevsky over 6 years ago
Ivan Muzalevsky wrote:
Vratislav Chudoba wrote:
После оценки переднего мертвого слоя в тонком детекторе можно предположить, что задний мертвый слой будет такой же.С энергиями поправленными на толщину мертвого слоя (под углом будет эффективная толщина отличаться) получить окончательные калибровочные параметры тонкого детектора на основе измерений под углом (там, где видны три пика).
Делаем следующие упрощения:
- толщину мертвого слоя берем однорудную по всей площади детектора
- поправку энергий делаем одинаковую (в зависимости от среднего угла) для всех стрипов
Проверку калибровки можно сделать следующим образом: выводим картинку сыммы депозита в тонком и депозита в толстом детекторе и должны увидеть 4 пика соотвествующие 226Ra на своих местах.
Попытался проверить так, как вы описали (метод заполнения деревьев можно увидеть в макросе fill1mm.C), получил следующую картинку:
Всё верно, получал точно такую же. Вроде бы мы лично это обсуждали. Так как я не увидел нормальные 4 пика, не стал это сюда выкладывать.
#14 Updated by Ivan Muzalevsky over 6 years ago
- File SQ20_58(14).cal added
Файл с калибровочными параметрами был приведён к go4-овскому виду.
#15 Updated by Ivan Muzalevsky over 6 years ago
From the data obtained in meauserements described in task 197 dead layer was calculated. calculations was made https://docs.google.com/spreadsheets/d/1RjBBLo_H3SWNhnWMjvQQSxE3np8uQBSFwvwXzX3idto/edit#gid=38392088 on the page Si20_new. calculated dead layer is about 0,864403939 um. Data from strips 3-6 were taken into account for calculations.
Using dead layer calibration parameters were calculated. data from 60 deg meausrements was used for this.
Linear calibration was made using 3 points. 3 peak positions in raw spectra was determined using TScectrum methods (showpeaks, getX position of polimarker). This time usual method of fitting gaus function was not used.
For different strips different effective value of dead layer was used (different strips were at different angles).
linear fits were looked pretty nice.
calibration parameters are written into SQ20_may18.cal
Same algorithm of calibration was used for 45 deg data (see task 197). The difference between new 45 deg calibration parameters and 60 deg pars was less than 5%. This error is considered as acceptable at the moment.
#16 Updated by Ivan Muzalevsky over 6 years ago
- File SQ20_may18.cal added