TopEFT Gridpack Creation
Dear Madgraph experts,
I am trying to create some gridpacks for ttbar production using an eft model (dim6top: http://
I am running the process :
p p > t t~ FCNC=0 DIM6=1, (t > w+ b, w+ > wdecay wdecay), (t~ > w- b~ , w- > wdecay wdecay)
where wdecay = g u c d s u~ c~ d~ s~ e+ mu+ tau+ e- mu- tau- ve vm vt ve~ vm~ vt~
I am however running into problems. When trying to create the gridpacks i end up with:
=== Results Summary for run: run_01 tag: tag_1 ===
Cross-section : 0 +- 0 pb
Nb of events : 0
before failing with the error:
Command "generate_events " interrupted with error:
IndexError : list index out of range
and in the log file (run_01_
File "/scratch/
return self.onecmd_
File "/scratch/
return func(arg, **opt)
File "/scratch/
self.
File "/scratch/
stop = Cmd.onecmd_
File "/scratch/
return func(arg, **opt)
File "/scratch/
proc_
File "/scratch/
return super(MultiEven
File "/scratch/
max_wgt = all_wgt[-1]
IndexError: list index out of range
Looking into some of the SubProcesses logs, I can see that none of the events passed the phase space cuts and it says to loosen cuts or increase max_events. I do not know where to increase the max_events from and the phase space cuts i think are the default ones which work fine with SM ttbar production.
Cheers,
Doug
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply:
Revision history for this message
|
#1 |
Here is the full error log:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.1 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events
Traceback (most recent call last):
File "/scratch/
return self.onecmd_
File "/scratch/
return func(arg, **opt)
File "/scratch/
self.
File "/scratch/
stop = Cmd.onecmd_
File "/scratch/
return func(arg, **opt)
File "/scratch/
proc_
File "/scratch/
return super(MultiEven
File "/scratch/
max_wgt = all_wgt[-1]
IndexError: list index out of range
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.1 2017-12-12 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model dim6top_LO_UFO
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define l = l+ l- ta+ ta-
define nu = vl vl~
define wdecay = j l nu
generate p p > t t~ FCNC=0 DIM6=1, (t > w+ b, w+ > wdecay wdecay), (t~\
> w- b~ , w- > wdecay wdecay)
output process -nojpeg
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
#######
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
#######
#######
## INFORMATION FOR DIM6
#######
Block dim6
1 1.000000e+03 # Lambda
2 0.000000e+00 # ctp
3 0.000000e+00 # ctpI
4 0.000000e+00 # cpQM
5 0.000000e+00 # cpQ3
6 0.000000e+00 # cpt
7 0.000000e+00 # cpb
8 0.000000e+00 # cptb
9 0.000000e+00 # cptbI
10 0.000000e+00 # ctW
11 0.000000e+00 # ctZ
12 0.000000e+00 # ctWI
13 0.000000e+00 # ctZI
14 0.000000e+00 # cbW
15 0.000000e+00 # cbWI
16 0.000000e+00 # ctG
17 0.000000e+00 # ctGI
18 0.000000e+00 # cQlM1
19 0.000000e+00 # cQlM2
20 0.000000e+00 # cQlM3
21 0.000000e+00 # cQl31
22 0.000000e+00 # cQl32
23 0.000000e+00 # cQl33
24 0.000000e+00 # cQe1
25 0.000000e+00 # cQe2
26 0.000000e+00 # cQe3
27 0.000000e+00 # ctl1
28 0.000000e+00 # ctl2
29 0.000000e+00 # ctl3
30 0.000000e+00 # cte1
31 0.000000e+00 # cte2
32 0.000000e+00 # cte3
33 0.000000e+00 # ctlS1
34 0.000000e+00 # ctlSI1
35 0.000000e+00 # ctlS2
36 0.000000e+00 # ctlSI2
37 0.000000e+00 # ctlS3
38 0.000000e+00 # ctlSI3
39 0.000000e+00 # ctlT1
40 0.000000e+00 # ctlTI1
41 0.000000e+00 # ctlT2
42 0.000000e+00 # ctlTI2
43 0.000000e+00 # ctlT3
44 0.000000e+00 # ctlTI3
45 0.000000e+00 # cblS1
46 0.000000e+00 # cblSI1
47 0.000000e+00 # cblS2
48 0.000000e+00 # cblSI2
49 0.000000e+00 # cblS3
50 0.000000e+00 # cblSI3
51 0.000000e+00 # cQq83
52 0.000000e+00 # cQq81
53 0.000000e+00 # cQu8
54 0.000000e+00 # cQd8
55 0.000000e+00 # ctq8
56 0.000000e+00 # ctu8
57 0.000000e+00 # ctd8
58 0.000000e+00 # cQq13
59 0.000000e+00 # cQq11
60 0.000000e+00 # cQu1
61 0.000000e+00 # cQd1
62 0.000000e+00 # ctq1
63 0.000000e+00 # ctu1
64 0.000000e+00 # ctd1
65 0.000000e+00 # cQQ1
66 0.000000e+00 # cQQ8
67 0.000000e+00 # cQt1
68 0.000000e+00 # cQb1
69 0.000000e+00 # ctt1
70 0.000000e+00 # ctb1
71 0.000000e+00 # cQt8
72 0.000000e+00 # cQb8
73 0.000000e+00 # ctb8
74 0.000000e+00 # ctQqu1
75 0.000000e+00 # ctQqu1I
76 0.000000e+00 # ctQqu8
77 0.000000e+00 # ctQqu8I
78 0.000000e+00 # cbQqd1
79 0.000000e+00 # cbQqd1I
80 0.000000e+00 # cbQqd8
81 0.000000e+00 # cbQqd8I
82 0.000000e+00 # cQtqd1
83 0.000000e+00 # cQtqd1I
84 0.000000e+00 # cQtqd8
85 0.000000e+00 # cQtqd8I
86 0.000000e+00 # cQbqu1
87 0.000000e+00 # cQbqu1I
88 0.000000e+00 # cQbqu8
89 0.000000e+00 # cQbqu8I
90 0.000000e+00 # cbtud1
91 0.000000e+00 # cbtud1I
92 0.000000e+00 # cbtud8
93 0.000000e+00 # cbtud8I
94 0.000000e+00 # cQtQb1
95 0.000000e+00 # cQtQb8
96 0.000000e+00 # cQtQb1I
97 0.000000e+00 # cQtQb8I
#######
## INFORMATION FOR FCNC
#######
Block fcnc
1 0.000000e+00 # ctpx13
2 0.000000e+00 # ctpx23
3 0.000000e+00 # ctpx31
4 0.000000e+00 # ctpx32
5 0.000000e+00 # ctpIx13
6 0.000000e+00 # ctpIx23
7 0.000000e+00 # ctpIx31
8 0.000000e+00 # ctpIx32
9 0.000000e+00 # cpQMx31
10 0.000000e+00 # cpQMx32
11 0.000000e+00 # cpQMIx31
12 0.000000e+00 # cpQMIx32
13 0.000000e+00 # cpQ3x31
14 0.000000e+00 # cpQ3x32
15 0.000000e+00 # cpQ3Ix31
16 0.000000e+00 # cpQ3Ix32
17 0.000000e+00 # cptx31
18 0.000000e+00 # cptx32
19 0.000000e+00 # cptIx31
20 0.000000e+00 # cptIx32
21 0.000000e+00 # cptbx13
22 0.000000e+00 # cptbx23
23 0.000000e+00 # cptbx31
24 0.000000e+00 # cptbx32
25 0.000000e+00 # cptbIx13
26 0.000000e+00 # cptbIx23
27 0.000000e+00 # cptbIx31
28 0.000000e+00 # cptbIx32
29 0.000000e+00 # ctWx13
30 0.000000e+00 # ctWx23
31 0.000000e+00 # ctWx31
32 0.000000e+00 # ctWx32
33 0.000000e+00 # ctWIx13
34 0.000000e+00 # ctWIx23
35 0.000000e+00 # ctWIx31
36 0.000000e+00 # ctWIx32
37 0.000000e+00 # ctZx13
38 0.000000e+00 # ctZx23
39 0.000000e+00 # ctZx31
40 0.000000e+00 # ctZx32
41 0.000000e+00 # ctZIx13
42 0.000000e+00 # ctZIx23
43 0.000000e+00 # ctZIx31
44 0.000000e+00 # ctZIx32
45 0.000000e+00 # cbWx13
46 0.000000e+00 # cbWx23
47 0.000000e+00 # cbWx31
48 0.000000e+00 # cbWx32
49 0.000000e+00 # cbWIx13
50 0.000000e+00 # cbWIx23
51 0.000000e+00 # cbWIx31
52 0.000000e+00 # cbWIx32
53 0.000000e+00 # ctGx13
54 0.000000e+00 # ctGx23
55 0.000000e+00 # ctGx31
56 0.000000e+00 # ctGx32
57 0.000000e+00 # ctGIx13
58 0.000000e+00 # ctGIx23
59 0.000000e+00 # ctGIx31
60 0.000000e+00 # ctGIx32
61 0.000000e+00 # cQl3x1x31
62 0.000000e+00 # cQl3x1x32
63 0.000000e+00 # cQl3x2x31
64 0.000000e+00 # cQl3x2x32
65 0.000000e+00 # cQl3x3x31
66 0.000000e+00 # cQl3x3x32
67 0.000000e+00 # cQl3Ix1x31
68 0.000000e+00 # cQl3Ix1x32
69 0.000000e+00 # cQl3Ix2x31
70 0.000000e+00 # cQl3Ix2x32
71 0.000000e+00 # cQl3Ix3x31
72 0.000000e+00 # cQl3Ix3x32
73 0.000000e+00 # cQlMx1x31
74 0.000000e+00 # cQlMx1x32
75 0.000000e+00 # cQlMx2x31
76 0.000000e+00 # cQlMx2x32
77 0.000000e+00 # cQlMx3x31
78 0.000000e+00 # cQlMx3x32
79 0.000000e+00 # cQlMIx1x31
80 0.000000e+00 # cQlMIx1x32
81 0.000000e+00 # cQlMIx2x31
82 0.000000e+00 # cQlMIx2x32
83 0.000000e+00 # cQlMIx3x31
84 0.000000e+00 # cQlMIx3x32
85 0.000000e+00 # cQex1x31
86 0.000000e+00 # cQex1x32
87 0.000000e+00 # cQex2x31
88 0.000000e+00 # cQex2x32
89 0.000000e+00 # cQex3x31
90 0.000000e+00 # cQex3x32
91 0.000000e+00 # cQeIx1x31
92 0.000000e+00 # cQeIx1x32
93 0.000000e+00 # cQeIx2x31
94 0.000000e+00 # cQeIx2x32
95 0.000000e+00 # cQeIx3x31
96 0.000000e+00 # cQeIx3x32
97 0.000000e+00 # ctlx1x31
98 0.000000e+00 # ctlx1x32
99 0.000000e+00 # ctlx2x31
100 0.000000e+00 # ctlx2x32
101 0.000000e+00 # ctlx3x31
102 0.000000e+00 # ctlx3x32
103 0.000000e+00 # ctlIx1x31
104 0.000000e+00 # ctlIx1x32
105 0.000000e+00 # ctlIx2x31
106 0.000000e+00 # ctlIx2x32
107 0.000000e+00 # ctlIx3x31
108 0.000000e+00 # ctlIx3x32
109 0.000000e+00 # ctex1x31
110 0.000000e+00 # ctex1x32
111 0.000000e+00 # ctex2x31
112 0.000000e+00 # ctex2x32
113 0.000000e+00 # ctex3x31
114 0.000000e+00 # ctex3x32
115 0.000000e+00 # cteIx1x31
116 0.000000e+00 # cteIx1x32
117 0.000000e+00 # cteIx2x31
118 0.000000e+00 # cteIx2x32
119 0.000000e+00 # cteIx3x31
120 0.000000e+00 # cteIx3x32
121 0.000000e+00 # ctlSx1x13
122 0.000000e+00 # ctlSx1x23
123 0.000000e+00 # ctlSx1x31
124 0.000000e+00 # ctlSx1x32
125 0.000000e+00 # ctlSx2x13
126 0.000000e+00 # ctlSx2x23
127 0.000000e+00 # ctlSx2x31
128 0.000000e+00 # ctlSx2x32
129 0.000000e+00 # ctlSx3x13
130 0.000000e+00 # ctlSx3x23
131 0.000000e+00 # ctlSx3x31
132 0.000000e+00 # ctlSx3x32
133 0.000000e+00 # ctlSIx1x13
134 0.000000e+00 # ctlSIx1x23
135 0.000000e+00 # ctlSIx1x31
136 0.000000e+00 # ctlSIx1x32
137 0.000000e+00 # ctlSIx2x13
138 0.000000e+00 # ctlSIx2x23
139 0.000000e+00 # ctlSIx2x31
140 0.000000e+00 # ctlSIx2x32
141 0.000000e+00 # ctlSIx3x13
142 0.000000e+00 # ctlSIx3x23
143 0.000000e+00 # ctlSIx3x31
144 0.000000e+00 # ctlSIx3x32
145 0.000000e+00 # ctlTx1x13
146 0.000000e+00 # ctlTx1x23
147 0.000000e+00 # ctlTx1x31
148 0.000000e+00 # ctlTx1x32
149 0.000000e+00 # ctlTx2x13
150 0.000000e+00 # ctlTx2x23
151 0.000000e+00 # ctlTx2x31
152 0.000000e+00 # ctlTx2x32
153 0.000000e+00 # ctlTx3x13
154 0.000000e+00 # ctlTx3x23
155 0.000000e+00 # ctlTx3x31
156 0.000000e+00 # ctlTx3x32
157 0.000000e+00 # ctlTIx1x13
158 0.000000e+00 # ctlTIx1x23
159 0.000000e+00 # ctlTIx1x31
160 0.000000e+00 # ctlTIx1x32
161 0.000000e+00 # ctlTIx2x13
162 0.000000e+00 # ctlTIx2x23
163 0.000000e+00 # ctlTIx2x31
164 0.000000e+00 # ctlTIx2x32
165 0.000000e+00 # ctlTIx3x13
166 0.000000e+00 # ctlTIx3x23
167 0.000000e+00 # ctlTIx3x31
168 0.000000e+00 # ctlTIx3x32
169 0.000000e+00 # cqq11x3331
170 0.000000e+00 # cqq11x3332
171 0.000000e+00 # cqq11Ix3331
172 0.000000e+00 # cqq11Ix3332
173 0.000000e+00 # cqq13x3331
174 0.000000e+00 # cqq13x3332
175 0.000000e+00 # cqq13Ix3331
176 0.000000e+00 # cqq13Ix3332
177 0.000000e+00 # cuu1x3331
178 0.000000e+00 # cuu1x3332
179 0.000000e+00 # cuu1Ix3331
180 0.000000e+00 # cuu1Ix3332
181 0.000000e+00 # cqu1x3331
182 0.000000e+00 # cqu1x3332
183 0.000000e+00 # cqu1Ix3331
184 0.000000e+00 # cqu1Ix3332
185 0.000000e+00 # cqu1x3133
186 0.000000e+00 # cqu1x3233
187 0.000000e+00 # cqu1Ix3133
188 0.000000e+00 # cqu1Ix3233
189 0.000000e+00 # cqu8x3331
190 0.000000e+00 # cqu8x3332
191 0.000000e+00 # cqu8Ix3331
192 0.000000e+00 # cqu8Ix3332
193 0.000000e+00 # cqu8x3233
194 0.000000e+00 # cqu8x3133
195 0.000000e+00 # cqu8Ix3133
196 0.000000e+00 # cqu8Ix3233
197 0.000000e+00 # cqd1x3331
198 0.000000e+00 # cqd1x3332
199 0.000000e+00 # cqd1Ix3331
200 0.000000e+00 # cqd1Ix3332
201 0.000000e+00 # cqd1x3133
202 0.000000e+00 # cqd1x3233
203 0.000000e+00 # cqd1Ix3133
204 0.000000e+00 # cqd1Ix3233
205 0.000000e+00 # cqd8x3331
206 0.000000e+00 # cqd8x3332
207 0.000000e+00 # cqd8Ix3331
208 0.000000e+00 # cqd8Ix3332
209 0.000000e+00 # cqd8x3133
210 0.000000e+00 # cqd8x3233
211 0.000000e+00 # cqd8Ix3133
212 0.000000e+00 # cqd8Ix3233
213 0.000000e+00 # cud1x3331
214 0.000000e+00 # cud1x3332
215 0.000000e+00 # cud1Ix3331
216 0.000000e+00 # cud1Ix3332
217 0.000000e+00 # cud1x3133
218 0.000000e+00 # cud1x3233
219 0.000000e+00 # cud1Ix3133
220 0.000000e+00 # cud1Ix3233
221 0.000000e+00 # cud8x3331
222 0.000000e+00 # cud8x3332
223 0.000000e+00 # cud8Ix3331
224 0.000000e+00 # cud8Ix3332
225 0.000000e+00 # cud8x3133
226 0.000000e+00 # cud8x3233
227 0.000000e+00 # cud8Ix3133
228 0.000000e+00 # cud8Ix3233
229 0.000000e+00 # cquqd1x3331
230 0.000000e+00 # cquqd1x3332
231 0.000000e+00 # cquqd1x3313
232 0.000000e+00 # cquqd1x3323
233 0.000000e+00 # cquqd1x3133
234 0.000000e+00 # cquqd1x3233
235 0.000000e+00 # cquqd1x1333
236 0.000000e+00 # cquqd1x2333
237 0.000000e+00 # cquqd1Ix3331
238 0.000000e+00 # cquqd1Ix3332
239 0.000000e+00 # cquqd1Ix3313
240 0.000000e+00 # cquqd1Ix3323
241 0.000000e+00 # cquqd1Ix3133
242 0.000000e+00 # cquqd1Ix3233
243 0.000000e+00 # cquqd1Ix1333
244 0.000000e+00 # cquqd1Ix2333
245 0.000000e+00 # cquqd8x3331
246 0.000000e+00 # cquqd8x3332
247 0.000000e+00 # cquqd8x3313
248 0.000000e+00 # cquqd8x3323
249 0.000000e+00 # cquqd8x3133
250 0.000000e+00 # cquqd8x3233
251 0.000000e+00 # cquqd8x1333
252 0.000000e+00 # cquqd8x2333
253 0.000000e+00 # cquqd8Ix3331
254 0.000000e+00 # cquqd8Ix3332
255 0.000000e+00 # cquqd8Ix3313
256 0.000000e+00 # cquqd8Ix3323
257 0.000000e+00 # cquqd8Ix3133
258 0.000000e+00 # cquqd8Ix3233
259 0.000000e+00 # cquqd8Ix1333
260 0.000000e+00 # cquqd8Ix2333
261 0.000000e+00 # cqq11x31ii
262 0.000000e+00 # cqq11x32ii
263 0.000000e+00 # cqq11Ix31ii
264 0.000000e+00 # cqq11Ix32ii
265 0.000000e+00 # cqq13x31ii
266 0.000000e+00 # cqq13x32ii
267 0.000000e+00 # cqq13Ix31ii
268 0.000000e+00 # cqq13Ix32ii
269 0.000000e+00 # cqq81x31ii
270 0.000000e+00 # cqq81x32ii
271 0.000000e+00 # cqq81Ix31ii
272 0.000000e+00 # cqq81Ix32ii
273 0.000000e+00 # cqq83x31ii
274 0.000000e+00 # cqq83x32ii
275 0.000000e+00 # cqq83Ix31ii
276 0.000000e+00 # cqq83Ix32ii
277 0.000000e+00 # cuu1x31ii
278 0.000000e+00 # cuu1x32ii
279 0.000000e+00 # cuu1Ix31ii
280 0.000000e+00 # cuu1Ix32ii
281 0.000000e+00 # cuu8x31ii
282 0.000000e+00 # cuu8x32ii
283 0.000000e+00 # cuu8Ix31ii
284 0.000000e+00 # cuu8Ix32ii
285 0.000000e+00 # cud1x31ii
286 0.000000e+00 # cud1x32ii
287 0.000000e+00 # cud1Ix31ii
288 0.000000e+00 # cud1Ix32ii
289 0.000000e+00 # cud8x31ii
290 0.000000e+00 # cud8x32ii
291 0.000000e+00 # cud8Ix31ii
292 0.000000e+00 # cud8Ix32ii
293 0.000000e+00 # cqu1x31ii
294 0.000000e+00 # cqu1x32ii
295 0.000000e+00 # cqu1Ix31ii
296 0.000000e+00 # cqu1Ix32ii
297 0.000000e+00 # cqu8x31ii
298 0.000000e+00 # cqu8x32ii
299 0.000000e+00 # cqu8Ix31ii
300 0.000000e+00 # cqu8Ix32ii
301 0.000000e+00 # cqu1xii31
302 0.000000e+00 # cqu1xii32
303 0.000000e+00 # cqu1Ixii31
304 0.000000e+00 # cqu1Ixii32
305 0.000000e+00 # cqu8xii31
306 0.000000e+00 # cqu8xii32
307 0.000000e+00 # cqu8Ixii31
308 0.000000e+00 # cqu8Ixii32
309 0.000000e+00 # cqd1x31ii
310 0.000000e+00 # cqd1x32ii
311 0.000000e+00 # cqd1Ix31ii
312 0.000000e+00 # cqd1Ix32ii
313 0.000000e+00 # cqd8Ix31ii
314 0.000000e+00 # cqd8Ix32ii
315 0.000000e+00 # cqd8x31ii
316 0.000000e+00 # cqd8x32ii
317 0.000000e+00 # cpbIx31
318 0.000000e+00 # cpbIx32
319 0.000000e+00 # cpbx31
320 0.000000e+00 # cpbx32
321 0.000000e+00 # cblSx1x13
322 0.000000e+00 # cblSx1x23
323 0.000000e+00 # cblSx1x31
324 0.000000e+00 # cblSx1x32
325 0.000000e+00 # cblSx2x13
326 0.000000e+00 # cblSx2x23
327 0.000000e+00 # cblSx2x31
328 0.000000e+00 # cblSx2x32
329 0.000000e+00 # cblSx3x13
330 0.000000e+00 # cblSx3x23
331 0.000000e+00 # cblSx3x31
332 0.000000e+00 # cblSx3x32
333 0.000000e+00 # cblSIx1x13
334 0.000000e+00 # cblSIx1x23
335 0.000000e+00 # cblSIx1x31
336 0.000000e+00 # cblSIx1x32
337 0.000000e+00 # cblSIx2x13
338 0.000000e+00 # cblSIx2x23
339 0.000000e+00 # cblSIx2x31
340 0.000000e+00 # cblSIx2x32
341 0.000000e+00 # cblSIx3x13
342 0.000000e+00 # cblSIx3x23
343 0.000000e+00 # cblSIx3x31
344 0.000000e+00 # cblSIx3x32
#######
## INFORMATION FOR MASS
#######
Block mass
5 4.700000e+00 # MB
6 1.720000e+02 # MT
15 1.777000e+00 # MTA
23 9.118760e+01 # MZ
25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000 # d : 0.0
2 0.000000 # u : 0.0
3 0.000000 # s : 0.0
4 0.000000 # c : 0.0
11 0.000000 # e- : 0.0
12 0.000000 # ve : 0.0
13 0.000000 # mu- : 0.0
14 0.000000 # vm : 0.0
16 0.000000 # vt : 0.0
21 0.000000 # g : 0.0
22 0.000000 # a : 0.0
24 79.824360 # w+ : cmath.sqrt(
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.279000e+02 # aEWM1
2 1.166370e-05 # Gf
3 1.184000e-01 # aS
#######
## INFORMATION FOR YUKAWA
#######
Block yukawa
5 4.700000e+00 # ymb
6 1.720000e+02 # ymt
15 1.777000e+00 # ymtau
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.508336e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.085000e+00 # WW
DECAY 25 4.070000e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
#******
#
#******
# Running parameters
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
10000 = nevents ! Number of unweighted events requested
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
1 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
6500.0 = ebeam1 ! beam 1 total energy in GeV
6500.0 = ebeam2 ! beam 2 total energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0.0 = polbeam1 ! beam polarization for beam 1
0.0 = polbeam2 ! beam polarization for beam 2
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
nn23lo1 = pdlabel ! PDF set
230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#******
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
.true. = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
3.0 = lhe_version ! Change the way clustering information pass to shower.
True = clusinfo ! include clustering tag in output
average = event_norm ! average/sum. Normalization of the weight in the LHEF
#******
# Matching parameter (MLM only)
#******
1 = ickkw ! 0 no matching, 1 MLM
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
4 = asrwgtflavor ! highest quark flavor for a_s reweight
False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
30.0 = xqcut ! minimum kt jet measure between partons
#******
#
#******
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
False = cut_decays ! Cut decay products
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
20.0 = ptj ! minimum pt for the jets
0.0 = ptb ! minimum pt for the b
10.0 = pta ! minimum pt for the photons
10.0 = ptl ! minimum pt for the charged leptons
0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
-1.0 = ptjmax ! maximum pt for the jets
-1.0 = ptbmax ! maximum pt for the b
-1.0 = ptamax ! maximum pt for the photons
-1.0 = ptlmax ! maximum pt for the charged leptons
-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
{} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
{} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#******
# Minimum and maximum E's (in the center of mass frame) *
#******
0.0 = ej ! minimum E for the jets
0.0 = eb ! minimum E for the b
0.0 = ea ! minimum E for the photons
0.0 = el ! minimum E for the charged leptons
-1.0 = ejmax ! maximum E for the jets
-1.0 = ebmax ! maximum E for the b
-1.0 = eamax ! maximum E for the photons
-1.0 = elmax ! maximum E for the charged leptons
{} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
{} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
5.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
0.0 = etalmin ! main rap for the charged leptons
{} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
{} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#******
# Minimum and maximum DeltaR distance *
#******
0 = drjj ! min distance between jets
0.0 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.0 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0 = drjl ! min distance between jet and lepton
0.0 = drab ! min distance between gamma and b
0.0 = drbl ! min distance between b and lepton
0.4 = dral ! min distance between gamma and lepton
0 = drjjmax ! max distance between jets
-1.0 = drbbmax ! max distance between b's
-1.0 = drllmax ! max distance between leptons
-1.0 = draamax ! max distance between gammas
-1.0 = drbjmax ! max distance between b and jet
-1.0 = drajmax ! max distance between gamma and jet
0 = drjlmax ! max distance between jet and lepton
-1.0 = drabmax ! max distance between gamma and b
-1.0 = drblmax ! max distance between b and lepton
-1.0 = dralmax ! maxdistance between gamma and lepton
#******
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#******
0.0 = mmjj ! min invariant mass of a jet pair
0.0 = mmbb ! min invariant mass of a b pair
0.0 = mmaa ! min invariant mass of gamma gamma pair
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
-1.0 = mmbbmax ! max invariant mass of a b pair
-1.0 = mmaamax ! max invariant mass of gamma gamma pair
-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_
#******
# Minimum and maximum invariant mass for all letpons *
#******
0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#******
# Inclusive cuts *
#******
0.0 = ptheavy ! minimum pt for at least one heavy final state
0.0 = xptj ! minimum pt for at least one jet
0.0 = xptb ! minimum pt for at least one b
0.0 = xpta ! minimum pt for at least one photon
0.0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
0.0 = ptj3min ! minimum pt for the third jet in pt
0.0 = ptj4min ! minimum pt for the fourth jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
-1.0 = ptj3max ! maximum pt for the third jet in pt
-1.0 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#******
# Control the pt's of leptons sorted by pt *
#******
0.0 = ptl1min ! minimum pt for the leading lepton in pt
0.0 = ptl2min ! minimum pt for the second lepton in pt
0.0 = ptl3min ! minimum pt for the third lepton in pt
0.0 = ptl4min ! minimum pt for the fourth lepton in pt
-1.0 = ptl1max ! maximum pt for the leading lepton in pt
-1.0 = ptl2max ! maximum pt for the second lepton in pt
-1.0 = ptl3max ! maximum pt for the third lepton in pt
-1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
0.0 = ht4min ! minimum Ht for the four leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
-1.0 = ht4max ! maximum Ht for the four leading jets
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#******
0.0 = ptgmin ! Min photon transverse momentum
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# WBF cuts *
#******
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#******
-1.0 = ktdurham
0.4 = dparameter
-1.0 = ptlund
1, 2, 3, 4, 5, 6, 21 = pdgs_for_
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
True = use_syst ! Enable systematics studies
#
#******
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#******
#
0.5 1 2 = sys_scalefact # factorization/
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
NNPDF23_
# MSTW2008nlo68cl
#
Revision history for this message
|
#2 |
Hi,
Looks like that the model request you to do something like:
generate p p > t t~ FCNC=0 DIM6=1, (t > w+ b DIM6=0, w+ > wdecay wdecay DIM6=0), (t~ > w- b~ DIM6=0, w- > wdecay wdecay DIM6=0)
Cheers,
Olivier
> On 16 Apr 2018, at 14:13, Douglas Burns <email address hidden> wrote:
>
> New question #667989 on MadGraph5_aMC@NLO:
> https:/
>
> Dear Madgraph experts,
>
> I am trying to create some gridpacks for ttbar production using an eft model (dim6top: http://
>
> I am running the process :
>
> p p > t t~ FCNC=0 DIM6=1, (t > w+ b, w+ > wdecay wdecay), (t~ > w- b~ , w- > wdecay wdecay)
> where wdecay = g u c d s u~ c~ d~ s~ e+ mu+ tau+ e- mu- tau- ve vm vt ve~ vm~ vt~
>
> I am however running into problems. When trying to create the gridpacks i end up with:
>
> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 0 +- 0 pb
> Nb of events : 0
>
> before failing with the error:
>
> Command "generate_events " interrupted with error:
> IndexError : list index out of range
>
> and in the log file (run_01_
>
> File "/scratch/
> return self.onecmd_
> File "/scratch/
> return func(arg, **opt)
> File "/scratch/
> self.exec_
> File "/scratch/
> stop = Cmd.onecmd_
> File "/scratch/
> return func(arg, **opt)
> File "/scratch/
> proc_charac=
> File "/scratch/
> return super(MultiEven
> File "/scratch/
> max_wgt = all_wgt[-1]
> IndexError: list index out of range
>
> Looking into some of the SubProcesses logs, I can see that none of the events passed the phase space cuts and it says to loosen cuts or increase max_events. I do not know where to increase the max_events from and the phase space cuts i think are the default ones which work fine with SM ttbar production.
>
> Cheers,
> Doug
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#3 |
Thanks Olivier Mattelaer, that solved my question.