MadGraph5_aMC@NLO

SLHA file input as param_card.dat not working

Asked by Akshat Khanna

I am using MG5_aMC_v3_5_0. I wrote the model file on SARAH, and generated a UFO file from the same. I tried to use the the SPheno output file "SPheno.spc.$ModelName", which is an SLHA file, as my param_card.dat. But I get an error, whenever I launch a process -

Command "generate_events run_01" interrupted with error:
InvalidParam :

The first few lines of the ME5_debug file reads like this -

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.0 2023-05-12 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/models/check_param_card.py", line 82, in load_str
    self.lhacode = tuple([int(d) for d in data[:-1]])
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/models/check_param_card.py", line 82, in <listcomp>
    self.lhacode = tuple([int(d) for d in data[:-1]])
ValueError: invalid literal for int() with base 10: 'l'

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/extended_cmd.py", line 1543, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
    return func(arg, **opt)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/madevent_interface.py", line 2395, in do_generate_events
    switch_mode = self.ask_run_configuration(mode, args)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/madevent_interface.py", line 6441, in ask_run_configuration
    self.ask_edit_cards(cards, plot=False, first_cmd=first_cmd)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/common_run_interface.py", line 1039, in ask_edit_cards
    self.ask_edit_card_static(cards, mode, plot, self.options['timeout'],
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/common_run_interface.py", line 1111, in ask_edit_card_static
    out = ask(question, '0', possible_answer, timeout=int(1.5*timeout),
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/extended_cmd.py", line 1111, in ask
    question_instance = obj(question, allow_arg=choices, default=default,
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/common_run_interface.py", line 4925, in __init__
    new_vars = set(getattr(self, 'init_%s' % name)(cards))
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/madgraph/interface/common_run_interface.py", line 5011, in init_param
    self.param_card = param_card_mod.ParamCard(self.paths['param'])
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/models/check_param_card.py", line 339, in __init__
    self.read(input_path)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/models/check_param_card.py", line 373, in read
    param.load_str(line[6:])
  File "/home/iitgn/Applications/MG5_aMC_v3_5_0/models/check_param_card.py", line 88, in load_str
    raise InvalidParam
models.check_param_card.InvalidParam
                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : True
        notification_center : True
                   run_mode : 2
              cluster_queue : None (user set)
               cluster_time : None (user set)
               cluster_size : 100
             cluster_memory : 100 (user set)
                    nb_core : 8 (user set)
          cluster_temp_path : None

                      Configuration Options
                      ---------------------
               pythia8_path : /home/iitgn/Applications/MG5_aMC_v3_5_0/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
                thepeg_path : None (user set)
                 hepmc_path : None (user set)
           madanalysis_path : None (user set)
          madanalysis5_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /home/iitgn/Applications/MG5_aMC_v3_5_0/Delphes (user set)
        exrootanalysis_path : None (user set)
               syscalc_path : None (user set)
                 rivet_path : None
                  yoda_path : None
                     lhapdf : /home/iitgn/Applications/MG5_aMC_v3_5_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                 lhapdf_py2 : None
                 lhapdf_py3 : /home/iitgn/Applications/MG5_aMC_v3_5_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                    timeout : 60
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                web_browser : None
                 eps_viewer : None
                text_editor : None
           fortran_compiler : None
               cpp_compiler : None
                auto_update : 7 (user set)
               cluster_type : condor
      cluster_status_update : (600, 30)
           cluster_nb_retry : 1
         cluster_local_path : None
         cluster_retry_wait : 300
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.0 2023-05-12 *
#* *
#* 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 low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
import model sm
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 vl = ve vm vt
define vl~ = ve~ vm~ vt~
convert model /home/iitgn/Applications/MG5_aMC_v3_5_0/models/3HDM-UFO/
import model /home/iitgn/Applications/MG5_aMC_v3_5_0/models/3HDM-UFO
generate hm2c > u3 d3bar
output h2bbart
# SUSY Les Houches Accord 2 - 3HDM-Z Spectrum + Decays + Flavour Observables
# SPheno module generated by SARAH
# ----------------------------------------------------------------------
# SPheno v4.0.5
# W. Porod, Comput. Phys. Commun. 153 (2003) 275-315, hep-ph/0301101
# W. Porod, F.Staub, Comput.Phys.Commun.183 (2012) 2458-2469, arXiv:1104.1573
# SARAH: 4.15.1
# F. Staub; arXiv:0806.0538 (online manual)
# F. Staub; Comput. Phys. Commun. 181 (2010) 1077-1086; arXiv:0909.2863
# F. Staub; Comput. Phys. Commun. 182 (2011) 808-833; arXiv:1002.0840
# F. Staub; Comput. Phys. Commun. 184 (2013) 1792-1809; arXiv:1207.0906
# F. Staub; Comput. Phys. Commun. 185 (2014) 1773-1790; arXiv:1309.7223
# Including the calculation of flavor observables based on FlavorKit
# W. Porod, F. Staub, A. Vicente; Eur.Phys.J. C74 (2014) 8, 2992; arXiv:1405.1434
# Two-loop mass corrections to Higgs fields based on
# M. D. Goodsell, K. Nickel, F. Staub; Eur.Phys.J. C75 (2015) no.6, 290; arXiv:1411.0675
# M. D. Goodsell, K. Nickel, F. Staub; Eur.Phys.J. C75 (2015) no.1, 32; arXiv:1503.03098
# M. D. Goodsell, F. Staub; Eur.Phys.J.C 77 (2017) 1, 46, arXiv:1604.05335
# J. Braathen, M. D. Goodsell, F. Staub; Eur.Phys.J.C 77 (2017) 11, 758; arXiv:1703.09237
# Loop decays based on
# M. D. Goodsell, S. Liebler, F. Staub; Eur.Phys.J.C 77 (2017) 11, 757; arXiv:1706.05372
# Unitarity routines based on
# M. D. Goodsell, F. Staub; Eur.Phys.J.C 78 (2018) 8, 649; arXiv:1805.07306
# M. D. Goodsell, R. Moutafis; arXiv:2012.09022
#
# in case of problems send email to <email address hidden> and <email address hidden>, or consult the forum/wiki
# ----------------------------------------------------------------------
# Created: 10.07.2023, 12:37
Block SPINFO # Program information
     1 SPhenoSARAH # spectrum calculator
     2 v4.0.5 # version number of SPheno
     9 4.15.1 # version number of SARAH
Block MODSEL # Input parameters
     1 1 # GUT scale input
     2 1 # Boundary conditions
     6 1 # switching on flavour violation
Block MINPAR # Input parameters
    1 1.25000000E+02 # mH1input
    2 3.98391490E+02 # mH2input
    3 4.29147091E+02 # mH3input
    4 5.78685513E+02 # mA1input
    5 2.56769225E+02 # mA2input
    6 4.31861223E+02 # mC1input
    7 2.41032063E+02 # mC2input
    8 9.38309883E-01 # tb1
    9 1.01344177E+00 # tb2
   10 1.73235228E+01 # tg1
   11 1.51774291E+01 # tg2
   12 9.38309883E-01 # ta1
   13 1.01344177E+00 # ta2
   14 2.17681397E+01 # ta4
.....

It would be great if you could please address this issue.

Thank you

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

HI,

I would need at least the full content of the param_card in order to be able to see if the content of that file makes sense.
From the part that you have included, I could see that the following block:
Block SPINFO # Program information
     1 SPhenoSARAH # spectrum calculator
     2 v4.0.5 # version number of SPheno
     9 4.15.1 # version number of SARAH
does not follow the expected convention, but this block is likely not included in the UFO format and therefore it might be ignored (or not).

From the error reported this is not that block responsible for the crash so it is likely further away in your card where the (first?) issue is.

Cheers,

Olivier

Revision history for this message
Akshat Khanna (khanna-akshat) said : 		#2

Hello,

Thanks for your message. I am hereby adding the remaining content of the ME5_debug, which also consist of the param_card details. The content of the param_card.dat is just the SPheno output file,

convert model /home/iitgn/Applications/MG5_aMC_v3_5_0/models/3HDM-UFO/
import model /home/iitgn/Applications/MG5_aMC_v3_5_0/models/3HDM-UFO
generate hm2c > u3 d3bar
output h2bbart
# SUSY Les Houches Accord 2 - 3HDM-Z Spectrum + Decays + Flavour Observables
# SPheno module generated by SARAH
# ----------------------------------------------------------------------
# SPheno v4.0.5
# W. Porod, Comput. Phys. Commun. 153 (2003) 275-315, hep-ph/0301101
# W. Porod, F.Staub, Comput.Phys.Commun.183 (2012) 2458-2469, arXiv:1104.1573
# SARAH: 4.15.1
# F. Staub; arXiv:0806.0538 (online manual)
# F. Staub; Comput. Phys. Commun. 181 (2010) 1077-1086; arXiv:0909.2863
# F. Staub; Comput. Phys. Commun. 182 (2011) 808-833; arXiv:1002.0840
# F. Staub; Comput. Phys. Commun. 184 (2013) 1792-1809; arXiv:1207.0906
# F. Staub; Comput. Phys. Commun. 185 (2014) 1773-1790; arXiv:1309.7223
# Including the calculation of flavor observables based on FlavorKit
# W. Porod, F. Staub, A. Vicente; Eur.Phys.J. C74 (2014) 8, 2992; arXiv:1405.1434
# Two-loop mass corrections to Higgs fields based on
# M. D. Goodsell, K. Nickel, F. Staub; Eur.Phys.J. C75 (2015) no.6, 290; arXiv:1411.0675
# M. D. Goodsell, K. Nickel, F. Staub; Eur.Phys.J. C75 (2015) no.1, 32; arXiv:1503.03098
# M. D. Goodsell, F. Staub; Eur.Phys.J.C 77 (2017) 1, 46, arXiv:1604.05335
# J. Braathen, M. D. Goodsell, F. Staub; Eur.Phys.J.C 77 (2017) 11, 758; arXiv:1703.09237
# Loop decays based on
# M. D. Goodsell, S. Liebler, F. Staub; Eur.Phys.J.C 77 (2017) 11, 757; arXiv:1706.05372
# Unitarity routines based on
# M. D. Goodsell, F. Staub; Eur.Phys.J.C 78 (2018) 8, 649; arXiv:1805.07306
# M. D. Goodsell, R. Moutafis; arXiv:2012.09022
#
# in case of problems send email to <email address hidden> and <email address hidden>, or consult the forum/wiki
# ----------------------------------------------------------------------
# Created: 10.07.2023, 12:37
Block SPINFO # Program information
     1 SPhenoSARAH # spectrum calculator
     2 v4.0.5 # version number of SPheno
     9 4.15.1 # version number of SARAH
Block MODSEL # Input parameters
     1 1 # GUT scale input
     2 1 # Boundary conditions
     6 1 # switching on flavour violation
Block MINPAR # Input parameters
    1 1.25000000E+02 # mH1input
    2 3.98391490E+02 # mH2input
    3 4.29147091E+02 # mH3input
    4 5.78685513E+02 # mA1input
    5 2.56769225E+02 # mA2input
    6 4.31861223E+02 # mC1input
    7 2.41032063E+02 # mC2input
    8 9.38309883E-01 # tb1
    9 1.01344177E+00 # tb2
   10 1.73235228E+01 # tg1
   11 1.51774291E+01 # tg2
   12 9.38309883E-01 # ta1
   13 1.01344177E+00 # ta2
   14 2.17681397E+01 # ta4
Block gaugeGUT Q= -1.00000000E+00 # (GUT scale)
   1 0.00000000E+00 # g1(Q)
   2 0.00000000E+00 # g2(Q)
   3 0.00000000E+00 # g3(Q)
Block SMINPUTS # SM parameters
         1 1.27932000E+02 # alpha_em^-1(MZ)^MSbar
         2 1.16637000E-05 # G_mu [GeV^-2]
         3 1.18700000E-01 # alpha_s(MZ)^MSbar
         4 9.11887000E+01 # m_Z(pole)
         5 4.18000000E+00 # m_b(m_b), MSbar
         6 1.72890000E+02 # m_t(pole)
         7 1.77669000E+00 # m_tau(pole)
Block GAUGE Q= 1.60000000E+02 # (Renormalization Scale)
   1 3.55690247E-01 # g1
   2 6.62798598E-01 # g2
   3 1.22132231E+00 # g3
Block MixedBlock Q= 1.60000000E+02 # (Renormalization Scale)
   1 -6.56303023E-01 # Lam10
   2 -7.14345827E-01 # Lam11
   3 -5.28157910E+00 # Lam12
   4 1.54639165E+00 # Lam1
   5 -1.98907571E-01 # Lam4
   6 3.65860240E+00 # Lam5
   7 4.39416903E+00 # Lam7
   8 -3.69234895E-01 # Lam8
   9 1.63654801E+00 # Lam2
  10 4.41220759E+00 # Lam6
  11 -4.07687340E-01 # Lam9
  12 1.28786743E+00 # Lam3
  13 -7.81250000E+03 # M112
  14 -7.81250000E+03 # M222
  15 -7.81250000E+03 # M332
Block HMIX Q= 1.60000000E+02 # (Renormalization Scale)
 102 1.26113773E+02 # v1
 103 1.18333799E+02 # v2
 104 1.75262651E+02 # v3
 102 1.26113773E+02 # v1
 103 1.18333799E+02 # v2
 104 1.75262651E+02 # v3
   3 2.46220569E+02 # v
  10 -9.47868674E-01 # beta1
  11 4.99273160E-01 # beta2
  12 -6.45664934E-01 # gamma1
  13 -6.49652651E-01 # gamma2
  14 -8.20914169E-01 # alpha1
  15 7.20437654E-01 # alpha2
  16 -1.52789393E+00 # alpha4
Block Yu Q= 1.60000000E+02 # (Renormalization Scale)
  1 1 1.96537724E-05 # Real(Yu(1,1),dp)
  1 2 4.54590803E-06 # Real(Yu(1,2),dp)
  1 3 6.90916244E-08 # Real(Yu(1,3),dp)
  2 1 -2.30878381E-03 # Real(Yu(2,1),dp)
  2 2 9.97537813E-03 # Real(Yu(2,2),dp)
  2 3 4.22048901E-04 # Real(Yu(2,3),dp)
  3 1 8.29634516E-03 # Real(Yu(3,1),dp)
  3 2 -5.70534366E-02 # Real(Yu(3,2),dp)
  3 3 1.39387656E+00 # Real(Yu(3,3),dp)
Block IMYu Q= 1.60000000E+02 # (Renormalization Scale)
  1 1 0.00000000E+00 # Aimag(Yu(1,1))
  1 2 0.00000000E+00 # Aimag(Yu(1,2))
  1 3 0.00000000E+00 # Aimag(Yu(1,3))
  2 1 0.00000000E+00 # Aimag(Yu(2,1))
  2 2 0.00000000E+00 # Aimag(Yu(2,2))
  2 3 0.00000000E+00 # Aimag(Yu(2,3))
  3 1 0.00000000E+00 # Aimag(Yu(3,1))
  3 2 0.00000000E+00 # Aimag(Yu(3,2))
  3 3 0.00000000E+00 # Aimag(Yu(3,3))
Block Ye Q= 1.60000000E+02 # (Renormalization Scale)
  1 1 5.73023549E-06 # Real(Ye(1,1),dp)
  1 2 0.00000000E+00 # Real(Ye(1,2),dp)
  1 3 0.00000000E+00 # Real(Ye(1,3),dp)
  2 1 0.00000000E+00 # Real(Ye(2,1),dp)
  2 2 1.18483095E-03 # Real(Ye(2,2),dp)
  2 3 0.00000000E+00 # Real(Ye(2,3),dp)
  3 1 0.00000000E+00 # Real(Ye(3,1),dp)
  3 2 0.00000000E+00 # Real(Ye(3,2),dp)
  3 3 1.99234313E-02 # Real(Ye(3,3),dp)
Block IMYe Q= 1.60000000E+02 # (Renormalization Scale)
  1 1 0.00000000E+00 # Aimag(Ye(1,1))
  1 2 0.00000000E+00 # Aimag(Ye(1,2))
  1 3 0.00000000E+00 # Aimag(Ye(1,3))
  2 1 0.00000000E+00 # Aimag(Ye(2,1))
  2 2 0.00000000E+00 # Aimag(Ye(2,2))
  2 3 0.00000000E+00 # Aimag(Ye(2,3))
  3 1 0.00000000E+00 # Aimag(Ye(3,1))
  3 2 0.00000000E+00 # Aimag(Ye(3,2))
  3 3 0.00000000E+00 # Aimag(Ye(3,3))
Block Yd Q= 1.60000000E+02 # (Renormalization Scale)
  1 1 5.97552672E-05 # Real(Yd(1,1),dp)
  1 2 0.00000000E+00 # Real(Yd(1,2),dp)
  1 3 0.00000000E+00 # Real(Yd(1,3),dp)
  2 1 0.00000000E+00 # Real(Yd(2,1),dp)
  2 2 1.13535008E-03 # Real(Yd(2,2),dp)
  2 3 0.00000000E+00 # Real(Yd(2,3),dp)
  3 1 0.00000000E+00 # Real(Yd(3,1),dp)
  3 2 0.00000000E+00 # Real(Yd(3,2),dp)
  3 3 4.99554034E-02 # Real(Yd(3,3),dp)
Block IMYd Q= 1.60000000E+02 # (Renormalization Scale)
  1 1 0.00000000E+00 # Aimag(Yd(1,1))
  1 2 0.00000000E+00 # Aimag(Yd(1,2))
  1 3 0.00000000E+00 # Aimag(Yd(1,3))
  2 1 0.00000000E+00 # Aimag(Yd(2,1))
  2 2 0.00000000E+00 # Aimag(Yd(2,2))
  2 3 0.00000000E+00 # Aimag(Yd(2,3))
  3 1 0.00000000E+00 # Aimag(Yd(3,1))
  3 2 0.00000000E+00 # Aimag(Yd(3,2))
  3 3 0.00000000E+00 # Aimag(Yd(3,3))
Block MASS # Mass spectrum
# PDG code mass particle
        25 1.25000000E+02 # hh_1
        35 3.98391490E+02 # hh_2
        45 4.29147091E+02 # hh_3
        36 2.56769225E+02 # Ah_2
       257 5.78685513E+02 # Ah_3
        37 2.41032063E+02 # Hm_2
       256 4.31861223E+02 # Hm_3
        22 0.00000000E+00 # VP
        23 9.11887000E+01 # VZ
        21 0.00000000E+00 # VG
        24 8.03497269E+01 # VWm
         1 5.00000000E-03 # Fd_1
         3 9.50000000E-02 # Fd_2
         5 4.18000000E+00 # Fd_3
         2 2.50000000E-03 # Fu_1
         4 1.27000000E+00 # Fu_2
         6 1.72890000E+02 # Fu_3
        11 5.10998930E-04 # Fe_1
        13 1.05658372E-01 # Fe_2
        15 1.77669000E+00 # Fe_3
        12 0.00000000E+00 # Fv_1
        14 0.00000000E+00 # Fv_2
        16 0.00000000E+00 # Fv_3
Block SCALARMIX Q= 1.60000000E+02 # ()
  1 1 5.12198365E-01 # ZH(1,1)
  1 2 4.80600788E-01 # ZH(1,2)
  1 3 7.11811575E-01 # ZH(1,3)
  2 1 -5.49937041E-01 # ZH(2,1)
  2 2 -4.53082631E-01 # ZH(2,2)
  2 3 7.01630516E-01 # ZH(2,3)
  3 1 6.59713639E-01 # ZH(3,1)
  3 2 -7.50825554E-01 # ZH(3,2)
  3 3 3.22319925E-02 # ZH(3,3)
Block PSEUDOSCALARMIX Q= 1.60000000E+02 # ()
  1 1 5.12198365E-01 # ZA(1,1)
  1 2 4.80600788E-01 # ZA(1,2)
  1 3 7.11811575E-01 # ZA(1,3)
  2 1 -7.13032480E-01 # ZA(2,1)
  2 2 6.99961642E-01 # ZA(2,2)
  2 3 4.04769391E-02 # ZA(2,3)
  3 1 -4.78787550E-01 # ZA(3,1)
  3 2 -5.28276994E-01 # ZA(3,2)
  3 3 7.01203180E-01 # ZA(3,3)
Block CHARGEMIX Q= 1.60000000E+02 # ()
  1 1 5.12198365E-01 # ZP(1,1)
  1 2 4.80600788E-01 # ZP(1,2)
  1 3 7.11811575E-01 # ZP(1,3)
  2 1 -7.16901995E-01 # ZP(2,1)
  2 2 6.95643010E-01 # ZP(2,2)
  2 3 4.61771802E-02 # ZP(2,3)
  3 1 -4.72973957E-01 # ZP(3,1)
  3 2 -5.33951014E-01 # ZP(3,2)
  3 3 7.00850876E-01 # ZP(3,3)
Block UDLMIX Q= 1.60000000E+02 # ()
  1 1 1.00000000E+00 # Real(ZDL(1,1),dp)
  1 2 0.00000000E+00 # Real(ZDL(1,2),dp)
  1 3 0.00000000E+00 # Real(ZDL(1,3),dp)
  2 1 0.00000000E+00 # Real(ZDL(2,1),dp)
  2 2 1.00000000E+00 # Real(ZDL(2,2),dp)
  2 3 0.00000000E+00 # Real(ZDL(2,3),dp)
  3 1 0.00000000E+00 # Real(ZDL(3,1),dp)
  3 2 0.00000000E+00 # Real(ZDL(3,2),dp)
  3 3 1.00000000E+00 # Real(ZDL(3,3),dp)
Block IMUDLMIX Q= 1.60000000E+02 # ()
  1 1 0.00000000E+00 # Aimag(ZDL(1,1))
  1 2 0.00000000E+00 # Aimag(ZDL(1,2))
  1 3 0.00000000E+00 # Aimag(ZDL(1,3))
  2 1 0.00000000E+00 # Aimag(ZDL(2,1))
  2 2 0.00000000E+00 # Aimag(ZDL(2,2))
  2 3 0.00000000E+00 # Aimag(ZDL(2,3))
  3 1 0.00000000E+00 # Aimag(ZDL(3,1))
  3 2 0.00000000E+00 # Aimag(ZDL(3,2))
  3 3 0.00000000E+00 # Aimag(ZDL(3,3))
Block UDRMIX Q= 1.60000000E+02 # ()
  1 1 1.00000000E+00 # Real(ZDR(1,1),dp)
  1 2 0.00000000E+00 # Real(ZDR(1,2),dp)
  1 3 0.00000000E+00 # Real(ZDR(1,3),dp)
  2 1 0.00000000E+00 # Real(ZDR(2,1),dp)
  2 2 1.00000000E+00 # Real(ZDR(2,2),dp)
  2 3 0.00000000E+00 # Real(ZDR(2,3),dp)
  3 1 0.00000000E+00 # Real(ZDR(3,1),dp)
  3 2 0.00000000E+00 # Real(ZDR(3,2),dp)
  3 3 1.00000000E+00 # Real(ZDR(3,3),dp)
Block IMUDRMIX Q= 1.60000000E+02 # ()
  1 1 0.00000000E+00 # Aimag(ZDR(1,1))
  1 2 0.00000000E+00 # Aimag(ZDR(1,2))
  1 3 0.00000000E+00 # Aimag(ZDR(1,3))
  2 1 0.00000000E+00 # Aimag(ZDR(2,1))
  2 2 0.00000000E+00 # Aimag(ZDR(2,2))
  2 3 0.00000000E+00 # Aimag(ZDR(2,3))
  3 1 0.00000000E+00 # Aimag(ZDR(3,1))
  3 2 0.00000000E+00 # Aimag(ZDR(3,2))
  3 3 0.00000000E+00 # Aimag(ZDR(3,3))
Block UULMIX Q= 1.60000000E+02 # ()
  1 1 9.74272160E-01 # Real(ZUL(1,1),dp)
  1 2 2.25348678E-01 # Real(ZUL(1,2),dp)
  1 3 3.42499367E-03 # Real(ZUL(1,3),dp)
  2 1 -2.25296231E-01 # Real(ZUL(2,1),dp)
  2 2 9.73419462E-01 # Real(ZUL(2,2),dp)
  2 3 4.11844653E-02 # Real(ZUL(2,3),dp)
  3 1 5.94690932E-03 # Real(ZUL(3,1),dp)
  3 2 -4.08965161E-02 # Real(ZUL(3,2),dp)
  3 3 9.99145690E-01 # Real(ZUL(3,3),dp)
Block IMUULMIX Q= 1.60000000E+02 # ()
  1 1 -0.00000000E+00 # Aimag(ZUL(1,1))
  1 2 -0.00000000E+00 # Aimag(ZUL(1,2))
  1 3 -0.00000000E+00 # Aimag(ZUL(1,3))
  2 1 0.00000000E+00 # Aimag(ZUL(2,1))
  2 2 0.00000000E+00 # Aimag(ZUL(2,2))
  2 3 0.00000000E+00 # Aimag(ZUL(2,3))
  3 1 0.00000000E+00 # Aimag(ZUL(3,1))
  3 2 -0.00000000E+00 # Aimag(ZUL(3,2))
  3 3 0.00000000E+00 # Aimag(ZUL(3,3))
Block UURMIX Q= 1.60000000E+02 # ()
  1 1 1.00000000E+00 # Real(ZUR(1,1),dp)
  1 2 -0.00000000E+00 # Real(ZUR(1,2),dp)
  1 3 -0.00000000E+00 # Real(ZUR(1,3),dp)
  2 1 0.00000000E+00 # Real(ZUR(2,1),dp)
  2 2 1.00000000E+00 # Real(ZUR(2,2),dp)
  2 3 0.00000000E+00 # Real(ZUR(2,3),dp)
  3 1 0.00000000E+00 # Real(ZUR(3,1),dp)
  3 2 0.00000000E+00 # Real(ZUR(3,2),dp)
  3 3 1.00000000E+00 # Real(ZUR(3,3),dp)
Block IMUURMIX Q= 1.60000000E+02 # ()
  1 1 0.00000000E+00 # Aimag(ZUR(1,1))
  1 2 0.00000000E+00 # Aimag(ZUR(1,2))
  1 3 0.00000000E+00 # Aimag(ZUR(1,3))
  2 1 -0.00000000E+00 # Aimag(ZUR(2,1))
  2 2 0.00000000E+00 # Aimag(ZUR(2,2))
  2 3 -0.00000000E+00 # Aimag(ZUR(2,3))
  3 1 0.00000000E+00 # Aimag(ZUR(3,1))
  3 2 0.00000000E+00 # Aimag(ZUR(3,2))
  3 3 0.00000000E+00 # Aimag(ZUR(3,3))
Block UELMIX Q= 1.60000000E+02 # ()
  1 1 1.00000000E+00 # Real(ZEL(1,1),dp)
  1 2 0.00000000E+00 # Real(ZEL(1,2),dp)
  1 3 0.00000000E+00 # Real(ZEL(1,3),dp)
  2 1 0.00000000E+00 # Real(ZEL(2,1),dp)
  2 2 1.00000000E+00 # Real(ZEL(2,2),dp)
  2 3 0.00000000E+00 # Real(ZEL(2,3),dp)
  3 1 0.00000000E+00 # Real(ZEL(3,1),dp)
  3 2 0.00000000E+00 # Real(ZEL(3,2),dp)
  3 3 1.00000000E+00 # Real(ZEL(3,3),dp)
Block IMUELMIX Q= 1.60000000E+02 # ()
  1 1 0.00000000E+00 # Aimag(ZEL(1,1))
  1 2 0.00000000E+00 # Aimag(ZEL(1,2))
  1 3 0.00000000E+00 # Aimag(ZEL(1,3))
  2 1 0.00000000E+00 # Aimag(ZEL(2,1))
  2 2 0.00000000E+00 # Aimag(ZEL(2,2))
  2 3 0.00000000E+00 # Aimag(ZEL(2,3))
  3 1 0.00000000E+00 # Aimag(ZEL(3,1))
  3 2 0.00000000E+00 # Aimag(ZEL(3,2))
  3 3 0.00000000E+00 # Aimag(ZEL(3,3))
Block UERMIX Q= 1.60000000E+02 # ()
  1 1 1.00000000E+00 # Real(ZER(1,1),dp)
  1 2 0.00000000E+00 # Real(ZER(1,2),dp)
  1 3 0.00000000E+00 # Real(ZER(1,3),dp)
  2 1 0.00000000E+00 # Real(ZER(2,1),dp)
  2 2 1.00000000E+00 # Real(ZER(2,2),dp)
  2 3 0.00000000E+00 # Real(ZER(2,3),dp)
  3 1 0.00000000E+00 # Real(ZER(3,1),dp)
  3 2 0.00000000E+00 # Real(ZER(3,2),dp)
  3 3 1.00000000E+00 # Real(ZER(3,3),dp)
Block IMUERMIX Q= 1.60000000E+02 # ()
  1 1 0.00000000E+00 # Aimag(ZER(1,1))
  1 2 0.00000000E+00 # Aimag(ZER(1,2))
  1 3 0.00000000E+00 # Aimag(ZER(1,3))
  2 1 0.00000000E+00 # Aimag(ZER(2,1))
  2 2 0.00000000E+00 # Aimag(ZER(2,2))
  2 3 0.00000000E+00 # Aimag(ZER(2,3))
  3 1 0.00000000E+00 # Aimag(ZER(3,1))
  3 2 0.00000000E+00 # Aimag(ZER(3,2))
  3 3 0.00000000E+00 # Aimag(ZER(3,3))
Block SPheno # SPheno internal parameters
         1 -1.00000000E+00 # ErrorLevel
         2 0.00000000E+00 # SPA_conventions
        11 1.00000000E+00 # Branching ratios
        13 1.00000000E+00 # 3 Body decays
        31 -1.00000000E+00 # GUT scale
        33 1.60000000E+02 # Renormalization scale
        34 1.00000000E-04 # Precision
        35 4.00000000E+01 # Iterations
        38 2.00000000E+00 # RGE level
        40 7.29735257E-03 # Alpha
        41 2.49520000E+00 # Gamma_Z
        42 2.06000000E+00 # Gamma_W
        50 1.00000000E+00 # Rotate negative fermion masses
        51 0.00000000E+00 # Switch to SCKM matrix
        52 0.00000000E+00 # Ignore negative masses
        53 0.00000000E+00 # Ignore negative masses at MZ
        55 0.00000000E+00 # Calculate one loop masses
        56 1.00000000E+00 # Calculate two-loop Higgs masses
        57 1.00000000E+00 # Calculate low energy
        60 1.00000000E+00 # Include kinetic mixing
        65 1.00000000E+00 # Solution of tadpole equation
Block HiggsLHC7 # Higgs production cross section at LHC7 [pb] (rescaled SM values from HXSWG)
  1 25 1.56982389E+01 # Gluon fusion
  2 25 1.20912339E+00 # Vector boson fusion
  3 25 5.41279544E-01 # W-H production
  4 25 3.16425082E-01 # Z-H production
  5 25 8.61242763E-02 # t-t-H production
  1 35 2.20252682E+00 # Gluon fusion
  2 35 4.92686931E-17 # Vector boson fusion
  1 45 8.10038284E-02 # Gluon fusion
  2 45 1.06992806E-17 # Vector boson fusion
  1 36 3.81927651E-01 # Gluon fusion
  1 257 5.27196464E-01 # Gluon fusion
Block HiggsLHC8 # Higgs production cross section at LHC8 [pb] (rescaled SM values from HXSWG)
  1 25 1.99812100E+01 # Gluon fusion
  2 25 1.55902259E+00 # Vector boson fusion
  3 25 6.93886067E-01 # W-H production
  4 25 3.95204818E-01 # Z-H production
  5 25 1.30360858E-01 # t-t-H production
  1 35 3.17460634E+00 # Gluon fusion
  2 35 7.85000112E-17 # Vector boson fusion
  1 45 1.16788342E-01 # Gluon fusion
  2 45 1.74642111E-17 # Vector boson fusion
  1 36 5.17219426E-01 # Gluon fusion
  1 257 7.97543668E-01 # Gluon fusion
Block HiggsLHC13 # Higgs production cross section at LHC13 [pb] (rescaled SM values from SusHI 1.5.0)
  1 25 4.66824447E+01 # Gluon fusion
  1 35 9.11435589E+00 # Gluon fusion
  1 45 3.72306027E-01 # Gluon fusion
  1 36 1.33012485E+00 # Gluon fusion
  1 257 2.96039096E+00 # Gluon fusion
Block HiggsLHC14 # Higgs production cross section at LHC14 [pb] (rescaled SM values from SusHI 1.5.0)
  1 25 5.25687295E+01 # Gluon fusion
  1 35 1.06945929E+01 # Gluon fusion
  1 45 4.38960170E-01 # Gluon fusion
  1 36 1.53038147E+00 # Gluon fusion
  1 257 3.53482945E+00 # Gluon fusion
Block HiggsFCC100 # Higgs production cross section at FCC-pp [pb] (rescaled SM values from SusHI 1.5.0)
  1 25 7.74152077E+02 # Gluon fusion
  1 35 3.38327690E+02 # Gluon fusion
  1 45 1.47033189E+01 # Gluon fusion
  1 36 3.46597805E+01 # Gluon fusion
  1 257 1.59743523E+02 # Gluon fusion
Block HPPloops # Loop contributions to H-Photon-Photon coupling
        25 24 8.63392581E+00 # h( 1)-VWm( 1)-loop
        25 37 -3.92334613E-01 # h( 1)-Hm( 2)-loop
        25 256 -3.51228864E-01 # h( 1)-Hm( 3)-loop
        25 1 2.91064547E-07 # h( 1)-Fd( 1)-loop
        25 3 5.75221665E-05 # h( 1)-Fd( 2)-loop
        25 5 2.27857874E-02 # h( 1)-Fd( 3)-loop
        25 2 2.77935527E-07 # h( 1)-Fu( 1)-loop
        25 4 1.29657197E-02 # h( 1)-Fu( 2)-loop
        25 6 -1.76301936E+00 # h( 1)-Fu( 3)-loop
        25 11 1.94419260E-08 # h( 1)-Fe( 1)-loop
        25 13 2.63122590E-04 # h( 1)-Fe( 2)-loop
        25 15 2.35834863E-02 # h( 1)-Fe( 3)-loop
        35 24 -8.55313260E-16 # h( 2)-VWm( 1)-loop
        35 37 -6.92261176E-01 # h( 2)-Hm( 2)-loop
        35 256 3.12629801E-02 # h( 2)-Hm( 3)-loop
        35 1 -2.82212780E-08 # h( 2)-Fd( 1)-loop
        35 3 -6.32151005E-06 # h( 2)-Fd( 2)-loop
        35 5 -3.97045335E-03 # h( 2)-Fd( 3)-loop
        35 2 2.74296749E-08 # h( 2)-Fu( 1)-loop
        35 4 1.85492949E-03 # h( 2)-Fu( 2)-loop
        35 6 -3.14392488E+00 # h( 2)-Fu( 3)-loop
        35 11 -2.47375313E-09 # h( 2)-Fe( 1)-loop
        35 13 -3.85221349E-05 # h( 2)-Fe( 2)-loop
        35 15 -4.41993666E-03 # h( 2)-Fe( 3)-loop
        45 24 6.67189647E-17 # h( 3)-VWm( 1)-loop
        45 37 -7.95579763E-03 # h( 3)-Hm( 2)-loop
        45 256 1.66916615E-02 # h( 3)-Hm( 3)-loop
        45 1 -4.04111664E-08 # h( 3)-Fd( 1)-loop
        45 3 -9.11788999E-06 # h( 3)-Fd( 2)-loop
        45 5 -5.85374538E-03 # h( 3)-Fd( 3)-loop
        45 2 1.08701139E-09 # h( 3)-Fu( 1)-loop
        45 4 7.50075849E-05 # h( 3)-Fu( 2)-loop
        45 6 -1.33824147E-01 # h( 3)-Fu( 3)-loop
        45 11 2.58667579E-09 # h( 3)-Fe( 1)-loop
        45 13 4.05959140E-05 # h( 3)-Fe( 2)-loop
        45 15 4.71397512E-03 # h( 3)-Fe( 3)-loop
Block EFFHIGGSCOUPLINGS # values of loop-induced couplings
        25 22 22 0.29184754E-04 # H-Photon-Photon
        25 21 21 0.66151673E-04 # H-Gluon-Gluon
        25 22 23 0.00000000E+00 # H-Photon-Z (not yet calculated by SPheno)
        35 22 22 0.18712749E-04 # H-Photon-Photon
        35 21 21 0.97173293E-04 # H-Gluon-Gluon
        35 22 23 0.00000000E+00 # H-Photon-Z (not yet calculated by SPheno)
        45 22 22 0.68562167E-06 # H-Photon-Photon
        45 21 21 0.47958187E-05 # H-Gluon-Gluon
        45 22 23 0.00000000E+00 # H-Photon-Z (not yet calculated by SPheno)
        36 22 22 0.89339079E-06 # A-Photon-Photon
        36 21 21 0.77337287E-05 # A-Gluon-Gluon
        36 22 23 0.00000000E+00 # A-Photon-Z (not yet calculated by SPheno)
       257 22 22 0.16436772E-04 # A-Photon-Photon
       257 21 21 0.10239533E-03 # A-Gluon-Gluon
       257 22 23 0.00000000E+00 # A-Photon-Z (not yet calculated by SPheno)
Block SPhenoLowEnergy # low energy observables
       1 -3.53301165E-02 # T-parameter (1-loop BSM)
       2 2.60525255E-02 # S-parameter (1-loop BSM)
       3 6.66835230E-04 # U-parameter (1-loop BSM)
      20 4.42178662E-25 # (g-2)_e
      21 1.64147667E-14 # (g-2)_mu
      22 7.45591141E-10 # (g-2)_tau
      23 0.00000000E+00 # EDM(e)
      24 0.00000000E+00 # EDM(mu)
      25 0.00000000E+00 # EDM(tau)
      39 -6.02010427E-04 # delta(rho)
Block FlavorKitQFV # quark flavor violating observables
     200 4.03861355E-04 # BR(B->X_s gamma)
     201 1.28209954E+00 # BR(B->X_s gamma)/BR(B->X_s gamma)_SM
     300 5.89318485E-04 # BR(D->mu nu)
     301 9.98996285E-01 # BR(D->mu nu)/BR(D->mu nu)_SM
     400 5.57280075E-03 # BR(Ds->mu nu)
     401 9.98898588E-01 # BR(Ds->mu nu)/BR(Ds->mu nu)_SM
     402 5.44634009E-02 # BR(Ds->tau nu)
     403 9.98898588E-01 # BR(Ds->tau nu)/BR(Ds->tau nu)_SM
     500 5.06715465E-07 # BR(B->mu nu)
     501 9.93275869E-01 # BR(B->mu nu)/BR(B->mu nu)_SM
     502 1.12750805E-04 # BR(B->tau nu)
     503 9.93275869E-01 # BR(B->tau nu)/BR(B->tau nu)_SM
     600 6.32338420E-01 # BR(K->mu nu)
     601 9.99940821E-01 # BR(K->mu nu)/BR(K->mu nu)_SM
     602 2.47647734E-05 # R_K = BR(K->e nu)/(K->mu nu)
     603 2.47618424E-05 # R_K^SM = BR(K->e nu)_SM/(K->mu nu)_SM
    1900 2.43286332E+01 # Delta(M_Bs)
    1901 1.29426358E+00 # Delta(M_Bs)/Delta(M_Bs)_SM
    1902 4.63887614E-01 # Delta(M_Bd)
    1903 1.10411095E+00 # Delta(M_Bd)/Delta(M_Bd)_SM
    4000 1.32312788E-15 # BR(B^0_d->e e)
    4001 1.26022751E+00 # BR(B^0_d->e e)/BR(B^0_d->e e)_SM
    4002 8.96080382E-14 # BR(B^0_s->e e)
    4003 1.26025639E+00 # BR(B^0_s->e e)/BR(B^0_s->e e)_SM
    4004 5.65225279E-11 # BR(B^0_d->mu mu)
    4005 1.26022736E+00 # BR(B^0_d->mu mu)/BR(B^0_d->mu mu)_SM
    4006 3.82805310E-09 # BR(B^0_s->mu mu)
    4007 1.26025624E+00 # BR(B^0_s->mu mu)/BR(B^0_s->mu mu)_SM
    4008 1.18304302E-08 # BR(B^0_d->tau tau)
    4009 1.26018553E+00 # BR(B^0_d->tau tau)/BR(B^0_d->tau tau)_SM
    4010 8.11829966E-07 # BR(B^0_s->tau tau)
    4011 1.26021440E+00 # BR(B^0_s->tau tau)/BR(B^0_s->tau tau)_SM
    5000 1.93490809E-06 # BR(B-> s e e)
    5001 1.16892393E+00 # BR(B-> s e e)/BR(B-> s e e)_SM
    5002 1.88127997E-06 # BR(B-> s mu mu)
    5003 1.17229043E+00 # BR(B-> s mu mu)/BR(B-> s mu mu)_SM
    6000 1.25213625E-07 # BR(B -> K mu mu)
    6001 1.12805068E+00 # BR(B -> K mu mu)/BR(B -> K mu mu)_SM
    6002 1.25213625E-07 # BR(B -> K mu mu)
    6003 1.12805068E+00 # BR(B -> K mu mu)/BR(B -> K mu mu)_SM
    7000 4.39778779E-05 # BR(B->s nu nu)
    7001 1.16479365E+00 # BR(B->s nu nu)/BR(B->s nu nu)_SM
    7002 9.32878435E-07 # BR(B->D nu nu)
    7003 1.16487858E+00 # BR(B->D nu nu)/BR(B->D nu nu)_SM
    8000 1.38414541E-10 # BR(K^+ -> pi^+ nu nu)
    8001 1.06111099E+00 # BR(K^+ -> pi^+ nu nu)/BR(K^+ -> pi^+ nu nu)_SM
    8002 1.62288751E-43 # BR(K_L -> pi^0 nu nu)
    8003 1.00000000E+00 # BR(K_L -> pi^0 nu nu)/BR(K_L -> pi^0 nu nu)_SM
    8004 0.00000000E+00 # BR(K^0_L -> e mu)
    8005 0.00000000E+00 # BR(K^0_L -> e mu)/BR(K^0_L -> e mu)_SM
    9100 1.95208577E-15 # Delta(M_K)
    9102 1.00291023E+00 # Delta(M_K)/Delta(M_K)_SM
    9103 1.83620120E-03 # epsilon_K
    9104 1.00000000E+00 # epsilon_K/epsilon_K^SM
Block FlavorKitLFV # lepton flavor violating observables
     701 0.00000000E+00 # BR(mu->e gamma)
     702 0.00000000E+00 # BR(tau->e gamma)
     703 0.00000000E+00 # BR(tau->mu gamma)
     800 0.00000000E+00 # CR(mu-e, Al)
     801 0.00000000E+00 # CR(mu-e, Ti)
     802 0.00000000E+00 # CR(mu-e, Sr)
     803 0.00000000E+00 # CR(mu-e, Sb)
     804 0.00000000E+00 # CR(mu-e, Au)
     805 0.00000000E+00 # CR(mu-e, Pb)
     901 0.00000000E+00 # BR(mu->3e)
     902 0.00000000E+00 # BR(tau->3e)
     903 0.00000000E+00 # BR(tau->3mu)
     904 0.00000000E+00 # BR(tau- -> e- mu+ mu-)
     905 0.00000000E+00 # BR(tau- -> mu- e+ e-)
     906 0.00000000E+00 # BR(tau- -> e+ mu- mu-)
     907 0.00000000E+00 # BR(tau- -> mu+ e- e-)
    1001 0.00000000E+00 # BR(Z->e mu)
    1002 0.00000000E+00 # BR(Z->e tau)
    1003 0.00000000E+00 # BR(Z->mu tau)
    1101 0.00000000E+00 # BR(h->e mu)
    1102 0.00000000E+00 # BR(h->e tau)
    1103 0.00000000E+00 # BR(h->mu tau)
    2001 0.00000000E+00 # BR(tau->e pi)
    2002 0.00000000E+00 # BR(tau->e eta)
    2003 0.00000000E+00 # BR(tau->e eta')
    2004 0.00000000E+00 # BR(tau->mu pi)
    2005 0.00000000E+00 # BR(tau->mu eta)
    2006 0.00000000E+00 # BR(tau->mu eta')
Block FWCOEF Q= 1.60000000E+02 # Wilson coefficients at scale Q
     0305 4422 00 0 -0.16426520E-08 # coeffC7sm
     0305 4422 00 2 -0.23932721E-08 # coeffC7
     0305 4322 00 2 -0.54343583E-10 # coeffC7p
     0305 4422 00 1 -0.75062014E-09 # coeffC7NP
     0305 4322 00 1 -0.54343583E-10 # coeffC7pNP
     0305 6421 00 0 -0.85359897E-09 # coeffC8sm
     0305 6421 00 2 -0.16425894E-08 # coeffC8
     0305 6321 00 2 -0.37283971E-10 # coeffC8p
     0305 6421 00 1 -0.78899041E-09 # coeffC8NP
     0305 6321 00 1 -0.37283971E-10 # coeffC8pNP
 03051111 4133 00 0 0.10133674E-08 # coeffC9eeSM
 03051111 4133 00 2 0.10622705E-08 # coeffC9ee
 03051111 4233 00 2 -0.53955498E-14 # coeffC9Pee
 03051111 4133 00 1 0.48903106E-10 # coeffC9eeNP
 03051111 4233 00 1 -0.53955498E-14 # coeffC9PeeNP
 03051111 4137 00 0 -0.37904146E-08 # coeffC10eeSM
 03051111 4137 00 2 -0.42534735E-08 # coeffC10ee
 03051111 4237 00 2 -0.51089950E-13 # coeffC10Pee
 03051111 4137 00 1 -0.46305888E-09 # coeffC10eeNP
 03051111 4237 00 1 0.51089950E-13 # coeffC10PeeNP
 03051313 4133 00 0 0.10133669E-08 # coeffC9mumuSM
 03051313 4133 00 2 0.10622702E-08 # coeffC9mumu
 03051313 4233 00 2 -0.53955422E-14 # coeffC9Pmumu
 03051313 4133 00 1 0.48903300E-10 # coeffC9mumuNP
 03051313 4233 00 1 -0.53955422E-14 # coeffC9PmumuNP
 03051313 4137 00 0 -0.37904151E-08 # coeffC10mumuSM
 03051313 4137 00 2 -0.42534738E-08 # coeffC10mumu
 03051313 4237 00 2 -0.51089957E-13 # coeffC10Pmumu
 03051313 4137 00 1 -0.46305869E-09 # coeffC10mumuNP
 03051313 4237 00 1 -0.51089957E-13 # coeffC10PmumuNP
 03051212 4141 00 0 -0.11693567E-07 # coeffCLnu1nu1SM
 03051212 4141 00 2 -0.12619685E-07 # coeffCLnu1nu1
 03051212 4241 00 2 0.10217990E-12 # coeffCLPnu1nu1
 03051212 4141 00 1 -0.92611777E-09 # coeffCLnu1nu1NP
 03051212 4241 00 1 0.10217990E-12 # coeffCLPnu1nu1NP
 03051414 4141 00 0 -0.11693471E-07 # coeffCLnu2nu2SM
 03051414 4141 00 2 -0.12619590E-07 # coeffCLnu2nu2
 03051414 4241 00 2 0.10217988E-12 # coeffCLPnu2nu2
 03051414 4141 00 1 -0.92611858E-09 # coeffCLnu2nu2NP
 03051414 4241 00 1 0.10217988E-12 # coeffCLPnu2nu2NP
 03051616 4141 00 0 -0.11671336E-07 # coeffCLnu3nu3SM
 03051616 4141 00 2 -0.12597683E-07 # coeffCLnu3nu3
 03051616 4241 00 2 0.10217563E-12 # coeffCLPnu3nu3
 03051616 4141 00 1 -0.92634696E-09 # coeffCLnu3nu3NP
 03051616 4241 00 1 0.10217563E-12 # coeffCLPnu3nu3NP
 03051212 4142 00 0 0.00000000E+00 # coeffCRnu1nu1SM
 03051212 4142 00 2 0.00000000E+00 # coeffCRnu1nu1
 03051212 4242 00 2 0.00000000E+00 # coeffCRPnu1nu1
 03051212 4142 00 1 0.00000000E+00 # coeffCRnu1nu1NP
 03051212 4242 00 1 0.00000000E+00 # coeffCRPnu1nu1NP
 03051414 4142 00 0 0.00000000E+00 # coeffCRnu2nu2SM
 03051414 4142 00 2 0.00000000E+00 # coeffCRnu2nu2
 03051414 4242 00 2 0.00000000E+00 # coeffCRPnu2nu2
 03051414 4142 00 1 0.00000000E+00 # coeffCRnu2nu2NP
 03051414 4242 00 1 0.00000000E+00 # coeffCRPnu2nu2NP
 03051616 4142 00 0 0.00000000E+00 # coeffCRnu3nu3SM
 03051616 4142 00 2 0.00000000E+00 # coeffCRnu3nu3
 03051616 4242 00 2 0.00000000E+00 # coeffCRPnu3nu3
 03051616 4142 00 1 0.00000000E+00 # coeffCRnu3nu3NP
 03051616 4242 00 1 0.00000000E+00 # coeffCRPnu3nu3NP
 01030103 3131 00 2 -0.52526240E-23 # coeffKK_SLL
 01030103 3232 00 2 -0.18961973E-20 # coeffKK_SRR
 01030103 3132 00 2 -0.10046273E-18 # coeffKK_SLR
 01030103 4141 00 2 0.27709052E-12 # coeffKK_VLL
 01030103 4242 00 2 0.57095046E-27 # coeffKK_VRR
 01030103 4142 00 2 -0.99799856E-22 # coeffKK_VLR
 01030103 4343 00 2 0.00000000E+00 # coeffKK_TLL
 01030103 4444 00 2 0.00000000E+00 # coeffKK_TRR
 01050105 3131 00 2 -0.14238374E-20 # coeffBB_SLL
 01050105 3232 00 2 -0.99511430E-15 # coeffBB_SRR
 01050105 3132 00 2 -0.22241159E-16 # coeffBB_SLR
 01050105 4141 00 2 0.25882247E-11 # coeffBB_VLL
 01050105 4242 00 2 0.26787235E-23 # coeffBB_VRR
 01050105 4142 00 2 -0.11903281E-17 # coeffBB_VLR
 01050105 4343 00 2 0.00000000E+00 # coeffBB_TLL
 01050105 4444 00 2 0.00000000E+00 # coeffBB_TRR
 03050305 3131 00 2 -0.24317776E-16 # coeffBsBs_SLL
 03050305 3232 00 2 -0.47079215E-13 # coeffBsBs_SRR
 03050305 3132 00 2 -0.20007937E-13 # coeffBsBs_SLR
 03050305 4141 00 2 0.12248856E-09 # coeffBsBs_VLL
 03050305 4242 00 2 0.45830182E-19 # coeffBsBs_VRR
 03050305 4142 00 2 -0.10699822E-14 # coeffBsBs_VLR
 03050305 4343 00 2 0.00000000E+00 # coeffBsBs_TLL
 03050305 4444 00 2 0.00000000E+00 # coeffBsBs_TRR
 01030103 3131 00 1 0.97209159E-24 # coeffKK_SLLNP
 01030103 3232 00 1 0.35092506E-21 # coeffKK_SRRNP
 01030103 3132 00 1 -0.22521205E-19 # coeffKK_SLRNP
 01030103 4141 00 1 0.11975525E-14 # coeffKK_VLLNP
 01030103 4242 00 1 0.22553114E-27 # coeffKK_VRRNP
 01030103 4142 00 1 0.18469740E-22 # coeffKK_VLRNP
 01030103 4343 00 1 0.00000000E+00 # coeffKK_TLLNP
 01030103 4444 00 1 0.00000000E+00 # coeffKK_TRRNP
 01050105 3131 00 1 0.39929168E-21 # coeffBB_SLLNP
 01050105 3232 00 1 0.27906336E-15 # coeffBB_SRRNP
 01050105 3132 00 1 -0.10098952E-16 # coeffBB_SLRNP
 01050105 4141 00 1 0.59102972E-12 # coeffBB_VLLNP
 01050105 4242 00 1 0.14926315E-23 # coeffBB_VRRNP
 01050105 4142 00 1 0.33380784E-18 # coeffBB_VLRNP
 01050105 4343 00 1 0.00000000E+00 # coeffBB_TLLNP
 01050105 4444 00 1 0.00000000E+00 # coeffBB_TRRNP
 03050305 3131 00 1 0.68182500E-17 # coeffBsBs_SLLNP
 03050305 3232 00 1 0.13200132E-13 # coeffBsBs_SRRNP
 03050305 3132 00 1 -0.90804866E-14 # coeffBsBs_SLRNP
 03050305 4141 00 1 0.27954211E-10 # coeffBsBs_VLLNP
 03050305 4242 00 1 0.25524100E-19 # coeffBsBs_VRRNP
 03050305 4142 00 1 0.30000300E-15 # coeffBsBs_VLRNP
 03050305 4343 00 1 0.00000000E+00 # coeffBsBs_TLLNP
 03050305 4444 00 1 0.00000000E+00 # coeffBsBs_TRRNP
 01030103 3131 00 0 -0.62247156E-23 # coeffKK_SLLSM
 01030103 3232 00 0 -0.22471223E-20 # coeffKK_SRRSM
 01030103 3132 00 0 -0.77941523E-19 # coeffKK_SLRSM
 01030103 4141 00 0 0.27589297E-12 # coeffKK_VLLSM
 01030103 4242 00 0 0.34541933E-27 # coeffKK_VRRSM
 01030103 4142 00 0 -0.11826960E-21 # coeffKK_VLRSM
 01030103 4343 00 0 0.00000000E+00 # coeffKK_TLLSM
 01030103 4444 00 0 0.00000000E+00 # coeffKK_TRRSM
 01050105 3131 00 0 -0.18231291E-20 # coeffBB_SLLSM
 01050105 3232 00 0 -0.12741777E-14 # coeffBB_SRRSM
 01050105 3132 00 0 -0.12142207E-16 # coeffBB_SLRSM
 01050105 4141 00 0 0.19971950E-11 # coeffBB_VLLSM
 01050105 4242 00 0 0.11860921E-23 # coeffBB_VRRSM
 01050105 4142 00 0 -0.15241360E-17 # coeffBB_VLRSM
 01050105 4343 00 0 0.00000000E+00 # coeffBB_TLLSM
 01050105 4444 00 0 0.00000000E+00 # coeffBB_TRRSM
 03050305 3131 00 0 -0.31136026E-16 # coeffBsBs_SLLSM
 03050305 3232 00 0 -0.60279347E-13 # coeffBsBs_SRRSM
 03050305 3132 00 0 -0.10927451E-13 # coeffBsBs_SLRSM
 03050305 4141 00 0 0.94534346E-10 # coeffBsBs_VLLSM
 03050305 4242 00 0 0.20306081E-19 # coeffBsBs_VRRSM
 03050305 4142 00 0 -0.13699852E-14 # coeffBsBs_VLRSM
 03050305 4343 00 0 0.00000000E+00 # coeffBsBs_TLLSM
 03050305 4444 00 0 0.00000000E+00 # coeffBsBs_TRRSM
Block IMFWCOEF Q= 1.60000000E+02 # Im(Wilson coefficients) at scale Q
     0305 4422 00 0 -0.00000000E+00 # coeffC7sm
     0305 4422 00 2 -0.00000000E+00 # coeffC7
     0305 4322 00 2 -0.00000000E+00 # coeffC7p
     0305 4422 00 1 0.00000000E+00 # coeffC7NP
     0305 4322 00 1 -0.00000000E+00 # coeffC7pNP
     0305 6421 00 0 -0.00000000E+00 # coeffC8sm
     0305 6421 00 2 -0.00000000E+00 # coeffC8
     0305 6321 00 2 -0.00000000E+00 # coeffC8p
     0305 6421 00 1 0.00000000E+00 # coeffC8NP
     0305 6321 00 1 -0.00000000E+00 # coeffC8pNP
 03051111 4133 00 0 0.00000000E+00 # coeffC9eeSM
 03051111 4133 00 2 0.00000000E+00 # coeffC9ee
 03051111 4233 00 2 0.00000000E+00 # coeffC9Pee
 03051111 4133 00 1 0.00000000E+00 # coeffC9eeNP
 03051111 4233 00 1 0.00000000E+00 # coeffC9PeeNP
 03051111 4137 00 0 0.00000000E+00 # coeffC10eeSM
 03051111 4137 00 2 0.00000000E+00 # coeffC10ee
 03051111 4237 00 2 0.00000000E+00 # coeffC10Pee
 03051111 4137 00 1 0.00000000E+00 # coeffC10eeNP
 03051111 4237 00 1 0.00000000E+00 # coeffC10PeeNP
 03051313 4133 00 0 0.00000000E+00 # coeffC9mumuSM
 03051313 4133 00 2 0.00000000E+00 # coeffC9mumu
 03051313 4233 00 2 0.00000000E+00 # coeffC9Pmumu
 03051313 4133 00 1 0.00000000E+00 # coeffC9mumuNP
 03051313 4233 00 1 0.00000000E+00 # coeffC9PmumuNP
 03051313 4137 00 0 0.00000000E+00 # coeffC10mumuSM
 03051313 4137 00 2 0.00000000E+00 # coeffC10mumu
 03051313 4237 00 2 0.00000000E+00 # coeffC10Pmumu
 03051313 4137 00 1 0.00000000E+00 # coeffC10mumuNP
 03051313 4237 00 1 0.00000000E+00 # coeffC10PmumuNP
 03051212 4141 00 0 0.00000000E+00 # coeffCLnu1nu1SM
 03051212 4141 00 2 0.00000000E+00 # coeffCLnu1nu1
 03051212 4241 00 2 0.00000000E+00 # coeffCLPnu1nu1
 03051212 4141 00 1 0.00000000E+00 # coeffCLnu1nu1NP
 03051212 4241 00 1 0.00000000E+00 # coeffCLPnu1nu1NP
 03051414 4141 00 0 0.00000000E+00 # coeffCLnu2nu2SM
 03051414 4141 00 2 0.00000000E+00 # coeffCLnu2nu2
 03051414 4241 00 2 0.00000000E+00 # coeffCLPnu2nu2
 03051414 4141 00 1 0.00000000E+00 # coeffCLnu2nu2NP
 03051414 4241 00 1 0.00000000E+00 # coeffCLPnu2nu2NP
 03051616 4141 00 0 0.00000000E+00 # coeffCLnu3nu3SM
 03051616 4141 00 2 0.00000000E+00 # coeffCLnu3nu3
 03051616 4241 00 2 0.00000000E+00 # coeffCLPnu3nu3
 03051616 4141 00 1 0.00000000E+00 # coeffCLnu3nu3NP
 03051616 4241 00 1 0.00000000E+00 # coeffCLPnu3nu3NP
 03051212 4142 00 0 0.00000000E+00 # coeffCRnu1nu1SM
 03051212 4142 00 2 0.00000000E+00 # coeffCRnu1nu1
 03051212 4242 00 2 0.00000000E+00 # coeffCRPnu1nu1
 03051212 4142 00 1 0.00000000E+00 # coeffCRnu1nu1NP
 03051212 4242 00 1 0.00000000E+00 # coeffCRPnu1nu1NP
 03051414 4142 00 0 0.00000000E+00 # coeffCRnu2nu2SM
 03051414 4142 00 2 0.00000000E+00 # coeffCRnu2nu2
 03051414 4242 00 2 0.00000000E+00 # coeffCRPnu2nu2
 03051414 4142 00 1 0.00000000E+00 # coeffCRnu2nu2NP
 03051414 4242 00 1 0.00000000E+00 # coeffCRPnu2nu2NP
 03051616 4142 00 0 0.00000000E+00 # coeffCRnu3nu3SM
 03051616 4142 00 2 0.00000000E+00 # coeffCRnu3nu3
 03051616 4242 00 2 0.00000000E+00 # coeffCRPnu3nu3
 03051616 4142 00 1 0.00000000E+00 # coeffCRnu3nu3NP
 03051616 4242 00 1 0.00000000E+00 # coeffCRPnu3nu3NP
 01030103 3131 00 2 0.00000000E+00 # coeffKK_SLL
 01030103 3232 00 2 0.00000000E+00 # coeffKK_SRR
 01030103 3132 00 2 0.00000000E+00 # coeffKK_SLR
 01030103 4141 00 2 0.00000000E+00 # coeffKK_VLL
 01030103 4242 00 2 0.00000000E+00 # coeffKK_VRR
 01030103 4142 00 2 0.00000000E+00 # coeffKK_VLR
 01030103 4343 00 2 0.00000000E+00 # coeffKK_TLL
 01030103 4444 00 2 0.00000000E+00 # coeffKK_TRR
 01050105 3131 00 2 0.00000000E+00 # coeffBB_SLL
 01050105 3232 00 2 0.00000000E+00 # coeffBB_SRR
 01050105 3132 00 2 0.00000000E+00 # coeffBB_SLR
 01050105 4141 00 2 0.00000000E+00 # coeffBB_VLL
 01050105 4242 00 2 0.00000000E+00 # coeffBB_VRR
 01050105 4142 00 2 0.00000000E+00 # coeffBB_VLR
 01050105 4343 00 2 0.00000000E+00 # coeffBB_TLL
 01050105 4444 00 2 0.00000000E+00 # coeffBB_TRR
 03050305 3131 00 2 0.00000000E+00 # coeffBsBs_SLL
 03050305 3232 00 2 0.00000000E+00 # coeffBsBs_SRR
 03050305 3132 00 2 0.00000000E+00 # coeffBsBs_SLR
 03050305 4141 00 2 0.00000000E+00 # coeffBsBs_VLL
 03050305 4242 00 2 0.00000000E+00 # coeffBsBs_VRR
 03050305 4142 00 2 0.00000000E+00 # coeffBsBs_VLR
 03050305 4343 00 2 0.00000000E+00 # coeffBsBs_TLL
 03050305 4444 00 2 0.00000000E+00 # coeffBsBs_TRR
 01030103 3131 00 1 0.00000000E+00 # coeffKK_SLLNP
 01030103 3232 00 1 0.00000000E+00 # coeffKK_SRRNP
 01030103 3132 00 1 0.00000000E+00 # coeffKK_SLRNP
 01030103 4141 00 1 0.00000000E+00 # coeffKK_VLLNP
 01030103 4242 00 1 0.00000000E+00 # coeffKK_VRRNP
 01030103 4142 00 1 0.00000000E+00 # coeffKK_VLRNP
 01030103 4343 00 1 0.00000000E+00 # coeffKK_TLLNP
 01030103 4444 00 1 0.00000000E+00 # coeffKK_TRRNP
 01050105 3131 00 1 0.00000000E+00 # coeffBB_SLLNP
 01050105 3232 00 1 0.00000000E+00 # coeffBB_SRRNP
 01050105 3132 00 1 0.00000000E+00 # coeffBB_SLRNP
 01050105 4141 00 1 0.00000000E+00 # coeffBB_VLLNP
 01050105 4242 00 1 0.00000000E+00 # coeffBB_VRRNP
 01050105 4142 00 1 0.00000000E+00 # coeffBB_VLRNP
 01050105 4343 00 1 0.00000000E+00 # coeffBB_TLLNP
 01050105 4444 00 1 0.00000000E+00 # coeffBB_TRRNP
 03050305 3131 00 1 0.00000000E+00 # coeffBsBs_SLLNP
 03050305 3232 00 1 0.00000000E+00 # coeffBsBs_SRRNP
 03050305 3132 00 1 0.00000000E+00 # coeffBsBs_SLRNP
 03050305 4141 00 1 0.00000000E+00 # coeffBsBs_VLLNP
 03050305 4242 00 1 0.00000000E+00 # coeffBsBs_VRRNP
 03050305 4142 00 1 0.00000000E+00 # coeffBsBs_VLRNP
 03050305 4343 00 1 0.00000000E+00 # coeffBsBs_TLLNP
 03050305 4444 00 1 0.00000000E+00 # coeffBsBs_TRRNP
 01030103 3131 00 0 0.00000000E+00 # coeffKK_SLLSM
 01030103 3232 00 0 0.00000000E+00 # coeffKK_SRRSM
 01030103 3132 00 0 0.00000000E+00 # coeffKK_SLRSM
 01030103 4141 00 0 0.00000000E+00 # coeffKK_VLLSM
 01030103 4242 00 0 0.00000000E+00 # coeffKK_VRRSM
 01030103 4142 00 0 0.00000000E+00 # coeffKK_VLRSM
 01030103 4343 00 0 0.00000000E+00 # coeffKK_TLLSM
 01030103 4444 00 0 0.00000000E+00 # coeffKK_TRRSM
 01050105 3131 00 0 0.00000000E+00 # coeffBB_SLLSM
 01050105 3232 00 0 0.00000000E+00 # coeffBB_SRRSM
 01050105 3132 00 0 0.00000000E+00 # coeffBB_SLRSM
 01050105 4141 00 0 0.00000000E+00 # coeffBB_VLLSM
 01050105 4242 00 0 0.00000000E+00 # coeffBB_VRRSM
 01050105 4142 00 0 0.00000000E+00 # coeffBB_VLRSM
 01050105 4343 00 0 0.00000000E+00 # coeffBB_TLLSM
 01050105 4444 00 0 0.00000000E+00 # coeffBB_TRRSM
 03050305 3131 00 0 0.00000000E+00 # coeffBsBs_SLLSM
 03050305 3232 00 0 0.00000000E+00 # coeffBsBs_SRRSM
 03050305 3132 00 0 0.00000000E+00 # coeffBsBs_SLRSM
 03050305 4141 00 0 0.00000000E+00 # coeffBsBs_VLLSM
 03050305 4242 00 0 0.00000000E+00 # coeffBsBs_VRRSM
 03050305 4142 00 0 0.00000000E+00 # coeffBsBs_VLRSM
 03050305 4343 00 0 0.00000000E+00 # coeffBsBs_TLLSM
 03050305 4444 00 0 0.00000000E+00 # coeffBsBs_TRRSM
Block TREELEVELUNITARITY #
       0 1.00000000E+00 # Tree-level unitarity limits fulfilled or not
       1 4.36228263E-01 # Maximal scattering eigenvalue
Block TREELEVELUNITARITYwTRILINEARS #
       0 1.00000000E+00 # Tree-level unitarity limits fulfilled or not
       1 4.07830496E-01 # Maximal scattering eigenvalue
       2 5.00000000E+03 # best scattering energy
      11 5.00000000E+02 # min scattering energy
      12 5.00000000E+03 # max scattering energy
      13 5.00000000E+00 # steps
DECAY 2 0.00000000E+00 # Fu_1
# BR NDA ID1 ID2
DECAY 4 3.79878857E-13 # Fu_2
# BR NDA ID1 ID2
# BR NDA ID1 ID2 ID3
     3.05502575E-02 3 2 -1 1 # BR(Fu_2 -> Fu_1 Fd_1^* Fd_1 )
     5.45954987E-01 3 2 -1 3 # BR(Fu_2 -> Fu_1 Fd_1^* Fd_2 )
     1.56486313E-03 3 2 -3 1 # BR(Fu_2 -> Fu_1 Fd_2^* Fd_1 )
     2.79270156E-02 3 2 -3 3 # BR(Fu_2 -> Fu_1 Fd_2^* Fd_2 )
     1.07295183E-02 3 1 -11 12 # BR(Fu_2 -> Fd_1 Fe_1^* Fv_1 )
     1.01645236E-02 3 1 -13 14 # BR(Fu_2 -> Fd_1 Fe_2^* Fv_2 )
     1.91744771E-01 3 3 -11 12 # BR(Fu_2 -> Fd_2 Fe_1^* Fv_1 )
     1.81364064E-01 3 3 -13 14 # BR(Fu_2 -> Fd_2 Fe_2^* Fv_2 )
DECAY 6 1.53621095E+00 # Fu_3
# BR NDA ID1 ID2
     1.67601048E-03 2 3 24 # BR(Fu_3 -> Fd_2 VWm^* )
     9.98288545E-01 2 5 24 # BR(Fu_3 -> Fd_3 VWm^* )
DECAY 11 0.00000000E+00 # Fe_1
# BR NDA ID1 ID2
DECAY 13 1.21230600E-18 # Fe_2
# BR NDA ID1 ID2
# BR NDA ID1 ID2 ID3
     2.63969808E-01 3 11 -12 14 # BR(Fe_2 -> Fe_1 Fv_1^* Fv_2 )
     7.36026661E-01 3 14 -2 1 # BR(Fe_2 -> Fv_2 Fu_1^* Fd_1 )
DECAY 15 2.15561802E-12 # Fe_3
# BR NDA ID1 ID2
# BR NDA ID1 ID2 ID3
     1.99664865E-01 3 11 -12 16 # BR(Fe_3 -> Fe_1 Fv_1^* Fv_3 )
     1.94368619E-01 3 13 -14 16 # BR(Fe_3 -> Fe_2 Fv_2^* Fv_3 )
     5.68537963E-01 3 16 -2 1 # BR(Fe_3 -> Fv_3 Fu_1^* Fd_1 )
     2.97768354E-02 3 16 -2 3 # BR(Fe_3 -> Fv_3 Fu_1^* Fd_2 )
     4.57258377E-04 3 16 -4 1 # BR(Fe_3 -> Fv_3 Fu_2^* Fd_1 )
     7.19445872E-03 3 16 -4 3 # BR(Fe_3 -> Fv_3 Fu_2^* Fd_2 )
DECAY 1 9.39808485E-27 # Fd_1
# BR NDA ID1 ID2
# BR NDA ID1 ID2 ID3
     1.00000000E+00 3 2 -12 11 # BR(Fd_1 -> Fu_1 Fv_1^* Fe_1 )
DECAY 3 3.58320026E-20 # Fd_2
# BR NDA ID1 ID2
# BR NDA ID1 ID2 ID3
     7.34967004E-01 3 1 -2 2 # BR(Fd_2 -> Fd_1 Fu_1^* Fu_1 )
     2.65032996E-01 3 2 -12 11 # BR(Fd_2 -> Fu_1 Fv_1^* Fe_1 )
DECAY 5 1.69296025E-13 # Fd_3
# BR NDA ID1 ID2
# BR NDA ID1 ID2 ID3
     6.12948851E-03 3 1 -2 2 # BR(Fd_3 -> Fd_1 Fu_1^* Fu_1 )
     4.52961996E-01 3 1 -2 4 # BR(Fd_3 -> Fd_1 Fu_1^* Fu_2 )
     1.67716645E-04 3 1 -4 2 # BR(Fd_3 -> Fd_1 Fu_2^* Fu_1 )
     8.80820957E-03 3 1 -4 4 # BR(Fd_3 -> Fd_1 Fu_2^* Fu_2 )
     3.27012200E-04 3 3 -2 2 # BR(Fd_3 -> Fd_2 Fu_1^* Fu_1 )
     2.41185539E-02 3 3 -2 4 # BR(Fd_3 -> Fd_2 Fu_1^* Fu_2 )
     3.11229758E-03 3 3 -4 2 # BR(Fd_3 -> Fd_2 Fu_2^* Fu_1 )
     1.62630977E-01 3 3 -4 4 # BR(Fd_3 -> Fd_2 Fu_2^* Fu_2 )
     2.15251763E-03 3 2 -12 11 # BR(Fd_3 -> Fu_1 Fv_1^* Fe_1 )
     2.14489928E-03 3 2 -14 13 # BR(Fd_3 -> Fu_1 Fv_2^* Fe_2 )
     5.83955549E-04 3 2 -16 15 # BR(Fd_3 -> Fu_1 Fv_3^* Fe_3 )
     1.59069257E-01 3 4 -12 11 # BR(Fd_3 -> Fu_2 Fv_1^* Fe_1 )
     1.58103583E-01 3 4 -14 13 # BR(Fd_3 -> Fu_2 Fv_2^* Fe_2 )
     1.96895356E-02 3 4 -16 15 # BR(Fd_3 -> Fu_2 Fv_3^* Fe_3 )
DECAY 25 3.58125124E-03 # hh_1
# BR NDA ID1 ID2
     2.31034848E-03 2 22 22 # BR(hh_1 -> VP VP )
     9.49592003E-02 2 21 21 # BR(hh_1 -> VG VG )
     2.53405594E-02 2 23 23 # BR(hh_1 -> VZ VZ )
     2.34518976E-01 2 24 -24 # BR(hh_1 -> VWm^* VWm_virt )
     2.04186700E-04 2 -3 3 # BR(hh_1 -> Fd_2^* Fd_2 )
     5.44594035E-01 2 -5 5 # BR(hh_1 -> Fd_3^* Fd_3 )
     2.49561413E-04 2 -13 13 # BR(hh_1 -> Fe_2^* Fe_2 )
     7.20361528E-02 2 -15 15 # BR(hh_1 -> Fe_3^* Fe_3 )
     2.57862980E-02 2 -4 4 # BR(hh_1 -> Fu_2^* Fu_2 )
DECAY 35 2.93058900E+00 # hh_2
# BR NDA ID1 ID2
     8.10642230E-03 2 21 21 # BR(hh_2 -> VG VG )
     6.89754706E-03 2 36 23 # BR(hh_2 -> Ah_2 VZ )
     1.62239432E-03 2 -5 5 # BR(hh_2 -> Fd_3^* Fd_3 )
     3.25778822E-04 2 -15 15 # BR(hh_2 -> Fe_3^* Fe_3 )
     9.53760479E-01 2 -6 6 # BR(hh_2 -> Fu_3^* Fu_3 )
     1.45826395E-02 2 -37 24 # BR(hh_2 -> Hm_2 VWm^* )
     1.45826395E-02 2 37 -24 # BR(hh_2 -> Hm_2^* VWm )
DECAY 45 1.67364991E+01 # hh_3
# BR NDA ID1 ID2
     2.53389785E-01 2 36 23 # BR(hh_3 -> Ah_2 VZ )
     8.32856417E-04 2 -5 5 # BR(hh_3 -> Fd_3^* Fd_3 )
     6.43785048E-04 2 -6 6 # BR(hh_3 -> Fu_3^* Fu_3 )
     3.72520061E-01 2 -37 24 # BR(hh_3 -> Hm_2 VWm^* )
     3.72520061E-01 2 37 -24 # BR(hh_3 -> Hm_2^* VWm )
DECAY 36 8.81162963E-03 # Ah_2
# BR NDA ID1 ID2
     4.57204689E-03 2 21 21 # BR(Ah_2 -> VG VG )
     3.26986527E-04 2 -3 3 # BR(Ah_2 -> Fd_2^* Fd_2 )
     8.77534065E-01 2 -5 5 # BR(Ah_2 -> Fd_3^* Fd_3 )
     4.05196147E-04 2 -13 13 # BR(Ah_2 -> Fe_2^* Fe_2 )
     1.17090515E-01 2 -15 15 # BR(Ah_2 -> Fe_3^* Fe_3 )
DECAY 257 3.95368291E+01 # Ah_3
# BR NDA ID1 ID2
     2.04477925E-03 2 21 21 # BR(Ah_3 -> VG VG )
     2.27248609E-04 2 -5 5 # BR(Ah_3 -> Fd_3^* Fd_3 )
     6.71211251E-01 2 -6 6 # BR(Ah_3 -> Fu_3^* Fu_3 )
     1.54582209E-01 2 35 23 # BR(Ah_3 -> hh_2 VZ )
     8.11556658E-04 2 45 23 # BR(Ah_3 -> hh_3 VZ )
     5.70553377E-05 2 -37 24 # BR(Ah_3 -> Hm_2 VWm^* )
     5.70553377E-05 2 37 -24 # BR(Ah_3 -> Hm_2^* VWm )
     8.54834704E-02 2 -256 24 # BR(Ah_3 -> Hm_3 VWm^* )
     8.54834704E-02 2 256 -24 # BR(Ah_3 -> Hm_3^* VWm )
DECAY 37 1.99116194E-02 # Hm_2
# BR NDA ID1 ID2
     5.80290399E-04 2 4 -3 # BR(Hm_2^* -> Fu_2 Fd_2^* )
     1.47918189E-03 2 4 -5 # BR(Hm_2^* -> Fu_2 Fd_3^* )
     1.18217999E-03 2 6 -3 # BR(Hm_2^* -> Fu_3 Fd_2^* )
     9.47392423E-01 2 6 -5 # BR(Hm_2^* -> Fu_3 Fd_3^* )
     1.73752305E-04 2 14 -13 # BR(Hm_2^* -> Fv_2 Fe_2^* )
     4.91246011E-02 2 16 -15 # BR(Hm_2^* -> Fv_3 Fe_3^* )
DECAY 256 1.73767435E+01 # Hm_3
# BR NDA ID1 ID2
     1.67232914E-03 2 6 -3 # BR(Hm_3^* -> Fu_3 Fd_2^* )
     9.98150136E-01 2 6 -5 # BR(Hm_3^* -> Fu_3 Fd_3^* )
DECAY 21 0.00000000E+00 # VG
DECAY 22 0.00000000E+00 # VP
DECAY 23 2.49520000E+00 # VZ
DECAY 24 2.14100000E+00 # VWm
DECAY 12 0.00000000E+00 # Fv_1
DECAY 14 0.00000000E+00 # Fv_2
DECAY 16 0.00000000E+00 # Fv_3
DECAY1L 4 1.02568589E-23 # Fu_2
# BR NDA ID1 ID2
     9.77032512E-01 2 2 21 # BR(Fu_2 -> Fu_1 VG )
     2.29674885E-02 2 2 22 # BR(Fu_2 -> Fu_1 VP )
DECAY1L 6 1.38640396E+00 # Fu_3
# BR NDA ID1 ID2
     1.67437423E-03 2 3 24 # BR(Fu_3 -> Fd_2 VWm^* )
     9.98290221E-01 2 5 24 # BR(Fu_3 -> Fd_3 VWm^* )
DECAY1L 3 1.32817627E-20 # Fd_2
# BR NDA ID1 ID2
     9.94550004E-01 2 1 21 # BR(Fd_2 -> Fd_1 VG )
     5.44999648E-03 2 1 22 # BR(Fd_2 -> Fd_1 VP )
DECAY1L 5 5.44948073E-14 # Fd_3
# BR NDA ID1 ID2
     2.06993325E-02 2 1 21 # BR(Fd_3 -> Fd_1 VG )
     9.78274777E-01 2 3 21 # BR(Fd_3 -> Fd_2 VG )
     1.00473476E-03 2 3 22 # BR(Fd_3 -> Fd_2 VP )
DECAY1L 25 6.35979581E-03 # hh_1
# BR NDA ID1 ID2
     4.30862243E-04 2 -3 3 # BR(hh_1 -> Fd_2^* Fd_2 )
     8.55589232E-01 2 -5 5 # BR(hh_1 -> Fd_3^* Fd_3 )
     1.43856830E-04 2 -13 13 # BR(hh_1 -> Fe_2^* Fe_2 )
     4.06294113E-02 2 -15 15 # BR(hh_1 -> Fe_3^* Fe_3 )
     7.62126344E-02 2 -4 4 # BR(hh_1 -> Fu_2^* Fu_2 )
     3.01202873E-02 2 21 21 # BR(hh_1 -> VG VG )
     2.45420977E-03 2 22 22 # BR(hh_1 -> VP VP )
DECAY1L 35 4.36520429E+00 # hh_2
# BR NDA ID1 ID2
     4.36633886E-03 2 36 23 # BR(hh_2 -> Ah_2 VZ )
     2.67357397E-03 2 -5 5 # BR(hh_2 -> Fd_3^* Fd_3 )
     1.97404872E-04 2 -15 15 # BR(hh_2 -> Fe_3^* Fe_3 )
     2.68675014E-04 2 -4 4 # BR(hh_2 -> Fu_2^* Fu_2 )
     9.67685271E-01 2 -6 6 # BR(hh_2 -> Fu_3^* Fu_3 )
     9.88482625E-03 2 -37 24 # BR(hh_2 -> Hm_2 VWm^* )
     9.88482625E-03 2 37 -24 # BR(hh_2 -> Hm_2^* VWm )
     1.12715793E-04 2 25 23 # BR(hh_2 -> hh_1 VZ )
     4.77956982E-03 2 21 21 # BR(hh_2 -> VG VG )
DECAY1L 45 1.67267474E+01 # hh_3
# BR NDA ID1 ID2
     2.45486199E-01 2 36 23 # BR(hh_3 -> Ah_2 VZ )
     2.12392067E-03 2 -5 5 # BR(hh_3 -> Fd_3^* Fd_3 )
     9.00288469E-04 2 -6 6 # BR(hh_3 -> Fu_3^* Fu_3 )
     3.75682998E-01 2 -37 24 # BR(hh_3 -> Hm_2 VWm^* )
     3.75682998E-01 2 37 -24 # BR(hh_3 -> Hm_2^* VWm )
DECAY1L 36 2.19666039E-02 # Ah_2
# BR NDA ID1 ID2
     4.85070984E-04 2 -3 3 # BR(Ah_2 -> Fd_2^* Fd_2 )
     9.54036410E-01 2 -5 5 # BR(Ah_2 -> Fd_3^* Fd_3 )
     1.59785971E-04 2 -13 13 # BR(Ah_2 -> Fe_2^* Fe_2 )
     4.51772254E-02 2 -15 15 # BR(Ah_2 -> Fe_3^* Fe_3 )
     1.27307125E-04 2 -4 4 # BR(Ah_2 -> Fu_2^* Fu_2 )
DECAY1L 257 4.19658275E+01 # Ah_3
# BR NDA ID1 ID2
     4.85699604E-04 2 -5 5 # BR(Ah_3 -> Fd_3^* Fd_3 )
     6.96637737E-01 2 -6 6 # BR(Ah_3 -> Fu_3^* Fu_3 )
     1.42131464E-01 2 35 23 # BR(Ah_3 -> hh_2 VZ )
     7.68768273E-04 2 45 23 # BR(Ah_3 -> hh_3 VZ )
     5.09354809E-05 2 -37 24 # BR(Ah_3 -> Hm_2 VWm^* )
     5.09354809E-05 2 37 -24 # BR(Ah_3 -> Hm_2^* VWm )
     7.97621410E-02 2 -256 24 # BR(Ah_3 -> Hm_3 VWm^* )
     7.97621410E-02 2 256 -24 # BR(Ah_3 -> Hm_3^* VWm )
     2.13454609E-04 2 23 22 # BR(Ah_3 -> VZ VP )
DECAY1L 37 2.67885201E-02 # Hm_2
# BR NDA ID1 ID2
     4.79029007E-04 2 4 -3 # BR(Hm_2^* -> Fu_2 Fd_2^* )
     1.22849655E-03 2 4 -5 # BR(Hm_2^* -> Fu_2 Fd_3^* )
     1.17999627E-03 2 6 -3 # BR(Hm_2^* -> Fu_3 Fd_2^* )
     9.43393154E-01 2 6 -5 # BR(Hm_2^* -> Fu_3 Fd_3^* )
     2.28553877E-04 2 14 -13 # BR(Hm_2^* -> Fv_2 Fe_2^* )
     5.33706461E-02 2 16 -15 # BR(Hm_2^* -> Fv_3 Fe_3^* )
DECAY1L 256 1.92278914E+01 # Hm_3
# BR NDA ID1 ID2
     1.67134414E-03 2 6 -3 # BR(Hm_3^* -> Fu_3 Fd_2^* )
     9.97990087E-01 2 6 -5 # BR(Hm_3^* -> Fu_3 Fd_3^* )
     1.59074904E-04 2 24 22 # BR(Hm_3^* -> VWm^* VP )
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update to_full *
#*********************************************************************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#*********************************************************************
  10000 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, *
# 2=elastic photon of proton/ion beam *
# +/-3=PDF of electron/positron beam *
# +/-4=PDF of muon/antimuon beam *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     6500.0 = ebeam1 ! beam 1 total energy in GeV
     6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"

#*********************************************************************
# PDF CHOICE: this automatically fixes alpha_s and its evol. *
# pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
# iww=Improved Weizsaecker-Williams Approx.[hep-ph/9310350] *
# eva=Effective W/Z/A Approx. [2111.02442] *
# edff=EDFF in gamma-UPC [eq.(11) in 2207.03012] *
# chff=ChFF in gamma-UPC [eq.(13) in 2207.03012] *
# none=No PDF, same as lhapdf with lppx=0 *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
     False = fixed_fac_scale ! if .true. use fixed fac scale
 91.188 = scale ! fixed ren scale
 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales

#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
                             ! 0: sum over helicity, 1: importance sampling
   2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
                             ! 1 is old strategy (using amp square)
        ! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#*********************************************************************
# Customization (custom cuts/scale/bias/...) *
# list of files containing fortran function that overwrite default *
#*********************************************************************
  = custom_fcts ! List of files containing user hook function
#*******************************
# Parton level cuts definition *
#*******************************
  0.0 = dsqrt_shat ! minimal shat for full process
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 0.0 = ptb ! minimum pt for the b
 -1.0 = ptbmax ! maximum pt for the b
 {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
 -1.0 = etab ! max rap for the b
 0.0 = etabmin ! min rap for the b
 {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
 {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = ptheavy ! minimum pt for at least one heavy final state
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   False = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule

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

So the issue is with the block DECAY1L which I have seen in your card for the first time, and that MG5aMC has issue to handle.

I can certainly patch MG5aMC to ignore such block but is that the expected behaviour, or should I use the sum of the two decay block in that case (or use the sum for NLO generation and only the first for LO)?

I do not find any official reference to such block, do you know them?
If you do not, I can try to contact the authors of that tools to try to have an idea, and in the mean time you can either
merge those two decay block into a single one and/or remove them according to what you think is the best behaviour.

Cheers,

Olivier

Revision history for this message
Akshat Khanna (khanna-akshat) said : 		#4

So these are the Branching Ratios calculated at one loop and yes even I didn't find any official reference to the same in the manual. Maybe we can ignore them and calculate the production cross section at the tree level. But, even if I delete the same from the param_card.dat and try to launch the process, I get the following error:

generate_events run_01
Traceback (most recent call last):
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1543, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
    return func(arg, **opt)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2395, in do_generate_events
    switch_mode = self.ask_run_configuration(mode, args)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 6441, in ask_run_configuration
    self.ask_edit_cards(cards, plot=False, first_cmd=first_cmd)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 1039, in ask_edit_cards
    self.ask_edit_card_static(cards, mode, plot, self.options['timeout'],
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 1111, in ask_edit_card_static
    out = ask(question, '0', possible_answer, timeout=int(1.5*timeout),
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1148, in ask
    value = Cmd.timed_input(question, default, timeout=timeout,
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1759, in timed_input
    result = fct(question)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 2130, in __call__
    return self.cmdloop()
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 2298, in cmdloop
    super(SmartQuestion,self).cmdloop(intro)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 178, in cmdloop
    stop = self.postcmd(stop, line)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 6584, in postcmd
    self.do_update('dependent', timer=20)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 6634, in do_update
    self.update_dependent(self.mother_interface, self.me_dir, self.param_card,
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 6800, in update_dependent
    restrict_modify = param_card.update_dependent(model, restrict_card, log_level)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/models/check_param_card.py", line 487, in update_dependent
    parameters = model.set_parameters_and_couplings(self)
  File "/home/iitgn/Applications/MG5_aMC_v3_5_1/models/model_reader.py", line 172, in set_parameters_and_couplings
    raise MadGraph5Error(msg)
madgraph.MadGraph5Error: Invalid restriction card (not same block)
    {'effhiggscouplings', 'udlmix', 'decay', 'sphenolowenergy', 'yd', 'gaugegut', 'higgsfcc100', 'imyu', 'treelevelunitarity', 'uulmix', 'udrmix', 'treelevelunitaritywtrilinears', 'mass', 'higgslhc7', 'flavorkitqfv', 'mixedblock', 'uermix', 'pseudoscalarmix', 'minpar', 'imuulmix', 'higgslhc8', 'modsel', 'ye', 'gauge', 'sminputs', 'imyd', 'higgslhc13', 'scalarmix', 'higgslhc14', 'imye', 'hmix', 'uelmix', 'imuurmix', 'imuelmix', 'chargemix', 'imudrmix', 'imudlmix', 'imuermix', 'flavorkitlfv', 'yu', 'uurmix', 'imfwcoef', 'fwcoef', 'spheno'} != {'hmix', 'effhiggscouplings', 'uelmix', 'none', 'udrmix', 'uulmix', 'imuurmix', 'imuelmix', 'sminputs', 'mass', 'imudrmix', 'udlmix', 'decay', 'imudlmix', 'imuermix', 'uermix', 'uurmix', 'imuulmix'}.
    Missing block: none
    Unknown block : sphenolowenergy,yd,gaugegut,higgsfcc100,treelevelunitarity,imyu,treelevelunitaritywtrilinears,higgslhc7,flavorkitqfv,mixedblock,pseudoscalarmix,higgslhc8,modsel,ye,gauge,imyd,higgslhc13,yu,scalarmix,higgslhc14,imye,chargemix,flavorkitlfv,minpar,imfwcoef,fwcoef,spheno

I suppose, I need to delete more such blocks from the param card. I was using the same as the SPheno manual says that the SLHA format file can be used directly as the param card in Madgraph.

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

Hi,

I have contacted Sarah/Spheno author, and they gave me this paper as reference:
https://arxiv.org/abs/1703.09237
and advise that by default to just ignore the DECAY1L block.

In future version of MG5aMC, I will try to include a function within MG5aMC to be able to switch easily from the LO to the NLO width (and back) easily.

I spot now, your second issue. It seems an issue of the restriction card, could you send me the UFO model such that I can test that?
Not fully sure to know what the issue is here. Maybe this should be turn to a simple warning not fully sure...

Cheers,

Olivier

Revision history for this message
Akshat Khanna (khanna-akshat) said : 		#6

Hello,

Thank you for acknowledging the Decay issue.

Regarding the second issue, I think I have found the error. So, the BLOCK in the param_card.dat, in which the Lambda parameter values of the potential are numerically defined, is named as "none". While in the SPheno output file it is named as "MixedBlock". So, by changing the Block name and then using it as a param card works.

Thank you

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

Thanks for following up,

Having a block named "none" sounds an issue of the UFO model, maybe you can contact Sarah authors to check with them why the block is named none.

This assume that the issue is not on the MG5aMC side. Maybe to check that hyppothesis, one can check that "none" block is indeed used within the UFO file. If "none" is not present, then it would indicate a MG5aMC bug.

Cheers,

Olivier

Revision history for this message
Akshat Khanna (khanna-akshat) said : 		#8

Thank you for your suggestions. So, this bug was on the SARAH end, which now I have found and fixed.

Thank you

Can you help with this problem?

Provide an answer of your own, or ask Akshat Khanna for more information if necessary.

To post a message you must log in.