TwinMSSM quark mass matrices and CH/UFO output
Posted: 28. Jan 2019, 16:22
Dear Florian,
I am trying to make a Twin MSSM model for use in micrOmegas (and eventually MadGraph). I am a bit stuck at the model checking and output stages. The model files are attached to this message. I have duplicated the existing MSSM model content, filling in the missing values from the SARAH/Models/particles.m and SARAH/Models/parameters.m files.
The first error I haven't been able to figure out is that the twin quark mass matrices are zero while the twin lepton matrix is nonzero. During CheckModel, I get the following:
No errors are generated during MakeVertexList[].
When MakeCHep[IncludeLoopDecyas -> False], I get the following error:
When attempting to import the model into MG5 (2.6.3.2), I get the following error:
It appears that sums are not being performed during the output for the four gluon vertex. The vertexListVVVV.m file lists the four-gluon vertex as identical to the MSSM output and no such errors pop up when running the same commands with SARAH's MSSM model file.
It takes me about twenty minutes to run through those commands, so I've attached the full worksheet with output as a PDF.
Any guidance you can give me would be greatly appreciated. I am using Mathematica 10.0.2 and SARAH 4.14.0.
Thanks,
Terry
I am trying to make a Twin MSSM model for use in micrOmegas (and eventually MadGraph). I am a bit stuck at the model checking and output stages. The model files are attached to this message. I have duplicated the existing MSSM model content, filling in the missing values from the SARAH/Models/particles.m and SARAH/Models/parameters.m files.
The first error I haven't been able to figure out is that the twin quark mass matrices are zero while the twin lepton matrix is nonzero. During CheckModel, I get the following:
Code: Select all
Checking matrix for {{FTwindL},{conj[FTwindR]}}
MassMatrix::OnlyZero: The mass matrix for {{FTwindL},{conj[FTwindR]}} contains only zero elements
Checking matrix for {{FTwinuL},{conj[FTwinuR]}}
MassMatrix::OnlyZero: The mass matrix for {{FTwinuL},{conj[FTwinuR]}} contains only zero elementsCode: Select all
In: MassMatrix[FTwind]
Out: {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
In: MassMatrix[FTwinu]
Out: {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}
In: MassMatrix[FTwine]
Out: {{(Twinvd TwinYe[1,1])/Sqrt[2],(Twinvd TwinYe[2,1])/Sqrt[2],(Twinvd TwinYe[3,1])/Sqrt[2]},{(Twinvd TwinYe[1,2])/Sqrt[2],(Twinvd TwinYe[2,2])/Sqrt[2],(Twinvd TwinYe[3,2])/Sqrt[2]},{(Twinvd TwinYe[1,3])/Sqrt[2],(Twinvd TwinYe[2,3])/Sqrt[2],(Twinvd TwinYe[3,3])/Sqrt[2]}}No errors are generated during MakeVertexList[].
When MakeCHep[IncludeLoopDecyas -> False], I get the following error:
Code: Select all
StringJoin::string: "String expected at position 2 in "|m1.m2*m3.m4"<>(-sum[j1,1,8,fSU3[j1,Tt1,Tt4] fSU3[j1,Tt2,Tt3]] + sum[j1,1,8,fSU3[j1,Tt1,Tt2] fSU3[j1,Tt3,Tt4]])/(sum[j1,1,8,fSU3[j1,Tt1,Tt4] fSU3[j1,Tt2,Tt3]] + sum[j1,1,8,fSU3[j1,Tt1,Tt3] fSU3[j1,Tt2,Tt4]])<>"m1.m3*m2.m4 "<>(1/16 (-sum[j1,1,8,fSU3[j1,Tt1,Tt3] fSU3[j1,Tt2,Tt4]] - sum[j1,1,8,fSU3[j1,Tt1,Tt2] fSU3[j1,Tt3,Tt4]]))<>"m1.m4*m2.m3" .Code: Select all
Command "import model TwinMSSM/" interrupted with error:
NameError : name 'Tt1' is not definedIt takes me about twenty minutes to run through those commands, so I've attached the full worksheet with output as a PDF.
Any guidance you can give me would be greatly appreciated. I am using Mathematica 10.0.2 and SARAH 4.14.0.
Thanks,
Terry