Образование заряженных π- и Κ- мезонов, протонов и антипротонов во взаимодействиях ρ + ρ, ρ + Pb и Pb + Pb при энергии 158 ГэВ на нуклон тема автореферата и диссертации по физике, 01.04.16 ВАК РФ

Баатар Батгэрэл АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Дубна МЕСТО ЗАЩИТЫ
2009 ГОД ЗАЩИТЫ
   
01.04.16 КОД ВАК РФ
Диссертация по физике на тему «Образование заряженных π- и Κ- мезонов, протонов и антипротонов во взаимодействиях ρ + ρ, ρ + Pb и Pb + Pb при энергии 158 ГэВ на нуклон»
 
Автореферат диссертации на тему "Образование заряженных π- и Κ- мезонов, протонов и антипротонов во взаимодействиях ρ + ρ, ρ + Pb и Pb + Pb при энергии 158 ГэВ на нуклон"

ОБЪЕДИНЕННЫМ ИНСТИТУТ ЯДЕРНЫХ ИССЛЕДОВАНИЙ

1-2009-151

На правах рукописи УДК 539.172.12+539.172.17

БААТАР Батгэрэл

ОБРАЗОВАНИЕ ЗАРЯЖЕННЫХ л- И К- МЕЗОНОВ, ПРОТОНОВ И АНТИПРОТОНОВ ВО ВЗАИМОДЕЙСТВИЯХ р +р,р + РЬ И РЬ + РЬ ПРИ ЭНЕРГИИ 158 ГэВ НА НУКЛОН

Г

Специальность: 01.04.16 — физика атомного ядра и элементарных частиц

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

12 ноя 2ВВЗ

Дубна 2009

003482920

Работа выполнена в Лаборатории физики высоких энергий Объединенного института ядерных исследований.

Научный руководитель:

доктор физико-математических наук

Георгий Левонович МЕЛКУМОВ

Официальные оппоненты:

доктор физико-математических наук,

профессор

Юрий Анатольевич ПАНЕБРАТЦЕВ

доктор физико-математических наук

Лариса Александровна ТИХОНОВА

Ведущее научно-исследовательское учреждение: Институт ядерных исследований РАН, г. Москва

Защита диссертации состоится "......."....................... 2009 года в "......." часов

на заседании диссертационного совета Д-720.001.02 в Лаборатории физики высоких энергий Объединенного института ядерных исследований, г. Дубна Московской области.

С диссертацией можно ознакомиться в библиотеке ЛФВЭ ОИЯИ. Автореферат разослан "......."....................... 2009 года

Ученый секретарь диссертационного совета кандидат физ.-мат. наук

В.А. Арефьев

Общая характеристика работы

Актуальность темы

Столкновение релятивистских тяжелых ионов представляет уникальную возможность для создания и изучения в лабораторных условиях предельно разогретой и плотной ядерной материи. В таком состоянии вещество кардинально меняет свои свойства по сравнению обычным ядерным веществом. Так, теория сильных взаимодействий квантовая хромодинамика (КХД) предсказывает, что при больших плотностях энергии может произойти переход от обычного состояния материи к новому состоянию, так называемой кварк-глюонной плазме (КГП), которая характеризуется как система, состоящая из квазисвободных кварков.

Эксперимент NA49 был спроектирован и создан для исследования рождения адронов во взаимодействиях Pb+Pb в пучках ускорителя SPS в ЦЕРН. Большой аксептанс установки для регистрации треков с высоким импульсным разрешением в больших времяпроекционных камерах ТРС и способностью к надежной идентификации частиц посредством измерения ионизационных потерь cLE/dx и времяпролетных измерений TOF позволили впоследствии существенно расширить программу исследований, дополнив ее экспериментами в нуклон-нуклонным и нуклон-ядерным столкновениях для изучения механизма образования частиц в относительно элементарных взаимодействиях р+р и р+Pb по сравнению со сложными процессами в столкновениях ядер Pb+Pb.

Настоящая работа посвящена, главным образом, определению спектров заряженных адронов 7г±, К±, р и р по поперечному импульсу (поперечной массе) частиц и анализу их характеристик (параметров). Это обусловлено тем что, что поперечные эффекты являются чистым проявлением динамики реакций, так как большая часть продольных импульсов рожденных частиц обусловлена внутренним движением сталкивающихся объектов, в то время как поперечные импульсы генерируются в самих столкновениях. Особое внимание уделено систематике параметров наклона поперечных спектров, которые являются главным инструментом в исследовании динамики взаимодействий и, в частности, изучении коллективных потоков в ядро-ядерных столкновениях. Основной интерес представляет анализ эволюции поперечных спектров частиц для различных сталкивающихся систем.

Кроме того, поперечные спектры барионов р и р и выходы (множественности) этих частиц во взаимодействиях Pb+Pb при энергии 158Л ГэВ измерены в широкой области центральности в столкновении ядер от периферических до самых центральных. Основные вопросы в анализе этих данных связаны с эффектами изменения барионной плотности (стоппинга), зависящей от степени центральности столкновения, а также аннигиляции антипротонов.

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

Цель работы

Целью работы являлось измерение и анализ спектров заряженных частиц 7г±, К±, р и р в р+р, р+Pb и центральных Pb+Pb взаимодействиях при энергии 158 ГэВ на нуклон. Одной из основных задач являлся анализ спектров заряженных протонов и антипротонов, рожденных в широкой области центральности Pb+Pb столкновений при энергии 158Л ГэВ.

Научная новизна работы

• Осуществлена калибровка времяпролетных измерений в экспериментах NA49 в пучках ускорителя SPS в ЦЕРН, позволившая провести надежную идентификацию заряженных частиц п±, К±, рирв широком интервале импульсов (2<р<12 ГэВ/с и 0<pt<2 ГэВ/с).

• Получены и проанализированы спектры частиц 7г, р и р по поперечному импульсу (поперечной массе) во взаимодействиях р+р, р+Pb и центральных столкновениях Pb+Pb, а также спектры р и р в широком интервале по центральности в столкновениях ядер Pb+Pb при энергии 158 ГэВ на нуклон. Определены основные характеристики - значения параметров наклона спектров частиц Т и выходы частиц (множественности) dn/dy в области центральной быстроты {уап = 2.9), которые были использованы для получения физических результатов в работе.

• Исследована систематика параметров наклона спектров в широкой области масс частиц, рожденных во взаимодействиях р+р, р+Pb и Pb+Pb с целью изучения коллективных радиальных потоков, возникающий при

столкновении релятивистских ядер.

• Исследована зависимость выходов частиц р и р и отношение р/р в зависимости от степени центральности в столкновении ядер Pb+Pb при энергии 158 ГэВ на нуклон для изучения стоппинга (торможения) барионов и аннигиляции антипротонов в плотной ядерной среде.

Научно-практическая значимость работы

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

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

Экспериментальные данные - значения параметров наклона и множествен-ностей тг±, К±, р и р во взаимодействиях р+р, р+Pb и Pb+Pb при энергии 158 ГэВ на нуклон, а также других величин, полученных из анализа спектров могут найти применение при планировании измерений на ускорителях частиц и ядер.

Автор защищает

1. Результаты работы по калибровке времяпролетных измерений в многоканальном детекторе в составе установки NA49 на ускорителе SPS в ЦЕРН и идентификации заряженных адронов используя данные о времени пролета и ионизационных потерях частиц, измеренных, соответственно, в детекторах TOF и ТРС.

2. Результаты измерения и анализа спектров заряженных адронов тг±, К±, р и р во взаимодействиях р+р, р+Pb и Pb+Pb, а также спектров р и р, полученных в широкой области центральности в столкновениях ядер Pb+Pb при энергии пучка 158Л ГэВ:

• выходы частиц dn/dy, параметры наклона Т и средние значения поперечных масс частиц < mt >.

• зависимость параметров наклона спектров от массы частиц.

• зависимость выходов dn/dy для рири отношения р/р от числа нуклонов Npart, участвующих во взаимодействии при столкновении ядер.

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

Апробация работы и публикации

Результаты, лежащие в основе диссертации, были доложены на совещаниях коллаборации ИА49, семинарах по релятивистской ядерной физике и международных конференциях: "Кварковая материя" Нант, Франция (2002) и Будапешт, Венгрия (2005), "Многочастичная динамика" Крым, Украина (2002), "Странная кварковая материя" Атлантик Бич, США (2003) и Кейптаун, ЮАР (2004), "Ядро-ядерные столкновения" Москва, Россия (2003).

Основные результаты диссертации опубликованы в работах [1-7].

Объем и структура работы

Диссертационная работа изложена на 90 страницах, состоит из введения, четырех глав и заключения, содержит 23 рисунка, 6 таблиц и список цитируемой литературы из 92 наименований.

Содержание диссертации

Во введении обоснована актуальность диссертации, а также сформулирована цель работы и приведено краткое содержание диссертации по главам.

В главе 1 приводится краткое описание экспериментальной установки КА49 и ее основных элементов, обеспечивших проведение экспериментов в условиях большого аксептанса, измерения треков частиц с высоким импульсным разрешением и надежной идентификацией частиц.

Схема установки КА49 представлена на рис. 1. Система мониторирова-ния пучка и триггера включает в себя сцинтилляционные счетчики 81, У0, Б2, 84 и газовые черенковские детекторы 82', 83. Счетчик 81 находится на расстоянии 27 м перед мишенью и представляет собой кварцевую пластину толщиной 200 мкм. Сигнал от которого используется в схеме триггера и служит также стартовым сигналом для времяпролетной системы. При работе с протонным пучком сигнал от черенковского излучения мал, поэтому кварцевая пластина заменяется на сцинтиллятор толщиной 5 мм. Для отбраковки "галло" пучка используется счетчик У0, который представляет собой сцинтиллятор толщиной 5 мм с круглым отверстием диаметром 1 см, который включен на антисовпадение в схеме формирования триггера. Счетчик 82 предназначен для работы в пучках протонов, а на пучках свинца используется система из двух газовых черенковских счетчиков 82' и 83, светоотражающие зеркала которых

VERTEX MAGNETS VTX-1 VTX-2

BPD-1 /

/ \ / \

/ Чч

4

VTPC-1

ГО

PESTOF-R \ ^ -'BPD 2 BPD-3 ^ -v. PESTOF-L

I a) -4—t—S-41--H-I1- 1

\ u VO 11 S? S3/

MTPC-L

-EHE3-

MTPC-R

RCAL COLL VCAL

-- BPD-2 BPD-3 ^ "BPD-2 BPD-3 CD

/ \ /

ib) — 11-Ы1—5-11--СП----o- I le)- "

\ u S2 VO u LH, S4 / \

J S2 VO

Рис 1. Схема установки NA49. Показаны детекторы пучка, сверхпроводящие магниты с времяпроекционными камерами в магнитах VTPC-1 и VTPC-2 и вне магнита MTPC-L и MTPC-R, времяпролетные детекторы TOF и калориметры RCAL и VCAL. Указаны состав пучковых камер, счетчиков и мишеней, используемых в экспериментах в реакциях Pb+Pb(a), p+p(b) и р+РЬ(с).

помещены в общий объем, заполненный гелием. Между ними установлена мишень - свинцовая фольга толщиной 224 мг/см2 0.5%A,ni).

Счетчик S2' служит для мониторирования пучка, а счетчик S3 для эффективного отбора в триггере событий взаимодействия ядер пучка с мишенью. Счетчик S4 расположен в промежутке между двумя время-проекционными камерами (VTPC-1 и VTPC-2) и используется в пучках протонов. Отсутствие сигнала с этого счетчика означает факт взаимодействия налетающего протона в мишени. При протон-протонных столкновениях используется жидкая водородная мишень размером 200 мм в длину и диаметром 30 мм.

Координаты траектории пучка определяются с помощью трех пропорциональных камер (Beam Position Detector - BPD 1,2,3), которые расположены на расстоянии 38, 16 и 6 метров от центра системы координат детектора NA49, соответственно. Считывание информации осуществляется с катодных стрипов. Камеры имеют размер 3x3 см2 и имеют координатную точность 40 мкм.

Трековая система установки состоит из четырех время-проекционных камер (Time Projection Chamber - ТРС) большого объема. Общее число электронных каналов во всех четырех ТРС равно 182000. Точность восстановления импульса частиц в камерах составляет Др/р2 = 3 х 10~5 GeV-1, а разрешение по измерению ионизационных потерь частиц в газе камер dE/dx o¡¿ ~ 4%.

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

Лшд-калориметр (RCAL) установки предназначен для измерения распределения поперечной энергии в ядро-ядерных взаимодействиях. Он расположен на расстоянии 18 м от мишени и имеет форму бочки с отверстием в центре.

На расстоянии 25 м от мишени расположен Vcio-калориметр (VCAL). Он является ключевым элементом при определении центральности взаимодействия. Калориметр состоит из двух последовательных секций. Первая (электромагнитная) состоит из слоев свинец/сцинтиллятор и имеет длину эквивалентную 16 радиционным длинам, вторая часть (адронная) длиной 7.5 ядерных длин состоит из слоев железо/сцинтиллятор. Поперечные размеры калориметра - 1x1 м2, сбор света осуществляется с помощью бокового световода и ФЭУ. Для ограничения аксептанса фрагментирующих вперед частиц перед калориметром установлен колиматор. Он представляет собой железный блок с толщиной по пучку 1 м и отверстием в центре размером 10x10 см2. Угловые размеры (аксептанс) колиматора составляют 0.3°. Вследствие этого, только непровзаимодействующие фрагменты (спекгаторы пучка) попадают в калориметр.

В этой же главе со сравнительно большой подробностью дано описание времяпролетной системы, с использованием которой в данной работе были получены результаты физического анализа. Это связано с тем, что создание 900-канального времяпролетного детектора явилось существенным вкладом Объединенного института ядерных исследований в эксперимент NA49.

Времяпролетная система включает детекторы (TOF-TL,TR) для расширения возможности идентификации частиц при относительно низких энергиях по сравнению с диапазоном идентификации в ТРС. Каждый детектор состоит из 891 счетчика. Сцинтилляционный счетчик времяпролетного детектора TOF представляет собой пластический сцинтиллятор прямоугольной формы, обернутый для светоизоляции алюминизированным майларом. В качестве фотодетектора используется фотоумножитель (ФЭУ). Сцинтиллятор приклеен к фотокатоду ФЭУ с помощью оптического клея Bicron ВС-600. К каждому ФЭУ присоединен делитель напряжения.

Конструктивным элементом сборки детектора являются алюминевые кассеты, в каждой из которых плотно укладываются 11 счетчиков. Конструкция кассеты обеспечивает необходимую жесткость и светоизоляцию. Все 81 кассеты устанавливаются компактно на трех полках механической передвижной ферме по 27 кассет на каждой, образуя из сцинтилляционных

пластиков сегмент сферической поверхности, расположенный непосредственно за большими камерами МТРС.

Основной параметр времяпролетного детектора - усредненное временное разрешение составляет а « 75 пикосекунд. Это позволяет надежно идентифицировать заряженные частицы п±, К±, р и р с импульсами до 14 ГэВ/с. Система высоковольтного питания TOF реализована в стандарте CAEN, а электроника регистрации, состоящая из модулей дискриминаторов и преобразователей TDC и QDC, выполнена в стандартах VME и FASTBUS.

Принципиальное значение для получения высокого разрешения по времени пролета частиц имеют вопросы калибровки детектора, изложенные в главе 2 вместе с методикой идентификации заряженных частиц, осуществляемых в эксперименте на основе измерений времени пролета в TOF детекторе и измерению ионизационных потерь частиц во времяпроекционных камерах ТРС. Основными этапами калибровки TOF являются: уточнение положения счетчиков, амплитудная нормировка сигналов со счетчиков, учет времени распространения света в сцинтилляторе, а также на зависимости момента срабатывания от амплитуды сигнала в счетчике.

При геометрической калибровке треки, восстановленные в ТРС, экстраполируются на TOF. Затем осуществляются сдвиги и вращения счетчиков в пространстве, добиваясь максимального количества попаданий в счетчик треков при наличии в нем сигнала. Точность определения местоположения счетчика в пространстве составляет 1 мм, что приводит к эквивалентной временной погрешности 10 пс. Так как амплитудные распределения сигналов со счетчиков TOF в каналах QDC имеют заметный разброс, то переход к Qnorm амплитуде, нормированной на значение сигнала в максимуме распределения, позволяет унифицировать критерий отбора треков для улучшения качества временной информации, используя условия 0.8 < Qnorm < 1-6.

В данной главе подробно описывается последовательность процедур коррекции временной информации, тип и степень влияния каждого фактора на временное разрешение. Для временной калибровки используются треки пионов, которые надежно идентифицируются по распределению ионизационных потерь частиц в газе ТРС. Величина (Т — Т,г), измеренного в TOF времени Т и вычисленного времени используется при определении коррекций. Время пролета Твычисляется по значению измеренного в ТРС импульса и длины трека. Коррекции проводятся на основе зависимости величины (Т - Т~) от координаты точки попадания трека в сцинтиллятор и амплитуды сигнала в данном счетчике, а также от амплитуды сигнала в стартовом счетчике.

В качестве иллюстрации на рис. 2 представлены зависимость (Т — Тж) от координаты точки попадания трека и от амплитуды сигнала в счетчике до и после применяемых коррекций.

Далее в главе изложена методика идентификации заряженных частиц и

Рис 2. Зависимость времени (Т — Т„) от координаты X и от амплитуды сигнала СЭпогт до (а, с) и после (6, сГ) коррекций.

процедура восстановления инклюзивных спектров частиц.

В результате методических исследований, были выработаны следукшдие критерии отбора треков частиц:

• трек должен быть реконструирован в камере за магнитом МТРС и, по крайней мере, в одной из камер в магните УТРС.

• если координата последней точки на треке находится внутри объема МТРС на расстоянии более 20 см от границы камеры, то частица считается распавшейся в данной точке и трек исключается из анализа.

• при попадании двух и более частиц в один счетчик треки отбрасываются из за невозможности прог>ссти их идентификации и^ орсмснп ирил^ха.

• точность определения положения отдельного счетчика в пространстве составляет порядка 1 мм, а потому в анализе используются треки, точка попадания которых в сцинтилляторе находится на расстоянии более 1 мм от края.

• так как треки со значением нормированной амплитуды Цпагт. < 0.8 и (Зпогт. > 1-6 ухудшают качество времяпролетной информации, отбираются треки, для которых 0.8 < <2погт. < 1.6.

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

Отобранные треки, попавшие в ТОР детектор, проходят затем процедуру идентификации, которая основывается на анализе двумерной зависимости

Рис 3. Характерные распределения по квадрату эффективных масс (а), ионизационным потерям частиц (<3), а также двумерные зависимости этих величин для пионов, каонов и протонов (Ь и с) в разных интервалах импульсов, используемые при идентификации частиц в эксперименте.

ионизационных потерь частиц dE/dx в газе ТРС и квадрата эффективной массы частиц m2, полученной из измерения времени пролета в TOF. Сорт частиц определяется путем отбора треков в заданных границах по т2 и dE/dx. В результате этого исключается доля частиц анализируемого сорта, а в отобранных событиях остается некоторая примесь фоновых частиц. Соответствующие коррекции на потери частиц и примесь фона производятся на основе экспериментальных распределений dE/dx и m2, используя для этого функциональное описание распределений (рис. 3). Приводится детальное описание процедуры определения параметров двумерных распределений для всех типов частиц. Суммарная эффективность в результате применения ограничений при идентификации и по критериям отбора треков составляет в среднем около 70%.

В завершении приводится описание процедуры построения двумерных распределений по быстроте у и поперечному импульсу pt (поперечной массе mt) идентифицированных частиц ж±, К±, р и р (рис. 4), на основе которых проводится анализ спектров частиц, обсуждаемых в диссертации.

В этой главе рассматриваются также вопросы моделирования установки и определения геометрической эффективности регистрации частиц (аксептанса) времяпролетного детектора с помощью пакета GEANT, которая также задается в виде двумерной зависимости pt(mt) от у. Геометрическая эффективность

Рис 4. Двумерные распределения по быстроте у и поперечному импульсу pt (поперечной масс 77?,£ — то) для 7г, К и р, идентифицированных в TOF детекторе.

изменяется от 90% для pt < 0.3 ГэВ/с до 10% при больших значениях поперечного импульса. Эффективность реконструкции треков в аксептансе TOF детектора близка к 100%. Результаты расчетов показали, что количество протонов от распадов составляет ~ 20% от полного числа частиц в спектрах протонов, а вклад антипротонов в соответствующих спектрах антипротонов равен « 35%. Примесь от вторичных пионов существенно меньше и составляет « 4%. Вклад вторичных каонов пренебрежимо мал (<1%).

Некоторые вопросы методики обработки и анализа данных, вкратце изложены в работах [5,7]

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

На рис. 5 представлены распределения по поперечной массе ггц — т. (mt = yjp\ + m2, m-масса частицы) для идентифицированных во вре-мяпролетном детекторе TOF частиц 7Г+, 7г_, К+, а также р и р, во взаимодействиях протонов р+р при энергии пучка 158 ГэВ на статистике около 3 млн. зарегистрированных событий. Известно, что при этой энергии пучка, соответствующей полной энергии в системе центра масс ,Jsññ = 17.3 ГэВ, рождается в среднем около 8 адронов.

Указанные спектры, также как и все поперечные спектры в данной работе, строились в интервале центральных быстрот, области наиболее

Рис 5. Распределения по поперечной массе частиц 7г+, 7г~, К+, К~, рирв р+рстолкновениях при энергии 158 ГэВ. Сплошными линиями показано фитирование экспериментальных точек по формуле 1.

интересной для исследования динамики взаимодействии при высоких энергиях. Обозначенные на спектрах ошибки статистические. Экспериментальные данные (рис. 5) фитировались функцией термального распределения Больцмана:

т( ¿т,1<1у \ Т

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

Значения параметров ¿п/йу, Т и ошибок этих параметров также можно получить фитируя спектры видоизмененной функцией:

1 (РИ ¿п/йу ( гщ-т\

- ехр--—— (2)

т1<1т1<1у Т{Т + т) \ Т

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

Из рис. 5 видно, что спектры заряженных пионов для выбранного интервала быстрот анализируются от значения т( — тп равного 0.3 ГэВ/с2, что определяется геометрическим аксептансом для регистрации частицы в данном интервале быстрот, а также тем обстоятельством, что при малых поперечных

Таблица 1. Значения параметров наклона спектров частиц по поперечной массе Т, плотности частиц на единицу быстроты dn/dy, средние значения поперечного импульса < pt > и поперечной массы <mt> —m в р+р взаимодействиях при энергии 158 ГэВ.

Частица У1 ■ ~У2 Т, МэВ dn/dy <pt>, МэВ <mt>-m, МэВ

7Г+ 2.6 - 3.5 176 ± 2 0.780 ± 0.030 440 ± 5 274 ± 3

7Г~ 2.6 - 3.5 174 ± 2 0.680 ± 0.040 418 ± 5 268 ±3

к+ 2.7 -3.2 177 ± 2 0.082 ± 0.001 601 ±4 229 ± 1

К- 2.7 -3.2 163 ± 2 0.058 ± 0.001 564 ±3 208 ± 1

р 2.4 -2.9 172 ± 2 0.094 ± 0.002 688 ± 8 196 ± 2

р 2.4 - 2.9 152 ±4 0.021 ± 0.001 655 ± 17 180 ± 5

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

В таблице 1 для каждой из частиц указаны интервалы быстрот у\ — у2 при построении спектров и приведены значения параметров наклона Т, выхода частиц «¿п/с?у, средних значений поперечных импульсов <рр- и поперечных масс <тг>-т.

В таблице 1 приведены также значения средних поперечных импульсов частиц <р(>. Вначале строились ргспектры, которые затем фитировались

выражением:

d?N ( mt - m .

С ■ Pt exp —-— , (3)

dptdy V T

Значения <pt> определялись интегрированием функций фитирования:

[PtTrdpt

■> apt -ть

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

Исследование протон-ядерных столкновений р+А позволяет изучить влияние ядерной среды на характеристики измеряемых частиц. В анализе использовались данные, относящиеся к так называемым "minimum bias" событиям р+РЬ взаимодействий при энергии пучка протонов 158 ГэВ. Это измерения, в которых не производился отбор по центральности столкновений. В таких взаимодействиях количество столкновений налетающего протона с нуклонами мишени (или число таких нуклонов мишени) определяется, согласно Глауберовскому представлению, величиной прицельного параметра

ю

s

a" '

I

e"

10 '' 10 г

0 0.25 0.5 0.75 1 О 0.25 0.5 0.75 1

■■1,-rn (CeV/c ) m,-m (GeV/c )

Рис 6. Распределения по поперечной массе частиц тт+, тг~, К+, Крирв р+Pb столкновениях при энергии 158 ГэВ. Сплошными линиями показано фитирование экспериментальных точек по формуле 1.

в столкновении. Оценка среднего числа столкновений в указанных событиях составляет величину V « 3.7.

На рис. 6 представлены распределения по поперечной массе, точнее по величине mt — m (m - масса частицы), для всех идентифицированных во времяпролетном детекторе TOF частиц ж+, 7Г~, К+, К~, р и р, рожденных в столкновениях р+Pb при энергии пучка 158 ГэВ. Соответствующие значения параметров наклона спектров частиц по поперечной массе Т, плотности частиц на единицу быстроты dn/dy, средние значения поперечного импульса < pt > и поперечной массы < mt > —m в столкновениях р+Pb при энергии пучка 158 ГэВ сведены в таблице 2.

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

Таблица 2. Значения параметров наклона спектров частиц по поперечной массе Т, плотности частиц на единицу быстроты dn/dy, средние значения поперечного импульса < pt > и поперечной массы < mt > —m в р+Pb столкновениях при энергии 158 ГэВ.

Частица У1 - У2 Т, МэВ dn/dy <pt>, МэВ <mt>-m, МэВ

7Г+ 2.6- 3.5 205 ± 3 1.070 ± 0.051 413 ± 6 325 ± 5

7Г~ 2.6- 3.5 191 ±3 0.915 ± 0.048 382 ± 6 300 ±4

К+ 2.7- 3.2 205 ± 3 0.151 ± 0.004 534 ± 8 298 ± 5

к~ 2.7 - 3.2 193 ± 3 0.093 ± 0.003 506 ± 9 274 ± 5

р 2.4- 2.9 203 ± 2 0.264 ± 0.005 638 ± 8 256 ±3

р 2.4- 2.9 165 ± 5 0.038 ± 0.002 560 ± 18 209 ± 6

IV ; чг

р (X 0.1)^

;

О 0.25 0.5 0.75 1 О 0.25 0.5 0.75 1

m,-m (GeV/c ) m,-ni (GeV/c1)

Рис 7. Распределения по поперечной массе частиц 7г+, К+, К~, р и р в Pb+Pb столкновениях при энергии 158Л ГэВ. Сплошными линиями показано фитирование экспериментальных точек по формуле 1.

сложным ядро-ядерным столкновениям в данной главе приводится результаты, полученные в центральных столкновениях Pb+Pb при энергии 158Л ГэВ [1, 5, 7], аналогично тем, которые были представлены выше для реакций р+р и р+Pb. Степень центральности в процессе измерения задавалась отбором 5%-ой доли наиболее центральных событий от полного сечения неупругих взаимодействий Pb+Pb при энергии 158Л ГэВ (ainei ~ 7.2 барн). Полученные результаты, распределения по поперечной массе частиц и их параметры, показаны, соответственно, на рис. 7 и в таблице 3.

Плотность частиц или выход частиц dn/dy можно также определять двумя способами, как указано выше, а именно, интегрированием фитирующей функцией и суммированием экспериментальных точек. Полученные результаты совпадают с точностью 1.5%.

Таблица 3. Значения параметров наклона спектров частиц по поперечной массе Т, плотности частиц на единицу быстроты dn/dy, средние значения поперечного импульса < pt > и поперечной массы <mt> —m в Pb+Pb столкновениях при энергии 158/1 ГэВ.

Частица У1 ~У2 Т, МэВ dn/dy <pt>, МэВ <m(>-m, МэВ

7Г+ 2.7 - 3.2 206 ± 1 130 ±2 440 ± 2 329 ± 2

7Г_ 2.7 -3.2 193 ± 1 139 ±2 407 ± 3 305 ±2

к+ 2.7 -3.2 233 ± 3 28 ± 1 610 ± 8 307 ±4

к~ 2.7 -3.2 229 ± 2 16 ± 3 580 ± 5 301 ± 3

р 2.4 -2.9 308 ±9 30 ± 1 884 ± 15 413 ± 13

р 2.4 -2.9 291 ± 15 1.7 ± 0.2 873 ± 27 384 ± 19

Приведенная оценка точности может рассматриваться в качестве одной из составляющих общей систематической ошибки в измерениях. Полная систематическая ошибка оценивается на уровне (6-8)%. Помимо указанной выше, она включает в себя также:

• ошибки в определении коэффициентов коррекции на критерии качества трека (2%)

• ошибки в определении параметров йЕ/йх и тг распределений и коэффициентов коррекции при идентификации частиц (1.5%)

• ошибка при определения геометрического аксептанса (1-3%)

• ошибка в определении вклада частиц от распадов (1-2.5%)

• ошибки в определении параметров спектров частиц при фитировании и абсолютной нормировки (2%).

Рис 8. Зависимость параметра наклона Т от массы покоя частиц (слева) в реакциях р+р, р+РЬ и центральных столкновениях ядер РЬ+РЬ, а также (справа) в реакции р+р и центральных столкновениях ядер С+С, и РЬ+РЬ при энергии 158/1 ГэВ.

Полученные результаты в форме зависимости параметров наклона Т от массы частиц тх+, К+ и р в реакциях р+р, р+РЬ и в центральных столкновениях РЬ+РЬ при энергии 158 ГэВ на нуклон [2,4] представлены на рис. 8. В р+р и р+РЬ взаимодействиях значения параметров наклона остаются практически постоянными, в то время как в центральных столкновениях РЬ+РЬ увеличение этого параметра с ростом массы частиц очевидно. Следует отметить также, что сама величина наклона для р+РЬ несколько выше, чем в р+р взаимодействиях. Это, возможно, связано с так называемым эффектом Кронина, который заключается в повышенном выходе частиц с большими поперечными импульсами в р+А столкновениях по сравнению с р+р взаимодействиями.

В эксперименте NA49 были проведены также измерения и выполнен анализ данных по взаимодействиям легких ядер С+С и Si+Si, в пучках, полученных в качестве фрагментов первичных ядер РЬ, выведенных из ускорителя SPS.

Систематика параметров наклона с включением данных в центральных столкновениях С+С и Si+Si представлена на рис. 8 [2,4]. Видно, что рост параметра наклона с увеличением массы частиц, отсутствующий в реакции р+р, становится заметным уже в С+С столкновениях и это может быть указанием на то, что формирование коллективных эффектов начинается уже в центральных столкновениях легких ядер и продолжает существенно возрастать с увеличением размеров сталкивающихся систем. Это наглядно показано на рис. 9, в котором представлена компиляция данных NA49 по систематике параметров наклона поперечных спектров в широкой области исследуемых частиц, начиная от пионов до дейтронов, включая странные частицы, образованные в центральных столкновениях Pb+Pb при энергии пучка 158 ГэВ на нуклон [5].

О 0.5 1 1Л 2

М (GeV/c2)

Рис 9. Зависимость параметра някппня т пг массы частиц (*) (с) с центральных

столкновениях ядер РЬ+РЬ при энергии 158/4 ГэВ.

Результаты показали, что параметр наклона распределений по поперечной массе постоянно растет с ростом массы рассматриваемых частиц и увеличением числа взаимодействующих нуклонов в сталкивающихся ядрах. Такое поведение объяснялось релятивистской суперпозицией локального термального распределения рожденных частиц с полем скоростей, созданным радиальным коллективным потоком при взаимодействии тяжелых ионов. Отклонение от такого поведения, линейного роста параметра наклона с увеличением массы, которое наблюдается для странных барионов, и в большой степени для мультистранных частиц Е и О, обсуждается в литературе и объясняется существенно меньшим сечением их взаимодействия с нуклонами

внутри ядра.

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

Измеряемый параметр наклона, заключает в себе две составляющие: актуальную температуру на стадии Ггееге-оШ'а Т/0 и компоненту поперечного потока. В таком представлении параметр наклона Т может быть выражен как:

где < /Зт > - средняя скорость поперечного потока, m - масса частицы. Однако, использование такого представления, в частности, предположение о фиксированной скорости потока является сильным упрощением. В настоящее время при расчетах по гидродинамической модели, многие исследования осуществляются с использованием так называемой "blast wave" параметризации (ударная волна, взрывная волна).

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

где р = atanh~lßT - поперечная быстрота, а ß-j- = ßs{r / Rmax)n зависит от выбранной формы (профиля) потока и скорости потока ßs на поверхности файрбола, а /о и К\ - модицифированные функции Бесселя.

Распределения по поперечной массе mt, полученные в эксперименте NA49 по взаимодействию ядер Pb+Pb для тг~, К~, К+, р, р, d, <р, Л и Л, а также для Б и Q при энергии налетающего ядра 158Л ГэВ [1,3], совместно фитировались по формуле (6). Результаты для тт~, КК+, p,piid получены с использованием данных по времени пролета в TOF, а для Л, Л, Е и Q по измерениям в ТРС, идентифицируя их по топологии распада и эффективной массе.

Полученные в результате анализа усредненные значения параметров составили Т/0 = 120 - 130 МэВ и ßr = 0.45 - 0.50 [1], что находится в разумном согласии с данными, полученными из совместного анализа корреляций тождественных пионов и спектров частиц. Аналогичные параметры были получены при энергии ускорителя RHIC и составили Tf0 110 МэВ, несколько меньшее значение, чем на SPS, но большее значение для скорости поперечных потоков ßx = 0.55 — 0.60, что может объясняться большим давлением в системе.

Т = Т/о + то < ßr >2,

(5)

rriTdmrdy

dN

(6)

В четвертой главе представлены результаты анализа данных по рождению протонов и антипротонов в широкой области центральности в столкновениях ядер Pb+Pb при максимальной энергии налетающих ядер ускорителя SPS равной 158Л ГэВ [5,7]. В обработке данных были использованы 320 тысяч событий с отбором центральных столкновений и 735 тысяч событий без ограничения по центральности (так называемые "minimum bias" события).

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

При наборе "minimum bias" событий Pb+Pb использовался черенковский счетчик , в котором газ заполнял область непосредственно за мишенью. Счетчик вырабатывал сигнал, с помощью которого осуществлялся запрет на запуск установки в большинстве случаев отсутствия взаимодействия ядер пучка с мишенью.

В процессе обработки сортировка событий по центральности осуществлялась разбиением их на 6 интервалов по распределению энергии Еуею в Уе£о-калориметре. Для каждого интервала определялось среднее значение количества нуклонов, участвующих во взаимодействии <Npart>, вычисляя экспериментально полное число барионов, уносимых заряженными частицами вне фазового пространства, занимаемого спектаторными частицами. Для этого использовались спектры заряженных частиц, измеренные в эксперименте NA49 в широком интервале телесного угла, близкого к 2тг, дополняя расчеты лишь небольшой экстраполяцией данных в область быстроты частиц в пучке. Моделирование с использованием генераторов событий VENUS и RQMD позволяло установить соотношение между <Npart> и диапазоном прицельного параметра столкновения Ь в каждом из 6 интервалов событий по центральности.

Во втором подходе оценивалось количество спектаторных частип в событии, используя непосредственно энергию, выделенную спектаторными фрагментами пучка в Veto-калориметре. Этот метод основывается на использовании генератора событий и моделировании отклика калориметра по программе GEANT. Оценки величины <Npart> с использованием измеренных спектров частиц и энергии в Veto хорошо согласуются и в дальнейшем при анализе используется их среднее значение, а величина разброса между ними служит в качестве оценки систематической неопределенности величины <Аграг«>- Оценка величины прицельного параметра b также основана на измерении энергии Eyct0, которая связана с b соотношением, полученным из моделирования. В дополнение, так называемое число раненых нуклонов <Nwmn,i> было получено по модели Глаубера для ядро-ядерных столкновений. Значение параметров для шести интервалов по центральности столкновений, используемых в анализе данных, приведены в таблице 4.

Таблица 4. Разбиение событий РЬ+РЬ столкновений на интервалы от центрального (№1) до периферийного (№6) при анализе данных. Приводятся отношения энерговыделения в Уе1о-калориметре к энергии пучка Еуно/^ьсат и соответствующая величина сечения(доля от Ст(0(), среднее число нуклонов, участвующих в реакции <А,раГ1> и "раненых" нуклонов <Л,11,ОТ1П^>, а также предельные значения прицельного параметра Ь.

№ Eyeto/Ebeam Доля от ош(%) <Npart> <NW0Und> Интервал b(fm)

1 0-0.25 0-5 366 357 ± 22 2.3 ± 0.9

2 0.25-0.40 5-12 309 288 ± 28 4.5 ± 0.8

3 0.40-0.58 12-23 242 211 ± 30 6.4 ± 0.8

4 0.58-0.71 23-33 178 146 ± 25 8.1 ± 0.8

5 0.71-0.80 33-43 132 99 ± 22 9.5 ± 0.8

6 0.80-1.00 43-100 85 42 ± 16 11.8 ± 1.0

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

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

Затем треки проходят процедуру идентификации, используя отбор по величине dE/dx и т2 (так называемая процедура "обрезания" по dE/dx и то2), чтобы отобрать антипротоны и протоны и исключить пионы и большую часть каонов. Эти "обрезания" выбираются с тем, чтобы оптимизировать долю выделенных частиц и минимизировать присутствие фона от других заряженных частиц. Величины коррекций, учитывающих долю исключенных из анализа антипротонов и протонов, и примесь фона оценивалась из параметризованного представления экспериментальных распределений dE/dx и т2.

Среди идентифицированных протонов присутствует примесь протонов от слабых распадов странных барионов (так называемые feeddown протоны), которые при реконструкции ошибочно восстанавливаются как треки, идущие от первичного взаимодействия ядер, и должны быть вычтены из протонных спектров. Коррекция указанного фона определялась моделированием по GEANT-y распадов Л и К гиперонов в условиях экспериментальной установки и последующей реконструкции распадных протонов по всей цепочке программ реконструкции реальных событий. При генерации гиперонов брались их распределения в фазовом пространстве, полученные в экспериментах WA97/NA57 и NA49 на ускорителе SPS по рождению Л частиц, включая Л

от распадов S0-гиперонов. Небольшой вклад фона протонов от распадов £+ вычислялся моделированием по программе RQMD.

Полный вклад распадных протонов оказался равным 20-25% и почти независящим от центральности столкновений, а спектр таких протонов по поперечной массе mt хорошо описывается экспоненциальной функцией с параметром наклона Т, медленно изменяющемся от 240 МэВ до 210 МэВ, соответственно, для наиболее центральных и периферических взаимодействий Pb+Pb при 158Л ГэВ. Оценка примеси антипротонов от слабых распадов антигиперонов была выполнена с использованием данных по рождению Л, полученных на том же экспериментальном материале в эксперименте NA49 и составила величину около 35%.

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

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

Эффективность реконструкции треков оценивалась с помощью так называемого метода "embedding tracks"- вложения моделированных треков в реальное событие, обработки таких событий и определения доли восстановленных треков из числа моделированных. Для треков, летящих в пределах аксептанса TOF детектора, эффективность реконструкции оказалась близкой к 100%.

После применения всех коррекций были получены спектры по поперечной массе (mt = + m2, m-масса частицы) для антипротонов и протонов (рис. 10). Экспериментальные точки распределения фитировались выражением, представляющим сумму двух экспоненциальных функций:

d'NJ = Сг ■ е-<">/г + С2 • (7)

mtdmtay

одна из которых содержит параметр наклона Т распределения для участка mt — m > 0.2 ГэВ/с2, где спектр хорошо описывается одиночной экспонентой, а другая функция учитывает дефицитную часть в спектре при малых mt, где наблюдается отклонение от экспоненты, вклад которой составляет около 10% в центральном столкновении и уменьшается, практически, до нуля в периферических столкновениях Pb+Pb.

Как показал анализ результатов, формула (7) хорошо воспроизводит данные, что особенно важно для центральных взаимодействий, где отклонение в экспериментальных распределениях от обычного термального распределения

Рис 10. Распределение по поперечной массе антипротонов и протонов в интервале быстрот 2.4 < у < 2.8 (ус.т,-2.9) для шести групп событий по центральности в столкновениях Pb+Pb при энергии 158Л ГэВ. Сплошными линиями обозначены результаты фитирования данных по формуле ( 7). Для наглядности каждый последующий спектр занижен в 2 раз.

Больцмановского типа выражено в большой степени.

Полученные на основе mt - спектров, выход частиц dn/dy, параметры наклона Т и средние значения поперечной массы <mt>-m, представлены в таблице 5. Выход частиц dn/dy, получен интегрированием параметризованного представления функций фитирования во всем интервале поперечной массы ггц. При определении величины <mt>—m, усредненной по всей кинематической области, осуществлялась небольшая экстраполяция спектров при больших mt, находящихся за пределами измерения. Для этого использовалась формула (7). Экстраполируемая часть спектра составила не более 10% от полного числа частиц в спектре.

Среди основных источников систематических ошибок при анализе полученных результатов рассматриваются неопределенности в величине коррекций на неэффективность регистрации при отборе одиночных заряженных треков в сцинтилляционных счетчиках TOF детектора, а также коррекций, связанных с процедурой идентификации частиц (обрезание по dE/dx и т2) и вычитанием фона. Дополнительная систематическая ошибка привносится коррекцией, учитывающей протоны (антипротоны) от распадов странных барионов (антибарионов). Кроме того, при вычислении выхода частиц dn/dy производится экстраполяция спектров за пределами области измерения, предполагая определенную форму описания mt- распределения, что также может внести вклад в систематическую ошибку измерения.

Таблица 5. Выход частиц с1п/<1у, параметр наклона Т и средняя поперечная масса <т(>-т для антипротонов и протонов (2.4 < у < 2.8) при различных центральностях в столкновениях РЬ+РЬ при 158/1 ГэВ. Ошибки приведены статистические.

Частица Центральность ¿п/йу Т <тп1>-т

(% от аш) (МэВ) (МэВ)

Р 0-5 1.66 ± 0.17 291 ± 15 384 ± 19

5-12 1.27 ± 0.11 299 ± 15 393 ± 16

12-23 1.05 ± 0.08 274 ± 22 370 ± 45

23-33 0.76 ± 0.06 269 ± 29 320 ± 35

33-43 0.55 ± 0.05 255 ± 28 309 ± 32

43-100 0.33 ± 0.04 236 ± 26 284 ± 28

Р 0-5 29.6 ± 0.9 308 ± 9 413 ± 13

5-12 22.2 ± 0.6 308 ±9 415 ± 14

12-23 14.5 ± 0.4 276 ± 9 362 ± 12

23-33 9.8 ± 0.3 273 ± 10 355 ± 12

33-43 5.7 ± 0.2 245 ± 10 315 ± 13

43-100 2.9 ± 0.1 216 ± 10 259 ± 12

Суммарная неопределенность в протонных данных складывается из приблизительно одинаковых по величине систематических и статистических ошибок и составляет около 15 МэВ для наклона Т, а для выхода частиц йп/йу варьируется от 5 до 10% в случае, соответственно, для центральных и периферических столкновений. Для выхода антипротонов оценка систематической ошибки составила около 10%, а для параметра наклона около 5% при слабой зависимости этих величин от степени центральности.

Из рассмотрения распределений по поперечной массе тг во взаимодействий РЬ+РЬ при энергии 158Л ГэВ (рис. 10) и результатов фитирова-ния (таблица 5), следует, что распределения становятся более "жесткими" по мере продвижения от периферических столкновений к центральным. Значения параметра наклона Т и средней поперечной массы <г?г(> - пг увеличиваются в центральных столкновениях, хотя следует отметить, что такое поведение несколько отличается для антипротонов и протонов, а именно, протонный наклон растет чуть быстрее, чем наклон антипротонов (рис. 11). Отклонение от экспоненциального распределения при малых значениях тп1 является другой характерной особенностью спектров, которая сильнее выражена в центральных столкновениях ядер. Такое поведение спектров в зависимости от центральности столкновений объясняется в рамках современных гидродинамических моделей коллективным расширением (радиальными потоками) как следствие давления, возникающего на начальной стадии взаимодействия релятивистских ядер и рассматривается в качестве еще

Рис 11. Зависимость средней поперечной массы <т(> - т для антипротонов и протонов от числа нуклонов Л^^, участвующих в РЬ+РЬ столкновениях при энергии 158А ГэВ.

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

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

Существенная разница между антипротонами и протонами обнаружена в характере зависимости выхода частиц (множественности) от центральности столкновений ядер РЬ+РЬ при энергии пучка 158А ГэВ [6,7]. Обе эти величины растут с увеличением центральности, но при этом выход протонов растет примерно в 2 раза быстрее. Это показано на рис. 12, где представлено отношение выхода частиц dn/dy для р и р к числу нуклонов Npart, участвующих во взаимодействии в зависимости от Npart для шести групп событий по центральности. Рост отношения (dn/dy)/Npart с увеличением центральности для р очевиден, в то время как для р это отношение остается почти постоянным.

Образование антипротонов во взаимодействиях РЬ+РЬ при максимальной энергии SPS изучалось с использованием микроскопической транспортной модели UrQMD. Оно показало, что присутствуем аннигиляции (поглощения)

Рис 12. Выход частиц dn/dy для антипротнов и протонов, нормированный на число участвующих во взаимодействии нуклонов NvaTt и отношение выходов р/р в зависимости от Npart в столкновении ядер Pb+Pb при энергии 158Л ГэВ.

антипротонов можно объяснить близкую к постоянной величину отношения (dn/dy)/Npart при изменении центральности, наблюдаемую в эксперименте. В то же время, при исключении из рассмотрения аннигиляции, а это -опция, при которой модель дает разумное описание большинства реакций в центральных столкновениях Pb+Pb при максимальной энергии SPS, расчеты показывают примерно двухкратный рост отношения при переходе от периферических к центральным столкновениям, именно то, что наблюдается в данном эксперименте для протонов, которые не подвержены аннигиляции. Существенное замечание в этих работах сводится к тому, что в адронной физике, основанной на теории струн, заложенной в модели UrQMD, имеет место повышенный выход антипротонов в центральных взаимодействиях ядер и процесс аннигиляций компенсирует повышенный выход р и приводит к почти постоянной зависимости отношения выхода р к NpaTt от центральности столкновения.

Указанные выше свойства выход протонов и антипротонов приводят к существенной зависимости отношения р/р от центральности, которая показана на рис. 12. Непрерывно возрастая при переходе от центральных столкновений к периферическим, это отношение достигает значения, примерно равного двум. Этот результат находится в хорошем количественном согласии с данными для Au+Au взаимодействии при 11Л А ГэВ на ускорителе AGS. Более того, данные подтверждают наблюдаемую тенденцию, демонстрируя еще более резкий рост отношения р/р при переходе к малым системам, и это видно также из результатов, приведенных в предыдущей

главе для реакций р+р и р+Pb. В противоположность этому, отношение р/р, измеренное на ускорителе RHIC, показывает слабую зависимость от центральности столкновения при энергии ^/swv = 130 ГэВ и, практически , отсутствие её при ^/s^n = 200 ГэВ. Таким образом, барионная плотность при взаимодействии ядер, значительная при низких энергиях, существенно влияет на выход антипротонов посредством аннигиляции барион-антибарионных пар в центральных столкновениях на AGS и SPS, но не на RHIC, где плотность барионов, образующаяся в процессе взаимодействия ядер, существенно меньше.

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

В заключении сформулированы основные выводы диссертации:

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

2. Впервые измерены инклюзивные спектры частиц 7Г±, К±, р и р в условиях одной экспериментальной установки (адронный детектор NA49 на ускорителе SPS в ЦЕРН), позволившие провести сравнительный анализ их характеристик во взаимодействиях р+р, р+Pb и центральных столкновениях ядер Pb+Pb при энергии 158 ГэВ/нуклон. Выполнены исследования образования частиц р и р в широкой области центральности взаимодействующих ядер Pb+Pb при энергии 158 ГэВ/нуклон.

3. Для широкого спектра частиц выполнен анализ зависимости параметров наклона распределений по поперечной массе mt (поперечному импульсу Pt). Показано, что значения параметров наклона возрастают с ростом массы частиц и размером сталкивающихся ядер, но остаются, практически, постоянными во взаимодействиях р+р и р+РЬ.

4. Показана сильная зависимость формы и параметров поперечных спектров р и р от степени центральности в столкновениях ядер Pb+Pb. Параметры наклона и средние значения поперечных масс

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

5. Полученные результаты (пп. 3 и 4) подтверждают предсказанные свойства инклюзивных поперечных спектров частиц, полученных в рамках современного представления о том, что в центральных ядро-ядерных взаимодействиях наряду с термальным движением частиц возникает сильное гидродинамическое расширение и коллективные радиальные потоки.

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

Основные результаты диссертации опубликованы в работах:

1. М. van Leeuwen, T.Anticic, B.Baatar et al., Recent results on spectra and yields from NA49. Nucl.Phys. A 715 (2003) 161.

2. C.Hohne, S.V.Afanasiev, T.Anticic, B.Baatar et al., System size dependence of strangeness production at 158 Л GeV. Nucl.Phys. A 715 (2003) 474.

3. V.Friese, C.Alt, T.Anticic, B.Baatar et al, Strangeness from 20A GeV to 158Л GeV. J.Phys. G: Nucl.Part.Phys. 30 (2004) 119-128.

4. i.Kraus, S.V.Afanasiev, T.Anticic, B.Baatar et al., System size dependence of strangeness production at 158Л GeV. J.Phys. G: Nucl.Part.Phys. 30 (2004) 583-588.

5. T.Anticic, B.Baatar et al., Energy and centrality dependence of deuteron and proton production in Pb+Pb collisions at relativistic energies. Phys.Rev. С 69, 024902 (2004).

6. C.Hohne, C.Alt, T.Anticic, B.Baatar et al., Results from NA49. Nucl.Phys. A 774 (2006) 35.

7. C.Alt, T.Anticic, B.Baatar et al., Energy and centrality dependence of p and p production and the A/p ratio in Pb+Pb collisions between 20Л GeV and 158Л GeV. Phys.Rev. С 73, 044910 (2006).

Получено 15 октября 2009 г.

Отпечатано методом прямого репродуцирования с оригинала, предоставленного автором.

Подписано в печать 16.10.2009. Формат 60 х 90/16. Бумага офсетная. Печать офсетная. Усл. печ. л. 1,81. Уч.-изд. л. 2,15. Тираж 100 экз. Заказ № 56741.

Издательский отдел Объединенного института ядерных исследований 141980. г. Дубна, Московская обл., ул. Жолио-Кюри, 6. E-mail: publ¡sh@jinr.ru www.jinr.ru/publish/

 
Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Баатар Батгэрэл

ВВЕДЕНИЕ

ГЛАВА 1. Эксперимент

1.1. Установка NA49 на ускорителе SPS (ЦЕРН).

1.2. Пучковые детекторы и условия триггера.

1.3. Магнитное поле.

1.4. Времяпроекционные камеры ТРС.

1.5. Калориметры.

1.6. Времяпролетный детектор TOF.

1.7. Электроника считывания и триггер TOF.

1.8. Система сбора данных TOF.

ГЛАВА 2. Обработка экспериментальных данных

2.1. Задачи реконструкции событий.

2.2. Реконструкция треков, вершины взаимодействия и времени пролета частиц.

2.3. Геометрическая калибровка времяпролетного детектора и амплитудная нормировка сигналов.

2.4. Коррекция измерений и определение времени пролета частиц.

2.5. Идентификация частиц.

2.6. Критерии качества треков.

2.7. Параметризация двумерных dE/dx — т2 распределений

2.8. Построение двумерных pt{mt) — у распределений.

2.9. Определение коррекций методом моделирования.

ГЛАВА 3. Спектры частиц тг*, /v ±, р и р в реакциях р+р, р+Pb и

Pb+Pb

3.1. Введение.

3.2. Построение инклюзивных распределений частиц.

3.3. Спектры частиц по поперечной массе.

3.4. Систематика параметров наклона спектров.

3.5. Поперечные спектры частиц и радиальные потоки.

ГЛАВА 4. Спектры частиц рирв реакции Pb+Pb

4.1. Введение.

4.2. Отбор событий по центральное ги столкновений.

4.3. Процедура обработки и коррекции данных.

4.4. Спектры частиц по поперечной массе.

4.5. Форма спектров и параметры наклона.

4.6. Выход частиц и отношение р/р.

 
Введение диссертация по физике, на тему "Образование заряженных π- и Κ- мезонов, протонов и антипротонов во взаимодействиях ρ + ρ, ρ + Pb и Pb + Pb при энергии 158 ГэВ на нуклон"

Столкновение релятивистских тяжелых ионов представляет уникальную возможность для создания и изучения в лабораторных условиях предельно разогретой и плотной ядерной материи. В таком состоянии вещество кардинально меняет свои свойства по сравнению обычным ядерным веществом. Так, теория сильных взаимодействий квантовая хромодинамика (КХД) предсказывает, что при больших плотностях энергии может произойти переход от обычного состояния материи к новому состоянию, так называемой кварк-глюоиной плазме (КГП), которая характеризуется как система, состоящая из квазисвободных кварков [1,2]. Главной целью исследований во многих ускорительных центрах было и остается поиск проявлений такого состояния в сильновзаимодействующей материи на ранней стадии ядро-ядерных взаимодействий.

Эксперимент NA49 [3] был спроектирован и создан для исследования рождения адронов во взаимодействиях Pb+Pb в пучках ускорителя SPS в ЦЕРН [4]. Большой аксептапс установки для регистрации треков с высоким импульсным разрешением в больших времяпроекционных камерах ТРС и надежной идентификации частиц посредством измерения ионизационных потерь dE/dx и времяпролетных измерений TOF позволили впоследствии существенно расширить программу исследований, дополнив ее экспериментами в нуклон-нуклонных и нуклон-ядерных столкновениях [5] с целыо изучения механизма образования частиц в сравнительно элементарных взаимодействиях и влияния холодной ядерной материи на выход частиц.

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

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

В диссертации представлены инклюзивные распределения по поперечному импульсу (поперечной массе) заряженных адронов (7г+, 7г~, К+, К~, р и р) во взаимодействиях р+р, р+РЬ и Pb+Pb при энергии 158 ГэВ на нуклон (158/1 ГэВ), измеренные в эксперименте NA49 на пучках ускорителя SPS в ЦЕРН. Спектры были использованы для сравнительного анализа данных, полученных в относительно простых взаимодействиях протонов и более сложных процессах ядерных столкновений. В частности, обсуждаются следующие основные вопросы:

• поведение (форма и параметры) спектров частиц по поперечному импульсу pt (поперечной массе mt),

• систематика параметров наклона спектров, их зависимость от типа (массы) частиц для различных сталкивающихся систем,

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

Поперечные спектры бариопов рнри выходы (множественности) этих частиц во взаимодействиях Pb+Pb при энергии 158А ГэВ измерены в широкой области центральности в столкновении ядер от периферических до самых центральных. Основные вопросы в аиализе этих данных связаны с эффектами изменения барионной плотности (стоппинга), зависящей от степени центральности столкновения, а также аннигиляции антипротонов.

Представленные в работе данные не направлены прямым образом на поиск кварк-глюонной плазмы (КГП) [6] и обнаружение ее основных признаков в столкновениях релятивистских ядер, таких как повышенный выход странности [7,8], подавление J/гр частиц [9] а также подавление частиц с большими поперечными импульсами [10,11], но их цель, скорее, в реакциях р+р, р+Pb и Pb+Pb изучить эволюцию в поведении поперечных импульсов, связанном с образованием высокой плотности энергии и барионной плотности (высокой плотности материи) па ранней стадии центральных столкновений релятивистских ионов, которые являются необходимым условием формирования нового состояния материи -кварк-глюонной плазмы.

Содержание диссертации представлено во Введении и последующих четырех главах.

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

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

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

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

Глава 3 посвящена изучению образования тт^, К±, р и р в реакциях р+р, р+Pb и центральных столкновениях ядер Pb+Pb при энергии 158Л ГэВ. Обсуждаются выходы частиц, параметры наклона и средние значения поперечных масс, полученных из спектров частиц.

Вопросы образования протонов и антипротонов в широкой области центральности в столкновениях ядер Pb+Pb при энергии 158А ГэВ рассматриваются в главе 4. Анализируются, главным образом, измеренные величины для указанных выше характеристик спектров частиц.

Большинство результатов, изложенных в диссертации, было представлено на международных конференциях: "Кварковая материя" Нант, Франция (2002) и Будапешт, Венгрия (2005), "Многочастичная динамика" Крым, Украина (2002), "Странная кварковая материя" Атлантик Бич, США (2003) и Кейптаун, ЮАР (2004), "Ядро-ядерные столкновения" Москва, Россия (2003).

Основные результаты диссертации опубликованы в работах [12-18].

 
Заключение диссертации по теме "Физика атомного ядра и элементарных частиц"

ЗАКЛЮЧЕНИЕ

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

2. Впервые измерены инклюзивные спектры частиц -к±, К±, р и р в условиях одной экспериментальной установки (адронный детектор NA49 на ускорителе SPS в ЦЕРН), позволившие провести сравнительный анализ их характеристик во взаимодействиях р+р, р+Pb и центральных столкновениях ядер Pb+Pb при энергии 158 ГэВ/нуклон. Выполнены исследования образования частиц р и р в широкой области центральности взаимодействующих ядер Pb+Pb при энергии 158 ГэВ/нуклон.

3. Для широкого спектра частиц выполнен анализ зависимости параметров наклона распределений по поперечной массе тп* (поперечному импульсу pt). Показано, что значения параметров наклона возрастают с ростом массы частиц и размером сталкивающихся ядер, но остаются, практически, постоянными во взаимодействиях р+р и р+РЬ.

4. Показана сильная зависимость формы и параметров поперечных спектров р и р от степени центральности в столкновениях ядер Pb+Pb. Параметры наклона и средние значения поперечных масс спектров р и р увеличиваются при переходе от периферических взаимодействий к центральным взаимодействиям ядер. В центральных столкновениях ядер при малых значениях поперечных масс mt в спектрах наблюдается заметное отклонение от экспоненциального Больцмановского распределения.

5. Полученные результаты (пп. 3 и 4) подтверждают предсказанные свойства инклюзивных поперечных спектров частиц, полученных в рамках современного представления о том, что в центральных ядро-ядерных взаимодействиях наряду с термальным движением частиц возникает сильное гидродинамическое расширение и коллективные радиальные потоки.

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

В заключении хочу поблагодарить профессора Малахова Александра Ивановича за постоянное внимание и поддержку на протяжении всего периода работы, основная часть которой была выполнена в Лаборатории высоких энергий ОИЯИ.

Я с большой благодарностью вспоминаю о той большой поддержке, которую оказывал Александр Михайлович Балдин участию Лаборатории в эксперименте NA49 в ЦЕРН.

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

Я благодарен моему научному руководителю Мелкумову Георгию Левоновичу за постановку задачи, постоянное внимание и большую помощь в работе.

Моя огромная признательность сотрудникам Лаборатории - коллегам, принимавшим участие в работе на разных ее этапах - Колесникову Вадиму Ивановичу за большую помощь в освоении программ и методики обработки, С.В. Афанасьеву, АЛО. Семенову, С.Г. Резникову и АЛО. Исупову за большой труд, который они внесли в создание времяпролетного детектора ОИЯИ для эксперимента NA49 в ЦЕРН, данные с которого положены в основу диссертации, а также Н.М. Агабабяну за полезные обсуждения.

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

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

 
Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Баатар Батгэрэл, Дубна

1. J.C. Collins and M.J. Perry, Phys. Rev. Lett. 34, 1353 (1975).

2. M. Kisslinger and P. Morley, Phys. Rev. D10, 2771 (1976).

3. S. Afanasiev et al, Nucl. Instrum. Methods Phys. Res. A430, 210 (1999).

4. Proposal NA35/NA49, CERN/SPSLC 91-31, SPSLC/P264 (1991).

5. H. Meyer-Ortmanns, Rev. Mod. Phys. 68, 473 (1996).

6. J. Rafelski and B. Miiller, Phys. Rev. Lett. 48, 1066 (1982).

7. J. Rafelski, Phys. Rept. 88, 331 (1982).

8. T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986).

9. M. Gyulassy and M. Pliimer, Phys. Lett. B243, 432 (1990).

10. M. Gyulassy and X.N. Wang, Nucl. Phys. B420, 583 (1994).

11. M. van Leeuwen, T. Anticic, B. Baatar et al, (NA49 Collaboration), Nucl. Phys. A715 (2003), 161c-171c.

12. C. Hohne, S.V. Afanasiev, T. Anticic, B. Baatar, et al., (NA49 Collaboration), Nucl. Phys. A715 (2003), 474c-477c.

13. V. Friese, C. Alt, T. Anticic, B. Baatar et al, (NA49 Collaboration), J. Phys. G: Nucl. Part. Phys. 30 (2004), SI 19—S128.

14. I. Kraus, S.V. Afanasiev, T. Anticic, B. Baatar et al, (NA49 Collaboration), J. Phys. G: Nucl. Part. Phys. 30 (2004), S583-S588.

15. T. Anticic, B. Baatar et al, (NA49 Collaboration), Phys. Rev. C69, 024902 (2004).

16. C. Holme, C. Alt, T. Anticic, B. Baatar et al, (NA49 Collaboration), Nucl. Phys. A774 (2006), 35-42.

17. C. Alt, T. Anticic, B. Baatar et al, (NA49 Collaboration), Phys. Rev. C73, 044910 (2006).

18. C. De Marzo et al, Nucl. Instr. and Meth. 217, 405 (1983).

19. J. Bachler et al, Phys. Lett. B184, 271 (1987).

20. K. Werner, Phys. Rep. 232, 87 (1993).

21. B.A. Григорьев и др, Энергоатомиздат, 1988

22. И.И. Анисимова и Б.М. Глуховский, Фотоэлектронные умножители, Москва, "Советское радио", 1974.

23. АЛО. Семенов, Кандидатская диссертация, ОИЯИ, Дубна, 1998.

24. С.В. Афанасьев и др., Краткие Сообщения ОИЯИ, Nl(81)-97, 1997, с.45.

25. С.В. Афанасьев и др., Краткие Сообщения ОИЯИ, N5(85)-97, 1997, с.69.

26. R. Zybert. DSPACK (version 1.0). NA49 Note 14.12.1993.

27. С. Borman et al, Nucl. Instrum. Methods Phys. Res. A374, 227 (1996).

28. Aplication Software Group, CN Division, CERN, DSPACK Data Manager, CERN, Geneva 1993.

29. В.И. Колесников, Кандидатская диссертация, ОИЯИ, Дубна, 2002.

30. R. Brun et al, CERN-DD/EE/84-1.

31. E. Fermi, Prog.Theor.Phys. 5, 570 (1950).

32. E.I. Finberg, Sov.Phys.Usp. 14, 455 (1972).

33. E.V. Shuryak, Sov.J.Phys. 20, 295 (1975).

34. R. Hagedorn, Sup. Nuovo Cimento, 3, 147 (1965).

35. G. Bearden et al, Phys. Rev. Lett. 78, 2080 (1997).

36. E. Schnederman, J. Solfrank and U. Heinz, Phys. Rev. C48, 2462, (1993).

37. N. Xu et al, Nucl. Phys. A610, 175c (1996).

38. R. Hagedorn and J. Ranft, Nucl. Phys. B48 (1972), 157.

39. F. Becaltini and G. Passavela, Eur. Phys. J. C23 (2002), 551.

40. P. Braun-Munzinger, K. Redlich and J. Stachel, e-print:nucl-th/0304013.

41. B. Nilsson-Almqvist and E.Stenlund, FRITIOF version 1.6, Comput. Phys. Commun 43 (1987) 387.

42. H. Pi, FRITIOF version 7.0 Comput. Phys. Commun. 71 (1992) 173.

43. R.J. Glauber and G. Matthiae, Nucl. Phys. B21 (1970), 135.

44. M. Kreps, NA49 Collaboration Meeting, April, 2001.

45. C. Alt, T. Anticic, B. Baatar et al, (NA49 Collaboration), Eur. Phys. J. C45, 343-381 (2006).

46. C. Alt, T. Anticic, B. Baatar etal., (NA49 Collaboration), Eur. Phys. J. C49, 897-917 (2007).

47. M.J. Tannenbaum, Int. J. Mod. Phys., A4 (1989) 3377.

48. A.Z. Mekjian, Phys. Rev. C17 (1978) 1051.

49. P.F. Kolb and U. Heinz, e-print:nucl-th/0305084.

50. J. Bachler et al, (NA49 Collaboration), Nucl. Phys. A661 (1999) 45.

51. H. Appelshauser et al, (NA49 Collaboration), Phys. Rev. Lett. 82, 2471 (1999).

52. G.E. Cooper et al, Ph. D. Thesis, University of California, Berkeley (2000).

53. I.G. Bearden et al, Phys. Rev. Lett. 78, 2080 (1997).

54. S.V. Afanasiev et al, (NA49 Collaboration), Phys. Lett. B486, 22 (2000).

55. J.W. Cronin et al., Phys. Rev. Dll, 3105 (1975).

56. C. Alt, T. Anticic, B. Baatar et al, (NA49 Collaboration), Phys. Rev. Lett. 94, 052301 (2005).

57. P. Braim-Munzinger, J. Stachel, J.P. Wessels and N. Xu. Phys. Lett. B344, 43, (1995).

58. H. Appelshauser et al, (NA49 Collaboration) Eur. Phys. J. C2, 661 (1998).

59. G. Roland, H. Appelshauser et al, (NA49 Collaboration), Nucl. Phys. A638, 91c, (1998).

60. H. van I-Iecke, H. Sorge and N. Xu, Nucl. Phys. A661, 493 (1999).

61. H. Sorge, J.L. Nagle and B.S. Kumar, Phys. Lett. B355, 27 (1995).

62. J. Soil frank et al, J. Phys. С 52, 593(1991).

63. U. Heinz, Nucl. Phys. A610, 264c (1996).

64. U. Heinz etal, Phys. Lett. B382, 181 (1996).

65. S. Chapman, J.R. Nix and U. Heinz, Phys. Rev. C52, 2694 (1995).

66. I.G. Bearden et al, Phys. Rev. Lett. 78, 2080 (1997).

67. J.Burward-Iioy and B.V.Jacak, nucl-ex/0111005 (2001).

68. Il.Sato and K.Yazaki, Phys. Lett. B98, 153 (1981).

69. E.Schnedermann and U.W. Heinz, Phys. Rev. C50, 1675 (1994).

70. G. Van Buren (STAR), Nucl. Phys. A715 (2003) 129-139.

71. J. Burward-Hoy (PHENIX), Nucl. Phys. A715 (2003) 498-501.

72. P. Koch, В. Muller, II. Stacker and W. Greiner, Mod. Phys. Lett. A3, 737 (1988).

73. J. Ellis, U. Heinz and H. Kowalski, Phys. Lett. B233, 223 (1989).

74. U. Heinz, P.R. Subramanian, II. Stacker and W. Greiner, J. Phys. G: Nucl. Part. Phys. 12, 1237(1986).

75. J. Rafelski and B. Muller, Phys. Rev. Lett. 48, 1066 (1982);

76. P. Koch, B. Muller and J. Rafelski, Phys. Rep. 142, 167 (1986).

77. M. Bleicher et al, Phys. Lett. B485, 133 (2000).

78. M. Bleicher et al, Phys. Rev. C59, 1844 (1999).

79. G.E. Cooper et al, (NA49 Collaboration), Nucl. Phys. A661, (1999) 362c.

80. F. Sikler et al., (NA49 Collaboration), Nucl. Phys. A661, (1999) 45c.

81. H. Sorge, Phys. Rev. C52, 3291 (1995).

82. R.A. Fini et al., J. Phys. G.: Nucl. Part. Phys. 27, 375 (2001); F. Antiniri et al, Nucl. Phys. A715, (2003) 140c.

83. A. Mischke et al., (NA49 Collaboration), Nucl. Phys. A715, 453c (2003).

84. T. Anticic, B. Baatar et al., (NA49 Collaboration), Phys. Rev. Lett. 93, 022302 (2004).

85. M. Gazdzicki, C. Alt., B. Baatar et al, (NA49 Collaboration), J. Phys. G: Nucl. Part. Phys. 30 (2004) 701-708.

86. Ch. Blume, S.V. Afanasiev., B. Baatar et al, (NA49 Collaboration), J. Phys. G: Nucl. Part. Phys. 31 (2005) 685.

87. I. Ahle et al, (E802), Phys. Rev. Lett. 81, 2650 (1998).

88. K. Guettler et al, Phys. Lett B64, 111 (1976).

89. I.G. Beardcn et al, Phys. Rev. C57, 837 (1998).

90. C. Adler et al, Phys. Rev. Lett. 86, 4778 (2001).

91. К. Adcox et al, Phys. Rev. C69, 024904 (2004).

92. J. Adams et al, Phys. Rev. Lett. 92, 112301 (2004).

93. S.S. Adler et al, Phys. Rev. C69, 034909 (2004).