Problem for generating generate p p > t t~ w j QED=3 QCD=1 [QCD]

Asked by Meng-Ju Tsai

Hello Madgraph experts,

We are trying to produce ttW EW sample with 1j@NLO (with QCD NLO correction).

Here are how we did to generate the samples. The command works to produce the diagrams. However, when we did the `launch` step, we have issues failing several tests. I added the model "import model loop_sm-no_b_mass" but this doesn't work for me.

Do you have ideas on how I should do with the errors?

Thank you very much!

Cheers,
Meng-Ju

---------- Command we used ----------

```
set gauge unitary

import model loop_sm-no_b_mass
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define w = w+ w-
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define uc = u c
define uc~ = u~ c~
define ds = d s
define ds~ = d~ s~

generate p p > t t~ w j QED=3 QCD=1 [QCD]

output -f

launch
```

---------- Logs ----------

```
INFO: Output of the failing test:
 Enter 0 to compute MC/MC(limit)
       1 to compute MC/ME(limit)
       2 to compute ME/ME(limit)
 Enter xi_i, y_ij to be used in coll/soft tests
  Enter -2 to generate them randomly
 Enter number of tests for soft and collinear limits
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11799999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model. Hence it will be set to: 5
 Give FKS configuration number ("0" loops over all)

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 7 j= 1
 with PDGs: i= 21 j= 21

 Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 1 , 1 , 0
  with seed 32
 Ranmar initialization seeds 10730 9407
tau_min 1 1 : 0.42642E+03 -- 0.43642E+03
tau_min 2 1 : 0.42642E+03 -- 0.43642E+03
tau_min 3 1 : 0.43642E+03 0.43642E+03 0.43642E+03
tau_min 4 1 : 0.43642E+03 0.43642E+03 0.43642E+03
tau_min 5 1 : 0.43642E+03 0.43642E+03 0.43642E+03
tau_min 6 1 : 0.42642E+03 -- 0.43642E+03
tau_min 7 1 : 0.42642E+03 -- 0.43642E+03
tau_min 8 1 : 0.42642E+03 -- 0.43642E+03
tau_min 9 1 : 0.42642E+03 -- 0.43642E+03
tau_min 10 1 : 0.42642E+03 -- 0.43642E+03
tau_min 11 1 : 0.42642E+03 -- 0.43642E+03
tau_min 12 1 : 0.42642E+03 -- 0.43642E+03
tau_min 13 1 : 0.42642E+03 -- 0.43642E+03
tau_min 14 1 : 0.42642E+03 -- 0.43642E+03
tau_min 15 1 : 0.42642E+03 -- 0.43642E+03
tau_min 16 1 : 0.42642E+03 -- 0.43642E+03

 Soft tests done for (Born) config 1
 Failures: 2
     Soft test 1 PASSED. Fraction of failures: 0.02

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 1 PASSED. Fraction of failures: 0.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 2
 FKS partons are: i= 7 j= 2
 with PDGs: i= 21 j= 2

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 3
     Soft test 2 PASSED. Fraction of failures: 0.03

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 2 PASSED. Fraction of failures: 0.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 3
 FKS partons are: i= 7 j= 4
 with PDGs: i= 21 j= 6

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 3 PASSED. Fraction of failures: 0.00

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 4
 FKS partons are: i= 7 j= 5
 with PDGs: i= 21 j= -6

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 4 PASSED. Fraction of failures: 0.00

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 5
 FKS partons are: i= 7 j= 6
 with PDGs: i= 21 j= 1

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 5 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 1
Collinear test 5 PASSED. Fraction of failures: 0.01

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 6
 FKS partons are: i= 7 j= 1
 with PDGs: i= -1 j= -1

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 6 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 96
Collinear test 6 FAILED. Fraction of failures: 0.96

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 7
 FKS partons are: i= 7 j= 1
 with PDGs: i= 1 j= 1

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 7 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 96
Collinear test 7 FAILED. Fraction of failures: 0.96

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 8
 FKS partons are: i= 7 j= 1
 with PDGs: i= -2 j= -2

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 8 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 100
Collinear test 8 FAILED. Fraction of failures: 1.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 9
 FKS partons are: i= 7 j= 1
 with PDGs: i= 2 j= 2

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 9 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 100
Collinear test 9 FAILED. Fraction of failures: 1.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 10
 FKS partons are: i= 7 j= 1
 with PDGs: i= -3 j= -3

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 10 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 97
Collinear test 10 FAILED. Fraction of failures: 0.97

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 11
 FKS partons are: i= 7 j= 1
 with PDGs: i= 3 j= 3

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 11 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 95
Collinear test 11 FAILED. Fraction of failures: 0.95

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 12
 FKS partons are: i= 7 j= 1
 with PDGs: i= -4 j= -4

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 12 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 100
Collinear test 12 FAILED. Fraction of failures: 1.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 13
 FKS partons are: i= 7 j= 1
 with PDGs: i= 4 j= 4

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 13 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 100
Collinear test 13 FAILED. Fraction of failures: 1.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 14
 FKS partons are: i= 7 j= 1
 with PDGs: i= -5 j= -5

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 14 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 96
Collinear test 14 FAILED. Fraction of failures: 0.96

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 15
 FKS partons are: i= 7 j= 1
 with PDGs: i= 5 j= 5

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 15 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 96
Collinear test 15 FAILED. Fraction of failures: 0.96

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 16
 FKS partons are: i= 7 j= 2
 with PDGs: i= -2 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 16 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 2
Collinear test 16 PASSED. Fraction of failures: 0.02
Error detected in "launch auto "
write debug file /atlas/data19/metsai/simulation/MG5_aMC_v2_9_5/ttWNLO/run_04_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str : Some tests failed, run cannot continue.
 Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.

```

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
davide.pagani.85 Edit question
Last query:
Last reply:
Revision history for this message
davide.pagani.85 (davide-pagani) said :
#1

Hi Meng-Ju,

sorry for the late answer.

Using the five-flavor scheme, it is necessary that proton and jets include bottom (anti)quarks.
Also, since the computation includes what is called the NLO3 for this process, 'p' and 'j' must include also photons and leptons for consistency (as discussed in appendix D of our EW automation paper (table 4)).
This is not done in your syntax.
Moreover, the model you are loading is not appropriate when EW calculations are involved. Although you calculate NLO QCD corrections, this are NLO QCD to LO3 therefore NLO3, which involves also NLO EW to LO2. The model should be

loop_qcd_qed_sm_Gmu

However, even with all of that correctly taken into account, I see another fundamental problem. Diagrams in some cases include in the final state both a stable top and a top decayed, ub->t tbar d t (t->W b). This is inconsistent from a computational point of view. This is not a problem of code, the calculation is simply not feasible.

You may try to perform the computation in the four-flavor scheme, where these problems are not present since p and j do not include the bottom.

You can load the model loop_qcd_qed_sm_Gmu_4FS-with_b_mass

I hope this can help

Cheers
Davide

Can you help with this problem?

Provide an answer of your own, or ask Meng-Ju Tsai for more information if necessary.

To post a message you must log in.