check_poles error on SMEFT@NLO

Asked by Eirini Kasimi on 2020-06-27

Dear experts.

I am trying to generate the decomposition terms using the SMEFT@NLO model,
I managed to generate all terms except the interference term between SM-BSM.

The syntax is the following:

import model SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO-NLO
generate p p > mu+ mu- e- ve~ QCD=99 QED=99 NP^2==2 [QCD]

and the error I get when I am trying to generate the events is :

Error detected in "launch auto "
write debug file /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/SMEFT_NLO/test/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Poles do not cancel, run cannot continue

The run_01_tag_1_debug.log file is the following:

launch auto
Traceback (most recent call last):
  File "/mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/madgraph/interface/extended_cmd.py", line 1515, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/madgraph/interface/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/madgraph/interface/amcatnlo_run_interface.py", line 1693, in do_launch
    self.compile(mode, options)
  File "/mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/madgraph/interface/amcatnlo_run_interface.py", line 5242, in compile
    self.check_tests(test, this_dir)
  File "/mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/madgraph/interface/amcatnlo_run_interface.py", line 5255, in check_tests
    return self.parse_check_poles_log(pjoin(dir, '%s.log' % test))
  File "/mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/madgraph/interface/amcatnlo_run_interface.py", line 5290, in parse_check_poles_log
    raise aMCatNLOError('Poles do not cancel, run cannot continue')
aMCatNLOError: Poles do not cancel, run cannot continue
aMCatNLOError: Poles do not cancel, run cannot continue
Value of current Options:
              text_editor : None
      notification_center : True
       cluster_local_path : None
  default_unset_couplings : 99
       group_subprocesses : Auto
ignore_six_quark_processes : False
    loop_optimized_output : True
    cluster_status_update : (600, 30)
         fortran_compiler : None
               hepmc_path : None
                  collier : None
              auto_update : 7
             pythia8_path : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/HEPTools/pythia8
                hwpp_path : None
low_mem_multicore_nlo_generation : False
                    golem : None
          pythia-pgs_path : None
                  td_path : None
             delphes_path : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/Delphes
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/HEPTools/madanalysis5/madanalysis5
      exrootanalysis_path : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/ExRootAnalysis
                      OLP : MadLoop
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet : None
                 run_mode : 2
              web_browser : None
   automatic_html_opening : False
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path : None
         madanalysis_path : None
                   lhapdf : lhapdf-config
             stdout_level : 20
                  nb_core : 8
            f2py_compiler : None
                    ninja : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/HEPTools/lib
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
           crash_on_error : False
mg5amc_py8_interface_path : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7/HEPTools/MG5aMC_PY8_interface
         loop_color_flows : False
                  samurai : None
         cluster_nb_retry : 1
                 mg5_path : /mnt/atlas_b/users/e/ekasimi/madgraph/MG5_aMC_v2_6_7
                  timeout : 60
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0

Do you have any idea what causes this error and how can I solve it?

Thanks in advance,
Eirini Kasimi

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Eleni Vryonidou Edit question
Solved by:
Eleni Vryonidou
Solved:
2020-06-29
Last query:
2020-06-29
Last reply:
2020-06-29
Best Eleni Vryonidou (evryonidou) said : #1

Hi Eirini,

the squared syntax doesn't work with [QCD]. So you can't do the decomposition like that at NLO. You need to do a fit by running for different values of the coefficients.

I tried

p p > mu+ mu- e- ve~ QCD=0 QED=4 NP=2 [QCD]

and this has no problems with pole cancellation.

I think it's important you specify the orders "QCD=0 QED=4 NP=2"

Cheers,

Eleni

Eirini Kasimi (eirinikasimi) said : #2

Thanks Eleni Vryonidou, that solved my question.