Another bug?
Hi! I am continuing to experience issues with maddump (see my previous question https:/
After resolving this bug (please see https:/
import model DM_mesons_2
import_events decay example_
decay pi0 > y1 a, y1 > xd xd~
define darkmatter xd
add process interaction @DIS
add process interaction @electron
output DP_electron
launch DP_electron
set flux_norm 2.0e20
set prod_xsec_in_norm false
set d_target_detector 5650.0
set detector_density 3.72
set off-axis false
set parallelepiped True
set x_side 187.0
set y_side 69.0
set depth 200.0
set ncores 16
set testplot True
set gvd11 -3.333333e-4
set gvu11 6.666666e-4
set gvd22 -3.333333e-4
set gvu22 6.666666e-4
set gvd33 -3.333333e-4
set gvu33 6.666666e-4
set gvl11 -1.000000e-3
set gvl22 -1.000000e-3
set gvl33 -1.000000e-3
set my1 scan1:[0.01*i for i in range(1,14)]
set mxd scan1:[0.01/3.*i for i in range(1,14)]
set wy1 auto
Exception :
Please report this bug to developers
I attach the log and debug files below:
_______
The log is
_______
INFO: Create the next param_card in the scan definition
INFO: change parameter mass with code (5000001,) to 0.01
INFO: change parameter mass with code (5000521,) to 0.00333333333333
INFO: Computing the width set on auto in the param_card.dat
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.3.0 2021-11-12 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https:/
* and *
* http://
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
*******
load MG5 configuration from ../input/
set fastjet to fastjet-config
lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/
Note that you can still compile and run aMC@NLO with the built-in PDFs
MG5_aMC> set lhapdf /PATH/TO/
None does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/
Note that you can still compile and run aMC@NLO with the built-in PDFs
MG5_aMC> set lhapdf /PATH/TO/
Using default text editor "vi". Set another one in ./input/
set automatic_
/home/name/
import model /home/name/
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt mu- ve~ vm~ vt~ mu+ t t~ pi0 pi+ eta omega rho etaprime z w+ ghz ghwp ghwm h xr xc y1 pi- w- ghz~ ghwp~ ghwm~ xc~ e- ta- xd e+ ta+ xd~
compute_widths 5000001 --path=
Please note that the automatic computation of the width is
only valid in narrow-width approximation and at tree-level.
INFO: get decay diagram for y1
Vertexlist of this model has not been searched.
Found 8 stable particles
Pass to numerical integration for computing the widths:
INFO: More info in temporary files:
/home/
INFO: set output information to level: 20
INFO: compile directory
INFO: Running Survey
INFO: P0_y1_xdxdx
INFO: P0_y1_ccx
INFO: P0_y1_uux
INFO: P0_y1_mummup
INFO: P0_y1_emep
INFO: P0_y1_bbx
INFO: P0_y1_ddx
INFO: P0_y1_ssx
INFO: P0_y1_bdx
INFO: P0_y1_vevex
INFO: P0_y1_vmvmx
INFO: P0_y1_vtvtx
INFO: Idle: 0, Running: 4, Completed: 8 [ current time: 17h03 ]
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
INFO: Idle: 0, Running: 3, Completed: 9 [ 4s ]
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
INFO: Idle: 0, Running: 0, Completed: 12 [ 5.5s ]
INFO: End survey
INFO: Combining Events
INFO: fail to reach target 10000
INFO: storing files of previous run
INFO: Done
INFO:
Results written to DP_electron/
INFO: Update the dependent parameter of the param_card.dat
/home/name/
import model /home/name/
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt mu- ve~ vm~ vt~ mu+ t t~ pi0 pi+ eta omega rho etaprime z w+ ghz ghwp ghwm h xr xc y1 pi- w- ghz~ ghwp~ ghwm~ xc~ e- ta- xd e+ ta+ xd~
*******
* *
* W E L C O M E to M A D S P I N *
* *
*******
set spinmode none
set cross_section {0:1.0}
set new_wgt BR
set input_format hepmc
import DP_electron/
INFO: Setup the code for pure decay mode
import model /home/name/
INFO: detected model: /home/name/
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = mu+ e+
Defined multiparticle l- = mu- e-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
Pass the definition of 'j' and 'p' to 5 flavour scheme.
Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt mu- ve~ vm~ vt~ mu+ t t~ pi0 pi+ eta omega rho etaprime z w+ ghz ghwp ghwm h xr xc y1 pi- w- ghz~ ghwp~ ghwm~ xc~ e- ta- xd e+ ta+ xd~
define all = 21 82 -82 2 4 1 3 5 -2 -4 -1 -3 -5 22 9000001 -9000001 12 14 16 13 -12 -14 -16 -13 6 -6 111 211 221 223 113 331 23 24 9000002 9000003 9000004 25 5000511 5000512 5000001 -211 -24 -9000002 -9000003 -9000004 -5000512 11 15 5000521 -11 -15 -5000521
define all = 21 82 -82 2 4 1 3 5 -2 -4 -1 -3 -5 22 9000001 -9000001 12 14 16 13 -12 -14 -16 -13 6 -6 111 211 221 223 113 331 23 24 9000002 9000003 9000004 25 5000511 5000512 5000001 -211 -24 -9000002 -9000003 -9000004 -5000512 11 15 5000521 -11 -15 -5000521
Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt mu- ve~ vm~ vt~ mu+ t t~ pi0 pi+ eta omega rho etaprime z w+ ghz ghwp ghwm h xr xc y1 pi- w- ghz~ ghwp~ ghwm~ xc~ e- ta- xd e+ ta+ xd~
define l- = 11 13
define l- = 11 13
Defined multiparticle l- = mu- e-
define j = 21 2 4 1 3 -2 -4 -1 -3 5 -5
define j = 21 2 4 1 3 -2 -4 -1 -3 5 -5
Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
define vl = 12 14 16
define vl = 12 14 16
Defined multiparticle vl = ve vm vt
define l+ = -11 -13
define l+ = -11 -13
Defined multiparticle l+ = mu+ e+
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5
Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
define vl~ = -12 -14 -16
define vl~ = -12 -14 -16
Defined multiparticle vl~ = ve~ vm~ vt~
decay pi0 > y1 a, y1 > xd xd~
launch
INFO: Parsing input event to know how many decay to generate. This can takes few minuts.
INFO: Will use seed 233264884
INFO: generate 6425 decay event for particle pi0
INFO: Will use seed 159831646
decaying event number 100 [0.149754047394 s]
decaying event number 200 [0.333559036255 s]
decaying event number 300 [0.505009174347 s]
decaying event number 400 [0.67059803009 s]
decaying event number 500 [0.836493015289 s]
decaying event number 600 [0.995557069778 s]
decaying event number 700 [1.17664813995 s]
decaying event number 800 [1.34749507904 s]
decaying event number 900 [1.52117109299 s]
Exception :
Please report this bug to developers
More information is found in 'debug'.
Please attach this file to your report.
_______
The debug is
_______
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
original_
File "/home/
'unweighted
File "/home/
self.
File "/home/
for event in lhe_evts:
File "/home/
return self.next_event()
File "/home/
raise Exception
Exception
Value of current Options:
notificat
ignore_
loop_
cluster_
low_mem_
default_
acknowledged_
include_
automatic_
output_
nlo_
mg5amc_
exrootana
complex_
max_
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## INFORMATION FOR CKMBLOCK
#######
BLOCK CKMBLOCK #
1 2.277360e-01 # cabi
#######
## INFORMATION FOR DMINPUTS
#######
BLOCK DMINPUTS #
1 0.000000e+00 # gvxc
2 1.000000e+00 # gvxd
3 0.000000e+00 # gaxd
4 -3.333333e-04 # gvd11
5 6.666666e-04 # gvu11
6 -3.333333e-04 # gvd22
7 6.666666e-04 # gvu22
8 -3.333333e-04 # gvd33
9 6.666666e-04 # gvu33
10 -1.000000e-03 # gvl11
11 -1.000000e-03 # gvl22
12 -1.000000e-03 # gvl33
13 0.000000e+00 # gad11
14 0.000000e+00 # gau11
15 0.000000e+00 # gad22
16 0.000000e+00 # gau22
17 0.000000e+00 # gad33
18 0.000000e+00 # gau33
19 0.000000e+00 # gal11
20 0.000000e+00 # gal22
21 0.000000e+00 # gal33
22 0.000000e+00 # gnu11
23 0.000000e+00 # gnu22
24 0.000000e+00 # gnu33
25 0.000000e+00 # gvu31
26 0.000000e+00 # gau31
27 0.000000e+00 # gvd31
28 0.000000e+00 # gad31
29 0.000000e+00 # gvh
#######
## INFORMATION FOR MASS
#######
BLOCK MASS #
6 1.720000e+02 # mt
11 5.110000e-04 # me
15 1.777000e+00 # mta
23 9.118760e+01 # mz
25 1.250000e+02 # mh
111 1.349800e-01 # mpi0
113 7.752600e-01 # mrho
211 1.395706e-01 # mpi
221 5.478620e-01 # meta
223 7.826500e-01 # momega
331 9.577800e-01 # metaprime
5000001 scan1:[0.01*i for i in range(1,14)] # my1
5000511 1.000000e+01 # mxr
5000512 1.000000e+01 # mxc
5000521 scan1:[0.01/3.*i for i in range(1,14)] # mxd
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
5 0.000000e+00 # b : 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 7.982436e+01 # w+ : cmath.sqrt(
9000002 9.118760e+01 # ghz : mz
9000003 7.982436e+01 # ghwp : mw
9000004 7.982436e+01 # ghwm : mw
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
1 1.279000e+02 # aewm1
2 1.166370e-05 # gf
3 1.184000e-01 # as (note that parameter not used if you use a pdf set)
#######
## INFORMATION FOR YUKAWA
#######
BLOCK YUKAWA #
6 1.720000e+02 # ymt
15 1.777000e+00 # ymtau
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.508336e+00 # wt
DECAY 23 2.495200e+00 # wz
DECAY 24 2.085000e+00 # ww
DECAY 25 4.070000e-03 # wh
DECAY 113 1.491000e-01 # wrho
DECAY 5000001 auto # wy1
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 111 0.000000e+00 # pi0 : 0.0
DECAY 211 0.000000e+00 # pi+ : 0.0
DECAY 221 0.000000e+00 # eta : 0.0
DECAY 223 0.000000e+00 # omega : 0.0
DECAY 331 0.000000e+00 # etaprime : 0.0
DECAY 5000511 0.000000e+00 # xr : 0.0
DECAY 5000512 0.000000e+00 # xc : 0.0
DECAY 5000521 0.000000e+00 # xd : 0.0
DECAY 9000002 2.495200e+00 # ghz : wz
DECAY 9000003 2.085000e+00 # ghwp : ww
DECAY 9000004 2.085000e+00 # ghwm : ww
#######
## INFORMATION FOR QNUMBERS 111
#######
BLOCK QNUMBERS 111 # pi0
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 211
#######
BLOCK QNUMBERS 211 # pi+
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 221
#######
BLOCK QNUMBERS 221 # eta
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 223
#######
BLOCK QNUMBERS 223 # omega
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 113
#######
BLOCK QNUMBERS 113 # rho
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 331
#######
BLOCK QNUMBERS 331 # etaprime
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 9000001
#######
BLOCK QNUMBERS 9000001 # gha
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 9000002
#######
BLOCK QNUMBERS 9000002 # ghz
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 9000003
#######
BLOCK QNUMBERS 9000003 # ghwp
1 3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 9000004
#######
BLOCK QNUMBERS 9000004 # ghwm
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 82
#######
BLOCK QNUMBERS 82 # ghg
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/
#######
## INFORMATION FOR QNUMBERS 5000511
#######
BLOCK QNUMBERS 5000511 # xr
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 5000512
#######
BLOCK QNUMBERS 5000512 # xc
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 5000521
#######
BLOCK QNUMBERS 5000521 # xd
1 0 # 3 times electric charge
2 2 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#######
## INFORMATION FOR QNUMBERS 5000001
#######
BLOCK QNUMBERS 5000001 # y1
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
What can be wrong here?
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- maddump Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Andrew McAdams for more information if necessary.