How does MagGraph5 calculate cross-section when mad-spin is used with NLO ?

Asked by Yue Xu on 2019-11-11

Hello MadGraph group,

I am generating events with a modified HC model at NLO and madspin is used to decay W and Z bosons, here is the process(labeled as A )(x0 is a heavy higgs):
import model HC_NLO_X0_UFO
generate p p > z x0 > z e+ e- z /a [QCD]
add process p p > z x0 > z mu+ mu- z /a [QCD]
add process p p > z x0 > w+ e+ e- w- /a [QCD]
add process p p > z x0 > w+ mu+ mu- w- /a [QCD]
add process p p > w- x0 > w- e+ e- z / a [QCD]
add process p p > w- x0 > w- mu+ mu- z / a [QCD]
add process p p > w+ x0 > w+ e+ e- z /a [QCD]
add process p p > w+ x0 > w+ mu+ mu- z /a [QCD]

and the set of madspin_card is:
decay w+ > j j
decay w- > j j
decay z > j j
launch

the model can be fount at:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_NLO_X0_UFO
the process is here:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0
With the set of parameter in the param_card, the cross-section is 0.01143 from crossx.html.

Since the final states of this process have different branch ratio, I generated the final states with same branch ratio separately.
Process 1:
import model HC_NLO_X0_UFO
generate p p > z x0 > z e+ e- z /a [QCD]
add process p p > z x0 > z mu+ mu- z /a [QCD]
output PROCNLO_HC_NLO_X0_UFO_0_1
madspin:
decay z > j j
launch
I got the cs before decay is 5.239e-03 and the branch ratio is 0.294154, so the final cross-section is 0.00154107. This process is here:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0_1

Process 2:
import model HC_NLO_X0_UFO
generate p p > z x0 > w+ e+ e- w- /a [QCD]
add process p p > z x0 > w+ mu+ mu- w- /a [QCD]
output PROCNLO_HC_NLO_X0_UFO_0_2
madspin:
decay w+ > j j
decay w- > j j
launch
The branch ratio is 0.444518 and the final cross-section is 0.00283000. The process is here:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0_2

Process 3:
import model HC_NLO_X0_UFO
generate p p > w- x0 > w- e+ e- z / a [QCD]
add process p p > w- x0 > w- mu+ mu- z / a [QCD]
add process p p > w+ x0 > w+ e+ e- z /a [QCD]
add process p p > w+ x0 > w+ mu+ mu- z /a [QCD]
output PROCNLO_HC_NLO_X0_UFO_0_3
madspin:
decay w+ > j j
decay w- > j j
decay z > j j
launch

The branch ratio is 0.361603 and the final cross-section is 0.005093. The process is here:
https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0_3

If I add the cross-section of process 1, 2 and 3 together, the total cross section is 0.00946407 which is much smaller than 0.01143.
But the final state of process 1,2 and 3 is the same as A, they should have the same cross section.
How does this happen?
How can we get the correct cross section when processes with different branch ratio are included?

Cheers,
Yue

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Yue Xu
Solved:
2019-12-03
Last query:
2019-12-03
Last reply:
2019-11-28

This should be related to statistical fluctuation no?
Having different BR means that we have to throw away events and this will be a new source of statistical error.

Cheers,

Olivier

> On 11 Nov 2019, at 09:23, Yue Xu <email address hidden> wrote:
>
> New question #685744 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Hello MadGraph group,
>
> I am generating events with a modified HC model at NLO and madspin is used to decay W and Z bosons, here is the process(labeled as A )(x0 is a heavy higgs):
> import model HC_NLO_X0_UFO
> generate p p > z x0 > z e+ e- z /a [QCD]
> add process p p > z x0 > z mu+ mu- z /a [QCD]
> add process p p > z x0 > w+ e+ e- w- /a [QCD]
> add process p p > z x0 > w+ mu+ mu- w- /a [QCD]
> add process p p > w- x0 > w- e+ e- z / a [QCD]
> add process p p > w- x0 > w- mu+ mu- z / a [QCD]
> add process p p > w+ x0 > w+ e+ e- z /a [QCD]
> add process p p > w+ x0 > w+ mu+ mu- z /a [QCD]
>
> and the set of madspin_card is:
> decay w+ > j j
> decay w- > j j
> decay z > j j
> launch
>
> the model can be fount at:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/HC_NLO_X0_UFO
> the process is here:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0
> With the set of parameter in the param_card, the cross-section is 0.01143 from crossx.html.
>
> Since the final states of this process have different branch ratio, I generated the final states with same branch ratio separately.
> Process 1:
> import model HC_NLO_X0_UFO
> generate p p > z x0 > z e+ e- z /a [QCD]
> add process p p > z x0 > z mu+ mu- z /a [QCD]
> output PROCNLO_HC_NLO_X0_UFO_0_1
> madspin:
> decay z > j j
> launch
> I got the cs before decay is 5.239e-03 and the branch ratio is 0.294154, so the final cross-section is 0.00154107. This process is here:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0_1
>
> Process 2:
> import model HC_NLO_X0_UFO
> generate p p > z x0 > w+ e+ e- w- /a [QCD]
> add process p p > z x0 > w+ mu+ mu- w- /a [QCD]
> output PROCNLO_HC_NLO_X0_UFO_0_2
> madspin:
> decay w+ > j j
> decay w- > j j
> launch
> The branch ratio is 0.444518 and the final cross-section is 0.00283000. The process is here:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0_2
>
> Process 3:
> import model HC_NLO_X0_UFO
> generate p p > w- x0 > w- e+ e- z / a [QCD]
> add process p p > w- x0 > w- mu+ mu- z / a [QCD]
> add process p p > w+ x0 > w+ e+ e- z /a [QCD]
> add process p p > w+ x0 > w+ mu+ mu- z /a [QCD]
> output PROCNLO_HC_NLO_X0_UFO_0_3
> madspin:
> decay w+ > j j
> decay w- > j j
> decay z > j j
> launch
>
> The branch ratio is 0.361603 and the final cross-section is 0.005093. The process is here:
> https://github.com/xuyue1231/MG5-heavy-Higgs/tree/master/PROCNLO_HC_NLO_X0_UFO_0_3
>
> If I add the cross-section of process 1, 2 and 3 together, the total cross section is 0.00946407 which is much smaller than 0.01143.
> But the final state of process 1,2 and 3 is the same as A, they should have the same cross section.
> How does this happen?
> How can we get the correct cross section when processes with different branch ratio are included?
>
> Cheers,
> Yue
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #2

In case A, the error is 0.01143 ± 1.4e-05. But the cross section difference between case A and the case adding process1,2,3 together reach 20%, it's so large.
If I want a precision cross section, is it better to generate events separately?

Cheers,
Yue

Hi,

I have run the following two scripts:
generate p p > z j
add process p p > w+ j
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z j
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
#generate p p > z z /a
generate p p > w+ j
#add process p p > w- z / a
#add process p p > w+ z /a
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j

and

generate p p > z z
add process p p > w+ w-
add process p p > w- z
add process p p > w+ z
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z z
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
#generate p p > z z /a
generate p p > w+ w-
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > w- z
add process p p > w+ z
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j

In both case, I have the expected result after madspin. (variation at the level of 1 per 1000)

Also this might be an effect related to the negative number of events.
In the case where I do not have a unique BR, I need to re-calculate the cross-section
from the generated events but since negative events are present this can increase statistical error.

I have tested this on

generate p p > z j [QCD]
add process p p > w+ j [QCD]
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z j [QCD]
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
#generate p p > z z /a
generate p p > w+ j [QCD]
#add process p p > w- z / a
#add process p p > w+ z /a
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j

and here I have a larger difference. But this expected since this is due to the amount of negative weights in the sample (in this case ~ 25%)

Cheers,

Olivier

> On 11 Nov 2019, at 14:57, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> In case A, the error is 0.01143 ± 1.4e-05. But the cross section difference between case A and the case adding process1,2,3 together reach 20%, it's so large.
> If I want a precision cross section, is it better to generate events separately?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #4

Hi Olivier,

We want to re-scale the LO cross-section to NLO cross-section.
But in the NLO case there is a large difference because of the negative weights.

Is it reliable to use such a cross-section? or it's better to calculate the cross-section separately?

Cheers,
Yue

To have a more accurate value, it is indeed better to not use this combine mode.

Cheers,

Olivier

> On 26 Nov 2019, at 03:43, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi Olivier,
>
> We want to re-scale the LO cross-section to NLO cross-section.
> But in the NLO case there is a large difference because of the negative weights.
>
> Is it reliable to use such a cross-section? or it's better to calculate
> the cross-section separately?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #6

Hi,

I'm considering the interference of the processes with same initial and final state.
For the process like:
generate p p > z x0 > l+ l- j j j j
add process p p > w+ x0 > l- l+ j j j j

the two sub-processes may have interference.
But for the process like:
generate p p > z x0 > w- l+ l- w+ /a [QCD]
add process p p > w+ x0 > z l+ l- z /a [QCD]
manspin=ON
decay z > j j
decay w->j j
decay w+>j j

Will the two sub-processes have interference ?

Cheers,
Yue

Hi,

None of the above syntax will have interference term. (even the one without using madspin)

Cheers,

Olivier

> On 28 Nov 2019, at 08:47, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi,
>
> I'm considering the interference of the processes with same initial and final state.
> For the process like:
> generate p p > z x0 > l+ l- j j j j
> add process p p > w+ x0 > l- l+ j j j j
>
> the two sub-processes may have interference.
> But for the process like:
> generate p p > z x0 > w- l+ l- w+ /a [QCD]
> add process p p > w+ x0 > z l+ l- z /a [QCD]
> manspin=ON
> decay z > j j
> decay w->j j
> decay w+>j j
>
> Will the two sub-processes have interference ?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #8

Hi Olivier,

Thanks a lot!
Then, separating the processes can be a choice to get an accurate cross section value.

I met another problem when use this HC NLO model to generate events.
The process is (tau is added to the lepton definition, b quark is added to the p and j definition):
import model HC_NLO_X0_UFO
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
generate p p > z x0 > z l+ l- z /a [QCD]
add process p p > z x0 > w+ l+ l- w- /a [QCD]
add process p p > w- x0 > w- l+ l- z / a [QCD]
add process p p > w+ x0 > w+ l+ l- z /a [QCD]

I got the warning like "WARNING: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computation "

I don't know how to make it work.
Could you do me a favour?

Thanks,
Yue

P.S.
When generating events with "./bin/generate_events", I got the following output:

INFO: Idle: 31, Running: 48, Completed: 59 [ 1m 23s ]
INFO: Idle: 30, Running: 48, Completed: 60 [ 1m 24s ]
INFO: Idle: 29, Running: 48, Completed: 61 [ 1m 24s ]
INFO: Idle: 28, Running: 48, Completed: 62 [ 1m 25s ]
INFO: Idle: 27, Running: 48, Completed: 63 [ 1m 25s ]
INFO: Idle: 26, Running: 48, Completed: 64 [ 1m 25s ]
INFO: Idle: 25, Running: 48, Completed: 65 [ 1m 25s ]
INFO: Idle: 24, Running: 48, Completed: 66 [ 1m 25s ]
INFO: Idle: 23, Running: 48, Completed: 67 [ 1m 31s ]
INFO: Idle: 22, Running: 48, Completed: 68 [ 2m 2s ]
WARNING: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computation
INFO: Idle: 22, Running: 47, Completed: 69 [ 2m 23s ]
INFO: Idle: 22, Running: 46, Completed: 70 [ 2m 23s ]
INFO: Idle: 22, Running: 45, Completed: 71 [ 2m 23s ]
INFO: Idle: 22, Running: 44, Completed: 72 [ 2m 23s ]
INFO: Idle: 22, Running: 43, Completed: 73 [ 2m 24s ]
INFO: Idle: 22, Running: 42, Completed: 74 [ 2m 24s ]
INFO: Idle: 22, Running: 41, Completed: 75 [ 2m 24s ]
INFO: Idle: 22, Running: 40, Completed: 76 [ 2m 24s ]
INFO: Idle: 22, Running: 39, Completed: 77 [ 2m 24s ]
INFO: Idle: 22, Running: 38, Completed: 78 [ 2m 24s ]
INFO: Idle: 22, Running: 37, Completed: 79 [ 2m 24s ]
INFO: Idle: 22, Running: 36, Completed: 80 [ 2m 24s ]
INFO: Idle: 22, Running: 35, Completed: 81 [ 2m 24s ]
INFO: Idle: 22, Running: 34, Completed: 82 [ 2m 24s ]
INFO: Idle: 22, Running: 33, Completed: 83 [ 2m 24s ]
INFO: Idle: 22, Running: 32, Completed: 84 [ 2m 24s ]
INFO: Idle: 22, Running: 31, Completed: 85 [ 2m 24s ]
INFO: Idle: 22, Running: 30, Completed: 86 [ 2m 24s ]
INFO: Idle: 22, Running: 29, Completed: 87 [ 2m 24s ]
INFO: Idle: 22, Running: 28, Completed: 88 [ 2m 24s ]
INFO: Idle: 22, Running: 27, Completed: 89 [ 2m 24s ]
INFO: Idle: 22, Running: 26, Completed: 90 [ 2m 24s ]
INFO: Idle: 22, Running: 25, Completed: 91 [ 2m 24s ]
INFO: Idle: 22, Running: 24, Completed: 92 [ 2m 24s ]
INFO: Idle: 22, Running: 23, Completed: 93 [ 2m 24s ]
INFO: Idle: 22, Running: 22, Completed: 94 [ 2m 24s ]
INFO: Idle: 22, Running: 21, Completed: 95 [ 2m 24s ]
INFO: Idle: 22, Running: 20, Completed: 96 [ 2m 25s ]
INFO: Idle: 22, Running: 19, Completed: 97 [ 2m 25s ]
INFO: Idle: 22, Running: 18, Completed: 98 [ 2m 25s ]
INFO: Idle: 22, Running: 17, Completed: 99 [ 2m 25s ]
INFO: Idle: 22, Running: 16, Completed: 100 [ 2m 25s ]
INFO: Idle: 22, Running: 15, Completed: 101 [ 2m 25s ]
INFO: Idle: 22, Running: 14, Completed: 102 [ 2m 25s ]
INFO: Idle: 22, Running: 13, Completed: 103 [ 2m 25s ]
INFO: Idle: 22, Running: 12, Completed: 104 [ 2m 25s ]
INFO: Idle: 22, Running: 11, Completed: 105 [ 2m 25s ]
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_dux_wmx0_wmtaptamz_no_a/ajob1: line 34: 27591 Terminated ../madevent_mintMC > log.txt < input_app.txt 2>&1
INFO: Idle: 22, Running: 10, Completed: 106 [ 2m 25s ]
INFO: Idle: 22, Running: 9, Completed: 107 [ 2m 25s ]
INFO: Idle: 22, Running: 8, Completed: 108 [ 2m 25s ]
INFO: Idle: 22, Running: 7, Completed: 109 [ 2m 25s ]
INFO: Idle: 22, Running: 6, Completed: 110 [ 2m 25s ]
INFO: Idle: 22, Running: 5, Completed: 111 [ 2m 25s ]
INFO: Idle: 22, Running: 4, Completed: 112 [ 2m 25s ]
INFO: Idle: 22, Running: 3, Completed: 113 [ 2m 25s ]
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/ajob1: line 34: 27792 Terminated ../madevent_mintMC > log.txt < input_app.txt 2>&1
INFO: Idle: 22, Running: 2, Completed: 114 [ 2m 25s ]

and the debug log is here:
launch
Traceback (most recent call last):
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1679, in do_launch
    evt_file = self.run(mode, options)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1910, in run
    self.run_all_jobs(jobs_to_run,mint_step,fixed_order=False)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 2116, in run_all_jobs
    self.wait_for_complete(run_type)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 4657, in wait_for_complete
    self.cluster.wait(self.me_dir, update_status)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/cluster.py", line 844, in wait
    raise Exception, self.fail_msg
Exception: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computationlaunch
Traceback (most recent call last):
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1679, in do_launch
    evt_file = self.run(mode, options)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1910, in run
    self.run_all_jobs(jobs_to_run,mint_step,fixed_order=False)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 2116, in run_all_jobs
    self.wait_for_complete(run_type)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 4657, in wait_for_complete
    self.cluster.wait(self.me_dir, update_status)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/cluster.py", line 844, in wait
    raise Exception, self.fail_msg
Exception: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computation
Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
                  timeout : 60
       cluster_local_path : None
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lhapdf6/bin/lhapdf-config
             cluster_size : 100
           cluster_memory : None
                    pjfry : None
    cluster_status_update : (600, 30)
             cluster_time : None
            f2py_compiler : None
                    ninja : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
               hepmc_path : None
mg5amc_py8_interface_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/MG5aMC_PY8_interface
             pythia8_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/pythia8
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : None
          pythia-pgs_path : None
                 mg5_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6
                  td_path : None
                  collier : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
             delphes_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Delphes
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
         fortran_compiler : None
                  nb_core : 48
                      OLP : MadLoop
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                    golem : None
             cpp_compiler : None
      notification_center : True
                 run_mode : 2
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.6 2018-06-28 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
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 HC_NLO_X0_UFO
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
generate p p > z x0 > z l+ l- z /a [QCD]
add process p p > z x0 > w+ l+ l- w- /a [QCD]
add process p p > w- x0 > w- l+ l- z / a [QCD]
add process p p > w+ x0 > w+ l+ l- z /a [QCD]
output ./Total/NLO/HC_NLO_LHC_HH_2lep
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR FRBLOCK
###################################
BLOCK FRBLOCK #
      1 5.000000e+03 # lambda
      2 1.000000e+00 # cosa
      3 5.000000e-02 # ksm
      4 0.000000e+00 # khtt
      5 0.000000e+00 # katt
      6 0.000000e+00 # khbb
      7 0.000000e+00 # kabb
      8 0.000000e+00 # khll
      9 0.000000e+00 # kall
      10 0.000000e+00 # khaa
      11 0.000000e+00 # kaaa
      12 0.000000e+00 # khza
      13 0.000000e+00 # kaza
      16 2.908000e-01 # khzz
      17 0.000000e+00 # kazz
      18 5.248000e-01 # khww
      19 0.000000e+00 # kaww
      20 0.000000e+00 # khda
      21 -2.624000e-01 # khdz
      22 -2.624000e-01 # khdwr (real part of khdw)
      23 0.000000e+00 # khdwi (imaginary part of khdw)
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.325070e+02 # aewm1
      2 1.166390e-05 # gf
      3 1.180000e-01 # as
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      5 4.700000e+00 # mb
      6 1.730000e+02 # mt
      15 1.777000e+00 # mta
      23 9.118800e+01 # mz
      25 6.000000e+02 # mx0
      1 0.000000e+00 # d : 0.0
      2 0.000000e+00 # u : 0.0
      3 0.000000e+00 # s : 0.0
      4 0.000000e+00 # c : 0.0
      11 0.000000e+00 # e- : 0.0
      12 0.000000e+00 # ve : 0.0
      13 0.000000e+00 # mu- : 0.0
      14 0.000000e+00 # vm : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
      82 0.000000e+00 # gh : 0.0
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      5 4.700000e+00 # ymb
      6 1.730000e+02 # ymt
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR QNUMBERS 82
###################################
BLOCK QNUMBERS 82 # gh
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.491500e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.441404e+00
#
# PDG Width
DECAY 24 2.047600e+00
#
# PDG Width
DECAY 25 1.634610e-01
# BR NDA ID1 ID2 ...
   7.066517e-01 2 -24 24 # 0.115509993534
   2.933483e-01 2 23 23 # 0.0479510064663
#
# PDG Width
DECAY 82 0.000000e+00
#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# Some of the values of variables can be list. These can either be *
# comma or space separated. *
# *
# To display additional parameter, you can use the command: *
# update to_full *
#***********************************************************************
#
#*******************
# Running parameters
#*******************
#
#***********************************************************************
# Tag name for the run (one word) *
#***********************************************************************
  tag_1 = run_tag ! name of the run
#***********************************************************************
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#***********************************************************************
  10000 = nevents ! Number of unweighted events requested
  -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
  -1 = nevt_job ! Max number of events per job in event generation.
                 ! (-1= no split).
#***********************************************************************
# Normalize the weights of LHE events such that they sum or average to *
# the total cross section *
#***********************************************************************
  average = event_norm ! valid settings: average, sum, bias
#***********************************************************************
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#***********************************************************************
  0.01 = req_acc_fo ! Required accuracy (-1=ignored, and use the
                           ! number of points and iter. below)
# These numbers are ignored except if req_acc_FO is equal to -1
  5000 = npoints_fo_grid ! number of points to setup grids
  4 = niters_fo_grid ! number of iter. to setup grids
  10000 = npoints_fo ! number of points to compute Xsec
  6 = niters_fo ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
  0 = iseed ! rnd seed (0=assigned automatically=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
  1 = lpp1 ! beam 1 type (0 = no PDF)
  1 = lpp2 ! beam 2 type (0 = no PDF)
  6500.0 = ebeam1 ! beam 1 energy in GeV
  6500.0 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
  nn23nlo = pdlabel ! PDF set
  244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
              ! numbers for central PDF sets are allowed. Can be a list;
              ! PDF sets beyond the first are included via reweighting.
#***********************************************************************
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#***********************************************************************
  PYTHIA8 = parton_shower
  1.0 = shower_scale_factor ! multiply default shower starting
                                  ! scale by this factor
#***********************************************************************
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses divided by two of all final state particles *
# and partons. This can be changed in SubProcesses/set_scales.f or via *
# dynamical_scale_choice option) *
#***********************************************************************
  False = fixed_ren_scale ! if .true. use fixed ren scale
  False = fixed_fac_scale ! if .true. use fixed fac scale
  91.118 = mur_ref_fixed ! fixed ren reference scale
  91.118 = muf_ref_fixed ! fixed fact reference scale
  3 = dynamical_scale_choice ! Choose one (or more) of the predefined
           ! dynamical choices. Can be a list; scale choices beyond the
           ! first are included via reweighting
  1.0 = mur_over_ref ! ratio of current muR over reference muR
  1.0 = muf_over_ref ! ratio of current muF over reference muF
#***********************************************************************
# Reweight variables for scale dependence and PDF uncertainty *
#***********************************************************************
  1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
  1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
  True = reweight_scale ! Reweight to get scale variation using the
            ! rw_rscale and rw_fscale factors. Should be a list of
            ! booleans of equal length to dynamical_scale_choice to
            ! specify for which choice to include scale dependence.
  False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a
            ! list booleans of equal length to lhaid to specify for
            ! which PDF set to include the uncertainties.
#***********************************************************************
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#***********************************************************************
  False = store_rwgt_info ! Store info for reweighting in LHE file
#***********************************************************************
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#***********************************************************************
  0 = ickkw
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
# written in the LHE event file *
#***********************************************************************
  15.0 = bwcutoff
#***********************************************************************
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
  0.7 = jetradius ! The radius parameter for the jet algorithm
  10.0 = ptj ! Min jet transverse momentum
  -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
#***********************************************************************
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  0.0 = ptl ! Min lepton transverse momentum
  -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
  0.0 = drll ! Min distance between opposite sign lepton pairs
  0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
  0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
  30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
# all the other parameters are ignored. *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  20.0 = ptgmin ! Min photon transverse momentum
  -1.0 = etagamma ! Max photon abs(pseudo-rap)
  0.4 = r0gamma ! Radius of isolation code
  1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
  1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
  True = isoem ! isolate photons from EM energy (photons and leptons)
#***********************************************************************
# Cuts associated to MASSIVE particles identified by their PDG codes. *
# All cuts are applied to both particles and anti-particles, so use *
# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
# {6:100, 25:200} for multiple particles *
#***********************************************************************
  {} = pt_min_pdg ! Min pT for a massive particle
  {} = pt_max_pdg ! Max pT for a massive particle
  {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
#***********************************************************************
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
#***********************************************************************
  0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
#***********************************************************************

Hi,

> b quark is added to the p and j definition):

This does not make sense.
Either your computation is 4 flavor (with b massive)
or 5 flavour (and b massless)

MG5amC, will automatically adapt the definition of p/j according to the fact that your model defines the b quark as massive or massless.

This being said, such computation is not going to work in the 5 flavour scheme due to the possibility to have a top onshell. (i.e. you need. to use diagram substraction for that which is available as a plugin but I do not know how that one work)

Cheers,

Olivier

> On 28 Nov 2019, at 13:47, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi Olivier,
>
> Thanks a lot!
> Then, separating the processes can be a choice to get an accurate cross section value.
>
> I met another problem when use this HC NLO model to generate events.
> The process is (tau is added to the lepton definition, b quark is added to the p and j definition):
> import model HC_NLO_X0_UFO
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = g u c d s b u~ c~ d~ s~ b~
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> generate p p > z x0 > z l+ l- z /a [QCD]
> add process p p > z x0 > w+ l+ l- w- /a [QCD]
> add process p p > w- x0 > w- l+ l- z / a [QCD]
> add process p p > w+ x0 > w+ l+ l- z /a [QCD]
>
> I got the warning like "WARNING: program
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1
> 2 F 0 0 launch ends with non zero status: 1. Stop all computation "
>
> I don't know how to make it work.
> Could you do me a favour?
>
> Thanks,
> Yue
>
> P.S.
> When generating events with "./bin/generate_events", I got the following output:
>
> INFO: Idle: 31, Running: 48, Completed: 59 [ 1m 23s ]
> INFO: Idle: 30, Running: 48, Completed: 60 [ 1m 24s ]
> INFO: Idle: 29, Running: 48, Completed: 61 [ 1m 24s ]
> INFO: Idle: 28, Running: 48, Completed: 62 [ 1m 25s ]
> INFO: Idle: 27, Running: 48, Completed: 63 [ 1m 25s ]
> INFO: Idle: 26, Running: 48, Completed: 64 [ 1m 25s ]
> INFO: Idle: 25, Running: 48, Completed: 65 [ 1m 25s ]
> INFO: Idle: 24, Running: 48, Completed: 66 [ 1m 25s ]
> INFO: Idle: 23, Running: 48, Completed: 67 [ 1m 31s ]
> INFO: Idle: 22, Running: 48, Completed: 68 [ 2m 2s ]
> WARNING: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computation
> INFO: Idle: 22, Running: 47, Completed: 69 [ 2m 23s ]
> INFO: Idle: 22, Running: 46, Completed: 70 [ 2m 23s ]
> INFO: Idle: 22, Running: 45, Completed: 71 [ 2m 23s ]
> INFO: Idle: 22, Running: 44, Completed: 72 [ 2m 23s ]
> INFO: Idle: 22, Running: 43, Completed: 73 [ 2m 24s ]
> INFO: Idle: 22, Running: 42, Completed: 74 [ 2m 24s ]
> INFO: Idle: 22, Running: 41, Completed: 75 [ 2m 24s ]
> INFO: Idle: 22, Running: 40, Completed: 76 [ 2m 24s ]
> INFO: Idle: 22, Running: 39, Completed: 77 [ 2m 24s ]
> INFO: Idle: 22, Running: 38, Completed: 78 [ 2m 24s ]
> INFO: Idle: 22, Running: 37, Completed: 79 [ 2m 24s ]
> INFO: Idle: 22, Running: 36, Completed: 80 [ 2m 24s ]
> INFO: Idle: 22, Running: 35, Completed: 81 [ 2m 24s ]
> INFO: Idle: 22, Running: 34, Completed: 82 [ 2m 24s ]
> INFO: Idle: 22, Running: 33, Completed: 83 [ 2m 24s ]
> INFO: Idle: 22, Running: 32, Completed: 84 [ 2m 24s ]
> INFO: Idle: 22, Running: 31, Completed: 85 [ 2m 24s ]
> INFO: Idle: 22, Running: 30, Completed: 86 [ 2m 24s ]
> INFO: Idle: 22, Running: 29, Completed: 87 [ 2m 24s ]
> INFO: Idle: 22, Running: 28, Completed: 88 [ 2m 24s ]
> INFO: Idle: 22, Running: 27, Completed: 89 [ 2m 24s ]
> INFO: Idle: 22, Running: 26, Completed: 90 [ 2m 24s ]
> INFO: Idle: 22, Running: 25, Completed: 91 [ 2m 24s ]
> INFO: Idle: 22, Running: 24, Completed: 92 [ 2m 24s ]
> INFO: Idle: 22, Running: 23, Completed: 93 [ 2m 24s ]
> INFO: Idle: 22, Running: 22, Completed: 94 [ 2m 24s ]
> INFO: Idle: 22, Running: 21, Completed: 95 [ 2m 24s ]
> INFO: Idle: 22, Running: 20, Completed: 96 [ 2m 25s ]
> INFO: Idle: 22, Running: 19, Completed: 97 [ 2m 25s ]
> INFO: Idle: 22, Running: 18, Completed: 98 [ 2m 25s ]
> INFO: Idle: 22, Running: 17, Completed: 99 [ 2m 25s ]
> INFO: Idle: 22, Running: 16, Completed: 100 [ 2m 25s ]
> INFO: Idle: 22, Running: 15, Completed: 101 [ 2m 25s ]
> INFO: Idle: 22, Running: 14, Completed: 102 [ 2m 25s ]
> INFO: Idle: 22, Running: 13, Completed: 103 [ 2m 25s ]
> INFO: Idle: 22, Running: 12, Completed: 104 [ 2m 25s ]
> INFO: Idle: 22, Running: 11, Completed: 105 [ 2m 25s ]
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_dux_wmx0_wmtaptamz_no_a/ajob1: line 34: 27591 Terminated ../madevent_mintMC > log.txt < input_app.txt 2>&1
> INFO: Idle: 22, Running: 10, Completed: 106 [ 2m 25s ]
> INFO: Idle: 22, Running: 9, Completed: 107 [ 2m 25s ]
> INFO: Idle: 22, Running: 8, Completed: 108 [ 2m 25s ]
> INFO: Idle: 22, Running: 7, Completed: 109 [ 2m 25s ]
> INFO: Idle: 22, Running: 6, Completed: 110 [ 2m 25s ]
> INFO: Idle: 22, Running: 5, Completed: 111 [ 2m 25s ]
> INFO: Idle: 22, Running: 4, Completed: 112 [ 2m 25s ]
> INFO: Idle: 22, Running: 3, Completed: 113 [ 2m 25s ]
> /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/ajob1: line 34: 27792 Terminated ../madevent_mintMC > log.txt < input_app.txt 2>&1
> INFO: Idle: 22, Running: 2, Completed: 114 [ 2m 25s ]
>
>
> and the debug log is here:
> launch
> Traceback (most recent call last):
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1514, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
> return func(arg, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1679, in do_launch
> evt_file = self.run(mode, options)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1910, in run
> self.run_all_jobs(jobs_to_run,mint_step,fixed_order=False)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 2116, in run_all_jobs
> self.wait_for_complete(run_type)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 4657, in wait_for_complete
> self.cluster.wait(self.me_dir, update_status)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/cluster.py", line 844, in wait
> raise Exception, self.fail_msg
> Exception: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computationlaunch
> Traceback (most recent call last):
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1514, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
> return func(arg, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1679, in do_launch
> evt_file = self.run(mode, options)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1910, in run
> self.run_all_jobs(jobs_to_run,mint_step,fixed_order=False)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 2116, in run_all_jobs
> self.wait_for_complete(run_type)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 4657, in wait_for_complete
> self.cluster.wait(self.me_dir, update_status)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/cluster.py", line 844, in wait
> raise Exception, self.fail_msg
> Exception: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_ddx_zx0_wptaptamwm_no_a/ajob1 2 F 0 0 launch ends with non zero status: 1. Stop all computation
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_local_path : None
> cluster_queue : None
> madanalysis_path : None
> lhapdf : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lhapdf6/bin/lhapdf-config
> cluster_size : 100
> cluster_memory : None
> pjfry : None
> cluster_status_update : (600, 30)
> cluster_time : None
> f2py_compiler : None
> ninja : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
> hepmc_path : None
> mg5amc_py8_interface_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/MG5aMC_PY8_interface
> pythia8_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/pythia8
> hwpp_path : None
> automatic_html_opening : False
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> mg5_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6
> td_path : None
> collier : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
> delphes_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Delphes
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
> fortran_compiler : None
> nb_core : 48
> OLP : MadLoop
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> notification_center : True
> run_mode : 2
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.6 2018-06-28 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set default_unset_couplings 99
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> 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 HC_NLO_X0_UFO
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = g u c d s b u~ c~ d~ s~ b~
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> generate p p > z x0 > z l+ l- z /a [QCD]
> add process p p > z x0 > w+ l+ l- w- /a [QCD]
> add process p p > w- x0 > w- l+ l- z / a [QCD]
> add process p p > w+ x0 > w+ l+ l- z /a [QCD]
> output ./Total/NLO/HC_NLO_LHC_HH_2lep
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR FRBLOCK
> ###################################
> BLOCK FRBLOCK #
> 1 5.000000e+03 # lambda
> 2 1.000000e+00 # cosa
> 3 5.000000e-02 # ksm
> 4 0.000000e+00 # khtt
> 5 0.000000e+00 # katt
> 6 0.000000e+00 # khbb
> 7 0.000000e+00 # kabb
> 8 0.000000e+00 # khll
> 9 0.000000e+00 # kall
> 10 0.000000e+00 # khaa
> 11 0.000000e+00 # kaaa
> 12 0.000000e+00 # khza
> 13 0.000000e+00 # kaza
> 16 2.908000e-01 # khzz
> 17 0.000000e+00 # kazz
> 18 5.248000e-01 # khww
> 19 0.000000e+00 # kaww
> 20 0.000000e+00 # khda
> 21 -2.624000e-01 # khdz
> 22 -2.624000e-01 # khdwr (real part of khdw)
> 23 0.000000e+00 # khdwi (imaginary part of khdw)
> ###################################
> ## INFORMATION FOR LOOP
> ###################################
> BLOCK LOOP #
> 1 9.118800e+01 # mu_r
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.325070e+02 # aewm1
> 2 1.166390e-05 # gf
> 3 1.180000e-01 # as
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 5 4.700000e+00 # mb
> 6 1.730000e+02 # mt
> 15 1.777000e+00 # mta
> 23 9.118800e+01 # mz
> 25 6.000000e+02 # mx0
> 1 0.000000e+00 # d : 0.0
> 2 0.000000e+00 # u : 0.0
> 3 0.000000e+00 # s : 0.0
> 4 0.000000e+00 # c : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 0.0
> 14 0.000000e+00 # vm : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
> 82 0.000000e+00 # gh : 0.0
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
> ###################################
> ## INFORMATION FOR QNUMBERS 82
> ###################################
> BLOCK QNUMBERS 82 # gh
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2s+1)
> 3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> #
> #*************************
> # Decay widths *
> #*************************
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.491500e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.441404e+00
> #
> # PDG Width
> DECAY 24 2.047600e+00
> #
> # PDG Width
> DECAY 25 1.634610e-01
> # BR NDA ID1 ID2 ...
> 7.066517e-01 2 -24 24 # 0.115509993534
> 2.933483e-01 2 23 23 # 0.0479510064663
> #
> # PDG Width
> DECAY 82 0.000000e+00
> #***********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat aMC@NLO *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a hash (#) are info or comments *
> # *
> # mind the format: value = variable ! comment *
> # *
> # Some of the values of variables can be list. These can either be *
> # comma or space separated. *
> # *
> # To display additional parameter, you can use the command: *
> # update to_full *
> #***********************************************************************
> #
> #*******************
> # Running parameters
> #*******************
> #
> #***********************************************************************
> # Tag name for the run (one word) *
> #***********************************************************************
> tag_1 = run_tag ! name of the run
> #***********************************************************************
> # Number of LHE events (and their normalization) and the required *
> # (relative) accuracy on the Xsec. *
> # These values are ignored for fixed order runs *
> #***********************************************************************
> 10000 = nevents ! Number of unweighted events requested
> -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
> -1 = nevt_job ! Max number of events per job in event generation.
> ! (-1= no split).
> #***********************************************************************
> # Normalize the weights of LHE events such that they sum or average to *
> # the total cross section *
> #***********************************************************************
> average = event_norm ! valid settings: average, sum, bias
> #***********************************************************************
> # Number of points per itegration channel (ignored for aMC@NLO runs) *
> #***********************************************************************
> 0.01 = req_acc_fo ! Required accuracy (-1=ignored, and use the
> ! number of points and iter. below)
> # These numbers are ignored except if req_acc_FO is equal to -1
> 5000 = npoints_fo_grid ! number of points to setup grids
> 4 = niters_fo_grid ! number of iter. to setup grids
> 10000 = npoints_fo ! number of points to compute Xsec
> 6 = niters_fo ! number of iter. to compute Xsec
> #***********************************************************************
> # Random number seed *
> #***********************************************************************
> 0 = iseed ! rnd seed (0=assigned automatically=default))
> #***********************************************************************
> # Collider type and energy *
> #***********************************************************************
> 1 = lpp1 ! beam 1 type (0 = no PDF)
> 1 = lpp2 ! beam 2 type (0 = no PDF)
> 6500.0 = ebeam1 ! beam 1 energy in GeV
> 6500.0 = ebeam2 ! beam 2 energy in GeV
> #***********************************************************************
> # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> #***********************************************************************
> nn23nlo = pdlabel ! PDF set
> 244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
> ! numbers for central PDF sets are allowed. Can be a list;
> ! PDF sets beyond the first are included via reweighting.
> #***********************************************************************
> # Include the NLO Monte Carlo subtr. terms for the following parton *
> # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
> # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
> #***********************************************************************
> PYTHIA8 = parton_shower
> 1.0 = shower_scale_factor ! multiply default shower starting
> ! scale by this factor
> #***********************************************************************
> # Renormalization and factorization scales *
> # (Default functional form for the non-fixed scales is the sum of *
> # the transverse masses divided by two of all final state particles *
> # and partons. This can be changed in SubProcesses/set_scales.f or via *
> # dynamical_scale_choice option) *
> #***********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.118 = mur_ref_fixed ! fixed ren reference scale
> 91.118 = muf_ref_fixed ! fixed fact reference scale
> 3 = dynamical_scale_choice ! Choose one (or more) of the predefined
> ! dynamical choices. Can be a list; scale choices beyond the
> ! first are included via reweighting
> 1.0 = mur_over_ref ! ratio of current muR over reference muR
> 1.0 = muf_over_ref ! ratio of current muF over reference muF
> #***********************************************************************
> # Reweight variables for scale dependence and PDF uncertainty *
> #***********************************************************************
> 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
> 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
> True = reweight_scale ! Reweight to get scale variation using the
> ! rw_rscale and rw_fscale factors. Should be a list of
> ! booleans of equal length to dynamical_scale_choice to
> ! specify for which choice to include scale dependence.
> False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a
> ! list booleans of equal length to lhaid to specify for
> ! which PDF set to include the uncertainties.
> #***********************************************************************
> # Store reweight information in the LHE file for off-line model- *
> # parameter reweighting at NLO+PS accuracy *
> #***********************************************************************
> False = store_rwgt_info ! Store info for reweighting in LHE file
> #***********************************************************************
> # ickkw parameter: *
> # 0: No merging *
> # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
> # level. After showering an MLM-type merging should be applied as *
> # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
> # 4: UNLOPS merging (with pythia8 only). No interface from within *
> # MG5_aMC available, but available in Pythia8. *
> # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
> #***********************************************************************
> 0 = ickkw
> #***********************************************************************
> #
> #***********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
> # written in the LHE event file *
> #***********************************************************************
> 15.0 = bwcutoff
> #***********************************************************************
> # Cuts on the jets. Jet clustering is performed by FastJet. *
> # - When matching to a parton shower, these generation cuts should be *
> # considerably softer than the analysis cuts. *
> # - More specific cuts can be specified in SubProcesses/cuts.f *
> #***********************************************************************
> 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 0.7 = jetradius ! The radius parameter for the jet algorithm
> 10.0 = ptj ! Min jet transverse momentum
> -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
> #***********************************************************************
> # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
> # More specific cuts can be specified in SubProcesses/cuts.f *
> #***********************************************************************
> 0.0 = ptl ! Min lepton transverse momentum
> -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 0.0 = drll ! Min distance between opposite sign lepton pairs
> 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
> #***********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
> # all the other parameters are ignored. *
> # More specific cuts can be specified in SubProcesses/cuts.f *
> #***********************************************************************
> 20.0 = ptgmin ! Min photon transverse momentum
> -1.0 = etagamma ! Max photon abs(pseudo-rap)
> 0.4 = r0gamma ! Radius of isolation code
> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> True = isoem ! isolate photons from EM energy (photons and leptons)
> #***********************************************************************
> # Cuts associated to MASSIVE particles identified by their PDG codes. *
> # All cuts are applied to both particles and anti-particles, so use *
> # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
> # {6:100, 25:200} for multiple particles *
> #***********************************************************************
> {} = pt_min_pdg ! Min pT for a massive particle
> {} = pt_max_pdg ! Max pT for a massive particle
> {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
> #***********************************************************************
> # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> #***********************************************************************
> 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> #***********************************************************************
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #10

Hi,

Even without the b quark, I still failed to generate the events.

The output is:
INFO: Compiling directories...
INFO: Compiling on 48 cores
INFO: Compiling P0_uux_zx0_ztaptamz_no_a...
INFO: Compiling P0_ddx_zx0_ztaptamz_no_a...
INFO: Compiling P0_uxu_zx0_ztaptamz_no_a...
Fontconfig warning: ignoring UTF-8: not a valid region tag
INFO: Compiling P0_dxd_zx0_ztaptamz_no_a...
INFO: Compiling P0_uux_zx0_wptaptamwm_no_a...
INFO: Compiling P0_ddx_zx0_wptaptamwm_no_a...
INFO: Compiling P0_uxu_zx0_wptaptamwm_no_a...
INFO: Compiling P0_dxd_zx0_wptaptamwm_no_a...
INFO: Compiling P0_dux_wmx0_wmtaptamz_no_a...
INFO: Compiling P0_uxd_wmx0_wmtaptamz_no_a...
INFO: Compiling P0_udx_wpx0_wptaptamz_no_a...
INFO: Compiling P0_dxu_wpx0_wptaptamz_no_a...
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
INFO: P0_uxu_zx0_wptaptamwm_no_a done.
INFO: P0_dxu_wpx0_wptaptamz_no_a done.
INFO: P0_uxd_wmx0_wmtaptamz_no_a done.
INFO: P0_ddx_zx0_ztaptamz_no_a done.
INFO: P0_dux_wmx0_wmtaptamz_no_a done.
INFO: P0_uux_zx0_wptaptamwm_no_a done.
INFO: P0_uux_zx0_ztaptamz_no_a done.
INFO: P0_udx_wpx0_wptaptamz_no_a done.
INFO: P0_ddx_zx0_wptaptamwm_no_a done.
INFO: P0_uxu_zx0_ztaptamz_no_a done.
INFO: P0_dxd_zx0_ztaptamz_no_a done.
INFO: P0_dxd_zx0_wptaptamwm_no_a done.
INFO: Checking test output:
INFO: P0_uux_zx0_ztaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_ddx_zx0_ztaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_uxu_zx0_ztaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_dxd_zx0_ztaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_uux_zx0_wptaptamwm_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_ddx_zx0_wptaptamwm_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_uxu_zx0_wptaptamwm_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_dxd_zx0_wptaptamwm_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_dux_wmx0_wmtaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_uxd_wmx0_wmtaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_udx_wpx0_wptaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_dxu_wpx0_wptaptamz_no_a
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: Starting run
INFO: Using 48 cores
INFO: Cleaning previous results
INFO: Doing NLO matched to parton shower
INFO: Setting up grids
INFO: Idle: 4, Running: 48, Completed: 0 [ current time: 21h20 ]
WARNING: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/ajob1 5 F 0 0 launch ends with non zero status: 1. Stop all computation
INFO: Idle: 4, Running: 47, Completed: 1 [ 3.3s ]
INFO: Idle: 4, Running: 46, Completed: 2 [ 3.3s ]
INFO: Idle: 4, Running: 43, Completed: 5 [ 3.4s ]
INFO: Idle: 4, Running: 42, Completed: 6 [ 3.4s ]
INFO: Idle: 4, Running: 40, Completed: 8 [ 3.5s ]
INFO: Idle: 4, Running: 39, Completed: 9 [ 3.5s ]
INFO: Idle: 4, Running: 37, Completed: 11 [ 3.6s ]
INFO: Idle: 4, Running: 36, Completed: 12 [ 3.6s ]
INFO: Idle: 4, Running: 35, Completed: 13 [ 3.6s ]
INFO: Idle: 4, Running: 34, Completed: 14 [ 3.7s ]
INFO: Idle: 4, Running: 33, Completed: 15 [ 3.7s ]
INFO: Idle: 4, Running: 32, Completed: 16 [ 3.8s ]
INFO: Idle: 4, Running: 31, Completed: 17 [ 3.8s ]
INFO: Idle: 4, Running: 30, Completed: 18 [ 3.8s ]
INFO: Idle: 4, Running: 29, Completed: 19 [ 3.9s ]
INFO: Idle: 4, Running: 28, Completed: 20 [ 3.9s ]
INFO: Idle: 4, Running: 27, Completed: 21 [ 3.9s ]
INFO: Idle: 4, Running: 26, Completed: 22 [ 4s ]
INFO: Idle: 4, Running: 25, Completed: 23 [ 4s ]
INFO: Idle: 4, Running: 24, Completed: 24 [ 4s ]

And the debug log file is:
launch
Traceback (most recent call last):
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1687, in do_launch
    evt_file = self.run(mode, options)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1918, in run
    self.run_all_jobs(jobs_to_run,mint_step,fixed_order=False)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 2124, in run_all_jobs
    self.wait_for_complete(run_type)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 4692, in wait_for_complete
    self.cluster.wait(self.me_dir, update_status)
  File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/cluster.py", line 844, in wait
    raise Exception, self.fail_msg
Exception: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/ajob1 5 F 0 0 launch ends with non zero status: 1. Stop all computation
Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
                  timeout : 60
       cluster_local_path : None
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lhapdf6/bin/lhapdf-config
             cluster_size : 100
           cluster_memory : None
                    pjfry : None
    cluster_status_update : (600, 30)
             cluster_time : None
            f2py_compiler : None
                    ninja : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
               hepmc_path : None
mg5amc_py8_interface_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/MG5aMC_PY8_interface
             pythia8_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/pythia8
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : None
          pythia-pgs_path : None
                 mg5_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6
                  td_path : None
                  collier : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
             delphes_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Delphes
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
         fortran_compiler : None
                  nb_core : 48
                      OLP : MadLoop
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                    golem : None
             cpp_compiler : None
      notification_center : True
                 run_mode : 2
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.7 2019-10-16 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
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 HC_NLO_X0_UFO
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
generate p p > z x0 > z l+ l- z /a [QCD]
add process p p > z x0 > w+ l+ l- w- /a [QCD]
add process p p > w- x0 > w- l+ l- z / a [QCD]
add process p p > w+ x0 > w+ l+ l- z /a [QCD]
output ./Total/NLO/HC_NLO_LHC_HH_2lep
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR FRBLOCK
###################################
BLOCK FRBLOCK #
      1 5.000000e+03 # lambda
      2 1.000000e+00 # cosa
      3 1.000000e+00 # ksm
      4 0.000000e+00 # khtt
      5 0.000000e+00 # katt
      6 0.000000e+00 # khbb
      7 0.000000e+00 # kabb
      8 0.000000e+00 # khll
      9 0.000000e+00 # kall
      10 0.000000e+00 # khaa
      11 0.000000e+00 # kaaa
      12 0.000000e+00 # khza
      13 0.000000e+00 # kaza
      16 2.908000e-01 # khzz
      17 0.000000e+00 # kazz
      18 5.248000e-01 # khww
      19 0.000000e+00 # kaww
      20 0.000000e+00 # khda
      21 -2.624000e-01 # khdz
      22 -2.624000e-01 # khdwr (real part of khdw)
      23 0.000000e+00 # khdwi (imaginary part of khdw)
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.325070e+02 # aewm1
      2 1.166390e-05 # gf
      3 1.180000e-01 # as
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      5 4.700000e+00 # mb
      6 1.730000e+02 # mt
      15 1.777000e+00 # mta
      23 9.118800e+01 # mz
      25 1.250000e+02 # mx0
      1 0.000000e+00 # d : 0.0
      2 0.000000e+00 # u : 0.0
      3 0.000000e+00 # s : 0.0
      4 0.000000e+00 # c : 0.0
      11 0.000000e+00 # e- : 0.0
      12 0.000000e+00 # ve : 0.0
      13 0.000000e+00 # mu- : 0.0
      14 0.000000e+00 # vm : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
      82 0.000000e+00 # gh : 0.0
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      5 4.700000e+00 # ymb
      6 1.730000e+02 # ymt
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
BLOCK QNUMBERS 82 # gh
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.491500e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.441404e+00
#
# PDG Width
DECAY 24 2.047600e+00
#
# PDG Width
DECAY 25 8.485861e-04
# BR NDA ID1 ID2 ...
   1.519292e-01 3 -2 1 24 # 0.000128925007304
   1.519292e-01 3 -4 3 24 # 0.000128925007304
   1.504738e-01 3 -24 -1 2 # 0.000127689975094
   1.504738e-01 3 -24 -3 4 # 0.000127689975094
   5.024122e-02 3 -12 11 24 # 4.2634000939e-05
   5.024122e-02 3 -14 13 24 # 4.2634000939e-05
   5.007860e-02 3 -24 -11 12 # 4.24960038675e-05
   5.007860e-02 3 -24 -13 14 # 4.24960038675e-05
   4.968264e-02 3 -16 15 24 # 4.21599977153e-05
   4.961429e-02 3 -24 -15 16 # 4.21019968554e-05
   1.488829e-02 3 -1 1 23 # 1.26339959468e-05
   1.488829e-02 3 -3 3 23 # 1.26339959468e-05
   1.194340e-02 3 -5 5 23 # 1.01350032267e-05
   1.176015e-02 3 -2 2 23 # 9.97949982391e-06
   1.176015e-02 3 -4 4 23 # 9.97949982391e-06
   6.692701e-03 3 -12 12 23 # 5.67933304006e-06
   6.692701e-03 3 -14 14 23 # 5.67933304006e-06
   6.692701e-03 3 -16 16 23 # 5.67933304006e-06
   3.333663e-03 3 -11 11 23 # 2.82890008388e-06
   3.333663e-03 3 -13 13 23 # 2.82890008388e-06
   3.271677e-03 3 -15 15 23 # 2.77629962589e-06
#
# PDG Width
DECAY 82 0.000000e+00
#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# Some of the values of variables can be list. These can either be *
# comma or space separated. *
# *
# To display additional parameter, you can use the command: *
# update to_full *
#***********************************************************************
#
#*******************
# Running parameters
#*******************
#
#***********************************************************************
# Tag name for the run (one word) *
#***********************************************************************
  tag_1 = run_tag ! name of the run
#***********************************************************************
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#***********************************************************************
  10000 = nevents ! Number of unweighted events requested
  -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
  -1 = nevt_job ! Max number of events per job in event generation.
                 ! (-1= no split).
#***********************************************************************
# Normalize the weights of LHE events such that they sum or average to *
# the total cross section *
#***********************************************************************
  average = event_norm ! valid settings: average, sum, bias
#***********************************************************************
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#***********************************************************************
  0.01 = req_acc_fo ! Required accuracy (-1=ignored, and use the
                     ! number of points and iter. below)
# These numbers are ignored except if req_acc_FO is equal to -1
  5000 = npoints_fo_grid ! number of points to setup grids
  4 = niters_fo_grid ! number of iter. to setup grids
  10000 = npoints_fo ! number of points to compute Xsec
  6 = niters_fo ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
  0 = iseed ! rnd seed (0=assigned automatically=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
  1 = lpp1 ! beam 1 type (0 = no PDF)
  1 = lpp2 ! beam 2 type (0 = no PDF)
  6500.0 = ebeam1 ! beam 1 energy in GeV
  6500.0 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
  nn23nlo = pdlabel ! PDF set
  244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
              ! numbers for central PDF sets are allowed. Can be a list;
              ! PDF sets beyond the first are included via reweighting.
#***********************************************************************
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#***********************************************************************
  PYTHIA8 = parton_shower
  1.0 = shower_scale_factor ! multiply default shower starting
                                  ! scale by this factor
#***********************************************************************
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses divided by two of all final state particles *
# and partons. This can be changed in SubProcesses/set_scales.f or via *
# dynamical_scale_choice option) *
#***********************************************************************
  False = fixed_ren_scale ! if .true. use fixed ren scale
  False = fixed_fac_scale ! if .true. use fixed fac scale
  91.118 = mur_ref_fixed ! fixed ren reference scale
  91.118 = muf_ref_fixed ! fixed fact reference scale
  3 = dynamical_scale_choice ! Choose one (or more) of the predefined
           ! dynamical choices. Can be a list; scale choices beyond the
           ! first are included via reweighting
  1.0 = mur_over_ref ! ratio of current muR over reference muR
  1.0 = muf_over_ref ! ratio of current muF over reference muF
#***********************************************************************
# Reweight variables for scale dependence and PDF uncertainty *
#***********************************************************************
  1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
  1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
  True = reweight_scale ! Reweight to get scale variation using the
            ! rw_rscale and rw_fscale factors. Should be a list of
            ! booleans of equal length to dynamical_scale_choice to
            ! specify for which choice to include scale dependence.
  False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a
            ! list booleans of equal length to lhaid to specify for
            ! which PDF set to include the uncertainties.
#***********************************************************************
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#***********************************************************************
  False = store_rwgt_info ! Store info for reweighting in LHE file
#***********************************************************************
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#***********************************************************************
  0 = ickkw
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
# written in the LHE event file *
#***********************************************************************
  15.0 = bwcutoff
#***********************************************************************
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
  0.7 = jetradius ! The radius parameter for the jet algorithm
  10.0 = ptj ! Min jet transverse momentum
  -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
#***********************************************************************
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  0.0 = ptl ! Min lepton transverse momentum
  -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
  0.0 = drll ! Min distance between opposite sign lepton pairs
  0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
  0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
  30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
# all the other parameters are ignored. *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  20.0 = ptgmin ! Min photon transverse momentum
  -1.0 = etagamma ! Max photon abs(pseudo-rap)
  0.4 = r0gamma ! Radius of isolation code
  1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
  1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
  True = isoem ! isolate photons from EM energy (photons and leptons)
#***********************************************************************
# Cuts associated to MASSIVE particles identified by their PDG codes. *
# All cuts are applied to both particles and anti-particles, so use *
# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
# {6:100, 25:200} for multiple particles *
#***********************************************************************
  {} = pt_min_pdg ! Min pT for a massive particle
  {} = pt_max_pdg ! Max pT for a massive particle
  {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
#***********************************************************************
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
#***********************************************************************
  0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
#***********************************************************************

You should have a directory
/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/G5_full

(not sure about the name of the last directory)
which contains a log file with usefull information

Cheers,

Olivier

> On 28 Nov 2019, at 14:28, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi,
>
> Even without the b quark, I still failed to generate the events.
>
> The output is:
> INFO: Compiling directories...
> INFO: Compiling on 48 cores
> INFO: Compiling P0_uux_zx0_ztaptamz_no_a...
> INFO: Compiling P0_ddx_zx0_ztaptamz_no_a...
> INFO: Compiling P0_uxu_zx0_ztaptamz_no_a...
> Fontconfig warning: ignoring UTF-8: not a valid region tag
> INFO: Compiling P0_dxd_zx0_ztaptamz_no_a...
> INFO: Compiling P0_uux_zx0_wptaptamwm_no_a...
> INFO: Compiling P0_ddx_zx0_wptaptamwm_no_a...
> INFO: Compiling P0_uxu_zx0_wptaptamwm_no_a...
> INFO: Compiling P0_dxd_zx0_wptaptamwm_no_a...
> INFO: Compiling P0_dux_wmx0_wmtaptamz_no_a...
> INFO: Compiling P0_uxd_wmx0_wmtaptamz_no_a...
> INFO: Compiling P0_udx_wpx0_wptaptamz_no_a...
> INFO: Compiling P0_dxu_wpx0_wptaptamz_no_a...
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> INFO: P0_uxu_zx0_wptaptamwm_no_a done.
> INFO: P0_dxu_wpx0_wptaptamz_no_a done.
> INFO: P0_uxd_wmx0_wmtaptamz_no_a done.
> INFO: P0_ddx_zx0_ztaptamz_no_a done.
> INFO: P0_dux_wmx0_wmtaptamz_no_a done.
> INFO: P0_uux_zx0_wptaptamwm_no_a done.
> INFO: P0_uux_zx0_ztaptamz_no_a done.
> INFO: P0_udx_wpx0_wptaptamz_no_a done.
> INFO: P0_ddx_zx0_wptaptamwm_no_a done.
> INFO: P0_uxu_zx0_ztaptamz_no_a done.
> INFO: P0_dxd_zx0_ztaptamz_no_a done.
> INFO: P0_dxd_zx0_wptaptamwm_no_a done.
> INFO: Checking test output:
> INFO: P0_uux_zx0_ztaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_ddx_zx0_ztaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_uxu_zx0_ztaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_dxd_zx0_ztaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_uux_zx0_wptaptamwm_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_ddx_zx0_wptaptamwm_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_uxu_zx0_wptaptamwm_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_dxd_zx0_wptaptamwm_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_dux_wmx0_wmtaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_uxd_wmx0_wmtaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_udx_wpx0_wptaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: P0_dxu_wpx0_wptaptamz_no_a
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> WARNING: 0 points have been tried
> INFO: Starting run
> INFO: Using 48 cores
> INFO: Cleaning previous results
> INFO: Doing NLO matched to parton shower
> INFO: Setting up grids
> INFO: Idle: 4, Running: 48, Completed: 0 [ current time: 21h20 ]
> WARNING: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/ajob1 5 F 0 0 launch ends with non zero status: 1. Stop all computation
> INFO: Idle: 4, Running: 47, Completed: 1 [ 3.3s ]
> INFO: Idle: 4, Running: 46, Completed: 2 [ 3.3s ]
> INFO: Idle: 4, Running: 43, Completed: 5 [ 3.4s ]
> INFO: Idle: 4, Running: 42, Completed: 6 [ 3.4s ]
> INFO: Idle: 4, Running: 40, Completed: 8 [ 3.5s ]
> INFO: Idle: 4, Running: 39, Completed: 9 [ 3.5s ]
> INFO: Idle: 4, Running: 37, Completed: 11 [ 3.6s ]
> INFO: Idle: 4, Running: 36, Completed: 12 [ 3.6s ]
> INFO: Idle: 4, Running: 35, Completed: 13 [ 3.6s ]
> INFO: Idle: 4, Running: 34, Completed: 14 [ 3.7s ]
> INFO: Idle: 4, Running: 33, Completed: 15 [ 3.7s ]
> INFO: Idle: 4, Running: 32, Completed: 16 [ 3.8s ]
> INFO: Idle: 4, Running: 31, Completed: 17 [ 3.8s ]
> INFO: Idle: 4, Running: 30, Completed: 18 [ 3.8s ]
> INFO: Idle: 4, Running: 29, Completed: 19 [ 3.9s ]
> INFO: Idle: 4, Running: 28, Completed: 20 [ 3.9s ]
> INFO: Idle: 4, Running: 27, Completed: 21 [ 3.9s ]
> INFO: Idle: 4, Running: 26, Completed: 22 [ 4s ]
> INFO: Idle: 4, Running: 25, Completed: 23 [ 4s ]
> INFO: Idle: 4, Running: 24, Completed: 24 [ 4s ]
>
>
> And the debug log file is:
> launch
> Traceback (most recent call last):
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1514, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
> return func(arg, **opt)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1687, in do_launch
> evt_file = self.run(mode, options)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 1918, in run
> self.run_all_jobs(jobs_to_run,mint_step,fixed_order=False)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 2124, in run_all_jobs
> self.wait_for_complete(run_type)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/amcatnlo_run_interface.py", line 4692, in wait_for_complete
> self.cluster.wait(self.me_dir, update_status)
> File "/home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/bin/internal/cluster.py", line 844, in wait
> raise Exception, self.fail_msg
> Exception: program /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Total/NLO/HC_NLO_LHC_HH_2lep/SubProcesses/P0_uxd_wmx0_wmtaptamz_no_a/ajob1 5 F 0 0 launch ends with non zero status: 1. Stop all computation
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_local_path : None
> cluster_queue : None
> madanalysis_path : None
> lhapdf : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lhapdf6/bin/lhapdf-config
> cluster_size : 100
> cluster_memory : None
> pjfry : None
> cluster_status_update : (600, 30)
> cluster_time : None
> f2py_compiler : None
> ninja : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
> hepmc_path : None
> mg5amc_py8_interface_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/MG5aMC_PY8_interface
> pythia8_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/pythia8
> hwpp_path : None
> automatic_html_opening : False
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> mg5_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6
> td_path : None
> collier : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/HEPTools/lib
> delphes_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/Delphes
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
> fortran_compiler : None
> nb_core : 48
> OLP : MadLoop
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> notification_center : True
> run_mode : 2
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.7 2019-10-16 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set default_unset_couplings 99
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> 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 HC_NLO_X0_UFO
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> generate p p > z x0 > z l+ l- z /a [QCD]
> add process p p > z x0 > w+ l+ l- w- /a [QCD]
> add process p p > w- x0 > w- l+ l- z / a [QCD]
> add process p p > w+ x0 > w+ l+ l- z /a [QCD]
> output ./Total/NLO/HC_NLO_LHC_HH_2lep
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR FRBLOCK
> ###################################
> BLOCK FRBLOCK #
> 1 5.000000e+03 # lambda
> 2 1.000000e+00 # cosa
> 3 1.000000e+00 # ksm
> 4 0.000000e+00 # khtt
> 5 0.000000e+00 # katt
> 6 0.000000e+00 # khbb
> 7 0.000000e+00 # kabb
> 8 0.000000e+00 # khll
> 9 0.000000e+00 # kall
> 10 0.000000e+00 # khaa
> 11 0.000000e+00 # kaaa
> 12 0.000000e+00 # khza
> 13 0.000000e+00 # kaza
> 16 2.908000e-01 # khzz
> 17 0.000000e+00 # kazz
> 18 5.248000e-01 # khww
> 19 0.000000e+00 # kaww
> 20 0.000000e+00 # khda
> 21 -2.624000e-01 # khdz
> 22 -2.624000e-01 # khdwr (real part of khdw)
> 23 0.000000e+00 # khdwi (imaginary part of khdw)
> ###################################
> ## INFORMATION FOR LOOP
> ###################################
> BLOCK LOOP #
> 1 9.118800e+01 # mu_r
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.325070e+02 # aewm1
> 2 1.166390e-05 # gf
> 3 1.180000e-01 # as
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 5 4.700000e+00 # mb
> 6 1.730000e+02 # mt
> 15 1.777000e+00 # mta
> 23 9.118800e+01 # mz
> 25 1.250000e+02 # mx0
> 1 0.000000e+00 # d : 0.0
> 2 0.000000e+00 # u : 0.0
> 3 0.000000e+00 # s : 0.0
> 4 0.000000e+00 # c : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 0.0
> 14 0.000000e+00 # vm : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
> 82 0.000000e+00 # gh : 0.0
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> BLOCK QNUMBERS 82 # gh
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2s+1)
> 3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # particle/antiparticle distinction (0=own anti)
> #
> #*************************
> # Decay widths *
> #*************************
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.491500e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.441404e+00
> #
> # PDG Width
> DECAY 24 2.047600e+00
> #
> # PDG Width
> DECAY 25 8.485861e-04
> # BR NDA ID1 ID2 ...
> 1.519292e-01 3 -2 1 24 # 0.000128925007304
> 1.519292e-01 3 -4 3 24 # 0.000128925007304
> 1.504738e-01 3 -24 -1 2 # 0.000127689975094
> 1.504738e-01 3 -24 -3 4 # 0.000127689975094
> 5.024122e-02 3 -12 11 24 # 4.2634000939e-05
> 5.024122e-02 3 -14 13 24 # 4.2634000939e-05
> 5.007860e-02 3 -24 -11 12 # 4.24960038675e-05
> 5.007860e-02 3 -24 -13 14 # 4.24960038675e-05
> 4.968264e-02 3 -16 15 24 # 4.21599977153e-05
> 4.961429e-02 3 -24 -15 16 # 4.21019968554e-05
> 1.488829e-02 3 -1 1 23 # 1.26339959468e-05
> 1.488829e-02 3 -3 3 23 # 1.26339959468e-05
> 1.194340e-02 3 -5 5 23 # 1.01350032267e-05
> 1.176015e-02 3 -2 2 23 # 9.97949982391e-06
> 1.176015e-02 3 -4 4 23 # 9.97949982391e-06
> 6.692701e-03 3 -12 12 23 # 5.67933304006e-06
> 6.692701e-03 3 -14 14 23 # 5.67933304006e-06
> 6.692701e-03 3 -16 16 23 # 5.67933304006e-06
> 3.333663e-03 3 -11 11 23 # 2.82890008388e-06
> 3.333663e-03 3 -13 13 23 # 2.82890008388e-06
> 3.271677e-03 3 -15 15 23 # 2.77629962589e-06
> #
> # PDG Width
> DECAY 82 0.000000e+00
> #***********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat aMC@NLO *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a hash (#) are info or comments *
> # *
> # mind the format: value = variable ! comment *
> # *
> # Some of the values of variables can be list. These can either be *
> # comma or space separated. *
> # *
> # To display additional parameter, you can use the command: *
> # update to_full *
> #***********************************************************************
> #
> #*******************
> # Running parameters
> #*******************
> #
> #***********************************************************************
> # Tag name for the run (one word) *
> #***********************************************************************
> tag_1 = run_tag ! name of the run
> #***********************************************************************
> # Number of LHE events (and their normalization) and the required *
> # (relative) accuracy on the Xsec. *
> # These values are ignored for fixed order runs *
> #***********************************************************************
> 10000 = nevents ! Number of unweighted events requested
> -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
> -1 = nevt_job ! Max number of events per job in event generation.
> ! (-1= no split).
> #***********************************************************************
> # Normalize the weights of LHE events such that they sum or average to *
> # the total cross section *
> #***********************************************************************
> average = event_norm ! valid settings: average, sum, bias
> #***********************************************************************
> # Number of points per itegration channel (ignored for aMC@NLO runs) *
> #***********************************************************************
> 0.01 = req_acc_fo ! Required accuracy (-1=ignored, and use the
> ! number of points and iter. below)
> # These numbers are ignored except if req_acc_FO is equal to -1
> 5000 = npoints_fo_grid ! number of points to setup grids
> 4 = niters_fo_grid ! number of iter. to setup grids
> 10000 = npoints_fo ! number of points to compute Xsec
> 6 = niters_fo ! number of iter. to compute Xsec
> #***********************************************************************
> # Random number seed *
> #***********************************************************************
> 0 = iseed ! rnd seed (0=assigned automatically=default))
> #***********************************************************************
> # Collider type and energy *
> #***********************************************************************
> 1 = lpp1 ! beam 1 type (0 = no PDF)
> 1 = lpp2 ! beam 2 type (0 = no PDF)
> 6500.0 = ebeam1 ! beam 1 energy in GeV
> 6500.0 = ebeam2 ! beam 2 energy in GeV
> #***********************************************************************
> # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> #***********************************************************************
> nn23nlo = pdlabel ! PDF set
> 244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
> ! numbers for central PDF sets are allowed. Can be a list;
> ! PDF sets beyond the first are included via reweighting.
> #***********************************************************************
> # Include the NLO Monte Carlo subtr. terms for the following parton *
> # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
> # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
> #***********************************************************************
> PYTHIA8 = parton_shower
> 1.0 = shower_scale_factor ! multiply default shower starting
> ! scale by this factor
> #***********************************************************************
> # Renormalization and factorization scales *
> # (Default functional form for the non-fixed scales is the sum of *
> # the transverse masses divided by two of all final state particles *
> # and partons. This can be changed in SubProcesses/set_scales.f or via *
> # dynamical_scale_choice option) *
> #***********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.118 = mur_ref_fixed ! fixed ren reference scale
> 91.118 = muf_ref_fixed ! fixed fact reference scale
> 3 = dynamical_scale_choice ! Choose one (or more) of the predefined
> ! dynamical choices. Can be a list; scale choices beyond the
> ! first are included via reweighting
> 1.0 = mur_over_ref ! ratio of current muR over reference muR
> 1.0 = muf_over_ref ! ratio of current muF over reference muF
> #***********************************************************************
> # Reweight variables for scale dependence and PDF uncertainty *
> #***********************************************************************
> 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
> 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
> True = reweight_scale ! Reweight to get scale variation using the
> ! rw_rscale and rw_fscale factors. Should be a list of
> ! booleans of equal length to dynamical_scale_choice to
> ! specify for which choice to include scale dependence.
> False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a
> ! list booleans of equal length to lhaid to specify for
> ! which PDF set to include the uncertainties.
> #***********************************************************************
> # Store reweight information in the LHE file for off-line model- *
> # parameter reweighting at NLO+PS accuracy *
> #***********************************************************************
> False = store_rwgt_info ! Store info for reweighting in LHE file
> #***********************************************************************
> # ickkw parameter: *
> # 0: No merging *
> # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
> # level. After showering an MLM-type merging should be applied as *
> # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
> # 4: UNLOPS merging (with pythia8 only). No interface from within *
> # MG5_aMC available, but available in Pythia8. *
> # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
> #***********************************************************************
> 0 = ickkw
> #***********************************************************************
> #
> #***********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
> # written in the LHE event file *
> #***********************************************************************
> 15.0 = bwcutoff
> #***********************************************************************
> # Cuts on the jets. Jet clustering is performed by FastJet. *
> # - When matching to a parton shower, these generation cuts should be *
> # considerably softer than the analysis cuts. *
> # - More specific cuts can be specified in SubProcesses/cuts.f *
> #***********************************************************************
> 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 0.7 = jetradius ! The radius parameter for the jet algorithm
> 10.0 = ptj ! Min jet transverse momentum
> -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
> #***********************************************************************
> # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
> # More specific cuts can be specified in SubProcesses/cuts.f *
> #***********************************************************************
> 0.0 = ptl ! Min lepton transverse momentum
> -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 0.0 = drll ! Min distance between opposite sign lepton pairs
> 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
> #***********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
> # all the other parameters are ignored. *
> # More specific cuts can be specified in SubProcesses/cuts.f *
> #***********************************************************************
> 20.0 = ptgmin ! Min photon transverse momentum
> -1.0 = etagamma ! Max photon abs(pseudo-rap)
> 0.4 = r0gamma ! Radius of isolation code
> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> True = isoem ! isolate photons from EM energy (photons and leptons)
> #***********************************************************************
> # Cuts associated to MASSIVE particles identified by their PDG codes. *
> # All cuts are applied to both particles and anti-particles, so use *
> # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
> # {6:100, 25:200} for multiple particles *
> #***********************************************************************
> {} = pt_min_pdg ! Min pT for a massive particle
> {} = pt_max_pdg ! Max pT for a massive particle
> {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
> #***********************************************************************
> # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> #***********************************************************************
> 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> #***********************************************************************
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #12

There are two log file in directory GF5. One is log.txt:

===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 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
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

  External Params
  ---------------------------------

 MU_R = 91.188000000000002
 aEWM1 = 132.50700000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 173.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MB = 4.7000000000000002
 mdl_MZ = 91.188000000000002
 mdl_MTA = 1.7769999999999999
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 2.0476000000000001
 mdl_MX0 = 125.00000000000000
 mdl_WX0 = 8.4858609999999995E-004
 mdl_Lambda = 5000.0000000000000
 mdl_cosa = 1.0000000000000000
 mdl_kSM = 1.0000000000000000
 mdl_kHtt = 0.0000000000000000
 mdl_kAtt = 0.0000000000000000
 mdl_kHbb = 0.0000000000000000
 mdl_kAbb = 0.0000000000000000
 mdl_kHll = 0.0000000000000000
 mdl_kAll = 0.0000000000000000
 mdl_kHaa = 0.0000000000000000
 mdl_kAaa = 0.0000000000000000
 mdl_kHza = 0.0000000000000000
 mdl_kAza = 0.0000000000000000
 mdl_kHzz = 0.29080000000000000
 mdl_kAzz = 0.0000000000000000
 mdl_kHww = 0.52480000000000004
 mdl_kAww = 0.0000000000000000
 mdl_kHda = 0.0000000000000000
 mdl_kHdz = -0.26240000000000002
 mdl_kHdwR = -0.26240000000000002
 mdl_kHdwI = 0.0000000000000000
  Internal Params
  ---------------------------------

 mdl_CKM11 = 1.0000000000000000
 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM11 = 1.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_MZ__exp__2 = 8315.2513440000002
 mdl_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MX0__exp__2 = 15625.000000000000
 mdl_cosa__exp__2 = 1.0000000000000000
 mdl_sina = 0.0000000000000000
 mdl_kHdw = (-0.26240000000000002 , 0.0000000000000000 )
 mdl_nb__2__exp__0_75 = 1.6817928305074290
 mdl_Ncol__exp__2 = 9.0000000000000000
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29929.000000000000
 mdl_conjg__kHdw = (-0.26240000000000002 , -0.0000000000000000 )
 mdl_aEW = 7.5467711139788835E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.6872153846781555E-002
 mdl_ee = 0.30795376724436879
 mdl_MW__exp__2 = 6467.2159543705357
 mdl_sw2 = 0.22224648578577766
 mdl_cw = 0.88190334743339216
 mdl_sqrt__sw2 = 0.47143025548407230
 mdl_sw = 0.47143025548407230
 mdl_g1 = 0.34919219678733299
 mdl_gw = 0.65323293034757990
 mdl_vev = 246.21845810181637
 mdl_vev__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yb = 2.6995554250465490E-002
 mdl_yt = 0.99366614581500623
 mdl_ytau = 1.0206617000654717E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.5430507588273299E-002
 mdl_VectorZDown = -0.13030376310310560
 mdl_VectorAUp = 0.20530251149624587
 mdl_VectorADown = -0.10265125574812294
 mdl_VectorWmDxU = 0.23095271737156670
 mdl_AxialWmDxU = -0.23095271737156670
 mdl_VectorWpUxD = 0.23095271737156670
 mdl_AxialWpUxD = -0.23095271737156670
 mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
 mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.96667059156454072 , 0.0000000000000000 )
 mdl_Axial_tbGp = ( -1.0206617000654716 , -0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.96667059156454072 , 0.0000000000000000 )
 mdl_Axial_tbGm = ( -1.0206617000654716 , -0.0000000000000000 )
 mdl_ee__exp__2 = 9.4835522759998875E-002
 mdl_gAaa = 1.3008566310666950E-005
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_gAza = 4.7794971072590281E-006
 mdl_gHaa = 2.5475109025056106E-005
 mdl_gHza = 3.9182129211851395E-005
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_sw__exp__2 = 0.22224648578577769
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_G__exp__4 = 2.1987899468922913
 mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
 mdl_GWcft_UV_b_1EPS_ = -3.1300472141406080E-003
 mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
 mdl_bWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_G__exp__3 = 1.8056676068262196
 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_GWcft_UV_b_FIN_ = -1.8563438626678915E-002
 mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
 mdl_bWcft_UV_FIN_ = -0.13642100947319838
 mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
 mdl_gAgg = 7.6274879753643885E-005
 mdl_gHgg = -5.0849919835762592E-005
  Couplings of HC_NLO_X0_UFO
  ---------------------------------

    R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
    R2_UUZ_V5 = 0.00000E+00 0.68702E-03
      R2_dxuW = -0.00000E+00 -0.11566E-01
         GC_5 = 0.00000E+00 0.12177E+01
       R2_QQq = 0.00000E+00 0.12520E-01
     GC_3005a = -0.00000E+00 -0.00000E+00
    GC_3005h2 = 0.00000E+00 0.52480E-04
    GC_3005h3 = -0.00000E+00 -0.10496E-03
    GC_3005h4 = -0.00000E+00 0.52480E-04
     GC_3007a = -0.00000E+00 -0.00000E+00
    GC_3007h2 = -0.00000E+00 -0.58160E-04
    GC_3007h3 = 0.00000E+00 0.52480E-04
         GC_7 = 0.00000E+00 0.57609E+00
        GC_22 = 0.00000E+00 0.28804E+00
        GC_23 = -0.00000E+00 -0.27437E-01
        GC_24 = 0.00000E+00 0.82310E-01
        GC_47 = 0.00000E+00 0.46191E+00
    GC_3005h1 = 0.00000E+00 0.52532E+02
    GC_3007h1 = 0.00000E+00 0.67544E+02
     GC_3008a = 0.00000E+00 0.00000E+00
     GC_3008h = -0.00000E+00 -0.00000E+00

 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23nlo
 alpha_s(Mz)= 0.1190 running at 2 loops.
 alpha_s(Mz)= 0.1190 running at 2 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.9999999999999999E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 5
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the all of this channel
 Normal integration (Sfunction != 1)
 about to integrate 13 -1 12 5
 imode is 0
channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
STOP 1
 ------- iteration 1
 Update # PS points (even): 1040 --> 1040
Using random seed offsets: 5 , 10 , 0
  with seed 35
 Ranmar initialization seeds 14391 9419
 Total number of FKS directories is 4
 FKS process map (sum= 3 ) :
           1 --> 2 : 1 3
           2 --> 2 : 2 4
 ================================
 process combination map (specified per FKS dir):
  1 map 1 2
  1 inv. map 1 2
  2 map 1 2
  2 inv. map 1 2
  3 map 1 2
  3 inv. map 1 2
  4 map 1 2
  4 inv. map 1 2
 ================================
tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
 bpower is 0.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
 muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
 muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
 QES, QES_reference: 0.198743D+03 0.198743D+03 1.00

 muR_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF1_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF2_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 QES_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state

 alpha_s= 0.10641796709416652
 ERROR: INTEGRAL APPEARS TO BE ZERO.
 TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
Time in seconds: 4

Another is log_MINT0.txt:

 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat
 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 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
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

  External Params
  ---------------------------------

 MU_R = 91.188000000000002
 aEWM1 = 132.50700000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 173.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MB = 4.7000000000000002
 mdl_MZ = 91.188000000000002
 mdl_MTA = 1.7769999999999999
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 2.0476000000000001
 mdl_MX0 = 125.00000000000000
 mdl_WX0 = 8.4858609999999995E-004
 mdl_Lambda = 5000.0000000000000
 mdl_cosa = 1.0000000000000000
 mdl_kSM = 1.0000000000000000
 mdl_kHtt = 0.0000000000000000
 mdl_kAtt = 0.0000000000000000
 mdl_kHbb = 0.0000000000000000
 mdl_kAbb = 0.0000000000000000
 mdl_kHll = 0.0000000000000000
 mdl_kAll = 0.0000000000000000
 mdl_kHaa = 0.0000000000000000
 mdl_kAaa = 0.0000000000000000
 mdl_kHza = 0.0000000000000000
 mdl_kAza = 0.0000000000000000
 mdl_kHzz = 0.29080000000000000
 mdl_kAzz = 0.0000000000000000
 mdl_kHww = 0.52480000000000004
 mdl_kAww = 0.0000000000000000
 mdl_kHda = 0.0000000000000000
 mdl_kHdz = -0.26240000000000002
 mdl_kHdwR = -0.26240000000000002
 mdl_kHdwI = 0.0000000000000000
  Internal Params
  ---------------------------------

 mdl_CKM11 = 1.0000000000000000
 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM11 = 1.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_MZ__exp__2 = 8315.2513440000002
 mdl_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MX0__exp__2 = 15625.000000000000
 mdl_cosa__exp__2 = 1.0000000000000000
 mdl_sina = 0.0000000000000000
 mdl_kHdw = (-0.26240000000000002 , 0.0000000000000000 )
 mdl_nb__2__exp__0_75 = 1.6817928305074290
 mdl_Ncol__exp__2 = 9.0000000000000000
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29929.000000000000
 mdl_conjg__kHdw = (-0.26240000000000002 , -0.0000000000000000 )
 mdl_aEW = 7.5467711139788835E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.6872153846781555E-002
 mdl_ee = 0.30795376724436879
 mdl_MW__exp__2 = 6467.2159543705357
 mdl_sw2 = 0.22224648578577766
 mdl_cw = 0.88190334743339216
 mdl_sqrt__sw2 = 0.47143025548407230
 mdl_sw = 0.47143025548407230
 mdl_g1 = 0.34919219678733299
 mdl_gw = 0.65323293034757990
 mdl_vev = 246.21845810181637
 mdl_vev__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yb = 2.6995554250465490E-002
 mdl_yt = 0.99366614581500623
 mdl_ytau = 1.0206617000654717E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.5430507588273299E-002
 mdl_VectorZDown = -0.13030376310310560
 mdl_VectorAUp = 0.20530251149624587
 mdl_VectorADown = -0.10265125574812294
 mdl_VectorWmDxU = 0.23095271737156670
 mdl_AxialWmDxU = -0.23095271737156670
 mdl_VectorWpUxD = 0.23095271737156670
 mdl_AxialWpUxD = -0.23095271737156670
 mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
 mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.96667059156454072 , 0.0000000000000000 )
 mdl_Axial_tbGp = ( -1.0206617000654716 , -0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.96667059156454072 , 0.0000000000000000 )
 mdl_Axial_tbGm = ( -1.0206617000654716 , -0.0000000000000000 )
 mdl_ee__exp__2 = 9.4835522759998875E-002
 mdl_gAaa = 1.3008566310666950E-005
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_gAza = 4.7794971072590281E-006
 mdl_gHaa = 2.5475109025056106E-005
 mdl_gHza = 3.9182129211851395E-005
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_sw__exp__2 = 0.22224648578577769
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_G__exp__4 = 2.1987899468922913
 mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
 mdl_GWcft_UV_b_1EPS_ = -3.1300472141406080E-003
 mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
 mdl_bWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
 mdl_G__exp__3 = 1.8056676068262196
 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_GWcft_UV_b_FIN_ = -1.8563438626678915E-002
 mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
 mdl_bWcft_UV_FIN_ = -0.13642100947319838
 mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
 mdl_gAgg = 7.6274879753643885E-005
 mdl_gHgg = -5.0849919835762592E-005
  Couplings of HC_NLO_X0_UFO
  ---------------------------------

    R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
    R2_UUZ_V5 = 0.00000E+00 0.68702E-03
      R2_dxuW = -0.00000E+00 -0.11566E-01
         GC_5 = 0.00000E+00 0.12177E+01
       R2_QQq = 0.00000E+00 0.12520E-01
     GC_3005a = -0.00000E+00 -0.00000E+00
    GC_3005h2 = 0.00000E+00 0.52480E-04
    GC_3005h3 = -0.00000E+00 -0.10496E-03
    GC_3005h4 = -0.00000E+00 0.52480E-04
     GC_3007a = -0.00000E+00 -0.00000E+00
    GC_3007h2 = -0.00000E+00 -0.58160E-04
    GC_3007h3 = 0.00000E+00 0.52480E-04
         GC_7 = 0.00000E+00 0.57609E+00
        GC_22 = 0.00000E+00 0.28804E+00
        GC_23 = -0.00000E+00 -0.27437E-01
        GC_24 = 0.00000E+00 0.82310E-01
        GC_47 = 0.00000E+00 0.46191E+00
    GC_3005h1 = 0.00000E+00 0.52532E+02
    GC_3007h1 = 0.00000E+00 0.67544E+02
     GC_3008a = 0.00000E+00 0.00000E+00
     GC_3008h = -0.00000E+00 -0.00000E+00

 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23nlo
 alpha_s(Mz)= 0.1190 running at 2 loops.
 alpha_s(Mz)= 0.1190 running at 2 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 Diagram information for clustering has been set-up for nFKSprocess 1
 Diagram information for clustering has been set-up for nFKSprocess 2
 Diagram information for clustering has been set-up for nFKSprocess 3
 Diagram information for clustering has been set-up for nFKSprocess 4
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.9999999999999999E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 5
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the all of this channel
 Normal integration (Sfunction != 1)
 about to integrate 13 -1 12 5
 imode is 0
channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
STOP 1
 ------- iteration 1
 Update # PS points (even): 1040 --> 1040
Using random seed offsets: 5 , 10 , 0
  with seed 35
 Ranmar initialization seeds 14391 9419
 Total number of FKS directories is 4
 FKS process map (sum= 3 ) :
           1 --> 2 : 1 3
           2 --> 2 : 2 4
 ================================
 process combination map (specified per FKS dir):
  1 map 1 2
  1 inv. map 1 2
  2 map 1 2
  2 inv. map 1 2
  3 map 1 2
  3 inv. map 1 2
  4 map 1 2
  4 inv. map 1 2
 ================================
tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
 bpower is 0.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
 muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
 muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
 QES, QES_reference: 0.198743D+03 0.198743D+03 1.00

 muR_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF1_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF2_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 QES_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state

 alpha_s= 0.10641796709416652
 ERROR: INTEGRAL APPEARS TO BE ZERO.
 TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
Time in seconds: 4

Hi,

Looks like that you have zero coupling in your model:
> GC_3008a = 0.00000E+00 0.00000E+00
> GC_3008h = -0.00000E+00 -0.00000E+00

@NLO the code is crashing as soon as one integral provides 0 result.
In this case this is likely due to the face that such coupling are zero.
I would advise to use a restriction card to remove such interaction from your model
and then retry.

Cheers,

Olivier

PS: https://answers.launchpad.net/mg5amcnlo/+faq/2312
> On 28 Nov 2019, at 16:02, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/685744
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> There are two log file in directory GF5. One is log.txt:
>
> ===============================================================
> INFO: MadFKS read these parameters from FKS_params.dat
> ===============================================================
>> IRPoleCheckThreshold = 1.0000000000000001E-005
>> PrecisionVirtualAtRunTime = 1.0000000000000000E-003
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000001E-003
> ===============================================================
> 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
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> MU_R = 91.188000000000002
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MZ = 91.188000000000002
> mdl_MTA = 1.7769999999999999
> mdl_WT = 1.4915000000000000
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_MX0 = 125.00000000000000
> mdl_WX0 = 8.4858609999999995E-004
> mdl_Lambda = 5000.0000000000000
> mdl_cosa = 1.0000000000000000
> mdl_kSM = 1.0000000000000000
> mdl_kHtt = 0.0000000000000000
> mdl_kAtt = 0.0000000000000000
> mdl_kHbb = 0.0000000000000000
> mdl_kAbb = 0.0000000000000000
> mdl_kHll = 0.0000000000000000
> mdl_kAll = 0.0000000000000000
> mdl_kHaa = 0.0000000000000000
> mdl_kAaa = 0.0000000000000000
> mdl_kHza = 0.0000000000000000
> mdl_kAza = 0.0000000000000000
> mdl_kHzz = 0.29080000000000000
> mdl_kAzz = 0.0000000000000000
> mdl_kHww = 0.52480000000000004
> mdl_kAww = 0.0000000000000000
> mdl_kHda = 0.0000000000000000
> mdl_kHdz = -0.26240000000000002
> mdl_kHdwR = -0.26240000000000002
> mdl_kHdwI = 0.0000000000000000
> Internal Params
> ---------------------------------
>
> mdl_CKM11 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM11 = 1.0000000000000000
> mdl_lhv = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM33 = 1.0000000000000000
> mdl_Ncol = 3.0000000000000000
> mdl_CA = 3.0000000000000000
> mdl_TF = 0.50000000000000000
> mdl_CF = 1.3333333333333333
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MX0__exp__2 = 15625.000000000000
> mdl_cosa__exp__2 = 1.0000000000000000
> mdl_sina = 0.0000000000000000
> mdl_kHdw = (-0.26240000000000002 , 0.0000000000000000 )
> mdl_nb__2__exp__0_75 = 1.6817928305074290
> mdl_Ncol__exp__2 = 9.0000000000000000
> mdl_MB__exp__2 = 22.090000000000003
> mdl_MT__exp__2 = 29929.000000000000
> mdl_conjg__kHdw = (-0.26240000000000002 , -0.0000000000000000 )
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_AxialZUp = -0.18517701861793787
> mdl_AxialZDown = 0.18517701861793787
> mdl_VectorZUp = 7.5430507588273299E-002
> mdl_VectorZDown = -0.13030376310310560
> mdl_VectorAUp = 0.20530251149624587
> mdl_VectorADown = -0.10265125574812294
> mdl_VectorWmDxU = 0.23095271737156670
> mdl_AxialWmDxU = -0.23095271737156670
> mdl_VectorWpUxD = 0.23095271737156670
> mdl_AxialWpUxD = -0.23095271737156670
> mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_Vector_tbGp = (-0.96667059156454072 , 0.0000000000000000 )
> mdl_Axial_tbGp = ( -1.0206617000654716 , -0.0000000000000000 )
> mdl_Vector_tbGm = ( 0.96667059156454072 , 0.0000000000000000 )
> mdl_Axial_tbGm = ( -1.0206617000654716 , -0.0000000000000000 )
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_gAaa = 1.3008566310666950E-005
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_gAza = 4.7794971072590281E-006
> mdl_gHaa = 2.5475109025056106E-005
> mdl_gHza = 3.9182129211851395E-005
> mdl_gw__exp__2 = 0.42671326129048615
> mdl_sw__exp__2 = 0.22224648578577769
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> mdl_G__exp__4 = 2.1987899468922913
> mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
> mdl_GWcft_UV_b_1EPS_ = -3.1300472141406080E-003
> mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
> mdl_bWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_G__exp__3 = 1.8056676068262196
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_GWcft_UV_b_FIN_ = -1.8563438626678915E-002
> mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
> mdl_bWcft_UV_FIN_ = -0.13642100947319838
> mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
> mdl_gAgg = 7.6274879753643885E-005
> mdl_gHgg = -5.0849919835762592E-005
> Couplings of HC_NLO_X0_UFO
> ---------------------------------
>
> R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
> R2_UUZ_V5 = 0.00000E+00 0.68702E-03
> R2_dxuW = -0.00000E+00 -0.11566E-01
> GC_5 = 0.00000E+00 0.12177E+01
> R2_QQq = 0.00000E+00 0.12520E-01
> GC_3005a = -0.00000E+00 -0.00000E+00
> GC_3005h2 = 0.00000E+00 0.52480E-04
> GC_3005h3 = -0.00000E+00 -0.10496E-03
> GC_3005h4 = -0.00000E+00 0.52480E-04
> GC_3007a = -0.00000E+00 -0.00000E+00
> GC_3007h2 = -0.00000E+00 -0.58160E-04
> GC_3007h3 = 0.00000E+00 0.52480E-04
> GC_7 = 0.00000E+00 0.57609E+00
> GC_22 = 0.00000E+00 0.28804E+00
> GC_23 = -0.00000E+00 -0.27437E-01
> GC_24 = 0.00000E+00 0.82310E-01
> GC_47 = 0.00000E+00 0.46191E+00
> GC_3005h1 = 0.00000E+00 0.52532E+02
> GC_3007h1 = 0.00000E+00 0.67544E+02
> GC_3008a = 0.00000E+00 0.00000E+00
> GC_3008h = -0.00000E+00 -0.00000E+00
>
> Collider parameters:
> --------------------
>
> Running at P P machine @ 13000.000000000000 GeV
> PDF set = nn23nlo
> alpha_s(Mz)= 0.1190 running at 2 loops.
> alpha_s(Mz)= 0.1190 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> Diagram information for clustering has been set-up for nFKSprocess 1
> Diagram information for clustering has been set-up for nFKSprocess 2
> Diagram information for clustering has been set-up for nFKSprocess 3
> Diagram information for clustering has been set-up for nFKSprocess 4
> getting user params
> Enter number of events and iterations:
> Number of events and iterations -1 12
> Enter desired fractional accuracy:
> Desired fractional accuracy: 2.9999999999999999E-002
> Enter alpha, beta for G_soft
> Enter alpha<0 to set G_soft=1 (no ME soft)
> for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
> Enter alpha, beta for G_azi
> Enter alpha>0 to set G_azi=0 (no azi corr)
> for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
> Doing the S and H events together
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Do MC over helicities for the virtuals
> Enter Configuration Number:
> Running Configuration Number: 5
> Enter running mode for MINT:
> 0 to set-up grids, 1 to integrate, 2 to generate events
> MINT running mode: 0
> Set the three folding parameters for MINT
> xi_i, phi_i, y_ij
> 1 1 1
> 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
> Enter 'born0' or 'virt0' to perform
> a pure n-body integration (no S functions)
> doing the all of this channel
> Normal integration (Sfunction != 1)
> about to integrate 13 -1 12 5
> imode is 0
> channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
> STOP 1
> ------- iteration 1
> Update # PS points (even): 1040 --> 1040
> Using random seed offsets: 5 , 10 , 0
> with seed 35
> Ranmar initialization seeds 14391 9419
> Total number of FKS directories is 4
> FKS process map (sum= 3 ) :
> 1 --> 2 : 1 3
> 2 --> 2 : 2 4
> ================================
> process combination map (specified per FKS dir):
> 1 map 1 2
> 1 inv. map 1 2
> 2 map 1 2
> 2 inv. map 1 2
> 3 map 1 2
> 3 inv. map 1 2
> 4 map 1 2
> 4 inv. map 1 2
> ================================
> tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
> bpower is 0.0000000000000000
> Scale values (may change event by event):
> muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
> muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
> muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
> QES, QES_reference: 0.198743D+03 0.198743D+03 1.00
>
> muR_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF1_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF2_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> QES_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
>
> alpha_s= 0.10641796709416652
> ERROR: INTEGRAL APPEARS TO BE ZERO.
> TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
> Time in seconds: 4
>
>
>
>
>
>
>
>
> Another is log_MINT0.txt:
>
> ===============================================================
> INFO: MadFKS read these parameters from FKS_params.dat
> ===============================================================
>> IRPoleCheckThreshold = 1.0000000000000001E-005
>> PrecisionVirtualAtRunTime = 1.0000000000000000E-003
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000001E-003
> ===============================================================
> 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
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> MU_R = 91.188000000000002
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MZ = 91.188000000000002
> mdl_MTA = 1.7769999999999999
> mdl_WT = 1.4915000000000000
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_MX0 = 125.00000000000000
> mdl_WX0 = 8.4858609999999995E-004
> mdl_Lambda = 5000.0000000000000
> mdl_cosa = 1.0000000000000000
> mdl_kSM = 1.0000000000000000
> mdl_kHtt = 0.0000000000000000
> mdl_kAtt = 0.0000000000000000
> mdl_kHbb = 0.0000000000000000
> mdl_kAbb = 0.0000000000000000
> mdl_kHll = 0.0000000000000000
> mdl_kAll = 0.0000000000000000
> mdl_kHaa = 0.0000000000000000
> mdl_kAaa = 0.0000000000000000
> mdl_kHza = 0.0000000000000000
> mdl_kAza = 0.0000000000000000
> mdl_kHzz = 0.29080000000000000
> mdl_kAzz = 0.0000000000000000
> mdl_kHww = 0.52480000000000004
> mdl_kAww = 0.0000000000000000
> mdl_kHda = 0.0000000000000000
> mdl_kHdz = -0.26240000000000002
> mdl_kHdwR = -0.26240000000000002
> mdl_kHdwI = 0.0000000000000000
> Internal Params
> ---------------------------------
>
> mdl_CKM11 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM11 = 1.0000000000000000
> mdl_lhv = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM33 = 1.0000000000000000
> mdl_Ncol = 3.0000000000000000
> mdl_CA = 3.0000000000000000
> mdl_TF = 0.50000000000000000
> mdl_CF = 1.3333333333333333
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MX0__exp__2 = 15625.000000000000
> mdl_cosa__exp__2 = 1.0000000000000000
> mdl_sina = 0.0000000000000000
> mdl_kHdw = (-0.26240000000000002 , 0.0000000000000000 )
> mdl_nb__2__exp__0_75 = 1.6817928305074290
> mdl_Ncol__exp__2 = 9.0000000000000000
> mdl_MB__exp__2 = 22.090000000000003
> mdl_MT__exp__2 = 29929.000000000000
> mdl_conjg__kHdw = (-0.26240000000000002 , -0.0000000000000000 )
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_AxialZUp = -0.18517701861793787
> mdl_AxialZDown = 0.18517701861793787
> mdl_VectorZUp = 7.5430507588273299E-002
> mdl_VectorZDown = -0.13030376310310560
> mdl_VectorAUp = 0.20530251149624587
> mdl_VectorADown = -0.10265125574812294
> mdl_VectorWmDxU = 0.23095271737156670
> mdl_AxialWmDxU = -0.23095271737156670
> mdl_VectorWpUxD = 0.23095271737156670
> mdl_AxialWpUxD = -0.23095271737156670
> mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_Vector_tbGp = (-0.96667059156454072 , 0.0000000000000000 )
> mdl_Axial_tbGp = ( -1.0206617000654716 , -0.0000000000000000 )
> mdl_Vector_tbGm = ( 0.96667059156454072 , 0.0000000000000000 )
> mdl_Axial_tbGm = ( -1.0206617000654716 , -0.0000000000000000 )
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_gAaa = 1.3008566310666950E-005
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_gAza = 4.7794971072590281E-006
> mdl_gHaa = 2.5475109025056106E-005
> mdl_gHza = 3.9182129211851395E-005
> mdl_gw__exp__2 = 0.42671326129048615
> mdl_sw__exp__2 = 0.22224648578577769
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> mdl_G__exp__4 = 2.1987899468922913
> mdl_R2MixedFactor_FIN_ = -2.5040377713124864E-002
> mdl_GWcft_UV_b_1EPS_ = -3.1300472141406080E-003
> mdl_GWcft_UV_t_1EPS_ = -3.1300472141406080E-003
> mdl_bWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_tWcft_UV_1EPS_ = -1.8780283284843650E-002
> mdl_G__exp__3 = 1.8056676068262196
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_GWcft_UV_b_FIN_ = -1.8563438626678915E-002
> mdl_GWcft_UV_t_FIN_ = 4.0087659331150384E-003
> mdl_bWcft_UV_FIN_ = -0.13642100947319838
> mdl_tWcft_UV_FIN_ = -9.8778211443463623E-004
> mdl_gAgg = 7.6274879753643885E-005
> mdl_gHgg = -5.0849919835762592E-005
> Couplings of HC_NLO_X0_UFO
> ---------------------------------
>
> R2_UUZ_V2 = -0.00000E+00 -0.72127E-02
> R2_UUZ_V5 = 0.00000E+00 0.68702E-03
> R2_dxuW = -0.00000E+00 -0.11566E-01
> GC_5 = 0.00000E+00 0.12177E+01
> R2_QQq = 0.00000E+00 0.12520E-01
> GC_3005a = -0.00000E+00 -0.00000E+00
> GC_3005h2 = 0.00000E+00 0.52480E-04
> GC_3005h3 = -0.00000E+00 -0.10496E-03
> GC_3005h4 = -0.00000E+00 0.52480E-04
> GC_3007a = -0.00000E+00 -0.00000E+00
> GC_3007h2 = -0.00000E+00 -0.58160E-04
> GC_3007h3 = 0.00000E+00 0.52480E-04
> GC_7 = 0.00000E+00 0.57609E+00
> GC_22 = 0.00000E+00 0.28804E+00
> GC_23 = -0.00000E+00 -0.27437E-01
> GC_24 = 0.00000E+00 0.82310E-01
> GC_47 = 0.00000E+00 0.46191E+00
> GC_3005h1 = 0.00000E+00 0.52532E+02
> GC_3007h1 = 0.00000E+00 0.67544E+02
> GC_3008a = 0.00000E+00 0.00000E+00
> GC_3008h = -0.00000E+00 -0.00000E+00
>
> Collider parameters:
> --------------------
>
> Running at P P machine @ 13000.000000000000 GeV
> PDF set = nn23nlo
> alpha_s(Mz)= 0.1190 running at 2 loops.
> alpha_s(Mz)= 0.1190 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> Diagram information for clustering has been set-up for nFKSprocess 1
> Diagram information for clustering has been set-up for nFKSprocess 2
> Diagram information for clustering has been set-up for nFKSprocess 3
> Diagram information for clustering has been set-up for nFKSprocess 4
> getting user params
> Enter number of events and iterations:
> Number of events and iterations -1 12
> Enter desired fractional accuracy:
> Desired fractional accuracy: 2.9999999999999999E-002
> Enter alpha, beta for G_soft
> Enter alpha<0 to set G_soft=1 (no ME soft)
> for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
> Enter alpha, beta for G_azi
> Enter alpha>0 to set G_azi=0 (no azi corr)
> for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
> Doing the S and H events together
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Do MC over helicities for the virtuals
> Enter Configuration Number:
> Running Configuration Number: 5
> Enter running mode for MINT:
> 0 to set-up grids, 1 to integrate, 2 to generate events
> MINT running mode: 0
> Set the three folding parameters for MINT
> xi_i, phi_i, y_ij
> 1 1 1
> 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
> Enter 'born0' or 'virt0' to perform
> a pure n-body integration (no S functions)
> doing the all of this channel
> Normal integration (Sfunction != 1)
> about to integrate 13 -1 12 5
> imode is 0
> channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
> STOP 1
> ------- iteration 1
> Update # PS points (even): 1040 --> 1040
> Using random seed offsets: 5 , 10 , 0
> with seed 35
> Ranmar initialization seeds 14391 9419
> Total number of FKS directories is 4
> FKS process map (sum= 3 ) :
> 1 --> 2 : 1 3
> 2 --> 2 : 2 4
> ================================
> process combination map (specified per FKS dir):
> 1 map 1 2
> 1 inv. map 1 2
> 2 map 1 2
> 2 inv. map 1 2
> 3 map 1 2
> 3 inv. map 1 2
> 4 map 1 2
> 4 inv. map 1 2
> ================================
> tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
> bpower is 0.0000000000000000
> Scale values (may change event by event):
> muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
> muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
> muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
> QES, QES_reference: 0.198743D+03 0.198743D+03 1.00
>
> muR_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF1_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF2_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> QES_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
>
> alpha_s= 0.10641796709416652
> ERROR: INTEGRAL APPEARS TO BE ZERO.
> TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
> Time in seconds: 4
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Yue Xu (yuexu) said : #14

Hi Olivier,

I added the restriction card and the couplings not interested are removed, then it works well.
Thanks a lot !

Cheers,
Yue