dilog 'silent' error at NLO

Asked by Hesham El Faham

Hi,

I am running one channel of three tops at NLO using the following command lines:
--> import model loop_qcd_qed_sm_Gmu-no_widths
--> generate p p > t t~ t~ t~ w+ [QCD QED]
I do fixed order computation aiming to look at the different contributions from each order. I treat the overlap from tttt, t>wb and the process converges and produces the expected LO1 and NLO1 cross-sections, but then beyond those orders, the other contributions are either zero or just random.

I noticed a long list of dilog errors in the log file of the run, one line of which looks like this:
------- iteration 3
 Update # PS points (even_rn): 29120 --> 24576
 ERROR in OneLOop dilog_r: |x|,iph = .1000000000000000000E+01, 2 , returning 0

The code never complained about this in run time, and so I am claiming the run failed silently and this relates to the random contributions I get beyond LO1 and NLO1.
I am assuming it is an error in the tensor reduction, but I am not sure how to fix or even debug this in mg5, may you please help?

Best,
Hesham

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
davide.pagani.85 Edit question
Solved by:
davide.pagani.85
Solved:
Last query:
Last reply:
Revision history for this message
Hesham El Faham (helfaham) said :
#1

Hi,

I am just adding more information that might be relevant: I tried changing the tensor reduction method in the MadLoop param card. Collier causes issues with poles cancellation. Ninja and CutTools work fine, but then they end up having the above-mentioned problem

Best,
Hesham

Revision history for this message
Best davide.pagani.85 (davide-pagani) said :
#2

Hi Hesham,

as you may already know the problem was due to the parameter OSthres, which was set equal to 1d-13 in Cards/MadLooopParams.dat for version 3.5.1.

From version 3.5.2 will be reverted to the value 1d-8, as was in all the versions before 3.5.1.

Cheers
Davide

Revision history for this message
Hesham El Faham (helfaham) said :
#3

Thanks davide.pagani.85, that solved my question.