Show pageOld revisionsBacklinksBack to top This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong. ==== 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''를 이용한다. tutorial/pythia.txt Last modified: 2018/12/17 06:33by minerva1993