MadGraph5_aMC@NLO

process generation for models other than SM

Asked by nabanita ganguly

Hi,

I am trying to generate some simple processes for inert doublet model using mg5. I am using the UFO file generated by Feynrule.

Now suppose my process is p p > ~H+ ~H- where ~H+/- is a new scalar. It can decay to W+/- and ~A0 where ~A0 is stable.
I have set WHch in param_card to "Auto" so that mg5 calculates it's decay channels depending on the allowed interactions.

1) I run the command p p > ~H+ ~H- , (~H+ > w+ ~A0 , ~H- > w- ~A0) - PYTHIA is set to "on". result is as per expectation.

2) Next I try with p p > ~H+ ~H- - again with PYTHIA "on". Now as in the param_card the possible decay modes of ~H+ are present, I would expect that PYTHIA will decay ~H+ automatically , but instead ~H+ is being declared as a stable particle (i.e. status code = 1). This is a bit confusing for me as I have generated SUSY processes using mg5 earlier and in that case one really doesn't need to specify decay mode explicitly - it is done automatically by reading off the param_card provided by some SUSY generator.

Can you please suggest why is it not working in this case? is it mandatory to mention decay mode explicitly or I am doing something wrong?

Thanks,
Nabanita

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
Best Olivier Mattelaer (olivier-mattelaer) said : 		#1

Dear Nabanita,

What you did sounds correct. Looks like that this is a pythia problem, so i can not really comment on this.
Unfortunately Pythia6 is not supported anymore by their author. You can still try to contact them.

Cheers,

Olivier

> On Jun 28, 2016, at 08:13, nabanita ganguly <email address hidden> wrote:
>
> New question #295732 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/295732
>
> Hi,
>
> I am trying to generate some simple processes for inert doublet model using mg5. I am using the UFO file generated by Feynrule.
>
> Now suppose my process is p p > ~H+ ~H- where ~H+/- is a new scalar. It can decay to W+/- and ~A0 where ~A0 is stable.
> I have set WHch in param_card to "Auto" so that mg5 calculates it's decay channels depending on the allowed interactions.
>
> 1) I run the command p p > ~H+ ~H- , (~H+ > w+ ~A0 , ~H- > w- ~A0) - PYTHIA is set to "on". result is as per expectation.
>
> 2) Next I try with p p > ~H+ ~H- - again with PYTHIA "on". Now as in the param_card the possible decay modes of ~H+ are present, I would expect that PYTHIA will decay ~H+ automatically , but instead ~H+ is being declared as a stable particle (i.e. status code = 1). This is a bit confusing for me as I have generated SUSY processes using mg5 earlier and in that case one really doesn't need to specify decay mode explicitly - it is done automatically by reading off the param_card provided by some SUSY generator.
>
> Can you please suggest why is it not working in this case? is it mandatory to mention decay mode explicitly or I am doing something wrong?
>
> Thanks,
> Nabanita
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
nabanita ganguly (nabanita-rimpi) said : 		#2

Thanks Olivier Mattelaer, that solved my question.