Madgraph를 이용해 대부분의 이벤트를 만들 수 있지만, 가끔 복잡한 경우가 있다. (예를 들어 VV: WW, WZ, ZZ) 특수한 경우 Pythia8에서 자체적으로 decay mode를 제공하여 편리하게 이벤트를 만들 수 있다. VV의 경우 CMS simulation도 pythia를 이용해 event를 만들고 있다. /MG5_aMC_v2_6_4/HEPTools/pythia8/share/Pythia8/examples에 필요한 파일들이 있긴 하지만 여기서는 별도로 Pythia를 설치하여 이용해보기로 한다.

cd work #some working folder
wget http://home.thep.lu.se/~torbjorn/pythia8/pythia8235.tgz
tar -xvf pythia8235.tgz
cd pythia8235

README를 보면 설치에 대한 전반적이 설명이 있다. 우리는 HEPMC포맷으로 이벤트를 만들고 싶으므로, configure에 추가 옵션을 주도록 한다.

./configure --with-hepmc2=/home/매드그래프_경로/HEPTools/hepmc
gmake

이제 필요한 event를 만들기 위해 example폴더에 있는 코드 중에서 하나를 골라야한다. CMSSW에서는 main89를 사용하는 것으로 보인다. ref

CMS CP5 Tune를 맞추기 위해 다음과 같이 ww.cmnd를 작성한다. !로 시작하는 행은 주석이다.

! 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

make main89
./main89 ww.cmnd ww_001.hepmc > log_001.txt 

PDF setting을 위해 다음이 (아마도) .bashrc에 포함되어야 할 것이다

export ROOTSYS=/매드그래프_경로
export PYTHONPATH=$ROOTSYS/HEPTools/lhapdf6/lib64/python2.7/site-packages:$PYTHONPATH
export LD_LIBRARY_PATH=$ROOTSYS/HEPTools/lhapdf6/lib:$LD_LIBRARY_PATH

추가한 후에는 source ~/.bashrc를 잊지 말자.

WW의 경우 WeakDoubleBoson:ffbar2WW=on, WZ는 WeakDoubleBoson:ffbar2ZW=on, ZZ는 WeakDoubleBoson:ffbar2gmZgmZ=on를 이용한다.