Error running madspin - KeyError : ((-5, 21), (-6, -5, 6))
Hello,
I have encountered a problem and I am not sure how to resolve this.
I am trying to do the following process in MadSpin:
p p > t t~, (t~ > b~ w- , w- > l- vl~ ), t > b h+
I am using the 2HDMtypeII model for this purpose.
I have an error as follows
Command "launch" interrupted with error:
KeyError : ((-5, 21), (-6, -5, 6))
Please report this bug to developers
More information is found in 'MS_debug'.
Please attach this file to your report.
It would be very greate if you would help me to resolve this problem.
Thank you!
Anna
My Madspin card:
set Nevents_
set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
set seed 765432
import run.lhe
set ms_dir MadSpin
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
decay t > b h+
decay t~ > b~ w-, w- > l- vl~
launch
My param_card:
#######
## INFORMATION FOR MASS
#######
Block mass
6 1.725000e+02 # MT
23 9.118760e+01 # MZ
25 125.0 # mh1 #
5 4.700000e+00 # MB
4 1.550000e+00 # MC
35 179.064232054 # mh2 #
36 179.064232054 # mh2 #
37 160 # mhc #
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000 # d : 0.0
2 0.000000 # u : 0.0
3 0.000000 # s : 0.0
# 4 0.000000 # c : 0.0
# 5 0.000000 # b : 0.0
11 0.000000 # e- : 0.0
12 0.000000 # ve : 0.0
13 0.000000 # mu- : 0.0
14 0.000000 # vm : 0.0
15 0.000000 # ta- : 0.0
16 0.000000 # vt : 0.0
21 0.000000 # g : 0.0
22 0.000000 # a : 0.0
24 80.39900 # w+ : cmath.sqrt(
82 0.000000 # ghg : 0.0
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.279000e+02 # aEWM1
2 1.166390e-05 # Gf
3 1.180000e-01 # aS
#######
## INFORMATION FOR YUKAWA
#######
Block yukawa
1 5.040000e-03 # ymdo
2 2.550000e-03 # ymup
3 1.010000e-01 # yms
4 1.270000e+00 # ymc
5 4.700000e+00 # ymb
6 1.725000e+02 # ymt
11 5.110000e-04 # yme
13 1.056600e-01 # ymm
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 1.000000e+00
DECAY 35 1.000000e+00
DECAY 36 1.000000e+00
DECAY 37 1.300000e+02 # whc
Block FRBlock
1 2.000000e+01 # beta
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
DECAY 82 0.000000 # ghg : 0.0
DECAY 9000001 0.000000 # gha : 0.0
DECAY 9000002 2.495200 # ghz : WZ
DECAY 9000003 2.085000 # ghwp : WW
DECAY 9000004 2.085000 # ghwm : WW
## DECAY ID WIDTH
## NDA ID1 ID2 ...
DECAY -24 2.085000e+00
3.377e-01 2 1 -2
3.377e-01 2 3 -4
1.082e-01 2 11 -12
1.082e-01 2 13 -14
1.082e-01 2 15 -16
DECAY 24 2.085000e+00
3.377e-01 2 -1 2
3.377e-01 2 -3 4
1.082e-01 2 -11 12
1.082e-01 2 -13 14
1.082e-01 2 -15 16
#With this, I can have small couplings to be consistent
# with the SM width and limit on BR (~ 0.1%), and yet generate the exotic width
#
DECAY -6 1.33
9.96e-01 2 -5 -24
0.10e-02 2 -4 25
0.10e-02 2 -2 25
0.05e-02 2 -4 23
0.05e-02 2 -2 23
0.05e-02 2 -4 22
0.05e-02 2 -2 22
0.05e-02 2 -5 -37
DECAY 6 1.33
9.96e-01 2 5 24
0.10e-02 2 4 25
0.10e-02 2 2 25
0.05e-02 2 4 23
0.05e-02 2 2 23
0.05e-02 2 4 22
0.05e-02 2 2 22
0.05e-02 2 5 37
#######
## INFORMATION FOR HIGGS
#######
Block Higgs
1 5.000000e-01 # l2
2 1.000000e+00 # l3
3 0.000000e+00 # lR7
5 3.000000e-01 # mixh
#######
## INFORMATION FOR LOOP
#######
Block LOOP
1 9.118800e+01 # MU_R
#======
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#======
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/
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/
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/
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/
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/
Block QNUMBERS 37 # h+
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/
Block QNUMBERS 35 # h2
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/
Block QNUMBERS 36 # h3
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/
My MS_debug report:
set Nevents_
set max_weight_ps_point 1000
set seed 765432
import run.lhe
set ms_dir MadSpin
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
decay t > b h+
decay t~ > b~ w-, w- > l- vl~
launch
Traceback (most recent call last):
File "/cvmfs/
return self.onecmd_
File "/cvmfs/
return func(arg, **opt)
File "/cvmfs/
out = f(self, *args, **opt)
File "/cvmfs/
generate_
File "/cvmfs/
self.
File "/cvmfs/
production_tag, event_map = self.load_event()
File "/cvmfs/
P_order = self.all_
KeyError: ((-5, 21), (-6, -5, 6))
Value of current Options:
max_
Nevents_
max_
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Anna Ivina for more information if necessary.