Error in GM NLO generation with MadSpin

Asked by Benjamin Freund on 2018-09-13

I am running MadGraph with Athena, the MadGraph version is 5.2.6.0.
I am trying to produce signal samples with the GM model with the following command:

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 l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model GM_UFO_nlo
generate p p > H5p > w+ z j j $$ w+ w- z [QCD] @1
add process p p > H5p~ > w- z j j $$ w+ w- z [QCD] @2

The event generation works without problem, the error occurs when trying to decay the W/Z bosons with MadSpin:

22:48:25 INFO: Running MadSpin
22:48:25 INFO: This functionality allows for the decay of resonances
22:48:25 INFO: in a .lhe file, keeping track of the spin correlation effets.
22:48:25 INFO: BE AWARE OF THE CURRENT LIMITATIONS:
22:48:25 INFO: (1) Only a succession of 2 body decay are currently allowed
22:48:25 ************************************************************
22:48:25 * *
22:48:25 * W E L C O M E to M A D S P I N *
22:48:25 * *
22:48:25 ************************************************************
22:48:26 INFO: Extracting the banner ...
22:48:28 INFO: process: p p > H5p > w+ z j j
22:48:28 INFO: options: $$ w+ w- z
22:48:29 INFO: detected model: GM_UFO_nlo. Loading...
22:48:29 set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
22:48:29 set seed 123456
22:48:29 decay t > w+ b, w+ > l+ vl
22:48:29 decay t~ > w- b~, w- > l- vl~
22:48:29 decay w+ > l+ vl
22:48:29 decay w- > l- vl~
22:48:29 decay z > l+ l-
22:48:29 launch
22:48:29 INFO: We need to recalculate the branching fractions for t~,w-,z,w+,t
22:48:29 INFO: using the FeynRules formula present in the model (arXiv:1402.1178)
22:48:29 INFO: load particles
22:48:29 INFO: load vertices
22:48:30 ^[[1;31mCommand "launch --name=GM_VBS_H5p_ew_lvl200ISHdot5_qcd0_nlo" interrupted with error:
22:48:30 KeyError : 'mdl_mh__2'
22:48:30 Please report this bug on https://bugs.launchpad.net/mg5amcnlo
22:48:30 More information is found in '/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/GM_VBS_H5p_\
ew_lvl200ISHdot5_qcd0_nlo_tag_1_debug.log'.
22:48:30 Please attach this file to your report.^[[0m
22:48:30 quit
22:48:30 INFO:
22:48:30
22:48:30 Py:MadGraphUtils ERROR MadSpin was run but can't find output folder Events/GM_VBS_H5p_ew_lvl200ISHdot5_qcd0_nlo_\
decayed_1/.
22:48:30 Shortened traceback (most recent user call last):
22:48:30 File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/AtlasProduction/19.2.5.34/InstallArea/jo\
bOptions/EvgenJobTransforms/skeleton.GENtoEVGEN.py", line 225, in <module>
22:48:30 include(jo)
22:48:30 File "./MC15.999999.MGPy8EG_A14NNPDF23LO_vbfGM_sH1_H5pWZ_lvll_m200.py", line 1, in <module>
22:48:30 include("./MadGraphPythia8EvtGenControl_GM_VBS.py")
22:48:30 File "./MadGraphPythia8EvtGenControl_GM_VBS.py", line 290, in <module>
22:48:30 isVBF
22:48:30 File "./MadGraphPythia8EvtGenControl_GM_VBS.py", line 214, in GM_Generation
22:48:30 grid_pack=False):
22:48:30 File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/AtlasProduction/19.2.5.34/InstallArea/py\
thon/MadGraphControl/MadGraphUtils.py", line 598, in generate
22:48:30 raise RuntimeError('MadSpin was run but can\'t find output folder %s.'%('Events/'+run_name+'_decayed_1/'))
22:48:30 RuntimeError: MadSpin was run but can't find output folder Events/GM_VBS_H5p_ew_lvl200ISHdot5_qcd0_nlo_decayed_1/.
22:48:30 Py:Athena INFO leaving with code 8: "an unknown exception occurred"

What is strange is that when I run a local copy of Madspin on the events.lhe I runs without problem, so I am not sure why the problem only occurs in Athena. Do you have any idea why MadSpin seems to have problems loading the Model and decaying the W and Z bosons?

Here is also the debug log with more information:
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.6.0 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch --name=GM_VBS_H5p_ew_lvl200ISHdot5_qcd0_nlo
Traceback (most recent call last):
  File "/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/bin/internal/extended_cmd.py", line 1438, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/bin/internal/amcatnlo_run_interface.py", line 1243, in do_launch
    self.exec_cmd('decay_events -from_cards', postcmd=False)
  File "/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/bin/internal/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/b/befreund/private/Athena/vbs_test/PROCNLO_GM_UFO_nlo_0/bin/internal/common_run_interface.py", line 3463, in do_decay_events
    madspin_cmd.import_command_file(path)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1580, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/various/misc.py", line 99, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/interface_madspin.py", line 579, in do_launch
    self.options)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/decay.py", line 2044, in __init__
    self.generate_all_matrix_element()
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/various/misc.py", line 99, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/various/misc.py", line 1785, in deco_f_set
    out = f(*args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/decay.py", line 2658, in generate_all_matrix_element
    self.get_branching_ratio()
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/decay.py", line 2943, in get_branching_ratio
    width.extract_br(self.list_branches, self.mgcmd)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/decay.py", line 1349, in extract_br
    self.launch_width_evaluation(to_decay, mgcmd)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/decay.py", line 1616, in launch_width_evaluation
    self.compute_widths(model, argument)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/MadSpin/decay.py", line 1654, in compute_widths
    cmd.exec_cmd(line)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/master_interface.py", line 334, in do_compute_widths
    return self.cmd.do_compute_widths(self, *args, **opts)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/madgraph_interface.py", line 7844, in do_compute_widths
    model = import_ufo.import_model(modelname, decay=True)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/models/import_ufo.py", line 127, in import_model
    model = import_full_model(model_path, decay, prefix)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/models/import_ufo.py", line 305, in import_full_model
    model.change_parameter_name_with_prefix()
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/core/base_objects.py", line 1580, in change_parameter_name_with_prefix
    part.partial_widths[key] = rep_pattern.sub(replace, value)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt/madgraph/core/base_objects.py", line 1538, in <lambda>
    replace = lambda match_pattern: change[match_pattern.groups()[0]]
KeyError: 'mdl_mh__2'
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 : /afs/cern.ch/work/b/befreund/private/Athena/vbs_test/lhapdf-config
      output_dependencies : internal
             cluster_size : 100
           cluster_memory : None
    cluster_status_update : (600, 30)
             cluster_time : None
            f2py_compiler : None
                    ninja : /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib
               hepmc_path : None
             pythia8_path : None
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : None
          pythia-pgs_path : None
                 mg5_path : /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.0.atlas/x86_64-slc6-gcc47-opt
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : None
         fortran_compiler : None
                  nb_core : 1
                  collier : /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//collier/1.1/x86_64-slc6-gcc47-opt
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : /afs/cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/syscalc/1.1.4/x86_64-slc6-gcc47-opt
                  fastjet : /afs/cern.ch/work/b/befreund/private/Athena/vbs_test/fastjet-config
             cpp_compiler : None
      notification_center : True
                 run_mode : 0
#************************************************************
#* 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
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 l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model GM_UFO_nlo
generate p p > H5p > w+ z j j $$ w+ w- z [QCD] @1
add process p p > H5p~ > w- z j j $$ w+ w- z [QCD] @2
output -f
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 8.038500e+01 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      6 1.725000e+02 # mt
      13 1.056580e-01 # mm
      15 1.776820e+00 # mta
      23 9.118760e+01 # mz
      24 8.038500e+01 # mw
      25 1.250000e+02 # mh
      252 2.793080e+02 # h
      253 2.726600e+02 # h3p
      255 2.000000e+02 # h5pp
      254 2.726600e+02 # h3z
      256 2.000000e+02 # h5p
      257 2.000000e+02 # h5z
      1 0.000000e+00 # md
      2 0.000000e+00 # mu
      3 0.000000e+00 # ms
      4 0.000000e+00 # mc
      5 0.000000e+00 # mb
      11 0.000000e+00 # me
      12 0.000000e+00 # ve
      14 0.000000e+00 # vm
      16 0.000000e+00 # vt
      21 0.000000e+00 # g
      22 0.000000e+00 # a
###################################
## INFORMATION FOR POTENTIALPARAM
###################################
BLOCK POTENTIALPARAM #
      1 8.000000e-02 # lam2
      2 -1.000000e-01 # lam3
      3 2.000000e-01 # lam4
      4 -8.663980e-01 # lam5
      5 2.889250e+02 # m1coeff
      6 4.815420e+01 # m2coeff
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.322330e+02 # aewm1
      2 1.166380e-05 # gf
      3 1.180000e-01 # as
###################################
## INFORMATION FOR VEV
###################################
BLOCK VEV #
      1 5.773500e-01 # tanth
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      1 5.040000e-03 # ymdo
      2 2.550000e-03 # ymup
      3 0.000000e+00 # yms
      4 1.643840e+00 # ymc
      5 4.579360e+00 # ymb
      6 1.725000e+02 # ymt
      11 5.110000e-04 # yme
      13 1.056580e-01 # ymm
      15 1.776820e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
BLOCK QNUMBERS 252 # h
      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)
###################################
## INFORMATION FOR QNUMBERS 253
###################################
BLOCK QNUMBERS 253 # h3p
      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)
###################################
## INFORMATION FOR QNUMBERS 254
###################################
BLOCK QNUMBERS 254 # h3z
      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)
###################################
## INFORMATION FOR QNUMBERS 255
###################################
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)
###################################
## INFORMATION FOR QNUMBERS 256
###################################
BLOCK QNUMBERS 256 # h5p
      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)
###################################
## INFORMATION FOR QNUMBERS 257
###################################
BLOCK QNUMBERS 257 # h5z
      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)
#
#*************************
# 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.508340e+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.495200e+00
#
# PDG Width
DECAY 24 2.085000e+00
#
# PDG Width
DECAY 25 7.382932e-03
# BR NDA ID1 ID2 ...
   8.192147e-01 2 5 -5 # 0.00604820617199
   4.106077e-02 2 15 -15 # 0.000303148839428
   2.103636e-02 3 -2 1 24 # 0.000155310003654
   2.103636e-02 3 -4 3 24 # 0.000155310003654
   2.091040e-02 3 -24 -1 2 # 0.000154380036989
   2.091040e-02 3 -24 -3 4 # 0.000154380036989
   7.103682e-03 3 -16 15 24 # 5.2446003942e-05
   7.053430e-03 3 -24 -11 12 # 5.20749972739e-05
   7.039614e-03 3 -12 11 24 # 5.19729949402e-05
   7.036229e-03 3 -14 13 24 # 5.19479998295e-05
   6.997897e-03 3 -24 -13 14 # 5.16649969953e-05
   6.962274e-03 3 -24 -15 16 # 5.14019964391e-05
   2.079391e-03 3 -5 5 23 # 1.53520011244e-05
   2.054875e-03 3 -1 1 23 # 1.5171000624e-05
   2.054875e-03 3 -3 3 23 # 1.5171000624e-05
   1.629230e-03 3 -2 2 23 # 1.20284972931e-05
   1.629230e-03 3 -4 4 23 # 1.20284972931e-05
   9.313373e-04 3 -12 12 23 # 6.87600021977e-06
   9.313373e-04 3 -14 14 23 # 6.87600021977e-06
   9.313373e-04 3 -16 16 23 # 6.87600021977e-06
   4.733757e-04 3 -13 13 23 # 3.49490034984e-06
   4.706802e-04 3 -11 11 23 # 3.47500019793e-06
   4.522729e-04 3 -15 15 23 # 3.33910024395e-06
#
# PDG Width
DECAY 252 1.487369e+00
# BR NDA ID1 ID2 ...
   6.358152e-01 2 -24 24 # 0.945692046732
   2.783653e-01 2 23 23 # 0.414031934031
   8.447669e-02 2 25 25 # 0.125648032896
   1.252127e-03 2 5 -5 # 0.00186237449026
   6.282013e-05 2 15 -15 # 9.34367220725e-05
   2.784783e-05 3 -24 23 24 # 4.142e-05
#
# PDG Width
DECAY 253 1.917892e+00
# BR NDA ID1 ID2 ...
   9.993819e-01 2 6 -5 # 1.91670626701
   5.199155e-04 2 25 24 # 0.000997141719398
   9.818403e-05 2 16 -15 # 0.00018830634528
   1.270052e-17 2 24 23 # 2.43582163289e-17
#
# PDG Width
DECAY 254 1.184878e-02
# BR NDA ID1 ID2 ...
   3.167177e-01 2 5 -5 # 0.00375271687669
   7.844694e-02 3 -2 1 256 # 0.000929500190708
   7.844694e-02 3 -4 3 256 # 0.000929500190708
   7.805873e-02 3 -256 -1 2 # 0.000924900321619
   7.805873e-02 3 -256 -3 4 # 0.000924900321619
   6.810798e-02 2 25 23 # 0.000806996118186
   3.008750e-02 3 -5 5 257 # 0.000356500008477
   2.981616e-02 3 -1 1 257 # 0.000353284999636
   2.981616e-02 3 -3 3 257 # 0.000353284999636
   2.624153e-02 3 -12 11 256 # 0.000310929990167
   2.618499e-02 3 -16 15 256 # 0.000310260033932
   2.567185e-02 3 -14 13 256 # 0.000304180012463
   2.346909e-02 3 -2 2 257 # 0.000278079980333
   2.346909e-02 3 -4 4 257 # 0.000278079980333
   1.589247e-02 2 15 -15 # 0.000188306344261
   1.335356e-02 3 -12 12 257 # 0.000158223354371
   1.335356e-02 3 -14 14 257 # 0.000158223354371
   1.335356e-02 3 -16 16 257 # 0.000158223354371
   1.010653e-02 3 -24 23 24 # 0.000119749972622
   6.766520e-03 3 -13 13 257 # 8.01749797428e-05
   6.665499e-03 3 -11 11 257 # 7.89780054724e-05
   6.574265e-03 3 -15 15 257 # 7.78969865312e-05
   1.340645e-03 3 -24 -5 6 # 1.58849995999e-05
#
# PDG Width
DECAY 255 2.603894e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 24 24 # 0.260389353461
#
# PDG Width
DECAY 256 2.243870e-01
# BR NDA ID1 ID2 ...
   1.000000e+00 2 24 23 # 0.224387043771
#
# PDG Width
DECAY 257 2.101888e-01
# BR NDA ID1 ID2 ...
   5.870547e-01 2 23 23 # 0.123392299937
   4.129453e-01 2 -24 24 # 0.0867964511536
#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
#***********************************************************************
#
#*******************
# 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 ! average or sum
#***********************************************************************
# 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=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
 1 = lpp1 ! beam 1 type (0 = no PDF)
 1 = lpp2 ! beam 2 type (0 = no PDF)
 6500 = ebeam1 ! beam 1 energy in GeV
 6500 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
 lhapdf = pdlabel ! PDF set
 247000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#***********************************************************************
# 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
#***********************************************************************
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses of all final state particles and partons. This *
# can be changed in SubProcesses/set_scales.f) *
#***********************************************************************
 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 = muF1_ref_fixed ! fixed fact reference scale for pdf1
 91.118 = muF2_ref_fixed ! fixed fact reference scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
#***********************************************************************
# Renormalization and factorization scales (advanced and NLO options) *
#***********************************************************************
 False = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
 91.118 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
 1.0 = muR_over_ref ! ratio of current muR over reference muR
 1.0 = muF1_over_ref ! ratio of current muF1 over reference muF1
 1.0 = muF2_over_ref ! ratio of current muF2 over reference muF2
 1.0 = QES_over_ref ! ratio of current QES over reference QES
#***********************************************************************
# Reweight flags to get scale dependence and PDF uncertainty *
# For scale dependence: factor rw_scale_up/down around central scale *
# For PDF uncertainty: use LHAPDF with supported set *
#***********************************************************************
 True = reweight_scale ! reweight to get scale dependence
 0.5 = rw_Rscale_down ! lower bound for ren scale variations
 2.0 = rw_Rscale_up ! upper bound for ren scale variations
 0.5 = rw_Fscale_down ! lower bound for fact scale variations
 2.0 = rw_Fscale_up ! upper bound for fact scale variations
 True = reweight_PDF ! reweight to get PDF uncertainty
 247000 = PDF_set_min ! First of the error PDF sets
 247101 = PDF_set_max ! Last of the error PDF sets
#***********************************************************************
# 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://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 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*Gamma) *
#***********************************************************************
 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
 0.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 gen 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
  0.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 *
#***********************************************************************
  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)
#***********************************************************************
# Maximal PDG code for quark to be considered a jet when applying cuts.*
# At least all massless quarks of the model should be included here. *
#***********************************************************************
 5 = maxjetflavor
#***********************************************************************
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
#***********************************************************************
 0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
#***********************************************************************

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2018-09-13
Last reply:
2018-09-28

Hi,

I actually fail to reproduce your error both in 2.6.0 and in the development version.
Maybe your model is not the same as the one online?

Can you check if your model is the same as the following one:
http://feynrules.irmp.ucl.ac.be/attachment/wiki/GeorgiMachacekModel/GM_UFO_nlo.tar.gz

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Benjamin Freund for more information if necessary.

To post a message you must log in.