aMCatNLOError : Poles do not cancel, run cannot continue

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.

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

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

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

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

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.
>