maddump

A bug?

Asked by Andrew McAdams

Hi!

I have installed mg5 v 3.3, and then installed maddump.

After downloading the model DM_mesons_2 I failed to launch the standard code with decays of pions into dark photons given in the appendix D.3 of https://arxiv.org/pdf/1812.06771.pdf.

Namely, I launched

import model DM_mesons_2
import_events decay example_MesonFlux.hepmc
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

The log is

AttributeError : 'HEPMC_EventFile' object has no attribute 'zip_mode'
Please report this bug to developers

Could you please tell me how can this be? The error occurs at the very beginning of the run, at the stage of the loading of .hepmc file... I have fixed this particular error by editing the file hepmc_parser.py (please look at the debug file listed below), but then a new one occurs... This suggests that either the overall code is broken.

The debug file is

Traceback (most recent call last):
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5_te/PLUGIN/maddump/new_interface/maddump_interface.py", line 482, in do_launch
    self._MDUMP.exec_cmd('launch')
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5_te/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 229, in do_launch
    self.run_launch()
  File "/home/name/Downloads/mg5_te/madgraph/interface/common_run_interface.py", line 7310, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/name/Downloads/mg5_te/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 267, in run_launch
    madspin_cmd.import_command_file(path)
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1692, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/name/Downloads/mg5_te/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5_te/madgraph/various/misc.py", line 110, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/home/name/Downloads/mg5_te/MadSpin/interface_madspin.py", line 606, in do_launch
    return self.run_bridge(line)
  File "/home/name/Downloads/mg5_te/MadSpin/interface_madspin.py", line 845, in run_bridge
    orig_lhe = hepmc_parser.HEPMC_EventFile(filename)
  File "/home/name/Downloads/mg5_te/madgraph/various/hepmc_parser.py", line 304, in __init__
    if 'b' in mode or self.zip_mode:
AttributeError: 'HEPMC_EventFile' object has no attribute 'zip_mode'
Value of current Options:
              text_editor : None
      notification_center : True
       auto_convert_model : False
       cluster_local_path : None
        f2py_compiler_py3 : None
       group_subprocesses : Auto
ignore_six_quark_processes : False
        max_t_for_channel : 99
    loop_optimized_output : True
    cluster_status_update : (600, 30)
           crash_on_error : False
                    ninja : ./HEPTools/lib
                  collier : ./HEPTools/lib
              auto_update : 7
             pythia8_path : None
                hwpp_path : None
low_mem_multicore_nlo_generation : False
                    golem : None
          pythia-pgs_path : ./pythia-pgs
        f2py_compiler_py2 : None
             delphes_path : ./Delphes
              thepeg_path : None
             cluster_type : condor
  default_unset_couplings : 99
        madanalysis5_path : None
         cluster_nb_retry : 1
                      OLP : MadLoop
 acknowledged_v3.1_syntax : False
                  fastjet : fastjet-config
include_lepton_initiated_processes : False
               eps_viewer : None
              web_browser : None
                 pineappl : pineappl
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path : ./SysCalc
         madanalysis_path : ./MadAnalysis
                   lhapdf : lhapdf-config
             stdout_level : 20
                  nb_core : None
            f2py_compiler : None
               hepmc_path : None
       zerowidth_tchannel : True
                  td_path : ./td
   automatic_html_opening : True
       cluster_retry_wait : 300
      output_dependencies : external
      nlo_mixed_expansion : True
               lhapdf_py3 : None
               lhapdf_py2 : None
         fortran_compiler : None
mg5amc_py8_interface_path : None
         loop_color_flows : False
                  samurai : None
      exrootanalysis_path : ./ExRootAnalysis
                  timeout : 60
                    gauge : unitary
                 run_mode : 2
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
######################################################################
## 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(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
      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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)

Question information

Language:
English Edit question
Status:
Open
For:
maddump Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Andrew McAdams (mcadams1994) said : 		#1

Please also find the log of the run:

load MG5 configuration from ../input/mg5_configuration.txt
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/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

None does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
set automatic_html_opening False --no_save
/home/name/Downloads/mg5_te/models/DM_mesons_2
import model /home/name/Downloads/mg5_te/models/DM_mesons_2
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=DP_electron/Cards/param_card.dat
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.Automatically run the model.find_vertexlist()
Found 8 stable particles
Pass to numerical integration for computing the widths:
INFO: More info in temporary files:
    /home/name/Downloads/mg5_te/bin/tmpl6auW8/temp_decay/index.html
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: 01h54 ]
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
INFO: Idle: 0, Running: 3, Completed: 9 [ 4.1s ]
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.1s ]
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/Cards/param_card.dat
INFO: Update the dependent parameter of the param_card.dat
/home/name/Downloads/mg5_te/models/DM_mesons_2
import model /home/name/Downloads/mg5_te/models/DM_mesons_2
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/Events_to_decay/example_MesonFlux.hepmc
INFO: Setup the code for pure decay mode
import model /home/name/Downloads/mg5_te/models/DM_mesons_2 DP_electron/Cards/param_card.dat
INFO: detected model: /home/name/Downloads/mg5_te/models/DM_mesons_2. Loading...
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
AttributeError : 'HEPMC_EventFile' object has no attribute 'zip_mode'
Please report this bug to developers

           More information is found in 'debug'.

           Please attach this file to your report.

Can you help with this problem?

Provide an answer of your own, or ask Andrew McAdams for more information if necessary.

To post a message you must log in.