Error using MakeVevacious for THDMCPV

Questions concerning the interface to Vevacious

Moderator: benoleary

Post Reply
AWlotzka

Error using MakeVevacious for THDMCPV

Post by AWlotzka » 29. Apr 2016, 09:07

Dear Florian Staub,

using the command

Code: Select all

MakeVevacious[];
for the THDM-CPV model results in Mathematica-errors

Code: Select all

Part::partw: Part 1 of {} does not exist. >>

Part::partw: Part 1 of {} does not exist. >>

Part::partw: Part 2 of {} does not exist. >>

General::stop: Further output of Part::partw will be suppressed during this calculation. >>
Nevertheless, the THDMCPV.vin is generated. For other models this error does not occur.
Is it of any importance?

Thanks,
Alexander Wlotzka
Top
FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: Error using MakeVevacious for THDMCPV

Post by FStaub » 29. Apr 2016, 12:13

Hi,

this problem is related to the issue with HOM4PS discussed in http://stauby.de/sarah_userforum/viewtopic.php?f=8&t=18 and actually, I think there might be a solution by changing the model file to

Code: Select all

(*
 DEFINITION[GaugeES][Phases]= 
{    {H1p, Exp[I eta]},
     {H10, Exp[I eta]}
    }; 
 *)

(* ----- VEVs ---- *)

DEFINITION[EWSB][VEVs]= 
{    {H10, {v1R, 1/Sqrt[2]}, {v1I, 1/Sqrt[2]}, {sigma1, \[ImaginaryI]/Sqrt[2]},{phi1, 1/Sqrt[2]}},
     {H20, {v2R, 1/Sqrt[2]}, {v2I, 1/Sqrt[2]}, {sigma2, \[ImaginaryI]/Sqrt[2]},{phi2, 1/Sqrt[2]}}     };
i.e. defining a real and imaginary part for the VEVs only and forgetting about the overall phase. In addition, it is necessary to write in parameters.m

Code: Select all

{v1R,        { LesHouches->{HMIX,1000} }},
{v1I,        { LesHouches->{HMIX,1001} }}, 
{v2R,        { LesHouches->{HMIX,1002} }},
{v2I,        { LesHouches->{HMIX,1003} }},


The only subtlety is that SPheno won't work with that ansatz correctly because the phases in the Higgs vertices (which are not needed for Vevacious because it uses only tadpoles & mass matrices) would be missing. This, it would be necessary to translate the SPheno conventions to the Vevacious ones before giving the SLHA file to Vevacious. That should be possible by editing InputOutput_THDMCPV.f90 written by SARAH for the standard THDMCPV as

Code: Select all

Write(io_L,106) "Block HMIX Q=",Q,"# (Renormalization Scale)" 
...
Write(io_L,104) 102,Real(vd,dp), "# vd" 
Write(io_L,104) 103,Real(vu,dp), "# vu" 
Write(io_L,104) 1001,Abs(vd)*Cos(eta), "# v1R" 
Write(io_L,104) 1002,Abs(vd)*Sin(eta), "# v1I" 
Write(io_L,104) 1003,Real(vu,dp), "# v2R" 
Write(io_L,104) 1004,0._dp, "# v2I" 
..
At least, I hope so, because I didn't test anything of that above myself right now. So, please be a bit careful.

Cheers,
Florian
Top
AWlotzka

Re: Error using MakeVevacious for THDMCPV

Post by AWlotzka » 29. Apr 2016, 15:40

Hello,

that sounds good, I'll try it out and let you know whether it works.

Thanks for the moment!
Alexander
Top
AWlotzka

Re: Error using MakeVevacious for THDMCPV

Post by AWlotzka » 2. May 2016, 10:50

Hi again,

I performed the steps you suggested, but unfortunately there is still something going wrong. I still get the same error as before from HOM4PS2, here is the content of the current_tadpoles.in:

Code: Select all

{ (a*(1.96666740E+06)) + (b*(1.96666740E+06)) + (-1.*c*(4.00000000E+04)) + (-1.*d*(4.00000000E+04)) + (a^3*(1.00000000E-01)) + (3.*a^2*b*(1.00000000E-01)) + (3.*a*b^2*(1.00000000E-01)) + (b^3*(1.00000000E-01)) + (0.5*a*c^2*(1.10000000E+00)) + (0.5*b*c^2*(1.10000000E+00)) + (a*c*d*(1.10000000E+00)) + (b*c*d*(1.10000000E+00)) + (0.5*a*d^2*(1.10000000E+00)) + (0.5*b*d^2*(1.10000000E+00)) + (0.5*a*c^2*(-5.00000000E-01)) + (0.5*b*c^2*(-5.00000000E-01)) + (a*c*d*(-5.00000000E-01)) + (b*c*d*(-5.00000000E-01)) + (0.5*a*d^2*(-5.00000000E-01)) + (0.5*b*d^2*(-5.00000000E-01)) + (0.5*a*c^2*(5.00000000E-01)) + (0.5*b*c^2*(5.00000000E-01)) + (a*c*d*(5.00000000E-01)) + (b*c*d*(5.00000000E-01)) + (0.5*a*d^2*(5.00000000E-01)) + (0.5*b*d^2*(5.00000000E-01));
(c*(-7.09137478E+03)) + (d*(-7.09137478E+03)) + (-1.*a*(4.00000000E+04)) + (-1.*b*(4.00000000E+04)) + (c^3*(1.30000000E-01)) + (3.*c^2*d*(1.30000000E-01)) + (3.*c*d^2*(1.30000000E-01)) + (d^3*(1.30000000E-01)) + (0.5*a^2*c*(1.10000000E+00)) + (a*b*c*(1.10000000E+00)) + (0.5*b^2*c*(1.10000000E+00)) + (0.5*a^2*d*(1.10000000E+00)) + (a*b*d*(1.10000000E+00)) + (0.5*b^2*d*(1.10000000E+00)) + (0.5*a^2*c*(-5.00000000E-01)) + (a*b*c*(-5.00000000E-01)) + (0.5*b^2*c*(-5.00000000E-01)) + (0.5*a^2*d*(-5.00000000E-01)) + (a*b*d*(-5.00000000E-01)) + (0.5*b^2*d*(-5.00000000E-01)) + (0.5*a^2*c*(5.00000000E-01)) + (a*b*c*(5.00000000E-01)) + (0.5*b^2*c*(5.00000000E-01)) + (0.5*a^2*d*(5.00000000E-01)) + (a*b*d*(5.00000000E-01)) + (0.5*b^2*d*(5.00000000E-01)); }
Moreover, now I get another error:

Code: Select all

Traceback (most recent call last):
  File "./Vevacious.py", line 29, in <module>
    EffectivePotential = VPD.LoopAndThermalCorrectedPotential )
  File "/itp/scratch/wlotzka/project_baryogenesis/Vevacious/bin/VevaciousParameterDependent.py", line 816, in __init__
    temperatureValue = currentTemperature )
  File "/itp/scratch/wlotzka/project_baryogenesis/Vevacious/bin/VevaciousParameterDependent.py", line 719, in __init__
    self.functionAtOrigin = FunctionToScale( 0.0, 0.0, 0.0 )
TypeError: LoopAndThermalCorrectedPotential() takes exactly 5 arguments (3 given)
The definition of the Function which caused the error in the file VevaciousParameterDependent.py is

Code: Select all

def LoopAndThermalCorrectedPotential( a, b, c, d, temperatureValue ):
If needed, here is the THDMCPV.vin which I generated with SARAH and the modified model and parameter files:

Code: Select all

<!--  
 This model file was automatically created by SARAH version4.8.5  
 SARAH References: arXiv:0806.0538, 0909.2863, 1002.0840, 1207.0906, 1309.7223 
 (c) Florian Staub, 2013   
 
 File created at 9:54 on 2.5.2016  
 
                                      -->
 
 
<Vevacious_stuff> 
<input_vevs v1I="SLHA::HMIX[1001]" v1R="SLHA::HMIX[1000]" v2I="SLHA::HMIX[1003]" v2R="SLHA::HMIX[1002]" /> 

<polynomial_part> 
(0.5*v1I^2*SLHA::HMIX[20.]) 
  + (v1I*v1R*SLHA::HMIX[20.]) 
  + (0.5*v1R^2*SLHA::HMIX[20.]) 
  + (0.5*v2I^2*SLHA::HMIX[21.]) 
  + (v2I*v2R*SLHA::HMIX[21.]) 
  + (0.5*v2R^2*SLHA::HMIX[21.]) 
  + (-1.*v1I*v2I*SLHA::HMIX[22.]) 
  + (-1.*v1R*v2I*SLHA::HMIX[22.]) 
  + (-1.*v1I*v2R*SLHA::HMIX[22.]) 
  + (-1.*v1R*v2R*SLHA::HMIX[22.]) 
  + (0.25*v1I^4*SLHA::HMIX[31.]) 
  + (v1I^3*v1R*SLHA::HMIX[31.]) 
  + (1.5*v1I^2*v1R^2*SLHA::HMIX[31.]) 
  + (v1I*v1R^3*SLHA::HMIX[31.]) 
  + (0.25*v1R^4*SLHA::HMIX[31.]) 
  + (0.25*v2I^4*SLHA::HMIX[32.]) 
  + (v2I^3*v2R*SLHA::HMIX[32.]) 
  + (1.5*v2I^2*v2R^2*SLHA::HMIX[32.]) 
  + (v2I*v2R^3*SLHA::HMIX[32.]) 
  + (0.25*v2R^4*SLHA::HMIX[32.]) 
  + (0.25*v1I^2*v2I^2*SLHA::HMIX[33.]) 
  + (0.5*v1I*v1R*v2I^2*SLHA::HMIX[33.]) 
  + (0.25*v1R^2*v2I^2*SLHA::HMIX[33.]) 
  + (0.5*v1I^2*v2I*v2R*SLHA::HMIX[33.]) 
  + (v1I*v1R*v2I*v2R*SLHA::HMIX[33.]) 
  + (0.5*v1R^2*v2I*v2R*SLHA::HMIX[33.]) 
  + (0.25*v1I^2*v2R^2*SLHA::HMIX[33.]) 
  + (0.5*v1I*v1R*v2R^2*SLHA::HMIX[33.]) 
  + (0.25*v1R^2*v2R^2*SLHA::HMIX[33.]) 
  + (0.25*v1I^2*v2I^2*SLHA::HMIX[34.]) 
  + (0.5*v1I*v1R*v2I^2*SLHA::HMIX[34.]) 
  + (0.25*v1R^2*v2I^2*SLHA::HMIX[34.]) 
  + (0.5*v1I^2*v2I*v2R*SLHA::HMIX[34.]) 
  + (v1I*v1R*v2I*v2R*SLHA::HMIX[34.]) 
  + (0.5*v1R^2*v2I*v2R*SLHA::HMIX[34.]) 
  + (0.25*v1I^2*v2R^2*SLHA::HMIX[34.]) 
  + (0.5*v1I*v1R*v2R^2*SLHA::HMIX[34.]) 
  + (0.25*v1R^2*v2R^2*SLHA::HMIX[34.]) 
  + (0.25*v1I^2*v2I^2*SLHA::HMIX[35.]) 
  + (0.5*v1I*v1R*v2I^2*SLHA::HMIX[35.]) 
  + (0.25*v1R^2*v2I^2*SLHA::HMIX[35.]) 
  + (0.5*v1I^2*v2I*v2R*SLHA::HMIX[35.]) 
  + (v1I*v1R*v2I*v2R*SLHA::HMIX[35.]) 
  + (0.5*v1R^2*v2I*v2R*SLHA::HMIX[35.]) 
  + (0.25*v1I^2*v2R^2*SLHA::HMIX[35.]) 
  + (0.5*v1I*v1R*v2R^2*SLHA::HMIX[35.]) 
  + (0.25*v1R^2*v2R^2*SLHA::HMIX[35.]) 
  + (0.5*v1I^3*v2I*SLHA::HMIX[36.]) 
  + (1.5*v1I^2*v1R*v2I*SLHA::HMIX[36.]) 
  + (1.5*v1I*v1R^2*v2I*SLHA::HMIX[36.]) 
  + (0.5*v1R^3*v2I*SLHA::HMIX[36.]) 
  + (0.5*v1I^3*v2R*SLHA::HMIX[36.]) 
  + (1.5*v1I^2*v1R*v2R*SLHA::HMIX[36.]) 
  + (1.5*v1I*v1R^2*v2R*SLHA::HMIX[36.]) 
  + (0.5*v1R^3*v2R*SLHA::HMIX[36.]) 
  + (0.5*v1I*v2I^3*SLHA::HMIX[37.]) 
  + (0.5*v1R*v2I^3*SLHA::HMIX[37.]) 
  + (1.5*v1I*v2I^2*v2R*SLHA::HMIX[37.]) 
  + (1.5*v1R*v2I^2*v2R*SLHA::HMIX[37.]) 
  + (1.5*v1I*v2I*v2R^2*SLHA::HMIX[37.]) 
  + (1.5*v1R*v2I*v2R^2*SLHA::HMIX[37.]) 
  + (0.5*v1I*v2R^3*SLHA::HMIX[37.]) 
  + (0.5*v1R*v2R^3*SLHA::HMIX[37.]) 
 
</polynomial_part>

<tadpoles> 
{ 
(v1I*SLHA::HMIX[20.]) 
+(v1R*SLHA::HMIX[20.]) 
+(-1.*v2I*SLHA::HMIX[22.]) 
+(-1.*v2R*SLHA::HMIX[22.]) 
+(v1I^3*SLHA::HMIX[31.]) 
+(3.*v1I^2*v1R*SLHA::HMIX[31.]) 
+(3.*v1I*v1R^2*SLHA::HMIX[31.]) 
+(v1R^3*SLHA::HMIX[31.]) 
+(0.5*v1I*v2I^2*SLHA::HMIX[33.]) 
+(0.5*v1R*v2I^2*SLHA::HMIX[33.]) 
+(v1I*v2I*v2R*SLHA::HMIX[33.]) 
+(v1R*v2I*v2R*SLHA::HMIX[33.]) 
+(0.5*v1I*v2R^2*SLHA::HMIX[33.]) 
+(0.5*v1R*v2R^2*SLHA::HMIX[33.]) 
+(0.5*v1I*v2I^2*SLHA::HMIX[34.]) 
+(0.5*v1R*v2I^2*SLHA::HMIX[34.]) 
+(v1I*v2I*v2R*SLHA::HMIX[34.]) 
+(v1R*v2I*v2R*SLHA::HMIX[34.]) 
+(0.5*v1I*v2R^2*SLHA::HMIX[34.]) 
+(0.5*v1R*v2R^2*SLHA::HMIX[34.]) 
+(0.5*v1I*v2I^2*SLHA::HMIX[35.]) 
+(0.5*v1R*v2I^2*SLHA::HMIX[35.]) 
+(v1I*v2I*v2R*SLHA::HMIX[35.]) 
+(v1R*v2I*v2R*SLHA::HMIX[35.]) 
+(0.5*v1I*v2R^2*SLHA::HMIX[35.]) 
+(0.5*v1R*v2R^2*SLHA::HMIX[35.]) 
+(1.5*v1I^2*v2I*SLHA::HMIX[36.]) 
+(3.*v1I*v1R*v2I*SLHA::HMIX[36.]) 
+(1.5*v1R^2*v2I*SLHA::HMIX[36.]) 
+(1.5*v1I^2*v2R*SLHA::HMIX[36.]) 
+(3.*v1I*v1R*v2R*SLHA::HMIX[36.]) 
+(1.5*v1R^2*v2R*SLHA::HMIX[36.]) 
+(0.5*v2I^3*SLHA::HMIX[37.]) 
+(1.5*v2I^2*v2R*SLHA::HMIX[37.]) 
+(1.5*v2I*v2R^2*SLHA::HMIX[37.]) 
+(0.5*v2R^3*SLHA::HMIX[37.]) 
; 
(v2I*SLHA::HMIX[21.]) 
+(v2R*SLHA::HMIX[21.]) 
+(-1.*v1I*SLHA::HMIX[22.]) 
+(-1.*v1R*SLHA::HMIX[22.]) 
+(v2I^3*SLHA::HMIX[32.]) 
+(3.*v2I^2*v2R*SLHA::HMIX[32.]) 
+(3.*v2I*v2R^2*SLHA::HMIX[32.]) 
+(v2R^3*SLHA::HMIX[32.]) 
+(0.5*v1I^2*v2I*SLHA::HMIX[33.]) 
+(v1I*v1R*v2I*SLHA::HMIX[33.]) 
+(0.5*v1R^2*v2I*SLHA::HMIX[33.]) 
+(0.5*v1I^2*v2R*SLHA::HMIX[33.]) 
+(v1I*v1R*v2R*SLHA::HMIX[33.]) 
+(0.5*v1R^2*v2R*SLHA::HMIX[33.]) 
+(0.5*v1I^2*v2I*SLHA::HMIX[34.]) 
+(v1I*v1R*v2I*SLHA::HMIX[34.]) 
+(0.5*v1R^2*v2I*SLHA::HMIX[34.]) 
+(0.5*v1I^2*v2R*SLHA::HMIX[34.]) 
+(v1I*v1R*v2R*SLHA::HMIX[34.]) 
+(0.5*v1R^2*v2R*SLHA::HMIX[34.]) 
+(0.5*v1I^2*v2I*SLHA::HMIX[35.]) 
+(v1I*v1R*v2I*SLHA::HMIX[35.]) 
+(0.5*v1R^2*v2I*SLHA::HMIX[35.]) 
+(0.5*v1I^2*v2R*SLHA::HMIX[35.]) 
+(v1I*v1R*v2R*SLHA::HMIX[35.]) 
+(0.5*v1R^2*v2R*SLHA::HMIX[35.]) 
+(0.5*v1I^3*SLHA::HMIX[36.]) 
+(1.5*v1I^2*v1R*SLHA::HMIX[36.]) 
+(1.5*v1I*v1R^2*SLHA::HMIX[36.]) 
+(0.5*v1R^3*SLHA::HMIX[36.]) 
+(1.5*v1I*v2I^2*SLHA::HMIX[37.]) 
+(1.5*v1R*v2I^2*SLHA::HMIX[37.]) 
+(3.*v1I*v2I*v2R*SLHA::HMIX[37.]) 
+(3.*v1R*v2I*v2R*SLHA::HMIX[37.]) 
+(1.5*v1I*v2R^2*SLHA::HMIX[37.]) 
+(1.5*v1R*v2R^2*SLHA::HMIX[37.]) 
; 
} 
</tadpoles> 
 
 
<mass-squared_matrix 
particle="hh"  rotationmatrix="ZH"  spin="scalar" factor="1" >     
(SLHA::HMIX[20.]+3.*v1I^2*SLHA::HMIX[31.]+6.*v1I*v1R*SLHA::HMIX[31.]+3.*v1R^2*SLHA::HMIX[31.]+0.5*v2I^2*SLHA::HMIX[33.]+v2I*v2R*SLHA::HMIX[33.]+0.5*v2R^2*SLHA::HMIX[33.]+0.5*v2I^2*SLHA::HMIX[34.]+v2I*v2R*SLHA::HMIX[34.]+0.5*v2R^2*SLHA::HMIX[34.]+0.5*v2I^2*SLHA::HMIX[35.]+v2I*v2R*SLHA::HMIX[35.]+0.5*v2R^2*SLHA::HMIX[35.]+3.*v1I*v2I*SLHA::HMIX[36.]+3.*v1R*v2I*SLHA::HMIX[36.]+3.*v1I*v2R*SLHA::HMIX[36.]+3.*v1R*v2R*SLHA::HMIX[36.]), 
(-1.*SLHA::HMIX[22.]+v1I*v2I*SLHA::HMIX[33.]+v1R*v2I*SLHA::HMIX[33.]+v1I*v2R*SLHA::HMIX[33.]+v1R*v2R*SLHA::HMIX[33.]+v1I*v2I*SLHA::HMIX[34.]+v1R*v2I*SLHA::HMIX[34.]+v1I*v2R*SLHA::HMIX[34.]+v1R*v2R*SLHA::HMIX[34.]+v1I*v2I*SLHA::HMIX[35.]+v1R*v2I*SLHA::HMIX[35.]+v1I*v2R*SLHA::HMIX[35.]+v1R*v2R*SLHA::HMIX[35.]+1.5*v1I^2*SLHA::HMIX[36.]+3.*v1I*v1R*SLHA::HMIX[36.]+1.5*v1R^2*SLHA::HMIX[36.]+1.5*v2I^2*SLHA::HMIX[37.]+3.*v2I*v2R*SLHA::HMIX[37.]+1.5*v2R^2*SLHA::HMIX[37.]), 
(0.), 
(0.), 
(-1.*SLHA::HMIX[22.]+v1I*v2I*SLHA::HMIX[33.]+v1R*v2I*SLHA::HMIX[33.]+v1I*v2R*SLHA::HMIX[33.]+v1R*v2R*SLHA::HMIX[33.]+v1I*v2I*SLHA::HMIX[34.]+v1R*v2I*SLHA::HMIX[34.]+v1I*v2R*SLHA::HMIX[34.]+v1R*v2R*SLHA::HMIX[34.]+v1I*v2I*SLHA::HMIX[35.]+v1R*v2I*SLHA::HMIX[35.]+v1I*v2R*SLHA::HMIX[35.]+v1R*v2R*SLHA::HMIX[35.]+1.5*v1I^2*SLHA::HMIX[36.]+3.*v1I*v1R*SLHA::HMIX[36.]+1.5*v1R^2*SLHA::HMIX[36.]+1.5*v2I^2*SLHA::HMIX[37.]+3.*v2I*v2R*SLHA::HMIX[37.]+1.5*v2R^2*SLHA::HMIX[37.]), 
(SLHA::HMIX[21.]+3.*v2I^2*SLHA::HMIX[32.]+6.*v2I*v2R*SLHA::HMIX[32.]+3.*v2R^2*SLHA::HMIX[32.]+0.5*v1I^2*SLHA::HMIX[33.]+v1I*v1R*SLHA::HMIX[33.]+0.5*v1R^2*SLHA::HMIX[33.]+0.5*v1I^2*SLHA::HMIX[34.]+v1I*v1R*SLHA::HMIX[34.]+0.5*v1R^2*SLHA::HMIX[34.]+0.5*v1I^2*SLHA::HMIX[35.]+v1I*v1R*SLHA::HMIX[35.]+0.5*v1R^2*SLHA::HMIX[35.]+3.*v1I*v2I*SLHA::HMIX[37.]+3.*v1R*v2I*SLHA::HMIX[37.]+3.*v1I*v2R*SLHA::HMIX[37.]+3.*v1R*v2R*SLHA::HMIX[37.]), 
(0.), 
(0.), 
(0.), 
(0.), 
(SLHA::HMIX[20.]+v1I^2*SLHA::HMIX[31.]+2.*v1I*v1R*SLHA::HMIX[31.]+v1R^2*SLHA::HMIX[31.]+0.5*v2I^2*SLHA::HMIX[33.]+v2I*v2R*SLHA::HMIX[33.]+0.5*v2R^2*SLHA::HMIX[33.]+0.5*v2I^2*SLHA::HMIX[34.]+v2I*v2R*SLHA::HMIX[34.]+0.5*v2R^2*SLHA::HMIX[34.]-0.5*v2I^2*SLHA::HMIX[35.]-1.*v2I*v2R*SLHA::HMIX[35.]-0.5*v2R^2*SLHA::HMIX[35.]+v1I*v2I*SLHA::HMIX[36.]+v1R*v2I*SLHA::HMIX[36.]+v1I*v2R*SLHA::HMIX[36.]+v1R*v2R*SLHA::HMIX[36.]), 
(-1.*SLHA::HMIX[22.]+v1I*v2I*SLHA::HMIX[35.]+v1R*v2I*SLHA::HMIX[35.]+v1I*v2R*SLHA::HMIX[35.]+v1R*v2R*SLHA::HMIX[35.]+0.5*v1I^2*SLHA::HMIX[36.]+v1I*v1R*SLHA::HMIX[36.]+0.5*v1R^2*SLHA::HMIX[36.]+0.5*v2I^2*SLHA::HMIX[37.]+v2I*v2R*SLHA::HMIX[37.]+0.5*v2R^2*SLHA::HMIX[37.]), 
(0.), 
(0.), 
(-1.*SLHA::HMIX[22.]+v1I*v2I*SLHA::HMIX[35.]+v1R*v2I*SLHA::HMIX[35.]+v1I*v2R*SLHA::HMIX[35.]+v1R*v2R*SLHA::HMIX[35.]+0.5*v1I^2*SLHA::HMIX[36.]+v1I*v1R*SLHA::HMIX[36.]+0.5*v1R^2*SLHA::HMIX[36.]+0.5*v2I^2*SLHA::HMIX[37.]+v2I*v2R*SLHA::HMIX[37.]+0.5*v2R^2*SLHA::HMIX[37.]), 
(SLHA::HMIX[21.]+v2I^2*SLHA::HMIX[32.]+2.*v2I*v2R*SLHA::HMIX[32.]+v2R^2*SLHA::HMIX[32.]+0.5*v1I^2*SLHA::HMIX[33.]+v1I*v1R*SLHA::HMIX[33.]+0.5*v1R^2*SLHA::HMIX[33.]+0.5*v1I^2*SLHA::HMIX[34.]+v1I*v1R*SLHA::HMIX[34.]+0.5*v1R^2*SLHA::HMIX[34.]-0.5*v1I^2*SLHA::HMIX[35.]-1.*v1I*v1R*SLHA::HMIX[35.]-0.5*v1R^2*SLHA::HMIX[35.]+v1I*v2I*SLHA::HMIX[37.]+v1R*v2I*SLHA::HMIX[37.]+v1I*v2R*SLHA::HMIX[37.]+v1R*v2R*SLHA::HMIX[37.]) 
</mass-squared_matrix> 
 
<mass-squared_matrix 
particle="Hm"  rotationmatrix="ZP"  spin="scalar" factor="2" >     
(SLHA::HMIX[20.]+v1I^2*SLHA::HMIX[31.]+2.*v1I*v1R*SLHA::HMIX[31.]+v1R^2*SLHA::HMIX[31.]+0.5*v2I^2*SLHA::HMIX[33.]+v2I*v2R*SLHA::HMIX[33.]+0.5*v2R^2*SLHA::HMIX[33.]+v1I*v2I*SLHA::HMIX[36.]+v1R*v2I*SLHA::HMIX[36.]+v1I*v2R*SLHA::HMIX[36.]+v1R*v2R*SLHA::HMIX[36.]), 
(-1.*SLHA::HMIX[22.]+0.5*v1I*v2I*SLHA::HMIX[34.]+0.5*v1R*v2I*SLHA::HMIX[34.]+0.5*v1I*v2R*SLHA::HMIX[34.]+0.5*v1R*v2R*SLHA::HMIX[34.]+0.5*v1I*v2I*SLHA::HMIX[35.]+0.5*v1R*v2I*SLHA::HMIX[35.]+0.5*v1I*v2R*SLHA::HMIX[35.]+0.5*v1R*v2R*SLHA::HMIX[35.]+0.5*v1I^2*SLHA::HMIX[36.]+v1I*v1R*SLHA::HMIX[36.]+0.5*v1R^2*SLHA::HMIX[36.]+0.5*v2I^2*SLHA::HMIX[37.]+v2I*v2R*SLHA::HMIX[37.]+0.5*v2R^2*SLHA::HMIX[37.]), 
(-1.*SLHA::HMIX[22.]+0.5*v1I*v2I*SLHA::HMIX[34.]+0.5*v1R*v2I*SLHA::HMIX[34.]+0.5*v1I*v2R*SLHA::HMIX[34.]+0.5*v1R*v2R*SLHA::HMIX[34.]+0.5*v1I*v2I*SLHA::HMIX[35.]+0.5*v1R*v2I*SLHA::HMIX[35.]+0.5*v1I*v2R*SLHA::HMIX[35.]+0.5*v1R*v2R*SLHA::HMIX[35.]+0.5*v1I^2*SLHA::HMIX[36.]+v1I*v1R*SLHA::HMIX[36.]+0.5*v1R^2*SLHA::HMIX[36.]+0.5*v2I^2*SLHA::HMIX[37.]+v2I*v2R*SLHA::HMIX[37.]+0.5*v2R^2*SLHA::HMIX[37.]), 
(SLHA::HMIX[21.]+v2I^2*SLHA::HMIX[32.]+2.*v2I*v2R*SLHA::HMIX[32.]+v2R^2*SLHA::HMIX[32.]+0.5*v1I^2*SLHA::HMIX[33.]+v1I*v1R*SLHA::HMIX[33.]+0.5*v1R^2*SLHA::HMIX[33.]+v1I*v2I*SLHA::HMIX[37.]+v1R*v2I*SLHA::HMIX[37.]+v1I*v2R*SLHA::HMIX[37.]+v1R*v2R*SLHA::HMIX[37.]) 
</mass-squared_matrix> 
 
<mass-squared_matrix 
particle="Fd"  rotationmatrix1="Vd"  rotationmatrix2="Ud"  spin="Weyl" factor="-12" >     
(0.5*v2I^2*SLHA::EPSD[1.,1.]^2+v2I*v2R*SLHA::EPSD[1.,1.]^2+0.5*v2R^2*SLHA::EPSD[1.,1.]^2+0.5*v2I^2*SLHA::EPSD[2.,1.]^2+v2I*v2R*SLHA::EPSD[2.,1.]^2+0.5*v2R^2*SLHA::EPSD[2.,1.]^2+0.5*v2I^2*SLHA::EPSD[3.,1.]^2+v2I*v2R*SLHA::EPSD[3.,1.]^2+0.5*v2R^2*SLHA::EPSD[3.,1.]^2+v1I*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,1.]+v1R*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,1.]+v1I*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,1.]+v1R*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,1.]+0.5*v1I^2*SLHA::YD[1.,1.]^2+v1I*v1R*SLHA::YD[1.,1.]^2+0.5*v1R^2*SLHA::YD[1.,1.]^2+v1I*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,1.]+v1R*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,1.]+v1I*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,1.]+v1R*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,1.]+0.5*v1I^2*SLHA::YD[2.,1.]^2+v1I*v1R*SLHA::YD[2.,1.]^2+0.5*v1R^2*SLHA::YD[2.,1.]^2+v1I*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,1.]+v1R*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,1.]+v1I*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,1.]+v1R*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,1.]+0.5*v1I^2*SLHA::YD[3.,1.]^2+v1I*v1R*SLHA::YD[3.,1.]^2+0.5*v1R^2*SLHA::YD[3.,1.]^2), 
(0.5*v2I^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,2.]+v2I*v2R*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,2.]+0.5*v2R^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,2.]+0.5*v2I^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,2.]+v2I*v2R*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,2.]+0.5*v2R^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,2.]+0.5*v2I^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,2.]+v2I*v2R*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,2.]+0.5*v2R^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,2.]+0.5*v1I*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1R*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1I*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1R*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1I*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1R*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1I*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1R*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1I^2*SLHA::YD[1.,1.]*SLHA::YD[1.,2.]+v1I*v1R*SLHA::YD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1R^2*SLHA::YD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1I*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1R*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1I*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1R*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1I*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1R*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1I*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1R*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1I^2*SLHA::YD[2.,1.]*SLHA::YD[2.,2.]+v1I*v1R*SLHA::YD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1R^2*SLHA::YD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1I*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1R*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1I*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1R*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1I*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1R*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1I*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1R*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1I^2*SLHA::YD[3.,1.]*SLHA::YD[3.,2.]+v1I*v1R*SLHA::YD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1R^2*SLHA::YD[3.,1.]*SLHA::YD[3.,2.]), 
(0.5*v2I^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,3.]+v2I*v2R*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,3.]+0.5*v2R^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,3.]+0.5*v2I^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,3.]+v2I*v2R*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,3.]+0.5*v2R^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,3.]+0.5*v2I^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,3.]+v2I*v2R*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,3.]+0.5*v2R^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,3.]+0.5*v1I*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1R*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1I*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1R*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1I*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1R*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1I*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1R*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1I^2*SLHA::YD[1.,1.]*SLHA::YD[1.,3.]+v1I*v1R*SLHA::YD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1R^2*SLHA::YD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1I*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1R*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1I*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1R*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1I*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1R*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1I*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1R*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1I^2*SLHA::YD[2.,1.]*SLHA::YD[2.,3.]+v1I*v1R*SLHA::YD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1R^2*SLHA::YD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1I*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1R*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1I*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1R*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1I*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1R*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1I*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1R*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1I^2*SLHA::YD[3.,1.]*SLHA::YD[3.,3.]+v1I*v1R*SLHA::YD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1R^2*SLHA::YD[3.,1.]*SLHA::YD[3.,3.]), 
(0.5*v2I^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,2.]+v2I*v2R*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,2.]+0.5*v2R^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,2.]+0.5*v2I^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,2.]+v2I*v2R*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,2.]+0.5*v2R^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,2.]+0.5*v2I^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,2.]+v2I*v2R*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,2.]+0.5*v2R^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,2.]+0.5*v1I*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1R*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1I*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1R*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,1.]+0.5*v1I*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1R*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1I*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1R*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1I^2*SLHA::YD[1.,1.]*SLHA::YD[1.,2.]+v1I*v1R*SLHA::YD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1R^2*SLHA::YD[1.,1.]*SLHA::YD[1.,2.]+0.5*v1I*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1R*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1I*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1R*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,1.]+0.5*v1I*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1R*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1I*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1R*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1I^2*SLHA::YD[2.,1.]*SLHA::YD[2.,2.]+v1I*v1R*SLHA::YD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1R^2*SLHA::YD[2.,1.]*SLHA::YD[2.,2.]+0.5*v1I*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1R*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1I*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1R*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,1.]+0.5*v1I*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1R*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1I*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1R*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1I^2*SLHA::YD[3.,1.]*SLHA::YD[3.,2.]+v1I*v1R*SLHA::YD[3.,1.]*SLHA::YD[3.,2.]+0.5*v1R^2*SLHA::YD[3.,1.]*SLHA::YD[3.,2.]), 
(0.5*v2I^2*SLHA::EPSD[1.,2.]^2+v2I*v2R*SLHA::EPSD[1.,2.]^2+0.5*v2R^2*SLHA::EPSD[1.,2.]^2+0.5*v2I^2*SLHA::EPSD[2.,2.]^2+v2I*v2R*SLHA::EPSD[2.,2.]^2+0.5*v2R^2*SLHA::EPSD[2.,2.]^2+0.5*v2I^2*SLHA::EPSD[3.,2.]^2+v2I*v2R*SLHA::EPSD[3.,2.]^2+0.5*v2R^2*SLHA::EPSD[3.,2.]^2+v1I*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,2.]+v1R*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,2.]+v1I*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,2.]+v1R*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,2.]+0.5*v1I^2*SLHA::YD[1.,2.]^2+v1I*v1R*SLHA::YD[1.,2.]^2+0.5*v1R^2*SLHA::YD[1.,2.]^2+v1I*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,2.]+v1R*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,2.]+v1I*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,2.]+v1R*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,2.]+0.5*v1I^2*SLHA::YD[2.,2.]^2+v1I*v1R*SLHA::YD[2.,2.]^2+0.5*v1R^2*SLHA::YD[2.,2.]^2+v1I*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,2.]+v1R*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,2.]+v1I*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,2.]+v1R*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,2.]+0.5*v1I^2*SLHA::YD[3.,2.]^2+v1I*v1R*SLHA::YD[3.,2.]^2+0.5*v1R^2*SLHA::YD[3.,2.]^2), 
(0.5*v2I^2*SLHA::EPSD[1.,2.]*SLHA::EPSD[1.,3.]+v2I*v2R*SLHA::EPSD[1.,2.]*SLHA::EPSD[1.,3.]+0.5*v2R^2*SLHA::EPSD[1.,2.]*SLHA::EPSD[1.,3.]+0.5*v2I^2*SLHA::EPSD[2.,2.]*SLHA::EPSD[2.,3.]+v2I*v2R*SLHA::EPSD[2.,2.]*SLHA::EPSD[2.,3.]+0.5*v2R^2*SLHA::EPSD[2.,2.]*SLHA::EPSD[2.,3.]+0.5*v2I^2*SLHA::EPSD[3.,2.]*SLHA::EPSD[3.,3.]+v2I*v2R*SLHA::EPSD[3.,2.]*SLHA::EPSD[3.,3.]+0.5*v2R^2*SLHA::EPSD[3.,2.]*SLHA::EPSD[3.,3.]+0.5*v1I*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1R*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1I*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1R*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1I*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1R*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1I*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1R*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1I^2*SLHA::YD[1.,2.]*SLHA::YD[1.,3.]+v1I*v1R*SLHA::YD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1R^2*SLHA::YD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1I*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1R*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1I*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1R*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1I*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1R*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1I*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1R*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1I^2*SLHA::YD[2.,2.]*SLHA::YD[2.,3.]+v1I*v1R*SLHA::YD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1R^2*SLHA::YD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1I*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1R*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1I*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1R*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1I*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1R*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1I*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1R*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1I^2*SLHA::YD[3.,2.]*SLHA::YD[3.,3.]+v1I*v1R*SLHA::YD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1R^2*SLHA::YD[3.,2.]*SLHA::YD[3.,3.]), 
(0.5*v2I^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,3.]+v2I*v2R*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,3.]+0.5*v2R^2*SLHA::EPSD[1.,1.]*SLHA::EPSD[1.,3.]+0.5*v2I^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,3.]+v2I*v2R*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,3.]+0.5*v2R^2*SLHA::EPSD[2.,1.]*SLHA::EPSD[2.,3.]+0.5*v2I^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,3.]+v2I*v2R*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,3.]+0.5*v2R^2*SLHA::EPSD[3.,1.]*SLHA::EPSD[3.,3.]+0.5*v1I*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1R*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1I*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1R*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,1.]+0.5*v1I*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1R*v2I*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1I*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1R*v2R*SLHA::EPSD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1I^2*SLHA::YD[1.,1.]*SLHA::YD[1.,3.]+v1I*v1R*SLHA::YD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1R^2*SLHA::YD[1.,1.]*SLHA::YD[1.,3.]+0.5*v1I*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1R*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1I*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1R*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,1.]+0.5*v1I*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1R*v2I*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1I*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1R*v2R*SLHA::EPSD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1I^2*SLHA::YD[2.,1.]*SLHA::YD[2.,3.]+v1I*v1R*SLHA::YD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1R^2*SLHA::YD[2.,1.]*SLHA::YD[2.,3.]+0.5*v1I*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1R*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1I*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1R*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,1.]+0.5*v1I*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1R*v2I*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1I*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1R*v2R*SLHA::EPSD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1I^2*SLHA::YD[3.,1.]*SLHA::YD[3.,3.]+v1I*v1R*SLHA::YD[3.,1.]*SLHA::YD[3.,3.]+0.5*v1R^2*SLHA::YD[3.,1.]*SLHA::YD[3.,3.]), 
(0.5*v2I^2*SLHA::EPSD[1.,2.]*SLHA::EPSD[1.,3.]+v2I*v2R*SLHA::EPSD[1.,2.]*SLHA::EPSD[1.,3.]+0.5*v2R^2*SLHA::EPSD[1.,2.]*SLHA::EPSD[1.,3.]+0.5*v2I^2*SLHA::EPSD[2.,2.]*SLHA::EPSD[2.,3.]+v2I*v2R*SLHA::EPSD[2.,2.]*SLHA::EPSD[2.,3.]+0.5*v2R^2*SLHA::EPSD[2.,2.]*SLHA::EPSD[2.,3.]+0.5*v2I^2*SLHA::EPSD[3.,2.]*SLHA::EPSD[3.,3.]+v2I*v2R*SLHA::EPSD[3.,2.]*SLHA::EPSD[3.,3.]+0.5*v2R^2*SLHA::EPSD[3.,2.]*SLHA::EPSD[3.,3.]+0.5*v1I*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1R*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1I*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1R*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,2.]+0.5*v1I*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1R*v2I*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1I*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1R*v2R*SLHA::EPSD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1I^2*SLHA::YD[1.,2.]*SLHA::YD[1.,3.]+v1I*v1R*SLHA::YD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1R^2*SLHA::YD[1.,2.]*SLHA::YD[1.,3.]+0.5*v1I*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1R*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1I*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1R*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,2.]+0.5*v1I*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1R*v2I*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1I*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1R*v2R*SLHA::EPSD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1I^2*SLHA::YD[2.,2.]*SLHA::YD[2.,3.]+v1I*v1R*SLHA::YD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1R^2*SLHA::YD[2.,2.]*SLHA::YD[2.,3.]+0.5*v1I*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1R*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1I*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1R*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,2.]+0.5*v1I*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1R*v2I*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1I*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1R*v2R*SLHA::EPSD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1I^2*SLHA::YD[3.,2.]*SLHA::YD[3.,3.]+v1I*v1R*SLHA::YD[3.,2.]*SLHA::YD[3.,3.]+0.5*v1R^2*SLHA::YD[3.,2.]*SLHA::YD[3.,3.]), 
(0.5*v2I^2*SLHA::EPSD[1.,3.]^2+v2I*v2R*SLHA::EPSD[1.,3.]^2+0.5*v2R^2*SLHA::EPSD[1.,3.]^2+0.5*v2I^2*SLHA::EPSD[2.,3.]^2+v2I*v2R*SLHA::EPSD[2.,3.]^2+0.5*v2R^2*SLHA::EPSD[2.,3.]^2+0.5*v2I^2*SLHA::EPSD[3.,3.]^2+v2I*v2R*SLHA::EPSD[3.,3.]^2+0.5*v2R^2*SLHA::EPSD[3.,3.]^2+v1I*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,3.]+v1R*v2I*SLHA::EPSD[1.,3.]*SLHA::YD[1.,3.]+v1I*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,3.]+v1R*v2R*SLHA::EPSD[1.,3.]*SLHA::YD[1.,3.]+0.5*v1I^2*SLHA::YD[1.,3.]^2+v1I*v1R*SLHA::YD[1.,3.]^2+0.5*v1R^2*SLHA::YD[1.,3.]^2+v1I*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,3.]+v1R*v2I*SLHA::EPSD[2.,3.]*SLHA::YD[2.,3.]+v1I*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,3.]+v1R*v2R*SLHA::EPSD[2.,3.]*SLHA::YD[2.,3.]+0.5*v1I^2*SLHA::YD[2.,3.]^2+v1I*v1R*SLHA::YD[2.,3.]^2+0.5*v1R^2*SLHA::YD[2.,3.]^2+v1I*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,3.]+v1R*v2I*SLHA::EPSD[3.,3.]*SLHA::YD[3.,3.]+v1I*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,3.]+v1R*v2R*SLHA::EPSD[3.,3.]*SLHA::YD[3.,3.]+0.5*v1I^2*SLHA::YD[3.,3.]^2+v1I*v1R*SLHA::YD[3.,3.]^2+0.5*v1R^2*SLHA::YD[3.,3.]^2) 
</mass-squared_matrix> 
 
<mass-squared_matrix 
particle="Fu"  rotationmatrix1="Vu"  rotationmatrix2="Uu"  spin="Weyl" factor="-12" >     
(0.5*v1I^2*SLHA::EPSU[1.,1.]^2+v1I*v1R*SLHA::EPSU[1.,1.]^2+0.5*v1R^2*SLHA::EPSU[1.,1.]^2+0.5*v1I^2*SLHA::EPSU[2.,1.]^2+v1I*v1R*SLHA::EPSU[2.,1.]^2+0.5*v1R^2*SLHA::EPSU[2.,1.]^2+0.5*v1I^2*SLHA::EPSU[3.,1.]^2+v1I*v1R*SLHA::EPSU[3.,1.]^2+0.5*v1R^2*SLHA::EPSU[3.,1.]^2+v1I*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,1.]+v1R*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,1.]+v1I*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,1.]+v1R*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,1.]+0.5*v2I^2*SLHA::YU[1.,1.]^2+v2I*v2R*SLHA::YU[1.,1.]^2+0.5*v2R^2*SLHA::YU[1.,1.]^2+v1I*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,1.]+v1R*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,1.]+v1I*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,1.]+v1R*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,1.]+0.5*v2I^2*SLHA::YU[2.,1.]^2+v2I*v2R*SLHA::YU[2.,1.]^2+0.5*v2R^2*SLHA::YU[2.,1.]^2+v1I*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,1.]+v1R*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,1.]+v1I*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,1.]+v1R*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,1.]+0.5*v2I^2*SLHA::YU[3.,1.]^2+v2I*v2R*SLHA::YU[3.,1.]^2+0.5*v2R^2*SLHA::YU[3.,1.]^2), 
(0.5*v1I^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,2.]+v1I*v1R*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,2.]+0.5*v1R^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,2.]+0.5*v1I^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,2.]+v1I*v1R*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,2.]+0.5*v1R^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,2.]+0.5*v1I^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,2.]+v1I*v1R*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,2.]+0.5*v1R^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,2.]+0.5*v1I*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1R*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1I*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1R*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1I*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1R*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1I*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1R*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v2I^2*SLHA::YU[1.,1.]*SLHA::YU[1.,2.]+v2I*v2R*SLHA::YU[1.,1.]*SLHA::YU[1.,2.]+0.5*v2R^2*SLHA::YU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1I*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1R*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1I*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1R*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1I*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1R*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1I*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1R*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v2I^2*SLHA::YU[2.,1.]*SLHA::YU[2.,2.]+v2I*v2R*SLHA::YU[2.,1.]*SLHA::YU[2.,2.]+0.5*v2R^2*SLHA::YU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1I*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1R*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1I*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1R*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1I*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v1R*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v1I*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v1R*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v2I^2*SLHA::YU[3.,1.]*SLHA::YU[3.,2.]+v2I*v2R*SLHA::YU[3.,1.]*SLHA::YU[3.,2.]+0.5*v2R^2*SLHA::YU[3.,1.]*SLHA::YU[3.,2.]), 
(0.5*v1I^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,3.]+v1I*v1R*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,3.]+0.5*v1R^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,3.]+0.5*v1I^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,3.]+v1I*v1R*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,3.]+0.5*v1R^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,3.]+0.5*v1I^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,3.]+v1I*v1R*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,3.]+0.5*v1R^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,3.]+0.5*v1I*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1R*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1I*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1R*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1I*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1R*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1I*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1R*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v2I^2*SLHA::YU[1.,1.]*SLHA::YU[1.,3.]+v2I*v2R*SLHA::YU[1.,1.]*SLHA::YU[1.,3.]+0.5*v2R^2*SLHA::YU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1I*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1R*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1I*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1R*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1I*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1R*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1I*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1R*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v2I^2*SLHA::YU[2.,1.]*SLHA::YU[2.,3.]+v2I*v2R*SLHA::YU[2.,1.]*SLHA::YU[2.,3.]+0.5*v2R^2*SLHA::YU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1I*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1R*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1I*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1R*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1I*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v1R*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v1I*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v1R*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v2I^2*SLHA::YU[3.,1.]*SLHA::YU[3.,3.]+v2I*v2R*SLHA::YU[3.,1.]*SLHA::YU[3.,3.]+0.5*v2R^2*SLHA::YU[3.,1.]*SLHA::YU[3.,3.]), 
(0.5*v1I^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,2.]+v1I*v1R*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,2.]+0.5*v1R^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,2.]+0.5*v1I^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,2.]+v1I*v1R*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,2.]+0.5*v1R^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,2.]+0.5*v1I^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,2.]+v1I*v1R*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,2.]+0.5*v1R^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,2.]+0.5*v1I*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1R*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1I*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1R*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,1.]+0.5*v1I*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1R*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1I*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1R*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,2.]+0.5*v2I^2*SLHA::YU[1.,1.]*SLHA::YU[1.,2.]+v2I*v2R*SLHA::YU[1.,1.]*SLHA::YU[1.,2.]+0.5*v2R^2*SLHA::YU[1.,1.]*SLHA::YU[1.,2.]+0.5*v1I*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1R*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1I*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1R*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,1.]+0.5*v1I*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1R*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1I*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1R*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,2.]+0.5*v2I^2*SLHA::YU[2.,1.]*SLHA::YU[2.,2.]+v2I*v2R*SLHA::YU[2.,1.]*SLHA::YU[2.,2.]+0.5*v2R^2*SLHA::YU[2.,1.]*SLHA::YU[2.,2.]+0.5*v1I*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1R*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1I*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1R*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,1.]+0.5*v1I*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v1R*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v1I*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v1R*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,2.]+0.5*v2I^2*SLHA::YU[3.,1.]*SLHA::YU[3.,2.]+v2I*v2R*SLHA::YU[3.,1.]*SLHA::YU[3.,2.]+0.5*v2R^2*SLHA::YU[3.,1.]*SLHA::YU[3.,2.]), 
(0.5*v1I^2*SLHA::EPSU[1.,2.]^2+v1I*v1R*SLHA::EPSU[1.,2.]^2+0.5*v1R^2*SLHA::EPSU[1.,2.]^2+0.5*v1I^2*SLHA::EPSU[2.,2.]^2+v1I*v1R*SLHA::EPSU[2.,2.]^2+0.5*v1R^2*SLHA::EPSU[2.,2.]^2+0.5*v1I^2*SLHA::EPSU[3.,2.]^2+v1I*v1R*SLHA::EPSU[3.,2.]^2+0.5*v1R^2*SLHA::EPSU[3.,2.]^2+v1I*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,2.]+v1R*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,2.]+v1I*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,2.]+v1R*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,2.]+0.5*v2I^2*SLHA::YU[1.,2.]^2+v2I*v2R*SLHA::YU[1.,2.]^2+0.5*v2R^2*SLHA::YU[1.,2.]^2+v1I*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,2.]+v1R*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,2.]+v1I*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,2.]+v1R*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,2.]+0.5*v2I^2*SLHA::YU[2.,2.]^2+v2I*v2R*SLHA::YU[2.,2.]^2+0.5*v2R^2*SLHA::YU[2.,2.]^2+v1I*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,2.]+v1R*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,2.]+v1I*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,2.]+v1R*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,2.]+0.5*v2I^2*SLHA::YU[3.,2.]^2+v2I*v2R*SLHA::YU[3.,2.]^2+0.5*v2R^2*SLHA::YU[3.,2.]^2), 
(0.5*v1I^2*SLHA::EPSU[1.,2.]*SLHA::EPSU[1.,3.]+v1I*v1R*SLHA::EPSU[1.,2.]*SLHA::EPSU[1.,3.]+0.5*v1R^2*SLHA::EPSU[1.,2.]*SLHA::EPSU[1.,3.]+0.5*v1I^2*SLHA::EPSU[2.,2.]*SLHA::EPSU[2.,3.]+v1I*v1R*SLHA::EPSU[2.,2.]*SLHA::EPSU[2.,3.]+0.5*v1R^2*SLHA::EPSU[2.,2.]*SLHA::EPSU[2.,3.]+0.5*v1I^2*SLHA::EPSU[3.,2.]*SLHA::EPSU[3.,3.]+v1I*v1R*SLHA::EPSU[3.,2.]*SLHA::EPSU[3.,3.]+0.5*v1R^2*SLHA::EPSU[3.,2.]*SLHA::EPSU[3.,3.]+0.5*v1I*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1R*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1I*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1R*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1I*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1R*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1I*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1R*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v2I^2*SLHA::YU[1.,2.]*SLHA::YU[1.,3.]+v2I*v2R*SLHA::YU[1.,2.]*SLHA::YU[1.,3.]+0.5*v2R^2*SLHA::YU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1I*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1R*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1I*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1R*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1I*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1R*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1I*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1R*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v2I^2*SLHA::YU[2.,2.]*SLHA::YU[2.,3.]+v2I*v2R*SLHA::YU[2.,2.]*SLHA::YU[2.,3.]+0.5*v2R^2*SLHA::YU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1I*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1R*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1I*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1R*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1I*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v1R*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v1I*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v1R*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v2I^2*SLHA::YU[3.,2.]*SLHA::YU[3.,3.]+v2I*v2R*SLHA::YU[3.,2.]*SLHA::YU[3.,3.]+0.5*v2R^2*SLHA::YU[3.,2.]*SLHA::YU[3.,3.]), 
(0.5*v1I^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,3.]+v1I*v1R*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,3.]+0.5*v1R^2*SLHA::EPSU[1.,1.]*SLHA::EPSU[1.,3.]+0.5*v1I^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,3.]+v1I*v1R*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,3.]+0.5*v1R^2*SLHA::EPSU[2.,1.]*SLHA::EPSU[2.,3.]+0.5*v1I^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,3.]+v1I*v1R*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,3.]+0.5*v1R^2*SLHA::EPSU[3.,1.]*SLHA::EPSU[3.,3.]+0.5*v1I*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1R*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1I*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1R*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,1.]+0.5*v1I*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1R*v2I*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1I*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1R*v2R*SLHA::EPSU[1.,1.]*SLHA::YU[1.,3.]+0.5*v2I^2*SLHA::YU[1.,1.]*SLHA::YU[1.,3.]+v2I*v2R*SLHA::YU[1.,1.]*SLHA::YU[1.,3.]+0.5*v2R^2*SLHA::YU[1.,1.]*SLHA::YU[1.,3.]+0.5*v1I*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1R*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1I*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1R*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,1.]+0.5*v1I*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1R*v2I*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1I*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1R*v2R*SLHA::EPSU[2.,1.]*SLHA::YU[2.,3.]+0.5*v2I^2*SLHA::YU[2.,1.]*SLHA::YU[2.,3.]+v2I*v2R*SLHA::YU[2.,1.]*SLHA::YU[2.,3.]+0.5*v2R^2*SLHA::YU[2.,1.]*SLHA::YU[2.,3.]+0.5*v1I*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1R*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1I*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1R*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,1.]+0.5*v1I*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v1R*v2I*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v1I*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v1R*v2R*SLHA::EPSU[3.,1.]*SLHA::YU[3.,3.]+0.5*v2I^2*SLHA::YU[3.,1.]*SLHA::YU[3.,3.]+v2I*v2R*SLHA::YU[3.,1.]*SLHA::YU[3.,3.]+0.5*v2R^2*SLHA::YU[3.,1.]*SLHA::YU[3.,3.]), 
(0.5*v1I^2*SLHA::EPSU[1.,2.]*SLHA::EPSU[1.,3.]+v1I*v1R*SLHA::EPSU[1.,2.]*SLHA::EPSU[1.,3.]+0.5*v1R^2*SLHA::EPSU[1.,2.]*SLHA::EPSU[1.,3.]+0.5*v1I^2*SLHA::EPSU[2.,2.]*SLHA::EPSU[2.,3.]+v1I*v1R*SLHA::EPSU[2.,2.]*SLHA::EPSU[2.,3.]+0.5*v1R^2*SLHA::EPSU[2.,2.]*SLHA::EPSU[2.,3.]+0.5*v1I^2*SLHA::EPSU[3.,2.]*SLHA::EPSU[3.,3.]+v1I*v1R*SLHA::EPSU[3.,2.]*SLHA::EPSU[3.,3.]+0.5*v1R^2*SLHA::EPSU[3.,2.]*SLHA::EPSU[3.,3.]+0.5*v1I*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1R*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1I*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1R*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,2.]+0.5*v1I*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1R*v2I*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1I*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1R*v2R*SLHA::EPSU[1.,2.]*SLHA::YU[1.,3.]+0.5*v2I^2*SLHA::YU[1.,2.]*SLHA::YU[1.,3.]+v2I*v2R*SLHA::YU[1.,2.]*SLHA::YU[1.,3.]+0.5*v2R^2*SLHA::YU[1.,2.]*SLHA::YU[1.,3.]+0.5*v1I*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1R*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1I*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1R*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,2.]+0.5*v1I*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1R*v2I*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1I*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1R*v2R*SLHA::EPSU[2.,2.]*SLHA::YU[2.,3.]+0.5*v2I^2*SLHA::YU[2.,2.]*SLHA::YU[2.,3.]+v2I*v2R*SLHA::YU[2.,2.]*SLHA::YU[2.,3.]+0.5*v2R^2*SLHA::YU[2.,2.]*SLHA::YU[2.,3.]+0.5*v1I*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1R*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1I*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1R*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,2.]+0.5*v1I*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v1R*v2I*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v1I*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v1R*v2R*SLHA::EPSU[3.,2.]*SLHA::YU[3.,3.]+0.5*v2I^2*SLHA::YU[3.,2.]*SLHA::YU[3.,3.]+v2I*v2R*SLHA::YU[3.,2.]*SLHA::YU[3.,3.]+0.5*v2R^2*SLHA::YU[3.,2.]*SLHA::YU[3.,3.]), 
(0.5*v1I^2*SLHA::EPSU[1.,3.]^2+v1I*v1R*SLHA::EPSU[1.,3.]^2+0.5*v1R^2*SLHA::EPSU[1.,3.]^2+0.5*v1I^2*SLHA::EPSU[2.,3.]^2+v1I*v1R*SLHA::EPSU[2.,3.]^2+0.5*v1R^2*SLHA::EPSU[2.,3.]^2+0.5*v1I^2*SLHA::EPSU[3.,3.]^2+v1I*v1R*SLHA::EPSU[3.,3.]^2+0.5*v1R^2*SLHA::EPSU[3.,3.]^2+v1I*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,3.]+v1R*v2I*SLHA::EPSU[1.,3.]*SLHA::YU[1.,3.]+v1I*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,3.]+v1R*v2R*SLHA::EPSU[1.,3.]*SLHA::YU[1.,3.]+0.5*v2I^2*SLHA::YU[1.,3.]^2+v2I*v2R*SLHA::YU[1.,3.]^2+0.5*v2R^2*SLHA::YU[1.,3.]^2+v1I*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,3.]+v1R*v2I*SLHA::EPSU[2.,3.]*SLHA::YU[2.,3.]+v1I*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,3.]+v1R*v2R*SLHA::EPSU[2.,3.]*SLHA::YU[2.,3.]+0.5*v2I^2*SLHA::YU[2.,3.]^2+v2I*v2R*SLHA::YU[2.,3.]^2+0.5*v2R^2*SLHA::YU[2.,3.]^2+v1I*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,3.]+v1R*v2I*SLHA::EPSU[3.,3.]*SLHA::YU[3.,3.]+v1I*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,3.]+v1R*v2R*SLHA::EPSU[3.,3.]*SLHA::YU[3.,3.]+0.5*v2I^2*SLHA::YU[3.,3.]^2+v2I*v2R*SLHA::YU[3.,3.]^2+0.5*v2R^2*SLHA::YU[3.,3.]^2) 
</mass-squared_matrix> 
 
<mass-squared_matrix 
particle="Fe"  rotationmatrix1="Ve"  rotationmatrix2="Ue"  spin="Weyl" factor="-4" >     
(0.5*v2I^2*SLHA::EPSE[1.,1.]^2+v2I*v2R*SLHA::EPSE[1.,1.]^2+0.5*v2R^2*SLHA::EPSE[1.,1.]^2+0.5*v2I^2*SLHA::EPSE[2.,1.]^2+v2I*v2R*SLHA::EPSE[2.,1.]^2+0.5*v2R^2*SLHA::EPSE[2.,1.]^2+0.5*v2I^2*SLHA::EPSE[3.,1.]^2+v2I*v2R*SLHA::EPSE[3.,1.]^2+0.5*v2R^2*SLHA::EPSE[3.,1.]^2+v1I*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,1.]+v1R*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,1.]+v1I*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,1.]+v1R*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,1.]+0.5*v1I^2*SLHA::YE[1.,1.]^2+v1I*v1R*SLHA::YE[1.,1.]^2+0.5*v1R^2*SLHA::YE[1.,1.]^2+v1I*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,1.]+v1R*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,1.]+v1I*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,1.]+v1R*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,1.]+0.5*v1I^2*SLHA::YE[2.,1.]^2+v1I*v1R*SLHA::YE[2.,1.]^2+0.5*v1R^2*SLHA::YE[2.,1.]^2+v1I*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,1.]+v1R*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,1.]+v1I*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,1.]+v1R*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,1.]+0.5*v1I^2*SLHA::YE[3.,1.]^2+v1I*v1R*SLHA::YE[3.,1.]^2+0.5*v1R^2*SLHA::YE[3.,1.]^2), 
(0.5*v2I^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,2.]+v2I*v2R*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,2.]+0.5*v2R^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,2.]+0.5*v2I^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,2.]+v2I*v2R*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,2.]+0.5*v2R^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,2.]+0.5*v2I^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,2.]+v2I*v2R*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,2.]+0.5*v2R^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,2.]+0.5*v1I*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1R*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1I*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1R*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1I*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1R*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1I*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1R*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1I^2*SLHA::YE[1.,1.]*SLHA::YE[1.,2.]+v1I*v1R*SLHA::YE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1R^2*SLHA::YE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1I*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1R*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1I*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1R*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1I*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1R*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1I*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1R*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1I^2*SLHA::YE[2.,1.]*SLHA::YE[2.,2.]+v1I*v1R*SLHA::YE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1R^2*SLHA::YE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1I*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1R*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1I*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1R*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1I*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1R*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1I*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1R*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1I^2*SLHA::YE[3.,1.]*SLHA::YE[3.,2.]+v1I*v1R*SLHA::YE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1R^2*SLHA::YE[3.,1.]*SLHA::YE[3.,2.]), 
(0.5*v2I^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,3.]+v2I*v2R*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,3.]+0.5*v2R^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,3.]+0.5*v2I^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,3.]+v2I*v2R*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,3.]+0.5*v2R^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,3.]+0.5*v2I^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,3.]+v2I*v2R*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,3.]+0.5*v2R^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,3.]+0.5*v1I*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1R*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1I*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1R*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1I*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1R*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1I*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1R*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1I^2*SLHA::YE[1.,1.]*SLHA::YE[1.,3.]+v1I*v1R*SLHA::YE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1R^2*SLHA::YE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1I*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1R*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1I*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1R*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1I*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1R*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1I*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1R*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1I^2*SLHA::YE[2.,1.]*SLHA::YE[2.,3.]+v1I*v1R*SLHA::YE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1R^2*SLHA::YE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1I*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1R*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1I*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1R*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1I*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1R*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1I*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1R*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1I^2*SLHA::YE[3.,1.]*SLHA::YE[3.,3.]+v1I*v1R*SLHA::YE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1R^2*SLHA::YE[3.,1.]*SLHA::YE[3.,3.]), 
(0.5*v2I^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,2.]+v2I*v2R*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,2.]+0.5*v2R^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,2.]+0.5*v2I^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,2.]+v2I*v2R*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,2.]+0.5*v2R^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,2.]+0.5*v2I^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,2.]+v2I*v2R*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,2.]+0.5*v2R^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,2.]+0.5*v1I*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1R*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1I*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1R*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,1.]+0.5*v1I*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1R*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1I*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1R*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1I^2*SLHA::YE[1.,1.]*SLHA::YE[1.,2.]+v1I*v1R*SLHA::YE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1R^2*SLHA::YE[1.,1.]*SLHA::YE[1.,2.]+0.5*v1I*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1R*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1I*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1R*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,1.]+0.5*v1I*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1R*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1I*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1R*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1I^2*SLHA::YE[2.,1.]*SLHA::YE[2.,2.]+v1I*v1R*SLHA::YE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1R^2*SLHA::YE[2.,1.]*SLHA::YE[2.,2.]+0.5*v1I*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1R*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1I*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1R*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,1.]+0.5*v1I*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1R*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1I*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1R*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1I^2*SLHA::YE[3.,1.]*SLHA::YE[3.,2.]+v1I*v1R*SLHA::YE[3.,1.]*SLHA::YE[3.,2.]+0.5*v1R^2*SLHA::YE[3.,1.]*SLHA::YE[3.,2.]), 
(0.5*v2I^2*SLHA::EPSE[1.,2.]^2+v2I*v2R*SLHA::EPSE[1.,2.]^2+0.5*v2R^2*SLHA::EPSE[1.,2.]^2+0.5*v2I^2*SLHA::EPSE[2.,2.]^2+v2I*v2R*SLHA::EPSE[2.,2.]^2+0.5*v2R^2*SLHA::EPSE[2.,2.]^2+0.5*v2I^2*SLHA::EPSE[3.,2.]^2+v2I*v2R*SLHA::EPSE[3.,2.]^2+0.5*v2R^2*SLHA::EPSE[3.,2.]^2+v1I*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,2.]+v1R*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,2.]+v1I*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,2.]+v1R*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,2.]+0.5*v1I^2*SLHA::YE[1.,2.]^2+v1I*v1R*SLHA::YE[1.,2.]^2+0.5*v1R^2*SLHA::YE[1.,2.]^2+v1I*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,2.]+v1R*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,2.]+v1I*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,2.]+v1R*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,2.]+0.5*v1I^2*SLHA::YE[2.,2.]^2+v1I*v1R*SLHA::YE[2.,2.]^2+0.5*v1R^2*SLHA::YE[2.,2.]^2+v1I*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,2.]+v1R*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,2.]+v1I*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,2.]+v1R*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,2.]+0.5*v1I^2*SLHA::YE[3.,2.]^2+v1I*v1R*SLHA::YE[3.,2.]^2+0.5*v1R^2*SLHA::YE[3.,2.]^2), 
(0.5*v2I^2*SLHA::EPSE[1.,2.]*SLHA::EPSE[1.,3.]+v2I*v2R*SLHA::EPSE[1.,2.]*SLHA::EPSE[1.,3.]+0.5*v2R^2*SLHA::EPSE[1.,2.]*SLHA::EPSE[1.,3.]+0.5*v2I^2*SLHA::EPSE[2.,2.]*SLHA::EPSE[2.,3.]+v2I*v2R*SLHA::EPSE[2.,2.]*SLHA::EPSE[2.,3.]+0.5*v2R^2*SLHA::EPSE[2.,2.]*SLHA::EPSE[2.,3.]+0.5*v2I^2*SLHA::EPSE[3.,2.]*SLHA::EPSE[3.,3.]+v2I*v2R*SLHA::EPSE[3.,2.]*SLHA::EPSE[3.,3.]+0.5*v2R^2*SLHA::EPSE[3.,2.]*SLHA::EPSE[3.,3.]+0.5*v1I*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1R*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1I*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1R*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1I*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1R*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1I*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1R*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1I^2*SLHA::YE[1.,2.]*SLHA::YE[1.,3.]+v1I*v1R*SLHA::YE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1R^2*SLHA::YE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1I*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1R*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1I*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1R*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1I*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1R*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1I*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1R*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1I^2*SLHA::YE[2.,2.]*SLHA::YE[2.,3.]+v1I*v1R*SLHA::YE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1R^2*SLHA::YE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1I*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1R*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1I*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1R*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1I*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1R*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1I*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1R*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1I^2*SLHA::YE[3.,2.]*SLHA::YE[3.,3.]+v1I*v1R*SLHA::YE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1R^2*SLHA::YE[3.,2.]*SLHA::YE[3.,3.]), 
(0.5*v2I^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,3.]+v2I*v2R*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,3.]+0.5*v2R^2*SLHA::EPSE[1.,1.]*SLHA::EPSE[1.,3.]+0.5*v2I^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,3.]+v2I*v2R*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,3.]+0.5*v2R^2*SLHA::EPSE[2.,1.]*SLHA::EPSE[2.,3.]+0.5*v2I^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,3.]+v2I*v2R*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,3.]+0.5*v2R^2*SLHA::EPSE[3.,1.]*SLHA::EPSE[3.,3.]+0.5*v1I*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1R*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1I*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1R*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,1.]+0.5*v1I*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1R*v2I*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1I*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1R*v2R*SLHA::EPSE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1I^2*SLHA::YE[1.,1.]*SLHA::YE[1.,3.]+v1I*v1R*SLHA::YE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1R^2*SLHA::YE[1.,1.]*SLHA::YE[1.,3.]+0.5*v1I*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1R*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1I*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1R*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,1.]+0.5*v1I*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1R*v2I*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1I*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1R*v2R*SLHA::EPSE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1I^2*SLHA::YE[2.,1.]*SLHA::YE[2.,3.]+v1I*v1R*SLHA::YE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1R^2*SLHA::YE[2.,1.]*SLHA::YE[2.,3.]+0.5*v1I*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1R*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1I*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1R*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,1.]+0.5*v1I*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1R*v2I*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1I*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1R*v2R*SLHA::EPSE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1I^2*SLHA::YE[3.,1.]*SLHA::YE[3.,3.]+v1I*v1R*SLHA::YE[3.,1.]*SLHA::YE[3.,3.]+0.5*v1R^2*SLHA::YE[3.,1.]*SLHA::YE[3.,3.]), 
(0.5*v2I^2*SLHA::EPSE[1.,2.]*SLHA::EPSE[1.,3.]+v2I*v2R*SLHA::EPSE[1.,2.]*SLHA::EPSE[1.,3.]+0.5*v2R^2*SLHA::EPSE[1.,2.]*SLHA::EPSE[1.,3.]+0.5*v2I^2*SLHA::EPSE[2.,2.]*SLHA::EPSE[2.,3.]+v2I*v2R*SLHA::EPSE[2.,2.]*SLHA::EPSE[2.,3.]+0.5*v2R^2*SLHA::EPSE[2.,2.]*SLHA::EPSE[2.,3.]+0.5*v2I^2*SLHA::EPSE[3.,2.]*SLHA::EPSE[3.,3.]+v2I*v2R*SLHA::EPSE[3.,2.]*SLHA::EPSE[3.,3.]+0.5*v2R^2*SLHA::EPSE[3.,2.]*SLHA::EPSE[3.,3.]+0.5*v1I*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1R*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1I*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1R*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,2.]+0.5*v1I*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1R*v2I*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1I*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1R*v2R*SLHA::EPSE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1I^2*SLHA::YE[1.,2.]*SLHA::YE[1.,3.]+v1I*v1R*SLHA::YE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1R^2*SLHA::YE[1.,2.]*SLHA::YE[1.,3.]+0.5*v1I*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1R*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1I*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1R*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,2.]+0.5*v1I*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1R*v2I*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1I*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1R*v2R*SLHA::EPSE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1I^2*SLHA::YE[2.,2.]*SLHA::YE[2.,3.]+v1I*v1R*SLHA::YE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1R^2*SLHA::YE[2.,2.]*SLHA::YE[2.,3.]+0.5*v1I*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1R*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1I*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1R*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,2.]+0.5*v1I*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1R*v2I*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1I*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1R*v2R*SLHA::EPSE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1I^2*SLHA::YE[3.,2.]*SLHA::YE[3.,3.]+v1I*v1R*SLHA::YE[3.,2.]*SLHA::YE[3.,3.]+0.5*v1R^2*SLHA::YE[3.,2.]*SLHA::YE[3.,3.]), 
(0.5*v2I^2*SLHA::EPSE[1.,3.]^2+v2I*v2R*SLHA::EPSE[1.,3.]^2+0.5*v2R^2*SLHA::EPSE[1.,3.]^2+0.5*v2I^2*SLHA::EPSE[2.,3.]^2+v2I*v2R*SLHA::EPSE[2.,3.]^2+0.5*v2R^2*SLHA::EPSE[2.,3.]^2+0.5*v2I^2*SLHA::EPSE[3.,3.]^2+v2I*v2R*SLHA::EPSE[3.,3.]^2+0.5*v2R^2*SLHA::EPSE[3.,3.]^2+v1I*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,3.]+v1R*v2I*SLHA::EPSE[1.,3.]*SLHA::YE[1.,3.]+v1I*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,3.]+v1R*v2R*SLHA::EPSE[1.,3.]*SLHA::YE[1.,3.]+0.5*v1I^2*SLHA::YE[1.,3.]^2+v1I*v1R*SLHA::YE[1.,3.]^2+0.5*v1R^2*SLHA::YE[1.,3.]^2+v1I*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,3.]+v1R*v2I*SLHA::EPSE[2.,3.]*SLHA::YE[2.,3.]+v1I*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,3.]+v1R*v2R*SLHA::EPSE[2.,3.]*SLHA::YE[2.,3.]+0.5*v1I^2*SLHA::YE[2.,3.]^2+v1I*v1R*SLHA::YE[2.,3.]^2+0.5*v1R^2*SLHA::YE[2.,3.]^2+v1I*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,3.]+v1R*v2I*SLHA::EPSE[3.,3.]*SLHA::YE[3.,3.]+v1I*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,3.]+v1R*v2R*SLHA::EPSE[3.,3.]*SLHA::YE[3.,3.]+0.5*v1I^2*SLHA::YE[3.,3.]^2+v1I*v1R*SLHA::YE[3.,3.]^2+0.5*v1R^2*SLHA::YE[3.,3.]^2) 
</mass-squared_matrix> 
 
<mass-squared_matrix 
particle="VP, VZ"  rotationmatrix="ZZ"  spin="vector" factor="3" constant="5/6" >
(0.25*v1I^2*SLHA::GAUGE[1.]^2+0.5*v1I*v1R*SLHA::GAUGE[1.]^2+0.25*v1R^2*SLHA::GAUGE[1.]^2+0.25*v2I^2*SLHA::GAUGE[1.]^2+0.5*v2I*v2R*SLHA::GAUGE[1.]^2+0.25*v2R^2*SLHA::GAUGE[1.]^2), 
(-0.25*v1I^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.5*v1I*v1R*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.25*v1R^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.25*v2I^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.5*v2I*v2R*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.25*v2R^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]), 
(-0.25*v1I^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.5*v1I*v1R*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.25*v1R^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.25*v2I^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.5*v2I*v2R*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]-0.25*v2R^2*SLHA::GAUGE[1.]*SLHA::GAUGE[2.]), 
(0.25*v1I^2*SLHA::GAUGE[2.]^2+0.5*v1I*v1R*SLHA::GAUGE[2.]^2+0.25*v1R^2*SLHA::GAUGE[2.]^2+0.25*v2I^2*SLHA::GAUGE[2.]^2+0.5*v2I*v2R*SLHA::GAUGE[2.]^2+0.25*v2R^2*SLHA::GAUGE[2.]^2) 
</mass-squared_matrix> 
 
<mass-squared_matrix 
particle="VWm"  rotationmatrix="ZW"  spin="vector" factor="3" constant="5/6" >
(0.25*v1I^2*SLHA::GAUGE[2.]^2+0.5*v1I*v1R*SLHA::GAUGE[2.]^2+0.25*v1R^2*SLHA::GAUGE[2.]^2+0.25*v2I^2*SLHA::GAUGE[2.]^2+0.5*v2I*v2R*SLHA::GAUGE[2.]^2+0.25*v2R^2*SLHA::GAUGE[2.]^2), 
(0.), 
(0.), 
(0.25*v1I^2*SLHA::GAUGE[2.]^2+0.5*v1I*v1R*SLHA::GAUGE[2.]^2+0.25*v1R^2*SLHA::GAUGE[2.]^2+0.25*v2I^2*SLHA::GAUGE[2.]^2+0.5*v2I*v2R*SLHA::GAUGE[2.]^2+0.25*v2R^2*SLHA::GAUGE[2.]^2) 
</mass-squared_matrix> 
 
</Vevacious_stuff>
Sorry for the inconvenience, I hope you can figure out what's going wrong.
Cheers,
Alexander
Top
FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: Error using MakeVevacious for THDMCPV

Post by FStaub » 3. May 2016, 14:28

Hi,

I see that you just have two tadpole equations in the files, but four VEVs. That's most likely the origin your problem.
I think the reason why you get only two equations is that you only run MakeVevacious[]: By default, all parameters are taken to be real to speed up the running time of Vevacious. However, in that limit two equations become just 0=0. Therefore, you need to run something like

Code: Select all

MakeVevacious[ComplexParameters -> {Lambda6, Lambda7}]
Cheers,
Florian
Top
AWlotzka

Re: Error using MakeVevacious for THDMCPV

Post by AWlotzka » 4. May 2016, 11:53

Hello,

thanks, I tried it out. Now in the THDMCPV.vin file there are four tadpoles instead of only two, but in the current_tadpoles.in still only two of them appear, the others are not passed to the file.
Seems like Vevacious is not detecting them.
The error message I get is still the same.
The exerpt of the .vin-file and the HOM4PS-file are shown below.

Thanks again,
cheers,
Alexander

THDMCPV.vin

Code: Select all

<tadpoles> 
{ 
(v1I*SLHA::HMIX[20.]) 
+(v1R*SLHA::HMIX[20.]) 
+(-1.*v2I*SLHA::HMIX[22.]) 
+(-1.*v2R*SLHA::HMIX[22.]) 
+(v1I^3*SLHA::HMIX[31.]) 
+(3.*v1I^2*v1R*SLHA::HMIX[31.]) 
+(3.*v1I*v1R^2*SLHA::HMIX[31.]) 
+(v1R^3*SLHA::HMIX[31.]) 
+(0.5*v1I*v2I^2*SLHA::HMIX[33.]) 
+(0.5*v1R*v2I^2*SLHA::HMIX[33.]) 
+(v1I*v2I*v2R*SLHA::HMIX[33.]) 
+(v1R*v2I*v2R*SLHA::HMIX[33.]) 
+(0.5*v1I*v2R^2*SLHA::HMIX[33.]) 
+(0.5*v1R*v2R^2*SLHA::HMIX[33.]) 
+(0.5*v1I*v2I^2*SLHA::HMIX[34.]) 
+(0.5*v1R*v2I^2*SLHA::HMIX[34.]) 
+(v1I*v2I*v2R*SLHA::HMIX[34.]) 
+(v1R*v2I*v2R*SLHA::HMIX[34.]) 
+(0.5*v1I*v2R^2*SLHA::HMIX[34.]) 
+(0.5*v1R*v2R^2*SLHA::HMIX[34.]) 
+(0.5*v1I*v2I^2*SLHA::HMIX[35.]) 
+(0.5*v1R*v2I^2*SLHA::HMIX[35.]) 
+(v1I*v2I*v2R*SLHA::HMIX[35.]) 
+(v1R*v2I*v2R*SLHA::HMIX[35.]) 
+(0.5*v1I*v2R^2*SLHA::HMIX[35.]) 
+(0.5*v1R*v2R^2*SLHA::HMIX[35.]) 
+(1.5*v1I^2*v2I*SLHA::HMIX[36.]) 
+(3.*v1I*v1R*v2I*SLHA::HMIX[36.]) 
+(1.5*v1R^2*v2I*SLHA::HMIX[36.]) 
+(1.5*v1I^2*v2R*SLHA::HMIX[36.]) 
+(3.*v1I*v1R*v2R*SLHA::HMIX[36.]) 
+(1.5*v1R^2*v2R*SLHA::HMIX[36.]) 
+(0.5*v2I^3*SLHA::HMIX[37.]) 
+(1.5*v2I^2*v2R*SLHA::HMIX[37.]) 
+(1.5*v2I*v2R^2*SLHA::HMIX[37.]) 
+(0.5*v2R^3*SLHA::HMIX[37.]) 
; 
(v2I*SLHA::HMIX[21.]) 
+(v2R*SLHA::HMIX[21.]) 
+(-1.*v1I*SLHA::HMIX[22.]) 
+(-1.*v1R*SLHA::HMIX[22.]) 
+(v2I^3*SLHA::HMIX[32.]) 
+(3.*v2I^2*v2R*SLHA::HMIX[32.]) 
+(3.*v2I*v2R^2*SLHA::HMIX[32.]) 
+(v2R^3*SLHA::HMIX[32.]) 
+(0.5*v1I^2*v2I*SLHA::HMIX[33.]) 
+(v1I*v1R*v2I*SLHA::HMIX[33.]) 
+(0.5*v1R^2*v2I*SLHA::HMIX[33.]) 
+(0.5*v1I^2*v2R*SLHA::HMIX[33.]) 
+(v1I*v1R*v2R*SLHA::HMIX[33.]) 
+(0.5*v1R^2*v2R*SLHA::HMIX[33.]) 
+(0.5*v1I^2*v2I*SLHA::HMIX[34.]) 
+(v1I*v1R*v2I*SLHA::HMIX[34.]) 
+(0.5*v1R^2*v2I*SLHA::HMIX[34.]) 
+(0.5*v1I^2*v2R*SLHA::HMIX[34.]) 
+(v1I*v1R*v2R*SLHA::HMIX[34.]) 
+(0.5*v1R^2*v2R*SLHA::HMIX[34.]) 
+(0.5*v1I^2*v2I*SLHA::HMIX[35.]) 
+(v1I*v1R*v2I*SLHA::HMIX[35.]) 
+(0.5*v1R^2*v2I*SLHA::HMIX[35.]) 
+(0.5*v1I^2*v2R*SLHA::HMIX[35.]) 
+(v1I*v1R*v2R*SLHA::HMIX[35.]) 
+(0.5*v1R^2*v2R*SLHA::HMIX[35.]) 
+(0.5*v1I^3*SLHA::HMIX[36.]) 
+(1.5*v1I^2*v1R*SLHA::HMIX[36.]) 
+(1.5*v1I*v1R^2*SLHA::HMIX[36.]) 
+(0.5*v1R^3*SLHA::HMIX[36.]) 
+(1.5*v1I*v2I^2*SLHA::HMIX[37.]) 
+(1.5*v1R*v2I^2*SLHA::HMIX[37.]) 
+(3.*v1I*v2I*v2R*SLHA::HMIX[37.]) 
+(3.*v1R*v2I*v2R*SLHA::HMIX[37.]) 
+(1.5*v1I*v2R^2*SLHA::HMIX[37.]) 
+(1.5*v1R*v2R^2*SLHA::HMIX[37.]) 
; 
(0.5*v1I^2*v2I*SLHA::IMHMIX[36.]) 
+(v1I*v1R*v2I*SLHA::IMHMIX[36.]) 
+(0.5*v1R^2*v2I*SLHA::IMHMIX[36.]) 
+(0.5*v1I^2*v2R*SLHA::IMHMIX[36.]) 
+(v1I*v1R*v2R*SLHA::IMHMIX[36.]) 
+(0.5*v1R^2*v2R*SLHA::IMHMIX[36.]) 
+(0.5*v2I^3*SLHA::IMHMIX[37.]) 
+(1.5*v2I^2*v2R*SLHA::IMHMIX[37.]) 
+(1.5*v2I*v2R^2*SLHA::IMHMIX[37.]) 
+(0.5*v2R^3*SLHA::IMHMIX[37.]) 
; 
(-0.5*v1I^3*SLHA::IMHMIX[36.]) 
+(-1.5*v1I^2*v1R*SLHA::IMHMIX[36.]) 
+(-1.5*v1I*v1R^2*SLHA::IMHMIX[36.]) 
+(-0.5*v1R^3*SLHA::IMHMIX[36.]) 
+(-0.5*v1I*v2I^2*SLHA::IMHMIX[37.]) 
+(-0.5*v1R*v2I^2*SLHA::IMHMIX[37.]) 
+(-1.*v1I*v2I*v2R*SLHA::IMHMIX[37.]) 
+(-1.*v1R*v2I*v2R*SLHA::IMHMIX[37.]) 
+(-0.5*v1I*v2R^2*SLHA::IMHMIX[37.]) 
+(-0.5*v1R*v2R^2*SLHA::IMHMIX[37.]) 
; 
} 
</tadpoles>
current_tadpoles.in

Code: Select all

{ (a*(1.81498477E+06)) + (b*(1.81498477E+06)) + (-1.*c*(4.00000000E+04)) + (-1.*d*(4.00000000E+04)) + (a^3*(1.00000000E-01)) + (3.*a^2*b*(1.00000000E-01)) + (3.*a*b^2*(1.00000000E-01)) + (b^3*(1.00000000E-01)) + (0.5*a*c^2*(1.10000000E+00)) + (0.5*b*c^2*(1.10000000E+00)) + (a*c*d*(1.10000000E+00)) + (b*c*d*(1.10000000E+00)) + (0.5*a*d^2*(1.10000000E+00)) + (0.5*b*d^2*(1.10000000E+00)) + (0.5*a*c^2*(-5.00000000E-01)) + (0.5*b*c^2*(-5.00000000E-01)) + (a*c*d*(-5.00000000E-01)) + (b*c*d*(-5.00000000E-01)) + (0.5*a*d^2*(-5.00000000E-01)) + (0.5*b*d^2*(-5.00000000E-01)) + (0.5*a*c^2*(5.00000000E-01)) + (0.5*b*c^2*(5.00000000E-01)) + (a*c*d*(5.00000000E-01)) + (b*c*d*(5.00000000E-01)) + (0.5*a*d^2*(5.00000000E-01)) + (0.5*b*d^2*(5.00000000E-01)) + (1.5*a^2*c*(1.00000000E-01)) + (3.*a*b*c*(1.00000000E-01)) + (1.5*b^2*c*(1.00000000E-01)) + (1.5*a^2*d*(1.00000000E-01)) + (3.*a*b*d*(1.00000000E-01)) + (1.5*b^2*d*(1.00000000E-01)) + (0.5*c^3*(1.00000000E-01)) + (1.5*c^2*d*(1.00000000E-01)) + (1.5*c*d^2*(1.00000000E-01)) + (0.5*d^3*(1.00000000E-01));
(c*(-7.27320000E+03)) + (d*(-7.27320000E+03)) + (-1.*a*(4.00000000E+04)) + (-1.*b*(4.00000000E+04)) + (c^3*(1.30000000E-01)) + (3.*c^2*d*(1.30000000E-01)) + (3.*c*d^2*(1.30000000E-01)) + (d^3*(1.30000000E-01)) + (0.5*a^2*c*(1.10000000E+00)) + (a*b*c*(1.10000000E+00)) + (0.5*b^2*c*(1.10000000E+00)) + (0.5*a^2*d*(1.10000000E+00)) + (a*b*d*(1.10000000E+00)) + (0.5*b^2*d*(1.10000000E+00)) + (0.5*a^2*c*(-5.00000000E-01)) + (a*b*c*(-5.00000000E-01)) + (0.5*b^2*c*(-5.00000000E-01)) + (0.5*a^2*d*(-5.00000000E-01)) + (a*b*d*(-5.00000000E-01)) + (0.5*b^2*d*(-5.00000000E-01)) + (0.5*a^2*c*(5.00000000E-01)) + (a*b*c*(5.00000000E-01)) + (0.5*b^2*c*(5.00000000E-01)) + (0.5*a^2*d*(5.00000000E-01)) + (a*b*d*(5.00000000E-01)) + (0.5*b^2*d*(5.00000000E-01)) + (0.5*a^3*(1.00000000E-01)) + (1.5*a^2*b*(1.00000000E-01)) + (1.5*a*b^2*(1.00000000E-01)) + (0.5*b^3*(1.00000000E-01)) + (1.5*a*c^2*(1.00000000E-01)) + (1.5*b*c^2*(1.00000000E-01)) + (3.*a*c*d*(1.00000000E-01)) + (3.*b*c*d*(1.00000000E-01)) + (1.5*a*d^2*(1.00000000E-01)) + (1.5*b*d^2*(1.00000000E-01)); }
Top
benoleary
Posts: 45
Joined: 3. May 2016, 10:49

Re: Error using MakeVevacious for THDMCPV

Post by benoleary » 4. May 2016, 11:59

Dear Alexander,

if you are using the same SPheno file as that which you sent to me, the problem is that the "missing" tadpoles only have terms which are dependent on the IMHMIX block, which is missing in the SPheno output file. In that case, Vevacious does detect them, but because it assumes that missing terms are zero, the tadpoles really do evaluate to zero equals zero.

Regards,
Ben
Top
AWlotzka

Re: Error using MakeVevacious for THDMCPV

Post by AWlotzka » 10. May 2016, 09:10

Hello again,

I tried it, but I think it's not working properly, yet.
HOM4PS2 never finds any roots and there is a MINUIT error:

Code: Select all

Setting up the VevaciousRunner object took 0.006822 seconds.
 Please enter the option of homotopy method
   1. The polyhedral homotopy
   2. The classical linear homotopy
 The classical linear homotopy
 The # of paths followed  =         81
 
 The # of roots =               0
 The # of real roots =          0
 The # of paths followed  =         81
 The # of blow_up =          0  (before endgame)
 The running time :          0.04400 s
 The solutions are listed in file "data.roots"
 
 Department of Mathematics, Michigan State University
 Contact Info: 
          Tsung-Lin Lee     : leetsung@msu.edu
          Tien-Yien Li      : li@math.msu.edu
          Chih-Hsiung Tsai  : tsaichih@msu.edu
 

Finding all the tree-level extrema and parsing the results took 0.078651 seconds.
rm: cannot remove ‘./Vevacious.pyc’: No such file or directory
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:580: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:587: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:594: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:601: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:608: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:615: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:622: VisibleDeprecationWarning: using a non-integer number instead of an integer will result in an error in the future
  massSquaredElementArray.size ), -1 )
trying to minimize { { v1I -> ( 157.381849 ), v1R -> ( 0.0 ), v2I -> ( 189.355545 ), v2R -> ( 0.0 ) }, TreeLevelPotentialValue -> (-199902404.64+0j), EffectivePotentialValue -> (380432708.034+0j) } (T = 0.0)
/itp/scratch/wlotzka/project/Vevacious/bin/VevaciousParameterDependent.py:991: RuntimeWarning: Cannot convert fcn(*arg) to double 
fcn is called with following arguments:
    a = +0.983637
    b = +0.000000
    c = +1.183472
    d = +0.000000

  minuitObject.migrad()
The MINUIT wrapper threw an exception:
can't convert complex to float
while trying to minimize the potential starting at { v1I -> ( 157.381849 ), v1R -> ( 0.0 ), v2I -> ( 189.355545 ), v2R -> ( 0.0 ) }. MINUIT stopped at { v1I -> ( 157.381849 ), v1R -> ( 0.0 ), v2I -> ( 189.355545 ), v2R -> ( 0.0 ) }.  MINUIT's estimate of how much deeper it should go is 655360000.0 GeV^4.
Taking the starting point as a vacuum candidate.
DSB vacuum: { { v1I -> ( 157.381849 ), v1R -> ( 0.0 ), v2I -> ( 189.355545 ), v2R -> ( 0.0 ) }, TreeLevelPotentialValue -> (-199902404.64+0j), EffectivePotentialValue -> (380432708.034+0j) }
DSB vacuum is stable (as far as the model file allows).
Result summary (not recapping warnings):
  <reference version="1.2.00" citation="arXiv:1307.1477, arXiv:1405.7376 (hep-ph)" />
  <stability slha_summary_code="1" > stable </stability>
  <quantum_stability> stable </quantum_stability>
  <thermal_stability> stable </thermal_stability>
  <global_minimum   relative_depth="(380432708.034+0j)" v1I="157.381849" v1R="0.0" v2I="189.355545" v2R="0.0" />
  <input_minimum   relative_depth="(380432708.034+0j)" v1I="157.381849" v1R="0.0" v2I="189.355545" v2R="0.0" />
  <lifetime  action_calculation="unnecessary" > -1.0 </lifetime>
 <tunneling_temperature survival_probability="1.0" > unnecessary </tunneling_temperature>

Minimizing the one-loop potential, calculating the tunneling time or its upper bound, and printing the results took 0.216761 seconds.

Vevacious finished. Total running time was 0.304178 seconds.
Here is the input:

Code: Select all

# SUSY Les Houches Accord 2 - THDM-CPV Spectrum + Decays + Flavor Observables
# SPheno module generated by SARAH
# ----------------------------------------------------------------------
# SPheno v3.3.8  
#   W. Porod, Comput. Phys. Commun. 153 (2003) 275-315, hep-ph/0301101
#   W. Porod, F.Staub, Comput.Phys.Commun.183 (2012) 2458-2469, arXiv:1104.1573
# SARAH: 4.8.5   
#   F. Staub; arXiv:0806.0538 (online manual)
#   F. Staub; Comput. Phys. Commun. 181 (2010) 1077-1086; arXiv:0909.2863
#   F. Staub; Comput. Phys. Commun. 182 (2011)  808-833; arXiv:1002.0840
#   F. Staub; Comput. Phys. Commun. 184 (2013)  1792-1809; arXiv:1207.0906
#   F. Staub; Comput. Phys. Commun. 185 (2014)  1773-1790; arXiv:1309.7223 
# Including the calculation of flavor observables based on the FlavorKit 
#   W. Porod, F. Staub, A. Vicente; Eur.Phys.J. C74 (2014) 8, 2992; arXiv:1405.1434 
# Two-loop masss corrections to Higgs fields based on 
#   M. D. Goodsell, K. Nickel, F. Staub; arXiv:1411.0675 
#   M. D. Goodsell, K. Nickel, F. Staub; arXiv:1503.03098
#  
# in case of problems send email to florian.staub@cern.ch and goodsell@lpthe.jussieu.fr 
# ----------------------------------------------------------------------
# Created: 10.05.2016,  08:47
Block SPINFO         # Program information
     1   SPhenoSARAH      # spectrum calculator
     2   v3.3.8      # version number of SPheno
     9   4.8.5       # version number of SARAH
Block MODSEL  # Input parameters
     1  1   #  GUT scale input
     2  1   #  Boundary conditions 
Block MINPAR  # Input parameters
    1    1.20000000E+00  # Lambda1Input
    2    1.50000000E+00  # Lambda2Input
    3    1.60000000E+00  # Lambda3Input
    4   -1.20000000E+00  # Lambda4Input
    5   -2.50000000E+00  # Lambda5Input
    6    0.00000000E+00  # Lambda6Input
    7    0.00000000E+00  # Lambda7Input
    9    3.29328000E+04  # M12input
   10    1.20316000E+00  # TanBeta
   11    0.00000000E+00  # etaInput
Block IMMINPAR  # Input parameters
    5    8.00000000E-01  # Lambda5Input
    6    1.00000000E-02  # Lambda6Input
    7    3.00000000E-01  # Lambda7Input
Block gaugeGUT Q= -1.00000000E+00  # (GUT scale)
   1    0.00000000E+00  # g1(Q)^DRbar
   2    0.00000000E+00  # g2(Q)^DRbar
   3    0.00000000E+00  # g3(Q)^DRbar
Block SMINPUTS  # SM parameters
         2     1.16637000E-05  # G_mu [GeV^-2]
         3     1.18700000E-01  # alpha_s(MZ)^MSbar
         4     9.11887000E+01  # m_Z(pole)
         5     4.18000000E+00  # m_b(m_b), MSbar
         6     1.73500000E+02  # m_t(pole)
         7     1.77669000E+00  # m_tau(pole)
Block GAUGE Q=  1.60000000E+02  # (Renormalization Scale)
   1    3.55734740E-01  # g1
   2    6.62881508E-01  # g2
   3    1.22132231E+00  # g3
Block HMIX Q=  1.60000000E+02  # (Renormalization Scale)
  36    0.00000000E+00  # Lam6
  35   -2.50000000E+00  # Lam5
  37    0.00000000E+00  # Lam7
  31    1.20000000E+00  # Lam1
  34   -1.20000000E+00  # Lam4
  33    1.60000000E+00  # Lam3
  32    1.50000000E+00  # Lam2
  22    3.29328008E+04  # M12
  20    4.75488716E+04  # M112
  21   -4.03863752E+02  # M222
 102    1.57381849E+02  # vd
 103    1.89355545E+02  # vu
1001    1.57381849E+02  # v1R
1002    0.00000000E+00  # v1I
1003    1.89355545E+02  # v2R
1004    0.00000000E+00  # v2I
   3    2.46220569E+02  # v
  10    8.77351122E-01  # betaH
Block IMHMIX Q=  1.60000000E+02  # (Renormalization Scale)
  36    1.00000000E-02  # Lam6
  35    8.00000000E-01  # Lam5
  37    3.00000000E-01  # Lam7
  31    0.00000000E+00  # Lam1
  34    0.00000000E+00  # Lam4
  33    0.00000000E+00  # Lam3
  32    0.00000000E+00  # Lam2
  22   -6.41827669E+03  # M12
  20   -0.00000000E+00  # M112
  21   -1.24676704E-12  # M222
Block TREEHMIX Q=  1.60000000E+02  # (Renormalization Scale)
  22    3.29328008E+04  # M12
  20    4.75488716E+04  # M112
  21   -4.03863752E+02  # M222
Block TREEIMHMIX Q=  1.60000000E+02  # (Renormalization Scale)
  22   -6.41827669E+03  # M12
  20   -0.00000000E+00  # M112
  21   -1.24676704E-12  # M222
Block LOOPHMIX Q=  1.60000000E+02  # (Renormalization Scale)
  22    0.00000000E+00  # M12
  20    0.00000000E+00  # M112
  21    0.00000000E+00  # M222
Block epsU Q=  1.60000000E+02  # (Renormalization Scale)
  1  1     0.00000000E+00   # Real(epsU(1,1),dp)
  1  2     0.00000000E+00   # Real(epsU(1,2),dp)
  1  3     0.00000000E+00   # Real(epsU(1,3),dp)
  2  1     0.00000000E+00   # Real(epsU(2,1),dp)
  2  2     0.00000000E+00   # Real(epsU(2,2),dp)
  2  3     0.00000000E+00   # Real(epsU(2,3),dp)
  3  1     0.00000000E+00   # Real(epsU(3,1),dp)
  3  2     0.00000000E+00   # Real(epsU(3,2),dp)
  3  3     0.00000000E+00   # Real(epsU(3,3),dp)
Block Yu Q=  1.60000000E+02  # (Renormalization Scale)
  1  1     1.86714041E-05   # Real(Yu(1,1),dp)
  1  2     0.00000000E+00   # Real(Yu(1,2),dp)
  1  3     0.00000000E+00   # Real(Yu(1,3),dp)
  2  1     0.00000000E+00   # Real(Yu(2,1),dp)
  2  2     9.48507331E-03   # Real(Yu(2,2),dp)
  2  3     0.00000000E+00   # Real(Yu(2,3),dp)
  3  1     0.00000000E+00   # Real(Yu(3,1),dp)
  3  2     0.00000000E+00   # Real(Yu(3,2),dp)
  3  3     1.29579545E+00   # Real(Yu(3,3),dp)
Block Yd Q=  1.60000000E+02  # (Renormalization Scale)
  1  1     4.49293732E-05   # Real(Yd(1,1),dp)
  1  2     0.00000000E+00   # Real(Yd(1,2),dp)
  1  3     0.00000000E+00   # Real(Yd(1,3),dp)
  2  1     0.00000000E+00   # Real(Yd(2,1),dp)
  2  2     8.53658091E-04   # Real(Yd(2,2),dp)
  2  3     0.00000000E+00   # Real(Yd(2,3),dp)
  3  1     0.00000000E+00   # Real(Yd(3,1),dp)
  3  2     0.00000000E+00   # Real(Yd(3,2),dp)
  3  3     3.75609560E-02   # Real(Yd(3,3),dp)
Block Ye Q=  1.60000000E+02  # (Renormalization Scale)
  1  1     4.59177233E-06   # Real(Ye(1,1),dp)
  1  2     0.00000000E+00   # Real(Ye(1,2),dp)
  1  3     0.00000000E+00   # Real(Ye(1,3),dp)
  2  1     0.00000000E+00   # Real(Ye(2,1),dp)
  2  2     9.49432881E-04   # Real(Ye(2,2),dp)
  2  3     0.00000000E+00   # Real(Ye(2,3),dp)
  3  1     0.00000000E+00   # Real(Ye(3,1),dp)
  3  2     0.00000000E+00   # Real(Ye(3,2),dp)
  3  3     1.59651136E-02   # Real(Ye(3,3),dp)
Block epsD Q=  1.60000000E+02  # (Renormalization Scale)
  1  1     0.00000000E+00   # Real(epsD(1,1),dp)
  1  2     0.00000000E+00   # Real(epsD(1,2),dp)
  1  3     0.00000000E+00   # Real(epsD(1,3),dp)
  2  1     0.00000000E+00   # Real(epsD(2,1),dp)
  2  2     0.00000000E+00   # Real(epsD(2,2),dp)
  2  3     0.00000000E+00   # Real(epsD(2,3),dp)
  3  1     0.00000000E+00   # Real(epsD(3,1),dp)
  3  2     0.00000000E+00   # Real(epsD(3,2),dp)
  3  3     0.00000000E+00   # Real(epsD(3,3),dp)
Block epsE Q=  1.60000000E+02  # (Renormalization Scale)
  1  1     0.00000000E+00   # Real(epsE(1,1),dp)
  1  2     0.00000000E+00   # Real(epsE(1,2),dp)
  1  3     0.00000000E+00   # Real(epsE(1,3),dp)
  2  1     0.00000000E+00   # Real(epsE(2,1),dp)
  2  2     0.00000000E+00   # Real(epsE(2,2),dp)
  2  3     0.00000000E+00   # Real(epsE(2,3),dp)
  3  1     0.00000000E+00   # Real(epsE(3,1),dp)
  3  2     0.00000000E+00   # Real(epsE(3,2),dp)
  3  3     0.00000000E+00   # Real(epsE(3,3),dp)
Block AddPars Q=1.0000E00 # Additional parameters not appearing elsewhere; scale Q just a dummy
   1    0.00000000E+00  # eta 
Block MASS  # Mass spectrum
#   PDG code      mass          particle
        25     1.36702100E+02  # hh_2
        35     4.54001498E+02  # hh_3
        36     4.77241632E+02  # hh_4
        37     4.23262096E+02  # Hm_2
        23     9.11887000E+01  # VZ
        24     8.03497269E+01  # VWm
         1     5.00000000E-03  # Fd_1
         3     9.50000000E-02  # Fd_2
         5     4.18000000E+00  # Fd_3
         2     2.50000000E-03  # Fu_1
         4     1.27000000E+00  # Fu_2
         6     1.73500000E+02  # Fu_3
        11     5.10998930E-04  # Fe_1
        13     1.05658372E-01  # Fe_2
        15     1.77669000E+00  # Fe_3
Block SCALARMIX Q=  1.60000000E+02  # ()
  1  1    -4.12227385E-16   # ZH(1,1)
  1  2    -4.09515075E-16   # ZH(1,2)
  1  3    -6.39190500E-01   # ZH(1,3)
  1  4    -7.69048441E-01   # ZH(1,4)
  2  1    -7.69985524E-01   # ZH(2,1)
  2  2    -6.33620377E-01   # ZH(2,2)
  2  3    -5.77939274E-02   # ZH(2,3)
  2  4     4.80351137E-02   # ZH(2,4)
  3  1     4.67806314E-01   # ZH(3,1)
  3  2    -6.40694701E-01   # ZH(3,2)
  3  3     4.68215710E-01   # ZH(3,3)
  3  4    -3.89154984E-01   # ZH(3,4)
  4  1     4.33911910E-01   # ZH(4,1)
  4  2    -4.33630623E-01   # ZH(4,2)
  4  3    -6.07346208E-01   # ZH(4,3)
  4  4     5.04792553E-01   # ZH(4,4)
Block CHARGEMIX Q=  1.60000000E+02  # ()
  1  1     6.39190500E-01   # Real(ZP(1,1),dp)
  1  2     7.69048441E-01   # Real(ZP(1,2),dp)
  2  1    -7.69048441E-01   # Real(ZP(2,1),dp)
  2  2     6.39190500E-01   # Real(ZP(2,2),dp)
Block IMCHARGEMIX Q=  1.60000000E+02  # ()
  1  1    -2.05684355E-17   # Aimag(ZP(1,1))
  1  2     0.00000000E+00   # Aimag(ZP(1,2))
  2  1     2.47471189E-17   # Aimag(ZP(2,1))
  2  2     0.00000000E+00   # Aimag(ZP(2,2))
Block UDLMIX Q=  1.60000000E+02  # ()
  1  1     1.00000000E+00   # Real(ZDL(1,1),dp)
  1  2     0.00000000E+00   # Real(ZDL(1,2),dp)
  1  3     0.00000000E+00   # Real(ZDL(1,3),dp)
  2  1     0.00000000E+00   # Real(ZDL(2,1),dp)
  2  2     1.00000000E+00   # Real(ZDL(2,2),dp)
  2  3     0.00000000E+00   # Real(ZDL(2,3),dp)
  3  1     0.00000000E+00   # Real(ZDL(3,1),dp)
  3  2     0.00000000E+00   # Real(ZDL(3,2),dp)
  3  3     1.00000000E+00   # Real(ZDL(3,3),dp)
Block UDRMIX Q=  1.60000000E+02  # ()
  1  1     1.00000000E+00   # Real(ZDR(1,1),dp)
  1  2     0.00000000E+00   # Real(ZDR(1,2),dp)
  1  3     0.00000000E+00   # Real(ZDR(1,3),dp)
  2  1     0.00000000E+00   # Real(ZDR(2,1),dp)
  2  2     1.00000000E+00   # Real(ZDR(2,2),dp)
  2  3     0.00000000E+00   # Real(ZDR(2,3),dp)
  3  1     0.00000000E+00   # Real(ZDR(3,1),dp)
  3  2     0.00000000E+00   # Real(ZDR(3,2),dp)
  3  3     1.00000000E+00   # Real(ZDR(3,3),dp)
Block UULMIX Q=  1.60000000E+02  # ()
  1  1    -1.00000000E+00   # Real(ZUL(1,1),dp)
  1  2     0.00000000E+00   # Real(ZUL(1,2),dp)
  1  3     0.00000000E+00   # Real(ZUL(1,3),dp)
  2  1     0.00000000E+00   # Real(ZUL(2,1),dp)
  2  2    -1.00000000E+00   # Real(ZUL(2,2),dp)
  2  3     0.00000000E+00   # Real(ZUL(2,3),dp)
  3  1     0.00000000E+00   # Real(ZUL(3,1),dp)
  3  2     0.00000000E+00   # Real(ZUL(3,2),dp)
  3  3    -1.00000000E+00   # Real(ZUL(3,3),dp)
Block UURMIX Q=  1.60000000E+02  # ()
  1  1     1.00000000E+00   # Real(ZUR(1,1),dp)
  1  2     0.00000000E+00   # Real(ZUR(1,2),dp)
  1  3     0.00000000E+00   # Real(ZUR(1,3),dp)
  2  1     0.00000000E+00   # Real(ZUR(2,1),dp)
  2  2     1.00000000E+00   # Real(ZUR(2,2),dp)
  2  3     0.00000000E+00   # Real(ZUR(2,3),dp)
  3  1     0.00000000E+00   # Real(ZUR(3,1),dp)
  3  2     0.00000000E+00   # Real(ZUR(3,2),dp)
  3  3     1.00000000E+00   # Real(ZUR(3,3),dp)
Block UELMIX Q=  1.60000000E+02  # ()
  1  1     1.00000000E+00   # Real(ZEL(1,1),dp)
  1  2     0.00000000E+00   # Real(ZEL(1,2),dp)
  1  3     0.00000000E+00   # Real(ZEL(1,3),dp)
  2  1     0.00000000E+00   # Real(ZEL(2,1),dp)
  2  2     1.00000000E+00   # Real(ZEL(2,2),dp)
  2  3     0.00000000E+00   # Real(ZEL(2,3),dp)
  3  1     0.00000000E+00   # Real(ZEL(3,1),dp)
  3  2     0.00000000E+00   # Real(ZEL(3,2),dp)
  3  3     1.00000000E+00   # Real(ZEL(3,3),dp)
Block UERMIX Q=  1.60000000E+02  # ()
  1  1     1.00000000E+00   # Real(ZER(1,1),dp)
  1  2     0.00000000E+00   # Real(ZER(1,2),dp)
  1  3     0.00000000E+00   # Real(ZER(1,3),dp)
  2  1     0.00000000E+00   # Real(ZER(2,1),dp)
  2  2     1.00000000E+00   # Real(ZER(2,2),dp)
  2  3     0.00000000E+00   # Real(ZER(2,3),dp)
  3  1     0.00000000E+00   # Real(ZER(3,1),dp)
  3  2     0.00000000E+00   # Real(ZER(3,2),dp)
  3  3     1.00000000E+00   # Real(ZER(3,3),dp)
Block SPheno # SPheno internal parameters 
         1    -1.00000000E+00  # ErrorLevel
         2     0.00000000E+00  # SPA_conventions
        11     1.00000000E+00  # Branching ratios
        13     1.00000000E+00  # 3 Body decays
        31    -1.00000000E+00  # GUT scale
        33     1.60000000E+02  # Renormalization scale
        34     1.00000000E-04  # Precision
        35     4.00000000E+01  # Iterations
        38     2.00000000E+00  # RGE level
        40     7.29735257E-03  # Alpha
        41     2.49520000E+00  # Gamma_Z
        42     2.06000000E+00  # Gamma_W
        50     1.00000000E+00  # Rotate negative fermion masses
        51     0.00000000E+00  # Switch to SCKM matrix
        52     0.00000000E+00  # Ignore negative masses
        53     0.00000000E+00  # Ignore negative masses at MZ
        55     0.00000000E+00  # Calculate one loop masses
        56     1.00000000E+00  # Calculate two-loop Higgs masses
        57     1.00000000E+00  # Calculate low energy
        60     1.00000000E+00  # Include kinetic mixing
        65     1.00000000E+00  # Solution of tadpole equation
And the Vevacious output blocks written to the spectrum file:

Code: Select all

BLOCK VEVACIOUSRESULTS # results from Vevacious version 1.2.00, documented in arXiv:1307.1477, arXiv:1405.7376 (hep-ph)
    0    0    1.00000000E+000    stable    # stability of input
    0    1    -1.00000000E+000    unnecessary    # tunneling time in Universe ages / calculation type
    0    2    0.00000000E+000    1.0    # estimated best tunneling temperature / survival probability at this temperature
    1    0    0.00000000E+000    relative_depth    # DSB vacuum potential energy
    1    1    1.57381849E+002    v1I    # DSB vacuum VEV
    1    2    0.00000000E+000    v1R    # DSB vacuum VEV
    1    3    1.89355545E+002    v2I    # DSB vacuum VEV
    1    4    0.00000000E+000    v2R    # DSB vacuum VEV
    2    0    0.00000000E+000    relative_depth    # panic vacuum potential
    2    1    1.57381849E+002    v1I    # panic vacuum VEV
    2    2    0.00000000E+000    v1R    # panic vacuum VEV
    2    3    1.89355545E+002    v2I    # panic vacuum VEV
    2    4    0.00000000E+000    v2R    # panic vacuum VEV
BLOCK VEVACIOUSWARNINGS # warnings from Vevacious
   1    The MINUIT wrapper threw an exception:can't convert complex to floatwhile trying to minimize the potential starting at { v1I -> ( 157.381849 ), v1R -> ( 0.0 ), v2I -> ( 189.355545 ), v2R -> ( 0.0 ) }. MINUIT stopped at { v1I -> ( 157.381849 ), v1R -> ( 0.0 ), v2I -> ( 189.355545 ), v2R -> ( 0.0 ) }.  MINUIT's estimate of how much deeper it should go is 655360000.0 GeV^4.Taking the starting point as a vacuum candidate.
Maybe there is a problem finding the roots and the MINUIT error is simply a follow up? Is the Vevacious runtime-warning important? I also realized that the vevs end up in v1I and v2I although they were in the v1R and v2R variables in the block HMIX...
For other models, this problem doesn't occur. Other input values didn't solve the problem.
I hope you can figure out whats going wrong. I can also provide you with the complete input files etc, if needed.

Thanks and cheers,
Alexander
Top
FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: Error using MakeVevacious for THDMCPV

Post by FStaub » 10. May 2016, 09:38

Hi Alexander,

I think it would be useful if you could attach a tar-file including
- the SARAH model file
- the Les Houches input for SPheno
- the Vevacious init file

I'll then try to run things the next days myself.

Cheers,
Florian
Top
benoleary
Posts: 45
Joined: 3. May 2016, 10:49

Re: Error using MakeVevacious for THDMCPV

Post by benoleary » 10. May 2016, 10:16

Hm, maybe also the code which I wrote to cope with spurious complex potential values when it should be a real-valued function isn't working (as the output mentions such values as "380432708.034+0j" which is obviously real but PyMinuit doesn't realize it).

However, there is also a problem if HOM4PS2 isn't producing output. Maybe Vevacious is producing incorrect HOM4PS2 input files, possibly because the complex parts of the tadpoles still have imaginary units in them, rather than being purely real constraints corresponding to the imaginary parts of the tadpoles.
Top
Post Reply

Return to “Vevacious”