madgraph学习小结

MADANALYSIS学习总结K7V2017.2.22目录•1.WhatisMadanalysis•2.Madanalysis5comeswithtwomodes•3.Madanalysis的一般使用•4.Madanalysis深入研究•5.Anexample•6.Conclusions1.WhatisMadanalysisAframeworkforphenomenologicalinvestigationsatparticlecolliders.BasedonaC++kernel.Madanalysisallowsustoefficientlyperform,inastraightforwardanduser-friendlyfashion,sophisticatedphysicsanalysesofeventfiles.2.TwomodesofMadanalysis•1.Thefirstone,easiertohandle,usesthestrengthsofapowerfulPythoninterfaceinordertoimplementphysicsanalysesbymeansofasetofintuitivecommands.•2.ThesecondonerequiresonetoimplementtheanalysesintheC++programminglanguage,directlywithinthecoreoftheanalysisframework.3.Madanalysis的一般使用installsamplesttbar_sl_1.lhe.gzttbar_sl_2.lhe.gzttbar_fh.lhe.gzzz.lhe.gzbin/ma5(一旦开始，文件就被保存在lib文件夹中）display_particlesdisplay_multiparticlesdefinemu=mu+mu-（定义新的粒子）displaybdisplayl+Theparticle’b’isdefinedbythePDG-id5.Themultiparticle’l+’isdefinedbythePDG-ids-11-13.//输入事例：importsamples/*.lhe.gz该命令等价于importsamples/ttbar_sl_1.lhe.gzimportsamples/ttbar_sl_2.lhe.gzimportsamples/ttbar_fh.lhe.gzimportsamples/zz.lhe.gz//CUT的选择以及画图的命令plotMETplotPT(mu)200100rejectMET100reject(mu)PT20plotM(mu+mu-)200100submitdirname//展示结果generate_latexdirnamegenerate_pdflatexdirnamegenerate_htmldirname3.Madanalysis的一般使用bin/ma5[level]bin/ma5--recolevelfilenamebin/ma5--script--recolevelfilenamebin/ma5-s--recolevelfilenamebin/ma5-s-Rfilename1filename2filename3bin/ma5–expertbin/ma5-ebin/ma5–Ebin/ma5-vbin/ma5--versionbin/ma5–releasebin/ma5-fbin/ma5–forcedbin/ma5--helpbin/ma5-h开始程序3.Madanalysis的一般使用3.Madanalysis的一般使用bin/ma5[options]help[action]（存储在.ma5history文件夹中）displayobjectdisplayobject.optionsetobject.option=xremoveobject（analysis正在运行的对象main,hadronicandinvisible是不能被删除的）!command等价于shellcommandreset3.Madanalysis的一般使用installsamplesimportpath-to-sample（输入MCSample）importdirectory/*（输入存储在directory中的Sample中，这里用*或者？表示输入一系列的Samples）对应的4种文件格式.lheor.lhe.gz.hepor.hep.gz.hepmcor.hepmc.gz.lhcoor.lhco.gz.importpath-to-signal-events/*assignalsetimportpath-to-background-events/*asbackgroundsetimportpath-to-events/ttbar*asttbarimportpath-to-events/zz*aszzimportpath-to-events/drellyan*asdrellyansetdataset.type=backgroundsetdataset.type=signal(默认的情况)setdataset.weight=weightsetdataset.xsection=value（包括linestyle,linewidth,linecolor,backstyleandbackcolor）setdataset1.backcolor=red+1setdataset2.backcolor=red+2setdataset3.backcolor=red+3setdataset4.backcolor=red+4setdataset.linestyle=valuesetdataset.linewidth=valuesetdataset.backstyle=valuesetdataset.title=stringdisplayttbarNameofthedataset=ttbar(signal)Title=’ttbar’User-imposedcrosssection=0.0User-imposedweightoftheset=1.0Linecolorinhistograms=autoLinestyleinhistograms=solidLinewidthinhistograms=1Backgroundcolorinhistograms=autoBackgroundstyleinhistograms=solidListofeventfilesincludedinthisdataset:-ttbar.lhe.gz数据设置程序3.Madanalysis的一般使用Particlesandmultiparticlesimportpath-to-UFO-files（输入UFO文件）display_particlesdisplay_multiparticlesdisplaylabel（thePDG-idlinkedtothelabellabel）definelabel=identifiers（Thecommandabovecreatesalabeldenotedbylabelandassociatestoitthecontentoftheparameteridentifiers.Ifthevalueofidentifiersconsistsinonesingleintegernumber,anewparticlelabelrelatedtothecorrespondingPDGidiscreated.）definelabel=labelidentifiersremovelabel3.Madanalysis的一般使用plotobservablenbinsminmaxdisplayselectiondisplayselection[i]removeselection[i]swapselection[i]selection[j]（修改selection）setselection[i].xmax=100setselection[i].logY=truesetselection[j].logX=falsesetselection[i].stacking_method=valuesetselection[i].titleX=stringsetselection[i].titleY=string’plotobservablenbinsminmax[options]CreatinghistogramsplotTHT[logYsuperimpose]plotSQRTSplotNPIDplotNAPIDplotobservable(label)plotPT(mu+)plotDELTAR(label1,label2)definemulti=particle1particle2plotobservable(multi)plotobservable(label[i])setselection[i].rank=valueplotobservable(type1type2)plotobservable(type1type2type3)setselection[i].statuscode=interstateplotobservable(prtcl1prtcl2)plotvobservable(prtcl1prtcl2)3.Madanalysis的一般使用plotM(mu+mu-)plotobservable(multimulti)plotsobservable(prtcl1...)（s表示sum，先计算出每个粒子的观测量，然后求和）plotobservable(allprtcl)plotdvobservable(prtcl1prtcl2)plotvdobservable(prtcl1prtcl2)plotdobservable(prtcl1prtcl2)plotdsobservable(prtcl1prtcl2)plotsdobservable(prtcl1prtcl2)Plotrobservable(prtcl1prtcl2)plotobservable(p1p2andp3p4)（1和2求和，3和4求和，再总的求和）plotobservable(label)可观测的矢量之差可观测的标量之差3.Madanalysis的一般使用3.Madanalysis的一般使用3.Madanalysis的一般使用selectcondition[options]rejectcondition[options]rejectPT(mu)50selectM(e+e-)100observablelogical-operatorvalue（condition的来源）cond1connector1cond2connector2cond3（连接词有and和or）select(particle)condition[options]reject(particle)condition[options]Selectioncuts3.Madanalysis的一般