mg@NLO with SMEFTatNLO

HI,

This is a physics limitation of the Long Term Stable version. You need to use the latest development version for such type of computation (so 3.3.2).

Cheers,

Olivier

> On 1 May 2022, at 06:20, Rameswar Sahu <email address hidden> wrote:
>
> New question #701608 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/701608
>
> Hi,
> I am using the model file SMEFTatNLO for a project work. But I am facing some issues. Here is what is happening: I have only kept the operator O_tz non-zero and have generated the process p p > t t~. The new physics contribution starts at QED=2. Still, when I am generating the above process with order QED=2 QCD=2 NP=2 [QCD], I am getting the warning "WARNING: All Born diagrams do not factorize the same sum of power(s) of 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." and the lunch is getting terminated. But when I am generating the process with order QED=2 QCD=0 NP=2 [QCD], everything works fine. But as you can clearly see, in the second method I am losing many standard model diagrams (QED=1 QCD=2 NP=0), which will be present in the background events giving weird results. Can you please suggest the proper way to generate that process?
> Note:
> I tried generating that process with the following set of commands:
> generate p p > t t~ QED=2 QCD=0 NP=2 [QCD]
> add process p p > t t~ QED=1 QCD=2 NP=0 [QCD]
>
> But, I am concerned that may be the interference term between the two sets of processes are not getting included. What are you thoughts.
>
> Best,
> Rameswar
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

actually I am using "MG5_aMC_v3_3_2" and the version "SMEFTatNLO_v1.0.3" of the model. Still I am getting the above mentioned problem. Any suggestions ....

Here are the contents of the log file :
"launch auto
Traceback (most recent call last):
  File "/home/rameswar.s/Desktop/softwares/MG5_aMC_v3_3_2/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/rameswar.s/Desktop/softwares/MG5_aMC_v3_3_2/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/rameswar.s/Desktop/softwares/MG5_aMC_v3_3_2/madgraph/interface/amcatnlo_run_interface.py", line 1783, in do_launch
    self.compile(mode, options)
  File "/home/rameswar.s/Desktop/softwares/MG5_aMC_v3_3_2/madgraph/interface/amcatnlo_run_interface.py", line 5399, in compile
    self.check_tests(test, this_dir)
  File "/home/rameswar.s/Desktop/softwares/MG5_aMC_v3_3_2/madgraph/interface/amcatnlo_run_interface.py", line 5410, in check_tests
    return self.parse_test_mx_log(pjoin(dir, '%s.log' % test))
  File "/home/rameswar.s/Desktop/softwares/MG5_aMC_v3_3_2/madgraph/interface/amcatnlo_run_interface.py", line 5420, in parse_test_mx_log
    raise 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.')
madgraph.interface.amcatnlo_run_interface.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."

Best,
Rameswar

Hi,
Forgot to mention that, I tried running the same process with model "loop_sm-no_b-mass" :
generate p p > t t~ QED=3 QCD=2 [QCD]

I am getting the same error i.e. "WARNING: All Born diagrams do not factorize the same sum of power(s) of 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."

Thanks,
Rameswar

Then here i would suggest to contact smeftatnlo author for guidance.

Cheers,

Olivier

Hi Olivier,

I don't think the problem is with the model file because, as I have mentioned in the last comment, I am getting the same error with the model file "loop_sm." I think the problem arises because of the presence of different QCD orders in a given subprocess (i.e., the appearance of both QED=2 QCD=2 diagrams with QED=0 QCD=4 diagrams together for a subprocess). I am not sure tho. Can you suggest any alternative methods to generate such processes?

Best,
Rameswar

Thanks Olivier Mattelaer, that solved my question.

Thanks Olivier,

I will try to contact the authors by email. Thanks for your suggestions.

Best,

Rameswar