1 of 19

SHINCHON school

2020.1.7. - 2020.1.11

Pusan National University

2 of 19

General information

  • 식사�https://docs.google.com/spreadsheets/d/1ytmI9aFHLihj-goPkXCowhTstf6qDl8fCLvROd6H-Qo/edit?usp=sharing�
  • Wifi�eduroam�PNU-Physics�
  • Indico page�https://n-ext.inha.ac.kr/event/382/�
  • Git repository�https://github.com/sh-lim/SHINCHON�
  • Slack workspace�mcshinchon.slack.com�

​

3 of 19

General information

  • Linux Server�164.125.40.30�164.125.40.31�윈도우 운영체제를 사용하는 경우 원활한 프로그램 사용이 어려울 수 있어서 위의 linux 서버를 사용하시면 됩니다.�test00-09 계정이 각 머신 마다 생성되어 있습니다. (패스워드는 계정명과 동일)�계정이 겹치지 않도록 아래 표에 이름을 적어두시고 사용하시면 됩니다. ��https://docs.google.com/spreadsheets/d/1ytmI9aFHLihj-goPkXCowhTstf6qDl8fCLvROd6H-Qo/edit?usp=sharing

ROOT 사용�source /opt/root-6.18.04-build/bin/thisroot.sh

​

​

4 of 19

Overview

  • Initial condition
    • MC Glauber�2D distribution of smeared wounded nucleons in transverse direction�Detailed format can be different for different hydro simulations �example grid size: 200x200 grid in 20x20 fm2 (0.1x0.1 fm2)�Convert into energy density by multiplying a scale factor�The scale factor should be tuned before to match charged particle multiplicity in data�N events for different centrality (geometry) can be prepared for hydro simulation
    • IP-Glasma�Initially from MC Glauber ⇒ additional fluctuation for gluons�Time evolution with weak coupling�
  • Hydro simulation
    • Read initial energy density profile
    • Time evolution until reaching a defined average temperature (150~170 MeV)
    • During the time evolution, energy density, temperature, and fluid velocity profiles for several time steps

5 of 19

Overview (Hydro simulation)

  • SONIC
    • 2+1D hydro
    • https://sites.google.com/site/revihy/home
    • Option for pre-hydrodynamic evolution based on AdS/CFT (superSONIC)
  • MUSIC
    • 2+1D or 3+1D hydro
    • http://www.physics.mcgill.ca/music/
  • For a given initial condition
    • Energy density, temperature, and fluid velocity profiles can be obtained

​

​

​

6 of 19

Heavy quark

  • Heavy quarks
    • Generate heavy quark pairs
    • Energy density profile from MC Glauber
    • pT distributions from FONLL
    • Langevin calculation

​

7 of 19

MC Glauber

  • Official page�https://tglaubermc.hepforge.org/�
  • Related papers�https://arxiv.org/abs/nucl-ex/0701025�
  • 1/7 할 일�충돌 시스템 및 impact parameter b 범위를 자유롭게 선택하여 100개의 이벤트를 만들어본다�대표 이벤트 하나의 2D 히스토그램을 저장하여 slack에 올린다.

8 of 19

IP-Glasma

  • Source file�https://github.com/schenke/ipglasma
  • Additional software�https://www.open-mpi.org/software/hwloc/v2.1/�https://www.open-mpi.org/software/ompi/v4.0/�http://www.fftw.org/download.html�./configure --prefix=XXX --enable-mpi --enable-threads�https://www.gnu.org/software/gsl/
  • src/GNUmakefile 추가로 설치해준 header, library 경로 추가�
  • openmpi installation for mac os�./configure --prefix=/xxx/openmpi-4.0.2-instal --enable-mpi-cxx�
  • Related papers�https://arxiv.org/abs/1202.6646�https://arxiv.org/abs/1807.05205�
  • 1/7 할 일�IP-Glasma compile

9 of 19

IP-Glasma

  • Configuration file�input�
  • How to run�./ipglasma input�
  • Macro to create 2D histogram�>ls -t epsilon-u* > file.lst�>root -l -b -q make_hist.C �

10 of 19

SONIC (hydro simulation)

  • Official source�https://sites.google.com/site/revihy/download
  • Required library�https://www.hdfgroup.org/downloads/hdf5/�Ubuntu (libhdf5-dev)�Macos (hdf5 from brew)
  • Compile�Mac의 경우 g++-6 대신에 g++-9도 사용 가능�xxx/uvh2-1>make initE�xxx/uvh2-1>make generate�xxx/uvh2-1>make vh2�xxx/b3d>make
  • Initial geometry for hydro�Hist2Txt.C macro in xxx/input�Require ROOT file from MC Glauber�Convert a 2D histogram of MC Glauber (Ncoll) to a text file of energy density�const int scaler = 2; //for RHIC�const int scaler = 4; //for LHC

11 of 19

SONIC (hydro simulation)

  • Copy files to template directory�>./copy_sonic.sh
  • Run hydro simulation�XXX/template>./generate�XXX/template>./initE�copy a text file from MC Glauber to input/inited.dat
  • template/data/param.txt 파일을 아래 파일을 참조하여 수정�https://github.com/sh-lim/SHINCHON/blob/master/Software/SONIC-Macos/template/data/params.txt�충분한 hadron의 수를 모으기 위해서 B3DEVENTS 100000 으로 변경
  • XXX/template>./vh2-2.1

​

12 of 19

SONIC (example output from hydro simulation)

This is VH2-2.1

IC = -2

RNUC = 6.4

ANUC = 0.54

SIGMANN = 40

TANORM = 197

NUMT = 100

B = 0

AT = 1.014

EPS = 0.001

TINIT = 0.5

TSTART = 0.370

STEPS = 2000

UPDATE = 250

SNAPUPDATE = 20000

TF = 0.170

L1COEF = 0.0

L2COEF = 0.0

FREEZE = 2

PTASIZE = 25

PHIPASIZE = 12

EOSNAME = BMWEOS

ETANAME = etaos-dyn

BETANAME = beta2-2

LAMBDANAME = lambda1-0

ZETANAME = zetaos-dyn

PTMAX = 4.0

BIEPS = 0.00

preeqflow = 0

FULL = 1

PCE = 0

MONOEPS = 0.0

SCAL = 1.0

SMOOTHING = 1

SMOOTH = 0.001

B3DEVENTS = 100

--------------------------------------------------------------------------------------------------------

==> Allocating memory

mi 50

TF=0.17

TSTART=0.37

===> Info: Successfully loaded EoS file of length 1961

===> Info: Successfully loaded eta/s file of length 451

===> Info: Successfully loaded zeta/s file of length 451

===> Info: Successfully loaded beta2 file of length 451

===> Info: Successfully loaded lambda1 file of length 451

Programmstart

Using SMOOTHING to be able to run lumpy conditions

Time[fm/c] T_cent [GeV] Ecc(spatial) Ecc(momentum)

0.5 0.24098 -0.44926 0

0.55002 0.23511 -0.24751 2.3295e-05

​

13 of 19

SONIC (hadronization)

  • Hadronization�XXX/template>./b3d default�XXX/template>./analyze default
  • Output files�analysis/default/cent0to5/details/
  • 각자 선택한 geometry에 대해서 100000 b3d 이벤트에 대한 v2(pT), v3(pT), v2/e2(pT), v3/e3(pT)를 그려보자
  • 예제 매크로 (git에 업데이트)�https://github.com/sh-lim/SHINCHON/blob/master/Software/SONIC-Macos/macro/Draw.C

​

14 of 19

SONIC (Preparation for HQ diffusion)

  • HQ diffusion을 위해서는 time step에 따라서 temperature, fluid vector의 profile이 필요!
  • data/param.txt 파일에서 SNAPUPDATE 100 으로 변경!
  • Hydro simulation(vh2-2.1)을 한 번 더 실행하면 짧
  • 은 시간 간격에서의 profile이 생성된다.�template/data/snapshot
  • git:Diffusion/macro/make_hist_profile.C 를 사용하여 텍스트 파일을 2D histogram으로 변환�(약간의 매크로 수정 필요!)

15 of 19

HQ Diffusion

  • Related paper�https://arxiv.org/abs/1307.2188

16 of 19

MUSIC (hydro simulation)

  • 기본 매뉴얼�MUSIC-XXX/notes.pdf�https://github.com/sh-lim/SHINCHON/blob/master/Software/MUSIC-Ubuntu/notes.pdf
  • Build MUSIC�MUSIC-XXX/build>rm -rf *�MUSIC-XXX/build>cmake ../�MUSIC-XXX/build>make install
  • MC Glauber input �SONIC과 MUSIC은 다른 형식의 input 파일이 필요�MUSIC-XXX/input/Hist2Txt.C 매크로를 사용�파일 이름, x-y 범위, grid 설정 등을 확인!
  • MUSIC 설정파일�MUSIC-XXX/music_input_mode_1�설정 파일 안에 “Initial_Distribution_input_filename”을 각자의 텍스트파일로 설정
  • MUSIC 실행�./mpihydro music_input_mode_1

17 of 19

Configuration in KIAF

  • 환경 설정�ALIPHYSICS_VERSION="vAN-20200216_ROOT6-1"�ALICE_ROOT_VERSION="v5-09-52b_ROOT6-1"
  • Pre-install�fftw (--enable-openmp --enable-threads --enable-shared=yes)�hdf5 (--enable-cxx --enable-hl)�hwloc ()�openmpi ()
  • MC Glauber�문제 없음
  • IPGlasma�git clone https://github.com/sh-lim/ipglasma.git�comment out lines for cxx compiler version check in CMakeLists.txt��cmake -DCMAKE_CXX_FLAGS="-std=c++11" -DGSL_INCLUDE_DIR=$GSL_DIR/include -DGSL_LIBRARY=$GSL_DIR/lib/libgsl.so -DGSL_CBLAS_LIBRARY=$GSL_DIR/lib/libgslcblas.so ../ipglasma���

���

18 of 19

Configuration in KIAF

  • SONIC�git clone -b b3d_ppbar https://@bitbucket.org/mluzum/uvh2-1.git uvh2-1�add include and lib path for GSL in Makefile��modify MADAI_GSLPATH and MADAI_HDF5_HOME in madai_repo/makefile_defs.mk�
  • MUSIC�git clone https://github.com/MUSIC-fluid/MUSIC.git�comment out lines for cxx compiler version check in CMakeLists.txt��Modify src/data.h for c++11�

���

19 of 19

https://github.com/sh-lim/SHINCHON.git