How does MagGraph5 calculate cross-section when mad-spin is used with NLO ?
Hello MadGraph group,
I am generating events with a modified HC model at NLO and madspin is used to decay W and Z bosons, here is the process(labeled as A )(x0 is a heavy higgs):
import model HC_NLO_X0_UFO
generate p p > z x0 > z e+ e- z /a [QCD]
add process p p > z x0 > z mu+ mu- z /a [QCD]
add process p p > z x0 > w+ e+ e- w- /a [QCD]
add process p p > z x0 > w+ mu+ mu- w- /a [QCD]
add process p p > w- x0 > w- e+ e- z / a [QCD]
add process p p > w- x0 > w- mu+ mu- z / a [QCD]
add process p p > w+ x0 > w+ e+ e- z /a [QCD]
add process p p > w+ x0 > w+ mu+ mu- z /a [QCD]
and the set of madspin_card is:
decay w+ > j j
decay w- > j j
decay z > j j
launch
the model can be fount at:
https:/
the process is here:
https:/
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_
madspin:
decay z > j j
launch
I got the cs before decay is 5.239e-03 and the branch ratio is 0.294154, so the final cross-section is 0.00154107. This process is here:
https:/
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_
madspin:
decay w+ > j j
decay w- > j j
launch
The branch ratio is 0.444518 and the final cross-section is 0.00283000. The process is here:
https:/
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_
madspin:
decay w+ > j j
decay w- > j j
decay z > j j
launch
The branch ratio is 0.361603 and the final cross-section is 0.005093. The process is here:
https:/
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
- Assignee:
- No assignee Edit question
- Solved by:
- Yue Xu
- Solved:
- 2019-12-03
- Last query:
- 2019-12-03
- Last reply:
- 2019-11-28
This should be related to statistical fluctuation no?
Having different BR means that we have to throw away events and this will be a new source of statistical error.
Cheers,
Olivier
> On 11 Nov 2019, at 09:23, Yue Xu <email address hidden> wrote:
>
> New question #685744 on MadGraph5_aMC@NLO:
> https:/
>
> Hello MadGraph group,
>
> I am generating events with a modified HC model at NLO and madspin is used to decay W and Z bosons, here is the process(labeled as A )(x0 is a heavy higgs):
> import model HC_NLO_X0_UFO
> generate p p > z x0 > z e+ e- z /a [QCD]
> add process p p > z x0 > z mu+ mu- z /a [QCD]
> add process p p > z x0 > w+ e+ e- w- /a [QCD]
> add process p p > z x0 > w+ mu+ mu- w- /a [QCD]
> add process p p > w- x0 > w- e+ e- z / a [QCD]
> add process p p > w- x0 > w- mu+ mu- z / a [QCD]
> add process p p > w+ x0 > w+ e+ e- z /a [QCD]
> add process p p > w+ x0 > w+ mu+ mu- z /a [QCD]
>
> and the set of madspin_card is:
> decay w+ > j j
> decay w- > j j
> decay z > j j
> launch
>
> the model can be fount at:
> https:/
> the process is here:
> https:/
> 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_
> madspin:
> decay z > j j
> launch
> I got the cs before decay is 5.239e-03 and the branch ratio is 0.294154, so the final cross-section is 0.00154107. This process is here:
> https:/
>
> 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_
> madspin:
> decay w+ > j j
> decay w- > j j
> launch
> The branch ratio is 0.444518 and the final cross-section is 0.00283000. The process is here:
> https:/
>
> 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_
> madspin:
> decay w+ > j j
> decay w- > j j
> decay z > j j
> launch
>
> The branch ratio is 0.361603 and the final cross-section is 0.005093. The process is here:
> https:/
>
> If I add the cross-section of process 1, 2 and 3 together, the total cross section is 0.00946407 which is much smaller than 0.01143.
> But the final state of process 1,2 and 3 is the same as A, they should have the same cross section.
> How does this happen?
> How can we get the correct cross section when processes with different branch ratio are included?
>
> Cheers,
> Yue
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #2 |
In case A, the error is 0.01143 ± 1.4e-05. But the cross section difference between case A and the case adding process1,2,3 together reach 20%, it's so large.
If I want a precision cross section, is it better to generate events separately?
Cheers,
Yue
Hi,
I have run the following two scripts:
generate p p > z j
add process p p > w+ j
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z j
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
#generate p p > z z /a
generate p p > w+ j
#add process p p > w- z / a
#add process p p > w+ z /a
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
and
generate p p > z z
add process p p > w+ w-
add process p p > w- z
add process p p > w+ z
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z z
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
#generate p p > z z /a
generate p p > w+ w-
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > w- z
add process p p > w+ z
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
In both case, I have the expected result after madspin. (variation at the level of 1 per 1000)
Also this might be an effect related to the negative number of events.
In the case where I do not have a unique BR, I need to re-calculate the cross-section
from the generated events but since negative events are present this can increase statistical error.
I have tested this on
generate p p > z j [QCD]
add process p p > w+ j [QCD]
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
generate p p > z j [QCD]
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
#generate p p > z z /a
generate p p > w+ j [QCD]
#add process p p > w- z / a
#add process p p > w+ z /a
output
launch
madspin=ON
decay z > j j
decay w+ > j j
decay w- > j j
and here I have a larger difference. But this expected since this is due to the amount of negative weights in the sample (in this case ~ 25%)
Cheers,
Olivier
> On 11 Nov 2019, at 14:57, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> In case A, the error is 0.01143 ± 1.4e-05. But the cross section difference between case A and the case adding process1,2,3 together reach 20%, it's so large.
> If I want a precision cross section, is it better to generate events separately?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #4 |
Hi Olivier,
We want to re-scale the LO cross-section to NLO cross-section.
But in the NLO case there is a large difference because of the negative weights.
Is it reliable to use such a cross-section? or it's better to calculate the cross-section separately?
Cheers,
Yue
To have a more accurate value, it is indeed better to not use this combine mode.
Cheers,
Olivier
> On 26 Nov 2019, at 03:43, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi Olivier,
>
> We want to re-scale the LO cross-section to NLO cross-section.
> But in the NLO case there is a large difference because of the negative weights.
>
> Is it reliable to use such a cross-section? or it's better to calculate
> the cross-section separately?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #6 |
Hi,
I'm considering the interference of the processes with same initial and final state.
For the process like:
generate p p > z x0 > l+ l- j j j j
add process p p > w+ x0 > l- l+ j j j j
the two sub-processes may have interference.
But for the process like:
generate p p > z x0 > w- l+ l- w+ /a [QCD]
add process p p > w+ x0 > z l+ l- z /a [QCD]
manspin=ON
decay z > j j
decay w->j j
decay w+>j j
Will the two sub-processes have interference ?
Cheers,
Yue
Hi,
None of the above syntax will have interference term. (even the one without using madspin)
Cheers,
Olivier
> On 28 Nov 2019, at 08:47, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi,
>
> I'm considering the interference of the processes with same initial and final state.
> For the process like:
> generate p p > z x0 > l+ l- j j j j
> add process p p > w+ x0 > l- l+ j j j j
>
> the two sub-processes may have interference.
> But for the process like:
> generate p p > z x0 > w- l+ l- w+ /a [QCD]
> add process p p > w+ x0 > z l+ l- z /a [QCD]
> manspin=ON
> decay z > j j
> decay w->j j
> decay w+>j j
>
> Will the two sub-processes have interference ?
>
> Cheers,
> Yue
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #8 |
Hi Olivier,
Thanks a lot!
Then, separating the processes can be a choice to get an accurate cross section value.
I met another problem when use this HC NLO model to generate events.
The process is (tau is added to the lepton definition, b quark is added to the p and j definition):
import model HC_NLO_X0_UFO
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
generate p p > z x0 > z l+ l- z /a [QCD]
add process p p > z x0 > w+ l+ l- w- /a [QCD]
add process p p > w- x0 > w- l+ l- z / a [QCD]
add process p p > w+ x0 > w+ l+ l- z /a [QCD]
I got the warning like "WARNING: program /home/storage/
I don't know how to make it work.
Could you do me a favour?
Thanks,
Yue
P.S.
When generating events with "./bin/
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/
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/
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/
INFO: Idle: 22, Running: 2, Completed: 114 [ 2m 25s ]
and the debug log is here:
launch
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
evt_file = self.run(mode, options)
File "/home/
self.
File "/home/
self.
File "/home/
self.
File "/home/
raise Exception, self.fail_msg
Exception: program /home/storage/
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
evt_file = self.run(mode, options)
File "/home/
self.
File "/home/
self.
File "/home/
self.
File "/home/
raise Exception, self.fail_msg
Exception: program /home/storage/
Value of current Options:
cluster_
mg5amc_
automatic_
exrootana
notificat
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.6 2018-06-28 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
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/
#######
## 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(
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/
#
#******
# Decay widths *
#******
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.491500e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.441404e+00
#
# PDG Width
DECAY 24 2.047600e+00
#
# PDG Width
DECAY 25 1.634610e-01
# BR NDA ID1 ID2 ...
7.066517e-01 2 -24 24 # 0.115509993534
2.933483e-01 2 23 23 # 0.0479510064663
#
# PDG Width
DECAY 82 0.000000e+00
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# 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
# 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=
#******
# 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
#******
# 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/
# dynamical_
#******
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_
! 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_
! 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://
# 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*
# 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://
#******
0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
#******
Hi,
> b quark is added to the p and j definition):
This does not make sense.
Either your computation is 4 flavor (with b massive)
or 5 flavour (and b massless)
MG5amC, will automatically adapt the definition of p/j according to the fact that your model defines the b quark as massive or massless.
This being said, such computation is not going to work in the 5 flavour scheme due to the possibility to have a top onshell. (i.e. you need. to use diagram substraction for that which is available as a plugin but I do not know how that one work)
Cheers,
Olivier
> On 28 Nov 2019, at 13:47, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi Olivier,
>
> Thanks a lot!
> Then, separating the processes can be a choice to get an accurate cross section value.
>
> I met another problem when use this HC NLO model to generate events.
> The process is (tau is added to the lepton definition, b quark is added to the p and j definition):
> import model HC_NLO_X0_UFO
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = g u c d s b u~ c~ d~ s~ b~
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> generate p p > z x0 > z l+ l- z /a [QCD]
> add process p p > z x0 > w+ l+ l- w- /a [QCD]
> add process p p > w- x0 > w- l+ l- z / a [QCD]
> add process p p > w+ x0 > w+ l+ l- z /a [QCD]
>
> I got the warning like "WARNING: program
> /home/storage/
> 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/
>
> 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/
> 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/
> 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/
> INFO: Idle: 22, Running: 2, Completed: 114 [ 2m 25s ]
>
>
> and the debug log is here:
> launch
> Traceback (most recent call last):
> File "/home/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> evt_file = self.run(mode, options)
> File "/home/
> self.run_
> File "/home/
> self.wait_
> File "/home/
> self.cluster.
> File "/home/
> raise Exception, self.fail_msg
> Exception: program /home/storage/
> Traceback (most recent call last):
> File "/home/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> evt_file = self.run(mode, options)
> File "/home/
> self.run_
> File "/home/
> self.wait_
> File "/home/
> self.cluster.
> File "/home/
> raise Exception, self.fail_msg
> Exception: program /home/storage/
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_local_path : None
> cluster_queue : None
> madanalysis_path : None
> lhapdf : /home/storage/
> cluster_size : 100
> cluster_memory : None
> pjfry : None
> cluster_
> cluster_time : None
> f2py_compiler : None
> ninja : /home/storage/
> hepmc_path : None
> mg5amc_
> pythia8_path : /home/storage/
> hwpp_path : None
> automatic_
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> mg5_path : /home/storage/
> td_path : None
> collier : /home/storage/
> delphes_path : /home/storage/
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : /home/storage/
> fortran_compiler : None
> nb_core : 48
> OLP : MadLoop
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> notification_center : True
> run_mode : 2
> #******
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.6 2018-06-28 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https:/
> #* *
> #******
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #******
> set default_
> set group_subprocesses Auto
> set ignore_
> set loop_optimized_
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_
> 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/
> #######
> ## 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(
> 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/
> #
> #******
> # Decay widths *
> #******
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.491500e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.441404e+00
> #
> # PDG Width
> DECAY 24 2.047600e+00
> #
> # PDG Width
> DECAY 25 1.634610e-01
> # BR NDA ID1 ID2 ...
> 7.066517e-01 2 -24 24 # 0.115509993534
> 2.933483e-01 2 23 23 # 0.0479510064663
> #
> # PDG Width
> DECAY 82 0.000000e+00
> #******
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat aMC@NLO *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/
> # *
> # 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=
> #******
> # 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/
> # dynamical_
> #******
> 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_
> ! 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_
> ! 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://
> # 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*
> # 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://
> #******
> 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> #******
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #10 |
Hi,
Even without the b quark, I still failed to generate the events.
The output is:
INFO: Compiling directories...
INFO: Compiling on 48 cores
INFO: Compiling P0_uux_
INFO: Compiling P0_ddx_
INFO: Compiling P0_uxu_
Fontconfig warning: ignoring UTF-8: not a valid region tag
INFO: Compiling P0_dxd_
INFO: Compiling P0_uux_
INFO: Compiling P0_ddx_
INFO: Compiling P0_uxu_
INFO: Compiling P0_dxd_
INFO: Compiling P0_dux_
INFO: Compiling P0_uxd_
INFO: Compiling P0_udx_
INFO: Compiling P0_dxu_
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_
INFO: P0_dxu_
INFO: P0_uxd_
INFO: P0_ddx_
INFO: P0_dux_
INFO: P0_uux_
INFO: P0_uux_
INFO: P0_udx_
INFO: P0_ddx_
INFO: P0_uxu_
INFO: P0_dxd_
INFO: P0_dxd_
INFO: Checking test output:
INFO: P0_uux_
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_
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_
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_
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_
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_
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_
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_
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_
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_
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_
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_
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/
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/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
evt_file = self.run(mode, options)
File "/home/
self.
File "/home/
self.
File "/home/
self.
File "/home/
raise Exception, self.fail_msg
Exception: program /home/storage/
Value of current Options:
cluster_
mg5amc_
automatic_
exrootana
notificat
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.7 2019-10-16 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
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/
#######
## 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(
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/
#
#******
# Decay widths *
#******
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.491500e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.441404e+00
#
# PDG Width
DECAY 24 2.047600e+00
#
# PDG Width
DECAY 25 8.485861e-04
# BR NDA ID1 ID2 ...
1.519292e-01 3 -2 1 24 # 0.000128925007304
1.519292e-01 3 -4 3 24 # 0.000128925007304
1.504738e-01 3 -24 -1 2 # 0.000127689975094
1.504738e-01 3 -24 -3 4 # 0.000127689975094
5.024122e-02 3 -12 11 24 # 4.2634000939e-05
5.024122e-02 3 -14 13 24 # 4.2634000939e-05
5.007860e-02 3 -24 -11 12 # 4.24960038675e-05
5.007860e-02 3 -24 -13 14 # 4.24960038675e-05
4.968264e-02 3 -16 15 24 # 4.21599977153e-05
4.961429e-02 3 -24 -15 16 # 4.21019968554e-05
1.488829e-02 3 -1 1 23 # 1.26339959468e-05
1.488829e-02 3 -3 3 23 # 1.26339959468e-05
1.194340e-02 3 -5 5 23 # 1.01350032267e-05
1.176015e-02 3 -2 2 23 # 9.97949982391e-06
1.176015e-02 3 -4 4 23 # 9.97949982391e-06
6.692701e-03 3 -12 12 23 # 5.67933304006e-06
6.692701e-03 3 -14 14 23 # 5.67933304006e-06
6.692701e-03 3 -16 16 23 # 5.67933304006e-06
3.333663e-03 3 -11 11 23 # 2.82890008388e-06
3.333663e-03 3 -13 13 23 # 2.82890008388e-06
3.271677e-03 3 -15 15 23 # 2.77629962589e-06
#
# PDG Width
DECAY 82 0.000000e+00
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# 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
# 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=
#******
# 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
#******
# 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/
# dynamical_
#******
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_
! 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_
! 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://
# 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*
# 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://
#******
0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
#******
You should have a directory
/home/storage/
(not sure about the name of the last directory)
which contains a log file with usefull information
Cheers,
Olivier
> On 28 Nov 2019, at 14:28, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> Hi,
>
> Even without the b quark, I still failed to generate the events.
>
> The output is:
> INFO: Compiling directories...
> INFO: Compiling on 48 cores
> INFO: Compiling P0_uux_
> INFO: Compiling P0_ddx_
> INFO: Compiling P0_uxu_
> Fontconfig warning: ignoring UTF-8: not a valid region tag
> INFO: Compiling P0_dxd_
> INFO: Compiling P0_uux_
> INFO: Compiling P0_ddx_
> INFO: Compiling P0_uxu_
> INFO: Compiling P0_dxd_
> INFO: Compiling P0_dux_
> INFO: Compiling P0_uxd_
> INFO: Compiling P0_udx_
> INFO: Compiling P0_dxu_
> 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_
> INFO: P0_dxu_
> INFO: P0_uxd_
> INFO: P0_ddx_
> INFO: P0_dux_
> INFO: P0_uux_
> INFO: P0_uux_
> INFO: P0_udx_
> INFO: P0_ddx_
> INFO: P0_uxu_
> INFO: P0_dxd_
> INFO: P0_dxd_
> INFO: Checking test output:
> INFO: P0_uux_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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_
> 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/
> 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/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> evt_file = self.run(mode, options)
> File "/home/
> self.run_
> File "/home/
> self.wait_
> File "/home/
> self.cluster.
> File "/home/
> raise Exception, self.fail_msg
> Exception: program /home/storage/
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_local_path : None
> cluster_queue : None
> madanalysis_path : None
> lhapdf : /home/storage/
> cluster_size : 100
> cluster_memory : None
> pjfry : None
> cluster_
> cluster_time : None
> f2py_compiler : None
> ninja : /home/storage/
> hepmc_path : None
> mg5amc_
> pythia8_path : /home/storage/
> hwpp_path : None
> automatic_
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> mg5_path : /home/storage/
> td_path : None
> collier : /home/storage/
> delphes_path : /home/storage/
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : /home/storage/
> fortran_compiler : None
> nb_core : 48
> OLP : MadLoop
> auto_update : 7
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> golem : None
> cpp_compiler : None
> notification_center : True
> run_mode : 2
> #******
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.7 2019-10-16 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https:/
> #* *
> #******
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #******
> set default_
> set group_subprocesses Auto
> set ignore_
> set loop_optimized_
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_
> 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/
> #######
> ## 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(
> 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/
> #
> #******
> # Decay widths *
> #******
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.491500e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.441404e+00
> #
> # PDG Width
> DECAY 24 2.047600e+00
> #
> # PDG Width
> DECAY 25 8.485861e-04
> # BR NDA ID1 ID2 ...
> 1.519292e-01 3 -2 1 24 # 0.000128925007304
> 1.519292e-01 3 -4 3 24 # 0.000128925007304
> 1.504738e-01 3 -24 -1 2 # 0.000127689975094
> 1.504738e-01 3 -24 -3 4 # 0.000127689975094
> 5.024122e-02 3 -12 11 24 # 4.2634000939e-05
> 5.024122e-02 3 -14 13 24 # 4.2634000939e-05
> 5.007860e-02 3 -24 -11 12 # 4.24960038675e-05
> 5.007860e-02 3 -24 -13 14 # 4.24960038675e-05
> 4.968264e-02 3 -16 15 24 # 4.21599977153e-05
> 4.961429e-02 3 -24 -15 16 # 4.21019968554e-05
> 1.488829e-02 3 -1 1 23 # 1.26339959468e-05
> 1.488829e-02 3 -3 3 23 # 1.26339959468e-05
> 1.194340e-02 3 -5 5 23 # 1.01350032267e-05
> 1.176015e-02 3 -2 2 23 # 9.97949982391e-06
> 1.176015e-02 3 -4 4 23 # 9.97949982391e-06
> 6.692701e-03 3 -12 12 23 # 5.67933304006e-06
> 6.692701e-03 3 -14 14 23 # 5.67933304006e-06
> 6.692701e-03 3 -16 16 23 # 5.67933304006e-06
> 3.333663e-03 3 -11 11 23 # 2.82890008388e-06
> 3.333663e-03 3 -13 13 23 # 2.82890008388e-06
> 3.271677e-03 3 -15 15 23 # 2.77629962589e-06
> #
> # PDG Width
> DECAY 82 0.000000e+00
> #******
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat aMC@NLO *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/
> # *
> # 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=
> #******
> # 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/
> # dynamical_
> #******
> 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_
> ! 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_
> ! 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://
> # 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*
> # 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://
> #******
> 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> #******
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #12 |
There are two log file in directory GF5. One is log.txt:
=======
INFO: MadFKS read these parameters from FKS_params.dat
======
> IRPoleCheckThre
> PrecisionVirtua
> NHelForMCoverHels = 4
> VirtualFraction = 1.0000000000000000
> MinVirtualFraction = 5.0000000000000
======
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_
*****
New value of alpha_s from PDF nn23nlo: 0.11899999999999999
******
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
******
External Params
-----
MU_R = 91.188000000000002
aEWM1 = 132.50700000000001
mdl_Gf = 1.1663900000000
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.4858609999999
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.262400000000
mdl_kHdwR = -0.262400000000
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.26240000000
mdl_nb_
mdl_Ncol__exp__2 = 9.0000000000000000
mdl_MB__exp__2 = 22.090000000000003
mdl_MT__exp__2 = 29929.000000000000
mdl_conjg__kHdw = (-0.26240000000
mdl_aEW = 7.5467711139788
mdl_MW = 80.419002445756163
mdl_sqrt__aEW = 8.6872153846781
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.6995554250465
mdl_yt = 0.99366614581500623
mdl_ytau = 1.0206617000654
mdl_muH = 88.388347648318430
mdl_AxialZUp = -0.185177018617
mdl_AxialZDown = 0.18517701861793787
mdl_VectorZUp = 7.5430507588273
mdl_VectorZDown = -0.130303763103
mdl_VectorAUp = 0.20530251149624587
mdl_VectorADown = -0.102651255748
mdl_VectorWmDxU = 0.23095271737156670
mdl_AxialWmDxU = -0.230952717371
mdl_VectorWpUxD = 0.23095271737156670
mdl_AxialWpUxD = -0.230952717371
mdl_I1x33 = ( 2.6995554250465
mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
mdl_I4x33 = ( 2.6995554250465
mdl_Vector_tbGp = (-0.96667059156
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.4835522759998
mdl_gAaa = 1.3008566310666
mdl_cw__exp__2 = 0.77775351421422245
mdl_gAza = 4.7794971072590
mdl_gHaa = 2.5475109025056
mdl_gHza = 3.9182129211851
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_R2MixedFac
mdl_GWcft_
mdl_GWcft_
mdl_bWcft_UV_1EPS_ = -1.878028328484
mdl_tWcft_UV_1EPS_ = -1.878028328484
mdl_G__exp__3 = 1.8056676068262196
mdl_MU_R__exp__2 = 8315.2513440000002
mdl_GWcft_
mdl_GWcft_
mdl_bWcft_UV_FIN_ = -0.136421009473
mdl_tWcft_UV_FIN_ = -9.877821144346
mdl_gAgg = 7.6274879753643
mdl_gHgg = -5.084991983576
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.9999999999999
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.100000000000
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.100000000000
Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 5
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
1 1 1
'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
Enter 'born0' or 'virt0' to perform
a pure n-body integration (no S functions)
doing the all of this channel
Normal integration (Sfunction != 1)
about to integrate 13 -1 12 5
imode is 0
channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
STOP 1
------- iteration 1
Update # PS points (even): 1040 --> 1040
Using random seed offsets: 5 , 10 , 0
with seed 35
Ranmar initialization seeds 14391 9419
Total number of FKS directories is 4
FKS process map (sum= 3 ) :
1 --> 2 : 1 3
2 --> 2 : 2 4
======
process combination map (specified per FKS dir):
1 map 1 2
1 inv. map 1 2
2 map 1 2
2 inv. map 1 2
3 map 1 2
3 inv. map 1 2
4 map 1 2
4 inv. map 1 2
======
tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
bpower is 0.0000000000000000
Scale values (may change event by event):
muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
QES, QES_reference: 0.198743D+03 0.198743D+03 1.00
muR_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
muF1_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
muF2_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
QES_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
alpha_s= 0.10641796709416652
ERROR: INTEGRAL APPEARS TO BE ZERO.
TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
Time in seconds: 4
Another is log_MINT0.txt:
======
INFO: MadFKS read these parameters from FKS_params.dat
======
> IRPoleCheckThre
> PrecisionVirtua
> NHelForMCoverHels = 4
> VirtualFraction = 1.0000000000000000
> MinVirtualFraction = 5.0000000000000
======
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_
*****
New value of alpha_s from PDF nn23nlo: 0.11899999999999999
******
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
******
External Params
-----
MU_R = 91.188000000000002
aEWM1 = 132.50700000000001
mdl_Gf = 1.1663900000000
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.4858609999999
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.262400000000
mdl_kHdwR = -0.262400000000
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.26240000000
mdl_nb_
mdl_Ncol__exp__2 = 9.0000000000000000
mdl_MB__exp__2 = 22.090000000000003
mdl_MT__exp__2 = 29929.000000000000
mdl_conjg__kHdw = (-0.26240000000
mdl_aEW = 7.5467711139788
mdl_MW = 80.419002445756163
mdl_sqrt__aEW = 8.6872153846781
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.6995554250465
mdl_yt = 0.99366614581500623
mdl_ytau = 1.0206617000654
mdl_muH = 88.388347648318430
mdl_AxialZUp = -0.185177018617
mdl_AxialZDown = 0.18517701861793787
mdl_VectorZUp = 7.5430507588273
mdl_VectorZDown = -0.130303763103
mdl_VectorAUp = 0.20530251149624587
mdl_VectorADown = -0.102651255748
mdl_VectorWmDxU = 0.23095271737156670
mdl_AxialWmDxU = -0.230952717371
mdl_VectorWpUxD = 0.23095271737156670
mdl_AxialWpUxD = -0.230952717371
mdl_I1x33 = ( 2.6995554250465
mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
mdl_I4x33 = ( 2.6995554250465
mdl_Vector_tbGp = (-0.96667059156
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.4835522759998
mdl_gAaa = 1.3008566310666
mdl_cw__exp__2 = 0.77775351421422245
mdl_gAza = 4.7794971072590
mdl_gHaa = 2.5475109025056
mdl_gHza = 3.9182129211851
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_R2MixedFac
mdl_GWcft_
mdl_GWcft_
mdl_bWcft_UV_1EPS_ = -1.878028328484
mdl_tWcft_UV_1EPS_ = -1.878028328484
mdl_G__exp__3 = 1.8056676068262196
mdl_MU_R__exp__2 = 8315.2513440000002
mdl_GWcft_
mdl_GWcft_
mdl_bWcft_UV_FIN_ = -0.136421009473
mdl_tWcft_UV_FIN_ = -9.877821144346
mdl_gAgg = 7.6274879753643
mdl_gHgg = -5.084991983576
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.9999999999999
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.100000000000
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.100000000000
Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 5
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
1 1 1
'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
Enter 'born0' or 'virt0' to perform
a pure n-body integration (no S functions)
doing the all of this channel
Normal integration (Sfunction != 1)
about to integrate 13 -1 12 5
imode is 0
channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
STOP 1
------- iteration 1
Update # PS points (even): 1040 --> 1040
Using random seed offsets: 5 , 10 , 0
with seed 35
Ranmar initialization seeds 14391 9419
Total number of FKS directories is 4
FKS process map (sum= 3 ) :
1 --> 2 : 1 3
2 --> 2 : 2 4
======
process combination map (specified per FKS dir):
1 map 1 2
1 inv. map 1 2
2 map 1 2
2 inv. map 1 2
3 map 1 2
3 inv. map 1 2
4 map 1 2
4 inv. map 1 2
======
tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
bpower is 0.0000000000000000
Scale values (may change event by event):
muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
QES, QES_reference: 0.198743D+03 0.198743D+03 1.00
muR_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
muF1_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
muF2_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
QES_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
alpha_s= 0.10641796709416652
ERROR: INTEGRAL APPEARS TO BE ZERO.
TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
Time in seconds: 4
Hi,
Looks like that you have zero coupling in your model:
> GC_3008a = 0.00000E+00 0.00000E+00
> GC_3008h = -0.00000E+00 -0.00000E+00
@NLO the code is crashing as soon as one integral provides 0 result.
In this case this is likely due to the face that such coupling are zero.
I would advise to use a restriction card to remove such interaction from your model
and then retry.
Cheers,
Olivier
PS: https:/
> On 28 Nov 2019, at 16:02, Yue Xu <email address hidden> wrote:
>
> Question #685744 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Yue Xu is still having a problem:
> There are two log file in directory GF5. One is log.txt:
>
> =======
> INFO: MadFKS read these parameters from FKS_params.dat
> =======
>> IRPoleCheckThre
>> PrecisionVirtua
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000
> =======
> 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_
> *******
> New value of alpha_s from PDF nn23nlo: 0.11899999999999999
> *******
> * MadGraph/MadEvent *
> * -------
> * http://
> * http://
> * http://
> * -------
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *******
>
> External Params
> -------
>
> MU_R = 91.188000000000002
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000
> 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.4858609999999
> 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.262400000000
> mdl_kHdwR = -0.262400000000
> 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.26240000000
> mdl_nb_
> mdl_Ncol__exp__2 = 9.0000000000000000
> mdl_MB__exp__2 = 22.090000000000003
> mdl_MT__exp__2 = 29929.000000000000
> mdl_conjg__kHdw = (-0.26240000000
> mdl_aEW = 7.5467711139788
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781
> 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.6995554250465
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654
> mdl_muH = 88.388347648318430
> mdl_AxialZUp = -0.185177018617
> mdl_AxialZDown = 0.18517701861793787
> mdl_VectorZUp = 7.5430507588273
> mdl_VectorZDown = -0.130303763103
> mdl_VectorAUp = 0.20530251149624587
> mdl_VectorADown = -0.102651255748
> mdl_VectorWmDxU = 0.23095271737156670
> mdl_AxialWmDxU = -0.230952717371
> mdl_VectorWpUxD = 0.23095271737156670
> mdl_AxialWpUxD = -0.230952717371
> mdl_I1x33 = ( 2.6995554250465
> mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995554250465
> mdl_Vector_tbGp = (-0.96667059156
> 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.4835522759998
> mdl_gAaa = 1.3008566310666
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_gAza = 4.7794971072590
> mdl_gHaa = 2.5475109025056
> mdl_gHza = 3.9182129211851
> 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_R2MixedFact
> mdl_GWcft_
> mdl_GWcft_
> mdl_bWcft_UV_1EPS_ = -1.878028328484
> mdl_tWcft_UV_1EPS_ = -1.878028328484
> mdl_G__exp__3 = 1.8056676068262196
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_GWcft_UV_b_FIN_ = -1.856343862667
> mdl_GWcft_UV_t_FIN_ = 4.0087659331150
> mdl_bWcft_UV_FIN_ = -0.136421009473
> mdl_tWcft_UV_FIN_ = -9.877821144346
> mdl_gAgg = 7.6274879753643
> mdl_gHgg = -5.084991983576
> 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.9999999999999
> 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.100000000000
> 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.100000000000
> Doing the S and H events together
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Do MC over helicities for the virtuals
> Enter Configuration Number:
> Running Configuration Number: 5
> Enter running mode for MINT:
> 0 to set-up grids, 1 to integrate, 2 to generate events
> MINT running mode: 0
> Set the three folding parameters for MINT
> xi_i, phi_i, y_ij
> 1 1 1
> 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
> Enter 'born0' or 'virt0' to perform
> a pure n-body integration (no S functions)
> doing the all of this channel
> Normal integration (Sfunction != 1)
> about to integrate 13 -1 12 5
> imode is 0
> channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
> STOP 1
> ------- iteration 1
> Update # PS points (even): 1040 --> 1040
> Using random seed offsets: 5 , 10 , 0
> with seed 35
> Ranmar initialization seeds 14391 9419
> Total number of FKS directories is 4
> FKS process map (sum= 3 ) :
> 1 --> 2 : 1 3
> 2 --> 2 : 2 4
> =======
> process combination map (specified per FKS dir):
> 1 map 1 2
> 1 inv. map 1 2
> 2 map 1 2
> 2 inv. map 1 2
> 3 map 1 2
> 3 inv. map 1 2
> 4 map 1 2
> 4 inv. map 1 2
> =======
> tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
> bpower is 0.0000000000000000
> Scale values (may change event by event):
> muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
> muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
> muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
> QES, QES_reference: 0.198743D+03 0.198743D+03 1.00
>
> muR_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF1_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF2_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> QES_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
>
> alpha_s= 0.10641796709416652
> ERROR: INTEGRAL APPEARS TO BE ZERO.
> TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
> Time in seconds: 4
>
>
>
>
>
>
>
>
> Another is log_MINT0.txt:
>
> =======
> INFO: MadFKS read these parameters from FKS_params.dat
> =======
>> IRPoleCheckThre
>> PrecisionVirtua
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000
> =======
> 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_
> *******
> New value of alpha_s from PDF nn23nlo: 0.11899999999999999
> *******
> * MadGraph/MadEvent *
> * -------
> * http://
> * http://
> * http://
> * -------
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *******
>
> External Params
> -------
>
> MU_R = 91.188000000000002
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000
> 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.4858609999999
> 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.262400000000
> mdl_kHdwR = -0.262400000000
> 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.26240000000
> mdl_nb_
> mdl_Ncol__exp__2 = 9.0000000000000000
> mdl_MB__exp__2 = 22.090000000000003
> mdl_MT__exp__2 = 29929.000000000000
> mdl_conjg__kHdw = (-0.26240000000
> mdl_aEW = 7.5467711139788
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781
> 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.6995554250465
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654
> mdl_muH = 88.388347648318430
> mdl_AxialZUp = -0.185177018617
> mdl_AxialZDown = 0.18517701861793787
> mdl_VectorZUp = 7.5430507588273
> mdl_VectorZDown = -0.130303763103
> mdl_VectorAUp = 0.20530251149624587
> mdl_VectorADown = -0.102651255748
> mdl_VectorWmDxU = 0.23095271737156670
> mdl_AxialWmDxU = -0.230952717371
> mdl_VectorWpUxD = 0.23095271737156670
> mdl_AxialWpUxD = -0.230952717371
> mdl_I1x33 = ( 2.6995554250465
> mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995554250465
> mdl_Vector_tbGp = (-0.96667059156
> 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.4835522759998
> mdl_gAaa = 1.3008566310666
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_gAza = 4.7794971072590
> mdl_gHaa = 2.5475109025056
> mdl_gHza = 3.9182129211851
> 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_R2MixedFact
> mdl_GWcft_
> mdl_GWcft_
> mdl_bWcft_UV_1EPS_ = -1.878028328484
> mdl_tWcft_UV_1EPS_ = -1.878028328484
> mdl_G__exp__3 = 1.8056676068262196
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_GWcft_UV_b_FIN_ = -1.856343862667
> mdl_GWcft_UV_t_FIN_ = 4.0087659331150
> mdl_bWcft_UV_FIN_ = -0.136421009473
> mdl_tWcft_UV_FIN_ = -9.877821144346
> mdl_gAgg = 7.6274879753643
> mdl_gHgg = -5.084991983576
> 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.9999999999999
> 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.100000000000
> 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.100000000000
> Doing the S and H events together
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Do MC over helicities for the virtuals
> Enter Configuration Number:
> Running Configuration Number: 5
> Enter running mode for MINT:
> 0 to set-up grids, 1 to integrate, 2 to generate events
> MINT running mode: 0
> Set the three folding parameters for MINT
> xi_i, phi_i, y_ij
> 1 1 1
> 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
> Enter 'born0' or 'virt0' to perform
> a pure n-body integration (no S functions)
> doing the all of this channel
> Normal integration (Sfunction != 1)
> about to integrate 13 -1 12 5
> imode is 0
> channel 1 : 5 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
> STOP 1
> ------- iteration 1
> Update # PS points (even): 1040 --> 1040
> Using random seed offsets: 5 , 10 , 0
> with seed 35
> Ranmar initialization seeds 14391 9419
> Total number of FKS directories is 4
> FKS process map (sum= 3 ) :
> 1 --> 2 : 1 3
> 2 --> 2 : 2 4
> =======
> process combination map (specified per FKS dir):
> 1 map 1 2
> 1 inv. map 1 2
> 2 map 1 2
> 2 inv. map 1 2
> 3 map 1 2
> 3 inv. map 1 2
> 4 map 1 2
> 4 inv. map 1 2
> =======
> tau_min 1 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 2 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 3 1 : 0.17516E+03 -- 0.20542E+03
> tau_min 4 1 : 0.17516E+03 -- 0.20542E+03
> bpower is 0.0000000000000000
> Scale values (may change event by event):
> muR, muR_reference: 0.198743D+03 0.198743D+03 1.00
> muF1, muF1_reference: 0.198743D+03 0.198743D+03 1.00
> muF2, muF2_reference: 0.198743D+03 0.198743D+03 1.00
> QES, QES_reference: 0.198743D+03 0.198743D+03 1.00
>
> muR_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF1_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF2_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> QES_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
>
> alpha_s= 0.10641796709416652
> ERROR: INTEGRAL APPEARS TO BE ZERO.
> TRIED 100880 PS POINTS AND ONLY 0 GAVE A NON-ZERO INTEGRAND.
> Time in seconds: 4
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Yue Xu (yuexu) said : | #14 |
Hi Olivier,
I added the restriction card and the couplings not interested are removed, then it works well.
Thanks a lot !
Cheers,
Yue
