bugs #168

Updated by Ivan Muzalevsky over 6 years ago

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

<p>На данный момент симуляция ломается на событиях, в которых я вижу, что кинетическая энергия 6He = 0. </p>

<p>191 строка ERDecayEXP1803.cxx : cout << lv6He(0) << " " << lv6He(1) << " " << lv6He(2)<< " " <<lv6he(3)> </lv6he(3)></p> &lt;&lt;lv6He(3) &lt;&lt; endl;  даёт первые три числа нули&lt;/p&gt;

<p>Данная проблема пропадает, если закомментировать в этом файле строку 306, и раскоментировать 307. (разыгрывать угловое распределение из изотропного распределения.) </p>

<p>Необходимо понять, откуда эти события с скоростью 6He равной 0 и для таких событий написать следующее условие, которое в любом случае понадобится при симуляции:</p>

<p>Если (ECM - f3He->Mass() - f5HMass)<0, (где ECM - энергия реакции в системе центра масс) то не нужно пытаться произвести бинарную реакцию по получению водорода 5. Может быть на данном этапе такие события просто не нужно писать а перехождить к следующему.</p>

<p>Возможно косяк просто изза моей криворукости, но я пока что его не нашёл.</p>

Back