command in shower_card.dat to keep top stable

Dear Shankha,

the shower will complain if it finds an undecayed tops, as it will have to built top-flavoured
hadrons. I think no sensible thing can emerge from that, and probably showers in the Herwig family will directly forbid (i.e. simply ignore) this setting.

This said, the commands to modify the decay behaviour can be found for example in MCatNLO/Scripts/MCatNLO_MadFKS_**.Script, where ** refers to the shower you're using.

If you look for example for the string 'TAUP_STABLE' inside those files, you will find the
(shower specific) commands to use to set the tau+ stable, for example

set /Herwig/Particles/tau+:Stable Stable for Herwig++, or
-15:maydecay = false ! stable tau+ for Pythia8,

and you can easily adapt those commands to the top case.
I reiterate once more not to use this for physics purposes, but just for checking, as you said.

Cheers.
Paolo

Dear Paolo,

Thank you very much for the detailed explanation. I will just use this for a check and will not use this for physics purposes.

Cheers,

Shankha.

Dear Paolo,

I took this file :

~/MG5_aMC_v2_3_3/Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIG6.Script

added the lines :

  $TP_STABLE_HW ! .TRUE. if one wants t stable
  $TM_STABLE_HW ! .TRUE. if one wants t~ stable

analogous to the tau lines

and

~/MG5_aMC_v2_3_3/Template/NLO/MCatNLO/Scripts/MCatNLO_MadFKS_HERWIGPP.Script

added the lines :

if [ $TP_STABLE_HWPP = ".TRUE." ]
then
cat <<EOF >> ./$ifile

#
# Prevent particles from decaying
# Particle names can be found in Herwig++/src/default/mesons.in or similar
#

# tau+
set /Herwig/Particles/t:Stable Stable

EOF
fi

if [ $TM_STABLE_HWPP = ".TRUE." ]
then
cat <<EOF >> ./$ifile

#
# Prevent particles from decaying
# Particle names can be found in Herwig++/src/default/mesons.in or similar
#

# tau-
set /Herwig/Particles/tbar:Stable Stable

EOF
fi

Then in the shower_card.dat file, I am adding these lines :

tp_stable = T # set t's stable
tm_stable = T # set tbar's stable

However, still I am getting this error :

Command "launch " interrupted with error:
ShowerCardError : Unknown entry: tp_stable = F
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/home/shankha/Packages_new/new_mg2/MG5_aMC_v2_3_3/tth/pp2tth_NLO2/run_04_tag_1_debug.log'.
Please attach this file to your report.

Can you please let me know what I am missing?

Thank you very much.

Cheers,

Shankha.

Dear Shankha,

sorry, I meant something simpler.

For Herwig6, just use the default MCatNLO_MadFKS_HERWIG6.Script. Go in srcHerwig/madfks_hwdriver.f, around line 223, and add

CALL HWUSTA('TQRK ')
CALL HWUSTA('TBAR ')

For Herwig++, just start from the default MCatNLO_MadFKS_HERWIGPP.Script
and simply add the lines

set /Herwig/Particles/t:Stable Stable
set /Herwig/Particles/tbar:Stable Stable

around line 785.
Namely, do not create if clauses if[ TP_STABLE_HWPP=’.TRUE.’ ], since
this would require modifying other parts of the code which is what
is creating the errors you see.

However, I doubt the showers will let you shower the sample with stable
top: most probably these settings will simply be ignored, or make the
run crash.

Cheers.
Paolo

Dear Paolo,

Thank you so much for your detailed message. I changed the files according to what you wrote but still I find that the tops are decaying. In my events_HERWIG6_0.lhe, I can see that the t and tbar have status code 2 and not 1. I am generating a t\bar{t}h event sample and I don't even see the Higgs as stable. I have put the line :

h_stable = T # set Higgs' stable

in the shower_card.dat file. An event in the shower lhe file looks like this :

<event>
 16 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
       6 2 44 44 0 0 0.86469805E+01 0.93705046E+02 0.49686618E+03 0.53447209E+03 0.17300000E+03 0. 0.
      -6 2 46 46 0 0 -0.77668565E+01 -0.23215025E+03 0.11714786E+03 0.31242069E+03 0.17300000E+03 0. 0.
       5 1 0 0 0 0 0.19494890E+03 -0.85966004E+03 0.14858672E+03 0.91343073E+03 0.18776908E+03 0. 0.
      21 1 0 0 0 0 -0.10464783E+03 0.61415785E+03 0.10047545E+04 0.12029151E+04 0.22210807E+03 0. 0.
       5 1 0 0 0 0 0.56191868E+03 0.16025025E+03 0.11747893E+04 0.13130479E+04 0.50319483E+02 0. 0.
       5 1 0 0 0 0 -0.29010313E+03 -0.34552728E+03 0.62685072E+03 0.77637597E+03 0.79176335E+02 0. 0.
      21 1 0 0 0 0 -0.13545610E+03 0.25184460E+02 0.10490260E+04 0.10586092E+04 0.34862052E+02 0. 0.
      21 1 0 0 0 0 -0.53676168E+02 0.11379287E+03 0.56907855E+03 0.58404589E+03 0.37805596E+02 0. 0.
      21 1 0 0 0 0 0.36717296E+02 0.23596684E+02 0.36091116E+03 0.36367269E+03 0.97977541E+01 0. 0.
       5 1 0 0 0 0 0.79795673E+01 0.25450546E+02 0.27793107E+02 0.38777975E+02 0.44576625E+01 0. 0.
      21 1 0 0 0 0 0.11732058E+02 0.21531422E+02 0.36894152E+01 0.25598315E+02 0.63575618E+01 0. 0.
      21 1 0 0 0 0 -0.81197647E-01 -0.17381681E+02 0.34545423E+02 0.38771592E+02 0.27786049E+01 0. 0.
      21 1 0 0 0 0 -0.14394840E+02 -0.78274017E+01 -0.51734072E+02 0.54671225E+02 0.66369423E+01 0. 0.
      21 1 0 0 0 0 -0.10050047E+02 -0.93661937E+01 0.17431262E+03 0.17489650E+03 0.38943620E+01 0. 0.
      21 1 0 0 0 0 -0.86842527E+01 0.81693277E+01 -0.23228096E+03 0.23261598E+03 0.36870730E+01 0. 0.
      12 1 0 0 0 0 0.40611439E+01 0.77583442E+01 0.00000000E+00 0.87569855E+01 0.00000000E+00 0. 0.
</event>

Am I still missing something?

I apologise for these naive questions.

Thank you very much.

Best regards,

Shankha.

Dear Shankha,

as I said, Herwig6 does not allow the top (or the Higgs) to stay undecayed.
I thought it would crash, but instead, most probably, it has just output a
warning message saying that it could not avoid decaying them.

You should find some message like that in file MCatNLO/RUN_HERWIG6**/
mcatnlo_run.log, where ** is an integer that labels your run.

I have actually a question on the sample lhe file you sent: what file is this?
Is it directly created by your analysis after shower, or by the stdhep interface,
or something else?

Cheers.
Paolo

On 03 Nov 2016, at 11:08, Shankha Banerjee <email address hidden> wrote:

> Question #403543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/403543
>
> Status: Answered => Open
>
> Shankha Banerjee is still having a problem:
> Dear Paolo,
>
> Thank you so much for your detailed message. I changed the files
> according to what you wrote but still I find that the tops are decaying.
> In my events_HERWIG6_0.lhe, I can see that the t and tbar have status
> code 2 and not 1. I am generating a t\bar{t}h event sample and I don't
> even see the Higgs as stable. I have put the line :
>
> h_stable = T # set Higgs' stable
>
> in the shower_card.dat file. An event in the shower lhe file looks like
> this :
>
> <event>
> 16 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
> 6 2 44 44 0 0 0.86469805E+01 0.93705046E+02 0.49686618E+03 0.53447209E+03 0.17300000E+03 0. 0.
> -6 2 46 46 0 0 -0.77668565E+01 -0.23215025E+03 0.11714786E+03 0.31242069E+03 0.17300000E+03 0. 0.
> 5 1 0 0 0 0 0.19494890E+03 -0.85966004E+03 0.14858672E+03 0.91343073E+03 0.18776908E+03 0. 0.
> 21 1 0 0 0 0 -0.10464783E+03 0.61415785E+03 0.10047545E+04 0.12029151E+04 0.22210807E+03 0. 0.
> 5 1 0 0 0 0 0.56191868E+03 0.16025025E+03 0.11747893E+04 0.13130479E+04 0.50319483E+02 0. 0.
> 5 1 0 0 0 0 -0.29010313E+03 -0.34552728E+03 0.62685072E+03 0.77637597E+03 0.79176335E+02 0. 0.
> 21 1 0 0 0 0 -0.13545610E+03 0.25184460E+02 0.10490260E+04 0.10586092E+04 0.34862052E+02 0. 0.
> 21 1 0 0 0 0 -0.53676168E+02 0.11379287E+03 0.56907855E+03 0.58404589E+03 0.37805596E+02 0. 0.
> 21 1 0 0 0 0 0.36717296E+02 0.23596684E+02 0.36091116E+03 0.36367269E+03 0.97977541E+01 0. 0.
> 5 1 0 0 0 0 0.79795673E+01 0.25450546E+02 0.27793107E+02 0.38777975E+02 0.44576625E+01 0. 0.
> 21 1 0 0 0 0 0.11732058E+02 0.21531422E+02 0.36894152E+01 0.25598315E+02 0.63575618E+01 0. 0.
> 21 1 0 0 0 0 -0.81197647E-01 -0.17381681E+02 0.34545423E+02 0.38771592E+02 0.27786049E+01 0. 0.
> 21 1 0 0 0 0 -0.14394840E+02 -0.78274017E+01 -0.51734072E+02 0.54671225E+02 0.66369423E+01 0. 0.
> 21 1 0 0 0 0 -0.10050047E+02 -0.93661937E+01 0.17431262E+03 0.17489650E+03 0.38943620E+01 0. 0.
> 21 1 0 0 0 0 -0.86842527E+01 0.81693277E+01 -0.23228096E+03 0.23261598E+03 0.36870730E+01 0. 0.
> 12 1 0 0 0 0 0.40611439E+01 0.77583442E+01 0.00000000E+00 0.87569855E+01 0.00000000E+00 0. 0.
> </event>
>
> Am I still missing something?
>
> I apologise for these naive questions.
>
> Thank you very much.
>
> Best regards,
>
> Shankha.
>
> --
> You received this question notification because you are subscribed to
> the question.

Hi,

> I have actually a question on the sample lhe file you sent: what file is this?
> Is it directly created by your analysis after shower, or by the stdhep interface,
> or something else?

Looks like the lhe file created by MadAnalysis to make default plot.
That file/plot should be used with a lot of caution since :
1) this does not contain the full information
2) the clustering algorithm used by MadAnalysis is a simple cone algorithm (so not infrared safe)
3) this type of file can not be used for detector simulation.
But for a quick look for checking sanity of the generation this is ok (and actually convenient).
I would never use such plot/output for publication.

Cheers,

Olivier

> On Nov 3, 2016, at 12:09, Paolo Torrielli <email address hidden> wrote:
>
> Question #403543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/403543
>
> Status: Open => Answered
>
> Paolo Torrielli proposed the following answer:
> Dear Shankha,
>
> as I said, Herwig6 does not allow the top (or the Higgs) to stay undecayed.
> I thought it would crash, but instead, most probably, it has just output a
> warning message saying that it could not avoid decaying them.
>
> You should find some message like that in file MCatNLO/RUN_HERWIG6**/
> mcatnlo_run.log, where ** is an integer that labels your run.
>
> I have actually a question on the sample lhe file you sent: what file is this?
> Is it directly created by your analysis after shower, or by the stdhep interface,
> or something else?
>
> Cheers.
> Paolo
>
> On 03 Nov 2016, at 11:08, Shankha Banerjee
> <email address hidden> wrote:
>
>> Question #403543 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/403543
>>
>> Status: Answered => Open
>>
>> Shankha Banerjee is still having a problem:
>> Dear Paolo,
>>
>> Thank you so much for your detailed message. I changed the files
>> according to what you wrote but still I find that the tops are decaying.
>> In my events_HERWIG6_0.lhe, I can see that the t and tbar have status
>> code 2 and not 1. I am generating a t\bar{t}h event sample and I don't
>> even see the Higgs as stable. I have put the line :
>>
>> h_stable = T # set Higgs' stable
>>
>> in the shower_card.dat file. An event in the shower lhe file looks like
>> this :
>>
>> <event>
>> 16 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
>> 6 2 44 44 0 0 0.86469805E+01 0.93705046E+02 0.49686618E+03 0.53447209E+03 0.17300000E+03 0. 0.
>> -6 2 46 46 0 0 -0.77668565E+01 -0.23215025E+03 0.11714786E+03 0.31242069E+03 0.17300000E+03 0. 0.
>> 5 1 0 0 0 0 0.19494890E+03 -0.85966004E+03 0.14858672E+03 0.91343073E+03 0.18776908E+03 0. 0.
>> 21 1 0 0 0 0 -0.10464783E+03 0.61415785E+03 0.10047545E+04 0.12029151E+04 0.22210807E+03 0. 0.
>> 5 1 0 0 0 0 0.56191868E+03 0.16025025E+03 0.11747893E+04 0.13130479E+04 0.50319483E+02 0. 0.
>> 5 1 0 0 0 0 -0.29010313E+03 -0.34552728E+03 0.62685072E+03 0.77637597E+03 0.79176335E+02 0. 0.
>> 21 1 0 0 0 0 -0.13545610E+03 0.25184460E+02 0.10490260E+04 0.10586092E+04 0.34862052E+02 0. 0.
>> 21 1 0 0 0 0 -0.53676168E+02 0.11379287E+03 0.56907855E+03 0.58404589E+03 0.37805596E+02 0. 0.
>> 21 1 0 0 0 0 0.36717296E+02 0.23596684E+02 0.36091116E+03 0.36367269E+03 0.97977541E+01 0. 0.
>> 5 1 0 0 0 0 0.79795673E+01 0.25450546E+02 0.27793107E+02 0.38777975E+02 0.44576625E+01 0. 0.
>> 21 1 0 0 0 0 0.11732058E+02 0.21531422E+02 0.36894152E+01 0.25598315E+02 0.63575618E+01 0. 0.
>> 21 1 0 0 0 0 -0.81197647E-01 -0.17381681E+02 0.34545423E+02 0.38771592E+02 0.27786049E+01 0. 0.
>> 21 1 0 0 0 0 -0.14394840E+02 -0.78274017E+01 -0.51734072E+02 0.54671225E+02 0.66369423E+01 0. 0.
>> 21 1 0 0 0 0 -0.10050047E+02 -0.93661937E+01 0.17431262E+03 0.17489650E+03 0.38943620E+01 0. 0.
>> 21 1 0 0 0 0 -0.86842527E+01 0.81693277E+01 -0.23228096E+03 0.23261598E+03 0.36870730E+01 0. 0.
>> 12 1 0 0 0 0 0.40611439E+01 0.77583442E+01 0.00000000E+00 0.87569855E+01 0.00000000E+00 0. 0.
>> </event>
>>
>> Am I still missing something?
>>
>> I apologise for these naive questions.
>>
>> Thank you very much.
>>
>> Best regards,
>>
>> Shankha.
>>
>> --
>> You received this question notification because you are subscribed to
>> the question.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Dear Paolo and Olivier,

Thank you very much for your replies. So as far as I gather, it is then not possible to have a situation with the higgs and the tops undecayed. I was not going to use this for any physics purposes but wanted to check some of the jets not coming from the top decays.

I am attaching my mcatnlo_run.log file. I do not see any warnings pertaining to the decays.

I took an event from the events_HERWIG6_0.lhe file which of course can not be used for detector simulation as Olivier said.

So, what I gather from these is that it is not possible to keep the tops and Higgs stable and have the showering at the same time.

Please correct me if I am still missing something.

Thank you very much for your time.

----- Original Message -----
From: "Olivier Mattelaer" <email address hidden>
To: <email address hidden>
Sent: Thursday, 3 November, 2016 17:23:22
Subject: Re: [Question #403543]: command in shower_card.dat to keep top stable

Your question #403543 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/403543

Olivier Mattelaer proposed the following answer:
Hi,

> I have actually a question on the sample lhe file you sent: what file is this?
> Is it directly created by your analysis after shower, or by the stdhep interface,
> or something else?

Looks like the lhe file created by MadAnalysis to make default plot.
That file/plot should be used with a lot of caution since :
1) this does not contain the full information
2) the clustering algorithm used by MadAnalysis is a simple cone algorithm (so not infrared safe)
3) this type of file can not be used for detector simulation.
But for a quick look for checking sanity of the generation this is ok (and actually convenient).
I would never use such plot/output for publication.

Cheers,

Olivier

> On Nov 3, 2016, at 12:09, Paolo Torrielli <email address hidden> wrote:
>
> Question #403543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/403543
>
> Status: Open => Answered
>
> Paolo Torrielli proposed the following answer:
> Dear Shankha,
>
> as I said, Herwig6 does not allow the top (or the Higgs) to stay undecayed.
> I thought it would crash, but instead, most probably, it has just output a
> warning message saying that it could not avoid decaying them.
>
> You should find some message like that in file MCatNLO/RUN_HERWIG6**/
> mcatnlo_run.log, where ** is an integer that labels your run.
>
> I have actually a question on the sample lhe file you sent: what file is this?
> Is it directly created by your analysis after shower, or by the stdhep interface,
> or something else?
>
> Cheers.
> Paolo
>
> On 03 Nov 2016, at 11:08, Shankha Banerjee
> <email address hidden> wrote:
>
>> Question #403543 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/403543
>>
>> Status: Answered => Open
>>
>> Shankha Banerjee is still having a problem:
>> Dear Paolo,
>>
>> Thank you so much for your detailed message. I changed the files
>> according to what you wrote but still I find that the tops are decaying.
>> In my events_HERWIG6_0.lhe, I can see that the t and tbar have status
>> code 2 and not 1. I am generating a t\bar{t}h event sample and I don't
>> even see the Higgs as stable. I have put the line :
>>
>> h_stable = T # set Higgs' stable
>>
>> in the shower_card.dat file. An event in the shower lhe file looks like
>> this :
>>
>> <event>
>> 16 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
>> 6 2 44 44 0 0 0.86469805E+01 0.93705046E+02 0.49686618E+03 0.53447209E+03 0.17300000E+03 0. 0.
>> -6 2 46 46 0 0 -0.77668565E+01 -0.23215025E+03 0.11714786E+03 0.31242069E+03 0.17300000E+03 0. 0.
>> 5 1 0 0 0 0 0.19494890E+03 -0.85966004E+03 0.14858672E+03 0.91343073E+03 0.18776908E+03 0. 0.
>> 21 1 0 0 0 0 -0.10464783E+03 0.61415785E+03 0.10047545E+04 0.12029151E+04 0.22210807E+03 0. 0.
>> 5 1 0 0 0 0 0.56191868E+03 0.16025025E+03 0.11747893E+04 0.13130479E+04 0.50319483E+02 0. 0.
>> 5 1 0 0 0 0 -0.29010313E+03 -0.34552728E+03 0.62685072E+03 0.77637597E+03 0.79176335E+02 0. 0.
>> 21 1 0 0 0 0 -0.13545610E+03 0.25184460E+02 0.10490260E+04 0.10586092E+04 0.34862052E+02 0. 0.
>> 21 1 0 0 0 0 -0.53676168E+02 0.11379287E+03 0.56907855E+03 0.58404589E+03 0.37805596E+02 0. 0.
>> 21 1 0 0 0 0 0.36717296E+02 0.23596684E+02 0.36091116E+03 0.36367269E+03 0.97977541E+01 0. 0.
>> 5 1 0 0 0 0 0.79795673E+01 0.25450546E+02 0.27793107E+02 0.38777975E+02 0.44576625E+01 0. 0.
>> 21 1 0 0 0 0 0.11732058E+02 0.21531422E+02 0.36894152E+01 0.25598315E+02 0.63575618E+01 0. 0.
>> 21 1 0 0 0 0 -0.81197647E-01 -0.17381681E+02 0.34545423E+02 0.38771592E+02 0.27786049E+01 0. 0.
>> 21 1 0 0 0 0 -0.14394840E+02 -0.78274017E+01 -0.51734072E+02 0.54671225E+02 0.66369423E+01 0. 0.
>> 21 1 0 0 0 0 -0.10050047E+02 -0.93661937E+01 0.17431262E+03 0.17489650E+03 0.38943620E+01 0. 0.
>> 21 1 0 0 0 0 -0.86842527E+01 0.81693277E+01 -0.23228096E+03 0.23261598E+03 0.36870730E+01 0. 0.
>> 12 1 0 0 0 0 0.40611439E+01 0.77583442E+01 0.00000000E+00 0.87569855E+01 0.00000000E+00 0. 0.
>> </event>
>>
>> Am I still missing something?
>>
>> I apologise for these naive questions.
>>
>> Thank you very much.
>>
>> Best regards,
>>
>> Shankha.
>>
>> --
>> You received this question notification because you are subscribed to
>> the question.
>
> --
> 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/403543/+confirm?answer_id=6

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/403543

You received this question notification because you asked the question.

Dear Shankha,

I cannot see the attachment, but however I trust you that no warning
is issued. Indeed, I run a quick p p > H sample requiring the Higgs
to be stable, and this setting was simply ignored by Herwig6, even
though the event record displays the message

          PARTICLE TYPE 201=HIGGS SET STABLE

So yes, with Herwig6 you cannot keep tops, H’s, W’s and Z’s stable,
as far as I can tell (but maybe the authors of Herwig6 would know how
to bypass this). In Pythia you would be allowed to keep H’s, W’s and
Z’s stable, and also tops, even though in the latter case it would
complain about the formation of top-flavoured hadrons.

Cheers.
Paolo

On 03 Nov 2016, at 13:27, Shankha Banerjee <email address hidden> wrote:

> Question #403543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/403543
>
> Status: Answered => Open
>
> Shankha Banerjee is still having a problem:
> Dear Paolo and Olivier,
>
> Thank you very much for your replies. So as far as I gather, it is then
> not possible to have a situation with the higgs and the tops undecayed.
> I was not going to use this for any physics purposes but wanted to check
> some of the jets not coming from the top decays.
>
> I am attaching my mcatnlo_run.log file. I do not see any warnings
> pertaining to the decays.
>
> I took an event from the events_HERWIG6_0.lhe file which of course can
> not be used for detector simulation as Olivier said.
>
> So, what I gather from these is that it is not possible to keep the tops
> and Higgs stable and have the showering at the same time.
>
> Please correct me if I am still missing something.
>
> Thank you very much for your time.
>
> ----- Original Message -----
> From: "Olivier Mattelaer" <email address hidden>
> To: <email address hidden>
> Sent: Thursday, 3 November, 2016 17:23:22
> Subject: Re: [Question #403543]: command in shower_card.dat to keep top stable
>
> Your question #403543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/403543
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
>> I have actually a question on the sample lhe file you sent: what file is this?
>> Is it directly created by your analysis after shower, or by the stdhep interface,
>> or something else?
>
> Looks like the lhe file created by MadAnalysis to make default plot.
> That file/plot should be used with a lot of caution since :
> 1) this does not contain the full information
> 2) the clustering algorithm used by MadAnalysis is a simple cone algorithm (so not infrared safe)
> 3) this type of file can not be used for detector simulation.
> But for a quick look for checking sanity of the generation this is ok (and actually convenient).
> I would never use such plot/output for publication.
>
> Cheers,
>
> Olivier
>
>> On Nov 3, 2016, at 12:09, Paolo Torrielli <email address hidden> wrote:
>>
>> Question #403543 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/403543
>>
>> Status: Open => Answered
>>
>> Paolo Torrielli proposed the following answer:
>> Dear Shankha,
>>
>> as I said, Herwig6 does not allow the top (or the Higgs) to stay undecayed.
>> I thought it would crash, but instead, most probably, it has just output a
>> warning message saying that it could not avoid decaying them.
>>
>> You should find some message like that in file MCatNLO/RUN_HERWIG6**/
>> mcatnlo_run.log, where ** is an integer that labels your run.
>>
>> I have actually a question on the sample lhe file you sent: what file is this?
>> Is it directly created by your analysis after shower, or by the stdhep interface,
>> or something else?
>>
>> Cheers.
>> Paolo
>>
>> On 03 Nov 2016, at 11:08, Shankha Banerjee
>> <email address hidden> wrote:
>>
>>> Question #403543 on MadGraph5_aMC@NLO changed:
>>> https://answers.launchpad.net/mg5amcnlo/+question/403543
>>>
>>> Status: Answered => Open
>>>
>>> Shankha Banerjee is still having a problem:
>>> Dear Paolo,
>>>
>>> Thank you so much for your detailed message. I changed the files
>>> according to what you wrote but still I find that the tops are decaying.
>>> In my events_HERWIG6_0.lhe, I can see that the t and tbar have status
>>> code 2 and not 1. I am generating a t\bar{t}h event sample and I don't
>>> even see the Higgs as stable. I have put the line :
>>>
>>> h_stable = T # set Higgs' stable
>>>
>>> in the shower_card.dat file. An event in the shower lhe file looks like
>>> this :
>>>
>>> <event>
>>> 16 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
>>> 6 2 44 44 0 0 0.86469805E+01 0.93705046E+02 0.49686618E+03 0.53447209E+03 0.17300000E+03 0. 0.
>>> -6 2 46 46 0 0 -0.77668565E+01 -0.23215025E+03 0.11714786E+03 0.31242069E+03 0.17300000E+03 0. 0.
>>> 5 1 0 0 0 0 0.19494890E+03 -0.85966004E+03 0.14858672E+03 0.91343073E+03 0.18776908E+03 0. 0.
>>> 21 1 0 0 0 0 -0.10464783E+03 0.61415785E+03 0.10047545E+04 0.12029151E+04 0.22210807E+03 0. 0.
>>> 5 1 0 0 0 0 0.56191868E+03 0.16025025E+03 0.11747893E+04 0.13130479E+04 0.50319483E+02 0. 0.
>>> 5 1 0 0 0 0 -0.29010313E+03 -0.34552728E+03 0.62685072E+03 0.77637597E+03 0.79176335E+02 0. 0.
>>> 21 1 0 0 0 0 -0.13545610E+03 0.25184460E+02 0.10490260E+04 0.10586092E+04 0.34862052E+02 0. 0.
>>> 21 1 0 0 0 0 -0.53676168E+02 0.11379287E+03 0.56907855E+03 0.58404589E+03 0.37805596E+02 0. 0.
>>> 21 1 0 0 0 0 0.36717296E+02 0.23596684E+02 0.36091116E+03 0.36367269E+03 0.97977541E+01 0. 0.
>>> 5 1 0 0 0 0 0.79795673E+01 0.25450546E+02 0.27793107E+02 0.38777975E+02 0.44576625E+01 0. 0.
>>> 21 1 0 0 0 0 0.11732058E+02 0.21531422E+02 0.36894152E+01 0.25598315E+02 0.63575618E+01 0. 0.
>>> 21 1 0 0 0 0 -0.81197647E-01 -0.17381681E+02 0.34545423E+02 0.38771592E+02 0.27786049E+01 0. 0.
>>> 21 1 0 0 0 0 -0.14394840E+02 -0.78274017E+01 -0.51734072E+02 0.54671225E+02 0.66369423E+01 0. 0.
>>> 21 1 0 0 0 0 -0.10050047E+02 -0.93661937E+01 0.17431262E+03 0.17489650E+03 0.38943620E+01 0. 0.
>>> 21 1 0 0 0 0 -0.86842527E+01 0.81693277E+01 -0.23228096E+03 0.23261598E+03 0.36870730E+01 0. 0.
>>> 12 1 0 0 0 0 0.40611439E+01 0.77583442E+01 0.00000000E+00 0.87569855E+01 0.00000000E+00 0. 0.
>>> </event>
>>>
>>> Am I still missing something?
>>>
>>> I apologise for these naive questions.
>>>
>>> Thank you very much.
>>>
>>> Best regards,
>>>
>>> Shankha.
>>>
>>> --
>>> You received this question notification because you are subscribed to
>>> the question.
>>
>> --
>> 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/403543/+confirm?answer_id=6
>
> 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/403543
>
> You received this question notification because you asked the question.
>
> You received this question notification because you are subscribed to
> the question.