ATLAS Offline Software
jobOptions.WhizardPythia.py
Go to the documentation of this file.
1 
13 
14 # IMPORTANT!!! The following line adds the correct meta data to Pythia output!
15 import AthenaCommon.AtlasUnixGeneratorJob
16 #--------------------------------------------------------------
17 
18 # Private Application Configuration option
19 from AthenaCommon.AppMgr import ServiceMgr
20 ServiceMgr.MessageSvc.OutputLevel = DEBUG
21 
22 #--------------------------------------------------------------
23 # Event related parameters
24 #--------------------------------------------------------------
25 # Number of events to be processed (default is 10)
26 theApp.EvtMax = 5
27 #--------------------------------------------------------------
28 # Algorithms Private Options
29 #--------------------------------------------------------------
30 from AthenaCommon.AlgSequence import AlgSequence
31 job=AlgSequence()
32 from Pythia_i.Pythia_iConf import Pythia
33 job +=Pythia()
34 
35 # Channel switches
36 Wtojj = "0"
37 Wtoln = "1"
38 Ztojj = "0"
39 Ztoll = "1"
40 #---------------------------------------------------------------
41 job.Pythia.PythiaCommand = [
42  # Use Madgraph interface
43  "pyinit user lhef",
44 
45  # Set process type and W' mass
46  "pysubs msel 0",
47 
48  "pydat2 pmas 25 1 1.E9", # Higgs mass
49 
50 
54 
55  # Z -> to quarks
56  "pydat3 mdme 174 1 "+Ztojj,
57  "pydat3 mdme 175 1 "+Ztojj,
58  "pydat3 mdme 176 1 "+Ztojj,
59  "pydat3 mdme 177 1 "+Ztojj,
60  "pydat3 mdme 178 1 "+Ztojj,
61  "pydat3 mdme 179 1 "+Ztojj,
62 
63  # Z -> ll
64  "pydat3 mdme 182 1 "+Ztoll,
65  "pydat3 mdme 183 1 0", # to nu_e
66  "pydat3 mdme 184 1 "+Ztoll,
67  "pydat3 mdme 185 1 0", # to nu_mu
68 
69  # Close Z -> tau channels
70  "pydat3 mdme 186 1 0",
71  "pydat3 mdme 187 1 0",
72  "pydat3 mdme 188 1 0",
73  "pydat3 mdme 189 1 0",
74 
75 
76 
78  "pydat3 mdme 206 1 "+Wtoln,
79  "pydat3 mdme 207 1 "+Wtoln,
80 
81  # Close W -> tau channels
82  "pydat3 mdme 208 1 0", # no tau
83  "pydat3 mdme 209 1 0",
84 
85 
86  # W -> quarks channels -------------------
87  "pydat3 mdme 190 1 "+Wtojj,
88  "pydat3 mdme 191 1 "+Wtojj,
89  "pydat3 mdme 192 1 "+Wtojj,
90  # "pydat3 mdme 193 1 "+Wtojj, # to t' : keep default value
91  "pydat3 mdme 194 1 "+Wtojj,
92  "pydat3 mdme 195 1 "+Wtojj,
93  "pydat3 mdme 196 1 "+Wtojj,
94  #"pydat3 mdme 197 1 "+Wtojj, # to t' : keep default value
95  "pydat3 mdme 198 1 "+Wtojj,
96  "pydat3 mdme 199 1 "+Wtojj,
97  "pydat3 mdme 200 1 "+Wtojj,
98 
103 
104  # No tau decays (Tauola)
105  "pydat3 mdcy 15 1 0",
106  # No FSR (Photos)
107  "pydat1 parj 90 20000",
108 
109 
110  "pysubs ckin 1 300", # minimun for the invariant mass
111  # of WZ
112  "pyinit pylistf 1",
113 
114  # list 1st 10 events
115  "pyinit dumpr 1 10",
116 
117  # print the decay info
118  "pyinit pylisti 12"
119 
120 ]
121 
122 #---------------------------------------------------------------
123 # Pool Persistency
124 #---------------------------------------------------------------
125 from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream
126 Stream1 = AthenaPoolOutputStream( "StreamEVGEN" )
127 Stream1.OutputFile = "pythia.pool.root"
128 Stream1.ItemList += [ 'EventInfo#*', 'McEventCollection#*' ]
129 
python.AlgSequence.AlgSequence
AlgSequence
Definition: PhysicsAnalysis/D3PDTools/AnaAlgorithm/python/AlgSequence.py:7
jobOptions.Pythia
Pythia
Definition: jobOptions.LhefMadGraphPythia.py:50
python.CreateOutputStreams.AthenaPoolOutputStream
def AthenaPoolOutputStream
backward compat
Definition: CreateOutputStreams.py:144