heavy higgs gluon fusion production with NMSSM specific decays

Since those are scalar, you can just use pythia for the decay.

Cheers,

Olivier
> On Jan 14, 2016, at 17:36, Matias Vazquez <email address hidden> wrote:
>
> New question #280980 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/280980
>
> Dear MG5_aMC team:
>
> I want to generate processes of the kind p p > h03 > h3 h3 > SM particles (h03 is the heaviest cp even higgs and h3 is the lightest pseudoscalar in the nmssm model). Of course, ggF production is not yet implemented in the nmssm model. Moreover, the higgs that I want to produce is heavy so i don't know how reliable an effective coupling would be... How could I generate such events in madgraph?
>
> At the moment I am trying with aMCSusHi (thanks to Marius Wiesemann!). It consists in ascript which generates a folder with the process p p > h using heft where the amplitudes in MadGraph are replaced by the ones from the external code SusHi. But, at this stage, I have no idea on how to decay the ggF produced Higgs through some specific channel, say h03 > h01 h01 or h03 > h3 h3.
>
> Thanks for your help!
> Matias.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

thanks for your answer. So is it possible to make a NMSSM-like Higgs to higgs decay at the Pythia level, even though madgraph didnt know about other scalars in the model?

So at this point I am not sure how to proceed. I have two approaches in mind:

a) I suitably modify the param_card and shower_card so they pass the information about QN, masses and BRs of the new states (not present in aMCSusHi) to Pythia6. I tried, without succeed, the following:

I add
Block mass
  25 6.20369408E+01 # lightest neutral scalar
 9000006 6.20369408E+01 # lightest neutral scalar
        35 1.22246098E+02 # second neutral scalar
 9000007 1.22246098E+02 # second neutral scalar
        45 2.56398486E+03 # third neutral scalar
 9000008 2.56398486E+03 # third neutral scalar
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.38202936E+00 # Top Quark
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 5.67548338E-05 # Lightest neutral Higgs scalar
DECAY 9000006 5.67548338E-05 # Lightest neutral Higgs scalar
DECAY 35 4.09421269E-03 # 2nd neutral Higgs scalar
DECAY 9000007 4.09421269E-03 # 2nd neutral Higgs scalar
DECAY 45 6.60335016E+01 # 3rd neutral Higgs scalar
DECAY 9000008 6.60335016E+01 # 3rd neutral Higgs scalar
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================
Block QNUMBERS 9000006 # h1
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000007 # 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)
Block QNUMBERS 9000008 # h03
        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)

in the param_card, and in the shower_card I simply add a line indicating the decaying of the heavy Higgs H03 > h2 h2. Note that I am producing the heavy Higgs H03, but since the model heft only has one Higgs (which amplitude is computed by SusHi), I indicate the decay with the PDGID 25:
 #***********************************************************************
# MadGraph5_aMC@NLO *
# *
# shower_card.dat aMC@NLO *
# *
# This file is used to set the parameters for the shower. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: variable = value # comment *
#***********************************************************************
#
#****************
# Shower settings
#****************
#
#***********************************************************************
# Number of events, jobs, errors, and random seeds *
#***********************************************************************
nevents = -1 # N evts to shower (< 0 = all)
nsplit_jobs = 1 # N jobs to run in parallel (< 100!!)
combine_td = T # combine the topdrawer files if nsplit_jobs > 1
maxprint = 5 # N evts to print in the log
maxerrs = 0.1 # max fraction of errors
rnd_seed = 0 # 1st random seed (0 = default)
rnd_seed2 = 0 # 2nd random seed (0 = default) !ONLY FOR HWERIG6!
#***********************************************************************
# PDFs and non-perturbative modelling *
#***********************************************************************
pdfcode = 1 # 0 = internal, 1 = same as NLO, other = lhaglue
ue_enabled = F # underlying event
hadronize = T # hadronisation on/off !IGNORED BY HERWIG6!
lambda_5 = -1 # Lambda_5 (< 0 = default) !IGNORED BY PYTHIA8!
#***********************************************************************
# Stable or unstable particles *
#***********************************************************************
b_stable = F # set B hadrons stable
pi_stable = T # set pi0's stable
wp_stable = F # set w+'s stable
wm_stable = F # set w-'s stable
z_stable = F # set z0's stable
h_stable = F # set Higgs' stable
tap_stable = T # set tau+'s stable
tam_stable = T # set tau-'s stable
mup_stable = F # set mu+'s stable
mum_stable = F # set mu-'s stable
#***********************************************************************
# Mass of the b quark *
#***********************************************************************
b_mass = -1 # if < 0 = read from SubProcesses/MCmasses_*.inc
#***********************************************************************
# Special settings *
#***********************************************************************
is_4lep = F # T if 4-lepton production !ONLY FOR PYTHIA6!
is_bbar = F # T if bb~ production !ONLY FOR HERWIG6!
#***********************************************************************
# FxFx merging parameters !ONLY FOR PYTHIA8!
#***********************************************************************
Qcut = -1.0 # Merging scale
njmax = 0 # Maximal multiplicity in the merging
#***********************************************************************
# Decay channels *
# Write down the decay channels for the resonances, to be performed by *
# the shower. *
# The syntax (for a two-body decay) is *
# DM_I = M > D1 D2 @ BR @ ME *
# where I < 100, M is the decaying resonance, D1, D2 are the decay *
# products (up to D5 if such a decay is supported by the shower), BR *
# is the branching ratio (only used by the HERWIG6 shower, ignored *
# otherwise) and ME is the type of matrix element to be used in the *
# decay (only used by HERWIG6, ignored otherwise). *
# BR's are correctly understood by HERWIG6 only if they add up to 1 *
# and only if no more than three modes are required for a given *
# resonance. *
# ME corresponds to the third entry of subroutine HWMODK, see the *
# relevant manual. *
# *
# WARNING: in HERWIG6, the order of decay products in > 2-body decays *
# IS RELEVANT. *
# WARNING: in PYTHIA6, turning hadronisation off disables top decays *
# WARNING: in PYTHIA6 and PYTHIA8, 1 -> n decays (with n > 2) are *
# handled through a sequence of 1 -> 2 decays. *
# *
# Examples of syntax: *
# Z -> e+ e- or mu+ mu- with BR = 0.5 each *
# DM_1 = 23 > -11 11 @ 0.5d0 @ 100 *
# DM_2 = 23 > -13 13 @ 0.5d0 @ 100 *
# H -> tau+ tau- with BR = 1 *
# DM_3 = 25 > -15 15 @ 1.0d0 @ 0
# t -> nu_e e+ b with BR = 1 (HERWIG) *
# DM_4 = 6 > 12 -11 5 @ 1d0 @ 100 *
# t -> nu_e e+ b with BR = 1 (PYTHIA) *
# DM_5 = 6 > 24 5 @ 1d0 @ 100 *
# DM_6 = 24 > 12 -11 @ 1d0 @ 100
#####################################
DM_7 = 25 > 9000007 9000007 @ 1d0 @ 0 *
#***********************************************************************
# Extra Libraries/analyses *
# The following lines need to be changed if the user does not want to *
# create a StdHEP/HepMC file, but to directly run an own analysis (to *
# be placed in HWAnalyzer or analogous MCatNLO subfolders). *
# Please use files in those folders as examples. *
#***********************************************************************
EXTRALIBS = stdhep Fmcfio # Needed extra-libraries (not LHAPDF).
                                    # Default: "stdhep Fmcfio"
EXTRAPATHS = ../lib
                                    # Path to the extra-libraries.
                                    # Default: "../lib"
INCLUDEPATHS=
                                    # Path to the dirs containing header files neede by C++.
                                    # Directory names are separated by white spaces
ANALYSE =
                                    # User's analysis and histogramming routines
                                    # (please use .o as extension and use spaces to separate files)

#
Of course, this naiv approach did not work. Is there any way to use the cards for instructing Pyhia6?

b) I generate the Higgs events through aMCSusHi, so I will obtain a events.lhe file. I write some Pythia6 code such that it teaches Pythia about the new scalars and its QN, BR and masses, and it showers the events.lhe file. Question: Is it possible to do the showering using Pythia6 with new states which were not present when generating the events.lhe file with madgraph(aMCSuHi)?

Cheers,
Matias.

Hi,

I’m not a pythia expert so I can only give you my opinion on this.
I can not ensure that what I will wrote is 100% correct.

Your method a) should work if you include the Branching ratio table inside the param_card such that pythia6/8 knows how the decay should be done.
(you can put fake one in principle to force the decay that you want but then you will need to change the cross-section accordingly.)

For your method b. This is going to work as well.
I do not understand your question
> Question: Is it possible to do the showering using
> Pythia6 with new states which were not present when generating the
> events.he file with madgraph(aMCSuHi)?

If they are no particle in the events file, then the shower will not apply on those particles (obviously).
If you mean that you need to change the identity of the particle (pdg code) between your aMCSuHi.
Then I do not see the problem.

Cheers,

Olivier

> On Jan 19, 2016, at 12:42, Matias Vazquez <email address hidden> wrote:
>
> Question #280980 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/280980
>
> Status: Answered => Open
>
> Matias Vazquez is still having a problem:
> Hi Olivier,
>
> thanks for your answer. So is it possible to make a NMSSM-like Higgs to
> higgs decay at the Pythia level, even though madgraph didnt know about
> other scalars in the model?
>
> So at this point I am not sure how to proceed. I have two approaches in
> mind:
>
> a) I suitably modify the param_card and shower_card so they pass the
> information about QN, masses and BRs of the new states (not present in
> aMCSusHi) to Pythia6. I tried, without succeed, the following:
>
> I add
> Block mass
> 25 6.20369408E+01 # lightest neutral scalar
> 9000006 6.20369408E+01 # lightest neutral scalar
> 35 1.22246098E+02 # second neutral scalar
> 9000007 1.22246098E+02 # second neutral scalar
> 45 2.56398486E+03 # third neutral scalar
> 9000008 2.56398486E+03 # third neutral scalar
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.38202936E+00 # Top Quark
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 5.67548338E-05 # Lightest neutral Higgs scalar
> DECAY 9000006 5.67548338E-05 # Lightest neutral Higgs scalar
> DECAY 35 4.09421269E-03 # 2nd neutral Higgs scalar
> DECAY 9000007 4.09421269E-03 # 2nd neutral Higgs scalar
> DECAY 45 6.60335016E+01 # 3rd neutral Higgs scalar
> DECAY 9000008 6.60335016E+01 # 3rd neutral Higgs scalar
> #===========================================================
> # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
> #===========================================================
> Block QNUMBERS 9000006 # h1
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2S+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # Particle/Antiparticle distinction (0=own anti)
> Block QNUMBERS 9000007 # 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)
> Block QNUMBERS 9000008 # h03
> 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)
>
> in the param_card, and in the shower_card I simply add a line indicating the decaying of the heavy Higgs H03 > h2 h2. Note that I am producing the heavy Higgs H03, but since the model heft only has one Higgs (which amplitude is computed by SusHi), I indicate the decay with the PDGID 25:
> #***********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # shower_card.dat aMC@NLO *
> # *
> # This file is used to set the parameters for the shower. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a hash (#) are info or comments *
> # *
> # mind the format: variable = value # comment *
> #***********************************************************************
> #
> #****************
> # Shower settings
> #****************
> #
> #***********************************************************************
> # Number of events, jobs, errors, and random seeds *
> #***********************************************************************
> nevents = -1 # N evts to shower (< 0 = all)
> nsplit_jobs = 1 # N jobs to run in parallel (< 100!!)
> combine_td = T # combine the topdrawer files if nsplit_jobs > 1
> maxprint = 5 # N evts to print in the log
> maxerrs = 0.1 # max fraction of errors
> rnd_seed = 0 # 1st random seed (0 = default)
> rnd_seed2 = 0 # 2nd random seed (0 = default) !ONLY FOR HWERIG6!
> #***********************************************************************
> # PDFs and non-perturbative modelling *
> #***********************************************************************
> pdfcode = 1 # 0 = internal, 1 = same as NLO, other = lhaglue
> ue_enabled = F # underlying event
> hadronize = T # hadronisation on/off !IGNORED BY HERWIG6!
> lambda_5 = -1 # Lambda_5 (< 0 = default) !IGNORED BY PYTHIA8!
> #***********************************************************************
> # Stable or unstable particles *
> #***********************************************************************
> b_stable = F # set B hadrons stable
> pi_stable = T # set pi0's stable
> wp_stable = F # set w+'s stable
> wm_stable = F # set w-'s stable
> z_stable = F # set z0's stable
> h_stable = F # set Higgs' stable
> tap_stable = T # set tau+'s stable
> tam_stable = T # set tau-'s stable
> mup_stable = F # set mu+'s stable
> mum_stable = F # set mu-'s stable
> #***********************************************************************
> # Mass of the b quark *
> #***********************************************************************
> b_mass = -1 # if < 0 = read from SubProcesses/MCmasses_*.inc
> #***********************************************************************
> # Special settings *
> #***********************************************************************
> is_4lep = F # T if 4-lepton production !ONLY FOR PYTHIA6!
> is_bbar = F # T if bb~ production !ONLY FOR HERWIG6!
> #***********************************************************************
> # FxFx merging parameters !ONLY FOR PYTHIA8!
> #***********************************************************************
> Qcut = -1.0 # Merging scale
> njmax = 0 # Maximal multiplicity in the merging
> #***********************************************************************
> # Decay channels *
> # Write down the decay channels for the resonances, to be performed by *
> # the shower. *
> # The syntax (for a two-body decay) is *
> # DM_I = M > D1 D2 @ BR @ ME *
> # where I < 100, M is the decaying resonance, D1, D2 are the decay *
> # products (up to D5 if such a decay is supported by the shower), BR *
> # is the branching ratio (only used by the HERWIG6 shower, ignored *
> # otherwise) and ME is the type of matrix element to be used in the *
> # decay (only used by HERWIG6, ignored otherwise). *
> # BR's are correctly understood by HERWIG6 only if they add up to 1 *
> # and only if no more than three modes are required for a given *
> # resonance. *
> # ME corresponds to the third entry of subroutine HWMODK, see the *
> # relevant manual. *
> # *
> # WARNING: in HERWIG6, the order of decay products in > 2-body decays *
> # IS RELEVANT. *
> # WARNING: in PYTHIA6, turning hadronisation off disables top decays *
> # WARNING: in PYTHIA6 and PYTHIA8, 1 -> n decays (with n > 2) are *
> # handled through a sequence of 1 -> 2 decays. *
> # *
> # Examples of syntax: *
> # Z -> e+ e- or mu+ mu- with BR = 0.5 each *
> # DM_1 = 23 > -11 11 @ 0.5d0 @ 100 *
> # DM_2 = 23 > -13 13 @ 0.5d0 @ 100 *
> # H -> tau+ tau- with BR = 1 *
> # DM_3 = 25 > -15 15 @ 1.0d0 @ 0
> # t -> nu_e e+ b with BR = 1 (HERWIG) *
> # DM_4 = 6 > 12 -11 5 @ 1d0 @ 100 *
> # t -> nu_e e+ b with BR = 1 (PYTHIA) *
> # DM_5 = 6 > 24 5 @ 1d0 @ 100 *
> # DM_6 = 24 > 12 -11 @ 1d0 @ 100
> #####################################
> DM_7 = 25 > 9000007 9000007 @ 1d0 @ 0 *
> #***********************************************************************
> # Extra Libraries/analyses *
> # The following lines need to be changed if the user does not want to *
> # create a StdHEP/HepMC file, but to directly run an own analysis (to *
> # be placed in HWAnalyzer or analogous MCatNLO subfolders). *
> # Please use files in those folders as examples. *
> #***********************************************************************
> EXTRALIBS = stdhep Fmcfio # Needed extra-libraries (not LHAPDF).
> # Default: "stdhep Fmcfio"
> EXTRAPATHS = ../lib
> # Path to the extra-libraries.
> # Default: "../lib"
> INCLUDEPATHS=
> # Path to the dirs containing header files neede by C++.
> # Directory names are separated by white spaces
> ANALYSE =
> # User's analysis and histogramming routines
> # (please use .o as extension and use spaces to separate files)
>
> #
> Of course, this naiv approach did not work. Is there any way to use the cards for instructing Pyhia6?
>
> b) I generate the Higgs events through aMCSusHi, so I will obtain a
> events.lhe file. I write some Pythia6 code such that it teaches Pythia
> about the new scalars and its QN, BR and masses, and it showers the
> events.lhe file. Question: Is it possible to do the showering using
> Pythia6 with new states which were not present when generating the
> events.lhe file with madgraph(aMCSuHi)?
>
> Cheers,
> Matias.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

The problem of modifying the cards is that when I run in the ggF process folder ./bin/generate_events, pythia seems to ignore the decays indicated in the param_card, but only reads those that I write in shower_card using the format specified in there. So if, for instance, I write

DM_1 = 25 > 9000007 9000007 @ 1d0 @ 0

where 9000007 was defined in the param_card, I find the following error in the pyhia log file:

 Error type 8 has occured after 0 PYEXEC calls:
     (PYONOF:) could not interpret particle code ****
 How are the events normalized ("average" or "sum")?
 stdhep initialized

     Advisory warning type 3 given after 0 PYEXEC calls:
     (PYRESD:) Failed to decay particle 25 with mass 2563.985

So clearly the inputs that I set in the param_card are being ignored. I tried setting a decay of the like
DECAY 25 5.67548338E-05 # Lightest neutral Higgs scalar
     0.50000000E+00 2 15 -15 # BR(H_3 -> H_2 H_2)
     0.50000000E+00 2 15 -15 # BR(H_3 -> H_2 H_2)

in the param_card and setting off any decay in the shower card: then Pythia ignores such decays and it might use some default ones since it is decaying my heavy Higgs into t tbar. On the other hand, if I set a certain SM-like decay in the shower_card:
DM_3 = 25 > -15 15 @ 1.0d0 @ 0

Any idea on how to tell the program to use the values in param_card instead of those in shower_card? or, alternatively, How can I instruct Pythia about the existence of other scalars?
It does work. I am just trying to avoid option B because it seems rather co

sorry, forgot the last line saying "It does work. I am just trying to avoid option B because it seems rather co"

HI,

You can not use the PDG=25 since this is reserved for the SM Higgs.
This means that Pythia will use the internal decay by decay and not the one specify in the cards.

Cheers,

Olivier

> On Jan 19, 2016, at 14:52, Matias Vazquez <email address hidden> wrote:
>
> Question #280980 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/280980
>
> Status: Answered => Open
>
> Matias Vazquez is still having a problem:
> Hi,
>
> The problem of modifying the cards is that when I run in the ggF process
> folder ./bin/generate_events, pythia seems to ignore the decays
> indicated in the param_card, but only reads those that I write in
> shower_card using the format specified in there. So if, for instance, I
> write
>
> DM_1 = 25 > 9000007 9000007 @ 1d0 @ 0
>
> where 9000007 was defined in the param_card, I find the following error in the pyhia log file:
>
> Error type 8 has occured after 0 PYEXEC calls:
> (PYONOF:) could not interpret particle code ****
> How are the events normalized ("average" or "sum")?
> stdhep initialized
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle 25 with mass 2563.985
>
> So clearly the inputs that I set in the param_card are being ignored. I tried setting a decay of the like
> DECAY 25 5.67548338E-05 # Lightest neutral Higgs scalar
> 0.50000000E+00 2 15 -15 # BR(H_3 -> H_2 H_2)
> 0.50000000E+00 2 15 -15 # BR(H_3 -> H_2 H_2)
>
> in the param_card and setting off any decay in the shower card: then Pythia ignores such decays and it might use some default ones since it is decaying my heavy Higgs into t tbar. On the other hand, if I set a certain SM-like decay in the shower_card:
> DM_3 = 25 > -15 15 @ 1.0d0 @ 0
>
> Any idea on how to tell the program to use the values in param_card instead of those in shower_card? or, alternatively, How can I instruct Pythia about the existence of other scalars?
> It does work. I am just trying to avoid option B because it seems rather co
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.