MadGraph5_aMC@NLO

Soft test Failed in ttbar production (EFT)

Asked by Andres Vasquez

Dear MadGraph team,

I am studying ttbar production at NLO in the SMEFT. For this I want to take advantage of the features of Madgraph5 v3, specifically I am using MG5_aMC_v3_0_1. I am copying below the content of the proc_card_mg5.dat

set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set include_lepton_initiated_processes False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
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~
import model SMEFTtopNLO
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate p p > t t~ NP=1 QCD=2 QED=0 [QCD]
output

The issue comes when I try to generate events, as the result of typing 'launch' is

INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
INFO: Compiling the code
INFO: For gauge cancellation, the width of 't' has been set to zero.
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 8 cores
INFO: Compiling P0_gg_ttx...
INFO: Compiling P0_uux_ttx...
INFO: Compiling P0_ccx_ttx...
INFO: Compiling P0_ddx_ttx...
INFO: Compiling P0_uxu_ttx...
INFO: Compiling P0_cxc_ttx...
INFO: Compiling P0_dxd_ttx...
INFO: Compiling P0_bbx_ttx...
INFO: P0_uxu_ttx done.
INFO: Compiling P0_bxb_ttx...
INFO: P0_bbx_ttx done.
INFO: P0_ccx_ttx done.
INFO: P0_ddx_ttx done.
INFO: P0_dxd_ttx done.
INFO: P0_cxc_ttx done.
INFO: P0_uux_ttx done.
INFO: P0_gg_ttx done.
INFO: P0_bxb_ttx done.
INFO: Checking test output:
INFO: P0_gg_ttx
INFO: Result for test_ME:
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
 SPLIT TYPE USED: F T F
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11839999999999999
  ****************************************

       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= 5 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.34600E+03 -- 0.34600E+03
tau_min 2 1 : 0.34600E+03 -- 0.34600E+03
tau_min 3 1 : 0.34600E+03 0.34600E+03 0.34600E+03
tau_min 4 1 : 0.34600E+03 0.34600E+03 0.34600E+03
tau_min 5 1 : 0.34600E+03 -- 0.34600E+03
tau_min 6 1 : 0.34600E+03 -- 0.34600E+03
tau_min 7 1 : 0.34600E+03 -- 0.34600E+03
tau_min 8 1 : 0.34600E+03 -- 0.34600E+03
tau_min 9 1 : 0.34600E+03 -- 0.34600E+03
tau_min 10 1 : 0.34600E+03 -- 0.34600E+03
tau_min 11 1 : 0.34600E+03 -- 0.34600E+03
tau_min 12 1 : 0.34600E+03 -- 0.34600E+03
tau_min 13 1 : 0.34600E+03 -- 0.34600E+03
tau_min 14 1 : 0.34600E+03 -- 0.34600E+03
tau_min 15 1 : 0.34600E+03 -- 0.34600E+03
tau_min 16 1 : 0.34600E+03 -- 0.34600E+03
tau_min 17 1 : 0.34600E+03 -- 0.34600E+03
tau_min 18 1 : 0.34600E+03 -- 0.34600E+03
tau_min 19 1 : 0.34600E+03 -- 0.34600E+03
tau_min 20 1 : 0.34600E+03 -- 0.34600E+03
 BORN: keeping split order 1
 BORN: keeping split order 2
 BORN: not keeping split order 3
 counterterm S.O 1 NP
 BORN: keeping split order 1
 BORN: not keeping split order 2
 BORN: not keeping split order 3
 counterterm S.O 2 QCD
 BORN: keeping split order 1
 BORN: keeping split order 2
 BORN: not keeping split order 3
 counterterm S.O 3 QED
 BORN: not keeping split order 1
 BORN: not keeping split order 2
 BORN: not keeping split order 3

 REAL 1: keeping split order 1
 REAL 1: keeping split order 2
 REAL 1: not keeping split order 3
 Soft tests done for (Born) config 1
 Failures: 103
     Soft test 1 FAILED. Fraction of failures: 1.03

 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= 5 j= 2
 with PDGs: i= 21 j= 21

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

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

 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= 5 j= 3
 with PDGs: i= 21 j= 6

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

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

 No collinear test for massive j_fks

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 4
 FKS partons are: i= 5 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: 100
     Soft test 4 FAILED. Fraction of failures: 1.00

 No collinear test for massive j_fks

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

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

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

 REAL 2: not keeping split order 1
 REAL 2: keeping split order 2
 REAL 2: keeping split order 3
 REAL 2: not keeping split order 4
 REAL 2: keeping split order 5
 REAL 2: not keeping split order 6
 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: 2
Collinear test 5 PASSED. Fraction of failures: 0.02

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

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

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

 REAL 3: not keeping split order 1
 REAL 3: keeping split order 2
 REAL 3: keeping split order 3
 REAL 3: not keeping split order 4
 REAL 3: keeping split order 5
 REAL 3: not keeping split order 6
 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: 2
Collinear test 6 PASSED. Fraction of failures: 0.02

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

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

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

 REAL 4: not keeping split order 1
 REAL 4: keeping split order 2
 REAL 4: keeping split order 3
 REAL 4: not keeping split order 4
 REAL 4: keeping split order 5
 REAL 4: not keeping split order 6
 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: 2
Collinear test 7 PASSED. Fraction of failures: 0.02

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

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

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

 REAL 5: not keeping split order 1
 REAL 5: keeping split order 2
 REAL 5: keeping split order 3
 REAL 5: not keeping split order 4
 REAL 5: keeping split order 5
 REAL 5: not keeping split order 6
 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: 2
Collinear test 8 PASSED. Fraction of failures: 0.02

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

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

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

 REAL 6: not keeping split order 1
 REAL 6: keeping split order 2
 REAL 6: keeping split order 3
 REAL 6: not keeping split order 4
 REAL 6: keeping split order 5
 REAL 6: not keeping split order 6
 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: 2
Collinear test 9 PASSED. Fraction of failures: 0.02

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

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

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

 REAL 7: not keeping split order 1
 REAL 7: keeping split order 2
 REAL 7: keeping split order 3
 REAL 7: not keeping split order 4
 REAL 7: keeping split order 5
 REAL 7: not keeping split order 6
 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: 0
Collinear test 10 PASSED. Fraction of failures: 0.00

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

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

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

 REAL 8: not keeping split order 1
 REAL 8: keeping split order 2
 REAL 8: keeping split order 3
 REAL 8: not keeping split order 4
 REAL 8: keeping split order 5
 REAL 8: not keeping split order 6
 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: 2
Collinear test 11 PASSED. Fraction of failures: 0.02

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

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

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

 REAL 9: not keeping split order 1
 REAL 9: keeping split order 2
 REAL 9: keeping split order 3
 REAL 9: not keeping split order 4
 REAL 9: keeping split order 5
 REAL 9: not keeping split order 6
 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: 0
Collinear test 12 PASSED. Fraction of failures: 0.00

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 13
 FKS partons are: i= 5 j= 2
 with PDGs: i= -2 j= -2

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

 REAL 10: not keeping split order 1
 REAL 10: keeping split order 2
 REAL 10: keeping split order 3
 REAL 10: not keeping split order 4
 REAL 10: keeping split order 5
 REAL 10: not keeping split order 6
 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: 0
Collinear test 13 PASSED. Fraction of failures: 0.00

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

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

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

 REAL 11: not keeping split order 1
 REAL 11: keeping split order 2
 REAL 11: keeping split order 3
 REAL 11: not keeping split order 4
 REAL 11: keeping split order 5
 REAL 11: not keeping split order 6
 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: 0
Collinear test 14 PASSED. Fraction of failures: 0.00

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 15
 FKS partons are: i= 5 j= 2
 with PDGs: i= -4 j= -4

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

 REAL 12: not keeping split order 1
 REAL 12: keeping split order 2
 REAL 12: keeping split order 3
 REAL 12: not keeping split order 4
 REAL 12: keeping split order 5
 REAL 12: not keeping split order 6
 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: 2
Collinear test 15 PASSED. Fraction of failures: 0.02

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 16
 FKS partons are: i= 5 j= 2
 with PDGs: i= 4 j= 4

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

 REAL 13: not keeping split order 1
 REAL 13: keeping split order 2
 REAL 13: keeping split order 3
 REAL 13: not keeping split order 4
 REAL 13: keeping split order 5
 REAL 13: not keeping split order 6
 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: 0
Collinear test 16 PASSED. Fraction of failures: 0.00

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

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

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

 REAL 14: not keeping split order 1
 REAL 14: keeping split order 2
 REAL 14: keeping split order 3
 REAL 14: not keeping split order 4
 REAL 14: keeping split order 5
 REAL 14: not keeping split order 6
 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 17 PASSED. Fraction of failures: 0.00

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

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

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

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

 REAL 15: not keeping split order 1
 REAL 15: keeping split order 2
 REAL 15: keeping split order 3
 REAL 15: not keeping split order 4
 REAL 15: keeping split order 5
 REAL 15: not keeping split order 6
 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 18 PASSED. Fraction of failures: 0.00

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

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 19
 FKS partons are: i= 5 j= 2
 with PDGs: i= -5 j= -5

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

 REAL 16: not keeping split order 1
 REAL 16: keeping split order 2
 REAL 16: keeping split order 3
 REAL 16: not keeping split order 4
 REAL 16: keeping split order 5
 REAL 16: not keeping split order 6
 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 19 PASSED. Fraction of failures: 0.00

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

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 20
 FKS partons are: i= 5 j= 2
 with PDGs: i= 5 j= 5

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

 REAL 17: not keeping split order 1
 REAL 17: keeping split order 2
 REAL 17: keeping split order 3
 REAL 17: not keeping split order 4
 REAL 17: keeping split order 5
 REAL 17: not keeping split order 6
 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 20 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 20 PASSED. Fraction of failures: 0.00
Error detected in "launch auto "
write debug file /home/vasqueztocor/MG5_aMC_v3_0_1/PROCNLO_SMEFTtopNLO_11/run_02_tag_1_debug.log
If you need help with this issue, please, contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Some tests failed, run cannot continue. Please search on https://answers.launchpad.net/mg5amcnlo for more information, and in case there is none, report the problem there.
quit

This doesn't look like an issue from the UFO model I am using, as all the Wilson coefficients are set to zero in the param_card.dat, besides I have tested the model in MG5 v2 and I could generate events.

I would be thankful if you could shed some light about this issue.

If you require the UFO model I could send it by email.

Andres.

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Rikkert Frederix Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi,
there are at least two sources of problems with the EFT computations. One is rather trivial, and is related to the very error you are showing. Basically it is due to fake non-zero results appearing in some amplitudes, which are not correctly recognised. It has been fixed in 3.0.2, wich you can get via bzr

https://bazaar.launchpad.net/~maddevelopers/mg5amcnlo/3.0.2

(this is the fix https://bazaar.launchpad.net/~maddevelopers/mg5amcnlo/3.0.2/revision/967)

The other problem require few hacks in the code, and prevents to have double EFT insertions in the virtuals (so EFT-like corrections to the pure-QCD born). This is a known fact, which I discussed in details with Eleni Vryonidou. In order to have it working, you have to comment the following lines (line numbers refer to the latest 3.0.2 revision). The changes amount in commenting 4 lines (those commented with ###MZ )

Please let me know if this helps (I checked and everything seems to go fine)

Cheers,

Marco

=== modified file 'madgraph/fks/fks_base.py'
--- madgraph/fks/fks_base.py 2018-11-12 10:58:08 +0000
+++ madgraph/fks/fks_base.py 2019-09-19 11:28:01 +0000
@@ -331,7 +331,7 @@
             myproc = copy.copy(born.born_amp['process'])
             # include all particles in the loops
             # i.e. allow all orders to be perturbed
- myproc['perturbation_couplings'] = myproc['model']['coupling_orders']
+ ###MZ myproc['perturbation_couplings'] = myproc['model']['coupling_orders']
             # take the orders that are actually used bu the matrix element
 #MZ myproc['orders'] = loop_orders
             myproc['legs'] = fks_common.to_legs(copy.copy(myproc['legs']))
@@ -393,7 +393,7 @@
                 self.process['perturbation_couplings'].append(o)
         # set the orders to empty, to force the use of the squared_orders
         self.process['orders'] = copy.copy(born_proc['orders'])
- self.process['orders'] = {}
+ ###MZ self.process['orders'] = {}

         legs = [(leg.get('id'), leg) for leg in leglist]
         self.pdgs = array.array('i',[s[0] for s in legs])

=== modified file 'madgraph/interface/amcatnlo_interface.py'
--- madgraph/interface/amcatnlo_interface.py 2018-12-17 09:00:23 +0000
+++ madgraph/interface/amcatnlo_interface.py 2019-09-19 11:49:27 +0000
@@ -542,10 +542,10 @@
                     ord in myprocdef['perturbation_couplings']])
         # finally set perturbation_couplings to **all** the coupling orders
         # avaliable in the model
- myprocdef['perturbation_couplings'] = list(myprocdef['model']['coupling_orders'])
-
-
- myprocdef['orders'] = {}
+ ###MZ myprocdef['perturbation_couplings'] = list(myprocdef['model']['coupling_orders'])
+
+
+ ###MZ myprocdef['orders'] = {}
         self._curr_proc_defs.append(myprocdef)

 # if myprocdef['perturbation_couplings']!=['QCD']:

Revision history for this message
Andres Vasquez (avasquezt) said : 		#3

Hi Marco, thanks for your reply.

I made all the changes you suggested, then used the same syntax as I wrote in the question above:

 'generate p p > t t~ NP=1 QCD=2 QED=0 [QCD]'

 With the changes this syntax works because it doesn't generate any diagram involving the EFT op.

Then, to get a cross-section sensitive to EFT op. I tried the syntax

 'generate p p > t t~ NP=2 QCD=2 QED=0 [QCD]'

but I couldn't get events, again the same Soft test failures happen when I launch with fixed_order = ON.

Cheers,

Andres.

Revision history for this message
marco zaro (marco-zaro) said : 		#4

Hi Andres,
I think the problem is related to the fact that in the model you have some qqphoton interactions which have just QNP=2 as coupling order, so they cannot be filtered out. I suggest to change the, adding QED=1, so that they will not be included.

Let me know.

Cheers,

Marco

Revision history for this message
Andres Vasquez (avasquezt) said : 		#5

Hi Marco,

Thanks for pointing out the issue with the qqphoton interactions. I got rid of them by calling a restrict model, where I set the coefficients c3pl1, c3pl2 and c3pl3 to zero. With this the real emission diagrams involving the photon are not generated anymore.

Finally I am getting the same diagrams in MG5 v3 and in MG5 v2, however in v3 I am still getting the message

aMCatNLOError : Poles do not cancel, run cannot continue

This happens with the syntax

 'generate p p > t t~ NP=2 QCD=2 QED=0 [QCD]'

Do you have any hint about this?

Thanks in advance,

Andres.

Revision history for this message
marco zaro (marco-zaro) said : 		#6

Can you send me the restrict card?
(mzaro at nikhef dot nl
Thanks,

Marco

> On 23 Sep 2019, at 17:10, Andres Vasquez <email address hidden> wrote:
>
> Question #683957 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/683957
>
> Andres Vasquez posted a new comment:
> Hi Marco,
>
> Thanks for pointing out the issue with the qqphoton interactions. I got
> rid of them by calling a restrict model, where I set the coefficients
> c3pl1, c3pl2 and c3pl3 to zero. With this the real emission diagrams
> involving the photon are not generated anymore.
>
> Finally I am getting the same diagrams in MG5 v3 and in MG5 v2, however
> in v3 I am still getting the message
>
> aMCatNLOError : Poles do not cancel, run cannot continue
>
> This happens with the syntax
>
> 'generate p p > t t~ NP=2 QCD=2 QED=0 [QCD]'
>
> Do you have any hint about this?
>
> Thanks in advance,
>
> Andres.
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
Fernando Cornet (fcornetgomez) said : 		#7

Hi Marco,
does this still apply for version 3.4.2? We are generating p p > t t~ QED=0 QCD=2 NP=2 [QCD] and applying some cuts in the invariant mass.

If we do set mxx_min_pdg {6:640} it fails in the soft test but if we do set mxx_min_pdg {6:690} it works.

Thanks,
Fernando

Revision history for this message
marco zaro (marco-zaro) said : 		#8

Hi Fernando,
how did we solve the problem long ago, or how did the discussion end up?

Marco

Revision history for this message
Reza Goldouzian (rgoldouz) said (last edit ): 		#9

Dear Marco and Fernando,

I am also facing the same issue using SMEFT@NLO with restricted NLO cards (all WC are set to zero except 2-light-2-heavy and CtG).
I am using MG VERSION 3.5.1. I receive similar error when I "set mxx_min_pdg {6:600}".
Please see the attached log file.
Is this issue fixed? How can I fix it?

Thanks a lot,
-Reza

Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt e- mu- ta- ve~ vm~ vt~ e+ mu+ ta+ t t~ z w+ ghz ghwp ghwm h w- ghz~ ghwp~ ghwm~
define p = g u c d s u~ c~ d~ s~ b b~
Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s u~ c~ d~ s~ b b~
Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
set acknowledged_v3.1_syntax True --global
save options acknowledged_v3.1_syntax
save configuration file to /scratch365/rgoldouz/MG33/new/genproductions/bin/MadGraph5_aMCatNLO/tt_NLO_SMEFT_mttMIN600/tt_NLO_SMEFT_mttMIN600_gridpack/work/MG5_aMC_v3_5_1/input/mg5_configuration.txt
generate p p > t t~ NP=2 QCD=2 QED=0 [QCD]@0
^[[1;34mInterpreting 'QCD=2' as 'QCD<=2'^[[0m
^[[1;34mInterpreting 'NP=2' as 'NP<=2'^[[0m
^[[1;34mOrder QED is not constrained as squared_orders. Using: QED^2=0^[[0m
^[[1;34mOrder QCD is not constrained as squared_orders. Using: QCD^2=4^[[0m
^[[1;34mOrder NP is not constrained as squared_orders. Using: NP^2=4^[[0m
INFO: Generating FKS-subtracted matrix elements for born process: g g > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (1 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: u u~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (2 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: c c~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (3 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: d d~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (4 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: s s~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (5 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: b b~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (6 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: u~ u > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (7 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: c~ c > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (8 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: d~ d > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (9 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: s~ s > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (10 / 11)
INFO: Generating FKS-subtracted matrix elements for born process: b~ b > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (11 / 11)
INFO: Generating virtual matrix element with MadLoop for process: g g > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (1 / 11)
^[[1;34mWARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them. ^[[0m
INFO: Generating virtual matrix element with MadLoop for process: u u~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (2 / 11)
^[[1;34mWARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions. ^[[0m
INFO: Generating virtual matrix element with MadLoop for process: c c~ > t t~ NP<=2 QCD<=3 NP^2=4 QCD^2=6 QED^2=0 (3 / 11)
^[[1;34mWARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions. ^[[0m
.
.
.
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 64 cores
INFO: Compiling P0_gg_ttx...
INFO: Compiling P0_uux_ttx...
INFO: Compiling P0_ccx_ttx...
INFO: Compiling P0_ddx_ttx...
INFO: Compiling P0_bbx_ttx...
INFO: Compiling P0_uxu_ttx...
INFO: Compiling P0_cxc_ttx...
INFO: Compiling P0_dxd_ttx...
INFO: Compiling P0_bxb_ttx...
INFO: P0_uxu_ttx done.
INFO: P0_uux_ttx done.
INFO: P0_cxc_ttx done.
INFO: P0_ddx_ttx done.
INFO: P0_ccx_ttx done.
INFO: P0_bxb_ttx done.
INFO: P0_dxd_ttx done.
INFO: P0_bbx_ttx done.
INFO: P0_gg_ttx done.
INFO: Checking test output:
INFO: P0_gg_ttx
INFO: Result for test_ME:
^[[1mINFO: Output of the failing test:
==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.3.0 loading /cvmfs/cms.cern.ch/slc7_amd64_gcc900/external/lhapdf/6.3.0/share/LHAPDF/NNPDF31_nnlo_as_0118_mc_hessian_pdfas/NNPDF31_nnlo_as_0118_mc_hessian_pdfas_0000.dat
NNPDF31_nnlo_as_0118_mc_hessian_pdfas PDF set, member #0, version 1; LHAPDF ID = 325300
 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
 SPLIT TYPE USED: F T F
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11839999999999999
 New value of alpha_s from PDF lhapdf : 0.11800208008122040
 Give FKS configuration number ("0" loops over all)

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 5 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.34500E+03 -- 0.60000E+03
tau_min 2 1 : 0.34500E+03 -- 0.60000E+03
tau_min 3 1 : 0.60000E+03 0.60000E+03 0.60000E+03
tau_min 4 1 : 0.60000E+03 0.60000E+03 0.60000E+03
tau_min 5 1 : 0.34500E+03 -- 0.60000E+03
tau_min 6 1 : 0.34500E+03 -- 0.60000E+03
tau_min 7 1 : 0.34500E+03 -- 0.60000E+03
tau_min 8 1 : 0.34500E+03 -- 0.60000E+03
tau_min 9 1 : 0.34500E+03 -- 0.60000E+03
tau_min 10 1 : 0.34500E+03 -- 0.60000E+03
tau_min 11 1 : 0.34500E+03 -- 0.60000E+03
tau_min 12 1 : 0.34500E+03 -- 0.60000E+03
tau_min 13 1 : 0.34500E+03 -- 0.60000E+03
tau_min 14 1 : 0.34500E+03 -- 0.60000E+03
tau_min 15 1 : 0.34500E+03 -- 0.60000E+03
tau_min 16 1 : 0.34500E+03 -- 0.60000E+03
tau_min 17 1 : 0.34500E+03 -- 0.60000E+03
tau_min 18 1 : 0.34500E+03 -- 0.60000E+03
tau_min 19 1 : 0.34500E+03 -- 0.60000E+03
tau_min 20 1 : 0.34500E+03 -- 0.60000E+03
 BORN: keeping split order 1
 BORN: keeping split order 2
 BORN: keeping split order 3
 counterterm S.O 1 NP
 BORN: keeping split order 1
 BORN: keeping split order 2
 BORN: not keeping split order 3
 counterterm S.O 2 QCD
 BORN: keeping split order 1
 BORN: keeping split order 2
 BORN: keeping split order 3
 counterterm S.O 3 QED
 BORN: not keeping split order 1
 BORN: not keeping split order 2
 BORN: not keeping split order 3

 REAL 1: keeping split order 1
 REAL 1: keeping split order 2
 REAL 1: keeping split order 3
 Soft tests done for (Born) config 1
 Failures: 100 0 0 0 100 100 100 0
Sum of all orders: Soft test 1 FAILED. Fraction of failures: 1.00
Split-order 4: NP: 0; QCD: 6; QED: 0; Soft test 1 FAILED. Fraction of failures: 1.00
Split-order 5: NP: 2; QCD: 6; QED: 0; Soft test 1 FAILED. Fraction of failures: 1.00
Split-order 6: NP: 4; QCD: 6; QED: 0; Soft test 1 FAILED. Fraction of failures: 1.00
Split-order 7: NP: 4; QCD: 2; QED: 0; Soft test 1 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 0 0 0 0 0 0 0 0
Sum of all orders: Collinear test 1 PASSED. Fraction of failures: 0.00
Split-order 4: NP: 0; QCD: 6; QED: 0; Collinear test 1 PASSED. Fraction of failures: 0.00
Split-order 5: NP: 2; QCD: 6; QED: 0; Collinear test 1 PASSED. Fraction of failures: 0.00
Split-order 6: NP: 4; QCD: 6; QED: 0; Collinear test 1 PASSED. Fraction of failures: 0.00
Split-order 7: NP: 4; QCD: 2; QED: 0; Collinear test 1 PASSED. Fraction of failures: 0.00

 =================================================
.
.
.
 Soft tests done for (Born) config 1
 Failures: 0 0 100 100 100 100 100 100
Thanks for using LHAPDF 6.3.0. Please make sure to cite the paper:
  Eur.Phys.J. C75 (2015) 3, 132 (http://arxiv.org/abs/1412.7420)
Sum of all orders: Soft test 19 PASSED. Fraction of failures: 0.00
Split-order 4: NP: 0; QCD: 6; QED: 0; Soft test 19 FAILED. Fraction of failures: 1.00
Split-order 5: NP: 2; QCD: 6; QED: 0; Soft test 19 FAILED. Fraction of failures: 1.00
Split-order 6: NP: 4; QCD: 6; QED: 0; Soft test 19 FAILED. Fraction of failures: 1.00
Split-order 7: NP: 4; QCD: 2; QED: 0; Soft test 19 FAILED. Fraction of failures: 1.00

 Collinear tests done for (Born) config 1
 Failures: 0 0 0 0 0 3 4 0
Sum of all orders: Collinear test 19 PASSED. Fraction of failures: 0.00
Split-order 4: NP: 0; QCD: 6; QED: 0; Collinear test 19 PASSED. Fraction of failures: 0.00
Split-order 5: NP: 2; QCD: 6; QED: 0; Collinear test 19 PASSED. Fraction of failures: 0.03
Split-order 6: NP: 4; QCD: 6; QED: 0; Collinear test 19 PASSED. Fraction of failures: 0.04
Split-order 7: NP: 4; QCD: 2; QED: 0; Collinear test 19 PASSED. Fraction of failures: 0.00

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

 NEW FKS CONFIGURATION:
 FKS configuration number is 20
 FKS partons are: i= 5 j= 2
 with PDGs: i= 5 j= 5

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

 REAL 17: keeping split order 1
 REAL 17: keeping split order 2
 REAL 17: keeping split order 3
 REAL 17: keeping split order 4
 REAL 17: keeping split order 5
 REAL 17: keeping split order 6
 Soft tests done for (Born) config 1
 Failures: 0 0 100 100 100 100 100 100
Sum of all orders: Soft test 20 PASSED. Fraction of failures: 0.00
Split-order 4: NP: 0; QCD: 6; QED: 0; Soft test 20 FAILED. Fraction of failures: 1.00
Split-order 5: NP: 2; QCD: 6; QED: 0; Soft test 20 FAILED. Fraction of failures: 1.00
Split-order 6: NP: 4; QCD: 6; QED: 0; Soft test 20 FAILED. Fraction of failures: 1.00
Split-order 7: NP: 4; QCD: 2; QED: 0; Soft test 20 FAILED. Fraction of failures: 1.00

 Collinear tests done for (Born) config 1
 Failures: 0 0 0 0 0 2 5 0
Sum of all orders: Collinear test 20 PASSED. Fraction of failures: 0.00
Split-order 4: NP: 0; QCD: 6; QED: 0; Collinear test 20 PASSED. Fraction of failures: 0.00
Split-order 5: NP: 2; QCD: 6; QED: 0; Collinear test 20 PASSED. Fraction of failures: 0.02
Split-order 6: NP: 4; QCD: 6; QED: 0; Collinear test 20 PASSED. Fraction of failures: 0.05
Split-order 7: NP: 4; QCD: 2; QED: 0; Collinear test 20 PASSED. Fraction of failures: 0.00 ^[[0m

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#10

Hi Rik,

Is this not the same issue as the one on the pt cut of the photon that you fix with this patch:
https://github.com/mg5amcnlo/mg5amcnlo/commit/a7fcab7a06c53d506a1ee03280d019cf7dd57715

If so can we do similar fix in this case?

Olivier

Revision history for this message
Reza Goldouzian (rgoldouz) said : 		#11

Dear Marco and Fernando,

Do you have any update/news about this issue?

Thanks a lot,
-Reza

Revision history for this message
Rikkert Frederix (frederix) said : 		#12

Hi Olivier,

You are right. It should be enough to add the PDG-level cuts also to the test_soft_col_limits.f. (Note that https://github.com/mg5amcnlo/mg5amcnlo/commit/a7fcab7a06c53d506a1ee03280d019cf7dd57715 is not enough to fix the issue. That whole block should be moved inside the loop over FKS directories, as is done in this fix: https://github.com/mg5amcnlo/mg5amcnlo/commit/ccca2744d753ee0dcc2b2f1d1d20fc03bb80664e )

I'm a bit busy right now -- could you quickly implement this? Otherwise I might have a look in about 3 weeks.

Best,
Rikkert

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#13

Hi,

Actually your fix for 3.5.2 should then fix it already...
(It just that I have to implement the same fix in 2.9.18)

@Reza can you try 3.5.2?

Cheers,

Olivier

Revision history for this message
Rikkert Frederix (frederix) said : 		#14

Hi Olivier,

I agree. In 3.5.2 this should already be fixed.

Best,
Rikkert

Can you help with this problem?

Provide an answer of your own, or ask Andres Vasquez for more information if necessary.

To post a message you must log in.