Skip to content

Commit

Permalink
rename symbols referring to the fermion DR-bar masses in the model
Browse files Browse the repository at this point in the history 
The new names have been proposed by Peter to emphasize that they are
of a 3x3 matrix type.
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Jul 7, 2016
1 parent 9fa11d3 commit 857c51b
Show file tree
Hide file tree
Showing 9 changed files with 118 additions and 106 deletions.
6 changes: 3 additions & 3 deletions meta/Constraint.m
View file
Expand Up @@ -46,17 +46,17 @@
];

ApplyConstraint[{parameter_ | parameter_[__] /; parameter === SARAH`UpYukawa,
value_ /; (!FreeQ[value, Global`topDRbar] || value === Automatic)},
value_ /; (!FreeQ[value, Global`upQuarksDRbar] || value === Automatic)},
modelName_String] :=
"calculate_" <> CConversion`ToValidCSymbolString[parameter] <> "_DRbar();\n";

ApplyConstraint[{parameter_ | parameter_[__] /; parameter === SARAH`DownYukawa,
value_ /; (!FreeQ[value, Global`bottomDRbar] || value === Automatic)},
value_ /; (!FreeQ[value, Global`downQuarksDRbar] || value === Automatic)},
modelName_String] :=
"calculate_" <> CConversion`ToValidCSymbolString[parameter] <> "_DRbar();\n";

ApplyConstraint[{parameter_ | parameter_[__] /; parameter === SARAH`ElectronYukawa,
value_ /; (!FreeQ[value, Global`electronDRbar] || value === Automatic)},
value_ /; (!FreeQ[value, Global`downLeptonsDRbar] || value === Automatic)},
modelName_String] :=
"calculate_" <> CConversion`ToValidCSymbolString[parameter] <> "_DRbar();\n";

Expand Down
36 changes: 24 additions & 12 deletions meta/FlexibleSUSY.m
View file
Expand Up @@ -642,15 +642,15 @@ FlexibleSUSY model file (FlexibleSUSY.m).
ThresholdCorrections`SetDRbarYukawaCouplingElectron[settings]
};
calculateDRbarMasses = {
LoopMasses`CallCalculateDRbarMass["Up Quark" , "Up-Quarks" , 1, "topDRbar" , "qedqcd.displayMass(softsusy::mUp)" ],
LoopMasses`CallCalculateDRbarMass["Charmed Quark" , "Up-Quarks" , 2, "topDRbar" , "qedqcd.displayMass(softsusy::mCharm)" ],
LoopMasses`CallCalculateDRbarMass["Top Quark" , "Up-Quarks" , 3, "topDRbar" , "qedqcd.displayPoleMt()" ],
LoopMasses`CallCalculateDRbarMass["Down Quark" , "Down-Quarks", 1, "bottomDRbar" , "qedqcd.displayMass(softsusy::mDown)" ],
LoopMasses`CallCalculateDRbarMass["Strange Quark" , "Down-Quarks", 2, "bottomDRbar" , "qedqcd.displayMass(softsusy::mStrange)" ],
LoopMasses`CallCalculateDRbarMass["Bottom Quark" , "Down-Quarks", 3, "bottomDRbar" , "qedqcd.displayMass(softsusy::mBottom)" ],
LoopMasses`CallCalculateDRbarMass["Electron" , "Leptons" , 1, "electronDRbar", "qedqcd.displayMass(softsusy::mElectron)"],
LoopMasses`CallCalculateDRbarMass["Muon" , "Leptons" , 2, "electronDRbar", "qedqcd.displayMass(softsusy::mMuon)" ],
LoopMasses`CallCalculateDRbarMass["Tau" , "Leptons" , 3, "electronDRbar", "qedqcd.displayMass(softsusy::mTau)" ],
LoopMasses`CallCalculateDRbarMass["Up Quark" , "Up-Quarks" , 1, "upQuarksDRbar", "qedqcd.displayMass(softsusy::mUp)" ],
LoopMasses`CallCalculateDRbarMass["Charmed Quark" , "Up-Quarks" , 2, "upQuarksDRbar", "qedqcd.displayMass(softsusy::mCharm)" ],
LoopMasses`CallCalculateDRbarMass["Top Quark" , "Up-Quarks" , 3, "upQuarksDRbar", "qedqcd.displayPoleMt()" ],
LoopMasses`CallCalculateDRbarMass["Down Quark" , "Down-Quarks", 1, "downQuarksDRbar", "qedqcd.displayMass(softsusy::mDown)" ],
LoopMasses`CallCalculateDRbarMass["Strange Quark" , "Down-Quarks", 2, "downQuarksDRbar", "qedqcd.displayMass(softsusy::mStrange)" ],
LoopMasses`CallCalculateDRbarMass["Bottom Quark" , "Down-Quarks", 3, "downQuarksDRbar", "qedqcd.displayMass(softsusy::mBottom)" ],
LoopMasses`CallCalculateDRbarMass["Electron" , "Leptons" , 1, "downLeptonsDRbar", "qedqcd.displayMass(softsusy::mElectron)"],
LoopMasses`CallCalculateDRbarMass["Muon" , "Leptons" , 2, "downLeptonsDRbar", "qedqcd.displayMass(softsusy::mMuon)" ],
LoopMasses`CallCalculateDRbarMass["Tau" , "Leptons" , 3, "downLeptonsDRbar", "qedqcd.displayMass(softsusy::mTau)" ],
LoopMasses`CallCalculateDRbarMass["Electron Neutrino", "Neutrinos" , 1, "neutrinoDRbar", "qedqcd.displayNeutrinoPoleMass(1)" ],
LoopMasses`CallCalculateDRbarMass["Muon Neutrino" , "Neutrinos" , 2, "neutrinoDRbar", "qedqcd.displayNeutrinoPoleMass(2)" ],
LoopMasses`CallCalculateDRbarMass["Tau Neutrino" , "Neutrinos" , 3, "neutrinoDRbar", "qedqcd.displayNeutrinoPoleMass(3)" ]
Expand Down Expand Up @@ -1844,6 +1844,18 @@ corresponding tadpole is real or imaginary (only in models with CP
susyBetaFunctions = DeleteBuggyBetaFunctions /@ susyBetaFunctions;
susyBreakingBetaFunctions = DeleteBuggyBetaFunctions /@ susyBreakingBetaFunctions;

(* backwards compatibility replacements in constraints *)
backwardsCompatRules = {
Global`topDRbar -> Global`upQuarksDRbar,
Global`bottomDRbar -> Global`downQuarksDRbar,
Global`electronDRbar -> Global`downLeptonsDRbar
};
FlexibleSUSY`LowScaleInput = FlexibleSUSY`LowScaleInput /. backwardsCompatRules;
FlexibleSUSY`SUSYScaleInput = FlexibleSUSY`SUSYScaleInput /. backwardsCompatRules;
FlexibleSUSY`HighScaleInput = FlexibleSUSY`HighScaleInput /. backwardsCompatRules;
FlexibleSUSY`InitialGuessAtLowScale = FlexibleSUSY`InitialGuessAtLowScale /. backwardsCompatRules;
FlexibleSUSY`InitialGuessAtHighScale = FlexibleSUSY`InitialGuessAtHighScale /. backwardsCompatRules;

(* store all model parameters *)
allParameters = ((#[[1]])& /@ Join[Join @@ susyBetaFunctions, Join @@ susyBreakingBetaFunctions]) /.
Parameters`StripSARAHIndicesRules[1] /.
Expand All @@ -1852,9 +1864,9 @@ corresponding tadpole is real or imaginary (only in models with CP
Parameters`StripSARAHIndicesRules[4];
allIndexReplacementRules = Join[
Parameters`CreateIndexReplacementRules[allParameters],
{Global`topDRbar[i_,j_] :> Global`topDRbar[i-1,j-1],
Global`bottomDRbar[i_,j_] :> Global`bottomDRbar[i-1,j-1],
Global`electronDRbar[i_,j_] :> Global`electronDRbar[i-1,j-1]}
{Global`upQuarksDRbar[i_,j_] :> Global`upQuarksDRbar[i-1,j-1],
Global`downQuarksDRbar[i_,j_] :> Global`downQuarksDRbar[i-1,j-1],
Global`downLeptonsDRbar[i_,j_] :> Global`downLeptonsDRbar[i-1,j-1]}
];
Parameters`SetModelParameters[allParameters];
DebugPrint["model parameters: ", allParameters];
Expand Down
6 changes: 3 additions & 3 deletions meta/ThresholdCorrections.m
View file
Expand Up @@ -235,13 +235,13 @@
];

SetDRbarYukawaCouplingTop[settings_] :=
SetDRbarYukawaCouplingFermion[SARAH`TopQuark, SARAH`UpYukawa, Global`topDRbar, settings];
SetDRbarYukawaCouplingFermion[SARAH`TopQuark, SARAH`UpYukawa, Global`upQuarksDRbar, settings];

SetDRbarYukawaCouplingBottom[settings_] :=
SetDRbarYukawaCouplingFermion[SARAH`BottomQuark, SARAH`DownYukawa, Global`bottomDRbar, settings];
SetDRbarYukawaCouplingFermion[SARAH`BottomQuark, SARAH`DownYukawa, Global`downQuarksDRbar, settings];

SetDRbarYukawaCouplingElectron[settings_] :=
SetDRbarYukawaCouplingFermion[SARAH`Electron, SARAH`ElectronYukawa, Global`electronDRbar, settings];
SetDRbarYukawaCouplingFermion[SARAH`Electron, SARAH`ElectronYukawa, Global`downLeptonsDRbar, settings];

SetDRbarYukawaCouplingFermionMatrix[fermion_, yukawa_, mass_, setting_] :=
Module[{y, f = setting},
Expand Down
6 changes: 3 additions & 3 deletions model_files/CMSSMCKM/FlexibleSUSY.m.in
View file
Expand Up @@ -50,9 +50,9 @@ LowScale = LowEnergyConstant[MZ];
LowScaleFirstGuess = LowEnergyConstant[MZ];

LowScaleInput = {
{Yu, Sqrt[2] MatMul[Tp[CKM], topDRbar] / vu},
{Yd, Sqrt[2] bottomDRbar / vd},
{Ye, Sqrt[2] electronDRbar / vd},
{Yu, Sqrt[2] MatMul[Tp[CKM], upQuarksDRbar] / vu},
{Yd, Sqrt[2] downQuarksDRbar / vd},
{Ye, Sqrt[2] downLeptonsDRbar / vd},
{vd, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
{vu, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
};
Expand Down
24 changes: 12 additions & 12 deletions model_files/U1xMSSM3G/FlexibleSUSY.m.in
View file
Expand Up @@ -138,13 +138,13 @@ LowScaleInput = {
{Temporary[gX1], 0},
{vd, vDRbar Cos[ArcTan[TanBeta]]},
{vu, vDRbar Sin[ArcTan[TanBeta]]},
{Yu[1,1], Sqrt[2] topDRbar[1,1] / vu},
{Yu[2,2], Sqrt[2] topDRbar[2,2] / vu},
{Yu[3,3], Sqrt[2] topDRbar[3,3] / vu},
{Yu[1,1], Sqrt[2] upQuarksDRbar[1,1] / vu},
{Yu[2,2], Sqrt[2] upQuarksDRbar[2,2] / vu},
{Yu[3,3], Sqrt[2] upQuarksDRbar[3,3] / vu},
{Yd, Automatic},
{Ye[1,1], Sqrt[2] electronDRbar[1,1] / vd},
{Ye[2,2], Sqrt[2] electronDRbar[2,2] / vd},
{Ye[3,3], Sqrt[2] electronDRbar[3,3] / vd}
{Ye[1,1], Sqrt[2] downLeptonsDRbar[1,1] / vd},
{Ye[2,2], Sqrt[2] downLeptonsDRbar[2,2] / vd},
{Ye[3,3], Sqrt[2] downLeptonsDRbar[3,3] / vd}
};


Expand All @@ -164,13 +164,13 @@ InitialGuessAtLowScale = {
{L[lw], 0},
{vd, LowEnergyConstant[vev] Cos[ArcTan[TanBeta]]},
{vu, LowEnergyConstant[vev] Sin[ArcTan[TanBeta]]},
{Yu[1,1], Sqrt[2] topDRbar[1,1] / vu},
{Yu[2,2], Sqrt[2] topDRbar[2,2] / vu},
{Yu[3,3], Sqrt[2] topDRbar[3,3] / vu},
{Yu[1,1], Sqrt[2] upQuarksDRbar[1,1] / vu},
{Yu[2,2], Sqrt[2] upQuarksDRbar[2,2] / vu},
{Yu[3,3], Sqrt[2] upQuarksDRbar[3,3] / vu},
{Yd, Automatic},
{Ye[1,1], Sqrt[2] electronDRbar[1,1] / vd},
{Ye[2,2], Sqrt[2] electronDRbar[2,2] / vd},
{Ye[3,3], Sqrt[2] electronDRbar[3,3] / vd}
{Ye[1,1], Sqrt[2] downLeptonsDRbar[1,1] / vd},
{Ye[2,2], Sqrt[2] downLeptonsDRbar[2,2] / vd},
{Ye[3,3], Sqrt[2] downLeptonsDRbar[3,3] / vd}
};


Expand Down
74 changes: 37 additions & 37 deletions src/standard_model.cpp
View file
Expand Up @@ -927,23 +927,23 @@ void Standard_model::initial_guess_for_parameters()

v = Re(Electroweak_constants::vev);

Eigen::Matrix<double,3,3> topDRbar(Eigen::Matrix<double,3,3>::Zero());
topDRbar(0,0) = mu_guess;
topDRbar(1,1) = mc_guess;
topDRbar(2,2) = mt_guess;
Yu = -((1.4142135623730951*topDRbar) / v).transpose();

Eigen::Matrix<double,3,3> bottomDRbar(Eigen::Matrix<double,3,3>::Zero());
bottomDRbar(0,0) = md_guess;
bottomDRbar(1,1) = ms_guess;
bottomDRbar(2,2) = mb_guess;
Yd = ((1.4142135623730951*bottomDRbar)/v).transpose();

Eigen::Matrix<double,3,3> electronDRbar(Eigen::Matrix<double,3,3>::Zero());
electronDRbar(0,0) = me_guess;
electronDRbar(1,1) = mm_guess;
electronDRbar(2,2) = mtau_guess;
Ye = ((1.4142135623730951*electronDRbar)/v).transpose();
Eigen::Matrix<double,3,3> upQuarksDRbar(Eigen::Matrix<double,3,3>::Zero());
upQuarksDRbar(0,0) = mu_guess;
upQuarksDRbar(1,1) = mc_guess;
upQuarksDRbar(2,2) = mt_guess;
Yu = -((1.4142135623730951*upQuarksDRbar) / v).transpose();

Eigen::Matrix<double,3,3> downQuarksDRbar(Eigen::Matrix<double,3,3>::Zero());
downQuarksDRbar(0,0) = md_guess;
downQuarksDRbar(1,1) = ms_guess;
downQuarksDRbar(2,2) = mb_guess;
Yd = ((1.4142135623730951*downQuarksDRbar)/v).transpose();

Eigen::Matrix<double,3,3> downLeptonsDRbar(Eigen::Matrix<double,3,3>::Zero());
downLeptonsDRbar(0,0) = me_guess;
downLeptonsDRbar(1,1) = mm_guess;
downLeptonsDRbar(2,2) = mtau_guess;
Ye = ((1.4142135623730951*downLeptonsDRbar)/v).transpose();

Lambdax = Sqr(MH) / Sqr(v);

Expand Down Expand Up @@ -1045,46 +1045,46 @@ double Standard_model::calculate_theta_w(double alpha_em_drbar)

void Standard_model::calculate_Yu_DRbar()
{
Eigen::Matrix<double,3,3> topDRbar(Eigen::Matrix<double,3,3>::Zero());
topDRbar(0,0) = qedqcd.displayMass(softsusy::mUp);
topDRbar(1,1) = qedqcd.displayMass(softsusy::mCharm);
topDRbar(2,2) = qedqcd.displayPoleMt();
Eigen::Matrix<double,3,3> upQuarksDRbar(Eigen::Matrix<double,3,3>::Zero());
upQuarksDRbar(0,0) = qedqcd.displayMass(softsusy::mUp);
upQuarksDRbar(1,1) = qedqcd.displayMass(softsusy::mCharm);
upQuarksDRbar(2,2) = qedqcd.displayPoleMt();

if (get_thresholds()) {
topDRbar(2,2) = calculate_MFu_DRbar(qedqcd.displayPoleMt(), 2);
upQuarksDRbar(2,2) = calculate_MFu_DRbar(qedqcd.displayPoleMt(), 2);
}

Yu = -((1.4142135623730951*topDRbar)/v).transpose();
Yu = -((1.4142135623730951*upQuarksDRbar)/v).transpose();
}

void Standard_model::calculate_Yd_DRbar()
{
Eigen::Matrix<double,3,3> bottomDRbar(Eigen::Matrix<double,3,3>::Zero());
bottomDRbar(0,0) = qedqcd.displayMass(softsusy::mDown);
bottomDRbar(1,1) = qedqcd.displayMass(softsusy::mStrange);
bottomDRbar(2,2) = qedqcd.displayMass(softsusy::mBottom);
Eigen::Matrix<double,3,3> downQuarksDRbar(Eigen::Matrix<double,3,3>::Zero());
downQuarksDRbar(0,0) = qedqcd.displayMass(softsusy::mDown);
downQuarksDRbar(1,1) = qedqcd.displayMass(softsusy::mStrange);
downQuarksDRbar(2,2) = qedqcd.displayMass(softsusy::mBottom);

if (get_thresholds()) {
bottomDRbar(2,2) = calculate_MFd_DRbar(qedqcd.displayMass(softsusy::mBottom), 2);
downQuarksDRbar(2,2) = calculate_MFd_DRbar(qedqcd.displayMass(softsusy::mBottom), 2);
}

Yd = ((1.4142135623730951*bottomDRbar)/v).transpose();
Yd = ((1.4142135623730951*downQuarksDRbar)/v).transpose();
}

void Standard_model::calculate_Ye_DRbar()
{
Eigen::Matrix<double,3,3> electronDRbar(Eigen::Matrix<double,3,3>::Zero());
electronDRbar(0,0) = qedqcd.displayPoleMel();
electronDRbar(1,1) = qedqcd.displayPoleMmuon();
electronDRbar(2,2) = qedqcd.displayPoleMtau();
Eigen::Matrix<double,3,3> downLeptonsDRbar(Eigen::Matrix<double,3,3>::Zero());
downLeptonsDRbar(0,0) = qedqcd.displayPoleMel();
downLeptonsDRbar(1,1) = qedqcd.displayPoleMmuon();
downLeptonsDRbar(2,2) = qedqcd.displayPoleMtau();

if (get_thresholds()) {
electronDRbar(0,0) = calculate_MFe_DRbar(qedqcd.displayMass(softsusy::mElectron), 0);
electronDRbar(1,1) = calculate_MFe_DRbar(qedqcd.displayMass(softsusy::mMuon), 1);
electronDRbar(2,2) = calculate_MFe_DRbar(qedqcd.displayMass(softsusy::mTau), 2);
downLeptonsDRbar(0,0) = calculate_MFe_DRbar(qedqcd.displayMass(softsusy::mElectron), 0);
downLeptonsDRbar(1,1) = calculate_MFe_DRbar(qedqcd.displayMass(softsusy::mMuon), 1);
downLeptonsDRbar(2,2) = calculate_MFe_DRbar(qedqcd.displayMass(softsusy::mTau), 2);
}

Ye = ((1.4142135623730951*electronDRbar)/v).transpose();
Ye = ((1.4142135623730951*downLeptonsDRbar)/v).transpose();

}

Expand Down
24 changes: 12 additions & 12 deletions templates/two_scale_initial_guesser.cpp.in
View file
Expand Up @@ -141,10 +141,10 @@ void @ModelName@_initial_guesser<Two_scale>::calculate_DRbar_yukawa_couplings()
*/
void @ModelName@_initial_guesser<Two_scale>::calculate_@UpYukawa@_DRbar()
{
Eigen::Matrix<std::complex<double>,3,3> topDRbar(ZEROMATRIXCOMPLEX(3,3));
topDRbar(0,0) = mu_guess;
topDRbar(1,1) = mc_guess;
topDRbar(2,2) = mt_guess;
Eigen::Matrix<std::complex<double>,3,3> upQuarksDRbar(ZEROMATRIXCOMPLEX(3,3));
upQuarksDRbar(0,0) = mu_guess;
upQuarksDRbar(1,1) = mc_guess;
upQuarksDRbar(2,2) = mt_guess;

@setDRbarUpQuarkYukawaCouplings@
}
Expand All @@ -156,10 +156,10 @@ void @ModelName@_initial_guesser<Two_scale>::calculate_@UpYukawa@_DRbar()
*/
void @ModelName@_initial_guesser<Two_scale>::calculate_@DownYukawa@_DRbar()
{
Eigen::Matrix<std::complex<double>,3,3> bottomDRbar(ZEROMATRIXCOMPLEX(3,3));
bottomDRbar(0,0) = md_guess;
bottomDRbar(1,1) = ms_guess;
bottomDRbar(2,2) = mb_guess;
Eigen::Matrix<std::complex<double>,3,3> downQuarksDRbar(ZEROMATRIXCOMPLEX(3,3));
downQuarksDRbar(0,0) = md_guess;
downQuarksDRbar(1,1) = ms_guess;
downQuarksDRbar(2,2) = mb_guess;

@setDRbarDownQuarkYukawaCouplings@
}
Expand All @@ -171,10 +171,10 @@ void @ModelName@_initial_guesser<Two_scale>::calculate_@DownYukawa@_DRbar()
*/
void @ModelName@_initial_guesser<Two_scale>::calculate_@ElectronYukawa@_DRbar()
{
Eigen::Matrix<std::complex<double>,3,3> electronDRbar(ZEROMATRIXCOMPLEX(3,3));
electronDRbar(0,0) = me_guess;
electronDRbar(1,1) = mm_guess;
electronDRbar(2,2) = mtau_guess;
Eigen::Matrix<std::complex<double>,3,3> downLeptonsDRbar(ZEROMATRIXCOMPLEX(3,3));
downLeptonsDRbar(0,0) = me_guess;
downLeptonsDRbar(1,1) = mm_guess;
downLeptonsDRbar(2,2) = mtau_guess;

@setDRbarElectronYukawaCouplings@
}
Expand Down
24 changes: 12 additions & 12 deletions templates/two_scale_initial_guesser_low_scale_model.cpp.in
View file
Expand Up @@ -138,10 +138,10 @@ void @ModelName@_initial_guesser<Two_scale>::calculate_DRbar_yukawa_couplings()
*/
void @ModelName@_initial_guesser<Two_scale>::calculate_@UpYukawa@_DRbar()
{
Eigen::Matrix<std::complex<double>,3,3> topDRbar(ZEROMATRIXCOMPLEX(3,3));
topDRbar(0,0) = mu_guess;
topDRbar(1,1) = mc_guess;
topDRbar(2,2) = mt_guess;
Eigen::Matrix<std::complex<double>,3,3> upQuarksDRbar(ZEROMATRIXCOMPLEX(3,3));
upQuarksDRbar(0,0) = mu_guess;
upQuarksDRbar(1,1) = mc_guess;
upQuarksDRbar(2,2) = mt_guess;

@setDRbarUpQuarkYukawaCouplings@
}
Expand All @@ -153,10 +153,10 @@ void @ModelName@_initial_guesser<Two_scale>::calculate_@UpYukawa@_DRbar()
*/
void @ModelName@_initial_guesser<Two_scale>::calculate_@DownYukawa@_DRbar()
{
Eigen::Matrix<std::complex<double>,3,3> bottomDRbar(ZEROMATRIXCOMPLEX(3,3));
bottomDRbar(0,0) = md_guess;
bottomDRbar(1,1) = ms_guess;
bottomDRbar(2,2) = mb_guess;
Eigen::Matrix<std::complex<double>,3,3> downQuarksDRbar(ZEROMATRIXCOMPLEX(3,3));
downQuarksDRbar(0,0) = md_guess;
downQuarksDRbar(1,1) = ms_guess;
downQuarksDRbar(2,2) = mb_guess;

@setDRbarDownQuarkYukawaCouplings@
}
Expand All @@ -168,10 +168,10 @@ void @ModelName@_initial_guesser<Two_scale>::calculate_@DownYukawa@_DRbar()
*/
void @ModelName@_initial_guesser<Two_scale>::calculate_@ElectronYukawa@_DRbar()
{
Eigen::Matrix<std::complex<double>,3,3> electronDRbar(ZEROMATRIXCOMPLEX(3,3));
electronDRbar(0,0) = me_guess;
electronDRbar(1,1) = mm_guess;
electronDRbar(2,2) = mtau_guess;
Eigen::Matrix<std::complex<double>,3,3> downLeptonsDRbar(ZEROMATRIXCOMPLEX(3,3));
downLeptonsDRbar(0,0) = me_guess;
downLeptonsDRbar(1,1) = mm_guess;
downLeptonsDRbar(2,2) = mtau_guess;

@setDRbarElectronYukawaCouplings@
}
Expand Down

0 comments on commit 857c51b

Please sign in to comment.