SMEFT NLO ttbar - Inconsistent integrals

Dear Pim,
we are aware of this bug, we are fixing it in a development branch. As soon as possible we will let you know how to correct the code behaviour.

Cheers,

Marco

Dear Marco,
Thats great to hear. Hope to hear from you soon then.
Cheers
Pim

Dear Marco,
Just checking if there has been any progress on the matter. Any indication on when this might be fixed?
Kind regards,
Pim

Dear Pim,
apologies for not having updated you. You can try this branch
https://code.launchpad.net/~maddevelopers/mg5amcnlo/3.0.3-betternlops
which includes the new sintax from 3.0.3-neworders, and fixes the behaviour at NLOPS.

Just a question, how do you shower the events? with the built-in analyses/tools or within a separate framework?

Let me know

Best wishes,

Marco

Hi Marco,
Thanks! I will try the new branch then. I will start with the built-in tools and then look later on to other shower frameworks. Will this branch be merged with an official release anytime soon?
Cheers,
Pim

Hi Pim,
the plans are to merge this branch in the next version (3.0.4).
If you start with the built-in tools (and I suppose you are going to use pythia8), there are few extra changes needed (I did not implement them yet, because I would like to have them implemented in a more flexible way):
Basically, you need to change the length of some string from 50 to 80. Below are the changes to be done.

Best,

Marco

=======

Template/MCatNLO/include/LHEFRead.h

< void lhef_read_wgtsinfo_(int &cwgtinfo_nn, char (cwgtinfo_weights_info[1024][50])) {
---
> void lhef_read_wgtsinfo_(int &cwgtinfo_nn, char (cwgtinfo_weights_info[1024][80])) {
 include/LHEFRead.h:
< void lhef_read_wgtsinfo_(int &cwgtinfo_nn, char (cwgtinfo_weights_info[1024][50])) {
—--
> void lhef_read_wgtsinfo_(int &cwgtinfo_nn, char (cwgtinfo_weights_info[1024][80])) {

< sprintf(cwgtinfo_weights_info[0], "%50s","central value");
--—
> sprintf(cwgtinfo_weights_info[0], “%80s","central value");

<
sprintf(cwgtinfo_weights_info[cwgtinfo_nn],"%50s",sRWGT.c_str());
--—
sprintf(cwgtinfo_weights_info[cwgtinfo_nn],”%80s",sRWGT.c_str());

- If you use the HwU analyses we ship within MG5_aMC:
diff ../Template/NLO/FixedOrderAnalysis/HwU.f FixedOrderAnalysis/
22c22
< character(len=50), allocatable :: wgts_info(:)
---
> character(len=80), allocatable :: wgts_info(:)

 diff ../Template/NLO/MCatNLO/PY8Analyzer/py8an_HwU_rates.f MCatNLO/PY8Analyzer/ ### this has to be done in the actual analysis you use
44c44
< character*50 weights_info(max_weight_shower)
---
> character*80 weights_info(max_weight_shower)

diff -r ../Template/NLO/MCatNLO/srcPythia8/Pythia8.cc <http://pythia8.cc/> MCatNLO/srcPythia8/Pythia8.cc <http://pythia8.cc/>
24c24
< void pyabeg_(int&,char(*)[50]);
---
> void pyabeg_(int&,char(*)[80]);
33c33
< char cwgtinfo_weights_info[1024][50];
---
> char cwgtinfo_weights_info[1024][80];
diff -r ../Template/NLO/MCatNLO/srcPythia8/Pythia82.cc <http://pythia82.cc/> MCatNLO/srcPythia8/Pythia82.cc <http://pythia82.cc/>
25c25
< void pyabeg_(int&,char(*)[50]);
---
> void pyabeg_(int&,char(*)[80]);
34c34
< char cwgtinfo_weights_info[1024][50];
---
> char cwgtinfo_weights_info[1024][80];

> On 8 Oct 2020, at 10:55, Pim Verschuuren <email address hidden> wrote:
>
> Question #692808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692808
>
> Pim Verschuuren posted a new comment:
> Hi Marco,
> Thanks! I will try the new branch then. I will start with the built-in tools and then look later on to other shower frameworks. Will this branch be merged with an official release anytime soon?
> Cheers,
> Pim
>
> --
> You received this question notification because you are assigned to this
> question.

Hi Marco,
The NLO+PS generation was working but with the latest commit 979 to branch betternlops it seems to break down again in the following way:

MG5_aMC>import model SMEFTatNLO-NLO
INFO: Restrict model SMEFTatNLO with file ../3.0.3-betternlops/models/SMEFTatNLO/restrict_NLO.dat .
INFO: Change particles name to pass to MG5 convention
Pass the definition of 'j' and 'p' to 5 flavour scheme.
Kept definitions of multiparticles l- / vl / l+ / vl~ unchanged
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~
MG5_aMC>generate p p > t t~ QCD=2 QED=0 NP=2 [QCD]
Command "generate p p > t t~ QCD=2 QED=0 NP=2 [QCD]" interrupted with error:
TypeError : cannot concatenate 'str' and 'list' objects
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

The problem seems to be the [QCD] argument which when removed removes the error message.

Cheers
Pim

Hi Pim,
I cannot reproduce your error, I can generate the process you tried and
output it.
Note that changes in rev 979 do not affect code generation, only the nlo+PS
analyses..
In case, please include the content of MG5_debug

Cheers,

Marco

On Tue, 13 Oct 2020, 18:45 Pim Verschuuren, <
<email address hidden>> wrote:

> Question #692808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692808
>
> Pim Verschuuren posted a new comment:
> Hi Marco,
> The NLO+PS generation was working but with the latest commit 979 to branch
> betternlops it seems to break down again in the following way:
>
> MG5_aMC>import model SMEFTatNLO-NLO
> INFO: Restrict model SMEFTatNLO with file
> ../3.0.3-betternlops/models/SMEFTatNLO/restrict_NLO.dat .
> INFO: Change particles name to pass to MG5 convention
> Pass the definition of 'j' and 'p' to 5 flavour scheme.
> Kept definitions of multiparticles l- / vl / l+ / vl~ unchanged
> 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~
> MG5_aMC>generate p p > t t~ QCD=2 QED=0 NP=2 [QCD]
> Command "generate p p > t t~ QCD=2 QED=0 NP=2 [QCD]" interrupted with
> error:
> TypeError : cannot concatenate 'str' and 'list' objects
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
>
> The problem seems to be the [QCD] argument which when removed removes the
> error message.
>
> Cheers
> Pim
>
> --
> You received this question notification because you are assigned to this
> question.
>

Hi Marco,
Strangely enough it seems to occur when I run the above on Scientific Linux CERN SLC release 6.10. When I run on CentOs Linux 7 it works fine.
Cheers,
Pim

Hi Pim,
what are the python versions in the two OSs?
Thanks,

Marco

> On 13 Oct 2020, at 21:20, Pim Verschuuren <email address hidden> wrote:
>
> Question #692808 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692808
>
> Pim Verschuuren posted a new comment:
> Hi Marco,
> Strangely enough it seems to occur when I run the above on Scientific Linux CERN SLC release 6.10. When I run on CentOs Linux 7 it works fine.
> Cheers,
> Pim
>
> --
> You received this question notification because you are assigned to this
> question.

Hi Marco,
Thanks for the hint. I changed from Python 2.7.4 GCC 48 to 2.7.10 GCC 62 and the code runs now.
Cheers,
Pim

Thanks marco zaro, that solved my question.