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

Asked by Yue Xu on 2019-11-11

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
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
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
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:
Assignee:
No assignee Edit question
Solved by:
Yue Xu
Solved:
2019-12-03
Last query:
2019-12-03
2019-11-28
 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-11: #1

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

 Yue Xu (yuexu) said on 2019-11-11: #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

 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-16: #3

Hi,

I have run the following two scripts:
generate p p > z j
add process p p > w+ j
output
launch
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z j
output
launch
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
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
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z z
output
launch
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
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
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
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z j [QCD]
output
launch
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
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:
>
>
> 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
>
> --

 Yue Xu (yuexu) said on 2019-11-26: #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

 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-26: #5

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

 Yue Xu (yuexu) said on 2019-11-28: #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

 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-28: #7

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

 Yue Xu (yuexu) said on 2019-11-28: #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
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
exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
fortran_compiler : None
nb_core : 48
auto_update : 7
cluster_nb_retry : 1
eps_viewer : None
syscalc_path : None
golem : None
cpp_compiler : None
run_mode : 2
#************************************************************
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.6 2018-06-28 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* *
#************************************************************
#* *
#* 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
#***********************************************************************
# *
# 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)
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)
#***********************************************************************

 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-28: #9

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:
>
>
> 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
> 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
> exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
> fortran_compiler : None
> nb_core : 48
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> run_mode : 2
> #************************************************************
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.6 2018-06-28 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* *
> #************************************************************
> #* *
> #* 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
> #***********************************************************************
> # *
> # 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)
> #***********************************************************************
>
> --

 Yue Xu (yuexu) said on 2019-11-28: #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
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
exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
fortran_compiler : None
nb_core : 48
auto_update : 7
cluster_nb_retry : 1
eps_viewer : None
syscalc_path : None
golem : None
cpp_compiler : None
run_mode : 2
#************************************************************
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.7 2019-10-16 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* *
#************************************************************
#* *
#* 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
#***********************************************************************
# *
# 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)
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)
#***********************************************************************

 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-28: #11

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:
>
>
> 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
> 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
> exrootanalysis_path : /home/storage/Users/xuyue/HeavyHiggs/plot_heavy/MG5_aMC_v2_6_6/ExRootAnalysis
> fortran_compiler : None
> nb_core : 48
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> run_mode : 2
> #************************************************************
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.7 2019-10-16 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* *
> #************************************************************
> #* *
> #* 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
> #***********************************************************************
> # *
> # 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)
> #***********************************************************************
>
> --

 Yue Xu (yuexu) said on 2019-11-28: #12

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

===============================================================
===============================================================
> 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
*****************************************************
* -------------------------------- *
* -------------------------------- *
* *
* 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_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:

===============================================================
===============================================================
> 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
*****************************************************
* -------------------------------- *
* -------------------------------- *
* *
* 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_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

 Olivier Mattelaer (olivier-mattelaer) said on 2019-11-28: #13

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

> On 28 Nov 2019, at 16:02, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
>
>
> Yue Xu is still having a problem:
> There are two log file in directory GF5. One is log.txt:
>
> ===============================================================
> ===============================================================
>> 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
> *****************************************************
> * -------------------------------- *
> * -------------------------------- *
> * *
> * 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_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:
>
> ===============================================================
> ===============================================================
>> 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
> *****************************************************
> * -------------------------------- *
> * -------------------------------- *
> * *
> * 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_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
>
> --

 Yue Xu (yuexu) said on 2019-12-03: #14

Hi Olivier,

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

Cheers,
Yue