EXPERT ROOT

Генератор распада 6Be --> alpha + p + p имеет форму текстогово файла. В прикрепленном  файле 6be_chex_noal_0_v2.dat пример 10 событий. В типичном файле находится 2 000 000 событий и его размер около 230 МБ. В файле находится 11 столбцов:

1. - енергия распада Е_Т; эту величину можно воспринимать также как энергию возбуждения 6Be и масса ядра тогда определяется как m_alpha = m_alpha + 2m_p + E_T

2. - 4. компоненты импулься частицы альфа - p_x, p_y, p_z в системе центра масс 6Be

5. - 7. компоненты импулься одного из протонов - p_x, p_y, p_z в системе центра масс 6Be

8. - 10. компоненты импулься второво из протонов - p_x, p_y, p_z в системе центра масс 6Be

11. некоторый параметр, который говорит про отношение состояния 0⁺ и 2⁺ в конктретном событию. Его стоит оставить в данных.

Чтение генератора можно найти в функции BeWork::FillSimKinFile(const char* generator, Int_t INPUTS, Double_t beamtheta, Double_t beamphi, const char* outputfile, Option_t *opt) {}.

Что надо сделать со значениями из генератора (переходы между системами координат) (неуклюжим образом) описано в функции BeWork::MonteCarloState(...).

Значит, задача стоит следующей: первычный пучок ⁶Li с энергией (ее неточностью) и размером, которые можем задавать, влетает в мишень и взаймодействует с ней. В каком-то месте в мишени (а именно в газе) происходит реакция p(⁶Li,⁶Be)n. Угловое распределение реакции можно разыграть по фазовому объему. Далее в детекторы влетает нейтрон. ⁶Be сразу распадается (как описано в генераторе) и в наш сетап летят альфа и два протоны.