Developing #147

Идентификация частицы в QTelescope

Added by Vitaliy Schetinin almost 7 years ago. Updated over 6 years ago.

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 almost 7 years ago

Для нахождения модуля импульса надо знать массу частицы, в нашем случае He-3, и ее кинетическую энергию.  Кинетическая энергия - это сумма энерговыделений в двух слоях телескопа плюс средние потери частицы с такой скоорстью и с таким Z в мертвых слоях.

Идентификация частицы - это предеелние ее A,Z и как слдествие массы. В основе идентификации анализ результатаов дижитизации.

DeltaE - энерговыделение в тонком слое, E - сумма в толстом и тонком слоях . графический кат на диаграмму DeltaE-E должен каким-то образом задать пользователь.  Остается вопрос - где и как задать этот кат?  В дополнение к этому кату могут быть еще и обсужденные ранее каты на преодоление порогов и малое отличие энерговыделений, считанных с двух сторон.   

Получается, на основе грубой E и DeltaE мы идентифицируем частицу, а дальше, уже зная массу, вычисляем коррекцию на мертвые слои.

#2 Updated by Sergey Belogurov almost 7 years ago

  • Assignee changed from Sergey Belogurov to Vitaliy Schetinin

#3 Updated by Vratislav Chudoba almost 7 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

Also available in: Atom PDF