EXPERT ROOT

Комментарий к первому коммиту

1. root Файл геометрии мы заливаем в репозиторий. Да, он бинарный, но все таки очень лёгкий, поэтому не страшно.
2. В файле create_BeamDet_geo_v1.С параметры геометрии (такие как Int_t platesCount = 2; Double_t plateOffsetZ = 1500;) стоит вынести наверх макроса. Фактически обычно именно они изменяются от версии к версии.
   В нуле всегда находится мишень. Все части BeamDet должны быть расположены до мишени. У тебя полчучилось, что MWPC и второй TOF после нуля. Надо глобально сдвинуть назад весь детектор.
   Лучше завязать все на 4 переменные: глобальное расположение по Z двух TOF и двух MWPC, а уже служебные параметры: что как соориентировано внутри Assembly считать из них
3. В классе ERBeamDet. Сейчас в репозитории творится путаница с тем как надо оформлять класс симуляции. Тут я виноват. Правильно так как сейчас написано в документации. То есть текущее состояние поинта хранят переменные класса. Статическими переменными мы не пользуемся. У тебя возникли две переменные eloss и fEloss, которые по смыслу - одно и тоже. Сделай это одной переменной и у тебя будут нормально писать поинты.
4. Увидел особенность, связанную с файлом gconfig/g4Config.C. Мы используем многослойный пирог из библиотек, где каждый норовит добавить что-то свое. Тут долго объяснять пока просто прими на веру, что создание класса должно выглядеть так:
   TG4RunConfiguration* runConfiguration
   = new TG4RunConfiguration("geomRoot", "QBBC", "stepLimiter+stackPopper");

1. Следует также создать не один, а два класса поинтов: ERBeamDetTofPoint и ERBeamDetMwpcPoint. И заполнять соответственно две коллекции. Пример можно посмотреть в ERGadast

2. В поинтах MWPC записывать в какую именно станцию и объём проволоки попал поинт. В Tof в какой по счету пластик. Как доставть номера текущих объёмов при транспорте можно смотерть в ERNeuRad

   Int_t curVolId =  gMC->CurrentVolID(fiberInModuleNb); - номер копии текущего объема
   Int_t corOffVolId = gMC->CurrentVolOffID(1, moduleNb); - номер копии материнского объема

Добавить мишень в симуляцию

1) Пробовал разные параметры в TG4RunConfiguration* runConfiguration = new TG4RunConfiguration("geomRoot", "QBBC", "stepLimiter+stackPopper"). Не понимаю, правильно ли сейчас работает симуляция?

2) Если мы убрали статические переменные из processHits(), то, может, следует преписать метод AddPoint() без параметров?

3) В документации сейчас в описании реализации закона Биркса статические переменные в параметрах AddPoint остались.

1) О правильности симуляции можно судить по двум вещам: eventDisplay.C - визуально и по поинтам в sim.root

На eventDisplay у тебя сейчас такая ситуация:

И ~40тысяч вторичных частиц.То есть ион отдает свою энергию на ионизацию воздуха. Потому что материал пещеры - Air. Мы забыли упомянуть, что все это должно происходить в вакуумной трубе. Пока для того, чтобы ее не вводить в геометрию давай глобально поменяем материал пещеры (geometry/cave.geo) на vacuum (проверить как правильно пишется в geometry/media.geo). После такого изменения видим следующую картину:

Тут видно, что он долетел до второго Tof и до MWPC. И родил уже в них вторичные частцы.

В gconfig/g4config.in нужно установить каты на рождение вторичных частиц:

/mcPhysics/rangeCutForGamma  1. mm
/mcPhysics/rangeCutForElectron 1. mm

В sim.root можно видеть такое энерговыделение на всех поинта(используй veiew->Editor для лог шкалы):

С помощью условия на атрибут fType, который ты ввел, можно разделить на eloss в Tof и MWPC. Нужно правой кнопкой на дереве Start Tree Viewer. Перетащить eloss мышкой на X, fType на Experssion. Потом Edit experssion и там c++ синтаксис предикатов.

Для симуляции нужно доправить: 1) материал пещеры, 2) Разбить на две ветки поинты 3) Теперь действительно мржно убрать аргументы из AddPoint(). 4) Лучше сделать AddTofPoint() и AddMWPCPoint() 5) Добавить водородную мишень

В остальном отлично! В документации сейчас поправлю, спасибо!

Разбил поинты на две разные ветки.

В Gadast и детекторе из документации несколько отличается инициализация в кострукторах. Является ли один из вариантов более приоритетным?

Метод CopyClone() в Gadast прописан только для ERGadastCsIPoint. Сделал также для ERBeamDetTOFPoint.

Да пусть CopyClones работает пока только для одного типа поинтов. Это не страшно. В другомслучае и просто пришлось бы наследовать от одного предка, который бы мы сдлеали только для этого. Инициализация пусть будет пока так, как в документации.

А зачем fStilbenNr в MWPC?

Тоже не нужен. Удаляй

Что такое stackPopper TG4RunConfiguration()?

Можно ли задвать разные каты для частиц, в зависимости от материала? 

Mikhail Kozlov писал(а):

Что такое stackPopper TG4RunConfiguration()?

В симуляциях мы не обращаемся к Geant4 напрямую, а используем некоторую прослойку абстрагирующую нас от конкретной библиотеки транспорта - Virtual Monte Cerlo - https://root.cern.ch/vmc.Идея в том, что в принципе все библиотеки транспорта работают примерно одинаково и можно иметь один и тот же интерфес к ним всем и без проблем переключаться.  Раньше это было более актуально, потому что в колаборациях одновременно использовали geant3(https://en.wikipedia.org/wiki/GEANT-3) и geant4(http://geant4.cern.ch/).Сейчас все практически используют geant4. Вvmc также есть поддержка fluka - http://www.fluka.org/fluka.php (в основном применяется для расчетов радиационной защиты) и http://garfield.web.cern.ch/garfield/(в основном для газовых детекторов).

TG4RunConfiguration это класс VMC для настройки Geant4. - http://ivana.home.cern.ch/ivana/g4vmc_html/classTG4RunConfiguration.html

Его параметры: тип используемой геометрии, Physic List - список используемых процессов(https://root.cern.ch/physics-list) и специальные процессы VMC. Под специальными процессами подразумевается то, что VMC сам может вторгаться в логику транспорта. К примеру SpecialCuts означает, что будут применены каты в стиле VMC. Они энергетические и задаются в макросе SetCuts с помощью объекта gMC. Если этого специального процесса нет, то работают каты по умолчанию или используется нативный способ задания катов библиотеки транспорта. Для нас сейчас удобен именно этот вариант. Сейчас стоит задавать каты с помощью g4config.in это макрос настройки geant4. В Geant4 введен свой способ задания катов - range cut(http://geant4.cern.ch/G4UsersDocuments/UsersGuides/ForApplicationDeveloper/html/TrackingAndPhysics/thresholdVScut.html). Это кат на пробег новой частицы в данном материале, То есть если она по расчету пробежит меньше ката, то не надо заниматься ее трекингом. ее энергия уйдет в edep в точке, где разыгрывается процесс.

Строки для задания катов:

/mcPhysics/rangeCuts 1 mm
/mcPhysics/rangeCutForGamma  1. mm
/mcPhysics/rangeCutForElectron 1 mm

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

Михаилу актуально задавать разные каты в разных материалах. В ТОФ надо пониже, а в газе MWPC и так ОК.

Для получения реальной координаты пучка на мишени, ее следует добавить как чувствительный объем в этот детектор (то есть прямо в геометрию). Ввести еще один тип поинта. Мишень как пассивный отдельный элемент убрать из симуляции.

Примитив мишени оставить Tube или как-то по-другому сделать?

Ровно также. Только она теперь в детекторе. Просто, чтобы не плодить отдельный класс мишени и  как то странно будет его несаледовать от ERDetector

Тест рассылки уведомлений на почту

Необходимо добавить в симуляцию запись такого объекта как fBeamDetMCProjectile класса ERBeamDetParticle. Это состояние иона на мишени в симуляции. Имя иона мы передаем в класс симуляции через интерфейс. С помощью этого при трекинге через мишень отбираем трек иона. Считываем его параметры при входе в мишень и строим объект  fBeamDetMCProjectile.

Также в симуляции необходимо сохранить объект fBeamDetMCTrack класс ERBeamDetTrack. В нем необходимо запомнить теже данные, что и в реконструкции. Только в данном случае мы по честному получаем координаты на мишени и на MWPC.