Using parton level lhe files for pythia 8

Asked by Sebastian Wild

Hello,

I'm a new user, and I'd like to use MadGraph5 in order to produce parton level lhe files, which then can be read in by pythia8 for the showering.
As a first (random) example, I used the 2HDM model, generated the process p p > h+ h- , and produced lhe files, via

import model 2HDM
generate p p > h+ h-
output
launch

,using the default cards. Then, I feed the lhe file "unweighted_events.lhe" into pythia 8, however it doesn't work: Pythia aborts by saying "no open decay channel for id = 37" (which is the h+ particle). Indeed, in the lhe file there is no information contained about the decay channels of h+.

What is the correct way of telling pythia the decay channels of the new particles? I guess via the slha interface, but I couldn't find out how...

Thanks a lot!

Question information

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

Hi,

For such simple process process, you can either make the decay of the h+ inside MadGraph.
via the following syntax:
generate p p > h + h-, h+ > ta+ vt, h- > ta- vt~

The other options is to have create the decay table for pythia.
the easiest way to do that, is to put the H+ width to “Auto” in the param_card (related paper: 1405.0301)
Which computes automatically the total width and the main partial width (see the paper for more details)

Cheers,

Olivier

On 14 Apr 2015, at 15:41, Sebastian Wild <email address hidden> wrote:

> New question #265178 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/265178
>
> Hello,
>
> I'm a new user, and I'd like to use MadGraph5 in order to produce parton level lhe files, which then can be read in by pythia8 for the showering.
> As a first (random) example, I used the 2HDM model, generated the process p p > h+ h- , and produced lhe files, via
>
> import model 2HDM
> generate p p > h+ h-
> output
> launch
>
> ,using the default cards. Then, I feed the lhe file "unweighted_events.lhe" into pythia 8, however it doesn't work: Pythia aborts by saying "no open decay channel for id = 37" (which is the h+ particle). Indeed, in the lhe file there is no information contained about the decay channels of h+.
>
> What is the correct way of telling pythia the decay channels of the new particles? I guess via the slha interface, but I couldn't find out how...
>
> Thanks a lot!
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Sebastian Wild (sebastian-wild-h) said :
#2

Thanks a lot! I tried both options, and both are working.