Using MadGraph/MadEvent/PYTHIA/PGS Ian-Woo Kim Seoul National University SNU, Mar 20, 2007 Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 1 / 18
Introduction MadGraph/MadEvent is a Feynman diagram generator, S-matrix calculator and Monte Carlo event generator. MadGraph : Feynman diagram generator, S-matrix calculator MadEvent : Event generator for hard process ( quark(parton) level ) To simulate a full LHC situation, we must have a hadron-level event generator and detector simulator. PYTHIA : Hadronization event generator ( it has more functionalities) PGS : Pretty Good (Detector) simulator A new model can be easily implemented in MadGraph : MSSM, 2HDM, etc unified framework for various models. parameter input must be set by another program. For SUSY, for example, SoftSUSY, SUSY-HIT ( SDECAY,HDECAY) Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 2 / 18
Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 3 / 18
My system Intel core 2-duo (Merom) T7200 2.0 GHz Memory 2G Scientific Linux v4.4 SUSY spectrum calculator : SoftSUSY 2.0.8 SUSY Decay rate calculator : SUSY-HIT MadGraph/MadEvent/PYTHIA/PGS Analysis tool : ROOT lhco-to-root file converter: ExRootAnalysis ( provide HEP definition for ROOT) Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 4 / 18
Web interface of MadGraph/MadEvent It has a web interface at http://madgraph.hep.uiuc.edu Need an id/password. It has a summary webpage for each process. To run MadEvent+..., we must have more previlaged id/password. impossible to use it. Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 5 / 18
Installation Install Scientific Linux v4.4 : http://www.scientificlinux.org I have the CDs. If you want, contact to me. Install ROOT : http://root.cern.ch download v5.14, compile or directly install binaries. Make sure you set environment variables correctly in.bashrc or.profile Install MadGraph/MadEvent/PYTHIA/PGS From MadGraph webpage, download MadGraph V4, Pythia and PGS package. Note that all the packages must be installed in the directory MG_ME_v4.1.19 Install ExRootAnalysis : download from MadGraph webpage. ExRootLHCO... will convert.lhco file to.root file. (.lhco = LHC Olympics format ) Install SUSY-HIT, SoftSUSY : They can be installed seperately. Their output is SLHA format which will be used for the param_card.dat file for MadGraph. Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 6 / 18
Running MadGraph/MadEvent/PYTHIA/PGS You can set your own model in MadGraph. Here I will only focus on SM and MSSM. copy Template with your own directory name( e.g. MyTestDir ) Prepare for proc_card.dat and param_card.dat. For SM, param_card.dat needs not be changed from Template. for MSSM param_card.dat is just SLHA format. Note that SoftSUSY generate only spectrum part of SLHA. We need decay rate input. (Using SUSY-HIT) newprocess : generate Feynman diagrams. build up numerical analysis codes for each process. You can see such processes in index.html generate_events : generate parton-level events. Output files are in the directory Events in.lhe.gz format. If PYTHIA/PGS are installed, then it automatically run them. cards MadGraph library MadEvent.lhe.gz PYTHIA/PGS.lhco Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 7 / 18
LHC Olympics format At collider, we can only identify the following elements for a given event. - A high energy jet (b-tagging is possible) - an electron e ± - a muon µ ± - a photon - Missing transverse momentum LHC olympics format contains such identification. # typ eta phi pt jmass ntrk btag had/em dummy dummy # is increasing simply in a given event. Next event starts with 0. Type : 0 = photon 1 = electron 2 = muon 3 = hadronically-decaying tau 4 = jet 6 = missing transverse energy Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 8 / 18
LHC Olympics format eta = pseudorapidity, phi = azimuthal angle, pt = transverse momentum jmass = invariant mass ntrk = number of track btag = b-tagging had/em = ratio of hadronic vs EM energy in the calorimeter Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 9 / 18
# typ eta phi pt jmas ntrk btag had/em dum1 0 1 2 1 2-0.629 4.513 23.54 0.11 1.0 5.0 0.10 0.0 2 4 0.208 5.953 165.26 9.24 5.0 0.0 2.63 0.0 3 4 2.610 1.793 55.56 4.93 7.0 0.0 2.66 0.0 4 4 0.960 0.523 88.96 12.02 5.0 0.0 1.48 0.0 5 4-0.406 3.685 173.10 28.34 22.0 1.0 1.44 0.0 6 4-2.456 1.476 51.38 4.75 7.0 0.0 3.54 0.0 7 6 0.000 3.362 53.11 0.00 0.0 0.0 0.00 0.0 0 2 1 1 1 1.683 2.868 40.96 0.00-1.0 0.0 0.01 0.0 2 6 0.000 5.774 48.36 0.00 0.0 0.0 0.00 0.0 0 3 1 1 1-1.634 5.182 36.88 0.00-1.0 0.0 0.00 0.0 2 6 0.000 2.111 36.23 0.00 0.0 0.0 0.00 0.0 0 4 0 1 4 1.052 2.529 26.70 0.00 1.0 0.0 0.01 0.0 2 6 0.000 5.918 33.03 0.00 0.0 0.0 0.00 0.0 0 5 0 1 4 1.903 4.242 48.42 6.39 16.0 0.0 0.41 0.0 2 4 2.428 5.696 78.36 5.64 16.0 0.0 2.04 0.0 Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 10 / 18
ExRootAnalysis ExRootAnalysis for interpreting.lhco format. run ExRootLHCOlympicsConverter, then obtain.root format file..root is a snapshot of ROOT object instances. an object LHCO which is an instance of TTree class. To browse LHCO in ROOT, gsystem->load("lib/libexrootanalysis.so"); TFile::Open("pgs_events.root"); LHCO->StartViewer(); Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 11 / 18
ExRootAnalysis LHCO contains several branches: Event, Photon, Electron, Muon, Tau, Jet, MisssingET They are TRootEvent, TRootPhoton, TRootElectron, TRootMuon, TRootTau, TRootJet, TRootMissingET class objects, respectively. Each object has data defined in ExRootAnalysis/doc/RotTreeDescription.html Using TTreeViewer, one can easily draw a histogram by drag and dropping. Use Draw with fields, cut and option, you can simply make a histogram as you want. For detail of histogram option, refer to ROOT User s guide. Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 12 / 18
Input files for MadGraph Now, we know how to deal with output fils, let s understand the input files for MadGraph. Every input file in Cards. proc_card.dat specify the model and process we analyze. To change models, just change the line in # Begin MODEL # This is TAG. Do not modify this line sm # End MODEL # This is TAG. Do not modify this line change sm to mssm, then we can analyze MSSM. Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 13 / 18
To specify processes, modify this # Begin PROCESS # This is TAG. Do not modify this line pp>e-ve~ @0 # First Process QCD=99 # Max QCD couplings QED=2 # Max QED couplings end_coup # End the couplings input pp>e-ve~j @1 QCD=99 QED=2 end_coup # Second Process # Max QCD couplings # Max QED couplings # End the couplings input pp>tt~ @2 # Third Process QCD=99 # Max QCD couplings end_coup # End the couplings input done # this tells MG there are no more procs # End PROCESS # This is TAG. Do not modify this line Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 14 / 18
For model description, we can refer to MG_ME_V4.1.19/Models/(model name) particles.dat contains the particle information. For SM, particle names are d, d d-quark u, u u-quark s, s s-quark c, c c-quark b, b b-quark t, t t-quark e-, e+ electron ve, ve e-neutrino mu-,mu+ muon vm, vm µ-neutrino ta-,ta+ tau vt, vt τ-neutrino g gluon z Z-boson a photon W-,W+ W-boson h higgs interaction.dat contains the interactions. Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 15 / 18
param_card.dat : LesHouches Accord format. This must contain mass spectrum and decay rates. For MSSM, use SUSY-HIT ( SUspect-SdecaY-Hdecay-InTerface ) SUSY-HIT has inputs from susyhit.in and suspect2.in It generates SLHA format file susyhit_slha.out It can directly be renamed to param_card.dat MadGraph has another input cards for Pythia and PGS. Not examined thoroughly yet. Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 16 / 18
After newprocess, MadGraph generates the info webpage. index.html For collective particle like protons and jets, in and out particle is one of the components, so a process is divided into several subprocesses. MadGraph generates a MadEvent code madevent.tar.gz After running generate_events, at the homepage, we can access data : Parton-level( LHE format ), Hadron-level(Pythia, STDHEP format), Recognized Objects( LHCO format ). Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 17 / 18
ENJOY YOUR MONTE CARLO SIMULATION! Ian-Woo Kim (SNU) MadGraph/MadEvent/PYTHIA/PGS SNU 03/2007 18 / 18