Developing #147
Идентификация частицы в QTelescope
| Status: | Закрыта | Start date: | 02/26/2018 | |
|---|---|---|---|---|
| Priority: | Низкий | Due date: | ||
| Assignee: | Vitaliy Schetinin | % Done: | 100% | |
| Category: | Software | |||
| Target version: | - |
Description
Под идентификацией частицы подразумевается востановление ее 4-х веткора на мишени.
Из ERQTelescopeTrackFinder знаем направление импульса частицы.
Как считаем модуль импульса?
Коллеги, прошу дописать данную задачу!
History
#1 Updated by Sergey Belogurov over 6 years ago
Для нахождения модуля импульса надо знать массу частицы, в нашем случае He-3, и ее кинетическую энергию. Кинетическая энергия - это сумма энерговыделений в двух слоях телескопа плюс средние потери частицы с такой скоорстью и с таким Z в мертвых слоях.
Идентификация частицы - это предеелние ее A,Z и как слдествие массы. В основе идентификации анализ результатаов дижитизации.
DeltaE - энерговыделение в тонком слое, E - сумма в толстом и тонком слоях . графический кат на диаграмму DeltaE-E должен каким-то образом задать пользователь. Остается вопрос - где и как задать этот кат? В дополнение к этому кату могут быть еще и обсужденные ранее каты на преодоление порогов и малое отличие энерговыделений, считанных с двух сторон.
Получается, на основе грубой E и DeltaE мы идентифицируем частицу, а дальше, уже зная массу, вычисляем коррекцию на мертвые слои.
#2 Updated by Sergey Belogurov over 6 years ago
- Assignee changed from Sergey Belogurov to Vitaliy Schetinin
#3 Updated by Vratislav Chudoba over 6 years ago
- Category set to Software
Если имеется информация о 4импульсе пучковой частицы в центре мишени (эта задача уже практически решена) и, допустим, нормированный вектор частицы попавшей в телескоп (задача #146), тогда можно сделать графический кат идентификации на основе картинок из макроса анализа (это макрос, которым выводим все возможные картинки из дерева с помощью TTree::Draw()) и реконструкцию частицы в том же макросе. Предлагаю эту задачу пока оставить на этап обработки данных, который покажет, что надо и какие подходы имеются. Считаю, что продумать все очень хорошо заранее может не получится.
#4 Updated by Vitaliy Schetinin over 6 years ago
- % Done changed from 0 to 90
В соответвующую ветку влита данная функциональность. Добавлен класс ERQTelescopePID со следующим интерфесом
ERQTelescopePID* qtelescopePID = new ERQTelescopePID(verbose);
qtelescopePID->SetStationParticle("DoubleSi_SD1_XY_0",1000020030);
qtelescopePID->SetStationParticle("DoubleSi_SD2_XY_1",1000020030);
qtelescopePID->SetStationParticle("DoubleSi_D1_XY_2",1000020030);
run->AddTask(qtelescopePID);
SetStationParticle("DoubleSi_SD1_XY_0",1000020030); - задает гипотезу о том, что частица с pdg 1000020030 прилетела в телескоп DoubleSi_SD1_XY_0
В результате данной задачи появляется ветки с названиями вида: ERQTelescope_DoubleSi_SD1_XY_0_1000020030 с объектами типа ERQTelescopeParticle, которые имеют атрибуты:
fDeadEloss - потери данной частицы в мертвых слоях телескопа и половине мишени
fLVTarget - четырех вектор частицы на мишени
fLVTelescope - четырех вектор частицы на телескопе.
Гипотез про типы прилетающих частиц в телескоп может быть произвольное количество. Под каждую гипотезу будет создана своя ветка.
#5 Updated by Vitaliy Schetinin over 6 years ago
Добавлена возможность добавлять пользовательский кат на дерево входного файла (в том числе и графический) . Только события которые прошли данный кат будут обработаны. Для все остальных будут записаны пустые коллеции.
Интерфейс такой:
qtelescopePID->SetUserCut("ERQTelescopeSiDigi_DoubleSi_D1_XY_2_X.fEdep > 0.00001");
Где "ERQTelescopeSiDigi_DoubleSi_D1_XY_2_X.fEdep > 0.00001" объект класса TCut.
#6 Updated by Vitaliy Schetinin over 6 years ago
Ветка влита в dev
#7 Updated by Vitaliy Schetinin over 6 years ago
- Status changed from Открыта to Закрыта
- % Done changed from 90 to 100