NeuRad¶

NeuRad neutron detector¶

Разрабатываемый детектор нейтронов NeuRad (аббревиатура от Neutron Radioactivity) предназначен для получения прецизионной информации об угловых корреляциях распадных нейтронов с заряженным фрагментом. Эта информация используется для определения энергии распада, происходящего вблизи мишени, и времени жизни экзотического ядра. Информация о направлении вылета тяжелого фрагмента поступает от кремниевого микрострипового детектора, а направления полета нейтронов должны быть получены с помощью специально изготовленного детектора NeuRad, обладающего высоким угловым разрешением. Необходимое угловое разрешение для нейтронов должно быть ~0,1-0.2 мрад. Покрываемый угловой диапазон ±6 мрад. Такой малый угловой диапазон отражает малый поперечный импульс, отвечающий энергии распада 0.1 кэВ – 100 кэВ.

Нейтрон - частица, не имеющая заряда и не способная ионизировать вещество. Поэтому нейтрон регистрируют косвенно. При попадании нейтрона в объём сцинтиллятора, происходит ядерная реакция выбивания заряженной частицы, в большинстве случаев протона, называемого протоном отдачи. Заряженная частица ионизирует сцинтиллятор, что приводит к появлению вспышки света. Порог чувствительности детектора к протонам отдачи должен быть не выше 1000 кэВ. Процедура реконструкции события должна позволять отличить взаимодействие двух распадных нейтронов от последовательных взаимодействий одного нейтрона и от фоновых взаимодействий гамма-излучения.

Детекторная система NeuRad имеет модульную структуру. Каждый модуль состоит из сцинтилляционных оптических волокон квадратного сечения 3мм*3мм и системы считывания. Волокна собраны в пучок квадратного сечения, состоящий из 256 штук, и ориентированы вдоль пучка ионов, падающих на мишень. Прототип пучка волокон, состоящий из 64 штук, показан на рис.2.1. совместно с технологической оснасткой, использующейся при его изготовлении.

Рисунок 1. Пучок оптических волокон в специально разработанной технологической оправке

Для оптического считывания пучка волокон используются многоанодные фотоэлектронные умножители (МА ФЭУ), расположенные с двух сторон. Фотокатод МА ФЭУ разбит на 64 пикселя, к каждому из которых подходит четыре волокна. Фотоэлектроны из каждого пикселя попадают в свою секцию динодной системы, где развивается лавина, приходящая, в конце концов, на соответствующий анод. Имеется возможность считывания с последнего динода сигнала, пропорционального сумме сигналов во всех анодах.

Принцип работы модуля упрощенно проиллюстрирован на рис. 2. Двустороннее считывание позволяет определить продольную координату места взаимодействия по разности задержек в приходе сигнала. В поперечном направлении точка прилета нейтрона определяется центром пикселя, в котором было зарегистрировано первое по времени взаимодействие.

Рисунок 2. Принцип действия NeuRad

Simulation¶

В задачи класса симуляции входит:

Транспорт всех частиц в событии через объём сцинтиллятора и разыгрывание происходящих с ними физических процессов, включающих ядерные и электромагнитные взаимодействия, в частности, потери энергии на ионизацию.
Преобразование энерговыделения на каждом шаге (dE_i, [MeV]) в величину той же размерности, но пропорциональную количеству рожденного света (dL_i, [MeV],), которое зависит не только от величины энерговыделения, но и от плотности ионизации, т.е. размера области пространства, в которой данное энерговыделение сосредоточено. Частицы с высокой плотностью ионизации (протоны, альфа-частицы, тяжёлые ионы, осколки деления) дают в большинстве сцинтилляторов меньшее количество фотонов, чем гамма-кванты, электроны или мюоны, потерявшие такую же энергию.

Геометрия¶

Модель состоит из 49 модулей, расположенных квадратным массивом 7x7. Каждый модуль состоит из 64 волокон, также собранных в квадратный массив 8x8. К модулю с каждой стороны подключен многоанодный ФЭУ, каждое волокно с каждой стороны просматривается одним пикселем. Волокно имеет длину 1 м и квадратное сечение 6*6 мм2. Таким образом, одно волокно в модели соответствует четырем волокнам в реальном детекторе. Конец модуля, обращенный к мишени, т.е. расположенный выше по пучку будем называть фронтальным, а противоположный конец - задним.

Материалы¶

Материал сцинтиллятора задан как стехиометрический эквивалент полистирола. Наличие МА ФЭУ на торцах прибора моделируется прямоугольными пластинами из алюминия толщиной 2.4 мм.

BC408               2  1.00794  12.0107  1.  6.  1.032  0.0764  0.916
                    1  1  20.  .001
                    0
aluminium           1  26.98 13. 2.7
                    0  1  20.  .001
                    0
vacuum              1  1.e-16  1.e-16  1.e-16
                    0  1  20.  .001
                    0

NeuRad.v1.geo.root¶

Геометрическая модель одного модуля детектора NeuRad.

Модель имеет следующее дерево объёмов:

cave
  NeuRad
    module
      fiber
        cladding 1
        cladding 2

NeuRad.v2.geo.root¶

Геометрическая модель всего детектора NeuRad.

Закон Биркса¶

Кроме стандартной характеристики Energy Loss, на этапе симуляции вычисляется Light Yield (световыход), величина которого определяется по закону Биркса. Вынос этих вычислений связан с необходимостью учёта неравномерности световыхода по длине отрезка энерговыдления. Параметризация закона Биркса взята из Geant4. Для имплементации закона Биркса введены следующие константы:

\[ \begin{align}\begin{aligned}dP = 1.032\\BirkC1 = 0.013/dP\\BirkC2 = \frac{9.6*10^{-6}}{dP^2}\end{aligned}\end{align} \]

dP - имеет смысл плотности вещесва. Световыход на данном шаге вычисляется так:

\[dL = \frac{dE}{1+BorkC1*\frac{dE}{dx} + BirkC2*\frac{dE}{dx}^2}\]

dL- величина светоыхода на данном шаге транспорта, dE - величина депозита энергии на данном шаге, dx - длина шага.

Предусмотрено изменение константы Биркса для случая высоко заряженных частиц: если энерговыделение оставляет частицца с зарядом больше 2 GeV, то

\[BirkC1 = BirkC1 * 7.2/12.6\]

Световыход поинта определяется как сумма световыходов на всех шагах VirtualMC, формирующих поинт

\[ \begin{align}\begin{aligned}E = \sum dE\ [GeV]\\L = \sum dL\ [GeV]\end{aligned}\end{align} \]

Для более точного расчёта времени прихода фотонов на катод, предусмотрен порог на длину поинта - 4 см. Когда накопленная длина поинта превысила порог, начинает формироваться новый поинт.

Реализация¶

Реализация в классах: ERNeuRad, ERNeuRadPoint, ERNeuRadStep

Алгоритм симуляции проведения частиц через детектор NeuRad реализован в методе ERNeuRad::ProcessHits(FairVolume* vol). Данный метод вызывается на каждом шаге работы VirtualMC, если шаг находится внутри активного объема детектора NeuRad. Активные объемы задаются в методе ERNeuRad::CheckIfSensetive(std::string name).

Результатом работы всех вызывов ERNeuRad::ProcessHits является коллекция объектов класса ERNeuRadPoint и, если указано, ERNeuRadStep. Поинт - прямолинейный отрезок энерговыделения в одном объёме, депозит энергии которого привысил порог. Порог на депозит энергии поинта в NeuRad равен нулю.

ERNeuRadStep - класс для хранения информации об одном шаге VirtualMC в NeuRad.

Добавление в макрос¶

Минимальный код добавления в макрос:

ERNeuRad* neuRad= new ERNeuRad("ERNeuRad", kTRUE,0);
neuRad->SetGeometryFileName("NeuRad.v2.geo.root");
run->AddModule(neuRad);

Третий параметр в конструкторе - verbosity - регулирует объём логов:

0 - только количество поинтов в событии
1 - информация о каждом поинте
2 - информация о каждом шаге Geant

Для анализа эффективности и работы ERNeuRadMatcher понадобилось хранить информацию о первом шаге Geant в NeuRad:

ERNeuRad* neuRad= new ERNeuRad("ERNeuRad", kTRUE,0);
neuRad->SetGeometryFileName("NeuRad.v2.geo.root");
neuRad->SetStorePrimarySteps();
run->AddModule(neuRad);

Для анализа процессов понадобилась возможность хранить информацию о всех шагах:

ERNeuRad* neuRad= new ERNeuRad("ERNeuRad", kTRUE,0);
neuRad->SetGeometryFileName("NeuRad.v2.geo.root");
neuRad->SetStoreAllSteps();
run->AddModule(neuRad);

Результаты симуляции¶

Результатом симуляции в общем случае является дерево со следующей структурой:

MCHeader

fNeuRadELoss - суммарный ELoss в событии

fNeuRadLY - суммарный LY в событии

NeuRadPoint

fEventID - ID события

fTrackID - ID трека

fMot0TrackID - ID материнского трека

fFiberInBundleNb - номер файбера в модуле

fBundleNb - номер модуля

fMass - масса частицы

fX, fY, fZ - координаты начала поинта

fX_out, fY_out, fZ_out - координаты конца поинта

fPx, fPy, fPz - импульс в начале поинта

fPx_out, fPy_out, fPz_out - импульс в конце поинта

fTimeIn, fTimeOut - время в начале и конце формирования поинта

fLightYield - световыход в поинте

fPID - PID частицы

fCharge - заряд частицы

fTrackLength - длина трека до начала поинта

fELoss - энерговыдление в поинте

fDetectorID - ID детектора из DetectorList

NeuRadStep

fEventID - ID события

fStepNr - номер шага в рамках поинта

fTrackID - ID трека

fMot0TrackID - ID материнского трека

fFiberInBundleNb - номер волокна в модуле

fX, fY, fZ - координаты шага

fPx, fPy, fPz - импульс на шаге

fTOF - время на шаге

fLength - длина шага ????

fPID - PID частицы

fMass - масса частицы

fTrackStatus - статус трека на шаге (начался, окончился, продолжает транспорт)

fEloss - энерговыделение на шаге

fCharge - заряд частицы

fProcessID - ID процесса из Geant4

Digitization¶

Задачей диджитизации является формирование формы сигнала на выходе фотокатода. Интегривание сигнала по времени вынесено в процедуры поиска хитов. Это связано с тем, что на данный момент система считвания не устоялась.

Моделируются следующие процессы:

Рождение фотонов:¶

\[Np = L * C_{SciFi}\]

где \(C_{SciFi}\) – световыход сцинтиллятора, т.е. количество фотонов, рожденных на один МэВ энергопотерь электрона. Для сцинтилляционных волокон Saint-Gobain BCF -12, используемых в NeuRad, \(C_{SciFi}=8000\ [\frac{photons}{MeV}]\).

Две причины затухания¶

Учитываются две причины зависимости количества фотонов, достигших фотокатода, от продольной координаты точки взаимодействия: телесный угол, под которым виден фотокатод из точки взаимодействия с учетом полного внутреннего отражения – быстро спадающая экспонента с характерной длиной 0.5 см, и поглощение света в материале волокна – экспонента с характерной длиной 2м. Соответствующие выражения имеют вид:

\[ \begin{align}\begin{aligned}N_p^f= N_p(k_1exp⁡(-Z_{fiber}/0.5)+k_2exp⁡(-Z_{fiber}/200));\\N_p^b= N_p(k_1exp⁡(-(l-Z_{fiber})/0.5)+k_2exp⁡(-(l-Z_{fiber}/200));\\k_2=C_{lftr};\\k_1=0.5-C_{lftr};\end{aligned}\end{align} \]

где \(C_{lftr}\)- константа захвата света в угол полного внутреннего отражения света (light fraction in total reflection), Z_fiber [cm]- расстояние от фотокатода до места энергвыделения, \(N_p^f,N_p^b\) – количество фотонов пришедших на фронтальный и задний пиксели соответственно.

Квантовая эффективность¶

Рассчитывается, учитывая квантовую эффективность, среднее число родившихся на каждом пикселе фотоэлектронов \(\hat{N_{pe}}\) для данного энерговыделения:

\[N_{pe} = N_p*C_{qe},\]

где \(C_{qe}\) –квантовая эффективность фотокатода. \(C_{qe}\) берётся из таблицы квантовых эффективностей по индексу волокна и модуля (\(C_{qe} \approx 20%\)).

Количество родившихся в событии фотоэлектронов разыгрывается по закону Пуассона со средним \(\hat{N_{pe}}\):

\[N_{pe}= P(\hat{N_{pe}}),\]

где P(∙)-функция распределения Пуассона.

Экспоненциальность высвечивания¶

Высвечивание фотонов происходит по экспоненциальному закону с характерным временем \(\tau_{sci}=3.2 нс\):

\[N_p(t)=1-exp⁡(\frac{-t}{\tau_{sci}}),\]

где \(N_p(t)\) – количество фотонов, высвеченных за время t. Для моделирования каждого фотоэлектрона необходимо знать время прихода соответствующего фотона на катод. Для этого найдем функцию обратную (2.6) и подадим ей в качестве аргумента равномерное распределение.

\[T_k=T_{edep}+(-1)*\tau_{sci}*log⁡(1-U[0:1])+Z_{fiber}/С_{sc},\]

где \(T_{edep}\) – время соответствующего энерговыделения, второе слагаемое – время до высвечивания данного фотона, определяемое обратным экспоненциальным распределением, U[0:1] – равномерное распределение, \(С_{sс}\) – скорость света в материале сцинтиллятора (n=1.58).

Динодная система¶

Амплитуда одноэлектронного сигнала \(A_{pe}\), рассчитывается как:

\[A_{pe}=\mid N(A,\sigma) \mid,\]

где \(N(A,\sigma)\) – функция нормального распределения c подобранными параметрами A и \(\sigma\). Cреднее значение распределения амплитуд \(A_{pe}\) нормировано на единицу для возможности выражения всех порогов в фотоэлектронах. Формула моделирует распределение, показанное на рис.3.

Рисунок 3. Распределение амплитуды одноэлектронного сигнала

Время прихода сигнала на анод рассчитывается как:

\[T_a=T_k+N(D_{PMT},J_{PMT}),\]

где \(D_{PMT}\) – задержка на динодной системе, \(J_{PMT}\) - флуктуация времени прохождения электронной лавины через динодную систему.

Форма сигнала¶

Получив параметры сигналов фотоэлектронов \(A_{pe}\ и\ T_a\), можно ввести функцию формы сигнала параметризованную на них. Была предложена следующая функция сигнала фотоэлектрона:

\[A(t)=T_a+ a*A_{pe}*T_aexp⁡(-T_a/b),\]

где a и b – коэффициенты для настраивания формы. Они определяются для ФЭУ из эксперимента. Характерный график такой зависимости в условных единицах без начального сдвига представлен на рис.

Рисунок 4. Форма одноэлектронного сигнала

Получив формы сигналов от каждого фотоэлектрона, можно получить суммарный токовый сигнал на аноде как:

\[Aa_j(t)=\sum_{i=1}^{N_{pe}} A_i(t)\]

где \(N_{pe}\) – количество фотоэлектронов, рожденных от света на одном пикселе, i – индекс фотоэлектрона, попавшего в динодную систему данного анода, j – индекс анода (волокна). Результат моделирования формы сигнала с двух сторон сцинтилляционного волокна представлен на рис.2.7.

Рисунок 5.Суммарные сигналы с двух сторон волокна

При расчёте сигнала с общего динода учитываются коэффициенты усиления:

\[Ad_j(t)=Aa_j(t)*Cg_j,\]

где \(Cg_j\)- коэффициент усиления динодной системы, который берётся из таблицы коэффициентов по индексу модуля и волокна.

Кросстолки¶

Также реализуется система кросстолков, моделирующая перетекание сигнала между анодами многоанодного ФЭУ.

Пример файла параметров¶

//-------------------------------------------------------------
//----------- Digitisation parameter for NeuRad    ------------
//-------------------------------------------------------------
[ERNeuRadDigiPar]
//-------------------------------------------------------------
ERNeuRadNofBundles: Int_t 49
ERNeuRadNofFibers: Int_t 64
ERNeuRadPMTQuantumEfficiency: Float_t \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 \
2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
ERNeuRadPMTGain: Float_t \ //[mV/p.e]
49947 4.78262 5.3909 4.98497 5.41213 4.97164 4.54956 4.96265 \
00396 4.79462 5.6956 4.50747 4.97553 4.27833 4.46966 4.30585 \
3837 4.63199 5.28986 4.80893 6.03045 4.38259 5.58267 4.7729 \
93261 4.75018 4.90878 5.92216 4.87858 5.99843 5.0024 4.78891 \
77027 5.04737 5.76235 5.60838 4.93185 4.90039 4.85312 4.94082 \
41985 5.00618 5.01103 4.96437 4.33919 5.87216 4.83381 4.87636 \
40989 4.40479 3.99512 5.66919 5.24409 5.14818 5.12196 4.45383 \
70045 4.60604 3.70993 4.39176 4.37836 4.88106 5.9002 4.76556
ERNeuRadPMTSigma: Float_t \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5. \
5. 5. 5. 5. 5. 5. 5.
ERNeuRadPMTCrosstalks: Float_t \
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.0 \
0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.1 0.1 0. 0.0 \
0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.

ERNeuRadNofModules и ERNeuRadNofFibers - количество молудей и волокон. ERNeuRadPMTQuantumEfficiency - матрица квантовых эффективностей. ERNeuRadPMTGain, ERNeuRadPMTSigma - матрица усилений и их сигм фотоумножителя. ERNeuRadPMTCrosstalks - матрица кросстолков.

Все перечисленные выше матрицы на данный момент применяются к каждому модулю. Поэтому количество элементов матрицы должно быть равно ERNeuRadNofFibers.

Матрица кросстолков для модуля состоит и матриц кросстолков для каждого волокна:

0. 0.
0. 0.1
0.1 0.1

Такая матрица означает, что вероятность перетекания фотоэлектрона к соседу справа, снизу и в правом нижнем углу равно 10%. Вероятность перетекания к остальным соседям равно нулю. В центре матрицы задан, условный ноль, чтобы пользователю не пришлось самостоятельно высчитывать оставшюся вероятность. Если волокно я вляется крайним вероятности перетекания к несуществующим соседям должны быть заданы нулём для общности.

Реализация¶

Реализация в классах: ERNeuRadDigitizer, ERNeuRadPMTSignal, ERNeuRadDigiPar, ERNeuRadSetup

ERNeuRadDigitizer - класс task процедуры диджитизации.

ERNeuRadPMTSignal - класс формы сигнала на выходе с фотоумножителя. Основной результат диджитизации.

Добавление в макрос¶

Int_t verbose = 1; // 1 - only standard log print, 2 - print digi information
ERNeuRadDigitizer* digitizer = new ERNeuRadDigitizer(verbose);
fRun->AddTask(digitizer);

Параметр в конструкторе - verbosity - регулирует объём логов:

0 - количество сигналов в событии
1 - информация о каждом сигнале

Reconstruction¶

Задача NeuRad – обнаружить точку входа нейтрона в детектор. Под реконструкцией события будем подразумевать идентификацию волокна, в которое влетел нейтрон. За координаты попадания нейтрона будем принимать координаты центра этого волокна. Задачу осложняет, то что ионизация, вызванная выбитыми из ядра сцинтиллятора частицами, происходит не только в объёме волокна, в который попал нейтрон, но и в соседних. В результате сигнал появляется на кластере пикселей.

Поэтому для NeuRad предусмотрены два режима работы:

Однофайберный режим. На каждый пиксель ставится порог на интеграл сигнала (\(T_p\)). Рассматриваются только те события, в которых этот порог с обоих концов преодолен только на одном из волокон. В этом режиме временная информация сигнала не учитывается. Порог вводится на интеграл импульса тока, т.е. заряд, как:

\[\begin{split}\begin{equation*} \begin{cases} Q_j^{front}=\int_{t}^{t+5[нс]}Aa_j^{front}(t)dt >T_p \\ Q_j^{back}=\int_{t}^{t+5[нс]}Aa_j^{back}(t)dt >T_p \end{cases} \end{equation*}\end{split}\]

где \(A_{PMT}^{front}(t),A_{PMT}^{back}(t)\) - сигналы на фронтальном и заднем пикселе соответственно, \(t \in [T_a,T_a+T_s], T_a=min⁡((T_a)_i),T_s\) – продолжительность сигнала, численно определяемая, как конечная точка сигнала последнего одноэлектрона. Порог \(T_p\) вводится в единицах фотоэлектронов, умноженных на интеграл сигнала одного фотоэлектрона.

Многофайберный режим. На каждый фронтальный пиксель ставится порог на мгновенный уровень сигнала. Момент времени, когда этот порог преодолен, служит временной отметкой (time stamp) данного пикселя. Отбираются события, в которых порог преодолен хотя бы на одном из фронтальных пикселей во всем детекторе. Затем отбираются события, в которых сумма интегралов динодных сигналов обоих МА ФЭУ какого-либо модуля преодолела порог \(T_M\). Считается, что нейтрон попал в волокно, фронтальный пиксель которого загорелся раньше остальных. Систему поставленных порогов в данном режиме можно описать так:

\[\begin{split}\begin{equation*} \begin{cases} Aa_j^{front}(t)>T_p \\ \sum_{j=1}^{M} \int_{t}^{t+5[нс]}Ad_j^{front}(t)dt+\sum_{j=1}^{M} \int_{t}^{t+5[нс]} Ad_j^{back}(t)dt >T_M, \end{cases} \end{equation*}\end{split}\]

где М – количество волокон в модуле, \(t \in [T_a,T_a+T_s]\), порог \(T_p\) в этом режиме вводится в единицах одноэлектронных амплитуд, порог \(T_M\) – в единицах фотоэлектронов, у множенных на интеграл сигнала одного фотоэлектрона. \(Ad_j\) посчитан по формуле 2.12. Результатом работы алгоритма реконструкции в детекторе NeuRad является хит – т.е предположительная точка, где нейтрон влетел в детектор.

В однофайберном режиме хитом служит центр единственного пикселя, где выполнено условие 2.12. В многофайберном режиме центр каждого переднего пикселя, в котором выполнено условие 2.13, является кандидатом в хиты, а хитом оказывается тот из них, в котором timestamp имеет наименьшее значение

Рисунок 6. Результаты работы однофайберного и многофйберного режима