MadGraph5_aMC@NLO

Decaying Higgs with Madspin

Asked by Aya Ismail Beshr

Hi,

I am using 2HDM typeI UFO model generated by SARAH. I have an issue with madspin, it is giving an error BranchingRatio is larger than one for particle 25, I did the following:

import model 2HDM --model name
define j = g u c d s b u~ c~ d~ s~ b~
define h = h1
generate p p > h2 > h1 h
output
launch

I edited the run card in part of max_jet_flavour to be 5 instead of 4.

In the param_card the decay width for sm like higgs h1 is 3.892474e-03.

In madspin card I did:

set spinmode none

decay h1 > b b~
decay h > z z, z > mu+ mu-, z > mu+ mu-

Is this true? bec. I guess I am making something wrong.

I have added the run, param and spin cards in the comment below.

Thanks a lot!
Aya

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Aya Ismail Beshr (ayabeshr) said : 		#1

I need that one of 2 sm like higgses to go to z bosons

Revision history for this message
Aya Ismail Beshr (ayabeshr) said : 		#2

************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#* Manual: *
#* cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
 set spinmode none # Use one of the madspin special mode
 set max_weight_ps_point 400 # number of PS to estimate the maximum for each event

# specify the decay for the final state particles
decay h1 > b b~
decay h > z z, z > mu+ mu-, z > mu+ mu-

# running the actual code
launch

Revision history for this message
Aya Ismail Beshr (ayabeshr) said : 		#3

######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR EFFHIGGSCOUPLINGS
###################################
BLOCK EFFHIGGSCOUPLINGS #
      25 21 21 6.614119e-05 # hgg1
      25 22 22 3.181640e-05 # hpp1
      35 21 21 1.923997e-07 # hgg2
      35 22 22 9.778655e-07 # hpp2
      36 21 21 2.109124e-06 # agg2
      36 22 22 2.786228e-07 # app2
###################################
## INFORMATION FOR HMIX
###################################
BLOCK HMIX #
      10 1.550799e+00 # betah
      11 3.138687e+00 # alphah
      31 1.000000e-01 # rlam1
      32 1.290650e-01 # rlam2
      33 1.100000e+00 # rlam3
      34 -5.000000e-01 # rlam4
      35 5.000000e-01 # rlam5
###################################
## INFORMATION FOR IMHMIX
###################################
BLOCK IMHMIX #
      31 0.000000e+00 # ilam1
      32 0.000000e+00 # ilam2
      33 0.000000e+00 # ilam3
      34 0.000000e+00 # ilam4
      35 0.000000e+00 # ilam5
###################################
## INFORMATION FOR IMUDLMIX
###################################
BLOCK IMUDLMIX #
      1 1 0.000000e+00 # izdl11
      1 2 0.000000e+00 # izdl12
      1 3 0.000000e+00 # izdl13
      2 1 0.000000e+00 # izdl21
      2 2 0.000000e+00 # izdl22
      2 3 0.000000e+00 # izdl23
      3 1 0.000000e+00 # izdl31
      3 2 0.000000e+00 # izdl32
      3 3 0.000000e+00 # izdl33
###################################
## INFORMATION FOR IMUDRMIX
###################################
BLOCK IMUDRMIX #
      1 1 0.000000e+00 # izdr11
      1 2 0.000000e+00 # izdr12
      1 3 0.000000e+00 # izdr13
      2 1 0.000000e+00 # izdr21
      2 2 0.000000e+00 # izdr22
      2 3 0.000000e+00 # izdr23
      3 1 0.000000e+00 # izdr31
      3 2 0.000000e+00 # izdr32
      3 3 0.000000e+00 # izdr33
###################################
## INFORMATION FOR IMUELMIX
###################################
BLOCK IMUELMIX #
      1 1 0.000000e+00 # izel11
      1 2 0.000000e+00 # izel12
      1 3 0.000000e+00 # izel13
      2 1 0.000000e+00 # izel21
      2 2 0.000000e+00 # izel22
      2 3 0.000000e+00 # izel23
      3 1 0.000000e+00 # izel31
      3 2 0.000000e+00 # izel32
      3 3 0.000000e+00 # izel33
###################################
## INFORMATION FOR IMUERMIX
###################################
BLOCK IMUERMIX #
      1 1 0.000000e+00 # izer11
      1 2 0.000000e+00 # izer12
      1 3 0.000000e+00 # izer13
      2 1 0.000000e+00 # izer21
      2 2 0.000000e+00 # izer22
      2 3 0.000000e+00 # izer23
      3 1 0.000000e+00 # izer31
      3 2 0.000000e+00 # izer32
      3 3 0.000000e+00 # izer33
###################################
## INFORMATION FOR IMUULMIX
###################################
BLOCK IMUULMIX #
      1 1 0.000000e+00 # izul11
      1 2 0.000000e+00 # izul12
      1 3 0.000000e+00 # izul13
      2 1 0.000000e+00 # izul21
      2 2 0.000000e+00 # izul22
      2 3 0.000000e+00 # izul23
      3 1 0.000000e+00 # izul31
      3 2 0.000000e+00 # izul32
      3 3 0.000000e+00 # izul33
###################################
## INFORMATION FOR IMUURMIX
###################################
BLOCK IMUURMIX #
      1 1 0.000000e+00 # izur11
      1 2 0.000000e+00 # izur12
      1 3 0.000000e+00 # izur13
      2 1 0.000000e+00 # izur21
      2 2 0.000000e+00 # izur22
      2 3 0.000000e+00 # izur23
      3 1 0.000000e+00 # izur31
      3 2 0.000000e+00 # izur32
      3 3 0.000000e+00 # izur33
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      1 5.000000e-03 # md1
      2 2.500000e-03 # mu1
      3 9.500000e-02 # md2
      4 1.270000e+00 # mu2
      5 4.180000e+00 # md3
      6 1.735000e+02 # mu3
      11 5.109989e-04 # me1
      13 1.056584e-01 # me2
      15 1.776690e+00 # me3
      23 9.118870e+01 # mz
      25 1.251658e+02 # mh1
      35 2.453508e+02 # mh2
      36 1.729265e+02 # mah2
      37 2.453892e+02 # mhm2
      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.094009e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - 8.0349727e+01
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.370360e+02 # aewm1
      2 1.166370e-05 # gf
      3 1.187000e-01 # as
###################################
## INFORMATION FOR UDLMIX
###################################
BLOCK UDLMIX #
      1 1 1.000000e+00 # rzdl11
      1 2 0.000000e+00 # rzdl12
      1 3 0.000000e+00 # rzdl13
      2 1 0.000000e+00 # rzdl21
      2 2 1.000000e+00 # rzdl22
      2 3 0.000000e+00 # rzdl23
      3 1 0.000000e+00 # rzdl31
      3 2 0.000000e+00 # rzdl32
      3 3 1.000000e+00 # rzdl33
###################################
## INFORMATION FOR UDRMIX
###################################
BLOCK UDRMIX #
      1 1 1.000000e+00 # rzdr11
      1 2 0.000000e+00 # rzdr12
      1 3 0.000000e+00 # rzdr13
      2 1 0.000000e+00 # rzdr21
      2 2 1.000000e+00 # rzdr22
      2 3 0.000000e+00 # rzdr23
      3 1 0.000000e+00 # rzdr31
      3 2 0.000000e+00 # rzdr32
      3 3 1.000000e+00 # rzdr33
###################################
## INFORMATION FOR UELMIX
###################################
BLOCK UELMIX #
      1 1 1.000000e+00 # rzel11
      1 2 0.000000e+00 # rzel12
      1 3 0.000000e+00 # rzel13
      2 1 0.000000e+00 # rzel21
      2 2 1.000000e+00 # rzel22
      2 3 0.000000e+00 # rzel23
      3 1 0.000000e+00 # rzel31
      3 2 0.000000e+00 # rzel32
      3 3 1.000000e+00 # rzel33
###################################
## INFORMATION FOR UERMIX
###################################
BLOCK UERMIX #
      1 1 1.000000e+00 # rzer11
      1 2 0.000000e+00 # rzer12
      1 3 0.000000e+00 # rzer13
      2 1 0.000000e+00 # rzer21
      2 2 1.000000e+00 # rzer22
      2 3 0.000000e+00 # rzer23
      3 1 0.000000e+00 # rzer31
      3 2 0.000000e+00 # rzer32
      3 3 1.000000e+00 # rzer33
###################################
## INFORMATION FOR UULMIX
###################################
BLOCK UULMIX #
      1 1 -1.000000e+00 # rzul11
      1 2 0.000000e+00 # rzul12
      1 3 0.000000e+00 # rzul13
      2 1 0.000000e+00 # rzul21
      2 2 -1.000000e+00 # rzul22
      2 3 0.000000e+00 # rzul23
      3 1 0.000000e+00 # rzul31
      3 2 0.000000e+00 # rzul32
      3 3 -1.000000e+00 # rzul33
###################################
## INFORMATION FOR UURMIX
###################################
BLOCK UURMIX #
      1 1 1.000000e+00 # rzur11
      1 2 0.000000e+00 # rzur12
      1 3 0.000000e+00 # rzur13
      2 1 0.000000e+00 # rzur21
      2 2 1.000000e+00 # rzur22
      2 3 0.000000e+00 # rzur23
      3 1 0.000000e+00 # rzur31
      3 2 0.000000e+00 # rzur32
      3 3 1.000000e+00 # rzur33
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.555653e+00 # wu3
DECAY 23 2.495200e+00 # wz
DECAY 24 2.141000e+00 # wwm
DECAY 25 3.892474e-03 # wh1
DECAY 35 2.109076e-03 # wh2
DECAY 36 2.286365e-06 # wah2
DECAY 37 1.863061e-03 # whm2
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
###################################
## 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)
###################################
## INFORMATION FOR QNUMBERS 36
###################################
BLOCK QNUMBERS 36 # h3
      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)

Revision history for this message
Aya Ismail Beshr (ayabeshr) said : 		#4

#*********************************************************************
# 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 full_run_card *
#*********************************************************************
#
#*********************************************************************
# 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 *
#*********************************************************************
  10 = 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=photon from proton, *
# 3=photon from electron *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     7000.0 = ebeam1 ! beam 1 total energy in GeV
     7000.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     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"

#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# 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) *
#*********************************************************************
 {} = 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) *
#*********************************************************************
 {} = 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) *
#*********************************************************************
 5 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = 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
# Syscalc is deprecated but to see the associate options type'update syscalc'

1e-12 = small_width_treatment

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

The first issue is the syntax:
decay h > z z, z > mu+ mu-, z > mu+ mu-
You can not ask an onshell Higgs to decay to two onshell Z.
decay h > z > mu+ mu- mu+ mu-

is already better.

Second,
define h = h1
generate p p > h2 > h1 h
is a weird syntax and you might create issue (or false hope) with that attempt.

generate p p > h2 > h h
and
decay h > b b~
decay h > z > mu+ mu- mu+ mu-

Should work better.

Third, The automatic computation of the branching ratio for the Higgs is not working (see 1402.1178)
So for the Higgs case it is advise to not use that tool (since it misses loop-induced decay)

Cheers,

Olivier

Revision history for this message
Aya Ismail Beshr (ayabeshr) said : 		#6

Hi OLivier,

Thanks a lot for the illustration!

I did as you kindly mentioned and everything works fine!

Thanks,
Aya

Revision history for this message
Aya Ismail Beshr (ayabeshr) said : 		#7

Thanks Olivier Mattelaer, that solved my question.