MadGraph5_aMC@NLO

wrong status for BSM particles after showering

Asked by Renjie Wang

Dear experts,

I have a quick question about particle status in LHE file after run with MG+pythia.

The question is I have a stable Dirac fermion particle (1009,-1009), when just run MG, the everything is good, but after parton shower with pythia, the particle goes to the intermediator status (2), not stable status (1).

How can I solve it?

Thanks ahead;
Ren-Jie

<event>
 13 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
      -2 -1 0 0 0 0 -0.85457660E+02 -0.12572625E+01 0.13608267E+03 0.16069563E+03 0.00000000E+00 0. 0.
       2 -1 0 0 0 0 0.66168770E+00 0.15730045E+01 -0.16887368E+04 0.16887377E+04 0.00000000E+00 0. 0.
      23 2 1 1 0 0 -0.41673835E+02 0.11808129E+03 0.30516587E+01 0.15537687E+03 0.91938886E+02 0. 0.
    1008 2 1 1 0 0 -0.43113524E+02 -0.11776541E+03 -0.15557334E+04 0.16940541E+04 0.65862396E+03 0. 0.
     -13 2 3 3 0 0 -0.60752329E+02 0.39187104E+02 -0.10603652E+02 0.73067936E+02 0.10566000E+00 0. 0.
      13 2 3 3 0 0 0.19078494E+02 0.78894190E+02 0.13655311E+02 0.82308935E+02 0.10566000E+00 0. 0.
    1009 2 4 4 0 0 0.16637392E+03 -0.12024735E+03 -0.34185607E+03 0.44610008E+03 0.20000000E+03 0. 0.
   -1009 2 4 4 0 0 -0.20948745E+03 0.24819364E+01 -0.12138773E+04 0.12479540E+04 0.20000000E+03 0. 0.
     -13 1 0 0 0 0 -0.60752276E+02 0.39187066E+02 -0.10603642E+02 0.73067870E+02 0.10566000E+00 0. 0.
      13 1 0 0 0 0 0.19078494E+02 0.78894190E+02 0.13655311E+02 0.82308935E+02 0.10566000E+00 0. 0.
      21 1 0 0 0 0 0.19881752E+03 -0.12887481E+03 -0.37975928E+03 0.44763579E+03 0.48546086E+01 0. 0.
      21 1 0 0 0 0 -0.21349587E+03 0.11188772E+02 -0.12295054E+04 0.12479540E+04 0.00000000E+00 0. 0.
      21 1 0 0 0 0 0.84554538E+02 0.18797700E+01 -0.22143626E+02 0.88237686E+02 0.11939250E+02 0. 0.
</event>

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

Hi,

You need to modify the file pythia-pgs/src/hep2lhe.f
and around line 291 adding a condition for stating that 1009 should be put in status 1
the status is saved in ISTUP variable.

This being said, I would strongly suggest to use the official pythia output (the .stdhep) instead of using the converted/simplified output that we create to make plot.

Cheers,

Olivier

On 23 Apr 2015, at 13:21, Renjie Wang <email address hidden> wrote:

> New question #265753 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/265753
>
> Dear experts,
>
> I have a quick question about particle status in LHE file after run with MG+pythia.
>
> The question is I have a stable Dirac fermion particle (1009,-1009), when just run MG, the everything is good, but after parton shower with pythia, the particle goes to the intermediator status (2), not stable status (1).
>
> How can I solve it?
>
> Thanks ahead;
> Ren-Jie
>
> <event>
> 13 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
> -2 -1 0 0 0 0 -0.85457660E+02 -0.12572625E+01 0.13608267E+03 0.16069563E+03 0.00000000E+00 0. 0.
> 2 -1 0 0 0 0 0.66168770E+00 0.15730045E+01 -0.16887368E+04 0.16887377E+04 0.00000000E+00 0. 0.
> 23 2 1 1 0 0 -0.41673835E+02 0.11808129E+03 0.30516587E+01 0.15537687E+03 0.91938886E+02 0. 0.
> 1008 2 1 1 0 0 -0.43113524E+02 -0.11776541E+03 -0.15557334E+04 0.16940541E+04 0.65862396E+03 0. 0.
> -13 2 3 3 0 0 -0.60752329E+02 0.39187104E+02 -0.10603652E+02 0.73067936E+02 0.10566000E+00 0. 0.
> 13 2 3 3 0 0 0.19078494E+02 0.78894190E+02 0.13655311E+02 0.82308935E+02 0.10566000E+00 0. 0.
> 1009 2 4 4 0 0 0.16637392E+03 -0.12024735E+03 -0.34185607E+03 0.44610008E+03 0.20000000E+03 0. 0.
> -1009 2 4 4 0 0 -0.20948745E+03 0.24819364E+01 -0.12138773E+04 0.12479540E+04 0.20000000E+03 0. 0.
> -13 1 0 0 0 0 -0.60752276E+02 0.39187066E+02 -0.10603642E+02 0.73067870E+02 0.10566000E+00 0. 0.
> 13 1 0 0 0 0 0.19078494E+02 0.78894190E+02 0.13655311E+02 0.82308935E+02 0.10566000E+00 0. 0.
> 21 1 0 0 0 0 0.19881752E+03 -0.12887481E+03 -0.37975928E+03 0.44763579E+03 0.48546086E+01 0. 0.
> 21 1 0 0 0 0 -0.21349587E+03 0.11188772E+02 -0.12295054E+04 0.12479540E+04 0.00000000E+00 0. 0.
> 21 1 0 0 0 0 0.84554538E+02 0.18797700E+01 -0.22143626E+02 0.88237686E+02 0.11939250E+02 0. 0.
> </event>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Renjie Wang (r-jiewang) said : 		#2

Thanks Olivier Mattelaer, that solved my question.

Revision history for this message
Renjie Wang (r-jiewang) said : 		#3

Hi Olivier,

Great, this is exactly what I need, thanks for the help.

Ren-Jie

> On Apr 23, 2015, at 2:36 PM, Olivier Mattelaer <email address hidden> wrote:
>
> Your question #265753 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/265753
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
>
> You need to modify the file pythia-pgs/src/hep2lhe.f
> and around line 291 adding a condition for stating that 1009 should be put in status 1
> the status is saved in ISTUP variable.
>
> This being said, I would strongly suggest to use the official pythia
> output (the .stdhep) instead of using the converted/simplified output
> that we create to make plot.
>
> Cheers,
>
> Olivier
>
>
> On 23 Apr 2015, at 13:21, Renjie Wang <email address hidden> wrote:
>
>> New question #265753 on MadGraph5_aMC@NLO:
>> https://answers.launchpad.net/mg5amcnlo/+question/265753
>>
>> Dear experts,
>>
>> I have a quick question about particle status in LHE file after run with MG+pythia.
>>
>> The question is I have a stable Dirac fermion particle (1009,-1009), when just run MG, the everything is good, but after parton shower with pythia, the particle goes to the intermediator status (2), not stable status (1).
>>
>> How can I solve it?
>>
>> Thanks ahead;
>> Ren-Jie
>>
>> <event>
>> 13 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
>> -2 -1 0 0 0 0 -0.85457660E+02 -0.12572625E+01 0.13608267E+03 0.16069563E+03 0.00000000E+00 0. 0.
>> 2 -1 0 0 0 0 0.66168770E+00 0.15730045E+01 -0.16887368E+04 0.16887377E+04 0.00000000E+00 0. 0.
>> 23 2 1 1 0 0 -0.41673835E+02 0.11808129E+03 0.30516587E+01 0.15537687E+03 0.91938886E+02 0. 0.
>> 1008 2 1 1 0 0 -0.43113524E+02 -0.11776541E+03 -0.15557334E+04 0.16940541E+04 0.65862396E+03 0. 0.
>> -13 2 3 3 0 0 -0.60752329E+02 0.39187104E+02 -0.10603652E+02 0.73067936E+02 0.10566000E+00 0. 0.
>> 13 2 3 3 0 0 0.19078494E+02 0.78894190E+02 0.13655311E+02 0.82308935E+02 0.10566000E+00 0. 0.
>> 1009 2 4 4 0 0 0.16637392E+03 -0.12024735E+03 -0.34185607E+03 0.44610008E+03 0.20000000E+03 0. 0.
>> -1009 2 4 4 0 0 -0.20948745E+03 0.24819364E+01 -0.12138773E+04 0.12479540E+04 0.20000000E+03 0. 0.
>> -13 1 0 0 0 0 -0.60752276E+02 0.39187066E+02 -0.10603642E+02 0.73067870E+02 0.10566000E+00 0. 0.
>> 13 1 0 0 0 0 0.19078494E+02 0.78894190E+02 0.13655311E+02 0.82308935E+02 0.10566000E+00 0. 0.
>> 21 1 0 0 0 0 0.19881752E+03 -0.12887481E+03 -0.37975928E+03 0.44763579E+03 0.48546086E+01 0. 0.
>> 21 1 0 0 0 0 -0.21349587E+03 0.11188772E+02 -0.12295054E+04 0.12479540E+04 0.00000000E+00 0. 0.
>> 21 1 0 0 0 0 0.84554538E+02 0.18797700E+01 -0.22143626E+02 0.88237686E+02 0.11939250E+02 0. 0.
>> </event>
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/mg5amcnlo/+question/265753/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/265753
>
> You received this question notification because you asked the question.