Differences

This shows you the differences between two versions of the page.

--- tutirial:pythia [2018/12/11 19:45] – minerva1993
+++ tutirial:pythia [2018/12/17 15:34] (current) – removed minerva1993
@@ Line 1: / Line 1: @@
-==== Event Generation Using Pythia ====
-Madgraph를 이용해 대부분의 이벤트를 만들 수 있지만, 가끔 복잡한 경우가 있다. (예를 들어 VV: WW, WZ, ZZ) 특수한 경우 Pythia8에서 자체적으로 decay mode를 제공하여 편리하게 이벤트를 만들 수 있다. VV의 경우 CMS simulation도 pythia를 이용해 event를 만들고 있다. ''/MG5_aMC_v2_6_4/HEPTools/pythia8/share/Pythia8/examples''에 필요한 파일들이 있긴 하지만 여기서는 별도로 Pythia를 설치하여 이용해보기로 한다.
-<code>
-cd work #some working folder
-wget http://home.thep.lu.se/~torbjorn/pythia8/pythia8235.tgz
-tar -xvf pythia8235.tgz
-cd pythia8235
-</code>
-''README''를 보면 설치에 대한 전반적이 설명이 있다. 우리는 HEPMC포맷으로 이벤트를 만들고 싶으므로, configure에 추가 옵션을 주도록 한다.
-<code>
-./configure --with-hepmc2=/home/매드그래프_경로/HEPTools/hepmc
-gmake
-</code>
-이제 필요한 event를 만들기 위해 ''example''폴더에 있는 코드 중에서 하나를 골라야한다. CMSSW에서는 main89를 사용하는 것으로 보인다. [[https://cmssdt.cern.ch/lxr/source/GeneratorInterface/Pythia8Interface/plugins/Pythia8Hadronizer.cc|ref]]
-CMS CP5 Tune를 맞추기 위해 다음과 같이 ww.cmnd를 작성한다. ''!''로 시작하는 행은 주석이다.
-<code>
-! File: main89.cmnd
-Main:numberOfEvents = 10000          ! number of events to generate
-! Settings related to output in init(), next() and stat().
-Init:showChangedSettings = on      ! list changed settings
-Init:showAllSettings = off         ! list all settings
-Init:showChangedParticleData = on  ! list changed particle data
-Init:showAllParticleData = off     ! list all particle data
-Next:numberCount = 1000            ! print message every n events
-Next:numberShowLHA = 1             ! print LHA information n times
-Next:numberShowInfo = 1            ! print event information n times
-Next:numberShowProcess = 1         ! print process record n times
-Next:numberShowEvent = 1           ! print event record n times
-Stat:showPartonLevel = on          ! additional statistics on MPI
-! Beam parameter settings. Values below agree with default ones.
-Beams:idA = 2212                   ! first beam, p = 2212, pbar = -2212
-Beams:idB = 2212                   ! second beam, p = 2212, pbar = -2212
-Beams:eCM = 13000.                 ! CM energy of collision
-WeakDoubleBoson:ffbar2WW=on
-Tune:pp=14
-Tune:ee=7
-PDF:pSet=20
-MultipartonInteractions:ecmPow=0.03344
-MultipartonInteractions:bProfile=2
-MultipartonInteractions:pT0Ref=1.41
-MultipartonInteractions:coreRadius=0.7634
-MultipartonInteractions:coreFraction=0.63
-ColourReconnection:range=5.176
-SigmaTotal:zeroAXB=off
-SpaceShower:alphaSorder=2
-SpaceShower:alphaSvalue=0.118
-SigmaProcess:alphaSvalue=0.118
-SigmaProcess:alphaSorder=2
-MultipartonInteractions:alphaSvalue=0.118
-MultipartonInteractions:alphaSorder=2
-TimeShower:alphaSorder=2
-TimeShower:alphaSvalue=0.118
-Tune:preferLHAPDF = 2
-Main:timesAllowErrors = 10000
-Check:epTolErr = 0.01
-Beams:setProductionScalesFromLHEF = off
-SLHA:keepSM = on
-SLHA:minMassSM = 1000.
-ParticleDecays:limitTau0 = on
-ParticleDecays:tau0Max = 10
-ParticleDecays:allowPhotonRadiation = on
-</code>
-<code>
-make main89
-./main89 ww.cmnd ww_001.hepmc > log_001.txt
-</code>
-PDF setting을 위해 다음이 (아마도) ''.bashrc''에 포함되어야 할 것이다
-<code>
-export ROOTSYS=/매드그래프_경로
-export PYTHONPATH=$ROOTSYS/HEPTools/lhapdf6/lib64/python2.7/site-packages:$PYTHONPATH
-export LD_LIBRARY_PATH=$ROOTSYS/HEPTools/lhapdf6/lib:$LD_LIBRARY_PATH
-</code>
-추가한 후에는 ''source ~/.bashrc''를 잊지 말자.
-WW의 경우 ''WeakDoubleBoson:ffbar2WW=on'', WZ는 ''WeakDoubleBoson:ffbar2ZW=on'', ZZ는 ''WeakDoubleBoson:ffbar2gmZgmZ=on''를 이용한다.