==== 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''를 이용한다.