maddump

Pion not decaying

Asked by Amalia Betancur

Dear Maddump team,

I am trying to run maddump using a pi0 decay to produce DM.

This is the debug file:
Traceback (most recent call last):
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1501, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/master_interface.py", line 274, in do_import
    self.cmd.do_import(self, *args, **opts)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/madgraph_interface.py", line 5268, in do_import
    self.import_command_file(args[1])
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1643, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/PLUGIN/maddump/new_interface/maddump_interface.py", line 482, in do_launch
    self._MDUMP.exec_cmd('launch')
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 229, in do_launch
    self.run_launch()
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/interface/common_run_interface.py", line 6801, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/PLUGIN/maddump/new_interface/maddump_run_interface.py", line 370, in run_launch
    'unweighted_events.lhe.gz',interaction_channel)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/PLUGIN/maddump/meshfitter2D.py", line 903, in __init__
    self.store_reweight(input_lhe_evts)
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/PLUGIN/maddump/meshfitter2D.py", line 1048, in store_reweight
    for event in lhe_evts:
  File "/home/amalia.betancur/heptools/madump_example-master/MG5_aMC_v2_6_4/madgraph/various/lhe_parser.py", line 288, in next
    raise Exception
Exception
Value of current Options:
              text_editor : None
      notification_center : True
                    pjfry : None
       cluster_local_path : None
  default_unset_couplings : 99
       group_subprocesses : Auto
ignore_six_quark_processes : False
    loop_optimized_output : True
    cluster_status_update : (600, 30)
         fortran_compiler : None
               hepmc_path : None
                  collier : ./HEPTools/lib
              auto_update : 0
             pythia8_path : None
                hwpp_path : None
low_mem_multicore_nlo_generation : False
                    golem : None
          pythia-pgs_path : ./pythia-pgs
                  td_path : ./td
             delphes_path : ./Delphes
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : ./ExRootAnalysis
                      OLP : MadLoop
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet : None
                 run_mode : 2
              web_browser : None
   automatic_html_opening : False
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path : ./SysCalc
         madanalysis_path : ./MadAnalysis
                   lhapdf : lhapdf-config
             stdout_level : 20
                  nb_core : 4
            f2py_compiler : None
                    ninja : ./HEPTools/lib
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
           crash_on_error : False
mg5amc_py8_interface_path : None
         loop_color_flows : False
                  samurai : None
         cluster_nb_retry : 1
                  timeout : 60
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR CKMBLOCK
###################################
BLOCK CKMBLOCK #
      1 4.880000e-01 # cabi
###################################
## INFORMATION FOR GAUGEMASS
###################################
BLOCK GAUGEMASS #
      1 9.118800e+01 # mz
###################################
## INFORMATION FOR HIDDEN
###################################
BLOCK HIDDEN #
      1 1.000000e-02 # epsilon
      2 1.000000e-05 # ttheta
      3 1.000000e+01 # adm1
      4 2.000000e+02 # mhs
###################################
## INFORMATION FOR HIGGS
###################################
BLOCK HIGGS #
      1 1.250000e+02 # mh
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      1 4.700000e-03 # md
      2 2.200000e-03 # mu
      3 9.600000e-02 # ms
      4 1.420000e+00 # mc
      5 4.700000e+00 # mb
      6 1.743000e+02 # mt
      11 5.110000e-04 # me
      13 1.060000e-01 # mm
      15 1.777000e+00 # mta
      52 5.000000e+01 # md1
      58 1.000000e+02 # md2
      111 1.349800e-01 # mpi0
      211 1.395700e-01 # mpit
      1023 1.500000e+02 # map
      12 0.000000e+00 # ve : 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 8.027720e+01 # w+ : cw*mz0
      23 9.118800e+01 # z : mz
      25 1.250000e+02 # h : mh
      35 2.000000e+02 # h2 : mhs
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 2.250000e-01 # swsq
      2 1.279000e+02 # aewm1
      3 1.166390e-05 # gf
      4 1.180000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      4 1.420000e+00 # ymc
      5 4.700000e+00 # ymb
      6 1.743000e+02 # ymt
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
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 # pit
      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 52
###################################
BLOCK QNUMBERS 52 # ~chi1
      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 58
###################################
BLOCK QNUMBERS 58 # ~chi2
      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 1023
###################################
BLOCK QNUMBERS 1023 # ap
      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 35
###################################
BLOCK QNUMBERS 35 # h2
      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.508336e+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.441404e+00
#
# PDG Width
DECAY 24 2.047600e+00
#
# PDG Width
DECAY 25 2.822990e-03
#
# PDG Width
DECAY 35 5.237950e+00
#
# PDG Width
DECAY 52 0.000000e+00
#
# PDG Width
DECAY 58 1.000000e+00
#
# PDG Width
DECAY 111 0.000000e+00
#
# PDG Width
DECAY 211 0.000000e+00
#
# PDG Width
DECAY 1023 1.119100e+04
# BR NDA ID1 ID2 ...
   1.000000e+00 2 22 111 # 11191.0
   4.981235e-09 2 -2 2 # 5.5745000885e-05
   4.979716e-09 2 -4 4 # 5.5728001756e-05
   4.060048e-09 2 -11 11 # 4.5435997168e-05
   4.060048e-09 2 -13 13 # 4.5435997168e-05
   4.059244e-09 2 -15 15 # 4.5426999604e-05
   2.750603e-09 2 -1 1 # 3.0781998173e-05
   2.750603e-09 2 -3 3 # 3.0781998173e-05
   2.734519e-09 2 -5 5 # 3.0602002129e-05
   1.829446e-09 2 -12 12 # 2.0473330186e-05
   1.829446e-09 2 -14 14 # 2.0473330186e-05
   1.829446e-09 2 -16 16 # 2.0473330186e-05
   4.070682e-24 2 -52 52 # 4.5555002262e-20

thank you

Question information

Language:
English Edit question
Status:
Open
For:
maddump Edit question
Assignee:
Luca Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

Hi,

Luca is the expert concerning maddump, so I hope he will answer you.
Waiting for that, could you explain your setup?
The question here is that you are speaking of hadron (a pion) that decay into a dark matter.
Such transition is not very common in simulation tools but at the same time can be done in multiple way but you need dedicated model/....

As far as I remember they are two ways it can be done in Maddump.
1) you provide an hepmc file of the pion spectra and ask madspin to decay the pion into dark matter (and then need a dedicated model where the lorentz structure of such decay is encoded)
2) You use an internal model of pythia to generate such decay from scratch.

From the above log, I'm not sure which of those two model you are trying (if you are in one of those two modes)

Cheers,

Olivier

Revision history for this message
Amalia Betancur (amaliab2) said : 		#2

Dear Oliver,

thanks for your quick reply.

I think I am using model 1. I have provided the hepmc file and have added the pion to the BSM model I am working with following the example provided in maddump. My understanding is that I have to modify the particles.py, params.py, vertices.py, and couplings.py. Not sure what else I need to modify so that the pion can decay.

Thanks again,
amalia

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#3

Hi,

Could you attach the model (I believe attachment are going trough now but if not, you will need to create a bug report, if you want to be sure, then convert to a bug report first).

Cheers,

Olivier

Revision history for this message
Amalia Betancur (amaliab2) said : 		#4

Hi Oliver,
I created a bug report and added the model files.

Thank you,
Amalia

Can you help with this problem?

Provide an answer of your own, or ask Amalia Betancur for more information if necessary.

To post a message you must log in.