Skip to content

Commit

Permalink
generalise MSSMtower
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Aug 24, 2016
1 parent 3f2d3ae commit ebcc1e6
Show file tree
Hide file tree
Showing 5 changed files with 177 additions and 56 deletions.
77 changes: 45 additions & 32 deletions model_files/MSSMMuBMu/FlexibleSUSY.m.in
View file
Expand Up @@ -8,32 +8,46 @@ MINPAR = {
};

EXTPAR = {
{0, Ms},
{14, Xtt},
{0, MSUSY},
{1, M1Input},
{2, M2Input},
{3, M3Input},
{4, MuInput},
{5, mAInput},
{25, TanBeta}
};

FSExtraInputParameters = {
{mq2Input, MSQ2IN, {3,3}}, (* 3x3 matrix *)
{mu2Input, MSU2IN, {3,3}}, (* 3x3 matrix *)
{md2Input, MSD2IN, {3,3}}, (* 3x3 matrix *)
{ml2Input, MSL2IN, {3,3}}, (* 3x3 matrix *)
{me2Input, MSE2IN, {3,3}}, (* 3x3 matrix *)
{AuInput, AUIN, {3,3}}, (* 3x3 matrix *)
{AdInput, ADIN, {3,3}}, (* 3x3 matrix *)
{AeInput, AEIN, {3,3}} (* 3x3 matrix *)
};

EWSBOutputParameters = { mHd2, mHu2 };

SUSYScale = Ms;
SUSYScale = MSUSY;

SUSYScaleFirstGuess = Ms;
SUSYScaleFirstGuess = MSUSY;

SUSYScaleInput = {
{MassB, Ms},
{MassWB, Ms},
{MassG, Ms},
{mq2, UNITMATRIX[3] Ms^2},
{mu2, UNITMATRIX[3] Ms^2},
{md2, UNITMATRIX[3] Ms^2},
{ml2, UNITMATRIX[3] Ms^2},
{me2, UNITMATRIX[3] Ms^2},
{\[Mu], Ms},
{B[\[Mu]], Sqr[Ms]/(TanBeta + 1/TanBeta)},
{T[Yu], \[Mu]/TanBeta Yu},
{T[Yd], \[Mu] TanBeta Yd},
{T[Ye], \[Mu] TanBeta Ye},
{T[Yu][3,3], (1/TanBeta + Xtt) Yu[3,3] Ms},
{MassB, M1Input},
{MassWB, M2Input},
{MassG, M3Input},
{mq2, mq2Input},
{mu2, mu2Input},
{md2, md2Input},
{ml2, ml2Input},
{me2, me2Input},
{\[Mu], MuInput},
{B[\[Mu]], mAInput^2/(TanBeta + 1/TanBeta)},
{T[Yu], AuInput Yu},
{T[Yd], AdInput Yd},
{T[Ye], AeInput Ye},
{vd, vu / TanBeta}
};

Expand All @@ -43,20 +57,19 @@ InitialGuessAtLowScale = {
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{MassB, Ms},
{MassWB, Ms},
{MassG, Ms},
{mq2, UNITMATRIX[3] Ms^2},
{mu2, UNITMATRIX[3] Ms^2},
{md2, UNITMATRIX[3] Ms^2},
{ml2, UNITMATRIX[3] Ms^2},
{me2, UNITMATRIX[3] Ms^2},
{\[Mu], Ms},
{B[\[Mu]], Sqr[Ms]/(TanBeta + 1./TanBeta)},
{T[Yu], MODELPARAMETER[Yu] Ms/TanBeta},
{T[Yd], MODELPARAMETER[Yd] Ms TanBeta},
{T[Ye], MODELPARAMETER[Ye] Ms TanBeta},
{T[Yu][3,3], (1/TanBeta + Xtt) Ms }
{MassB, M1Input},
{MassWB, M2Input},
{MassG, M3Input},
{mq2, mq2Input},
{mu2, mu2Input},
{md2, md2Input},
{ml2, ml2Input},
{me2, me2Input},
{\[Mu], MuInput},
{B[\[Mu]], mAInput^2/(TanBeta + 1/TanBeta)},
{T[Yu], AuInput Yu},
{T[Yd], AdInput Yd},
{T[Ye], AeInput Ye}
};

LowScale = LowEnergyConstant[MZ];
Expand Down
47 changes: 45 additions & 2 deletions model_files/MSSMMuBMu/LesHouches.in.MSSMMuBMu
View file
Expand Up @@ -17,6 +17,13 @@ Block FlexibleSUSY
13 1 # Top quark 2-loop corrections QCD
14 1.000000000e-11 # beta-function zero threshold
15 0 # calculate observables (a_muon, ...)
16 0 # force positive majorana masses
17 0 # pole mass renormalization scale (0 = SUSY scale)
18 0 # pole mass renormalization scale in the EFT (0 = min(SUSY scale, Mt))
19 0 # EFT matching scale (0 = SUSY scale)
20 2 # EFT loop order for upwards matching
21 1 # EFT loop order for downwards matching
22 0 # EFT index of SM-like Higgs in the BSM model
Block SMINPUTS # Standard Model inputs
1 1.279440000e+02 # alpha^(-1) SM MSbar(MZ)
2 1.166380000e-05 # G_Fermi
Expand All @@ -39,5 +46,41 @@ Block MINPAR # Input parameters
4 1 # SignMu
Block EXTPAR
0 2000 # Ms
14 0 # Xt / Ms
25 5 # tan(beta) at Ms
1 2000 # M1(MSUSY)
2 2000 # M2(MSUSY)
3 2000 # M3(MSUSY)
4 2000 # Mu(MSUSY)
5 2000 # mA(MSUSY)
25 5 # TanBeta(MSUSY)
Block MSQ2IN
1 1 4.00000000E+06 # mq2(1,1)
2 2 4.00000000E+06 # mq2(2,2)
3 3 4.00000000E+06 # mq2(3,3)
Block MSE2IN
1 1 4.00000000E+06 # me2(1,1)
2 2 4.00000000E+06 # me2(2,2)
3 3 4.00000000E+06 # me2(3,3)
Block MSL2IN
1 1 4.00000000E+06 # ml2(1,1)
2 2 4.00000000E+06 # ml2(2,2)
3 3 4.00000000E+06 # ml2(3,3)
Block MSU2IN
1 1 4.00000000E+06 # mu2(1,1)
2 2 4.00000000E+06 # mu2(2,2)
3 3 4.00000000E+06 # mu2(3,3)
Block MSD2IN
1 1 4.00000000E+06 # md2(1,1)
2 2 4.00000000E+06 # md2(2,2)
3 3 4.00000000E+06 # md2(3,3)
Block AUIN
1 1 400 # Au(1,1)
2 2 400 # Au(2,2)
3 3 400 # Au(3,3)
Block ADIN
1 1 1e4 # Ad(1,1)
2 2 1e4 # Ad(2,2)
3 3 1e4 # Ad(3,3)
Block AEIN
1 1 1e4 # Ad(1,1)
2 2 1e4 # Ad(2,2)
3 3 1e4 # Ad(3,3)
50 changes: 32 additions & 18 deletions model_files/MSSMtower/FlexibleSUSY.m.in
View file
Expand Up @@ -9,32 +9,46 @@ MINPAR = {
};

EXTPAR = {
{0, Ms},
{14, Xtt},
{0, MSUSY},
{1, M1Input},
{2, M2Input},
{3, M3Input},
{4, MuInput},
{5, mAInput},
{25, TanBeta}
};

FSExtraInputParameters = {
{mq2Input, MSQ2IN, {3,3}}, (* 3x3 matrix *)
{mu2Input, MSU2IN, {3,3}}, (* 3x3 matrix *)
{md2Input, MSD2IN, {3,3}}, (* 3x3 matrix *)
{ml2Input, MSL2IN, {3,3}}, (* 3x3 matrix *)
{me2Input, MSE2IN, {3,3}}, (* 3x3 matrix *)
{AuInput, AUIN, {3,3}}, (* 3x3 matrix *)
{AdInput, ADIN, {3,3}}, (* 3x3 matrix *)
{AeInput, AEIN, {3,3}} (* 3x3 matrix *)
};

EWSBOutputParameters = { mHd2, mHu2 };

SUSYScale = Ms;
SUSYScale = MSUSY;

SUSYScaleFirstGuess = Ms;
SUSYScaleFirstGuess = MSUSY;

SUSYScaleInput = {
{MassB, Ms},
{MassWB, Ms},
{MassG, Ms},
{mq2, UNITMATRIX[3] Ms^2},
{mu2, UNITMATRIX[3] Ms^2},
{md2, UNITMATRIX[3] Ms^2},
{ml2, UNITMATRIX[3] Ms^2},
{me2, UNITMATRIX[3] Ms^2},
{\[Mu], Ms},
{B[\[Mu]], Sqr[Ms]/(TanBeta + 1/TanBeta)},
{T[Yu], Ms/TanBeta Yu},
{T[Yd], Ms TanBeta Yd},
{T[Ye], Ms TanBeta Ye},
{T[Yu][3,3], (Ms/TanBeta + Xtt Ms) Yu[3,3]}
{MassB, M1Input},
{MassWB, M2Input},
{MassG, M3Input},
{mq2, mq2Input},
{mu2, mu2Input},
{md2, md2Input},
{ml2, ml2Input},
{me2, me2Input},
{\[Mu], MuInput},
{B[\[Mu]], mAInput^2/(TanBeta + 1/TanBeta)},
{T[Yu], AuInput Yu},
{T[Yd], AdInput Yd},
{T[Ye], AeInput Ye}
};

InitialGuessAtSUSYScale = SUSYScaleInput;
Expand Down
44 changes: 40 additions & 4 deletions model_files/MSSMtower/LesHouches.in.MSSMtower
View file
@@ -1,5 +1,5 @@
Block MODSEL # Select model
# 12 1000 # DRbar parameter output scale (GeV)
6 0 # flavour violation
Block FlexibleSUSY
0 1.000000000e-05 # precision goal
1 0 # max. iterations (0 = automatic)
Expand Down Expand Up @@ -27,7 +27,7 @@ Block FlexibleSUSY
Block SMINPUTS # Standard Model inputs
1 1.279440000e+02 # alpha^(-1) SM MSbar(MZ)
2 1.166380000e-05 # G_Fermi
3 1.185000000e-01 # alpha_s(MZ) SM MSbar
3 1.184000000e-01 # alpha_s(MZ) SM MSbar
4 9.118760000e+01 # MZ(pole)
5 4.180000000e+00 # mb(mb) SM MSbar
6 1.733400000e+02 # mtop(pole)
Expand All @@ -46,5 +46,41 @@ Block MINPAR # Input parameters
4 1 # SignMu
Block EXTPAR
0 2000 # Ms
14 0 # Xt / Ms
25 5 # tan(beta) at Ms
1 2000 # M1(MSUSY)
2 2000 # M2(MSUSY)
3 2000 # M3(MSUSY)
4 2000 # Mu(MSUSY)
5 2000 # mA(MSUSY)
25 5 # TanBeta(MSUSY)
Block MSQ2IN
1 1 4.00000000E+06 # mq2(1,1)
2 2 4.00000000E+06 # mq2(2,2)
3 3 4.00000000E+06 # mq2(3,3)
Block MSE2IN
1 1 4.00000000E+06 # me2(1,1)
2 2 4.00000000E+06 # me2(2,2)
3 3 4.00000000E+06 # me2(3,3)
Block MSL2IN
1 1 4.00000000E+06 # ml2(1,1)
2 2 4.00000000E+06 # ml2(2,2)
3 3 4.00000000E+06 # ml2(3,3)
Block MSU2IN
1 1 4.00000000E+06 # mu2(1,1)
2 2 4.00000000E+06 # mu2(2,2)
3 3 4.00000000E+06 # mu2(3,3)
Block MSD2IN
1 1 4.00000000E+06 # md2(1,1)
2 2 4.00000000E+06 # md2(2,2)
3 3 4.00000000E+06 # md2(3,3)
Block AUIN
1 1 400 # Au(1,1)
2 2 400 # Au(2,2)
3 3 400 # Au(3,3)
Block ADIN
1 1 1e4 # Ad(1,1)
2 2 1e4 # Ad(2,2)
3 3 1e4 # Ad(3,3)
Block AEIN
1 1 1e4 # Ad(1,1)
2 2 1e4 # Ad(2,2)
3 3 1e4 # Ad(3,3)
15 changes: 15 additions & 0 deletions test/test_MSSMtower.sh
View file
Expand Up @@ -113,6 +113,9 @@ run_sg() {
local SG="$1"
local MS2=$(echo "scale=5; ${MS}^2" | bc)
local At=$(echo "scale=10; (1./${TB} + ${Xt}) * ${MS}" | bc)
local Au=$(echo "scale=10; (1./${TB} + 0) * ${MS}" | bc)
local Ad=$(echo "scale=10; (${TB} + 0) * ${MS}" | bc)
local Ae="$Ad"
local slha_output=
local block=
local value=
Expand Down Expand Up @@ -154,6 +157,18 @@ Block MSD2IN
1 1 ${MS2} # md2(1,1)
2 2 ${MS2} # md2(2,2)
3 3 ${MS2} # md2(3,3)
Block AUIN
1 1 ${Au} # Au(1,1)
2 2 ${Au} # Au(2,2)
3 3 ${At} # Au(3,3)
Block ADIN
1 1 ${Ad} # Ad(1,1)
2 2 ${Ad} # Ad(2,2)
3 3 ${Ad} # Ad(3,3)
Block AEIN
1 1 ${Ae} # Ad(1,1)
2 2 ${Ae} # Ad(2,2)
3 3 ${Ae} # Ad(3,3)
EOF
})

Expand Down

0 comments on commit ebcc1e6

Please sign in to comment.