EXP1803

<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, &quot;Segoe UI&quot;, Helvetica, Arial, sans-serif, &quot;Apple Color Emoji&quot;, &quot;Segoe UI Emoji&quot;, &quot;Segoe UI Symbol&quot;; 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, &quot;Liberation Mono&quot;, Menlo, Courier, monospace; white-space: nowrap;">commit </span><span class="sha user-select-contain" style="font-family: SFMono-Regular, Consolas, &quot;Liberation Mono&quot;, 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, &quot;Liberation Mono&quot;, Menlo, Courier, monospace; color: rgb(68, 77, 86); white-space: nowrap;">макросы симуляции, обработки и  отрисовки прикреплены. Функциональность была доавлена в класс ERDecayEXP1803. Добавленные методы основаны на </span>http://old.pskgu.ru/ebooks/djekson/djekson_12.pdf .</p>

<p style="background-color: rgb(255, 255, 255);"></p> 255);&quot;&gt;&lt;img alt=&quot;&quot; src=&quot;http://er.jinr.ru/develop/system/rich/rich_files/rich_files/000/000/059/original/LISE%20output1.png&quot; style=&quot;vertical-align: middle; max-width: 100%;&quot; /&gt;&lt;/p&gt;

<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> получилось. Иван молодец.(шутка)&lt;/span&gt;&quot; &lt;/p&gt;

<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>