set DECAY in Madgraph2.5.2

Hi,

First the set command only modify the total width and do not touch to the Branching ratio.
When a width is set on Auto then the code will replace the "Auto" value and the various branching ratio associated.

So this explains that your SVN width still has BR information since we do not touch those.
The same for your gluino case, since it looks like that the width is still on Auto (i.e. the computation of the width is not done yet and the BR information is still the old one)

If you want to force the computation of the auto-width such that you can check the result before generating the events, you can type
"compute_widths" after have set all the width that you want to compute on Auto.
You can also do "compute_widths 1000021 1000016"
Such that you do not have to set the width to auto in the first place.

Cheers,

Olivier

> On 8 Jul 2017, at 16:44, yingpu xiahou <email address hidden> wrote:
>
> New question #647138 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/647138
>
> Hi,
>
> Sorry for disturbing.
>
> I am using MG5_2_5_2 and MSSM model. The script are:
>
> import model sm
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> import model mssm
> generate p p > go go, (go > n1 t t~)
> add process p p > go go j, (go > n1 t t~)
> add process p p > go go j j, (go > n1 t t~)
> output xxx
> launch
> set DECAY 1000016 0 #wsn
> set DECAY 1000021 Auto #wgo
>
> I only want 1000021 to happen 3-body decay (go > n1 t t~).
>
> But, when I check the param_card.dat, I found all 1000021 decay are 2-body:
>
> DECAY 1000021 Auto # gluino decays
> 2.084542e-02 2 -1 1000001 # br(~g -> ~d_l db)
> 2.084542e-02 2 -1000001 1 # br(~g -> ~d_l* d )
> 5.070753e-02 2 -1 2000001 # br(~g -> ~d_r db)
> 5.070753e-02 2 -2000001 1 # br(~g -> ~d_r* d )
> 2.897878e-02 2 -2 1000002 # br(~g -> ~u_l ub)
> 2.897878e-02 2 -1000002 2 # br(~g -> ~u_l* u )
> 4.468728e-02 2 -2 2000002 # br(~g -> ~u_r ub)
> 4.468728e-02 2 -2000002 2 # br(~g -> ~u_r* u )
> 2.084542e-02 2 -3 1000003 # br(~g -> ~s_l sb)
> 2.084542e-02 2 -1000003 3 # br(~g -> ~s_l* s )
> 5.070753e-02 2 -3 2000003 # br(~g -> ~s_r sb)
> 5.070753e-02 2 -2000003 3 # br(~g -> ~s_r* s )
> 2.897878e-02 2 -4 1000004 # br(~g -> ~c_l cb)
> 2.897878e-02 2 -1000004 4 # br(~g -> ~c_l* c )
> 4.468728e-02 2 -4 2000004 # br(~g -> ~c_r cb)
> 4.468728e-02 2 -2000004 4 # br(~g -> ~c_r* c )
> 1.058402e-01 2 -5 1000005 # br(~g -> ~b_1 bb)
> 1.058402e-01 2 -1000005 5 # br(~g -> ~b_1* b )
> 5.565748e-02 2 -5 2000005 # br(~g -> ~b_2 bb)
> 5.565748e-02 2 -2000005 5 # br(~g -> ~b_2* b )
> 4.806428e-02 2 -6 1000006 # br(~g -> ~t_1 tb)
> 4.806428e-02 2 -1000006 6 # br(~g -> ~t_1* t )
> 0.000000e+00 2 -6 2000006 # br(~g -> ~t_2 tb)
> 0.000000e+00 2 -2000006 6 # br(~g -> ~t_2* t )
>
> And although I set 1000016 DECAY to 0, there is still some non-zero branching ratio appears:
>
> DECAY 1000016 0.000000e+00 # snu_taul decays
> 9.859945e-01 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
> 6.251296e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
> 0.000000e+00 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
> 0.000000e+00 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
> 7.754175e-03 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
> 0.000000e+00 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
> 0.000000e+00 2 -1000015 -37 # br(~nu_taul -> ~tau_1+ h-)
> 0.000000e+00 2 -2000015 -37 # br(~nu_taul -> ~tau_2+ h-)
> 0.000000e+00 2 -1000015 -24 # br(~nu_taul -> ~tau_1+ w-)
> 0.000000e+00 2 -2000015 -24 # br(~nu_taul -> ~tau_2+ w-)
> So, is there something wrong in my command?
>
> Thanks!
>
> Yingpu
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

Sorry. I am more interested in how to manually set Branching ratio currently. For PGS=1000016 as an example, how can I set all the bracing ratio of 1000016 to be zero?

DECAY 1000016 0.000000e+00 # snu_taul decays
    9.859945e-01 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
    6.251296e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
    0.000000e+00 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
    0.000000e+00 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
    7.754175e-03 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
    0.000000e+00 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
    0.000000e+00 2 -1000015 -37 # br(~nu_taul -> ~tau_1+ h-)
    0.000000e+00 2 -2000015 -37 # br(~nu_taul -> ~tau_2+ h-)
    0.000000e+00 2 -1000015 -24 # br(~nu_taul -> ~tau_1+ w-)
          0.000000e+00 2 -2000015 -24 # br(~nu_taul -> ~tau_2+ w-)

As I don’t want the Branching ratio to be non-zero.

Thanks,

Yingpu

在 2017年7月8日，19:57，Olivier Mattelaer <<email address hidden><mailto:<email address hidden>>> 写道：

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

   Status: Open => Answered

Olivier Mattelaer proposed the following answer:
Hi,

First the set command only modify the total width and do not touch to the Branching ratio.
When a width is set on Auto then the code will replace the "Auto" value and the various branching ratio associated.

So this explains that your SVN width still has BR information since we do not touch those.
The same for your gluino case, since it looks like that the width is still on Auto (i.e. the computation of the width is not done yet and the BR information is still the old one)

If you want to force the computation of the auto-width such that you can check the result before generating the events, you can type
"compute_widths" after have set all the width that you want to compute on Auto.
You can also do "compute_widths 1000021 1000016"
Such that you do not have to set the width to auto in the first place.

Cheers,

Olivier

On 8 Jul 2017, at 16:44, yingpu xiahou <<email address hidden><mailto:<email address hidden>>> wrote:

New question #647138 on MadGraph5_aMC@NLO:
https://answers.launchpad.net/mg5amcnlo/+question/647138

Hi,

Sorry for disturbing.

I am using MG5_2_5_2 and MSSM model. The script are:

import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model mssm
generate p p > go go, (go > n1 t t~)
add process p p > go go j, (go > n1 t t~)
add process p p > go go j j, (go > n1 t t~)
output xxx
launch
set DECAY 1000016 0 #wsn
set DECAY 1000021 Auto #wgo

I only want 1000021 to happen 3-body decay (go > n1 t t~).

But, when I check the param_card.dat, I found all 1000021 decay are 2-body:

DECAY 1000021 Auto # gluino decays
    2.084542e-02 2 -1 1000001 # br(~g -> ~d_l db)
    2.084542e-02 2 -1000001 1 # br(~g -> ~d_l* d )
    5.070753e-02 2 -1 2000001 # br(~g -> ~d_r db)
    5.070753e-02 2 -2000001 1 # br(~g -> ~d_r* d )
    2.897878e-02 2 -2 1000002 # br(~g -> ~u_l ub)
    2.897878e-02 2 -1000002 2 # br(~g -> ~u_l* u )
    4.468728e-02 2 -2 2000002 # br(~g -> ~u_r ub)
    4.468728e-02 2 -2000002 2 # br(~g -> ~u_r* u )
    2.084542e-02 2 -3 1000003 # br(~g -> ~s_l sb)
    2.084542e-02 2 -1000003 3 # br(~g -> ~s_l* s )
    5.070753e-02 2 -3 2000003 # br(~g -> ~s_r sb)
    5.070753e-02 2 -2000003 3 # br(~g -> ~s_r* s )
    2.897878e-02 2 -4 1000004 # br(~g -> ~c_l cb)
    2.897878e-02 2 -1000004 4 # br(~g -> ~c_l* c )
    4.468728e-02 2 -4 2000004 # br(~g -> ~c_r cb)
    4.468728e-02 2 -2000004 4 # br(~g -> ~c_r* c )
    1.058402e-01 2 -5 1000005 # br(~g -> ~b_1 bb)
    1.058402e-01 2 -1000005 5 # br(~g -> ~b_1* b )
    5.565748e-02 2 -5 2000005 # br(~g -> ~b_2 bb)
    5.565748e-02 2 -2000005 5 # br(~g -> ~b_2* b )
    4.806428e-02 2 -6 1000006 # br(~g -> ~t_1 tb)
    4.806428e-02 2 -1000006 6 # br(~g -> ~t_1* t )
    0.000000e+00 2 -6 2000006 # br(~g -> ~t_2 tb)
    0.000000e+00 2 -2000006 6 # br(~g -> ~t_2* t )

And although I set 1000016 DECAY to 0, there is still some non-zero branching ratio appears:

DECAY 1000016 0.000000e+00 # snu_taul decays
    9.859945e-01 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
    6.251296e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
    0.000000e+00 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
    0.000000e+00 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
    7.754175e-03 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
    0.000000e+00 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
    0.000000e+00 2 -1000015 -37 # br(~nu_taul -> ~tau_1+ h-)
    0.000000e+00 2 -2000015 -37 # br(~nu_taul -> ~tau_2+ h-)
    0.000000e+00 2 -1000015 -24 # br(~nu_taul -> ~tau_1+ w-)
    0.000000e+00 2 -2000015 -24 # br(~nu_taul -> ~tau_2+ w-)
So, is there something wrong in my command?

Thanks!

Yingpu

--
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/647138/+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/647138

You received this question notification because you asked the question.

Hi,

Setting the width to Auto should do the trick (since the BR line will be removed)
The only other options is to remove them manually by opening the file and remove them.
We do not have any way to remove those via a nice command.

Cheers,

Olivier

> On 9 Jul 2017, at 03:09, yingpu xiahou <email address hidden> wrote:
>
> Question #647138 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/647138
>
> Status: Answered => Open
>
> yingpu xiahou is still having a problem:
> Hi,
>
> Sorry. I am more interested in how to manually set Branching ratio
> currently. For PGS=1000016 as an example, how can I set all the bracing
> ratio of 1000016 to be zero?
>
> DECAY 1000016 0.000000e+00 # snu_taul decays
> 9.859945e-01 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
> 6.251296e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
> 0.000000e+00 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
> 0.000000e+00 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
> 7.754175e-03 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
> 0.000000e+00 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
> 0.000000e+00 2 -1000015 -37 # br(~nu_taul -> ~tau_1+ h-)
> 0.000000e+00 2 -2000015 -37 # br(~nu_taul -> ~tau_2+ h-)
> 0.000000e+00 2 -1000015 -24 # br(~nu_taul -> ~tau_1+ w-)
> 0.000000e+00 2 -2000015 -24 # br(~nu_taul -> ~tau_2+ w-)
>
> As I don’t want the Branching ratio to be non-zero.
>
> Thanks,
>
> Yingpu
>
>
> 在 2017年7月8日，19:57，Olivier Mattelaer <<email address hidden><mailto:<email address hidden>>> 写道：
>
> Your question #647138 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/647138
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> First the set command only modify the total width and do not touch to the Branching ratio.
> When a width is set on Auto then the code will replace the "Auto" value and the various branching ratio associated.
>
> So this explains that your SVN width still has BR information since we do not touch those.
> The same for your gluino case, since it looks like that the width is still on Auto (i.e. the computation of the width is not done yet and the BR information is still the old one)
>
> If you want to force the computation of the auto-width such that you can check the result before generating the events, you can type
> "compute_widths" after have set all the width that you want to compute on Auto.
> You can also do "compute_widths 1000021 1000016"
> Such that you do not have to set the width to auto in the first place.
>
> Cheers,
>
> Olivier
>
>
> On 8 Jul 2017, at 16:44, yingpu xiahou <<email address hidden><mailto:<email address hidden>>> wrote:
>
> New question #647138 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/647138
>
> Hi,
>
> Sorry for disturbing.
>
> I am using MG5_2_5_2 and MSSM model. The script are:
>
> import model sm
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> import model mssm
> generate p p > go go, (go > n1 t t~)
> add process p p > go go j, (go > n1 t t~)
> add process p p > go go j j, (go > n1 t t~)
> output xxx
> launch
> set DECAY 1000016 0 #wsn
> set DECAY 1000021 Auto #wgo
>
> I only want 1000021 to happen 3-body decay (go > n1 t t~).
>
> But, when I check the param_card.dat, I found all 1000021 decay are
> 2-body:
>
> DECAY 1000021 Auto # gluino decays
> 2.084542e-02 2 -1 1000001 # br(~g -> ~d_l db)
> 2.084542e-02 2 -1000001 1 # br(~g -> ~d_l* d )
> 5.070753e-02 2 -1 2000001 # br(~g -> ~d_r db)
> 5.070753e-02 2 -2000001 1 # br(~g -> ~d_r* d )
> 2.897878e-02 2 -2 1000002 # br(~g -> ~u_l ub)
> 2.897878e-02 2 -1000002 2 # br(~g -> ~u_l* u )
> 4.468728e-02 2 -2 2000002 # br(~g -> ~u_r ub)
> 4.468728e-02 2 -2000002 2 # br(~g -> ~u_r* u )
> 2.084542e-02 2 -3 1000003 # br(~g -> ~s_l sb)
> 2.084542e-02 2 -1000003 3 # br(~g -> ~s_l* s )
> 5.070753e-02 2 -3 2000003 # br(~g -> ~s_r sb)
> 5.070753e-02 2 -2000003 3 # br(~g -> ~s_r* s )
> 2.897878e-02 2 -4 1000004 # br(~g -> ~c_l cb)
> 2.897878e-02 2 -1000004 4 # br(~g -> ~c_l* c )
> 4.468728e-02 2 -4 2000004 # br(~g -> ~c_r cb)
> 4.468728e-02 2 -2000004 4 # br(~g -> ~c_r* c )
> 1.058402e-01 2 -5 1000005 # br(~g -> ~b_1 bb)
> 1.058402e-01 2 -1000005 5 # br(~g -> ~b_1* b )
> 5.565748e-02 2 -5 2000005 # br(~g -> ~b_2 bb)
> 5.565748e-02 2 -2000005 5 # br(~g -> ~b_2* b )
> 4.806428e-02 2 -6 1000006 # br(~g -> ~t_1 tb)
> 4.806428e-02 2 -1000006 6 # br(~g -> ~t_1* t )
> 0.000000e+00 2 -6 2000006 # br(~g -> ~t_2 tb)
> 0.000000e+00 2 -2000006 6 # br(~g -> ~t_2* t )
>
> And although I set 1000016 DECAY to 0, there is still some non-zero
> branching ratio appears:
>
> DECAY 1000016 0.000000e+00 # snu_taul decays
> 9.859945e-01 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)
> 6.251296e-03 2 16 1000023 # br(~nu_taul -> ~chi_20 nu_tau)
> 0.000000e+00 2 16 1000025 # br(~nu_taul -> ~chi_30 nu_tau)
> 0.000000e+00 2 16 1000035 # br(~nu_taul -> ~chi_40 nu_tau)
> 7.754175e-03 2 15 1000024 # br(~nu_taul -> ~chi_1+ tau-)
> 0.000000e+00 2 15 1000037 # br(~nu_taul -> ~chi_2+ tau-)
> 0.000000e+00 2 -1000015 -37 # br(~nu_taul -> ~tau_1+ h-)
> 0.000000e+00 2 -2000015 -37 # br(~nu_taul -> ~tau_2+ h-)
> 0.000000e+00 2 -1000015 -24 # br(~nu_taul -> ~tau_1+ w-)
> 0.000000e+00 2 -2000015 -24 # br(~nu_taul -> ~tau_2+ w-)
> So, is there something wrong in my command?
>
> Thanks!
>
> Yingpu
>
>
> --
> 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/647138/+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/647138
>
> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.