Skip to content

Commit

Permalink
make helper functions private
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Dec 21, 2016
1 parent 5b669f9 commit fa44905
Show file tree
Hide file tree
Showing 2 changed files with 19 additions and 8 deletions.
19 changes: 19 additions & 0 deletions legacy/lowe_legacy.cpp
Expand Up @@ -16,6 +16,25 @@ namespace softsusy {

namespace {

// Given a value of mt, and alphas(MZ), find alphas(mt) to 1 loops in qcd:
// it's a very good approximation at these scales, better than 10^-3 accuracy
double getAsmt(double mtop, double alphasMz) {
using std::log;
return alphasMz /
(1.0 - 23.0 * alphasMz / (6.0 * M_PI) * log(MZ / mtop));
}

// Input pole mass of top and alphaS(mt), outputs running mass mt(mt)
// including one-loop standard model correction only
double getRunMt(double poleMt, double asmt) {
return poleMt / (1.0 + (4.0 / (3.0 * M_PI)) * asmt);
}

// Given pole mass and alphaS(MZ), returns running top mass -- one loop qcd
double getRunMtFromMz(double poleMt, double asMZ) {
return getRunMt(poleMt, getAsmt(poleMt, asMZ));
}

/// external object temp used to get objects into external routines, however:
/// don't use it!
QedQcd_legacy *tempLe;
Expand Down
8 changes: 0 additions & 8 deletions legacy/lowe_legacy.h
Expand Up @@ -185,14 +185,6 @@ istream & operator >>(istream &left, QedQcd_legacy &m);
/// Call with fname "" if you want it to come from standard input
/// "massIn" is an example of a data initialisation file:
void readIn(QedQcd_legacy & oneset, const char fname[80]);
/// Input pole mass of top and alphaS(mt), outputs running mass mt(mt)
/// including one-loop standard model correction only
double getRunMt(double poleMt, double asmt);
/// Given a value of mt, and alphas(MZ), find alphas(mt) to 1 loops in qcd:
/// it's a very good approximation at these scales, better than 10^-3 accuracy
double getAsmt(double mtop, double alphasMz);
/// Given pole mass and alphaS(MZ), returns running top mass -- one loop qcd
double getRunMtFromMz(double poleMt, double asMZ);

inline QedQcd_legacy::QedQcd_legacy(const QedQcd_legacy &m)
: RGE()
Expand Down

0 comments on commit fa44905

Please sign in to comment.