Two dark matter components

Asked by Y Zhong

Hi,

I am working on a model of excited DM. In the model, I have two dark matter components: chi and psi, with a scalar mediator s0 connecting to the SM. I would like to compute the relic density of chi and psi.

Under the maddm 2.1 MG5 interface, I have tried

import model my_model
define darkmatter chi psi
generate relic_density
output my_output

It seems only chi is identified as DM. I was wondering if there is something I have missed? The output is attached in the end of the email.

Thank you very much!

INFO: Update the dependent parameter of the param_card.dat
WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...)
INFO: Found Dark Matter candidate: chi
add model /project/splash/mg5_simulation/mg25/MG5_aMC_v2_5_5/PLUGIN/maddm/EffOperators/REAL ~fdm=chi --recreate --keep_decay
INFO: The two model defines the block 'MASS' with id '[3001]' with different parameter name 'fdmm', 'MChi'
     We will merge those two parameters in a single one
INFO: Adding interaction for the following particles: [-1, 1, 999000007]
INFO: Adding interaction for the following particles: [-3, 3, 999000007]
INFO: Adding interaction for the following particles: [-5, 5, 999000007]
INFO: Adding interaction for the following particles: [-1, 1, 999000006, 999000006]
INFO: Adding interaction for the following particles: [-3, 3, 999000006, 999000006]
INFO: Adding interaction for the following particles: [-5, 5, 999000006, 999000006]
INFO: Adding interaction for the following particles: [3001, 3001, 999000007]
INFO: Adding interaction for the following particles: [-2, 2, 999000007]
INFO: Adding interaction for the following particles: [-4, 4, 999000007]
INFO: Adding interaction for the following particles: [-6, 6, 999000007]
INFO: Adding interaction for the following particles: [-2, 2, 999000006, 999000006]
INFO: Adding interaction for the following particles: [-4, 4, 999000006, 999000006]
INFO: Adding interaction for the following particles: [-6, 6, 999000006, 999000006]
INFO: Adding interaction for the following particles: [-1, 1, 999000009, 999000009]
INFO: Adding interaction for the following particles: [-3, 3, 999000009, 999000009]
INFO: Adding interaction for the following particles: [-5, 5, 999000009, 999000009]
INFO: Adding interaction for the following particles: [-2, 2, 999000009, 999000009]
INFO: Adding interaction for the following particles: [-4, 4, 999000009, 999000009]
INFO: Adding interaction for the following particles: [-6, 6, 999000009, 999000009]
INFO: Adding interaction for the following particles: [-1, 1, 999000010]
INFO: Adding interaction for the following particles: [-3, 3, 999000010]
INFO: Adding interaction for the following particles: [-5, 5, 999000010]
INFO: Adding interaction for the following particles: [3001, 3001, 999000010]
INFO: Adding interaction for the following particles: [-2, 2, 999000010]
INFO: Adding interaction for the following particles: [-4, 4, 999000010]
INFO: Adding interaction for the following particles: [-6, 6, 999000010]
INFO: load particles
INFO: load vertices
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles l- / bsm / j / vl / l+ / p / vl~ unchanged
Defined multiparticle all = g a ve vm vt ve~ vm~ vt~ u c t d s b u~ c~ t~ d~ s~ b~ s0 z w+ h ~sdm fsdm ~vdm fvdm w- chi psi e- mu- ta- e+ mu+ ta+
INFO: No coannihilation partners found.
INFO: DM is allowed to annihilate into the following BSM particles: s0 ~sdm fsdm ~vdm fvdm
INFO: use generate relic_density / X to forbid the decay the annihilation to X
INFO: if you want to forbid DM to annihilate into BSM particles do relic_density / bsm
Multiparticle labels:
all = g a ve vm vt ve~ vm~ vt~ u c t d s b u~ c~ t~ d~ s~ b~ s0 z w+ h ~sdm fsdm ~vdm fvdm w- chi psi e- mu- ta- e+ mu+ ta+
l- = e- mu-
dm_particle_0 = chi
bsm = s0 chi psi
dm_particles = chi
j = g u c d s u~ c~ d~ s~
vl = ve vm vt
l+ = e+ mu+
p = g u c d s u~ c~ d~ s~
vl~ = ve~ vm~ vt~
INFO: Trying process: chi chi > g g QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g a QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g z QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g h QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g s0 QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g ~sdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > g fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a a QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a z QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a h QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a s0 QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a ~sdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > a fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ve ve~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ve vm~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ve vt~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > vm ve~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > vm vm~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > vm vt~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > vt ve~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > vt vm~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > vt vt~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > d d~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > d s~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > d b~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > u u~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > u c~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > u t~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s d~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s s~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s b~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > c u~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > c c~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > c t~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > b d~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > b s~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > b b~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > t u~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > t c~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > t t~ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z z QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z h QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z s0 QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z ~sdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > z fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > w+ w- QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > h h QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > h s0 QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > h ~sdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > h fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > h ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > h fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s0 s0 QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Process has 2 diagrams
INFO: Trying process: chi chi > s0 ~sdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s0 fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s0 ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > s0 fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ~sdm ~sdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ~sdm fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ~sdm ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ~sdm fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > fsdm fsdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > fsdm ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > fsdm fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ~vdm ~vdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ~vdm fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > fvdm fvdm QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > e- e+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > e- mu+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > e- ta+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > mu- e+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > mu- mu+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > mu- ta+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ta- e+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ta- mu+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
INFO: Trying process: chi chi > ta- ta+ QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1999
1 processes with 2 diagrams generated in 0.201 s
Total: 1 processes with 2 diagrams
INFO: Trying process: chi chi > chi chi QED<=4 SDEFFF=0 SDEFFV=0 SIEFFF=0 SIEFFS=0 SIEFFV=0 @1998

Question information

Language:
English Edit question
Status:
Answered
For:
MadDM Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

looks like the command define darkmatter chi psi is actually ignored.

Here is a patch which should solve the problem (hoping formatting --number of space-- will not be modifed):
=== modified file 'maddm_interface.py'
--- maddm_interface.py 2017-08-09 13:25:48 +0000
+++ maddm_interface.py 2017-08-24 12:34:37 +0000
@@ -117,11 +117,15 @@
                     self.search_dm_candidate([])
                 elif args[1].startswith('/'):
                     self.search_dm_candidate([a.replace('/', '') for a in args[1:]])
- elif len(args)==2:
- self._dm_candidate = [self._curr_model.get_particle(args[1])]
- if not self._dm_candidate[0]:
+ elif len(args)>=2:
+ self._dm_candidate = [self._curr_model.get_particle(dm) for dm in args[1:]]
+ #remove None
+ self._dm_candidate = [dm for dm in self._dm_candidate if dm]
+ if not self._dm_candidate:
                         raise DMError, '%s is not a valid particle for the model.' % args[1]
- self.update_model_with_EFT()
+ if len(self._dm_candidate) == 1:
+ # No update of the model if 2(or more) DM since DD is not possible
+ self.update_model_with_EFT()

             elif args[0] == 'coannihilator':
                 if not self._dm_candidate:

If you can send me your model on my email (<email address hidden>), I would be able to test this.
Otherwise, you can try the patch yourself.

Cheers,

Olivier

Revision history for this message
Y Zhong (ymzhong) said :
#2

Hi Olivier,

Thanks a lot for the patch! Now I found both psi and chi are identified as DM particle. However, Maddm fails when I try "generate relic_density".

Command "generate relic_density" interrupted with error:
KeyError : '1'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

 In below please find information from "MG5_debug"

*************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.1 2017-12-12 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set low_mem_multicore_nlo_generation False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model /projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/models/ExDarkLepton_UFO
define darkmatter chi psi
generate relic_density
Traceback (most recent call last):
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/extended_cmd.py", line 1467, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/master_interface.py", line 269, in do_generate
    return self.cmd.do_generate(self, line, *args, **opts)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/madgraph_interface.py", line 4446, in do_generate
    self.do_add(" ".join(args))
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/PLUGIN/maddm/maddm_interface.py", line 352, in do_add
    return self.generate_relic([])
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/PLUGIN/maddm/maddm_interface.py", line 546, in generate_relic
    self.do_generate(proc)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/master_interface.py", line 269, in do_generate
    return self.cmd.do_generate(self, line, *args, **opts)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/madgraph_interface.py", line 4446, in do_generate
    self.do_add(" ".join(args))
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/PLUGIN/maddm/maddm_interface.py", line 370, in do_add
    return super(MadDM_interface, self).do_add(line)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/master_interface.py", line 220, in do_add
    return self.cmd.do_add(self, line, *args, **opts)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/interface/madgraph_interface.py", line 3079, in do_add
    optimize=optimize, diagram_filter=diagram_filter)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/core/diagram_generation.py", line 1595, in __init__
    self.get('amplitudes')
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/core/diagram_generation.py", line 1643, in get
    diagram_filter=self['diagram_filter']))
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/core/diagram_generation.py", line 1817, in generate_multi_amplitudes
    result = amplitude.generate_diagrams(diagram_filter=diagram_filter)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/core/diagram_generation.py", line 783, in generate_diagrams
    diagram.calculate_orders(model)
  File "/projectnb/leptoq/mg5_simulation/mg90/MG5_aMC_v2_6_1/madgraph/core/base_objects.py", line 2461, in calculate_orders
    (c,n) in couplings.items()])
KeyError: '1'
                          MadGraph5_aMC@NLO Options
                          ----------------
        complex_mass_scheme : False
                      gauge : unitary
         group_subprocesses : Auto
  ignore_six_quark_processes : False
           loop_color_flows : False
      loop_optimized_output : True
  low_mem_multicore_nlo_generation : False
     max_npoint_for_channel : 0
               stdout_level : 20 (user set)

                         MadEvent Options
                          ----------------
     automatic_html_opening : True
                    nb_core : 16 (user set)
        notification_center : True
                   run_mode : 2

                      Configuration Options
                      ---------------------
                        OLP : MadLoop
                    amcfast : amcfast-config
                   applgrid : applgrid-config
                auto_update : 7
         cluster_local_path : None
           cluster_nb_retry : 1
              cluster_queue : None (user set)
         cluster_retry_wait : 300
               cluster_size : 100
      cluster_status_update : (600, 30)
          cluster_temp_path : None
               cluster_type : condor
                    collier : ./HEPTools/lib
               cpp_compiler : None
             crash_on_error : False
               delphes_path : ./Delphes
                 eps_viewer : None
        exrootanalysis_path : ./ExRootAnalysis
              f2py_compiler : None
                    fastjet : fastjet-config
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : /projectnb/leptoq/mg5_simulation/mg1/MG5_aMC_v2_6_1/HEPTools/lhapdf6/bin/lhapdf-config (user set)
          madanalysis5_path : None (user set)
           madanalysis_path : ./MadAnalysis
  mg5amc_py8_interface_path : None (user set)
                      ninja : ./HEPTools/lib
        output_dependencies : external
                      pjfry : None (user set)
            pythia-pgs_path : ./pythia-pgs
               pythia8_path : None (user set)
                    samurai : None
               syscalc_path : ./SysCalc
                    td_path : ./td
                text_editor : None
                thepeg_path : None (user set)
                    timeout : 60
                web_browser : None

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi,

That means that your UFO model does not follow the normal convention.
So you should contact the author of your model and ask them to fix it.
The problem is that some coupling are defined to a coupling order named "1"
but that name for coupling order is not defined in couplings_order.py
(actually that name itself is not valid since any name should start by a letter since those are used as python variable)

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Y Zhong for more information if necessary.

To post a message you must log in.