aMCatNLOError : Some tests failed, run cannot continue.

Asked by Liqing Zhang on 2018-09-25

Dear aMC@NLO experts,

I'm willing to calculate NLO EWK correction about pp > w+ w+ w- process.
The LO process of pp > w+ w+ w- has been done successfully. But NLO process pp > w+ w+ w- [QED] has a error at the event generation step (after typing "launch")

INFO: load configuration from /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/bin/TEST_WWW/Cards/amcatnlo_configuration.txt
INFO: load configuration from /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/input/mg5_configuration.txt
INFO: load configuration from /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/bin/TEST_WWW/Cards/amcatnlo_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
No valid web browser found. Please set in ./input/mg5_configuration.txt
launch auto
The following switches determine which programs are run:
/================== Description ==================|=========== values ===========|================ other options ================\
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = ON | No NLO+PS available for EW correction |
\================================================================================================================================/
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'madspin=onshell' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>
use 0
INFO: will run in mode: NLO
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. FO_analyse : FO_analyse_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, FO_analyse, enter path][90s to answer]
>
use 0
INFO: Starting run
Not able to open file /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/bin/TEST_WWW/crossx.html since no program configured.Please set one in ./input/mg5_configuration.txt
INFO: Compiling the code
INFO: For gauge cancellation, the width of 'W+' has been set to zero.
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 12 cores
INFO: Compiling P0_udx_wpwpwm...
INFO: Compiling P0_dxu_wpwpwm...
INFO: P0_udx_wpwpwm done.
INFO: P0_dxu_wpwpwm done.
INFO: Checking test output:
INFO: P0_udx_wpwpwm
INFO: Result for test_ME:
INFO: Output of the failing test:
 Enter 0 to compute MC/MC(limit)
       1 to compute MC/ME(limit)
       2 to compute ME/ME(limit)
 Enter xi_i, y_ij to be used in coll/soft tests
  Enter -2 to generate them randomly
 Enter number of tests for soft and collinear limits
 SPLIT TYPE USED: F T
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11799999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model. Hence it will be set to: 5
 Give FKS configuration number ("0" loops over all)

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 6 j= 1
 with PDGs: i= 22 j= 2

 Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 1 , 1 , 0
  with seed 32
 Ranmar initialization seeds 10730 9407
tau_min 1 1 : 0.24126E+03 -- 0.24126E+03
tau_min 2 1 : 0.24126E+03 -- 0.24126E+03
tau_min 3 1 : 0.24126E+03 0.24126E+03 0.24126E+03
tau_min 4 1 : 0.24126E+03 0.24126E+03 0.24126E+03
tau_min 5 1 : 0.24126E+03 -- 0.24126E+03
tau_min 6 1 : 0.24126E+03 -- 0.24126E+03
 BORN: keeping split order 1
 counterterm S.O 1 QCD
 BORN: not keeping split order 1
 counterterm S.O 2 QED
 BORN: keeping split order 1

 REAL 1: keeping split order 1
 Soft tests done for (Born) config 1
 Failures: 155
     Soft test 1 FAILED. Fraction of failures: 1.55

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 1 PASSED. Fraction of failures: 0.00

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 2
 FKS partons are: i= 6 j= 2
 with PDGs: i= 22 j= -1

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 163
     Soft test 2 FAILED. Fraction of failures: 1.63

 Collinear tests done for (Born) config 1
 Failures: 2
Collinear test 2 PASSED. Fraction of failures: 0.02

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 3
 FKS partons are: i= 6 j= 4
 with PDGs: i= 22 j= 24

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 167
     Soft test 3 FAILED. Fraction of failures: 1.67

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 4
 FKS partons are: i= 6 j= 5
 with PDGs: i= 22 j= -24

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 171
     Soft test 4 FAILED. Fraction of failures: 1.71

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 5
 FKS partons are: i= 6 j= 1
 with PDGs: i= -2 j= 22

 Enter graph number (iconfig), '0' loops over all graphs

 REAL 2: keeping split order 1
 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 5 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 2
Collinear test 5 PASSED. Fraction of failures: 0.02

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 6
 FKS partons are: i= 6 j= 2
 with PDGs: i= 1 j= 22

 Enter graph number (iconfig), '0' loops over all graphs

 REAL 3: keeping split order 1
 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 6 PASSED. Fraction of failures: 0.00

 Collinear tests done for (Born) config 1
 Failures: 0
Collinear test 6 PASSED. Fraction of failures: 0.00
Error detected in "launch auto "
write debug file /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/bin/TEST_WWW/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Some tests failed, run cannot continue.
 Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.

The contents of the run_01_tag_1_debug.log file are:

launch auto
Traceback (most recent call last):
  File "/lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/madgraph/interface/extended_cmd.py", line 1504, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/madgraph/interface/extended_cmd.py", line 1453, in onecmd_orig
    return func(arg, **opt)
  File "/lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/madgraph/interface/amcatnlo_run_interface.py", line 1744, in do_launch
    self.compile(mode, options)
  File "/lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/madgraph/interface/amcatnlo_run_interface.py", line 5293, in compile
    self.check_tests(test, this_dir)
  File "/lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/madgraph/interface/amcatnlo_run_interface.py", line 5304, in check_tests
    return self.parse_test_mx_log(pjoin(dir, '%s.log' % test))
  File "/lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/madgraph/interface/amcatnlo_run_interface.py", line 5316, in parse_test_mx_log
    ' set to 0 in the param_card.dat.')
aMCatNLOError: Some tests failed, run cannot continue.
Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.
Value of current Options:
              text_editor : None
      notification_center : True
                    pjfry : None
       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 : ./HEPTools/lib
              auto_update : 7
             pythia8_path : None
                   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 : None
         madanalysis_path : None
                   lhapdf : lhapdf-config
             stdout_level : 20
                  nb_core : 12
            f2py_compiler : None
                    ninja : /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/HEPTools/lib
include_lepton_initiated_processes : False
                  amcfast : amcfast-config
       cluster_retry_wait : 300
    output_dependencies : external
           crash_on_error : False
mg5amc_py8_interface_path : None
         loop_color_flows : False
                  samurai : None
         cluster_nb_retry : 1
                 mg5_path : /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0
                  timeout : 60
                    gauge : Feynman
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0

Any idea of what could be going wrong?

Thanks a lot in advance!
Liqing

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
davide.pagani.85 Edit question
Solved by:
davide.pagani.85
Solved:
2018-10-13
Last query:
2018-10-13
Last reply:
2018-09-26

Hi Liqing,

Can you copy the content of the file

Cards/proc_card_mg5.dat

inside the folder you have generated for www production?

Thanks
Davide

Liqing Zhang (liqing.zhang) said : #2

Hi Davide,

Thanks a lot for your reply.

The following is the content.

set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set include_lepton_initiated_processes False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model loop_qcd_qed_sm
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 22 # pass to 5 flavors
define j = p
generate p p > w+ w- w+ [QED]
output TEST_WWW

Hi Liqing,

the problem is that the widths of the particles are all different from zero in the default param_card.dat, but you did not turn on the complex_mass_scheme (as you can see in your proc_card).

The code set the W width to zero, but it is not enough.
You are performing the calculation in the on-shell scheme, which for WWW production is correct, but you have to set all the widths to zero.

Alternatively you can turn on the complex mass scheme with

set complex_mass_scheme True

and set only the W width to zero. The numerical difference for the two approaches should be negligible.

Unless you are particularly interested in the renormalization of the EW interactions in the alpha(MZ) scheme, I suggest you to import the model

loop_qcd_qed_sm_Gmu

which employs the Gmu scheme and therefore is better for this process.

I hope this can help you, let me know if you still have problems.

Davide

Liqing Zhang (liqing.zhang) said : #4

Hi Davide,

Thanks a lot. it works now.

But I still have two questions:

1> when I import the model of loop_qcd_qed_sm_Gmu, I get a error

import loop_qcd_qed_sm_Gmu
The import format was not given, so we guess it as proc_v4
ERROR: I/O error on file loop_qcd_qed_sm_Gmu (2): No such file or directory
Command "import loop_qcd_qed_sm_Gmu" interrupted in sub-command:
"#import proc_v4 loop_qcd_qed_sm_Gmu" with error:
InvalidCmd : "loop_qcd_qed_sm_Gmu" is not a valid path

2>It's been 10 hours since I type 'launch', but it's still running. Is this normal?

Best,
Liqing

Hi Liqing,

try

import model loop_qcd_qed_sm_Gmu

if you do not have the model you need the internet connection for downloading it automatically.

With which precision are you running it and on how many cores?

With the default run_card

 0.01 = req_acc_FO

I have

  Inclusive timing profile :
    Overall slowest channel 0:40:17 (P0_dxu_wpwpwm/all_G1 [step 0])
    Average channel running time 0:17:49
    Aggregated total running time 2:22:36

Let me know
Cheers

Davide

Liqing Zhang (liqing.zhang) said : #6

Hi Davide,

Sorry for the late feedback. I was off line for a long period.

Thanks a lot for your explanation.

Here is a example which is submited two days before. But it’s not finished yet.
The precision I used is 0.01, and the output shows there are 12 cores.

The commands I used are :
     set complex_mass_scheme True
     generate p p > w+ w- w+ [QED]
     output WWW
     launch

The following is the output:
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'madspin=onshell' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>
INFO: will run in mode: NLO
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. FO_analyse : FO_analyse_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, FO_analyse, enter path][90s to answer]
>
INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
Not able to open file /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/bin/10.8-NLO/crossx.html since no program configured.Please set one in ./input/mg5_configuration.txt
INFO: Compiling the code
INFO: For gauge cancellation, the width of 'W+' has been set to zero.
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 12 cores
INFO: Compiling P0_udx_wpwpwm...
INFO: Compiling P0_dxu_wpwpwm...
INFO: P0_udx_wpwpwm done.
INFO: P0_dxu_wpwpwm done.
INFO: Checking test output:
INFO: P0_udx_wpwpwm
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxu_wpwpwm
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 12 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
INFO: Setting up grids
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 11h30 ]

Do you know the reason why it's so slow?

Thanks,
Liqing

________________________________________
From: <email address hidden> [<email address hidden>] on behalf of davide.pagani.85 [<email address hidden>]
Sent: Wednesday, September 26, 2018 23:57
To: Liqing Zhang
Subject: Re: [Question #674413]: aMCatNLOError : Some tests failed, run cannot continue.

Your question #674413 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/674413

    Status: Open => Answered

davide.pagani.85 proposed the following answer:
Hi Liqing,

try

import model loop_qcd_qed_sm_Gmu

if you do not have the model you need the internet connection for
downloading it automatically.

With which precision are you running it and on how many cores?

With the default run_card

 0.01 = req_acc_FO

I have

  Inclusive timing profile :
    Overall slowest channel 0:40:17 (P0_dxu_wpwpwm/all_G1 [step 0])
    Average channel running time 0:17:49
    Aggregated total running time 2:22:36

Let me know
Cheers

Davide

--
If this answers your question, please go to the following page to let us
know that it is solved:
https://answers.launchpad.net/mg5amcnlo/+question/674413/+confirm?answer_id=4

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/mg5amcnlo/+question/674413

You received this question notification because you asked the question.

Liqing Zhang (liqing.zhang) said : #7

Hi Davide,

The problem is solved.

The reason is that I found the Gmu scheme in this file 'models/loop_qcd_qed_sm/SMWidth/scheme.dat'
And I thought it's equivalent to what you mentioned about the loop_qcd_qed_sm_Gmu.

After the command ‘import model loop_qcd_qed_sm_Gmu’, everything is ok .

Thanks,
Liqing
________________________________________
From: Liqing Zhang
Sent: Wednesday, October 10, 2018 21:43
To: <email address hidden>
Subject: RE: [Question #674413]: aMCatNLOError : Some tests failed, run cannot continue.

Hi Davide,

Sorry for the late feedback. I was off line for a long period.

Thanks a lot for your explanation.

Here is a example which is submited two days before. But it’s not finished yet.
The precision I used is 0.01, and the output shows there are 12 cores.

The commands I used are :
     set complex_mass_scheme True
     generate p p > w+ w- w+ [QED]
     output WWW
     launch

The following is the output:
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'madspin=onshell' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>
INFO: will run in mode: NLO
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. FO_analyse : FO_analyse_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, FO_analyse, enter path][90s to answer]
>
INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
Not able to open file /lustre/AtlUser/lqzhang/Workarea/WWW/NLOEWK/V3.0/MG5_aMC_v3_0_0/bin/10.8-NLO/crossx.html since no program configured.Please set one in ./input/mg5_configuration.txt
INFO: Compiling the code
INFO: For gauge cancellation, the width of 'W+' has been set to zero.
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 12 cores
INFO: Compiling P0_udx_wpwpwm...
INFO: Compiling P0_dxu_wpwpwm...
INFO: P0_udx_wpwpwm done.
INFO: P0_dxu_wpwpwm done.
INFO: Checking test output:
INFO: P0_udx_wpwpwm
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxu_wpwpwm
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 12 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
INFO: Setting up grids
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 11h30 ]

Do you know the reason why it's so slow?

Thanks,
Liqing

________________________________________
From: <email address hidden> [<email address hidden>] on behalf of davide.pagani.85 [<email address hidden>]
Sent: Wednesday, September 26, 2018 23:57
To: Liqing Zhang
Subject: Re: [Question #674413]: aMCatNLOError : Some tests failed, run cannot continue.

Your question #674413 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/674413

    Status: Open => Answered

davide.pagani.85 proposed the following answer:
Hi Liqing,

try

import model loop_qcd_qed_sm_Gmu

if you do not have the model you need the internet connection for
downloading it automatically.

With which precision are you running it and on how many cores?

With the default run_card

 0.01 = req_acc_FO

I have

  Inclusive timing profile :
    Overall slowest channel 0:40:17 (P0_dxu_wpwpwm/all_G1 [step 0])
    Average channel running time 0:17:49
    Aggregated total running time 2:22:36

Let me know
Cheers

Davide

--
If this answers your question, please go to the following page to let us
know that it is solved:
https://answers.launchpad.net/mg5amcnlo/+question/674413/+confirm?answer_id=4

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/mg5amcnlo/+question/674413

You received this question notification because you asked the question.

Liqing Zhang (liqing.zhang) said : #8

Thanks davide.pagani.85, that solved my question.