Developing #167
Updated by Ivan Muzalevsky over 6 years ago
<p>На данный момент в методах класса ERDecayEXP1803 для бинарной реакции используются методы TGenPhaseSpace, в которой заложено, что угловое распределение вылетающих частиц - изотропно. </p>
<p>Добавить методы, позволяющие разыгрывать реакцию, использую заданную пользователем функцию углового распределения.</p>
<h2>Результат:</h2>
<p style="background-color: rgb(255, 255, 255);"><span style="background-color: rgb(255, 255, 221);">Вся функциональность проверялась в ветке </span><a href="https://github.com/ExpertRootGroup/er/compare/168_EventsWith6HeZeroSpeed" style="color: inherit; word-wrap: break-word; background-color: rgb(234, 245, 255); font-family: -apple-system, BlinkMacSystemFont, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", "Segoe UI Symbol"; font-weight: 600;">168_EventsWith6HeZeroSpeed</a></p>
<p style="background-color: rgb(255, 255, 255);"><span style="color: rgb(88, 96, 105); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; white-space: nowrap;">commit </span><span class="sha user-select-contain" style="font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; color: rgb(68, 77, 86); white-space: nowrap;">2a576fbe7a34efde6a5cdf6ab3177af99fedf40f</span></p>
<p style="background-color: rgb(255, 255, 255);"><span class="sha user-select-contain" style="font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; color: rgb(68, 77, 86); white-space: nowrap;">макросы симуляции, обработки и отрисовки прикреплены. Функциональность была доавлена в класс ERDecayEXP1803. Добавленные методы основаны на </span>http://old.pskgu.ru/ebooks/djekson/djekson_12.pdf .</p> ERDecayEXP1803.</span></p>
<p style="background-color: rgb(255, 255, 255);"><img alt="" src="http://er.jinr.ru/develop/system/rich/rich_files/rich_files/000/000/059/original/LISE%20output1.png" style="vertical-align: middle; max-width: 100%;" /></p>
<p style="background-color: rgb(255, 255, 255);">В симуляции использовалось распределение энергии возбуждения 5H состоящее из двух гауссов: "mean"="0.4" "shape"="gauss" "strength"="1" "width"="0.02355"<br />
"mean"="1.2" "shape"="gauss" "strength"="1" "width"="0.2355". Значение, разыгрывающееся случайным образом из описанного распределения использовалось как добавка к массе оснровного состояния 5H, равного 4.69036244; // [GeV]. масса основного состояния задается пользователем в макросе симуляции. Исходя из заданного распределения, понятны раздвоения линий.</p>
<p style="background-color: rgb(255, 255, 255);">При разыгрывания бинарной реакции, в случае, если раыгранная энергия возбуждения превышает порог реакции, то энергия возбуждения разыгрывается ещё раз, пока бинарная реакция не сможет произойти. </p>
<p style="background-color: rgb(255, 255, 255);">Для разыгрывания бинарной реакции (6he + h2 - > 5H + 3He) использовалось распределение полярного угла вылета продуктов реакции в с.ц.м реакции. Распределение считывалось из прикрепленного файла.</p>
<p style="background-color: rgb(255, 255, 255);"><img alt="" src="http://er.jinr.ru/develop/system/rich/rich_files/rich_files/000/000/060/original/LISE%20output2.png" style="vertical-align: middle; max-width: 100%;" /></p>
<p style="background-color: rgb(255, 255, 255);">Рассматривались события, соответствующие второму максимуму углового распределения. (зелёные линии). Как и в рассчётах SIMONE, в данном случае, He3 в лабе вылетает под углами (21,33) в лабе.</p>
<p style="background-color: rgb(255, 255, 255);">Все результаты совпадают с полученными в SIMONE. </p>
<p style="background-color: rgb(255, 255, 255);">Используется нестандартный для er метод разыгрывания AD для бинарной реакции, не используются методы класса ERTextDecay. Цитата Виталика по этому поводу: "<span style="background-color: rgb(255, 255, 255);">Относительно ERTextDecay - давайте чтение из текстового файла применять, только тогда, когда не можем запрограммировать быстро и вычислительно просто результат распада. В этом случае получилось. Иван молодец.(шутка)</span>" </p>
<h2 style="padding-top: 2px; padding-left: 0px; background-color: rgb(255, 255, 255);"><span style="font-size: 14px;">Существующие проблемы</span></h2>
<p style="background-color: rgb(255, 255, 255);"><span style="background-color: rgb(255, 255, 255);">Возникла проблема, скорее всего связанная с использованием генератора случайных чисел, при разыгрывании угла вылета продуктов бинарной реакции в сцм. Double_t thetaCM = fADFunction->GetRandom(1.,150.)*TMath::DegToRad(); Было выяснено, что иногда (при большой статистике > 10000 событий) возникают события, в которых кинетическая энергия 6He равна 0, поэтому симуляция зацикливается на этапе разыгрывания бинарной реации. Задача перенесена </span>http://er.jinr.ru/develop/issues/169</p>
<p>Добавить методы, позволяющие разыгрывать реакцию, использую заданную пользователем функцию углового распределения.</p>
<h2>Результат:</h2>
<p style="background-color: rgb(255, 255, 255);"><span style="background-color: rgb(255, 255, 221);">Вся функциональность проверялась в ветке </span><a href="https://github.com/ExpertRootGroup/er/compare/168_EventsWith6HeZeroSpeed" style="color: inherit; word-wrap: break-word; background-color: rgb(234, 245, 255); font-family: -apple-system, BlinkMacSystemFont, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", "Segoe UI Symbol"; font-weight: 600;">168_EventsWith6HeZeroSpeed</a></p>
<p style="background-color: rgb(255, 255, 255);"><span style="color: rgb(88, 96, 105); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; white-space: nowrap;">commit </span><span class="sha user-select-contain" style="font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; color: rgb(68, 77, 86); white-space: nowrap;">2a576fbe7a34efde6a5cdf6ab3177af99fedf40f</span></p>
<p style="background-color: rgb(255, 255, 255);"><span class="sha user-select-contain" style="font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; color: rgb(68, 77, 86); white-space: nowrap;">макросы симуляции, обработки и отрисовки прикреплены. Функциональность была доавлена в класс ERDecayEXP1803. Добавленные методы основаны на </span>http://old.pskgu.ru/ebooks/djekson/djekson_12.pdf .</p> ERDecayEXP1803.</span></p>
<p style="background-color: rgb(255, 255, 255);"><img alt="" src="http://er.jinr.ru/develop/system/rich/rich_files/rich_files/000/000/059/original/LISE%20output1.png" style="vertical-align: middle; max-width: 100%;" /></p>
<p style="background-color: rgb(255, 255, 255);">В симуляции использовалось распределение энергии возбуждения 5H состоящее из двух гауссов: "mean"="0.4" "shape"="gauss" "strength"="1" "width"="0.02355"<br />
"mean"="1.2" "shape"="gauss" "strength"="1" "width"="0.2355". Значение, разыгрывающееся случайным образом из описанного распределения использовалось как добавка к массе оснровного состояния 5H, равного 4.69036244; // [GeV]. масса основного состояния задается пользователем в макросе симуляции. Исходя из заданного распределения, понятны раздвоения линий.</p>
<p style="background-color: rgb(255, 255, 255);">При разыгрывания бинарной реакции, в случае, если раыгранная энергия возбуждения превышает порог реакции, то энергия возбуждения разыгрывается ещё раз, пока бинарная реакция не сможет произойти. </p>
<p style="background-color: rgb(255, 255, 255);">Для разыгрывания бинарной реакции (6he + h2 - > 5H + 3He) использовалось распределение полярного угла вылета продуктов реакции в с.ц.м реакции. Распределение считывалось из прикрепленного файла.</p>
<p style="background-color: rgb(255, 255, 255);"><img alt="" src="http://er.jinr.ru/develop/system/rich/rich_files/rich_files/000/000/060/original/LISE%20output2.png" style="vertical-align: middle; max-width: 100%;" /></p>
<p style="background-color: rgb(255, 255, 255);">Рассматривались события, соответствующие второму максимуму углового распределения. (зелёные линии). Как и в рассчётах SIMONE, в данном случае, He3 в лабе вылетает под углами (21,33) в лабе.</p>
<p style="background-color: rgb(255, 255, 255);">Все результаты совпадают с полученными в SIMONE. </p>
<p style="background-color: rgb(255, 255, 255);">Используется нестандартный для er метод разыгрывания AD для бинарной реакции, не используются методы класса ERTextDecay. Цитата Виталика по этому поводу: "<span style="background-color: rgb(255, 255, 255);">Относительно ERTextDecay - давайте чтение из текстового файла применять, только тогда, когда не можем запрограммировать быстро и вычислительно просто результат распада. В этом случае получилось. Иван молодец.(шутка)</span>" </p>
<h2 style="padding-top: 2px; padding-left: 0px; background-color: rgb(255, 255, 255);"><span style="font-size: 14px;">Существующие проблемы</span></h2>
<p style="background-color: rgb(255, 255, 255);"><span style="background-color: rgb(255, 255, 255);">Возникла проблема, скорее всего связанная с использованием генератора случайных чисел, при разыгрывании угла вылета продуктов бинарной реакции в сцм. Double_t thetaCM = fADFunction->GetRandom(1.,150.)*TMath::DegToRad(); Было выяснено, что иногда (при большой статистике > 10000 событий) возникают события, в которых кинетическая энергия 6He равна 0, поэтому симуляция зацикливается на этапе разыгрывания бинарной реации. Задача перенесена </span>http://er.jinr.ru/develop/issues/169</p>