Error related to pdg code when simulating decays

Asked by Andrew McAdams

Hi!

I have a model in which a particle Ups4s decays into b+ b-. I am trying to simulate the decay events with an external .hepmc file:

import model modelN_UFO
import_events decay Upsilon-distr.hepmc
decay ups4s > b+ b-
add process displaced_decay ups4s
output test
launch test

The structure of the file Upsilon-distr.hepmc is as follows:

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
E 6 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 1 0 0 0 0

The run proceeded...
NFO: Parsing input event to know how many decay to generate. This can takes few minuts.
INFO: Will use seed 192343602
decaying event number 100 [0.00400590896606 s]
decaying event number 200 [0.00806593894958 s]
decaying event number 300 [0.0121018886566 s]
decaying event number 400 [0.0160369873047 s]
decaying event number 500 [0.0199959278107 s]
decaying event number 600 [0.0240159034729 s]
decaying event number 700 [0.0286679267883 s]
decaying event number 800 [0.032928943634 s]
decaying event number 900 [0.037801027298 s]
decaying event number 1000 [0.0425419807434 s]
decaying event number 2000 [0.0868349075317 s]
decaying event number 3000 [0.137882947922 s]
decaying event number 4000 [0.188770055771 s]
decaying event number 5000 [0.242496013641 s]
decaying event number 6000 [0.293478012085 s]
decaying event number 7000 [0.341568946838 s]
decaying event number 8000 [0.398047924042 s]
decaying event number 9000 [0.4542760849 s]
decaying event number 10000 [0.503362894058 s]
decaying event number 20000 [0.91261100769 s]
decaying event number 30000 [1.32661390305 s]
decaying event number 40000 [1.731539011 s]
decaying event number 50000 [2.14628195763 s]
decaying event number 60000 [2.55866098404 s]
decaying event number 70000 [2.96986198425 s]
decaying event number 80000 [3.37958192825 s]
decaying event number 90000 [3.77543592453 s]
decaying event number 100000 [4.17021298409 s]

Until the following error is displayed:

Traceback (most recent call last):
  File "/home/name/Downloads/mg5_new/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/name/Downloads/mg5_new/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5_new/PLUGIN/maddump/new_interface/maddump_interface.py", line 482, in do_launch
    self._MDUMP.exec_cmd('launch')
  File "/home/name/Downloads/mg5_new/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/name/Downloads/mg5_new/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/name/Downloads/mg5_new/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 229, in do_launch
    self.run_launch()
  File "/home/name/Downloads/mg5_new/madgraph/interface/common_run_interface.py", line 7295, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/name/Downloads/mg5_new/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 338, in run_launch
    int(self.proc_characteristics['pdg_mother']))
  File "/home/name/Downloads/mg5_new/PLUGIN/maddump/new_interface/displaced_decay.py", line 29, in __init__
    for BR in self.param_card['decay'].decay_table[self.pdg_code]:
KeyError: 300553

Could you please tell me what is the problem and how to solve it?

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
Luca (lbuono) said :
#1

Dear Andrew,

from a quick look at your question and output,
I guess that the BRs are not present in the param card.
You should ask for the computation of the width, which
is automatically performed (at tree-level) in mg5.
Add/give, after launching the run, the command

set ups4s auto

Cheers,
Luca

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

Dear Luca,

excuse me, what do you mean by "after launching the run"?

Revision history for this message
Luca (lbuono) said :
#3

Hi,

I mean after the command

launch test

To be more specific, if you are providing the commands in batch-mode passing a text file, you can add it
after "launch test". If you are instead using mg5 interactively, then after "launch test",
you prompt the above command when it is asked to edit the various input cards.

Cheers,
Luca

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

So after
launch test
it asks me to edit param card or madspin card or done, and I then choose the param card, and input
set ups4s auto
somewhere in the card?

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

Nevermind, please ignore my last question. Thanks!

Revision history for this message
Luca (lbuono) said :
#6

It is not necessary to enter the param card;
instead, just give the command
set ups4s auto

the shell should prompted a message saying the parameter card
has been modified accordingly.

(I recall you that you can use the 'set' command to change all parameters appearing
in all the cards listed using the syntax
set keyword value
where 'keyword' identified the parameter and 'value' the value you want
to assign. Tab completion is enabled; an error message should be prompted if you try to assign a numerical
value to a logical parameter, for example. )

Best,
Luca

Revision history for this message
Luca (lbuono) said :
#7