Issue with MadSpin decaying the Z in a BSM model

Asked by Emanuele Usai

Dear MadGraph Experts,

We are using MadGraph+MadSpin to generate the following process:
pp > Z’ > t T’ with T’ > t Z (t and Z to anything)
Where Z’ is a Z-like heavy resonance and T’ a heavy vector-like quark.
The MG model used can be found here [1] with the param [2] and run [3] cards.
We use this proc card [4], the relevant line is (rho0 is the Z’ and topo is the T'):
generate p p > rho0 > topo tt
The MadSpin card can be found here [5], the ony modification with respect to the standard madspin card is:
decay topo > z tt / rho+ rho- rho0 , z > all all, ( t > w+ b, w+ > all all ), (t~ > w- b~, w- > all all)
When I run MG+MadSpin, MadSpin crashes with this error [6] and log [7]. The compilation and generation of full square matrix element (with decay) steps are very slow as well (~30 min for both of them).
See the shell output [6].

I tried removing the decay of the Z from MadSpin, like this:
decay topo > z tt / rho+ rho- rho0 , ( t > w+ b, w+ > all all ), (t~ > w- b~, w- > all all)
and things go smooth and fast (but of course I’d like to decay the Z with MS as well).
If I let only the Z decay, like this:
decay topo > z tt / rho+ rho- rho0 , z > all all
MS crashes again.

So it seems there is a problem in the treatment of the Z.
Is there a problem in my MadSpin setup?

[1] https://github.com/emanueleusai/Vp_to_VLQ/tree/master/onerho
[2] https://github.com/emanueleusai/Vp_to_VLQ/blob/master/Single/Zp_Zt_t_madspin/param_card.dat
[3] https://github.com/emanueleusai/Vp_to_VLQ/blob/master/run_card.dat
[4] https://github.com/emanueleusai/Vp_to_VLQ/blob/master/Single/Zp_Zt_t_madspin/proc_card.dat
[5] https://github.com/emanueleusai/Vp_to_VLQ/blob/master/Single/Zp_Zt_t_madspin/madspin_card.dat
[6] https://github.com/emanueleusai/Vp_to_VLQ/blob/master/out2
[7] https://github.com/emanueleusai/Vp_to_VLQ/blob/master/out

Cheers,
Emanuele

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Pierre Artoisenet Edit question
Solved by:
Pierre Artoisenet
Solved:
Last query:
Last reply:
Revision history for this message
Pierre Artoisenet (partois) said :
#1

Hi Emanuele,

Can you load your model on Launchpad (perhaprs as a .tar.gz ?)

Pierre

Revision history for this message
Emanuele Usai (emanuele-usai) said :
#2

Hi Pierre,

Thanks for looking into my issue.
Here you can find the onerho model as a tar gzipped archive:
https://cernbox.cern.ch/public.php?service=files&t=e7edd6d6435a76d0f6d50bc548c6d3ed&download

Cheers,
Emanuele

Revision history for this message
Sébastien Brochet (blinkseb) said :
#3

Hi Emanuele,

You can look here for a possible patch for your issue : https://bugs.launchpad.net/mg5amcnlo/+bug/1435527

I'm having the exact same issue as you, and the attached patch should solve it.

Cheers,

Sébastien

Revision history for this message
Emanuele Usai (emanuele-usai) said :
#4

Hi Sébastien,

Awesome! The patch removes the error, now I can produce the events.

However, in the final LHE file, the decay products added by madspin have the wrong mother assignment.
For example here in this LHE file [1] you can see that the decay products of the Z are linked to the top, and the decay products of the top are linked to one of the decay products of the Z.

I'm applying the patch to version 2.2.3

[1] https://cernbox.cern.ch/public.php?service=files&t=b4ee53f808cfe983e8e8310bc66ee9a3&download

Cheers,
Emanuele

Revision history for this message
Pierre Artoisenet (partois) said :
#5

Hi Emanuele,

There is also a bug in the routine writing the event in the loco event file (decay_one_event_new)
I am having a look now

Pierre

Revision history for this message
Best Pierre Artoisenet (partois) said :
#6

Hi Emanuele,

I corrected the file decay.py in the branch lp:~maddevelopers/mg5amcnlo/MS_fix_bug

This should fix your problem. I did check that it solves the mother assignement in the case
(h> t t~). It would be great if you could do the check with your process as well.

Many thanks for your reports,

Pierre

Revision history for this message
Emanuele Usai (emanuele-usai) said :
#7

Hi Pierre,

I tested your fix with my model and it works correctly! The mother assignment is correct.
Thank you very much for addressing the issues.

Cheers,
Emanuele

Revision history for this message
Emanuele Usai (emanuele-usai) said :
#8

Thanks Pierre Artoisenet, that solved my question.