Skip to content

Commit

Permalink
rename variable for clarity
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Aug 11, 2016
1 parent 9f95c69 commit 9deb24e
Show file tree
Hide file tree
Showing 2 changed files with 19 additions and 18 deletions.
26 changes: 13 additions & 13 deletions src/mssm_twoloophiggs.cpp
View file
Expand Up @@ -81,15 +81,15 @@ Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_ab_as_mssm(
}

Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_atau_atau_mssm(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2)
{
Eigen::Matrix<double, 2, 1> result;

LOCK_MUTEX();

tausqtad_(&mtau2, &mA2, &msnusq, &mstau12, &mstau22, &sintau,
tausqtad_(&mtau2, &mA2, &msv2, &mstau12, &mstau22, &sintau,
&costau, &scale2, &mu, &tanb, &vev2, &result(0), &result(1));

return -result;
Expand Down Expand Up @@ -147,15 +147,15 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm_with_tadpoles(
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm_with_tadpoles(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2, int scheme)
{
Eigen::Matrix<double, 2, 2> result;

LOCK_MUTEX();

tausqhiggs_(&mtau2, &mA2, &msnusq, &mstau12, &mstau22, &sintau,
tausqhiggs_(&mtau2, &mA2, &msv2, &mstau12, &mstau22, &sintau,
&costau, &scale2, &mu, &tanb, &vev2, &scheme,
&result(0,0), &result(1,1), &result(0,1));

Expand Down Expand Up @@ -206,15 +206,15 @@ double self_energy_pseudoscalar_2loop_ab_as_mssm_with_tadpoles(
}

double self_energy_pseudoscalar_2loop_atau_atau_mssm_with_tadpoles(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2)
{
double result;

LOCK_MUTEX();

tausqodd_(&mtau2, &mA2, &msnusq, &mstau12, &mstau22, &sintau,
tausqodd_(&mtau2, &mA2, &msv2, &mstau12, &mstau22, &sintau,
&costau, &scale2, &mu, &tanb, &vev2, &result);

return -result;
Expand Down Expand Up @@ -324,22 +324,22 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm(
}

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2, int scheme)
{
const Eigen::Matrix<double, 2, 2> result =
self_energy_higgs_2loop_atau_atau_mssm_with_tadpoles(
mtau2, mA2, msnusq, mstau12, mstau22, sintau, costau, scale2,
mtau2, mA2, msv2, mstau12, mstau22, sintau, costau, scale2,
mu, tanb, vev2, scheme);

const double dMA = self_energy_pseudoscalar_2loop_atau_atau_mssm_with_tadpoles(
mtau2, mA2, msnusq, mstau12, mstau22, sintau, costau, scale2,
mtau2, mA2, msv2, mstau12, mstau22, sintau, costau, scale2,
mu, tanb, vev2);

const Eigen::Matrix<double, 2, 1> tadpoles =
tadpole_higgs_2loop_atau_atau_mssm(
mtau2, mA2, msnusq, mstau12, mstau22, sintau, costau, scale2,
mtau2, mA2, msv2, mstau12, mstau22, sintau, costau, scale2,
mu, tanb, vev2);

const Eigen::Matrix<double, 2, 2> tM =
Expand Down Expand Up @@ -422,17 +422,17 @@ Eigen::Matrix<double, 2, 2> self_energy_pseudoscalar_2loop_ab_as_mssm(
}

Eigen::Matrix<double, 2, 2> self_energy_pseudoscalar_2loop_atau_atau_mssm(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2)
{
const double se = self_energy_pseudoscalar_2loop_atau_atau_mssm_with_tadpoles(
mtau2, mA2, msnusq, mstau12, mstau22, sintau, costau, scale2,
mtau2, mA2, msv2, mstau12, mstau22, sintau, costau, scale2,
mu, tanb, vev2);

const Eigen::Matrix<double, 2, 1> tadpoles =
tadpole_higgs_2loop_atau_atau_mssm(
mtau2, mA2, msnusq, mstau12, mstau22, sintau, costau, scale2,
mtau2, mA2, msv2, mstau12, mstau22, sintau, costau, scale2,
mu, tanb, vev2);

return subtract_mssm_tadpoles_pseudoscalar(se, tadpoles, tanb);
Expand Down
11 changes: 6 additions & 5 deletions src/mssm_twoloophiggs.hpp
View file
Expand Up @@ -39,6 +39,7 @@
* msb22 : squared DR-bar heaviest sbottom mass
* mstau12: squared DR-bar lightest stau mass
* mstau22: squared DR-bar heaviest stau mass
* msv2 : squared DR-bar tau sneutrino mass
*
* sxt : sine of DR-bar stop mixing angle in the MSSM
* cxt : cosine of DR-bar stop mixing angle in the MSSM
Expand Down Expand Up @@ -78,7 +79,7 @@ Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_ab_as_mssm(
double mu, double cotb, double vev2, double gs);

Eigen::Matrix<double, 2, 1> tadpole_higgs_2loop_atau_atau_mssm(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2);

Expand All @@ -101,7 +102,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm(
double cotb, double vev2, double gs, int scheme = 0);

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2, int scheme = 0);

Expand All @@ -123,7 +124,7 @@ Eigen::Matrix<double, 2, 2> self_energy_pseudoscalar_2loop_ab_as_mssm(
double cotb, double vev2, double gs);

Eigen::Matrix<double, 2, 2> self_energy_pseudoscalar_2loop_atau_atau_mssm(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2);

Expand All @@ -146,7 +147,7 @@ Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_ab_as_mssm_with_tadpoles(
double cotb, double vev2, double gs, int scheme = 0);

Eigen::Matrix<double, 2, 2> self_energy_higgs_2loop_atau_atau_mssm_with_tadpoles(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2, int scheme = 0);

Expand All @@ -168,7 +169,7 @@ double self_energy_pseudoscalar_2loop_ab_as_mssm_with_tadpoles(
double cotb, double vev2, double gs);

double self_energy_pseudoscalar_2loop_atau_atau_mssm_with_tadpoles(
double mtau2, double mA2, double msnusq, double mstau12,
double mtau2, double mA2, double msv2, double mstau12,
double mstau22, double sintau, double costau, double scale2,
double mu, double tanb, double vev2);

Expand Down

0 comments on commit 9deb24e

Please sign in to comment.