maddump

A problem with maddump

Asked by Andrew McAdams

Hi! I have the following question. I simulate the decay
ups4s > b+ b-
using my customary model.

Using pure mg5, I write
generate ups4s > b+ b-
output test2
launch test2

Everything goes fine.

When launching the same process in maddump, I do the following:

import_events Upsilon-distr.hempc
decay ups4s > b+ b-
add process displaced_decay ups4s
output test
launch test

where Upsilon-distr.hempc has the form

HepMC::Version 2.06.09
HepMC::IO_GenEvent-START_EVENT_LISTING
E 0 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0
U GEV MM
V -1 0 0 0 0 0 0 1 0
P 1 300553 0.0000000000000000e+00 0.0000000000000000e+00 3.0000000000000000e+00 1.099710871092943e+01 1.02800000000000000e+01 0 0 0 0 0
E 1 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0
U GEV MM
V -1 0 0 0 0 0 0 1 0
P 1 300553 0.0000000000000000e+00 0.0000000000000000e+00 3.0000000000000000e+00 1.099710871092943e+01 1.02800000000000000e+01 0 0 0 0 0
E 2 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0
U GEV MM
V -1 0 0 0 0 0 0 1 0
P 1 300553 0.0000000000000000e+00 0.0000000000000000e+00 3.0000000000000000e+00 1.099710871092943e+01 1.02800000000000000e+01 0 0 0 0 0
E 3 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0
U GEV MM
V -1 0 0 0 0 0 0 1 0
P 1 300553 0.0000000000000000e+00 0.0000000000000000e+00 3.0000000000000000e+00 1.099710871092943e+01 1.02800000000000000e+01 0 0 0 0 0
E 4 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0
U GEV MM
V -1 0 0 0 0 0 0 1 0
P 1 300553 0.0000000000000000e+00 0.0000000000000000e+00 3.0000000000000000e+00 1.099710871092943e+01 1.02800000000000000e+01 0 0 0 0 0
E 5 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0
U GEV MM
V -1 0 0 0 0 0 0 1 0
P 1 300553 0.0000000000000000e+00 0.0000000000000000e+00 3.0000000000000000e+00 1.099710871092943e+01 1.02800000000000000e+01 0 0 0 0 0
...

The weird thing is that it does not properly generate the information in the output folder.

The output is
InvalidCmd : No such file or directory : test/Events_to_decay/unweighted_events.lhe.gz

What may be wrong with these commands?

Question information

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

The debug file:

Traceback (most recent call last):
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1541, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5/PLUGIN/maddump/new_interface/maddump_interface.py", line 482, in do_launch
    self._MDUMP.exec_cmd('launch')
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 229, in do_launch
    self.run_launch()
  File "/home/name/Downloads/mg5/madgraph/interface/common_run_interface.py", line 7111, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/name/Downloads/mg5/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 267, in run_launch
    madspin_cmd.import_command_file(path)
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1689, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/name/Downloads/mg5/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5/MadSpin/interface_madspin.py", line 198, in do_import
    raise self.InvalidCmd('No such file or directory : %s' % inputfile)
InvalidCmd: No such file or directory : test/Events_to_decay/unweighted_events.lhe.gz
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 : /home/name/Downloads/mg5/HEPTools/lib
                  collier : /home/name/Downloads/mg5/HEPTools/lib
              auto_update : 7
             pythia8_path : ./HEPTools/pythia8
                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
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet : fastjet-config
                 run_mode : 2
              web_browser : None
   automatic_html_opening : True
        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
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
               lhapdf_py3 : None
               lhapdf_py2 : None
         fortran_compiler : None
mg5amc_py8_interface_path : /home/name/Downloads/mg5/HEPTools/MG5aMC_PY8_interface
         loop_color_flows : False
                  samurai : None
      exrootanalysis_path : ./ExRootAnalysis
                  timeout : 60
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
######################################################################
## 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
    2 3.508000e-03 # th13
    3 4.153900e-02 # th23
    4 1.200000e+00 # del13

###################################
## 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
   25 1.251000e+02 # MH
  111 1.349800e-01 # MPi0
  113 7.752600e-01 # Mrho0
  211 1.395700e-01 # MPip
  213 7.751100e-01 # Mrho
  221 5.478600e-01 # Meta
  223 7.826500e-01 # Mom
  311 4.976110e-01 # MK0
  321 4.936770e-01 # MK
  323 8.917000e-01 # MKstar
  331 9.578000e-01 # Metap
  333 1.019461e+00 # Mphi
  411 1.869000e+00 # MDd
  421 1.865000e+00 # MD0
  431 1.968300e+00 # MDs
  511 5.280000e+00 # MB0
  521 5.279000e+00 # MBu
  523 5.324700e+00 # MBstar
  531 5.367000e+00 # MB0s
  535 5.839860e+00 # MBs2Star
  541 6.275000e+00 # MBc
  300553 1.058000e+01 # MUps4s
## 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.
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  24 7.983323e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt((-2.220315292925759*aEW*MZ__exp__2)/Gf + MZ__exp__4/4.))
  12000 0.000000e+00 # v1 : 0.0
  14000 0.000000e+00 # v2 : 0.0
  16000 0.000000e+00 # v3 : 0.0
  18000 1.000000e+00 # n4 : MN4

###################################
## INFORMATION FOR MESONBLOCK
###################################
Block mesonblock
    2 1.556000e-01 # fK
    3 1.302000e-01 # fpi
    4 8.159000e-02 # feta
    5 -9.460000e-02 # fetap
    6 2.120000e-01 # fD
    7 1.700000e-01 # frho
    8 2.490000e-01 # fDs
    9 1.550000e-01 # fomega
   10 2.320000e-01 # fphi
   11 1.900000e-01 # fB
   12 2.300000e-01 # fBs
   13 1.900000e-01 # fBc
   14 1.900000e-01 # fBs2BK
   16 1.900000e-01 # fD0Kpi
   17 1.900000e-01 # frho0pipi
   18 1.900000e-01 # fUps4sBB
   34 1.770000e-01 # fKstar

###################################
## 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
    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
   16 1.000000e-05 # modthetae
   17 0.000000e+00 # argthetae
   18 1.000000e-05 # modthetamu
   19 0.000000e+00 # argthetamu
   20 1.000000e-05 # modthetatau
   21 0.000000e+00 # argthetatau
   22 1.000000e+00 # MN4

###################################
## 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.478000e-01 # Wrho0
DECAY 213 1.491000e-01 # Wrho
DECAY 223 8.490000e-03 # Wom
DECAY 323 5.080000e-02 # WKstar
DECAY 333 4.249000e-03 # Wphi
DECAY 411 6.320000e-13 # WD
DECAY 421 1.620000e-12 # WD0
DECAY 431 1.300000e-12 # WDs
DECAY 511 4.330000e-13 # WB0
DECAY 521 4.000000e-13 # WBu
DECAY 523 1.000000e-13 # WBstar
DECAY 531 4.320000e-13 # WB0s
DECAY 535 1.490000e-03 # WBs2Star
DECAY 300553 2.050000e-03 # WUps4s
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.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 13 0.000000e+00 # mu- : 0.0
DECAY 15 0.000000e+00 # ta- : 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 311 0.000000e+00 # k0 : 0.0
DECAY 321 0.000000e+00 # k+ : 0.0
DECAY 331 0.000000e+00 # etap : 0.0
DECAY 541 0.000000e+00 # bc+ : 0.0
DECAY 12000 0.000000e+00 # v1 : 0.0
DECAY 14000 0.000000e+00 # v2 : 0.0
DECAY 16000 0.000000e+00 # v3 : 0.0
DECAY 18000 0.000000e+00 # n4 : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 12000 # v1
        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 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 14000 # v2
        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 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 16000 # v3
        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 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 18000 # n4
        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 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 321 # k+
        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 211 # pi+
        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 411 # d+
        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 421 # d0
        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)
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)
Block QNUMBERS 311 # k0
        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)
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)
Block QNUMBERS 331 # etap
        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 431 # ds+
        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 521 # b+
        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 541 # bc+
        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 511 # b0
        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)
Block QNUMBERS 531 # b0s
        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)
Block QNUMBERS 535 # bs2
        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)
Block QNUMBERS 113 # rho0
        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)
Block QNUMBERS 213 # rho+
        1 3 # 3 times electric charge
        2 3 # 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 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)
Block QNUMBERS 323 # kstar+
        1 3 # 3 times electric charge
        2 3 # 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 333 # phimeson
        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)
Block QNUMBERS 523 # bstar+
        1 3 # 3 times electric charge
        2 3 # 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 300553 # ups4s
        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)

Revision history for this message
Andrew McAdams (mcadams1994) said : 		#2

I get the same error when using the example code for the production of dark matter particle in decays of dark photon:

import model DM_mesons_2
import_events decay path/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

And then

Revision history for this message
Andrew McAdams (mcadams1994) said : 		#3

Excuse me, there was an explicit error in my run command.