bugs #132

Симуляция эксперимента по водороду-5

Added by Ivan Muzalevsky almost 7 years ago. Updated almost 7 years ago.

Status:ЗакрытаStart date:
Priority:НеотложныйDue date:
Assignee:Mikhail Kozlov% Done:

100%

Category:-
Target version:-

Description

1) при симуляции я вижу, что в больгшинстве событий 5h не рождается. По словам Вратислава, симуляция была настроена так, что рассматриваются только эти события. Другие нам в принципе не нужны. На данный момент, на примере файла, ссылку на гугл диске которого я тебе скинул, можно сказать, что тоьлко в 8-9 процентах событий я вижу интересный мне 5H. (из 10000 событий 5H я вижу в 847 событиях) 

2) Для проверки кинематики, я хочу получить 4 TLorentzVector'a (5h и его продукты распада: 3h,n1,n2). Так как в некоторых событиях нейтронов рождается больше, чем два, то я делал два метода отбора нужных мне нейтронов:   1 - сумма 4-векторов двух нейтронов с 4-вектором h3 должна равняться 4-вектору 5h.  2 - motherID у нейтронов должен быть такой же как у h3. Для обоих методов возникла проблема, которая заключается в том, что существуют события, в которых рождаются несколько h3. Поэтому возникает теперь вопрос, как выбрать нружный h3, котрый возник из распада h5? 

3) посылаю тебе макрос showData.C, на всяких случай, которым я читал свой файл full.root. а также output.txt - список pdg частиц, рождающихся в событиях, где я вижу несколько h3. Может тебе это как то поможет, не знаю.

showData.C Magnifier - reading full.root (3.66 KB) Ivan Muzalevsky, 02/12/2018 04:04 PM

output.txt Magnifier - events are separated by the string (next event) (951 Bytes) Ivan Muzalevsky, 02/12/2018 04:04 PM

sim_full.C.txt Magnifier (12.1 KB) Ivan Muzalevsky, 02/12/2018 11:10 PM

History

#1 Updated by Mikhail Kozlov almost 7 years ago

Мне бы еще макрос, который генерирует выходной файл, посмотреть. 
По-моему, пучок просто не попадает в мишень в большинстве случаев. Когда я запускаю свою симуляцию без BeamDet напрямую в мишень, у меня во всех событиях происходит реакция. С ним же происходит отклонение и реакция отсутсвует.

#2 Updated by Ivan Muzalevsky almost 7 years ago

Mikhail Kozlov wrote:

Мне бы еще макрос, который генерирует выходной файл, посмотреть. 
По-моему, пучок просто не попадает в мишень в большинстве случаев. Когда я запускаю свою симуляцию без BeamDet напрямую в мишень, у меня во всех событиях происходит реакция. С ним же происходит отклонение и реакция отсутсвует.

Макрос - sim_full.C в ветке, exp1803. я его не менял ни разу. ну скидываю тут на всякий случай.

#3 Updated by Vratislav Chudoba almost 7 years ago


У нас есть возможность посмотреть энергию пучковой частици в объеме мишени? У реакции есть порог и если энергия пучка ниже его, то реакция не может осуществиться. С другой стороны, если мы сильно не ложанулись в толщине слоев пучковых детекторов и окон мишени, тогда при 40 МэВ/А пучка, у нас должно оставаться достаточно энергии.

#4 Updated by Mikhail Kozlov almost 7 years ago

Вратислав, энергию пучка можно посмотреть путем ввода чувствительной мишени в BeamDet с помощью интерфейса: setupBeamDet->AddSensitiveTarget(). Притом пассивную мишень нужно убрать из симуляции. Я сейчас пропустил без отклонений пучок He6 40MeV/nucleon. После прохождения многопроволочных камер импульс пучка на входе в мишень 1.653 GeV.

#5 Updated by Ivan Muzalevsky almost 7 years ago

  • % Done changed from 0 to 70

Можно ли закрывать задачу? Проблему, которую я описывал при создании была решена.

#6 Updated by Ivan Muzalevsky almost 7 years ago

  • Status changed from Открыта to Закрыта
  • % Done changed from 0 to 100

Also available in: Atom PDF