Identified charged particle flow and unfolding event-by-event flow in heavy ion collisions

Abstract

Yksi raskasionitörmäyksen ominaisuuksista on syntyvän matalaenergisen hiukkasvirran epäisotrooppisuus. Tämä epäisotrooppisuus karakterisoidaan useimmiten Fourierin sarjakehitelmän ja sen kertoimien avulla.

Tässä työssä tutkin kahta puolta hiukkasvirrassa. Ensimmäinen on tunnistettujen hiukkasten virta ja niin sanottu kvarkkilukuskaalaus. RHIC-kiihdyttimellä on huomattu, että skaalaamalla hadronien ja mesonien havaitut isotrooppisuuskertoimet kvarkkiluvuilla erot eri hiukkasten välillä häviävät lähes kokonaan. Tämän tulkittiin tarkoittavan sitä, että hiukkasvirran kehitykseen vaikuttaa lähinnä törmäyksen partoninen vaihe ja hadronisaation jälkeisellä vaiheella ei ole vaikutusta.

Toinen puoli tutkimuksessani on hiukkasvirran muuttuminen törmäyksestä toiseen ja havaittujen Fourier-kertoimien jakauma. Kokeelliset menetelmät kuitenkin vääristävät jakaumia. Vääristymien poistamiseen käytän bayesilaista menetelmää, jota kokeilin ensin yksinkertaisella Monte Carlo -simulaatiolla.

Tämän työn olen suorittanut tutkimalla AMPT-mallilla tuotettua LHC-energiaa vastaavaa törmäysdataa keskeisyydeltään erilaisissa törmäyksissä. AMPT on hybridimalli, joka käyttää Monte Carlo -menetelmää raskasionitörmäyksen simulointiin. Havaitsin, että AMPT-malli LHC-energialla ei noudata kvarkkilukuskaalausta. Esitän myös tulokset Fourier-kertoimien jakaumista. Monte Carlo -testin perusteella uskon, että toisen asteen kertoimien jakauma vastaa todellista jakaumaa koko tutkitussa datassa. Kolmannen asteen kertoimien kohdalla bayesilainen menetelmä sen sijaan toimii ainoastaan keskeisimmissä törmäyksissä, joissa signaali on vahvin.

In this thesis I study two aspects related to anisotropic flow coefficients v2 and v3 in heavy-ion collisions. The study is done for Pb-Pb collisions at √sNN = 2.76 TeV using data simulated by A MultiPhase Transport (AMPT) model. First I study flow of identified charged particles, pions, kaons and protons. At RHIC in Au-Au collisions at √sNN = 200 GeV it has been observed that scaling v2 of identified hadrons with the number of quarks and plotting it as a function of transverse kinetic energy KET produces almost perfect scaling between different particle species. This was taken as an indication that flow is mainly generated in the partonic phase and is not strongly affected by the hadronic phase. However, in Pb-Pb col- lisions at √sNN = 2.76 TeV in LHC the scaling has been observed to break down. AMPT model uses a simple quark coalescence model, which was used to explain the scaling at RHIC energies. Because of the scaling breakdown at LHC the coalescence model has been challenged in the field. In my studies I have observed that AMPT does not produce perfect quark number scaling, even though it would be expected because of the coalescence model. Another aspect studied here is event-by-event flow. Event-by-event flow is connected to the fluctuations in the initial collisions. Only recently the field has started to study fluctuations and event-by-event flow. I will show distributions of event-by-event flow coefficients in the AMPT model. In addition to the true fluctuations the distributions have a significant smearing component from limited resolution resulting from finite multiplicity in a single event. I will use a data-driven unfolding method based on an iterative Bayesian procedure to remove the smearing effects. I will test the procedure in a toy Monte Carlo simulation to test its performance and apply it to AMPT data. I have observed that based on the Monte Carlo the procedure works for v2 in general and for v3 in central collisions.