No decayed particles in the .lhe file using dm_s_spin1 model

there a typo in the proc_card.dat above: the corrected is below:

set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
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 DMsimp_UFO ../models/DMsimp_UFO/param_card.dat --bypass_check
generate p p > Y1
# Output processes to MadEvent directory
output MSSMD -nojpeg

Hi,

Do you look at the correct output?
In this case, you should have two directories in the Events one.
If you do not, you shoud have a crash (or you did not ask MasSpin to run at all)

Actually, for such simple process, I would advise to not use MadSpin but rather do
p p > Y1, Y1 > xd xd~

Cheers,

Olivier

> On 13 Mar 2019, at 15:32, mehdirahmani <email address hidden> wrote:
>
> New question #679170 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/679170
>
> Hello,
>
> I am using the spin 1 simplified dark matter model on feynrules website to generate some events. In this model Y1 is a spin 1 mediator and Xd is a Dirac dark matter spinor. I want to do the following for now:
>
> generate p p > Y1
> decay Y1 > Xd Xd
>
> When I try to generate events in madgraph the output .lhe file will contain the Y1s but there is no sign of Xd.
>
> Can you help me with this?
>
> below is all the info you may need to answer my question:
>
> Thank you,
>
>
>
> ===========================================
>
>
>
> Here is a link to the model file: http://feynrules.irmp.ucl.ac.be/attachment/wiki/DMsimp/dm_s_spin1.fr
>
> The model page: https://feynrules.irmp.ucl.ac.be/wiki/DMsimp
>
>
> Here is how my proc_card.dat looks like:
>
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> import model sm
> define p = 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 DMsimp_UFO ../models/DMsimp_UFO/param_card.dat --bypass_check
> generate p p > Y1
> # Output processes to MadEvent directory
> output MSSMD -nojpeg
> # This will create a directory PROC_$MODELNAME_$X
> # If you want to specify the path/name of the directory use
> # output PATH
>
> Here is how my madspin_card.dat looks like:
>
> set spinmode none
> # param_card.dat directory relative to work directory
> import model DMsimp_UFO ../models/DMsimp_UFO/param_card.dat --bypass_check
> # specify the decay for the final state particles
> decay Y1 > xd xd~
> # running the actual code
> launch
>
> here are my parameters:
>
> ## INFORMATION FOR MASS
> 5000521 4.000000e+00 # MXd
> 5000001 1.000000e+01 # MY1
>
> ## INFORMATION FOR DECAY
> DECAY 5000001 1.000000e+01 #Y1
> DECAY 5000521 0.000000 # Xd
>
> Block QNUMBERS 5000521 # Xd
> 1 0 # 3 times electric charge
> 2 2 # number of spin states (2S+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # Particle/Antiparticle distinction (0=own anti)
>
> Block QNUMBERS 5000521 # Xd
> 1 0 # 3 times electric charge
> 2 2 # number of spin states (2S+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # Particle/Antiparticle distinction (0=own anti)
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thank you for your response Oliver. Your suggestion solved the problem. However I still get an error when I do:

generate p p > Y1, Y1 > Xd Xd1~, Xd1 > Xd2 Y1, Y1 > mu+ mu-

Y1: Spin 1 mediator
Mass: 90 GeV
Width: 2 GeV

Xd1: Dirac Dark matter
Mass: 30 GeV
Width: 0.001 GeV

Xd2: Dirac Dark matter
Mass: 4 GeV
Width: 0

This is the error that I get from MadGraph:

INFO: Trying process: fd1 > fd2 gammad WEIGHTED<=2 DMV<=4

Command "import /afs/cern.ch/user/m/merahman/MadGraph5_DMEditedModel/MG5_aMC_v2_4_2/proc_card.dat" interrupted in sub-command:
"generate p p > gammaD, gammaD > Fd1 Fd1~, Fd1 > Fd2 gammaD, gammaD > mu+ mu-" with error:
NoDiagramException : No amplitudes generated from process Process: fd1 > fd2 gammad DMV=4 WEIGHTED=2. Please enter a valid process

quit

Hi,

1) You can not have a Y1 onshell decaying to a Xd1 onshell, decaying to a Y1 onshell.
This is just not allowed. You are not using the correct syntax.

2) For sub-decay (an onshell particle decaying into an onshell particles which in turns decay) you have to put paranthesis.

so one possibility (assuming that I assume correctly which particles are not onshell in your above chain would be)

> generate p p > gammaD, (gammaD > Fd1 Fd1~, Fd1 > Fd2 mu+ mu)

3) Your error seems actually different:
> NoDiagramException : No amplitudes generated from process Process: fd1 > fd2 gammad DMV=4 WEIGHTED=2. Please enter a valid process

This seems to indicate that you model does not have such kind of vertex.

Cheers,

Olivier

> On 22 Mar 2019, at 15:47, mehdirahmani <email address hidden> wrote:
>
> Question #679170 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/679170
>
> Status: Answered => Open
>
> mehdirahmani is still having a problem:
> Thank you for your response Oliver. Your suggestion solved the problem.
> However I still get an error when I do:
>
> generate p p > Y1, Y1 > Xd Xd1~, Xd1 > Xd2 Y1, Y1 > mu+ mu-
>
> Y1: Spin 1 mediator
> Mass: 90 GeV
> Width: 2 GeV
>
> Xd1: Dirac Dark matter
> Mass: 30 GeV
> Width: 0.001 GeV
>
> Xd2: Dirac Dark matter
> Mass: 4 GeV
> Width: 0
>
> This is the error that I get from MadGraph:
>
> INFO: Trying process: fd1 > fd2 gammad WEIGHTED<=2 DMV<=4
>
> Command "import /afs/cern.ch/user/m/merahman/MadGraph5_DMEditedModel/MG5_aMC_v2_4_2/proc_card.dat" interrupted in sub-command:
> "generate p p > gammaD, gammaD > Fd1 Fd1~, Fd1 > Fd2 gammaD, gammaD > mu+ mu-" with error:
> NoDiagramException : No amplitudes generated from process Process: fd1 > fd2 gammad DMV=4 WEIGHTED=2. Please enter a valid process
>
> quit
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thank you for your response.

Second gammaD (Y1) can be offshell. I used your syntax and got the same error:

INFO: Trying process: fd1 > fd2 mu+ mu- WEIGHTED<=4 DMV<=4
Command "import /afs/cern.ch/user/m/merahman/MadGraph5_DMEditedModel/MG5_aMC_v2_4_2/proc_card.dat" interrupted in sub-command:
"generate p p > gammaD, (gammaD > Fd1 Fd1~, Fd1 > Fd2 mu+ mu-)" with error:
NoDiagramException : No amplitudes generated from process Process: fd1 > fd2 mu+ mu- DMV=4 WEIGHTED=4. Please enter a valid process
quit

I see your point about vertex but I have defined that vertex in my model:

L1Fd1Fd2gammaD := Fd1bar.Ga[mu].gFd1Fd2.Fd2 gammaD[mu];

Could it be that my parameters are off?

Y1 (gammaD): Spin 1 mediator
Mass: 90 GeV
Width: 2 GeV

Xd1 (Fd1): Dirac Dark matter
Mass: 30 GeV
Width: 0.001 GeV

Xd2 (Fd2): Dirac Dark matter
Mass: 4 GeV
Width: 0

gFd1Fd2 coupling constant: 0.25

I encountered the same problem at some point when I was trying to do
p p > gammaD, gammaD > Fd1 Fd1~
with gammaD mass of 10 GeV but increasing the mass solved the problem.

Thank you,

Hi,

> I see your point about vertex but I have defined that vertex in my
> model:

I do not know FeynRules (or other similar code), so I can not comment on that.
I can only claim that such error occurs when your diagram does not exists.

> I encountered the same problem at some point when I was trying to do
> p p > gammaD, gammaD > Fd1 Fd1~
> with gammaD mass of 10 GeV but increasing the mass solved the problem.

You can not have the SAME error since the generation of the Feynman Diagram does not depend
of the benchmark.
Now I can see a lot of issues at later stage of the computation (breaking of perturbativity and impossibiity for gammaD to be onshell due to the mass of Fd1 to only cite two)

Cheers,

Olivier

> On 22 Mar 2019, at 16:37, mehdirahmani <email address hidden> wrote:
>
> Question #679170 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/679170
>
> Status: Answered => Open
>
> mehdirahmani is still having a problem:
> Thank you for your response.
>
> Second gammaD (Y1) can be offshell. I used your syntax and got the same
> error:
>
> INFO: Trying process: fd1 > fd2 mu+ mu- WEIGHTED<=4 DMV<=4
> Command "import /afs/cern.ch/user/m/merahman/MadGraph5_DMEditedModel/MG5_aMC_v2_4_2/proc_card.dat" interrupted in sub-command:
> "generate p p > gammaD, (gammaD > Fd1 Fd1~, Fd1 > Fd2 mu+ mu-)" with error:
> NoDiagramException : No amplitudes generated from process Process: fd1 > fd2 mu+ mu- DMV=4 WEIGHTED=4. Please enter a valid process
> quit
>
> I see your point about vertex but I have defined that vertex in my
> model:
>
> L1Fd1Fd2gammaD := Fd1bar.Ga[mu].gFd1Fd2.Fd2 gammaD[mu];
>
> Could it be that my parameters are off?
>
> Y1 (gammaD): Spin 1 mediator
> Mass: 90 GeV
> Width: 2 GeV
>
> Xd1 (Fd1): Dirac Dark matter
> Mass: 30 GeV
> Width: 0.001 GeV
>
> Xd2 (Fd2): Dirac Dark matter
> Mass: 4 GeV
> Width: 0
>
> gFd1Fd2 coupling constant: 0.25
>
>
> I encountered the same problem at some point when I was trying to do
> p p > gammaD, gammaD > Fd1 Fd1~
> with gammaD mass of 10 GeV but increasing the mass solved the problem.
>
> Thank you,
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

I am not sure if was the same exact error but some error occurred when the dark photon mass was set to 10 GeV.
But can you offer some suggestion on how I should go about solving these issues? I am correct to understand that I have some problems in my model file?

Thanks,

> But can you offer some suggestion on how I should go about solving these issues? I am correct to understand that I have some problems in my model file?

That's my bet yes. Now I can not help more than that since I have not used FeynRules since 2010.

Cheers,

Olivier

> On 22 Mar 2019, at 17:02, mehdirahmani <email address hidden> wrote:
>
> Question #679170 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/679170
>
> Status: Answered => Open
>
> mehdirahmani is still having a problem:
> I am not sure if was the same exact error but some error occurred when the dark photon mass was set to 10 GeV.
> But can you offer some suggestion on how I should go about solving these issues? I am correct to understand that I have some problems in my model file?
>
> Thanks,
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.