No valid decay for particle 6 with mad spin

Asked by David Uzan

Hello,

I am trying to simulate an unphysical process as a test of p p > t t~, where the top mass has been changed to 4.7.
I also change the following:
aewm1=1e8
mb = 1
mz=1
ww=1e-3
wt =1e-3

I am using the dim6top_LO model, with the full process above written as:
p p > t t~ DIM6=1 FCNC=0 QED=0.

The problem arises when decaying the "top" particles using Madspin t > b w+, w+ > l+ vl.
I receive the error "no valid decay for 6. no 2 body decay for that particle. (three body are not supported by madspin)".
I have tried changing the values written above, and also tried writing DIM6=1 FCNC=0 after each decay but it does not seem to work.

The odd thing is that when doing this decay directly :
p p > t t~ DIM6=1 FCNC=0 QED=0, (t > b l+ vl DIM6=1 FCNC=0), (t~ > b~ l- vl~ DIM6=1 FCNC=0)
the simulation seems to work fine.

Is Madpsin limited in some way to make this process unavailable, or do I miss something when writing the Madpsin card?

Thank you in advance,
David

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

The issue is that the decay that you asked is unphysical.
Since you do not specify the branching ratio for your "top", madspin request madwidth to do the computation.
MadWidth detects that the decay is below the QCD scale and discard the decay channel since such computation is meaningless.

If you ask the decay directly, Madgraph does not test if the decay does makes sense or not and therefore return a (non physical) result. I would think to check if they are a way to easily detects such type of decay and forbide them as well within madgraph.

Cheers,

Olivier

Revision history for this message
David Uzan (david-uz) said :
#2

Thank you for your answer, but I don't understand your suggestion.
Do you suggest forbidding my decay within madgrpah, as in not looking at this decay at all?
If I misunderstood what you meant, I have a different question, Is there a way to change the BR directly?
I am changing the decay width of the "top", but does it affect Madspin aswell?

David

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi,

Yes, I would not use/trust any result for that syntax since you are going to be dominated by non perturbative physics.

For MG5aMC, the only input used is the total width (all BR information will be ignored).

MadSpin is only going to use the provided total width if you also have the Branching Ratio information.
Otherwise, it will compute the total width (and all the BR) and use that one.
So by providing the width/BR for all particle in madspin, you can bypass the issue that MadSpin detects that the decay is non perturbative.

Olivier

Revision history for this message
David Uzan (david-uz) said :
#4

Thank you again.

Lastly, I couldn't find a way to change the BR/width specifically in madspin.
Is it done through the param_card or is it based on the model I use?

I have changed the width of the "top" in the param_card and the error still persists, so I figured it is not the former.

David

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#5

Yes,

It is done via the param_card according to the community standard for that file (see SLHA1: Susy LesHouches Accord)
You have to have the BR information included for all particles that have a line decay in the madspin_card to avoid that madspin triggers the automatic computation of the width.

But can you try to generate a physical process as a test, running such unphysical process leads you to unphysical regime where all code are expected to either crash (best case scenario) or to return non-sense. Trying to avoid those crash is just a waste of time.
I do not really understand what you are hoping to achieve here.

Cheers,

Olivier

Revision history for this message
David Uzan (david-uz) said :
#6

Thank you very much, this does help me run the simulation.

For the reason for me doing this, I am using this process as a probe of the polarization and spin correlation of both of the "top" particles. From the mass of these "tops", I am trying to look at the decaying bottoms, and as far as I understand it, this process fits my needs.

I used this process for actual tops as well, and using madspin for the decay resulted in better output. This is why I am trying to make it work through madspin.

Best,
David

Can you help with this problem?

Provide an answer of your own, or ask David Uzan for more information if necessary.

To post a message you must log in.