TwinMSSM quark mass matrices and CH/UFO output

[terryp]

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:

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 elements

Code: 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" .

When attempting to import the model into MG5 (2.6.3.2), I get the following error:

Code: Select all

Command "import model TwinMSSM/" interrupted with error:
NameError : name 'Tt1' is not defined

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

[FStaub]

Sorry, I have no time to look at your fields. However, I attach a bunch of Twin models which I implemented some time ago. Those are not part of the public package because they aren't checked in all detail, but I hope they are helpful anyways.

Cheers,
Florian