|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
|
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
SLHA_TABLE = """
Block MODSEL # Model selection
1 0 # Generic MSSM
Block MINPAR # Input parameters
3 1.00000000E+00 # tanb at m_Z
#
Block SMINPUTS # SM parameters
1 1.27931277E+02 # alpha_em^-1(MZ)^MSbar
2 1.16639000E-05 # G_mu [GeV^-2]
3 1.17200000E-01 # alpha_s(MZ)^MSbar
4 9.11876000E+01 # m_Z(pole)
5 4.20000000E+00 # m_b(m_b), MSbar
6 1.74300000E+02 # m_t(pole)
7 1.77700000E+00 # m_tau(pole)
Block MASS # Mass spectrum
# PDG code mass particle
24 8.04009772E+01 # W+
25 10.0E+4 # h0
35 10.0E+4 # H0
36 10.0E+4 # A0
37 10.0E+4 # H+
1000001 10.0E+4 # ~d_L
2000001 10.0E+4 # ~d_R
1000002 10.0E+4 # ~u_L
2000002 10.0E+4 # ~u_R
1000003 10.0E+4 # ~s_L
2000003 10.0E+4 # ~s_R
1000004 10.0E+4 # ~c_L
2000004 10.0E+4 # ~c_R
1000005 10.0E+4 # ~b_1
2000005 10.0E+4 # ~b_2
1000006 10.0E+4 # ~t_1
2000006 10.0E+4 # ~t_2
1000011 10.0E+4 # ~e_L-
2000011 10.0E+4 # ~e_R-
1000012 10.0E+4 # ~nu_eL
1000013 10.0E+4 # ~mu_L-
2000013 10.0E+4 # ~mu_R-
1000014 10.0E+4 # ~nu_muL
1000015 10.0E+4 # ~tau_1-
2000015 10.0E+4 # ~tau_2-
1000016 10.0E+4 # ~nu_tauL
1000021 1.50E+3 # ~g
1000022 1.00E+2 # ~chi_10
1000023 10.0E+4 # ~chi_20
1000025 10.0E+4 # ~chi_30
1000035 10.0E+4 # ~chi_40
1000024 10.0E+4 # ~chi_1+
1000037 10.0E+4 # ~chi_2+
#
BLOCK NMIX # Neutralino Mixing Matrix
1 1 9.79183656E-01 # N_11
1 2 -8.70017948E-02 # N_12
1 3 1.75813037E-01 # N_13
1 4 -5.21520034E-02 # N_14
2 1 1.39174513E-01 # N_21
2 2 9.44472080E-01 # N_22
2 3 -2.71658234E-01 # N_23
2 4 1.21674770E-01 # N_24
3 1 -7.50233573E-02 # N_31
3 2 1.16844446E-01 # N_32
3 3 6.87186106E-01 # N_33
3 4 7.13087741E-01 # N_34
4 1 -1.27284400E-01 # N_41
4 2 2.94534470E-01 # N_42
4 3 6.50435881E-01 # N_43
4 4 -6.88462993E-01 # N_44
#
BLOCK UMIX # Chargino Mixing Matrix U
1 1 9.15480281E-01 # U_11
1 2 -4.02362840E-01 # U_12
2 1 4.02362840E-01 # U_21
2 2 9.15480281E-01 # U_22
#
BLOCK VMIX # Chargino Mixing Matrix V
1 1 9.82636204E-01 # V_11
1 2 -1.85542692E-01 # V_12
2 1 1.85542692E-01 # V_21
2 2 9.82636204E-01 # V_22
#
BLOCK STOPMIX # Stop Mixing Matrix
1 1 1.0 # cos(theta_t)
1 2 0.0 # sin(theta_t)
2 1 0.0 # -sin(theta_t)
2 2 1.0 # cos(theta_t)
#
BLOCK SBOTMIX # Sbottom Mixing Matrix
1 1 9.66726392E-01 # cos(theta_b)
1 2 2.55812594E-01 # sin(theta_b)
2 1 -2.55812594E-01 # -sin(theta_b)
2 2 9.66726392E-01 # cos(theta_b)
#
BLOCK STAUMIX # Stau Mixing Matrix
1 1 4.51419848E-01 # cos(theta_tau)
1 2 8.92311672E-01 # sin(theta_tau)
2 1 -8.92311672E-01 # -sin(theta_tau)
2 2 4.51419848E-01 # cos(theta_tau)
#
BLOCK ALPHA # Higgs mixing
-1.13676047E-01 # Mixing angle in the neutral Higgs boson sector
#
BLOCK HMIX Q= 2.90528802E+02 # DRbar Higgs Parameters
1 3.05599351E+02 # mu(Q)MSSM
#
BLOCK AU Q= 2.90528802E+02 # The trilinear couplings
1 1 0.00000000E+00 # A_u(Q) DRbar
2 2 0.00000000E+00 # A_c(Q) DRbar
3 3 -4.46245994E+02 # A_t(Q) DRbar
#
BLOCK AD Q= 2.90528802E+02 # The trilinear couplings
1 1 0.00000000E+00 # A_d(Q) DRbar
2 2 0.00000000E+00 # A_s(Q) DRbar
3 3 -8.28806503E+02 # A_b(Q) DRbar
#
BLOCK AE Q= 2.90528802E+02 # The trilinear couplings
1 1 0.00000000E+00 # A_e(Q) DRbar
2 2 0.00000000E+00 # A_mu(Q) DRbar
3 3 -4.92306701E+02 # A_tau(Q) DRbar
#
BLOCK MSOFT Q= 2.90528802E+02 # The soft SUSY breaking masses at the scale Q
1 6.39136864E+01 # M_1(Q)
2 1.22006983E+02 # M_2(Q)
3 3.90619532E+02 # M_3(Q)
21 4.42860395E+04 # mH1^2(Q)
22 -9.76585434E+04 # mH2^2(Q)
31 2.26648170E+02 # meL(Q)
32 2.26648170E+02 # mmuL(Q)
33 2.24355944E+02 # mtauL(Q)
34 2.08394096E+02 # meR(Q)
35 2.08394096E+02 # mmuR(Q)
36 2.03337218E+02 # mtauR(Q)
41 4.08594291E+02 # mqL1(Q)
42 4.08594291E+02 # mqL2(Q)
43 3.46134575E+02 # mqL3(Q)
44 3.98943379E+02 # muR(Q)
45 3.98943379E+02 # mcR(Q)
46 2.58021672E+02 # mtR(Q)
47 3.95211849E+02 # mdR(Q)
48 3.95211849E+02 # msR(Q)
49 3.90320031E+02 # mbR(Q)
#
#
#
# =================
# |The decay table|
# =================
#
# - The QCD corrections to the decays gluino -> squark + quark
# squark -> gaugino + quark_prime
# squark -> squark_prime + Higgs
# squark -> gluino + quark
# are included.
#
# - The multi-body decays for the inos, stops and sbottoms are included.
#
# - The loop induced decays for the gluino, neutralinos and stops
# are included.
#
# - The SUSY decays of the top quark are included.
#
#
#
# PDG Width
DECAY 1000022 0.00000000E+00 # neutralino1 decays
DECAY 1000021 1.00000000E+00 # gluino decays
# BR NDA ID1 ID2
0.000000E+00 3 1000022 1 -1
0.000000E+00 3 1000022 2 -2
0.000000E+00 3 1000022 3 -3
0.000000E+00 3 1000022 4 -4
0.000000E+00 3 1000022 5 -5
1.000000E+00 3 1000022 6 -6
#
# PDG Width
DECAY 1000006 0.00000000E+00 # stop1 decays
DECAY 2000006 0.00000000E+00 # stop2 decays
DECAY 1000005 0.00000000E+00 # sbottom1 decays
DECAY 2000005 0.00000000E+00 # sbottom2 decays
#
# PDG Width
DECAY 1000002 1.00000000E+00 # sup_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 2 # BR(~u_L -> ~chi_10 u)
#
# PDG Width
DECAY 2000002 1.00000000E+00 # sup_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 2 # BR(~u_R -> ~chi_10 u)
#
# PDG Width
DECAY 1000001 1.00000000E+00 # sdown_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 1 # BR(~d_L -> ~chi_10 d)
#
# PDG Width
DECAY 2000001 1.00000000E+00 # sdown_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 1 # BR(~d_R -> ~chi_10 d)
#
# PDG Width
DECAY 1000004 1.00000000E+00 # scharm_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 4 # BR(~c_L -> ~chi_10 c)
#
# PDG Width
DECAY 2000004 1.00000000E+00 # scharm_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 4 # BR(~c_R -> ~chi_10 c)
#
# PDG Width
DECAY 1000003 1.00000000E+00 # sstrange_L decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 3 # BR(~s_L -> ~chi_10 s)
#
# PDG Width
DECAY 2000003 1.00000000E+00 # sstrange_R decays
# BR NDA ID1 ID2
1.00000000E+00 2 1000022 3 # BR(~s_R -> ~chi_10 s)
#
# PDG Width
DECAY 1000011 0.00000000E+00 # selectron_L decays
DECAY 2000011 0.00000000E+00 # selectron_R decays
DECAY 1000013 0.00000000E+00 # smuon_L decays
DECAY 2000013 0.00000000E+00 # smuon_R decays
DECAY 1000015 0.00000000E+00 # stau_1 decays
DECAY 2000015 0.00000000E+00 # stau_2 decays
#
# PDG Width
DECAY 1000012 0.00000000E+00 # snu_elL decays
DECAY 1000014 0.00000000E+00 # snu_muL decays
DECAY 1000016 0.00000000E+00 # snu_tauL decays
# PDG Width
DECAY 1000024 2.02592183E-05 # chargino1+ decays
# BR NDA ID1 ID2 ID3
0.33333333E+00 3 1000022 -11 12 # BR(~chi_1+ -> ~chi_10 e+ nu_e)
0.33333333E+00 3 1000022 -13 14 # BR(~chi_1+ -> ~chi_10 mu+ nu_mu)
0.33333333E+00 3 1000022 -15 16 # BR(~chi_1+ -> ~chi_10 tau+ nu_tau)
# PDG Width
DECAY 1000037 0.00000000E+00 # chargino2+ decays
DECAY 1000023 0.00000000E+00 # neutralino2 decays
DECAY 1000025 0.00000000E+00 # neutralino3 decays
DECAY 1000035 0.00000000E+00 # neutralino4 decays
"""
generator = cms.EDFilter("Pythia8ConcurrentGeneratorFilter",
pythiaPylistVerbosity = cms.untracked.int32(0),
filterEfficiency = cms.untracked.double(1.0),
pythiaHepMCVerbosity = cms.untracked.bool(False),
comEnergy = cms.double(14000),
SLHATableForPythia8 = cms.string(SLHA_TABLE),
maxEventsToPrint = cms.untracked.int32(0),
PythiaParameters = cms.PSet(
pythia8CommonSettingsBlock,
pythia8CP5SettingsBlock,
processParameters = cms.vstring(
'SUSY:all = off',
'SUSY:gg2gluinogluino = on',
'SUSY:qqbar2gluinogluino = on',
),
parameterSets = cms.vstring(
'pythia8CommonSettings',
'pythia8CP5Settings',
'processParameters'
)
)
)
|