EXP1803

В ветку BeamDetParticle добавить набор переменных, которые описывают реконструированную частицу в плоскости центра мишени. На данный момент удобно работать со всеми текущими переменными (fFourMomentum, fPx, fPy, fPz, fP, ...), новые переменные могут называться например fPxt, fPxT или fPx_t.

Текущие переменные описывают пучковую частицу (в некотором достаточном приближении) на выходе из второго пластика (ToF2). Эта информация нужна, но надо добавить листы с энергией после прохождения всего вещества на пути до центра мишени, причем направление останется таким же. Для каждого события надо посчитать с помощью калькулятора потер энергии потери в (типично) MWPC1, MPWC2, передних окнах ячейки мишени и половине толщины газа.

Новая энергия будет отличаться на 5 - 15 МэВ, что чрезмерно важно для реконструкции спектра изучаемого ядра.

Результат:¶

Задача в первом приближении работает. Симмуляция и реконструкция проводилась с помощью макросов в прикрепленных файлах в ветке "137_BeamDetEloss"

commit 08edfe61a72b707d153b6ba68dfaf22e6bbfd64c
Author: SchetininVitaliy
Date:   Wed Feb 28 10:15:29 2018 +0300

С помощью макроса в прикрепленном файле видим:

В первом столбце наблюдается восстановленная кинетическая энергия пучка после прохождения детекторов TOF. На остальных картинках видим кинетическую энергию пучка восстановленную в центральной плоскости мишени. Наличие трех пиков можно объяснить тем, что путь частицы реконструированный с помощью MWPC очень сильно поврежден многократным рассеянием на всех плоскостях и при провождению частистицы обратно в плоскость TOF1 траектория не проходит через все детекторы. Тем не менее можем легко понять, какой из пиков правильный (самый низкоэнергетичный).

Известные проблеммы:

1. Нужна количественная проверка на потери в детекторах TOF.
2. Надо посмотреть, как сильно влияет "Особым образом учитываются энергопотери в мишени. Так как мы не знаем точку реакии в реконструкции, берутся энергопотери, как если бы ион пролетел мишень полностью и деляться на 2." (ВЩ). Такой расчет в принципе не правильный, но может быть, что из-за точности разницу не заметим.
3. Надо разобраться с подавлением многократного рассеяния на детекторах TOF, чтобы очистить вторую панель вышеуказанной картинки.
4. Есть несогласие в массе ⁶He. В симмуляции и реконструкции стоит 5605,54 MэВ, согласно данным из nrv.jinr.ru должна быть масса 5606,556095 МэВ. В макросе симмуляции не найдена возможность управлять массей пучкового иона.