USRMODERROR: Some inconsistency in the mass assignment in the model

[terryp]

Hi,
 I am trying to calculate relic abundance for a model using MadDM in MG5_aMC_v2.6.0. I receive the following error log when running ./maddm.py with my model. I have tried modifying the PDG code, but get the same error. I have included both the error log and the param card. Any help would be greatly appreciated.

Thanks,
Terry

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
MG5_debug start
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.0 2017-08-16 *
#* *
#* 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 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 GMIDM_00q
generate relic_density
Traceback (most recent call last):
  File "/Users/terryp/work/maddm/mg/madgraph/interface/extended_cmd.py", line 1438, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/master_interface.py", line 266, in do_generate
    return self.cmd.do_generate(self, line, *args, **opts)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/madgraph_interface.py", line 4342, in do_generate
    self.do_add(" ".join(args))
  File "/Users/terryp/work/maddm/mg/PLUGIN/maddm/maddm_interface.py", line 348, in do_add
    return self.generate_relic([])
  File "/Users/terryp/work/maddm/mg/PLUGIN/maddm/maddm_interface.py", line 491, in generate_relic
    self.search_dm_candidate(excluded_particles)
  File "/Users/terryp/work/maddm/mg/PLUGIN/maddm/maddm_interface.py", line 241, in search_dm_candidate
    self.update_model_with_EFT()
  File "/Users/terryp/work/maddm/mg/PLUGIN/maddm/maddm_interface.py", line 772, in update_model_with_EFT
    self.exec_cmd(mg5_command)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/Users/terryp/work/maddm/mg/PLUGIN/maddm/maddm_interface.py", line 366, in do_add
    return super(MadDM_interface, self).do_add(line)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/master_interface.py", line 217, in do_add
    return self.cmd.do_add(self, line, *args, **opts)
  File "/Users/terryp/work/maddm/mg/madgraph/interface/madgraph_interface.py", line 2968, in do_add
    return self.add_model(args[1:])
  File "/Users/terryp/work/maddm/mg/madgraph/interface/madgraph_interface.py", line 3125, in add_model
    base_model.add_model(path=model_path, identify_particles=identify)
  File "/Users/terryp/work/maddm/mg/models/usermod.py", line 969, in add_model
    self.add_particle(particle, identify=identify_particles[particle.name])
  File "/Users/terryp/work/maddm/mg/models/usermod.py", line 559, in add_particle
    return self.check_mass_width_of_particle(old_part, particle)
  File "/Users/terryp/work/maddm/mg/models/usermod.py", line 594, in check_mass_width_of_particle
    ' pdg_code: %s %s' % (p_base.pdg_code, p_plugin.pdg_code)
USRMODERROR: Some inconsistency in the mass assignment in the model
     Mass: MH0 and sdmm
     conflict name {'Gf': 'Gf__1', 'Wsdm': 'WH0'}
     pdg_code: 9900035 999000006
                          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 : False (user set)
                    nb_core : None
        notification_center : False (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                        OLP : MadLoop
                    amcfast : amcfast-config
                   applgrid : applgrid-config
                auto_update : 0 (user set)
         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 : None (user set)
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : lhapdf-config
          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 : bbedit (user set)
                thepeg_path : None (user set)
                    timeout : 60
                web_browser : None

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
MG5_debug end
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
param_card.dat start
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR CKMBLOCK
###################################
Block ckmblock
    1 2.277360e-01 # cabi

###################################
## INFORMATION FOR FRBLOCK
###################################
Block frblock
    1 1.166370e-05 # Gf
    2 1.570800e+00 # lam2
    3 1.570800e+00 # lam3
    4 1.000000e-01 # tth
    5 3.500000e+02 # MH5
    6 3.300000e+02 # MH3
    7 1.250000e+02 # MH
    8 4.000000e+02 # MH0
    9 5.000000e+02 # MA0
   10 4.500000e+02 # MHP

###################################
## INFORMATION FOR MASS
###################################
Block mass
    1 5.040000e-03 # MD
    2 2.550000e-03 # MU
    3 1.010000e-01 # MS
    4 1.270000e+00 # MC
    5 4.700000e+00 # MB
    6 1.720000e+02 # MT
   11 5.110000e-04 # Me
   13 1.056600e-01 # MMU
   15 1.777000e+00 # MTA
   23 9.118760e+01 # MZ
## 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.
  12 0.000000 # ve : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 79.824360 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
  252 318.802912 # hh : cmath.sqrt(M11sq + M22sq - MH__exp__2)
  25 125.000000 # h : MH
  9900035 400.000000 # h0 : MH0
  9000036 500.000000 # a0 : MA0
  9000037 450.000000 # h+ : MHP
  254 330.000000 # h30 : MH3
  253 330.000000 # h3+ : MH3
  257 350.000000 # h50 : MH5
  256 350.000000 # h5+ : MH5
  255 350.000000 # h5++ : MH5

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.279000e+02 # aEWM1
    3 1.184000e-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.720000e+02 # ymt
   11 5.110000e-04 # yme
   13 1.056600e-01 # ymm
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################

Block QNUMBERS 9900035 # h0
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000036 # a0
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000037 # 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/Antiparticle distinction (0=own anti)
Block QNUMBERS 252 # hh
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 254 # h30
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 253 # h3+
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 257 # h50
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 256 # h5+
        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/Antiparticle distinction (0=own anti)
Block QNUMBERS 255 # h5++
        1 6 # 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/Antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.466877e+00
# BR NDA ID1 ID2 ...
   1.000000e+00 2 24 5 # 1.46687686934
#
# 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.411596e+00
# BR NDA ID1 ID2 ...
   1.520727e-01 2 1 -1 # 0.366737836784
   1.520719e-01 2 3 -3 # 0.366736009828
   1.504304e-01 2 5 -5 # 0.362777280229
   1.178265e-01 2 2 -2 # 0.284149783088
   1.177064e-01 2 4 -4 # 0.283860247524
   6.878490e-02 2 16 -16 # 0.165881357446
   6.878490e-02 2 14 -14 # 0.165881357446
   6.878490e-02 2 12 -12 # 0.165881357446
   3.453867e-02 2 11 -11 # 0.0832933003408
   3.453839e-02 2 13 -13 # 0.0832926322158
   3.446033e-02 2 15 -15 # 0.0831043891058
#
# PDG Width
DECAY 24 2.002524e+00
# BR NDA ID1 ID2 ...
   3.164088e-01 2 2 -1 # 0.633616177565
   3.162879e-01 2 4 -3 # 0.633374083785
   1.111345e-01 2 12 -11 # 0.222549463313
   1.111342e-01 2 14 -13 # 0.222548878444
   1.110519e-01 2 16 -15 # 0.22238403048
   1.699462e-02 2 2 -3 # 0.0340321258784
   1.698820e-02 2 4 -1 # 0.0340192858031
#
# PDG Width
DECAY 25 6.882998e-03
# BR NDA ID1 ID2 ...
   7.804202e-01 2 5 -5 # 0.00537163044673
   5.743354e-02 2 4 -4 # 0.000395314914913
   3.745881e-02 2 15 -15 # 0.000257828885951
   1.807352e-02 3 -24 -1 2 # 0.00012440000251
   1.796020e-02 3 -2 1 24 # 0.000123619996324
   1.783380e-02 3 -4 3 24 # 0.000122749979561
   1.778731e-02 3 -24 -3 4 # 0.000122430027434
   6.335176e-03 3 -24 -11 12 # 4.36050042434e-05
   6.330672e-03 3 -14 13 24 # 4.35740006003e-05
   6.279677e-03 3 -16 15 24 # 4.32230020905e-05
   6.274301e-03 3 -12 11 24 # 4.31859977422e-05
   6.274010e-03 3 -24 -13 14 # 4.31839975071e-05
   6.269652e-03 3 -24 -15 16 # 4.31540025295e-05
   1.727736e-03 3 -1 1 23 # 1.18920042336e-05
   1.721198e-03 3 -3 3 23 # 1.18469989452e-05
   1.377205e-03 3 -5 5 23 # 9.47930192032e-06
   1.351476e-03 3 -2 2 23 # 9.30220418699e-06
   1.312568e-03 3 -4 4 23 # 9.03440348771e-06
   9.584922e-04 3 -24 -3 2 # 6.59730001651e-06
   9.568651e-04 3 -24 -1 4 # 6.58610040988e-06
   9.567489e-04 3 -2 3 24 # 6.58530031586e-06
   9.563420e-04 3 -4 1 24 # 6.58249998681e-06
   7.661875e-04 3 -12 12 23 # 5.2736667023e-06
   7.661875e-04 3 -14 14 23 # 5.2736667023e-06
   7.661875e-04 3 -16 16 23 # 5.2736667023e-06
   3.906728e-04 3 -11 11 23 # 2.68899973749e-06
   3.856459e-04 3 -13 13 23 # 2.65439994527e-06
   3.784398e-04 3 -15 15 23 # 2.60480009984e-06
   3.634697e-04 2 3 -3 # 2.50176113859e-06
   1.325944e-04 2 13 -13 # 9.12647128886e-07
   9.050826e-07 2 1 -1 # 6.22968180091e-09
   2.316901e-07 2 2 -2 # 1.59472262394e-09
   3.101328e-09 2 11 -11 # 2.13464343091e-11
#
# PDG Width
DECAY 252 2.558321e+00
# BR NDA ID1 ID2 ...
   6.825555e-01 2 25 25 # 1.74619587046
   2.191700e-01 2 -24 24 # 0.560707121926
   9.819238e-02 2 23 23 # 0.251207599077
   7.326345e-05 2 5 -5 # 0.000187431415462
   5.355800e-06 2 4 -4 # 1.37018546498e-05
   3.494869e-06 2 15 -15 # 8.94099608042e-06
   3.387667e-08 2 3 -3 # 8.66673971584e-08
   1.235828e-08 2 13 -13 # 3.16164380037e-08
   8.435665e-11 2 1 -1 # 2.15811379093e-10
   2.159428e-11 2 2 -2 # 5.52450827552e-11
   2.890539e-13 2 11 -11 # 7.39492567445e-13
#
# PDG Width
DECAY 253 5.570073e-01
# BR NDA ID1 ID2 ...
   8.167019e-01 2 25 24 # 0.454908918443
   1.832668e-01 2 6 -5 # 0.102080939327
   1.796789e-05 2 4 -3 # 1.00082440925e-05
   1.227762e-05 2 16 -15 # 6.83872417635e-06
   9.590242e-07 2 4 -1 # 5.34183423462e-07
   4.340959e-08 2 14 -13 # 2.41794569297e-08
   6.069425e-09 2 2 -3 # 3.38071377797e-09
   3.531919e-10 2 2 -1 # 1.9673046821e-10
   1.015328e-12 2 12 -11 # 5.65545067542e-13
#
# PDG Width
DECAY 254 4.134789e-01
# BR NDA ID1 ID2 ...
   9.968289e-01 2 25 23 # 0.41216773101
   1.475408e-03 3 -24 23 24 # 0.000610050055653
   6.234660e-04 3 -24 -5 6 # 0.000257790053889
   6.209265e-04 3 -6 5 24 # 0.000256740042706
   3.469866e-04 2 5 -5 # 0.000143471661394
   6.224984e-05 3 -24 24 25 # 2.57389972941e-05
   2.534474e-05 2 4 -4 # 1.04795166124e-05
   1.653947e-05 2 15 -15 # 6.8387241591e-06
   1.603009e-07 2 3 -3 # 6.62810438823e-08
   5.847809e-08 2 13 -13 # 2.41794569297e-08
   3.991668e-10 2 1 -1 # 1.65047042341e-10
   1.021818e-10 2 2 -2 # 4.22500312262e-11
   1.367772e-12 2 11 -11 # 5.65545067542e-13
#
# PDG Width
DECAY 255 1.022074e-01
# BR NDA ID1 ID2 ...
   9.999256e-01 2 24 24 # 0.102199793203
   7.438894e-05 3 -5 6 24 # 7.6031e-06
#
# PDG Width
DECAY 256 9.762803e-02
# BR NDA ID1 ID2 ...
   9.989710e-01 2 24 23 # 0.097527575747
   9.992929e-04 3 -24 24 24 # 9.75590030095e-05
   2.967590e-05 3 -5 6 23 # 2.89720000417e-06
#
# PDG Width
DECAY 257 9.641928e-02
# BR NDA ID1 ID2 ...
   6.443330e-01 2 23 23 # 0.0621261178892
   3.533173e-01 2 -24 24 # 0.0340665977342
   2.349738e-03 3 -24 23 24 # 0.00022656
#
# PDG Width
DECAY 9000036 1.303190e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 9900035 23 # 0.13031903242
#
# PDG Width
DECAY 9000037 4.517742e-04
# BR NDA ID1 ID2 ...
   3.214438e-01 3 -1 2 9900035 # 0.00014522001559
   3.140064e-01 3 -3 4 9900035 # 0.000141859990155
   1.102453e-01 3 -11 12 9900035 # 4.98059822113e-05
   1.101679e-01 3 -15 16 9900035 # 4.97710148882e-05
   1.101280e-01 3 -13 14 9900035 # 4.97529890976e-05
   1.715149e-02 3 -1 4 9900035 # 7.74860067356e-06
   1.685709e-02 3 -3 2 9900035 # 7.61559834908e-06
#
# PDG Width
DECAY 9900035 0.000000e+00

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
param_card.dat end
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Can you attach your model, such that I can try to reproduce the problem on my machine?

Thanks,

Olivier

[terryp (terrydpilkington)]

Here are both the FeynRules and UFO files.

[Olivier Mattelaer (olivier-mattelaer)]

Hi Terry,

Thanks I have succesfully reproduce the problem.
The reason is actually a wrongly formatted UFO model.

If I look at the particles H0:
H0 = Particle(pdg_code = 9900035,
              name = 'H0',
              antiname = 'H0',
              spin = 1,
              color = 1,
              mass = Param.MH0,
              width = Param.WH0,
              texname = 'H^{0}',
              antitexname = 'H^{0}',
              charge = 0,
              GhostNumber = 0,
              LeptonNumber = 0,
              Y = 0)

and If I look at the parameter MH0
MH0 = Parameter(name = 'MH0',
                nature = 'external',
                type = 'real',
                value = 400,
                texname = '\\text{MH0}',
                lhablock = 'FRBlock',
                lhacode = [ 8 ])

For Mass (same for width), it is conventional to have the lhablock set to 'MASS' (respectively 'DECAY').
Not following this convention is not a problem for MG5aMC in itself.
But can be a problem for MadSpin, Pythia and all post-processing tools.
(especially if you see the following message when running MG5aMC: WARNING: Failed to update dependent parameter

In the case of MadDM, we have to merge your model with a EFT model to be able to run the indirect detection if needed.
Such merging of model is quite tricky and can not support model where the mass (in the BLOCK MASS) is not a free parameter.
In your case MG5aMC consider it as an internal parameter (i.e. a parameter given by a formula).with the formula given by the identity operator of the eight entry of the FRBlock

If I modify your parameter MH0 to
MH0 = Parameter(name = 'MH0',
                nature = 'external',
                type = 'real',
                value = 400,
                texname = '\\text{MH0}',
                lhablock = 'MASS',
                lhacode = [ 9900035 ])

Then everything is working fine.
Fixing the merging of model to support those non standard model is (to my point of view) to risky/complicated.

Cheers,

Olivier

[terryp (terrydpilkington)]

Thanks for the help! After modifying the FeynRules model file to get rid of all the automatically-created FRBlock parameters, I can successfully do the 'generate relic_density' step, but when I attempt to output I get an error. The steps I follow are:
1. open maddm using "./maddm.py"
2. import the model using "import model GM_IDM_00a"
3. tell it which is the DM candidate "define darkmatter h0"
4. calculate relic abundance using "generate relic_density" (same error if I also have an "add direct_detection" after this step)
5. output the process "output" which results in the error:
----------
Command "output" interrupted with error:
IOError : [Errno 2] No such file or directory: '/Users/terryp/work/maddm/mg/bin/PROC_PLUGIN_GMIDM_00a__REAL_1/SubProcesses/procdef_mg5.dat'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.
----------
If I go into MG itself and do, for example, "generate h0 h0 > w+ w-; output" there is no error and the procdef_mg5.dat is written in the correct folder.

Thanks again,
Terry

===============
MG5_debug file
===============

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.0 2017-08-16 *
#* *
#* 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 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 GMIDM_00a
define darkmatter h0
import model /Users/terryp/work/maddm/mg/models/GMIDM_00a__REAL
generate relic_density
output
Traceback (most recent call last):
  File "/Users/terryp/work/maddm/mg/madgraph/interface/extended_cmd.py", line 1438, in onecmd
    return self.onecmd_orig(line, **opt)
  Fil...

[terryp (terrydpilkington)]

UFO model file

[Olivier Mattelaer (olivier-mattelaer)]

Thanks,

I will update this to fix this in maddm, either you can reinstall maddm in a couple of hours,
or use the patch in attachment.

Thanks a lot,

Olivier Mattelaer

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

We have some problem with our cluster right now, so re-installing maddm is still not working.

Cheers,

Olivier

[terryp (terrydpilkington)]

I patched the two files in the MadDM plugin folder, and everything appears to be working as intended. Thank-you very much for all your help.

[Olivier Mattelaer (olivier-mattelaer)]

Thanks to you for patience and your precious feed back.

Cheers,

Olivier
> On 13 Oct 2017, at 00:59, terryp <email address hidden> wrote:
>
> I patched the two files in the MadDM plugin folder, and everything
> appears to be working as intended. Thank-you very much for all your
> help.
>
> --
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1720194
>
> Title:
> USRMODERROR: Some inconsistency in the mass assignment in the model
>
> Status in MadDM:
> Invalid
> Status in MadGraph5_aMC@NLO:
> Won't Fix
>
> Bug description:
> Hi,
> I am trying to calculate relic abundance for a model using MadDM in MG5_aMC_v2.6.0. I receive the following error log when running ./maddm.py with my model. I have tried modifying the PDG code, but get the same error. I have included both the error log and the param card. Any help would be greatly appreciated.
>
> Thanks,
> Terry
>
>
> =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
> MG5_debug start
> =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.0 2017-08-16 *
> #* *
> #* 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 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 GMIDM_00q
> generate relic_density
> Traceback (most recent call last):
> File "/Users/terryp/work/maddm/mg/madgraph/interface/extended_cmd.py", line 1438, in onecmd
> r...