Анализ редкого распада K⁺ → e⁺νμ⁺μ⁻ в эксперименте NA62

Теоретические основы и мотивация изучения распада K⁺ → e⁺νμ⁺μ⁻

Добро пожаловать в курс, посвященный анализу одного из самых интригующих процессов в физике элементарных частиц — редкого распада заряженного каона. В этой первой статье мы заложим фундамент для понимания того, почему тысячи ученых в ЦЕРН (CERN) тратят годы на изучение событий, которые происходят с вероятностью менее одной миллиардной.

Мы разберем, что такое каоны, как Стандартная модель описывает их поведение, и почему именно канал распада с образованием четырех лептонов (позитрона, нейтрино и пары мюонов) так важен для современной науки.

1. Каоны и их место в микромире

Прежде чем говорить о распаде, нужно понять, что именно распадается. Каоны (или К-мезоны) — это частицы, которые сыграли ключевую роль в становлении физики частиц в середине XX века. Именно благодаря им мы узнали о существовании такого свойства материи, как «странность».

Структура K⁺-мезона

Каон — это адрон, то есть частица, состоящая из кварков. В частности, положительно заряженный каон () состоит из двух кварков:

!Кварковая структура положительного каона

С точки зрения кварковой модели, состав выглядит так:

Где:

— положительно заряженный каон.

— верхний (up) кварк.

— странный (strange) антикварк.

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

2. Стандартная модель и необходимость эффективных теорий

Стандартная модель (СМ) — это наша лучшая теория, описывающая фундаментальные частицы и их взаимодействия. Однако, когда дело доходит до описания адронов при низких энергиях (как в случае с распадом каона), математический аппарат СМ становится невероятно сложным.

Взаимодействие кварков описывается Квантовой Хромодинамикой (КХД). При высоких энергиях кварки ведут себя почти как свободные частицы, и мы можем легко считать их взаимодействия. Но при низких энергиях (масса каона всего около 493 МэВ/c²) взаимодействие становится настолько сильным, что методы теории возмущений перестают работать. Мы не можем просто сложить ряд диаграмм Фейнмана и получить точный ответ.

Киральная теория возмущений (ChPT)

Чтобы обойти эту проблему, физики используют Киральную теорию возмущений (Chiral Perturbation Theory, ChPT). Это эффективная теория поля.

> Эффективная теория похожа на карту местности: она не описывает каждый атом асфальта (как полная КХД), но отлично показывает дороги и маршруты (поведение мезонов), которые нам нужны для навигации.

В ChPT вместо кварков и глюонов мы работаем сразу с мезонами (пионами, каонами) как с фундаментальными степенями свободы. Это позволяет делать предсказания для распадов каонов с высокой точностью.

3. Распад K⁺ → e⁺νμ⁺μ⁻: Анатомия процесса

Теперь перейдем к главному герою нашего курса. Мы изучаем процесс, в котором положительный каон распадается на четыре частицы:

Где:

— распадающийся каон.

— позитрон (античастица электрона).

— электронное нейтрино.

и — пара мюонов (положительный и отрицательный).

Этот процесс является подвидом более общего класса распадов , где и — это лептоны (электроны или мюоны).

!Диаграмма процесса: рождение лептонной пары через виртуальный фотон

Механизм распада

Этот распад происходит следующим образом:

Внутри каона кварки аннигилируют через слабое взаимодействие, рождая виртуальный -бозон.

-бозон распадается на позитрон и нейтрино .

В процессе этого распада испускается виртуальный фотон (электромагнитное взаимодействие).

Виртуальный фотон превращается в пару мюонов .

Таким образом, этот распад интересен тем, что он объединяет в себе сразу два фундаментальных взаимодействия: слабое (распад каона) и электромагнитное (рождение пары мюонов).

4. Форм-факторы: Заглядывая внутрь «черного ящика»

Самая сложная часть теоретического описания — это то, как именно кварковая структура каона влияет на испускание фотона. В физике мы прячем это незнание (или сложное знание) в функции, называемые форм-факторами.

Амплитуда распада (вероятность процесса) зависит от форм-факторов, которые описывают распределение заряда и токов внутри каона. Для нашего распада ключевую роль играют:

* Векторный форм-фактор () * Аксиальный форм-фактор ()

Теория ChPT предсказывает значения этих форм-факторов. Эксперимент NA62 измеряет параметры распада, чтобы проверить, совпадают ли реальные значения и с предсказаниями теоретиков.

Дифференциальная ширина распада (скорость распада в зависимости от энергии частиц) может быть записана в упрощенном виде как:

Где:

— вероятность распада в зависимости от инвариантной массы.

— квадрат инвариантной массы мюонной пары (энергия, переданная фотону).

— векторный форм-фактор.

— аксиальный форм-фактор.

— знак пропорциональности.

Если эксперимент покажет отклонение от этой формулы, это будет означать, что либо наши расчеты в ChPT неверны, либо (что гораздо интереснее) существует Новая физика.

5. Мотивация: Зачем нам это нужно?

Изучение распада в эксперименте NA62 имеет три главные цели.

А. Проверка Киральной теории возмущений (ChPT)

Как уже упоминалось, ChPT — это эффективная теория. Она работает хорошо, но имеет свои границы применимости. Этот распад позволяет протестировать предсказания ChPT в области энергий, близких к массе каона. Это «стресс-тест» для теории: выдержит ли она проверку высокой точностью эксперимента NA62?

Б. Лептонная универсальность (Lepton Flavor Universality)

В Стандартной модели лептоны (электрон, мюон, тау) считаются идентичными во всем, кроме массы. Это свойство называется лептонной универсальностью. Взаимодействие должно быть одинаковым для электрона и мюона.

Сравнивая распад с похожими распадами (например, ), мы можем проверить, действительно ли природа относится к разным лептонам одинаково. Любое нарушение этого принципа стало бы сенсацией мирового масштаба.

В. Поиск Новой физики

Редкие распады — это идеальное место для поиска явлений, выходящих за рамки Стандартной модели. Если существуют новые, еще не открытые частицы (например, темные фотоны или тяжелые нейтрино), они могут вмешиваться в этот процесс на квантовом уровне, слегка изменяя вероятность распада или угловые распределения вылетающих частиц.

6. Эксперимент NA62

Эксперимент NA62 расположен в ЦЕРН и использует протонный пучок от супер-протонного синхротрона (SPS). Его главная цель — измерение ультра-редкого распада . Однако, благодаря высокой интенсивности пучка и отличным детекторам, NA62 является «фабрикой каонов», позволяющей изучать и другие редкие процессы, такие как наш , с беспрецедентной точностью.

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

Заключение

Распад — это сложный танец кварков и лептонов, управляемый законами Квантовой Хромодинамики и электрослабого взаимодействия. Его изучение позволяет нам:

Уточнять параметры эффективных теорий (ChPT).

Проверять фундаментальные симметрии природы (лептонная универсальность).

Искать следы Новой физики.

В следующей лекции мы перейдем от теории к практике и рассмотрим устройство экспериментальной установки NA62.

Экспериментальная установка NA62: детекторы и система сбора данных

В предыдущей лекции мы обсудили теоретическую мотивацию изучения редкого распада . Мы выяснили, что этот процесс является уникальной лабораторией для проверки Киральной теории возмущений и поиска Новой физики. Однако теория мертва без экспериментального подтверждения.

Сегодня мы перенесемся в Северную зону (North Area) ЦЕРН, чтобы изучить «микроскоп», который позволяет нам видеть эти редчайшие события. Эксперимент NA62 — это сложнейший комплекс детекторов, растянувшийся почти на 300 метров. Его задача — найти «иголку в стоге сена», где «стог» — это миллиарды обычных распадов, а «иголка» — наше искомое событие.

1. Пучок частиц: Откуда берутся каоны?

Всё начинается с Супер-протонного синхротрона (SPS) — второго по величине ускорителя в ЦЕРН после Большого адронного коллайдера (LHC). SPS разгоняет протоны до огромных энергий.

!Схема формирования вторичного пучка адронов в эксперименте NA62

Процесс рождения пучка выглядит так:

Первичный пучок: Протоны с импульсом 400 ГэВ/c (гигаэлектронвольт на скорость света) выводятся из SPS.

Мишень: Протоны ударяют в бериллиевую мишень (T10). При столкновении рождается ливень вторичных частиц: пионов, протонов, каонов и других.

Селекция: С помощью системы магнитов и коллиматоров мы отбираем только положительно заряженные частицы с импульсом .

Важно понимать, что полученный пучок — смешанный. На один каон () приходится примерно 15 других частиц (в основном пионов и протонов ). Доля каонов в пучке составляет всего около 6%.

2. Система трекинга: Глаза эксперимента

Чтобы восстановить картину распада, нам нужно знать траектории частиц до и после распада. Для этого используются трековые детекторы (трекеры).

GTK: Гига-трекер (GigaTracker)

Первым в игру вступает детектор GTK. Он установлен на пути пучка до зоны распада. Это кремниевый пиксельный детектор, который работает в условиях колоссальной радиационной нагрузки.

Его задачи: * Измерить импульс входящего каона. * Определить точное время пролета. * Определить угол влета частицы.

GTK позволяет нам знать 4-импульс «родителя» (каона) с высокой точностью.

Вакуумная камера распада

После прохождения GTK частицы попадают в огромную вакуумную трубу длиной около 65 метров. Вакуум необходим, чтобы каоны не сталкивались с молекулами воздуха, что создало бы паразитный фон. Именно здесь происходит магия: каон распадается.

STRAW: Строу-трекер

Если каон распался внутри вакуумной трубы, продукты распада (в нашем случае ) летят дальше и попадают в спектрометр STRAW. Он состоит из тысяч тонких трубок (соломинок), наполненных газом.

Посередине спектрометра установлен мощный дипольный магнит (MNP33). Магнитное поле искривляет траектории заряженных частиц. По радиусу кривизны мы вычисляем импульс частицы, используя формулу Лоренца:

Где: * — импульс частицы (в ГэВ/c). * — коэффициент пересчета размерностей. * — индукция магнитного поля (в Тесла). * — радиус кривизны траектории (в метрах).

Чем меньше искривляется траектория (больше ), тем выше импульс частицы.

3. Идентификация частиц (PID): Кто есть кто?

Знать импульс недостаточно. Нам нужно отличить электрон от мюона или пиона. Для нашего распада это критически важно, так как нам нужно точно идентифицировать три заряженные частицы: один позитрон и два мюона.

RICH: Детектор черенковского излучения

RICH (Ring Imaging CHerenkov detector) — это огромный бак, заполненный газом неоном. Когда частица летит быстрее скорости света в этой среде, она испускает конус света (черенковское излучение).

!Принцип работы RICH: конусы черенковского света превращаются в кольца на детекторах

Угол раствора этого конуса зависит от скорости частицы :

Где: * — черенковский угол. * — показатель преломления среды (для неона ). * — скорость частицы в единицах скорости света.

Зная импульс (из STRAW) и скорость (из RICH), мы можем вычислить массу частицы и понять, кто перед нами: легкий электрон (большое кольцо), средний мюон или тяжелый пион (маленькое кольцо).

LKr: Электромагнитный калориметр

Далее стоит калориметр на жидком криптоне (LKr). Это «ловушка» для частиц. Он измеряет полную энергию, оставленную частицей. * Электроны/позитроны полностью застревают в LKr, отдавая всю свою энергию. * Мюоны и пионы пролетают сквозь него, оставляя лишь малую часть энергии (ионизационные потери).

Для нашего анализа LKr — главный инструмент для подтверждения того, что кандидат в позитрон действительно является позитроном (отношение энергии к импульсу должно быть близко к единице: ).

MUV: Мюонная система вето

В самом конце установки находятся мюонные детекторы (MUV1, MUV2, MUV3). Они спрятаны за толстой железной стеной. * Электроны и адроны (пионы) застревают в железе. * Только мюоны обладают высокой проникающей способностью и долетают до детекторов MUV.

Если детектор MUV сработал — мы почти наверняка видим мюон. Для канала нам нужно зарегистрировать два сигнала в MUV (для и ).

4. Триггер и сбор данных (TDAQ)

Интенсивность пучка в NA62 колоссальна: на детекторы прилетает около 750 миллионов частиц в секунду (750 МГц). Записать информацию о каждой из них невозможно — не хватит никаких жестких дисков.

Для отбора интересных событий используется многоуровневая система триггеров:

Уровень 0 (L0): Аппаратный триггер. Специальные электронные платы (FPGA) анализируют сигналы с детекторов в реальном времени. Если событие похоже на интересное (например, есть сигнал в мюонном детекторе), дается команда на считывание. Это снижает поток данных с 750 МГц до 1 МГц.

Уровень 1 (L1): Программный триггер. Компьютерная ферма проводит быстрый анализ собранных фрагментов данных. Если событие проходит проверку, оно записывается на диск.

5. Специфика для распада

Для успешного анализа нашего редкого распада установка должна работать как единый организм. Вот чек-лист успешного события:

GTK: Видит один входящий трек ().

STRAW: Видит три исходящих трека, выходящих из одной точки (вершины распада).

Заряд: Суммарный заряд трех треков должен быть равен заряду каона (+1). То есть . Это соответствует составу .

PID:

* Один трек идентифицирован как (сигнал в LKr, кольцо в RICH). * Два трека идентифицированы как мюоны (сигналы в MUV, кольца в RICH).

Кинематика: Мы не видим нейтрино , но можем вычислить его присутствие через «недостающую массу».

Квадрат недостающей массы вычисляется как:

Где: * — квадрат недостающей массы. * — 4-импульс каона (из GTK). * — 4-импульсы продуктов распада (из STRAW и калориметров).

Если в распаде действительно улетело только нейтрино (которое почти не имеет массы), то величина должна быть близка к нулю.

Заключение

Эксперимент NA62 — это триумф инженерной мысли, позволяющий регистрировать частицы с невероятной точностью во времени и пространстве. Комбинация точного трекинга (GTK, STRAW) и мощной идентификации частиц (RICH, LKr, MUV) дает нам шанс выделить несколько десятков событий распада из триллионов фоновых процессов.

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

Реконструкция событий и критерии отбора сигнала

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

Но детекторы — это лишь «железо». Они выдают нам не готовые названия частиц («вот пролетел электрон»), а набор электрических сигналов: время срабатывания, номер сработавшего пикселя, количество собранного заряда. Этот сырой поток данных (raw data) похож на набор разрозненных пазлов.

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

1. От сигналов к частицам: Процесс реконструкции

Прежде чем искать конкретный распад, нам нужно превратить показания приборов в физические объекты. Этот процесс называется реконструкцией и выполняется специальным программным обеспечением на серверах ЦЕРН.

Построение треков (Tracking)

Представьте, что частица летит сквозь детекторы, оставляя точки срабатывания (хиты) в каждой плоскости. Задача алгоритма — соединить эти точки одной линией.

Трек пучка: Используя данные детектора GTK, мы строим прямую линию, показывающую траекторию каона до распада.

Треки продуктов распада: Используя данные спектрометра STRAW, мы ищем траектории вторичных частиц. Поскольку в STRAW есть магнит, эти линии будут искривлены. По кривизне линии мы вычисляем импульс частицы.

Поиск вершины (Vertexing)

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

Эта точка называется вершиной распада (decay vertex). Для нашего сигнала вершина должна находиться строго внутри вакуумной трубы.

!Схема реконструкции вершины распада: соединение трека каона и треков вторичных частиц.

2. Топология сигнала: Что мы ищем?

В физике частиц «топология» означает геометрическую конфигурацию события. Для нашего распада топология выглядит так:

* 1 входящий трек: Положительно заряженный каон (). * 3 исходящих трека: Три заряженные частицы, вылетающие из одной точки. * Суммарный заряд: Электрический заряд должен сохраняться. Заряд каона . Значит, сумма зарядов продуктов тоже должна быть . Нам нужна комбинация , то есть две положительные частицы и одна отрицательная.

Если событие не соответствует этой базовой геометрии, мы его сразу отбрасываем. Это первый и самый грубый фильтр.

3. Идентификация частиц (PID): Фейс-контроль

После того как мы нашли три трека, нам нужно понять: кто есть кто? Нам нужны один позитрон () и пара мюонов (). Главные враги здесь — пионы (), которых в эксперименте рождается огромное количество.

Мы применяем критерии отбора (cuts) к каждому треку.

Критерии для позитрона ()

Позитрон — это легкая частица, которая охотно взаимодействует с веществом. Для его идентификации мы используем калориметр LKr.

Ключевая переменная — отношение энергии к импульсу:

Где: * — отношение энергии к импульсу. * — энергия, измеренная калориметром (в ГэВ). * — импульс, измеренный спектрометром (в ГэВ/c).

Для электронов и позитронов это отношение должно быть близко к 1.0 (почти вся энергия поглощается калориметром). Для мюонов и пионов это число обычно намного меньше 1.

Наш критерий: .

Критерии для мюонов ()

Мюоны — это «тяжелые электроны», которые обладают высокой проникающей способностью.

Калориметр: Они оставляют мало энергии в LKr ().

Мюонный детектор (MUV): Они должны пробить железную стену и оставить сигнал в детекторах MUV3.

Если трек ведет себя как «призрак» (проходит сквозь калориметр и железо), мы помечаем его как мюон.

4. Кинематический анализ и недостающая масса

Даже если мы нашли правильные частицы (), это еще не значит, что мы нашли наш распад. Существуют фоновые процессы, которые могут имитировать этот набор.

Главная особенность нашего сигнала — наличие нейтрино (). Нейтрино невидимо для детекторов NA62. Оно улетает, унося с собой энергию и импульс.

Как нам «увидеть» невидимое? Мы используем закон сохранения энергии и импульса. Мы знаем, что было в начале (каон), и видим почти всё, что получилось в конце, кроме нейтрино.

Мы вычисляем квадрат недостающей массы ():

Где: * — квадрат недостающей массы системы. * — 4-импульс начального каона. * — 4-импульсы зарегистрированных позитрона и мюонов.

Физический смысл: Если в распаде действительно улетело только нейтрино (масса которого ничтожно мала, ), то величина должна быть близка к нулю.

Если же мы ошиблись и пропустили какую-то тяжелую частицу (например, ), то будет большим.

!Распределение недостающей массы: сигнал концентрируется около нуля.

5. Борьба с фоном

Главный враг нашего анализа — распад . Это очень частый процесс (около 6% всех распадов каона). Как он может притвориться нашим сигналом?

Представьте ситуацию:

Один ошибочно принимается за позитрон (ошибка PID).

Другой и ошибочно принимаются за мюоны (распадаются на лету).

Чтобы убить этот фон, мы используем тот факт, что масса трех пионов фиксирована ( МэВ). Если мы посчитаем инвариантную массу трех треков, предполагая, что это пионы, и получим массу каона — значит, это фоновое событие . Мы вырезаем такие события из данных.

6. Слепой анализ (Blind Analysis)

В современной физике частиц принят стандарт «слепого анализа». Это защита от человеческого фактора.

Ученый очень хочет найти редкий распад. Если он будет видеть данные в реальном времени, он может подсознательно подкручивать критерии отбора так, чтобы «натянуть» результат на желаемое открытие.

Чтобы этого избежать, область сигнала (region of interest) закрывается (masking). Мы видим данные вокруг сигнала (боковые полосы), настраиваем на них наши алгоритмы, проверяем уровень фона, и только когда мы на 100% уверены в своих методах, мы «открываем ящик» (unblinding).

> «Нельзя смотреть в ответ в конце учебника, пока не решил задачу сам».

Заключение

Реконструкция и отбор событий — это многоступенчатая воронка:

Триггер: Снижает поток с 750 МГц до 1 МГц.

Реконструкция: Превращает сигналы в треки.

Топология: Оставляет только события с 3 треками.

PID: Оставляет только события с составом .

Кинематика: Оставляет только события с недостающей массой .

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

Моделирование и методы подавления фоновых процессов

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

Но как физик может быть уверен, что отобранные события — это действительно искомый сигнал, а не случайный шум? И как узнать, какую долю реальных событий мы потеряли из-за несовершенства приборов?

В этой статье мы погрузимся в мир Монте-Карло моделирования и детально разберем стратегии борьбы с фоновыми процессами, которые пытаются замаскироваться под наш редкий распад.

1. Виртуальный эксперимент: Метод Монте-Карло

В физике высоких энергий невозможно провести анализ данных без создания точной цифровой копии эксперимента. Мы называем это MC-моделированием (Monte Carlo simulation).

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

Этапы моделирования

Процесс создания «виртуальных данных» состоит из двух шагов:

Генерация событий (Event Generation): Специальные программы (генераторы) создают список частиц, вылетающих из точки распада, основываясь на теоретических формулах Стандартной модели. Например, мы «просим» программу создать 10 миллионов распадов .

Транспорт частиц (Detector Simulation): Здесь в игру вступает пакет Geant4. Это библиотека, которая знает физику прохождения частиц через материю. Она рассчитывает, как частицы взаимодействуют с кремнием трекеров, газом в трубках STRAW и жидким криптоном в калориметре. Мы моделируем отклики электроники, шумы и даже «битые» пиксели.

!Схема потока данных: сравнение пути реальных данных и симуляции Монте-Карло.

2. Понятие аксептанса (Acceptance)

Ни один детектор не идеален. У него есть «слепые зоны», ограниченное разрешение и эффективность регистрации меньше 100%. Доля событий, которые мы способны зарегистрировать и правильно отобрать, называется аксептансом (или эффективностью регистрации).

Мы вычисляем аксептанс () с помощью MC-моделирования:

Где:

— аксептанс (безразмерная величина от 0 до 1).

— количество событий, которые успешно прошли реконструкцию и все критерии отбора в симуляции.

— полное количество сгенерированных событий.

Для редких распадов аксептанс обычно составляет от 1% до 10%. Это означает, что мы видим лишь малую часть того, что происходит на самом деле. Зная , мы можем восстановить истинное количество распадов.

3. Главные враги: Фоновые процессы

Фон — это любые процессы, которые имитируют сигнатуру нашего сигнала ( + недостающая энергия), но им не являются. Поскольку мы ищем редкий распад (вероятность ), даже обычные распады (вероятность ), происходящие часто, могут создать огромный фон, если детектор ошибется.

Основной фон: ()

Этот распад происходит в 6% случаев. В конечном состоянии мы имеем три заряженных пиона. Наш сигнал требует позитрон и два мюона. Как три пиона могут превратиться в ?

Неправильная идентификация (MisID): Один детектор ошибочно принимает за позитрон . Вероятность такой ошибки мала, но не нулевая.

Распад на лету (Decay in flight): Пионы нестабильны. и могут распасться прямо внутри детектора на мюон и нейтрино: . Детектор видит трек, который начинается как пион, а заканчивается как мюон, и регистрирует его как мюон.

В итоге мы получаем набор частиц, который выглядит как сигнал, но имеет другое происхождение.

4. Методы подавления фона

Чтобы избавиться от фона, мы используем комбинацию кинематики и идентификации частиц.

Кинематическое подавление

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

Когда мы рассчитываем для фонового события, мы подставляем в формулу массу мюона () вместо массы пиона (). Эта ошибка сдвигает значение далеко от нуля.

Где:

— вычисленная недостающая масса для фонового события.

— массы каона, электрона и мюона соответственно.

Благодаря этому сдвигу, сигнал и фон разделяются на графике.

!Разделение сигнала и фона по переменной недостающей массы.

Подавление с помощью PID (Particle ID)

Даже если кинематика не помогла, у нас есть детекторы RICH и LKr. Мы требуем строгих условий:

* Вероятность ошибки : Мы настраиваем алгоритмы так, чтобы вероятность принять пион за электрон была меньше . * Вероятность ошибки : Аналогично для мюонов.

Перемножение этих малых вероятностей дает колоссальное подавление фона. Если вероятность одной ошибки , то вероятность трех одновременных ошибок будет .

5. Нормировка: Как посчитать вероятность распада?

Конечная цель эксперимента — измерить вероятность распада (Branching Ratio, ). Мы не можем просто посчитать количество сигнальных событий , потому что не знаем точно, сколько всего каонов влетело в установку (счетчики могут сбиваться при потоке 750 МГц).

Поэтому мы используем метод нормировки. Мы выбираем частый и хорошо изученный распад, например (назовем его нормировочным каналом, ), и считаем отношение.

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

Где:

— искомая вероятность распада .

— количество отобранных событий сигнала и нормировки в реальных данных.

— геометрический аксептанс для сигнала и нормировки (из MC).

— эффективность триггера (срабатывания электроники).

— известная из справочников вероятность нормировочного распада.

Этот метод позволяет сократить многие систематические ошибки, так как неточности в измерении потока каонов влияют и на числитель, и на знаменатель одинаково.

6. Валидация фона: Контрольные регионы

Мы не можем слепо верить Монте-Карло. Программа может ошибаться в моделировании «хвостов» распределений. Чтобы проверить качество моделирования, физики используют контрольные регионы (Control Regions).

Контрольный регион — это набор критериев отбора, специально инвертированных так, чтобы запретить сигнал и выделить только фон.

Пример: Мы намеренно выбираем события с большой недостающей массой (где сигнала быть не должно). Если в этом регионе количество событий в реальных данных совпадает с предсказанием MC, значит, мы правильно понимаем физику фона. Если есть расхождение — нужно вводить поправки в симуляцию.

Заключение

Моделирование и подавление фона — это битва за чистоту данных. Используя мощь Geant4, мы рассчитываем аксептанс (), а применяя хитрые кинематические трюки и PID, мы отсекаем миллиарды паразитных событий .

Использование нормировочного канала позволяет нам превратить количество найденных событий в физическую константу — вероятность распада. Теперь, когда у нас есть все инструменты, в следующей, заключительной статье курса мы обсудим полученные результаты, оценку систематических ошибок и то, подтвердил ли эксперимент NA62 предсказания Стандартной модели.

Статистический анализ, нормализация и вычисление вероятности распада

Добро пожаловать на заключительную лекцию курса, посвященного анализу редкого распада в эксперименте NA62. Мы прошли долгий путь: от теоретических основ и кварковой структуры каона до сложнейших детекторов и алгоритмов фильтрации данных.

К этому моменту у нас на руках есть «чистый» набор данных: несколько событий-кандидатов, которые прошли все наши строгие фильтры, и оценка фонового шума. Но просто сказать «мы нашли 5 событий» — это не наука. Наука начинается там, где появляются числа с погрешностями и сравнение с теорией.

Сегодня мы займемся статистическим анализом. Мы узнаем, как превратить количество найденных событий в физическую величину — вероятность распада (Branching Ratio), как оценить точность нашего измерения и как правильно объявить об открытии (или его отсутствии).

1. Проблема абсолютного счета

Представьте, что вы хотите узнать, с какой вероятностью при подбрасывании монеты выпадает «орел». Вы бросили монету много раз и насчитали 50 «орлов». Можете ли вы вычислить вероятность? Нет, если вы не знаете, сколько всего раз вы бросали монету.

В эксперименте NA62 ситуация похожая. Мы знаем количество найденных распадов сигнала (), но нам нужно знать полное количество каонов (), пролетевших через детектор за время эксперимента.

Где:

(Branching Ratio) — вероятность (доля) распада по конкретному каналу.

— количество зарегистрированных сигнальных событий.

— полное количество распавшихся каонов.

Проблема в том, что посчитать напрямую с высокой точностью очень сложно. Интенсивность пучка колеблется, детекторы имеют «мертвое время», эффективность триггера меняется. Прямой подсчет дал бы огромную ошибку.

2. Метод нормализации

Чтобы обойти эту проблему, физики используют метод относительной нормализации. Мы не считаем полное число каонов. Вместо этого мы сравниваем наш редкий сигнал с другим, очень частым и хорошо изученным распадом того же каона.

В качестве такого «эталона» (нормалировочного канала) мы используем распад (обозначим его индексом ). Его вероятность известна с высокой точностью из прошлых экспериментов.

Мы составляем уравнение отношения:

Где:

— количество найденных событий сигнала и нормализации в наших данных.

— вероятности распадов (искомая и известная).

— геометрический аксептанс (доля событий, попадающих в геометрию детектора, вычисляется из Монте-Карло).

— эффективность регистрации (триггер, реконструкция, отбор).

Из этого уравнения мы выражаем искомую величину :

В чем магия этого метода? Обратите внимание, что количество каонов исчезло из формулы. Оно сократилось, так как оно одинаково и для сигнала, и для нормализации (они летят в одном пучке). Более того, многие систематические ошибки, связанные с работой детекторов, тоже сокращаются или сильно уменьшаются, так как они влияют на оба канала схожим образом.

[VISUALIZATION: схематичное изображение весов, на одной чаше которых лежит редкий золотой самородок (сигнал), а на другой — гора обычных камней (нормализация). Подпись: Метод нормализации позволяет взвесить редкий сигнал относительно известного процесса, исключая необходимость знать общий вес горы.]

3. Учет фона и статистика

В реальности количество событий, которое мы насчитали в сигнальном регионе (), состоит из реальных сигнальных событий и случайно просочившегося фона ().

Поэтому чистое число сигнальных событий вычисляется так:

Где:

— оценочное число истинных сигнальных событий.

— число событий, наблюдаемых в эксперименте после всех отборов.

— ожидаемое число фоновых событий (оценено с помощью Монте-Карло и контрольных регионов).

Статистическая погрешность

Поскольку радиоактивный распад — процесс случайный, он подчиняется статистике Пуассона. Если мы насчитали событий, то статистическая погрешность (ошибка) этого измерения составляет корень из этого числа:

Где:

— статистическое стандартное отклонение.

— количество событий.

Например, если мы нашли 100 событий, погрешность будет . Если нашли 4 события, погрешность (то есть 50%!). Именно поэтому для редких распадов так важно набрать как можно больше статистики.

4. Систематические ошибки: Главный враг экспериментатора

Если статистическая ошибка уменьшается просто по мере набора данных (чем дольше работает эксперимент, тем точнее результат), то систематическая ошибка никуда не исчезает сама по себе.

Систематика — это наше незнание о несовершенстве приборов или моделей. Она смещает результат в одну сторону.

Основные источники систематики в нашем анализе:

Триггер: Мы можем неточно знать эффективность срабатывания электроники. Ошибка в 1% здесь даст ошибку в 1% в ответе.

PID (Идентификация частиц): Насколько точно мы знаем эффективность разделения пионов и электронов? Различия между данными и Монте-Карло вносят вклад в систематику.

Модель фона: Мы оценили фон как 1.5 события. А что, если на самом деле их 2.0? Эта неопределенность тоже является систематической ошибкой.

Суммирование ошибок

В физике принято разделять эти два типа ошибок. Финальный результат записывается в виде:

Где:

— измеренная величина.

— центральное значение.

— статистическая ошибка.

— систематическая ошибка.

Чтобы получить полную ошибку, мы складываем их «в квадратуре» (как катеты в теореме Пифагора):

Где:

— полная (тотальная) ошибка измерения.

— статистическая ошибка.

— систематическая ошибка.

5. Сравнение с теорией и интерпретация

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

Предположим, наш результат:

Где:

— экспериментально полученное значение.

Теперь мы открываем теоретические статьи по Киральной теории возмущений (ChPT), о которых говорили в первой лекции. Допустим, теория предсказывает:

Где:

— теоретическое предсказание.

Мы видим, что числа совпадают в пределах погрешностей. Это означает:

Триумф Стандартной модели: Наши представления о взаимодействии кварков при низких энергиях верны.

Отсутствие Новой физики: В пределах точности нашего эксперимента мы не увидели влияния неизвестных частиц.

Если бы результат отличался на 3-5 стандартных отклонений (сигм), это было бы указанием на открытие.

[VISUALIZATION: График с двумя точками и планками погрешностей (error bars). Одна точка красная (Эксперимент NA62), другая синяя (Теория ChPT). Точки находятся на одном уровне по вертикали, их планки погрешностей перекрываются. Подпись: Сравнение экспериментальных данных с теоретическим предсказанием. Перекрытие планок погрешностей указывает на согласие с теорией.]

6. Заключение курса

Мы завершаем наш курс по анализу редкого распада в NA62. Давайте оглянемся назад и посмотрим на полный цикл научного исследования:

Мотивация: Мы начали с вопроса «почему?». Мы поняли, что — это тест для теории ChPT.

Инструмент: Мы изучили установку NA62 — гигантский «микроскоп» длиной 300 метров.

Реконструкция: Мы разобрали, как из электрических сигналов восстановить траектории частиц.

Отбор: Мы научились отсеивать миллиарды фоновых событий, используя кинематику и PID.

Результат: Сегодня мы превратили отобранные события в физическую константу, используя метод нормализации.

Именно так работает современная физика высоких энергий. Это коллективный труд сотен людей, объединяющий квантовую теорию, передовую инженерию, программирование и статистику.

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