EXP1803

На данный момент в методах класса ERDecayEXP1803 для бинарной реакции используются методы TGenPhaseSpace, в которой заложено, что угловое распределение вылетающих частиц - изотропно. 

Добавить методы, позволяющие разыгрывать реакцию, использую заданную пользователем функцию углового распределения.

Результат:¶

Вся функциональность проверялась в ветке 168_EventsWith6HeZeroSpeed

commit 2a576fbe7a34efde6a5cdf6ab3177af99fedf40f

макросы симуляции, обработки и  отрисовки прикреплены. Функциональность была доавлена в класс ERDecayEXP1803. Добавленные методы основаны на http://old.pskgu.ru/ebooks/djekson/djekson_12.pdf .



В симуляции использовалось распределение энергии возбуждения 5H состоящее из двух гауссов: "mean"="0.4" "shape"="gauss" "strength"="1" "width"="0.02355"
"mean"="1.2" "shape"="gauss" "strength"="1" "width"="0.2355". Значение, разыгрывающееся случайным образом из описанного распределения использовалось как добавка к массе оснровного состояния 5H, равного 4.69036244;  // [GeV]. масса основного состояния задается пользователем в макросе симуляции. Исходя из заданного распределения, понятны раздвоения линий.

При разыгрывания бинарной реакции, в случае, если раыгранная энергия возбуждения превышает порог реакции, то энергия возбуждения разыгрывается ещё раз, пока бинарная реакция не сможет произойти. 

Для разыгрывания бинарной реакции (6he + h2 - > 5H + 3He)  использовалось распределение полярного угла вылета продуктов реакции в с.ц.м реакции. Распределение считывалось из прикрепленного файла.

Рассматривались события, соответствующие второму максимуму углового распределения.  (зелёные линии). Как и в рассчётах SIMONE, в данном случае, He3 в лабе вылетает под углами (21,33) в лабе.

Все результаты совпадают с полученными в SIMONE. 

Используется нестандартный для er метод разыгрывания AD для бинарной реакции, не используются методы класса ERTextDecay. Цитата Виталика по этому поводу: "Относительно ERTextDecay - давайте чтение из текстового файла применять, только тогда, когда не можем запрограммировать быстро и вычислительно просто результат распада. В этом случае получилось." 

Существующие проблемы¶

Возникла проблема, скорее всего связанная с  использованием генератора случайных чисел, при разыгрывании угла вылета продуктов бинарной реакции в сцм. Double_t thetaCM = fADFunction->GetRandom(1.,150.)*TMath::DegToRad(); Было выяснено, что иногда (при большой статистике > 10000 событий) возникают события, в которых кинетическая энергия 6He равна 0, поэтому симуляция зацикливается на этапе разыгрывания бинарной реации.  Задача перенесена http://er.jinr.ru/develop/issues/169