-- ShionChen - 2015-01-05

基本的なsignal generationの方法

1. 適当なathena releaseをセット。例えば、

$ asetup 17.2.13.6,slc5,64

2. jobOption fileを取得

https://svnweb.cern.ch/trac/atlasoff/browser/Generators/MC12JobOptions/trunk 

例えば、

https://svnweb.cern.ch/trac/atlasoff/browser/Generators/MC12JobOptions/trunk/share/DSID178xxx/MC12.178223.Herwigpp_UEEE3_CTEQ6L1_direct_C1N2_200_N1_100_WZ_3lep.py

後は、Herwig++の場合はSLHA fileが必要。これの場合、

https://svnweb.cern.ch/trac/atlasoff/browser/Generators/MC12JobOptions/trunk/share/DSID178xxx/susy_direct_C1N2_200_N1_100_WZ.slha

3. 上を置いたdirectoryで、下記のコマンドを実行。

$ Generate_trf.py ecmEnergy=14000 firstEvent=1 randomSeed=1 jobConfig= [[https://svnweb.cern.ch/trac/atlasoff/browser/Generators/MC12JobOptions/trunk/share/DSID178xxx/MC12.178223.Herwigpp_UEEE3_CTEQ6L1_direct_C1N2_200_N1_100_WZ_3lep.py][MC12.178223.Herwigpp_UEEE3_CTEQ6L1_direct_C1N2_200_N1_100_WZ_3lep.py]] outputEVNTFile=EVNT.root evgenJobOpts=MC12JobOpts-00-12-75_v0.tar.gz runNumber=178223 maxEvents=-1

4. EVNT.rootができるので、これをinputにして、NTUP_TRUTHを作る。

今までp1032 tagのNTUP_TRUTHを使っていたようなので、このtagを作りたい場合は、asetup 17.2.2.1.2 を使う。

$ Reco_trf.py inputEVNTFile=EVNT.root outputNTUP_TRUTHFile=NTUP_TRUTH.root

Les Houches Accord format

Ref: hep-ph/0109068

■ Output columns

Pasted Graphic 1.pdf

3rd/4th columns: parents

5th/6th columns: color flow

Pasted Graphic 2.pdf

Parton showering by Pythia

MadGraph の場合は、

SMなら、

include ( "MC12JobOptions/Pythia_AUET2B_CTEQ6L1_Common.py" )

topAlg.Pythia.PythiaCommand += [ "pyinit user madgraph",

"pystat 1 3 4 5",

"pyinit dumpr 1 5",

"pyinit pylistf 1",

]

evgenConfig.inputfilecheck = 'madgraph.300410.WZ'

evgenConfig.generators += ["MadGraph", "Pythia"]

evgenConfig.description = 'WZ'

evgenConfig.keywords = ["WZ"]

evgenConfig.contact = ["somebody"]

phojf=open('./pythia_card.dat', 'w')

phojinp = """

!...Matching parameters...

IEXCFILE=0

"""%(runArgs)

phojf.write(phojinp)

phojf.close()

SUSYなら、

include ( "MC12JobOptions/Pythia_AUET2B_CTEQ6L1_Common.py" )

topAlg.Pythia.PythiaCommand += [ "pyinit user madgraph",

"pystat 1 3 4 5",

"pyinit dumpr 1 5",

"pyinit pylistf 1",

]

evgenConfig.inputfilecheck = 'madgraph.300411.C1N2'

evgenConfig.generators += ["MadGraph", "Pythia"]

evgenConfig.description = 'C1N2'

evgenConfig.keywords = ["SUSY"]

evgenConfig.contact = ["somebody"]

phojf=open('./pythia_card.dat', 'w')

phojinp = """

!...Matching parameters...

IEXCFILE=0

imss(21)=24

imss(22)=24

"""%(runArgs)

phojf.write(phojinp)

phojf.close()

で、inputGeneratorFileは、MadGraphで作ったLHE fileを

group.phys-gener.madgraph.300410.WZ_13TeV.TXT.mc12_v1._00001.events

みたいに名づけて、

group.phys-gener.madgraph.300410.WZ_13TeV.TXT.mc12_v1._00001.tar.gz

とすればいいはず。

Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2015-01-05 - ShionChen
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback