aMCatNLOError : Poles do not cancel, run cannot continue

Asked by Deepak Kumar on 2018-05-02

Dear Olivier,

I am trying to generate this process:

MG5_aMC>import model Pseudoscalar_2HDM
generate p p > b b~ xd xd~ [QCD]

but I am getting error[1]

Can you please let me know why I am getting this error.

Thanks,
Deepak

[1]:
Error detected in "launch auto "
write debug file /Users/dekumar/Documents/MG5_aMC_v2_5_5/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

---------------------------------------
LOGFILE:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.5.5 2017-05-26 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch auto
Traceback (most recent call last):
  File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/extended_cmd.py", line 1430, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/extended_cmd.py", line 1384, in onecmd_orig
    return func(arg, **opt)
  File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 1226, in do_launch
    self.compile(mode, options)
  File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 4712, in compile
    self.check_tests(test, this_dir)
  File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 4725, in check_tests
    return self.parse_check_poles_log(pjoin(dir, '%s.log' % test))
  File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 4760, in parse_check_poles_log
    raise 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
                    pjfry : None
       cluster_local_path : None
       group_subprocesses : Auto
ignore_six_quark_processes : False
    loop_optimized_output : True
    cluster_status_update : (600, 30)
         fortran_compiler : None
               hepmc_path : None
                  collier : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/lib
              auto_update : 7
             pythia8_path : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/pythia8
                hwpp_path : None
low_mem_multicore_nlo_generation : False
                    golem : None
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : None
                      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 : /Users/dekumar/Documents/MG5_aMC_v2_5_5/SysCalc
         madanalysis_path : None
                   lhapdf : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/lhapdf6/bin/lhapdf-config
             stdout_level : 20
                  nb_core : 4
            f2py_compiler : None
                    ninja : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/lib
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
           crash_on_error : False
mg5amc_py8_interface_path : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/MG5aMC_PY8_interface
         loop_color_flows : False
                  samurai : None
         cluster_nb_retry : 1
                 mg5_path : /Users/dekumar/Documents/MG5_aMC_v2_5_5
                  timeout : 60
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
Deepak Kumar
Solved:
2018-05-15
Last query:
2018-05-15
Last reply:
2018-05-02

Hi,

Your model is not in our database:

MG5_aMC>import model Pseudoscalar_2HDM
DEBUG: no model with that name found online
Error detected in "import model Pseudoscalar_2HDM"
write debug file MG5_debug
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
UFOImportError : Path Pseudoscalar_2HDM is not a valid pathname

So i can not reproduce / comment in a meaningful way.

This error means that the code realise that you do not have a proper cancelation of the pole between the real and the virtual contribution. In other word, it detects that something breaks the NLO accuracy of the computation.

The three main reasons for such type of problem are

1) you are using a 5 flavor model (then you can not specify "b" in the final state and be NLO accurate)

2) Your real contribution has some additional resonant diagram. Those should be subtracted to have a NLO accurate result and such subtraction method are not implemented so far. (In some case they are work around but not in all case)

3) You have some not pure QCD loop (loop with W/Z propagating in the loop) that are not included so far.
Such loop can makes the poles cancellation to fail. Sometimes it is fine to ignore that when you know that such loop are irrelevant.

Cheers,

Olivier

> On 2 May 2018, at 14:13, Deepak Kumar <email address hidden> wrote:
>
> New question #668490 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/668490
>
> Dear Olivier,
>
> I am trying to generate this process:
>
> MG5_aMC>import model Pseudoscalar_2HDM
> generate p p > b b~ xd xd~ [QCD]
>
> but I am getting error[1]
>
> Can you please let me know why I am getting this error.
>
> Thanks,
> Deepak
>
> [1]:
> Error detected in "launch auto "
> write debug file /Users/dekumar/Documents/MG5_aMC_v2_5_5/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
>
>
>
> ---------------------------------------
> LOGFILE:
>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.5.5 2017-05-26 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* and *
> #* http://amcatnlo.cern.ch *
> #* *
> #************************************************************
> #* *
> #* Command File for aMCatNLO *
> #* *
> #* run as ./bin/aMCatNLO.py filename *
> #* *
> #************************************************************
> launch auto
> Traceback (most recent call last):
> File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/extended_cmd.py", line 1430, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/extended_cmd.py", line 1384, in onecmd_orig
> return func(arg, **opt)
> File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 1226, in do_launch
> self.compile(mode, options)
> File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 4712, in compile
> self.check_tests(test, this_dir)
> File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 4725, in check_tests
> return self.parse_check_poles_log(pjoin(dir, '%s.log' % test))
> File "/Users/dekumar/Documents/MG5_aMC_v2_5_5/madgraph/interface/amcatnlo_run_interface.py", line 4760, in parse_check_poles_log
> raise 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
> pjfry : None
> cluster_local_path : None
> group_subprocesses : Auto
> ignore_six_quark_processes : False
> loop_optimized_output : True
> cluster_status_update : (600, 30)
> fortran_compiler : None
> hepmc_path : None
> collier : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/lib
> auto_update : 7
> pythia8_path : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/pythia8
> hwpp_path : None
> low_mem_multicore_nlo_generation : False
> golem : None
> pythia-pgs_path : None
> td_path : None
> delphes_path : None
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : None
> 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 : /Users/dekumar/Documents/MG5_aMC_v2_5_5/SysCalc
> madanalysis_path : None
> lhapdf : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/lhapdf6/bin/lhapdf-config
> stdout_level : 20
> nb_core : 4
> f2py_compiler : None
> ninja : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/lib
> amcfast : amcfast-config
> cluster_retry_wait : 300
> output_dependencies : external
> crash_on_error : False
> mg5amc_py8_interface_path : /Users/dekumar/Documents/MG5_aMC_v2_5_5/HEPTools/MG5aMC_PY8_interface
> loop_color_flows : False
> samurai : None
> cluster_nb_retry : 1
> mg5_path : /Users/dekumar/Documents/MG5_aMC_v2_5_5
> timeout : 60
> gauge : unitary
> complex_mass_scheme : False
> cpp_compiler : None
> max_npoint_for_channel : 0
>
>
>
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Deepak Kumar (deepak1kumar) said : #2

Dear Olivier,

This is the model which I am using :
https://gitlab.cern.ch/lhc-dmwg-material/model-repository/tree/master/models/Pseudoscalar_2HDM

Thanks,
Deepak

Hi,

Actually it does work on my laptop.

But if I look at the debug statement, you can see that the pole cancelation is not 100% perfect. So the fact that it crash on some machine and not some other is possible.

I will assign this to Marco who should be able to comment more on this:

Cheers,

Olivier

PS for marco:

INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxu_xdxdxbbx
INFO: Result for test_ME:
INFO: Passed.
DEBUG:
     Soft test 1 PASSED. Fraction of failures: 0.00
Collinear test 1 PASSED. Fraction of failures: 0.00
     Soft test 2 PASSED. Fraction of failures: 0.01
Collinear test 2 PASSED. Fraction of failures: 0.00
     Soft test 3 PASSED. Fraction of failures: 0.00
     Soft test 4 PASSED. Fraction of failures: 0.00
     Soft test 5 PASSED. Fraction of failures: 0.00
Collinear test 5 PASSED. Fraction of failures: 0.28
     Soft test 6 PASSED. Fraction of failures: 0.00
Collinear test 6 PASSED. Fraction of failures: 0.21
INFO: Result for test_MC:
INFO: Passed.
DEBUG:
     Soft test 1 PASSED. Fraction of failures: 0.00
Collinear test 1 PASSED. Fraction of failures: 0.00
     Soft test 2 PASSED. Fraction of failures: 0.00
Collinear test 2 PASSED. Fraction of failures: 0.00
     Soft test 3 PASSED. Fraction of failures: 0.00
     Soft test 4 PASSED. Fraction of failures: 0.00
     Soft test 5 PASSED. Fraction of failures: 0.00
Collinear test 5 PASSED. Fraction of failures: 0.00
     Soft test 6 PASSED. Fraction of failures: 0.00
Collinear test 6 PASSED. Fraction of failures: 0.00

marco zaro (marco-zaro) said : #4

Hi,
apologies for the delay in the answer.
Is the problem solved?
@Olivier, what you see as 'not perfect' (with ~20% failures) are not the tests for the poles cancelation, rather those for the correct numerical behaviour of the matrix elements close to the singular (soft/collinear) limit. For a 2->4 process (with colored particles in the final state?) similar fraction of failures are expected.

Best,

Marco

Deepak Kumar (deepak1kumar) said : #6

Hi Marco,

My problem has been solved I was not using correct model.

But now I want to set mediator pT cut in 2HDM+a model so I took help from here :
https://github.com/violatingcp/MadGraph5_aMCatNLO_Grid/blob/master/Cards/BBbarDMJets_pseudo_NLO_Mchi_Mphi_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/cuts.f#L76

and added it in my cut.f file, but it was not working.

Now the problem is that I am not able to generate any process.
If I am generating new process then getting this error[1]

[1]: Error detected in "launch auto "
write debug file /Users/dekumar/DATA/Software/MG5_aMC_v2_6_2/test/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : CutTools installation in /Users/dekumar/DATA/Software/MG5_aMC_v2_6_2/vendor/CutTools/includects/libcts.a seems to have been compiled with a different compiler than the one specified in MG5_aMC. Please recompile CutTools.

Please help

Thanks,
Deepak

marco zaro (marco-zaro) said : #7

Hi Deepak,
did you change anything related to the compiler you used?
 /Users/dekumar/DATA/Software/MG5_aMC_v2_6_2/vendor/CutTools/includects/libcts.a seems to have been compiled with a different compiler than the one specified in MG5_aMC. Please recompile CutTools.
Otherwise, you need to run make clean inside vendor/CutTools and vendor/IREGI/src

Let me know

Marco

Deepak Kumar (deepak1kumar) said : #8

Hi Marco,

I am using iMac desktop. In my iMac pyroot was not working so i uninstalled
it and install again with command.

brew install root --with-python@2 --without-pythonearlier I installed
root with the command brew install root

Now I have uninstalled root again.

I did make clean inside CutTools and IREGI/src. But still getting error[1].

Thanks,

Deepak

[1]:

Working on SubProcesses

INFO: P1_qq_ll

*Error detected in "generate_events run_01"*

*write debug file
/Users/dekumar/DATA/Software/MG5_aMC_v2_6_2/test/run_01_tag_1_debug.log
*

*If you need help with this issue please contact us on
https://answers.launchpad.net/mg5amcnlo
<https://answers.launchpad.net/mg5amcnlo>*

*MadGraph5Error : A compilation Error occurs when trying to compile
/Users/dekumar/DATA/Software/MG5_aMC_v2_6_2/test/SubProcesses/P1_qq_ll.*

* The compilation fails with the following output message:*

* gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w
-c driver.f -I../../Source/*

* gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w
-c myamp.f -I../../Source/*

* gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w
-c genps.f -I../../Source/*

* gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w
-c unwgt.f -I../../Source/*

* genps.f:328:25:*

* & nbranch,P,M,S,X,jac,pswgt)*

* 1*

* Error: Actual argument contains too few elements for dummy
argument 'x' (16/40) at (1)*

* genps.f:60:53:*

* call gen_mom(iconfig,mincfig,maxcfig,invar,wgt,x,p)*

* 1*

* Error: Actual argument contains too few elements for dummy
argument 'x' (5/16) at (1)*

* make: *** [genps.o] Error 1*

* make: *** Waiting for unfinished jobs....*

* Please try to fix this compilations issue and retry.*

* Help might be found at https://answers.launchpad.net/mg5amcnlo
<https://answers.launchpad.net/mg5amcnlo>.*

* If you think that this is a bug, you can report this at
https://bugs.launchpad.net/mg5amcnlo
<https://bugs.launchpad.net/mg5amcnlo>*

quit

On Wed, May 16, 2018 at 12:23 PM, marco zaro <
<email address hidden>> wrote:

> Your question #668490 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/668490
>
> marco zaro posted a new comment:
> Hi Deepak,
> did you change anything related to the compiler you used?
> /Users/dekumar/DATA/Software/MG5_aMC_v2_6_2/vendor/CutTools/includects/libcts.a
> seems to have been compiled with a different compiler than the one
> specified in MG5_aMC. Please recompile CutTools.
> Otherwise, you need to run make clean inside vendor/CutTools and
> vendor/IREGI/src
>
> Let me know
>
> Marco
>
> --
> You received this question notification because you asked the question.
>