Madgraph decay
Affects | Status | Importance | Assigned to | Milestone | |
---|---|---|---|---|---|
MadGraph5_aMC@NLO |
Fix Released
|
Undecided
|
Unassigned |
Bug Description
Dear MG team,
I am trying to generate events for process gg > h2 > h3 z* > h1 z* z* > bb jj mu+mu-, with h1 > b b~, h2 is. SM-like and mh1, mh3 < mh2 (2HDM type-1).
PS: z's are off-shell, where z* > mu+ mu- and z* > j j. The syntax suggested is : gg> h2 z > h3 mu+ mu-, ( h3 > z > h1 j j, h1 > b b~)
To generate g g > h2, the model is imported (2HDM type-I) , followed by adding the hgg_plugin (after changing the pdg of H from 25 to 35). I have checked the cross section computed by MG, it is \sigma (g g > h2) =17.63pb.
After that, I used MG to compute the cross-section of the process. I did some checks without/with decaying h1 to b b~ (Part 1/part2) ,
Part 1:
gg> h2 z > h3 mu+ mu-, ( h3 > z > h1 j j)
cross section = 3.57x10^-4 pb
Part 2 : after decaying h1 > b b~
gg> h2 z > h3 mu+ mu-, ( h3 > z > h1 j j, h1 > b b~)
cross section = 0.31pb
with, mh1 = 12.68 GeV, mh2 = 125.00 GeV, mh3 = 68.87 GeV, mhc = 135.00, tb = 6.3, sinbma = -0.16, wh1 = 4.680000e-09, wh2 = 4.520000e-03, wh3 = 1.130000e-04 , l2 = 0.262, l3 = -.592, lr7 = 0.00
The first value (part1) seems more reasonable since h2 is SM-like. After decaying h1 to b b~, the cross section become huge (part2),. It is worth to mention that the total width of h1 is set in the parameter card to 4.68x10^-9 and Br(h1 -> b b~) ~0.69.
Setting the total width of h1 to auto is not working, an error message is appearing. Please find below the error message:
INFO: failed to load model from pickle file. Try importing UFO from File
INFO: load particles
INFO: load vertices
INFO: Change particles name to pass to MG5 convention
INFO: syntax: compute_widths Particle [Particles] [OPTIONS]
INFO: -- Compute the widths for the particles specified.
INFO: By default, this takes the current param_card and overwrites it.
INFO: Precision allows to define when to include three/four/... body decays (LO).
INFO: If this number is an integer then all N-body decay will be included.
INFO: Various options:
INFO: --body_decay=X: Parameter to control the precision of the computation
INFO: if X is an integer, we compute all channels up to X-body decay.
INFO: if X <1, then we stop when the estimated error is lower than X.
INFO: if X >1 BUT not an integer, then we X = N + M, with M <1 and N an integer
INFO: We then either stop at the N-body decay or when the estimated error is lower than M.
INFO: default: 4.0025
INFO: --min_br=X: All channel which are estimated below this value will not be integrated numerically.
INFO: default: precision (decimal part of the body_decay options) divided by four
INFO: --precision_
INFO: default: 0.01
INFO: --path=X: path for param_card
INFO: default: take value from the model
INFO: --output=X: path where to write the resulting card.
INFO: default: overwrite input file. If no input file, write it in the model directory
INFO: --nlo: Compute NLO width [if the model support it]
ERROR: Invalid command: 25 is not a valid argument for compute_widths
INFO: Update the dependent parameter of the param_card.dat
Command "generate_events run_01" interrupted with error:
ValueError : invalid literal for int() with base 10: 'autoautoautoau
Please report this bug on https:/
More information is found in 'ME5_debug'.
Please attach this file to your report.
I did some checks:
- No problem after using the model file without hgg plugin and setting the width(s) to auto for h1 and/or h2 (good values)
No problem after using the model file with hgg plugin and setting the width to auto for h2 (good values)
Error message mentioned above, after using the model file with hgg plugin and setting the width to auto for h1.
Many thanks for any suggestions.
Kind regards
Changed in mg5amcnlo: | |
status: | Fix Committed → Fix Released |
Please find attached the model file