Problem with applying MadSpin for decay chains

Asked by Afiq Aize

Hi,

I have generated some events using the HC model: p p > x0 for a couple of points. I then tried to decay this resonance to ttbar using MadSpin; I have a card which goes like [1] which I run with [2]. So it ran fine, but I found that the events_decayed.lhe.gz contain only the x0 > t t~ decay and not the subsequent desired t > w b and w > l vl decays. Then I naively gave the decayed LHE file to MadSpin for the subsequent decays but it complained [3].

Then I tried modifying the card, specifying the whole chain in [4], but this only ends me up with a new complain [5].

May I know the proper way of get the desired final state? All the decays in the chain are two-body decays so as far as I understand MadSpin should have no issue handling them.

Cheers,
Afiq

[1]
decay x0 > t t~
decay t > w+ b, w+ > l+ vl
decay t~ > w- b~, w- > l- vl~

import x0-lo/Events/ma500/events.lhe.gz
launch

import x0-lo/Events/mh500/events.lhe.gz
launch

[...]

[2] ./MadSpin/madspin < x0tt_decay.txt

[3]
Command "import output/x0-lo/Events/mh500/events_decayed.lhe.gz" interrupted with error:
InvalidCmd : This event file was already decayed by MS. This is not possible to add to it a second decay
Command "launch" interrupted with error:
AttributeError : 'NoneType' object has no attribute 'get'
Please report this bug to developers

           More information is found in 'MS_debug'.

[4] decay x0 > t t~, t > w+ b, w+ > l+ vl, t~ > w- b~, w- > l- vl~

[5]
Command "launch" interrupted with error:
ValueError : invalid literal for float(): Error#13 in genps_madspin.f

Please report this bug to developers

           More information is found in 'MS_debug'.

           Please attach this file to your report.

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

you miss the parenthesis for sub-decay:

> [4] decay x0 > t t~, t > w+ b, w+ > l+ vl, t~ > w- b~, w- > l- vl~
should be
> [4] decay x0 > t t~, (t > w+ b, w+ > l+ vl),( t~ > w- b~, w- > l- vl~)

> On Jan 7, 2016, at 18:07, Afiq Aize <email address hidden> wrote:
>
> New question #280715 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/280715
>
> Hi,
>
> I have generated some events using the HC model: p p > x0 for a couple of points. I then tried to decay this resonance to ttbar using MadSpin; I have a card which goes like [1] which I run with [2]. So it ran fine, but I found that the events_decayed.lhe.gz contain only the x0 > t t~ decay and not the subsequent desired t > w b and w > l vl decays. Then I naively gave the decayed LHE file to MadSpin for the subsequent decays but it complained [3].
>
> Then I tried modifying the card, specifying the whole chain in [4], but this only ends me up with a new complain [5].
>
> May I know the proper way of get the desired final state? All the decays in the chain are two-body decays so as far as I understand MadSpin should have no issue handling them.
>
> Cheers,
> Afiq
>
> [1]
> decay x0 > t t~
> decay t > w+ b, w+ > l+ vl
> decay t~ > w- b~, w- > l- vl~
>
> import x0-lo/Events/ma500/events.lhe.gz
> launch
>
> import x0-lo/Events/mh500/events.lhe.gz
> launch
>
> [...]
>
> [2] ./MadSpin/madspin < x0tt_decay.txt
>
> [3]
> Command "import output/x0-lo/Events/mh500/events_decayed.lhe.gz" interrupted with error:
> InvalidCmd : This event file was already decayed by MS. This is not possible to add to it a second decay
> Command "launch" interrupted with error:
> AttributeError : 'NoneType' object has no attribute 'get'
> Please report this bug to developers
>
> More information is found in 'MS_debug'.
>
> [4] decay x0 > t t~, t > w+ b, w+ > l+ vl, t~ > w- b~, w- > l- vl~
>
> [5]
> Command "launch" interrupted with error:
> ValueError : invalid literal for float(): Error#13 in genps_madspin.f
>
> Please report this bug to developers
>
> More information is found in 'MS_debug'.
>
> Please attach this file to your report.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Afiq Aize (watashinokonoka) said :
#2

Thanks Olivier Mattelaer, that solved my question.