diff --git a/CRunDec3/CRunDec.3.0.cpp b/CRunDec3/CRunDec.3.0.cpp
deleted file mode 100644
index 9113734..0000000
--- a/CRunDec3/CRunDec.3.0.cpp
+++ /dev/null
@@ -1,3086 +0,0 @@
-/*
- CRunDec.cpp, v3.0
-
- Authors: Florian Herren and Matthias Steinhauser (Feb 2017)
-
- F. Herren, M. Steinhauser,
- "Version 3 of RunDec and CRunDec"
- arXiv:1703.03751, TTP17-011
-
- Changes since v2.1:
- Dec 2016: added five-loop beta function (see arXiv:1606.08659);
- extended LamExpl, LamImpl, AlphasLam to five loops;
- ported DecAsDownMS, DecMqDownMS and DecMq{Up,Down}SI from RunDec;
- extended Dec{As,Mq}{Up,Down}{MS,SI} to four loops;
- Jan 2017: fixed Segmentation fault in case null pointer is passed as "mq" in
- mass conversion routines, minor bugfixes in threshold mass routines
- added DecLambdaUp/Down to 5 loops
- extended mRGI2mMS/mMS2mRGI/mMS2mRGImod to 5 loops
- included light quark mass effects in MS-OS relation up to three loops (arXiv:0708.1729)
-*/
-
-/*
- CRunDec.cpp, v2.1
-
- Authors: Florian Herren and Matthias Steinhauser (Jun 2016)
-
- Changes since v2.0:
- Jun 2016: improved 4-loop coefficients in mMS2mOS, mOS2mSI, mOS2mMS
-*/
-
-/*
- Based on:
-
- CRunDec.cpp, v2.0
-
- Authors: Florian Herren and Matthias Steinhauser (Jan 2016)
-
- Changes since v1.1:
- Sep 2015: extended mOS2mMS, mMS2mOS to four loops (see arXiv:1502.01030);
- mMS2mMS to five loops (see arXiv:1402.6611);
- DecMqUp/DownOS to four loops (see hep-ph/0512058, hep-ph/0512060
- and arXiv:1502.04719);
- prepared AsmMSrunexact, AlphasExact for five loops
- Oct 2015: implemented PS, 1S and RS masses, now requires C++11 for tgamma,
- extended mOS2mSI to four loops (see arXiv:1502.01030)
- Dec 2015: bugfix: nl -> nlq in fas6to5ms, ported DecAsDownSI/DecAsUpSI from RunDec
- Jan 2016: Add RS' mass relations (see arXiv:1502.01030)
-*/
-
-/*
- Based on:
-
- CRunDec.cpp, v1.1
-
- Author: Barbara Schmidt
-
- For documentation see
-
- CRunDec: a C++ package for running and decoupling of the
- strong coupling and quark masses
-
- by
-
- Barbara Schmidt, Matthias Steinhauser
-
- Comput.Phys.Commun. 183 (2012) 1845-1848
- arXiv:1201.6149
- SFB/CPP-12-03
- TTP12-02
-
- See also:
- K.G. Chetyrkin, J.H. Kuhn and M. Steinhauser,
- "RunDec: A Mathematica package for running and decoupling of the strong
- coupling and quark masses"
- Comput. Phys. Commun. 133 (2000) 43,
- arXiv:hep-ph/0004189
-
- May 2013: minor modifications to avoid "-Wunused-parameter"
- during compilation with g++ -Wall -Wextra
- Sep 2013: bug fix in fRungeKuttaImpl (thanks to Stephen Jones)
- Oct 2014: bug fix in AsmMSrunexact
-*/
-
-/*
-License:
-
-CRunDec is free software: you can redistribute it and/or modify
-it under the terms of the GNU General Public License as published by
-the Free Software Foundation, either version 3 of the License, or
-(at your option) any later version.
-
-This program is distributed in the hope that it will be useful,
-but WITHOUT ANY WARRANTY; without even the implied warranty of
-MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-GNU General Public License for more details.
-
-You should have received a copy of the GNU General Public License
-along with this program. If not, see .
-*/
-
-#include
-#include
-#include
-
-#include "CRunDec.h"
-
-// Some constants:
-#define cf 4./3.
-#define ca 3.
-#define tr 1./2.
-#define B4 -1.762800087073770864061897634679818807215137274389016762629478603776
-#define A4 0.5174790616738993863307581618988629456223774751413792582443193479770
-#define A5 0.5084005792422687074591088492585899413195411256648216487244977963526
-#define Pi 3.14159265358979323846
-#define Zeta2 (Pi*Pi)/6.
-#define Zeta3 1.20205690315959428539973816151144999076498629234049888179227155534
-#define Zeta4 (Pi*Pi*Pi*Pi)/90.
-#define Zeta5 1.03692775514336992633136548645703416805708091950191281197419267790
-#define Zeta6 (Pi*Pi*Pi*Pi*Pi*Pi)/945.
-#define Zeta7 1.00834927738192282683979754984979675959986356056523870641728313657160147831735573534609696891385132
-#define ln2 log(2)
-
-using namespace std;
-
-// Default constructor:
-CRunDec::CRunDec(){
- for(int i=0; i<4; i++){
- mq[i].first = 0.;
- mq[i].second = 0.;
- nfMmu[i].Mth=0.;
- nfMmu[i].muth=0.;
- nfMmu[i].nf=0;
- }
- AM.Asexact=0.;
- AM.mMSexact=0.;
-}
-
-// Constructor called with number for active flavours nf
-CRunDec::CRunDec(int n){
- this->SetConstants(n);
- for(int i=0; i<4; i++) {
- mq[i].first = 0.;
- mq[i].second = 0.;
- nfMmu[i].Mth=0.;
- nfMmu[i].muth=0.;
- nfMmu[i].nf=0;
- }
- AM.Asexact=0.;
- AM.mMSexact=0.;
-}
-
-// Define constants (if not already done in constructor)
-void CRunDec::SetConstants(int n){
- double nf=(double) n;
- Nf=nf;
- double nl = nf-1;
- Beta[0]= (double)0.25*(11.0 - 2.*nf/3.0);
- Beta[1]= (double)(102.0 - 38.*nf/3.0)/16.0;
- Beta[2]= (double)(0.5*2857. - 5033.*nf/18.0 + 325.*nf*nf/54.0)/64.0;
- Beta[3]= (double)(149753./6.0 + 3564.*Zeta3 +
- (-1078361./162.0 - 6508.*Zeta3/27.0)*nf +
- (50065./162.0 + 6472.*Zeta3/81.0)*nf*nf +
- 1093.*nf*nf*nf/729.0)/256.0;
- Beta[4]= (double)(8157455.0/16.0 + 621885.0/2.0*Zeta3 - 88209.0/2.0*Zeta4 - 288090.0*Zeta5 +
- (-336460813./1944. - 4811164./81.*Zeta3 + 33935./6.*Zeta4 + 1358995./27.*Zeta5)*nf +
- (25960913./1944. + 698531./81.*Zeta3 - 10526./9.*Zeta4 - 381760./81.*Zeta5)*nf*nf +
- (-630559./5832. - 48722./243.*Zeta3 + 1618./27.*Zeta4 + 460./9.*Zeta5)*nf*nf*nf +
- (1205./2916. - 152./81.*Zeta3)*nf*nf*nf*nf)/1024.0;
- Betap[0]= (double)0.25*(11.0 - 2.*nl/3.0);
- Betap[1]= (double)(102.0 - 38.*nl/3.0)/16.0;
- Betap[2]= (double)(0.5*2857. - 5033.*nl/18.0 + 325.*nl*nl/54.0)/64.0;
- Betap[3]= (double)(149753./6.0 + 3564.*Zeta3 +
- (-1078361./162.0 - 6508.*Zeta3/27.0)*nl +
- (50065./162.0 + 6472.*Zeta3/81.0)*nl*nl +
- 1093.*nl*nl*nl/729.0)/256.0;
- Betap[4]= (double)(8157455.0/16.0 + 621885.0/2.0*Zeta3 - 88209.0/2.0*Zeta4 - 288090.0*Zeta5 +
- (-336460813./1944. - 4811164./81.*Zeta3 + 33935./6.*Zeta4 + 1358995./27.*Zeta5)*nl +
- (25960913./1944. + 698531./81.*Zeta3 - 10526./9.*Zeta4 - 381760./81.*Zeta5)*nl*nl +
- (-630559./5832. - 48722./243.*Zeta3 + 1618./27.*Zeta4 + 460./9.*Zeta5)*nl*nl*nl +
- (1205./2916. - 152./81.*Zeta3)*nl*nl*nl*nl)/1024.0;
-
- Gamma[0]=(double)1.0;
- Gamma[1]=(double)(202./3.-20.*nf/9.)/16.;
- Gamma[2]=(double)(1249. + (-2216./27. - 160.*Zeta3/3.)*nf-
- 140.*nf*nf/81.)/64.;
- Gamma[3]=(double)(4603055./162. + 135680.*Zeta3/27. - 8800.*Zeta5 +
- (-91723./27. - 34192.*Zeta3/9. +
- 880.*Zeta4 + 18400.*Zeta5/9.)*nf +
- (5242./243. + 800.*Zeta3/9. - 160.*Zeta4/3.)*nf*nf +
- (-332./243. + 64.*Zeta3/27.)*nf*nf*nf)/256.;
- Gamma[4]=(double)(99512327./162. + 46402466.*Zeta3/243. + 96800.*Zeta3*Zeta3 - 698126.*Zeta4/9.
- -231757160.*Zeta5/243. + 242000.*Zeta6 + 412720.*Zeta7
- +nf*(-150736283./1458. - 12538016.*Zeta3/81. - 75680.*Zeta3*Zeta3/9. + 2038742.*Zeta4/27.
- + 49876180.*Zeta5/243. - 638000.*Zeta6/9. - 1820000.*Zeta7/27.)
- +nf*nf*(1320742./729. + 2010824.*Zeta3/243. + 46400.*Zeta3*Zeta3/27. - 166300.*Zeta4/27. - 264040.*Zeta5/81. + 92000.*Zeta6/27.)
- +nf*nf*nf*(91865./1458. + 12848.*Zeta3/81. + 448.*Zeta4/9. - 5120.*Zeta5/27.)
- +nf*nf*nf*nf*(-260./243 - 320.*Zeta3/243. + 64.*Zeta4/27.)
- )/(4*4*4*4*4);
-
-
- for(int i=0; i<5; i++) {
- B[i]=Beta[i]/Beta[0];
- Bp[i]=Betap[i]/Betap[0];
- C[i]=Gamma[i]/Beta[0];
- }
-}
-
-
-// Function int CRunDec::GetNf()
-// Returns number of active flavours.
-int CRunDec::GetNf(){
- return (int)Nf;
-}
-
-// Function void CRunDec::SetNf(int nf)
-// Set the private component Nf to the number of active flavours.
-void CRunDec::SetNf(int nf){
- this->SetConstants(nf);
-}
-
-// Aux. function to exit function.
-int CRunDec::Abbruch(void){
- RETURN
-}
-
-// PolyLog[n,x] for -1 < x < 1 and positive integer n
-double CRunDec::PolyLog(unsigned int n, double x) {
- if(x >= 1 || x <= -1 || n == 0){
- cout << "PolyLog NOT IMPLEMENTED FOR n = " << n << " and x = " << x << endl;
- RETURN
- }
- double m = n;
- double ret = 0.0;
- double sav1 = 1.0;
- double sav2 = 0.0;
- double diff = 1.0;
- double k = 1.0;
- while(abs(diff) > 1e-16) {
- sav2 = sav1;
- sav1 = pow(x,k)*pow(k,-m);
- ret += sav1;
- diff = sav2 - sav1;
- k += 1.0;
- }
- return ret;
-}
-
-// Function double CRunDec::LamExpl(double AlphaS, double Mu, int nl)
-// Compute \Lambda using eq.(4) of [RunDec].
-double CRunDec::LamExpl(double AlphaS, double Mu, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fSetAsL(x0,Mu,nl,AlphaS);
- double f1= this->fSetAsL(x1,Mu,nl,AlphaS);
- if(f0*f1>0){
- cout<<"WARNING: No root can be calculatet!"<= epsilonX);
- double Lambda=xTest;
- return Lambda;
-}
-
-// Right-hand side of differential equation for \alpha_s times 1/\mu^2,
-// see eq.(1) of [1].
-double fSetdydx(CRunDec S, double A,int nl){
- double f=0.0;
- double sum[5];
- double B=A*A;
- sum[0]=-S.Beta[0]*B;
- sum[1]=-S.Beta[1]*B*A;
- sum[2]=-S.Beta[2]*B*B;
- sum[3]=-S.Beta[3]*B*B*A;
- sum[4]=-S.Beta[4]*B*B*B;
- for(int i=1; i<= nl; i++) {
- f+=sum[i-1];
- }
- return (f*2);
-}
-
-// Implicit Runge-Kutte step (4th order)
-// Call with x=\mu^2, y=\alpha_s(\mu)/\pi, step size h, number of loops nl
-// New y value is returned at x+h
-double CRunDec::fRungeKuttaImpl(double &x, double y, double &htry, int nl,
- double (*f)(CRunDec, double, int)){
- // Precision
- double eps=1e-10;
- double yerr,ytemp,htemp, hnext;
- double h=htry;
- double xnew; // new variable
- double k1,k2,k3,k4,k5,k6;
- for(;;){
- k1=h*f(*this,y,nl);
- k2=h*f(*this,y+b21*k1,nl);
- k3=h*f(*this,y+b31*k1+b32*k2,nl);
- k4=h*f(*this,y+b41*k1+b42*k2+b43*k3,nl);
- k5=h*f(*this,y+b51*k1+b52*k2+b53*k3+b54*k4,nl);
- k6=h*f(*this,y+b61*k1+b62*k2+b63*k3+b64*k4+b65*k5,nl);
- // y value at x+h as a sum of the previous value and the
- // correspondingly weighted function evaluations
- ytemp= y+ c1*k1+ c2*k2+ c3*k3+ c4*k4+ c5*k5+ c6*k6;
- // Estimate of uncertainty
- yerr=dc1*k1 + dc2*k2 + dc3*k3 + dc4*k4 + dc5*k5 + dc6*k6;
- double err=0.;
- err=fmax(err,fabs(yerr/eps));
-
- // Uncertainty too big? -> Discard result and reduce step size
- if(err>1.){
- htemp=0.9*h*pow(err,-0.25);
- if(h>=0.){h=fmax(htemp,0.1*h);}
- else{h=fmin(htemp,0.1*h);}
- xnew=x+h; // modification to previous code
- //decide whether reduced stepsize is still big enough
- //(in order to prevent a closed loop)
- if(xnew==x){cout<<"stepsize too small"<1.89e-4){
- hnext=0.9*h*pow(err,-0.2);
- }
- // Uncertainty OK? -> take y value, increase h
- else{
- hnext=5.*h;
- }
- x+=h;
-
- y=ytemp;
- htry=hnext;
- break;
- }
- }
- return y;
-}
-
-// Function: double CRunDec::AlphasExact(double AlphaS0, double Mu0,
-// double MuEnd, int nl)
-// Compute \alpha_s using eq.(1) of [RunDec]
-double CRunDec::AlphasExact(double AlphaS0, double Mu0, double MuEnd, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fRungeKuttaImpl(x,y,h,nl,fSetdydx);
- if(x+h>=log(MuEnd)){
- h=log(MuEnd)-x;
- }
- }
- return(y*Pi);
- }
- else{h=-1e-4;}
- while(x>(log(MuEnd))){
- y=this->fRungeKuttaImpl(x,y,h,nl,fSetdydx);
- if(x+h<=log(MuEnd)){
- h=log(MuEnd)-x;
- }
- }
- return(y*Pi);
-}
-
-// Eq.(10) of [RunDec]
-double CRunDec::fSetcx(double x, int nl){
- double sum[5];
- sum[0]=1;
- sum[1]=(C[1]-B[1]*C[0])*x;
- sum[2]=0.5*((C[1]-B[1]*C[0])*(C[1]-B[1]*C[0]) + C[2] - B[1]*C[1]+
- B[1]*B[1]*C[0] - B[2]*C[0])*x*x;
- sum[3]=((C[1]-B[1]*C[0])*(C[1]-B[1]*C[0])*(C[1]-B[1]*C[0])/6. +
- 0.5*(C[1]-B[1]*C[0])*(C[2]-B[1]*C[1]+B[1]*B[1]*C[0]-B[2]*C[0])+
- (C[3]-B[1]*C[2]+B[1]*B[1]*C[1]-B[2]*C[1]-B[1]*B[1]*B[1]*C[0] +
- 2.*B[1]*B[2]*C[0] - B[3]*C[0])/3.)*x*x*x;
- sum[4]=((C[1]-B[1]*C[0])*(C[1]-B[1]*C[0])*(C[1]-B[1]*C[0])*(C[1]-B[1]*C[0])/24.
- +(C[1]-B[1]*C[0])*(C[1]-B[1]*C[0])*(C[2]/2. - B[1]*C[1]/2.+ B[1]*B[1]*C[0]/2. - B[2]*C[0]/2.)/2.
- +(C[2]/2. - B[1]*C[1]/2.+ B[1]*B[1]*C[0]/2. - B[2]*C[0]/2.)*(C[2]/2. - B[1]*C[1]/2.+ B[1]*B[1]*C[0]/2. - B[2]*C[0]/2.)/2.
- +((C[1]-B[1]*C[0])*(C[3]-B[1]*C[2]+B[1]*B[1]*C[1]-B[2]*C[1]-B[1]*B[1]*B[1]*C[0] + 2.*B[1]*B[2]*C[0] - B[3]*C[0])/3.)
- +B[1]*B[1]*B[1]*B[1]*C[0]/4. - 3.*B[1]*B[1]*B[2]*C[0]/4. + B[2]*B[2]*C[0]/4. + B[1]*B[3]*C[0]/2. - B[4]*C[0]/4.
- -B[1]*B[1]*B[1]*C[1]/4. + B[1]*B[2]*C[1]/2. - B[3]*C[1]/4. + B[1]*B[1]*C[2]/4. - B[2]*C[2]/4. - B[1]*C[3]/4. + C[4]/4.)*x*x*x*x;
- double erg=0.0;
- for(int i=1; i<=nl; i++){
- erg+=sum[i-1];
- }
- return (pow(x,C[0])*erg);
-}
-
-// Function double CRunDec::mMS2mMS(double mmu0, double AlphaS0,
-// double AlphaSEnd, int nl)
-// Compute m_q(\mu) using eqs.(9) and (10) of [RunDec]
-double CRunDec::mMS2mMS(double mmu0, double AlphaS0, double AlphaSEnd, int nl){
- if(nl<0||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fSetcx(AlphaS0/Pi, nl);
- double cAlphaSEnd= this->fSetcx(AlphaSEnd/Pi, nl);
- return mmu0*cAlphaSEnd/cAlphaS0;
-}
-
-// Aux. functions (r.h.s of diff. eqs. for alpha_s and m_q)
-double fSetdydxM1(CRunDec S, double A, double M){
- return (M*(S.Gamma[0])/(S.Beta[0]*A));
-}
-
-double fSetdydxa1(CRunDec S, double x, double A){
- if (x == 1.) x=1.;
- return (-2.*(S.Beta[0]*A*A));
-}
-
-double fSetdydxM2(CRunDec S, double A, double M){
- return (M*(S.Gamma[0]+S.Gamma[1]*A)/(S.Beta[0]*A+S.Beta[1]*A*A));
-}
-
-double fSetdydxa2(CRunDec S, double x, double A){
- if (x == 1.) x=1.;
- return (-2.*(S.Beta[0]*A*A+S.Beta[1]*A*A*A));
-}
-
-double fSetdydxM3(CRunDec S, double A, double M){
- return (M*(S.Gamma[0]+S.Gamma[1]*A+S.Gamma[2]*A*A)/
- (S.Beta[0]*A+S.Beta[1]*A*A+S.Beta[2]*A*A*A));
-}
-
-double fSetdydxa3(CRunDec S, double x, double A){
- if (x == 1.) x=1.;
- return (-2.*(S.Beta[0]*A*A+S.Beta[1]*A*A*A+S.Beta[2]*A*A*A*A));
-}
-
-double fSetdydxM4(CRunDec S, double A, double M){
- return (M*(S.Gamma[0]+S.Gamma[1]*A+S.Gamma[2]*A*A+S.Gamma[3]*A*A*A)/
- (S.Beta[0]*A+S.Beta[1]*A*A+S.Beta[2]*A*A*A+S.Beta[3]*A*A*A*A));
-}
-
-double fSetdydxa4(CRunDec S, double x, double A){
- if (x == 1.) x=1.;
- return (-2.*(S.Beta[0]*A*A+S.Beta[1]*A*A*A+S.Beta[2]*A*A*A*A+
- S.Beta[3]*A*A*A*A*A));
-}
-
-double fSetdydxM5(CRunDec S, double A, double M){
- return (M*(S.Gamma[0]+S.Gamma[1]*A+S.Gamma[2]*A*A+S.Gamma[3]*A*A*A + S.Gamma[4]*A*A*A*A)/
- (S.Beta[0]*A+S.Beta[1]*A*A+S.Beta[2]*A*A*A+S.Beta[3]*A*A*A*A+S.Beta[4]*A*A*A*A*A));
-}
-
-double fSetdydxa5(CRunDec S, double x, double A){
- if (x == 1.) x=1.;
- return (-2.*(S.Beta[0]*A*A+S.Beta[1]*A*A*A+S.Beta[2]*A*A*A*A+
- S.Beta[3]*A*A*A*A*A+S.Beta[4]*A*A*A*A*A*A));
-}
-
-// Runge-Kutta step for implicit procedure
-double CRunDec::fRKSchritt(double x, double y, double h, double &yerr,
- double (*f)(CRunDec, double , double)){
- double k1,k2,k3,k4,k5,k6;
- k1=h*f(*this,x,y);
- k2=h*f(*this,x+a2*h,y+b21*k1);
- k3=h*f(*this,x+a3*h,y+b31*k1+b32*k2);
- k4=h*f(*this,x+a4*h,y+b41*k1+b42*k2+b43*k3);
- k5=h*f(*this,x+a5*h,y+b51*k1+b52*k2+b53*k3+b54*k4);
- k6=h*f(*this,x+a6*h,y+b61*k1+b62*k2+b63*k3+b64*k4+b65*k5);
-
- yerr=(dc1*k1 + dc2*k2 + dc3*k3 + dc4*k4 + dc5*k5 + dc6*k6);
- return (y+ c1*k1+ c2*k2+ c3*k3+ c4*k4+ c5*k5+ c6*k6);
-}
-
-// Function AsmMS CRunDec::AsmMSrunexact(double mMu, double AlphaS0, double Mu0,
-// double MuEnd, int nl)
-// Compute \alpha_s and m_q solving diff. eqs. simultaneously.
-AsmMS CRunDec::AsmMSrunexact(double mMu, double AlphaS0, double Mu0,
- double MuEnd, int nl){
- AsmMS Erg;
- if(nl<0||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<Abbruch();
- }
- if(nl==0){
- Erg.Asexact=AlphaS0;
- Erg.mMSexact=mMu;
- return Erg;
- }
- double yerr0,ytemp0,xnew,h=1e-3;
- float eps=1e-15;
- float errmax;
- double x0=log(Mu0);
- double y0=AlphaS0/Pi;
- double xEnd=log(MuEnd);
- double yscal0=abs(x0)+abs(h*y0);
-
- double (*falpha)(CRunDec, double , double);
- double (*fmMS)(CRunDec, double , double);
-
- if(nl==1){
- falpha=fSetdydxa1;
- fmMS=fSetdydxM1;
- }
- if(nl==2){
- falpha=fSetdydxa2;
- fmMS=fSetdydxM2;
- }
- if(nl==3){
- falpha=fSetdydxa3;
- fmMS=fSetdydxM3;
- }
- if(nl==4){
- falpha=fSetdydxa4;
- fmMS=fSetdydxM4;
- }
- if(nl==5){
- falpha=fSetdydxa5;
- fmMS=fSetdydxM5;
- }
-
- if(Mu01){
- h*=0.9;
- xnew=x0+h;
- if(xnew==x0){cout<<"stepsize too small!"<1.89e-4){h=0.9*h*pow((double)errmax,-0.2);}
- else{h=5.*h;}
- break;
- }//else
- }//for
- if(x0+h>=xEnd){
- h=xEnd-x0;
- }//if
- }//while
- Erg.Asexact=y0*Pi;
-
- x0=AlphaS0/Pi;
- xEnd=y0;
- y0=mMu;
- yscal0=abs(x0)+abs(h*y0);
- eps=1e-10;
- h=-1e-3;
-
- while(x0>xEnd){
- for(;;){
- ytemp0=fRKSchritt(x0,y0,h,yerr0,fmMS);
- errmax=0.;
- errmax=fmax(errmax,fabs((float)yerr0/yscal0));
- errmax/=eps;
- if(errmax>1){
- h*=0.9;
- xnew=x0+h;
- if(xnew==x0){cout<<"stepsize too small!"<1.89e-4){h=0.9*h*pow((double)errmax,-0.2);}
- else{h=5.*h;}
- break;
- }//else
- }//for
- if(x0+h<=xEnd){
- h=xEnd-x0;
- }//if
- }//while
- Erg.mMSexact=y0;
- return Erg;
- }//if
- else{
- h=-1e-2;
- while(x0>xEnd){
- for(;;){
- ytemp0=fRKSchritt(x0,y0,h,yerr0,falpha);
- errmax=0.;
- errmax=fmax(errmax,fabs((float)yerr0/yscal0));
- errmax/=eps;
- if(errmax>1){
- h*=0.9;
- xnew=x0+h;
- if(xnew==x0){cout<<"stepsize too small!"<1.89e-4){h=0.9*h*pow((double)errmax,-0.2);}
- else{h=5.*h;}
- break;
- }//else
- }//for
- if(x0+h<=xEnd){
- h=xEnd-x0;
- }//if
- }//while
- Erg.Asexact=y0*Pi;
-
- x0=AlphaS0/Pi;
- xEnd=y0;
- y0=mMu;
- yscal0=abs(x0)+abs(h*y0);
- eps=1e-10;
- h=1e-3;
- while(x01){
- h*=0.9;
- xnew=x0+h;
- if(xnew==x0){cout<<"stepsize too small!"<1.89e-4){h=0.9*h*pow((double)errmax,-0.2);}
- else{h=5.*h;}
- break;
- }//else
- }//for
- if(x0+h>=xEnd){
- h=xEnd-x0;
- }//if
- }//while
- Erg.mMSexact=y0;
- return Erg;
- }//else
-}
-
-// Coefficients of eq.(13) of [RunDec]
-double CRunDec::fMsFromOs1(double mu, double M){
- double lmM=log((mu*mu)/(M*M));
- double erg;
- erg= (double) (-cf - (3.*cf*lmM)/4.);
- return erg;
-}
-
-double CRunDec::fMsFromOs2(double mu, double M, double nl){
- double lmM=log((mu*mu)/(M*M));
- double erg;
- erg= (double) ((-1111.*ca*cf)/384. + (7.*cf*cf)/128. -
- (185.*ca*cf*lmM)/96. + (21.*cf*cf*lmM)/32. - (11.*ca*cf*lmM*lmM)/32. +
- (9.*cf*cf*lmM*lmM)/32. + (143.*cf*tr)/96. + (13.*cf*lmM*tr)/24. +
- (cf*lmM*lmM*tr)/8. +
- (71.*cf*nl*tr)/96. + (13.*cf*lmM*nl*tr)/24. + (cf*lmM*lmM*nl*tr)/8. +
- (ca*cf*Zeta2)/2 - (15.*cf*cf*Zeta2)/8. - (3.*ca*cf*log(2)*Zeta2)/2. +
- 3.*cf*cf*log(2)*Zeta2 -
- cf*tr*Zeta2 + (cf*nl*tr*Zeta2)/2. + (3.*ca*cf*Zeta3)/8. -
- (3.*cf*cf*Zeta3)/4.);
- return erg;
-}
-
-double CRunDec::fMsFromOs3(double mu, double M, double nl){
- double lmM=log((mu*mu)/(M*M));
- double erg;
- erg= (double) (lmM*lmM*(-2341.*ca*ca*cf + 1962.*ca*cf*cf - 243.*cf*cf*cf
- + 1492.*ca*cf*tr -
- 468.*cf*cf*tr + 1492.*ca*cf*nl*tr - 468.*cf*cf*nl*tr - 208.*cf*tr*tr -
- 416.*cf*nl*tr*tr - 208.*cf*nl*nl*tr*tr))/1152. +
- (lmM*lmM*lmM*(-242.*ca*ca*cf + 297.*ca*cf*cf - 81.*cf*cf*cf +
- 176.*ca*cf*tr - 108.*cf*cf*tr + 176.*ca*cf*nl*tr - 108.*cf*cf*nl*tr -
- 32.*cf*tr*tr - 64.*cf*nl*tr*tr - 32.*cf*nl*nl*tr*tr))/1152. +
- (lmM*(-105944.*ca*ca*cf + 52317.*ca*cf*cf - 13203.*cf*cf*cf +
- 74624.*ca*cf*tr -
- 5436.*cf*cf*tr + 55616.*ca*cf*nl*tr + 2340.*cf*cf*nl*tr -
- 12608.*cf*tr*tr -
- 18304.*cf*nl*tr*tr - 5696.*cf*nl*nl*tr*tr + 12672.*ca*ca*cf*Zeta2 -
- 52704.*ca*cf*cf*Zeta2 + 19440.*cf*cf*cf*Zeta2 -
- 38016.*ca*ca*cf*log(2)*Zeta2 +
- 91584.*ca*cf*cf*log(2)*Zeta2 - 31104.*cf*cf*cf*log(2)*Zeta2 -
- 29952.*ca*cf*tr*Zeta2 +
- 27648.*cf*cf*tr*Zeta2 + 13824.*ca*cf*log(2)*tr*Zeta2 -
- 27648.*cf*cf*log(2)*tr*Zeta2 +
- 8064.*ca*cf*nl*tr*Zeta2 + 12096.*cf*cf*nl*tr*Zeta2 +
- 13824.*ca*cf*log(2)*nl*tr*Zeta2 -
- 27648.*cf*cf*log(2)*nl*tr*Zeta2 + 9216.*cf*tr*tr*Zeta2 +
- 4608.*cf*nl*tr*tr*Zeta2 -
- 4608.*cf*nl*nl*tr*tr*Zeta2 + 9504.*ca*ca*cf*Zeta3 -
- 22896.*ca*cf*cf*Zeta3 +
- 7776.*cf*cf*cf*Zeta3 + 6912.*ca*cf*tr*Zeta3 - 3456.*cf*cf*tr*Zeta3 +
- 6912.*ca*cf*nl*tr*Zeta3 - 3456.*cf*cf*nl*tr*Zeta3))/13824.;
- return erg;
-}
-
-//z[3,m](M) according to eq.(15)
-double CRunDec::fZmM(double nl){
- double erg;
- erg= (double) -9478333./93312. + 55.*log(2)*log(2)*log(2)*log(2)/162. +
- (-644201./6480. + 587.*log(2)/27. + 44.*log(2)*log(2)/27.)*Zeta2 -
- 61.*Zeta3/27. + 3475*Zeta4/432. + 1439.*Zeta2*Zeta3/72. -
- 1975.*Zeta5/216. + 220.*A4/27. + nl*(246643./23328. -
- log(2)*log(2)*log(2)*log(2)/81. +(967./108. + 22.*log(2)/27. -
- 4.*log(2)*log(2)/27.)*Zeta2 + 241.*Zeta3/72. - 305.*Zeta4/108. -
- 8.*A4/27.) + nl*nl*(-2353./23328. - 13.*Zeta2/54 - 7.*Zeta3/54.);
- return erg;
-}
-
-double CRunDec::fMsFromOs4(double mu, double M, double nl, double err){
- double lmM=log((mu*mu)/(M*M));
- double erg;
-
- erg = - 3654.15040757339*err - 1524.2292266911543*lmM - 288.778291935394*lmM*lmM - 32.54735725308642*lmM*lmM*lmM -
- 1.85546875*lmM*lmM*lmM*lmM +
- nl*nl*nl*(0. + 0.678141025604516*err + 0.3205521521864135*lmM + 0.0800290327210927*lmM*lmM +
- 0.010030864197530864*lmM*lmM*lmM) +
- nl*nl*(0. - 43.48241924867489*err - 19.82672048099557*lmM - 4.482957520194182*lmM*lmM -
- 0.5270061728395061*lmM*lmM*lmM - 0.04108796296296297*lmM*lmM*lmM*lmM) +
- nl*(0. + 756.9421565599532*err + 330.1770776731065*lmM + 67.99849534415492*lmM*lmM +
- 7.595293209876542*lmM*lmM*lmM + 0.48119212962962954*lmM*lmM*lmM*lmM);
- return erg;
-}
-
-// Function: double CRunDec::mOS2mMS(double mOS, std::pair* mq, double asmu,
-// double Mu, int nl, double fdelm)
-double CRunDec::mOS2mMS(double mOS, std::pair* mq, double asmu, double Mu,int nl, double fdelm){
- if(nl<0||nl>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<6)){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR 4 LOOPS and "<< nl << " FLAVORS" <fMsFromOs1(Mu, mOS))/Pi;
- sum[2]=asmu*asmu*(this-> fMsFromOs2(Mu, mOS, Nf-1))/(Pi*Pi);
- sum[3]=asmu*asmu*asmu*(this-> fMsFromOs3(Mu, mOS, Nf-1)+this->fZmM(Nf-1))/(Pi*Pi*Pi);
- sum[4]=asmu*asmu*asmu*asmu*(this->fMsFromOs4(Mu, mOS, Nf-1, fdelm))/(Pi*Pi*Pi*Pi);
- double erg=0.0;
- if(nl==0){
- erg=1;
- }
- else{
- for(int i=0; i<=nl; i++){
- erg+=sum[i];
- }
- }
- erg += deltalight;
- return mOS*erg;
-}
-
-// Function: double CRunDec::mMS2mSI(double mMS, double asmu, double mu, int nl)
-double CRunDec::mMS2mSI(double mMS, double asmu, double mu, int nl){
- double epsilonX = 1e-8;
- AsmMS asmq;
- asmq.Asexact = asmu;
- asmq.mMSexact = mMS;
- for (;;) {
- double mbold = asmq.mMSexact;
- asmq = AsmMSrunexact(mMS, asmu, mu, mbold, nl);
- if (abs(asmq.mMSexact - mbold) < epsilonX) break;
- }
- return asmq.mMSexact;
-}
-
-// Coefficients of eq.(18) of [RunDec]
-double CRunDec::fMsFromRi1(void){
- return (double) -4./3.;
-}
-double CRunDec::fMsFromRi2(void){
- double erg= (double) -995./72. + 19.*Zeta3/6. + 89.*Nf/144.;
- return erg;
-}
-double CRunDec::fMsFromRi3(void){
- double erg= (double) -6663911./41472. + 408007.*Zeta3/6912. -185.*Zeta5/36.
- + (118325./7776. + 5*Zeta4/12. - 617.*Zeta3/216.)*Nf +
- (-4459./23328. - Zeta3/54.)*Nf*Nf;
- return erg;
-}
-
-// Function: double CRunDec::mRI2mMS(double mRI, double asmu, int nl)
-double CRunDec::mRI2mMS(double mRI, double asmu, int nl){
- if(nl<0||nl>3){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fMsFromRi1())/Pi;
- sum[2]=asmu*asmu*(this-> fMsFromRi2())/(Pi*Pi);
- sum[3]=asmu*asmu*asmu*(this-> fMsFromRi3())/(Pi*Pi*Pi);
- double erg=0.0;
- if(nl==0){
- erg=1;
- }
- else{
- for(int i=0; i<=nl; i++){
- erg+=sum[i];
- }
- }
- return mRI*erg;
-}
-
-// Function: double CRunDec::mMS2mRGI(double mMS, double asmu, int nl)
-double CRunDec::mMS2mRGI(double mMS, double asmu, int nl){
- if(nl<0||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fSetcx(asmu/Pi, nl);
- return (double) mMS/cAsmu;
- }
-}
-
-// Function: double CRunDec::mRGI2mMS(double mRGI, double asmu, int nl)
-double CRunDec::mRGI2mMS(double mRGI, double asmu, int nl){
- if(nl<0||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fSetcx(asmu/Pi, nl);
- return (double) mRGI*cAsmu;
-}
-
-// Coefficients of eq.(17) of [RunDec]
-double CRunDec::fOsFromMs1(double mu, double M){
- double lmM=log((mu*mu)/(M*M));
- double erg;
- erg= (double) (cf + (3.*cf*lmM)/4.);
- return erg;
-}
-
-double CRunDec::fOsFromMs2(double mu, double M, double nl){
- double lmM=log((mu*mu)/(M*M));
- double erg;
- erg= (double) ((1111.*ca*cf)/384. - (71.*cf*cf)/128. -
- (143.*cf*tr)/96. - (71.*cf*nl*tr)/96. +
- lmM*((185.*ca*cf)/96. - (9.*cf*cf)/32. - (13.*cf*tr)/24. -
- (13.*cf*nl*tr)/24.) +
- lmM*lmM*((11.*ca*cf)/32. + (9.*cf*cf)/32. - (cf*tr)/8. - (cf*nl*tr)/8.) -
- (ca*cf*Zeta2)/2. + (15.*cf*cf*Zeta2)/8. + (3.*ca*cf*log(2)*Zeta2)/2.
- - 3.*cf*cf*log(2)*Zeta2 +
- cf*tr*Zeta2 - (cf*nl*tr*Zeta2)/2. - (3.*ca*cf*Zeta3)/8. +
- (3.*cf*cf*Zeta3)/4.);
- return erg;
-}
-
-double CRunDec::fOsFromMs3(double mu, double M, double nl){
- double lmM=log((mu*mu)/(M*M));
- double erg;
-
- erg= (double) (lmM*lmM*lmM*((121.*ca*ca*cf)/576. + (33.*ca*cf*cf)/128. +
- (9.*cf*cf*cf)/128. - (11.*ca*cf*tr)/72. - (3.*cf*cf*tr)/32. -
- (11.*ca*cf*nl*tr)/72. -
- (3.*cf*cf*nl*tr)/32. + (cf*tr*tr)/36. + (cf*nl*tr*tr)/18. +
- (cf*nl*nl*tr*tr)/36.) + lmM*lmM*((2341.*ca*ca*cf)/1152.
- + (21.*ca*cf*cf)/64. -
- (63.*cf*cf*cf)/128. - (373.*ca*cf*tr)/288. - (3.*cf*cf*tr)/32. -
- (373.*ca*cf*nl*tr)/288. - (3.*cf*cf*nl*tr)/32. + (13.*cf*tr*tr)/72. +
- (13.*cf*nl*tr*tr)/36. + (13.*cf*nl*nl*tr*tr)/72.) +
- lmM*((13243.*ca*ca*cf)/1728. - (4219.*ca*cf*cf)/1536. +
- (495.*cf*cf*cf)/512. -
- (583.*ca*cf*tr)/108. - (307.*cf*cf*tr)/384. - (869.*ca*cf*nl*tr)/216. -
- (91.*cf*cf*nl*tr)/384. + (197.*cf*tr*tr)/216. + (143.*cf*nl*tr*tr)/108. +
- (89.*cf*nl*nl*tr*tr)/216. - (11.*ca*ca*cf*Zeta2)/12. +
- (49.*ca*cf*cf*Zeta2)/16. +
- (45.*cf*cf*cf*Zeta2)/32. + (11.*ca*ca*cf*log(2)*Zeta2)/4. -
- (35.*ca*cf*cf*log(2)*Zeta2)/8. -
- (9.*cf*cf*cf*log(2)*Zeta2)/4. + (13.*ca*cf*tr*Zeta2)/6. -
- (cf*cf*tr*Zeta2)/2. -
- ca*cf*log(2)*tr*Zeta2 + 2.*cf*cf*log(2)*tr*Zeta2 -
- (7.*ca*cf*nl*tr*Zeta2)/12. -
- (13.*cf*cf*nl*tr*Zeta2)/8. - ca*cf*log(2)*nl*tr*Zeta2 +
- 2.*cf*cf*log(2)*nl*tr*Zeta2 -
- (2.*cf*tr*tr*Zeta2)/3. - (cf*nl*tr*tr*Zeta2)/3. +
- (cf*nl*nl*tr*tr*Zeta2)/3. -
- (11.*ca*ca*cf*Zeta3)/16. + (35.*ca*cf*cf*Zeta3)/32. +
- (9.*cf*cf*cf*Zeta3)/16. -
- (ca*cf*tr*Zeta3)/2. + (cf*cf*tr*Zeta3)/4. - (ca*cf*nl*tr*Zeta3)/2. +
- (cf*cf*nl*tr*Zeta3)/4.));
- return erg;
-}
-
-double CRunDec::fOsFromMs4(double mu, double M, double nl, double err){
- double lmM=log((mu*mu)/(M*M));
- double erg;
- erg = 3567.602784989066*err + 1727.2260148986106*lmM + 409.2429990574718*lmM*lmM + 66.93663194444443*lmM*lmM*lmM +
- 8.056278935185185*lmM*lmM*lmM*lmM +
- nl*nl*nl*(-0.678141025604516*err - 0.3205521521864134*lmM - 0.0800290327210927*lmM*lmM -
- 0.010030864197530864*lmM*lmM*lmM) +
- nl*(-745.7207145811878*err - 358.29765085086774*lmM - 87.39262571554698*lmM*lmM -
- 11.883873456790122*lmM*lmM*lmM - 1.2705439814814814*lmM*lmM*lmM*lmM) +
- nl*nl*(43.396250117985666*err + 20.528466368867228*lmM + 4.971905254812516*lmM*lmM +
- 0.6304012345679011*lmM*lmM*lmM + 0.06655092592592593*lmM*lmM*lmM*lmM);
- return erg;
-}
-
-// Compute 2 Loop contribution to mOS2mMS at 3 loops for MS light quarks
-double CRunDec::deltamOS2mMS(double mOS, std::pair* mq, double asmu, double mu, int nlq, int nloops){
- double erg=0.0;
- if(!mq)
- return 0.0;
-
- for(int i = 0; i < 4; i++) {
- if(mq[i].first == 0.0)
- continue;
-
- double x = mq[i].first/mOS;
- double muf = mq[i].second;
- int nl = nlq-i;
- if(nloops >= 2) {
- erg += asmu*asmu*2./3.*(48.*x*x*x*x*log(x)*log(x) + 48.*x*x*log(x)+72.*x*x + 8.*Pi*Pi*(x*x*x*x-3.*x*x*x-3.*x)
- -48.*(x+1.)*(x+1.)*(x*x - x + 1.)*(log(x)*log(x+1) + PolyLog(2,-x))
- -48.*(x-1.)*(x-1.)*(x*x + x + 1.)*(log(x)*log(1-x) + PolyLog(2,x)))/96.;
- }
- if(nloops >= 3) {
- erg += asmu*asmu*asmu*(-21.8714*x - 4.348*x*x - 1.02211*x*x*x - 0.0493333*x*x*x*x + nl*(0.982667*x + 0.300333*x*x)
- + log(mu*mu/(mOS*mOS))*(-6.61056*x + 2.46511*x*x - 0.724333*x*x*x + nl*(0.534667*x - 0.22*x*x + 0.067*x*x*x))
- + log(x)*(16.9477*x - 1.10133*nl*x + 2.78756*x*x - 0.0343333*x*x*x)
- + 8./9.*(3./2.*log(muf*muf/(mq[i].first*mq[i].first))+2.)*x/24.*(24.*x*x*x*log(x)*log(x) + 12.*x*log(x)+ 24.*x
- + Pi*Pi*(4.*x*x*x-9.*x*x-3.) - 6.*(4.*x*x*x + 3.*x*x + 1.)*(log(x)*log(x+1) + PolyLog(2,-x))
- -6.*(x-1.)*(4.*x*x + x + 1.)*(log(x)*log(1-x) + PolyLog(2,x))));
- }
- }
-
- return erg;
-}
-
-// Compute 2 Loop contribution to mMS2mOS at 3 loops for MS light quarks
-double CRunDec::deltamMS2mOS(double mMS, std::pair* mq, double asmu, double mu, int nlq, int nloops){
- double erg=0.0;
- double lmu = log(mu*mu/(mMS*mMS));
- if(!mq)
- return 0.0;
-
- for(int i = 0; i < 4; i++) {
- if(mq[i].first == 0.0)
- continue;
-
- double x = mq[i].first/mMS;
- double muf = mq[i].second;
- int nl = nlq-i;
- if(nloops >= 2) {
- erg += -asmu*asmu*2./3.*(48.*x*x*x*x*log(x)*log(x) + 48.*x*x*log(x)+72.*x*x + 8.*Pi*Pi*(x*x*x*x-3.*x*x*x-3.*x)
- -48.*(x+1.)*(x+1.)*(x*x - x + 1.)*(log(x)*log(x+1) + PolyLog(2,-x))
- -48.*(x-1.)*(x-1.)*(x*x + x + 1.)*(log(x)*log(1-x) + PolyLog(2,x)))/96.;
- }
- if(nloops >= 3) {
- erg += -asmu*asmu*asmu*(-21.8714*x - 4.348*x*x - 1.02211*x*x*x - 0.0493333*x*x*x*x + nl*(0.982667*x + 0.300333*x*x)
- + lmu*(-6.61056*x + 2.46511*x*x - 0.724333*x*x*x + nl*(0.534667*x - 0.22*x*x + 0.067*x*x*x))
- + log(x)*(16.9477*x - 1.10133*nl*x + 2.78756*x*x - 0.0343333*x*x*x)
- + 8./9.*(3./2.*log(muf*muf/(mq[i].first*mq[i].first))+2.)*x/24.*(24.*x*x*x*log(x)*log(x) + 12.*x*log(x)+ 24.*x
- + Pi*Pi*(4.*x*x*x-9.*x*x-3.) - 6.*(4.*x*x*x + 3.*x*x + 1.)*(log(x)*log(x+1) + PolyLog(2,-x))
- -6.*(x-1.)*(4.*x*x + x + 1.)*(log(x)*log(1-x) + PolyLog(2,x))))
- + asmu*asmu*asmu*8./9.*(48.*(3.*lmu + 7.)*x*x*x*x*log(x)*log(x) + 144.*x*x*log(x)+312.*x*x
- + 8.*Pi*Pi*(7.*x*x*x*x-15.*x*x*x-3.*x) - lmu*(-72.*x*x +
- + 12.*Pi*Pi*(-2.*x*x*x*x+3.*x*x*x - 3*x))
- - 48.*((7.*x+5.)*x*x*x + x + 3.*lmu/2.*(2.*x*x*x*x+x*x*x-x-2.) -1.)*(log(x)*log(x+1) + PolyLog(2,-x))
- - 48.*(x-1.)*(3.*lmu/2.*(2.*x*x*x+x*x+x+2.) + 7.*x*x*x + 2.*x*x + 2.*x + 1.)*(log(x)*log(1-x) + PolyLog(2,x)))/192.;
- }
- }
-
- return erg;
-}
-
-// z_m^inv
-double CRunDec::fZmInvM(double nl){
- double erg;
- erg=(8481925./93312. +
- (137.*nl)/216. + (652841.*Pi*Pi)/38880. - (nl*Pi*Pi)/27. -
- (695.*Pi*Pi*Pi*Pi)/7776. - (575.*Pi*Pi*log(2))/162. -
- (22.*Pi*Pi*log(2)*log(2))/81. -
- (55.*log(2)*log(2)*log(2)*log(2))/162. - (220.*A4)/27. -
- nl*nl*(-2353./23328. - (13.*Pi*Pi)/324. - (7.*Zeta3)/54.) +
- (58.*Zeta3)/27. -
- (1439.*Pi*Pi*Zeta3)/432. - nl*(246643./23328. + (967.*Pi*Pi)/648. -
- (61.*Pi*Pi*Pi*Pi)/1944. + (11.*Pi*Pi*log(2))/81. -
- (2.*Pi*Pi*log(2)*log(2))/81. -
- log(2)*log(2)*log(2)*log(2)/81. - (8.*A4)/27. +
- (241.*Zeta3)/72.) +
- (1975.*Zeta5)/216.);
- return erg;
-}
-
-// Function: double CRunDec::mMS2mOS(double mMS, std::pair* mq, double asmu,
-// double mu,int nl)
-double CRunDec::mMS2mOS(double mMS, std::pair* mq, double asmu, double mu, int nl, double fdelm){
- if(nl<0||nl>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<6)){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR 4 LOOPS and "<< nl << " FLAVORS"<fOsFromMs1(mu, mMS))/Pi;
- sum[2]=asmu*asmu*((this-> fOsFromMs2(mu, mMS, Nf-1)))/(Pi*Pi);
- sum[3]=asmu*asmu*asmu*(this-> fOsFromMs3(mu, mMS, Nf-1)+
- this->fZmInvM(Nf-1))/(Pi*Pi*Pi);
- sum[4]=asmu*asmu*asmu*asmu*(this->fOsFromMs4(mu, mMS, Nf-1, fdelm))/(Pi*Pi*Pi*Pi);
- double erg=0.0;
- if(nl==0){
- erg=1;
- }
- else{
- for(int i=0; i<=nl; i++){
- erg+=sum[i];
- }
- }
- erg += deltalight;
- return mMS*erg;
-}
-
-// Function: double CRunDec::mMS2mOSmod(double mMS, std::pair* mq, double asmu,
-// double mu, int nf, int nloop, double fdelm), needed for several other mass relations
-double CRunDec::mMS2mOSmod(double mMS, std::pair* mq, double asmu, double mu,int nf, int nloop, double fdelm){
- if(nloop<0||nloop>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloop <<" LOOPS"<6)){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR 4 LOOPS and "<< nf << " FLAVORS"<fOsFromMs1(mu, mMS))/(Pi);
- sum[2]=asmu*asmu*((this-> fOsFromMs2(mu, mMS, Nf-1)))/(Pi*Pi) +
- asmu*asmu*(lmm)/(6.*Pi)*(this ->fOsFromMs1(mu, mMS))/(Pi);
- sum[3]=asmu*asmu*asmu*(this-> fOsFromMs3(mu, mMS, Nf-1)+this->fZmInvM(Nf-1))/(Pi*Pi*Pi)+
- asmu*asmu*asmu*2.*(lmm)/(6.*Pi)*((this-> fOsFromMs2(mu, mMS, Nf-1)))/(Pi*Pi)+
- asmu*asmu*asmu*(-11./72. + (11.*lmm)/24. + lmm*lmm/36.)/(Pi*Pi)*(this ->fOsFromMs1(mu, mMS))/(Pi);
- sum[4]=asmu*asmu*asmu*asmu*(this->fOsFromMs4(mu, mMS, Nf-1, fdelm))/(Pi*Pi*Pi*Pi) +
- asmu*asmu*asmu*asmu*3.*(lmm)/(6.*Pi)*(this-> fOsFromMs3(mu, mMS, Nf-1)+this->fZmInvM(Nf-1))/(Pi*Pi*Pi) +
- asmu*asmu*asmu*asmu*2.*(-11./72. + (11.*lmm)/24. + lmm*lmm/36.)/(Pi*Pi)*
- ((this-> fOsFromMs2(mu, mMS, Nf-1)))/(Pi*Pi) +
- asmu*asmu*asmu*asmu*(-564731./124416. + (2645.*lmm)/1728. + (167.*lmm*lmm)/576. + lmm*lmm*lmm/216. +
- (2633./31104. - (67.*lmm)/576. + lmm*lmm/36.)*(Nf-1) + (82043.*Zeta3)/27648.)/(Pi*Pi*Pi)
- *(this ->fOsFromMs1(mu, mMS))/(Pi) +
- asmu*asmu*asmu*asmu*(lmm)/(6.*Pi)*(lmm)/(6.*Pi)*((this-> fOsFromMs2(mu, mMS, Nf-1)))/(Pi*Pi);
-
- double erg=0.0;
- if(nloop==0){
- erg=1;
- }
- else{
- for(int i=0; i<=nloop; i++){
- erg+=sum[i];
- }
- }
- return mMS*erg;
-}
-
-// Coefficients of eq.(16) of [RunDec]
-double CRunDec::fMumFromOs1(void){
- return (double) -cf;
-}
-
-double CRunDec::fMumFromOs2(void){
- double erg;
- erg= (double) ((-1111.*ca*cf)/384. + (199.*cf*cf)/128. + (143.*cf*tr)/96. +
- (71.*cf*(Nf-1)*tr)/96. + (ca*cf*Zeta2)/2. - (15.*cf*cf*Zeta2)/8. -
- (3.*ca*cf*log(2)*Zeta2)/2. +
- 3.*cf*cf*log(2)*Zeta2 - cf*tr*Zeta2 + (cf*(Nf-1)*tr*Zeta2)/2. +
- (3.*ca*cf*Zeta3)/8. -
- (3.*cf*cf*Zeta3)/4.);
- return erg;
-}
-
-// z_m^SI
-double CRunDec::fMumFromOs3(void){
- double erg;
- erg= (double) -7172965./93312. -
- (293.*(Nf-1))/216. - (618281.*Pi*Pi)/38880. - ((Nf-1)*Pi*Pi)/9. +
- (695.*Pi*Pi*Pi*Pi)/7776. + (623.*Pi*Pi*log(2))/162. +
- (22.*Pi*Pi*log(2)*log(2))/81. +
- (55.*log(2)*log(2)*log(2)*log(2))/162. + (220.*A4)/27. +
- (Nf-1)*(Nf-1)*(-2353./23328. - (13.*Pi*Pi)/324. - (7.*Zeta3)/54.) -
- (70.*Zeta3)/27. +
- (1439.*Pi*Pi*Zeta3)/432. + (Nf-1)*(246643./23328. + (967.*Pi*Pi)/648. -
- (61.*Pi*Pi*Pi*Pi)/1944. + (11.*Pi*Pi*log(2))/81. -
- (2*Pi*Pi*log(2)*log(2))/81. -
- log(2)*log(2)*log(2)*log(2)/81. - (8.*A4)/27. + (241.*Zeta3)/72.) -
- (1975.*Zeta5)/216.;
- return erg;
-}
-
-double CRunDec::fMumFromOs4(double err){
- double erg;
- erg = err*(-3214.227044839041 + 692.4809215366435*(-1. + Nf) - 41.95978562498058*(-1. + Nf)*(-1. + Nf)
- + 0.678141025604516*(-1. + Nf)*(-1. + Nf)*(-1. + Nf));
- return erg;
-}
-
-// Function: double CRunDec::mOS2mSI(double mOS, std::pair* mq, double asM,
-// int nl)
-double CRunDec::mOS2mSI(double mOS, std::pair* mq, double asM, int nl, double fdelm){
- if(nl<0||nl>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<6)){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR 4 LOOPS and "<< Nf << " FLAVORS"<fMumFromOs1())/Pi;
- sum[2]=asM*asM*(this-> fMumFromOs2())/(Pi*Pi);
- sum[3]=asM*asM*asM*(this-> fMumFromOs3())/(Pi*Pi*Pi);
- sum[4]=asM*asM*asM*asM*(this->fMumFromOs4(fdelm))/(Pi*Pi*Pi*Pi);
-
-
- double erg=0.0;
- if(nl==0){
- erg=1;
- }
- else{
- for(int i=0; i<=nl; i++){
- erg+=sum[i];
- }
- }
- erg += deltalight;
- return mOS*erg;
-
-}
-
-// Function: double CRunDec::mOS2mMSrun(double mOS, std::pair* mq, double asmu,
-// double mu, int nl)
-double CRunDec::mOS2mMSrun(double mOS, std::pair* mq, double asmu, double mu,
- int nl){
- double asM=0.0;
- asM= this-> AlphasExact(asmu, mu, mOS, nl);
- double mum= this-> mOS2mSI(mOS, mq, asM, nl);
- double asmum= this-> AlphasExact(asmu, mu, mum, nl);
- double newM= this->mMS2mMS(mum, asmum, asmu, nl);
- return newM;
-}
-
-// Function: double CRunDec::mMS2mOSrun(double mMS, std::pair* mq, double asmu,
-// double mu, int nl)
-double CRunDec::mMS2mOSrun(double mMS, std::pair* mq, double asmu, double mu,
- int nl){
- double mNeu = mMS2mSI(mMS, asmu, mu, nl);
- double asmNeu = AlphasExact(asmu, mu, mNeu, nl);
- return mMS2mOS(mNeu, mq, asmNeu, mNeu, nl);
-}
-
-// Coefficients of eq.(19) of [RunDec]
-double CRunDec::fRiFromMs(double alpha, double nl){
- double sum[4];
- sum[0]= 1.;
- sum[1]= (4.*alpha)/3.;
- sum[2]= alpha*alpha*((1123./72. - (89.*Nf)/144. - (19.*Zeta3)/6.));
- sum[3]= alpha*alpha*alpha*(6663911./41472. - (118325.*Nf)/7776. +
- (4459.*Nf*Nf)/23328. +
- (4.*(1123./72. - (89.*Nf)/144. - (19.*Zeta3)/6.))/3. -
- (408007.*Zeta3)/6912. +
- (617.*Nf*Zeta3)/216. + (Nf*Nf*Zeta3)/54. - (4.*(-995./72. +
- (89.*Nf)/144. + (19.*Zeta3)/6.))/3. -
- (5.*Nf*Zeta4)/12. + (185.*Zeta5)/36.);
- double erg=0.0;
- if(nl==0){
- erg=1.;
- }
- else{
- erg=0.;
- for(int i=0; i<=nl; i++){
- erg+=sum[i];
- }
- }
- return (erg);
-}
-
-// Function: double CRunDec::mMS2mRI(double mMS, double asmu, int nl)
-double CRunDec::mMS2mRI(double mMS, double asmu, int nl){
- if(nl<0||nl>3){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fRiFromMs((asmu/Pi), nl));
-}
-
-// Coefficients needed for the transformation of mOS to mMSit
-double CRunDec::fHelpmOS2mMSit(double mMS, double mOS, std::pair* mq, double asmu,
- double mu, int nl){
- double sum[4];
- double deltalight = deltamMS2mOS(mMS,mq,asmu/Pi,mu,Nf-1,nl);
- sum[0]=(double) 1.;
- sum[1]=asmu*(this ->fOsFromMs1(mu, mMS))/Pi;
- sum[2]=asmu*asmu*((this-> fOsFromMs2(mu, mMS, Nf-1)))/(Pi*Pi);
- sum[3]=asmu*asmu*asmu*(this-> fOsFromMs3(mu, mMS,Nf-1)+
- +this->fZmInvM(Nf-1))/(Pi*Pi*Pi);
- double erg=0.0;
- if(nl==0){
- erg=1;
- }
- else{
- for(int i=0; i<=nl; i++){
- erg+=sum[i];
- }
- }
- erg += deltalight;
- return erg;
-}
-
-// Function: double CRunDec::mOS2mMSit(double mOS, std::pair* mq, double asmu,
-// double mu, int nl)
-double CRunDec::mOS2mMSit(double mOS, std::pair* mq, double asmu, double mu,
- int nl){
- if(nl<0||nl>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fHelpmOS2mMSit(x0,mOS, mq, asmu, mu, nl))-mOS);
- double f1= (x1*(this->fHelpmOS2mMSit(x1,mOS, mq, asmu, mu, nl))-mOS);
- if(f0*f1>0){
- cout<<"WARNING: No root can be calculatet!"<fHelpmOS2mMSit(xTest,mOS, mq, asmu, mu, nl))-mOS);
- if(f0*fTest<=0){x1= xTest;}
- else {x0= xTest;}
- }
- while(abs(x1-x0)>= epsilonX);
- double mNeu=xTest;
- return mNeu;
-}
-
-double CRunDec::PSdelta(double asmu, double muf, double mu, int nl, int nloops) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu, double mu, double muf, int nl, int nloops)
-double CRunDec::mOS2mPS(double mOS, std::pair* mq, double asmu, double mu, double muf, int nl, int nloops) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
-// double mu, double muf, int nl, int nloops, double fdelm)
-double CRunDec::mMS2mPS(double mMS, std::pair* mq, double asmu,
- double mu, double muf, int nl, int nloops, double fdelm) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<mq, asmu, mu, nl+1, nloops, fdelm);
- return (exmOS - muf*delmuf);
-}
-
-// Function: double CRunDec::mPS2mMS(double mPS, std::pair* mq, double asmu,
-// double mu, double muf, int nl, int nloops, double fdelm)
-double CRunDec::mPS2mMS(double mPS, std::pair* mq, double asmu,
- double mu, double muf, int nl, int nloops, double fdelm) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
-// double (*as)(double), double muf, int nl, int nloops, double fdelm)
-// The function pointer passed should contain the adress of a function computing alpha_s in dependence of mu
-double CRunDec::mPS2mSI(double mPS, std::pair* mq,
- double (*as)(double), double muf, int nl, int nloops, double fdelm) {
- if(as == NULL) {
- cout << "Pointer to as == NULL! Aborting..." << endl;
- RETURN
- }
- double mMS1 = 0;
- double mMS = mPS;
- double acc = 10e-6;
- while(fabs(mMS1 - mMS) > acc) {
- mMS1 = mMS;
- mMS = mPS2mMS(mPS, mq, as(mMS1), mMS1, muf, nl, nloops, fdelm);
- }
- return mMS;
-}
-
-
-double CRunDec::E1p(double mOS, double asmu, double mu, int nl, int nloops) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
-// double asmu, double mu, int nl, int nloops)
-double CRunDec::mOS2m1S(double mOS, std::pair* mq, double asmu, double mu, int nl, int nloops) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
-// double asmu, double mu, int nl, int nloops, double fdelm)
-double CRunDec::mMS2m1S(double mMS, std::pair* mq,
- double asmu, double mu, int nl, int nloops, double fdelm) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR 4 LOOPS and " << nl << " FLAVORS"<fOsFromMs4(mu,mMS,nl,fdelm)) +
- (asmu*asmu*mMS*(2763 - 142*nl + 96*Pi*Pi - 16*nl*Pi*Pi + 32*Pi*Pi*ln2 +
- 2036*lmm - 104*nl*lmm + 564*lmm*lmm - 24*nl*lmm*lmm - 48*Zeta3))/
- (288*Pi*Pi) + (-2763*asmu*asmu*asmu*asmu*mMS + 142*asmu*asmu*asmu*asmu*mMS*nl -
- 96*asmu*asmu*asmu*asmu*mMS*Pi*Pi + 16*asmu*asmu*asmu*asmu*mMS*nl*Pi*Pi -
- 32*asmu*asmu*asmu*asmu*mMS*Pi*Pi*ln2 - 2036*asmu*asmu*asmu*asmu*mMS*lmm +
- 104*asmu*asmu*asmu*asmu*mMS*nl*lmm - 564*asmu*asmu*asmu*asmu*mMS*lmm*lmm +
- 24*asmu*asmu*asmu*asmu*mMS*nl*lmm*lmm + 48*asmu*asmu*asmu*asmu*mMS*Zeta3)/(1296*Pi*Pi)
- + (-129096*asmu*asmu*asmu*asmu*mMS + 20316*asmu*asmu*asmu*asmu*mMS*nl - 616*asmu*asmu*asmu*asmu*mMS*nl*nl -
- 11668*asmu*asmu*asmu*asmu*mMS*Pi*Pi + 528*asmu*asmu*asmu*asmu*mMS*nl*Pi*Pi -
- 16*asmu*asmu*asmu*asmu*mMS*nl*nl*Pi*Pi + 243*asmu*asmu*asmu*asmu*mMS*Pi*Pi*Pi*Pi -
- 200232*asmu*asmu*asmu*asmu*mMS*log34 +
- 27792*asmu*asmu*asmu*asmu*mMS*nl*log34 -
- 864*asmu*asmu*asmu*asmu*mMS*nl*nl*log34 -
- 78408*asmu*asmu*asmu*asmu*mMS*log34*log34 +
- 9504*asmu*asmu*asmu*asmu*mMS*nl*log34*log34 -
- 288*asmu*asmu*asmu*asmu*mMS*nl*nl*log34*log34 -
- 59400*asmu*asmu*asmu*asmu*mMS*Zeta3 + 8208*asmu*asmu*asmu*asmu*mMS*nl*Zeta3 -
- 192*asmu*asmu*asmu*asmu*mMS*nl*nl*Zeta3)/(3888*Pi*Pi) +
- (asmu*asmu*asmu*mMS*(42314585 - 4636940*nl + 47060*nl*nl + 7834092*Pi*Pi -
- 713520*nl*Pi*Pi + 18720*nl*nl*Pi*Pi - 41700*Pi*Pi*Pi*Pi + 14640*nl*Pi*Pi*Pi*Pi -
- 1656000*Pi*Pi*ln2 - 63360*nl*Pi*Pi*ln2 - 126720*Pi*Pi*ln2*ln2 +
- 11520*nl*Pi*Pi*ln2*ln2 - 158400*ln2*ln2*ln2*ln2 + 5760*nl*ln2*ln2*ln2*ln2 +
- 33620760*lmm - 3723120*nl*lmm +
- 64080*nl*nl*lmm + 1010880*Pi*Pi*lmm -
- 220320*nl*Pi*Pi*lmm + 8640*nl*nl*Pi*Pi*lmm +
- 336960*Pi*Pi*ln2*lmm -
- 17280*nl*Pi*Pi*ln2*lmm + 11726100*lmm*lmm -
- 1247400*nl*lmm*lmm + 28080*nl*nl*lmm*lmm +
- 2009880*lmm*lmm*lmm - 185760*nl*lmm*lmm*lmm +
- 4320*nl*nl*lmm*lmm*lmm - 3801600*A4 +
- 138240*nl*A4 + 1002240*Zeta3 - 1561680*nl*Zeta3 +
- 60480*nl*nl*Zeta3 - 1554120*Pi*Pi*Zeta3 -
- 894240*lmm*Zeta3 - 362880*nl*lmm*Zeta3 +
- 4266000*Zeta5))/(466560*Pi*Pi*Pi);
- break;
- }
-
-
-
- return m1S;
-}
-
-// Function: double CRunDec::m1S2mMS(double m1S, std::pair* mq,
-// double asmu, double mu, int nl, int nloops, double fdelm)
-double CRunDec::m1S2mMS(double m1S, std::pair* mq,
- double asmu, double mu, int nl, int nloops, double fdelm) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
-// double (*as)(double), int nl, int nloops, double fdelm)
-// The function pointer passed should contain the adress of a function computing alpha_s in dependence of mu
-double CRunDec::m1S2mSI(double m1S, std::pair* mq,
- double (*as)(double), int nl, int nloops, double fdelm) {
- if(as == NULL) {
- cout << "Pointer to as == NULL! Aborting..." << endl;
- RETURN
- }
- double mMS1 = 0;
- double mMS = m1S;
- double acc = 10e-6;
- while(fabs(mMS1 - mMS) > acc) {
- mMS1 = mMS;
- mMS = m1S2mMS(m1S, mq, as(mMS1), mMS1, nl, nloops, fdelm);
- }
- return mMS;
-}
-
-double CRunDec::exOS2RS(double api, double mmu, double nnuf, int nnl, int nloops) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
-// double mu, double nuf, int nl, int nloops, bool prime)
-double CRunDec::mOS2mRS(double mOS, std::pair* mq, double asmu,
- double mu, double nuf, int nl, int nloops, bool prime) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
-// double mu, double nuf, int nl, int nloops, double fdelm, bool prime)
-double CRunDec::mMS2mRS(double mMS, std::pair* mq, double asmu,
- double mu, double nuf, int nl, int nloops, double fdelm, bool prime) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
-// double mu, double muf, int nl, int nloops, double fdelm, bool prime)
-double CRunDec::mRS2mMS(double mRS, std::pair* mq, double asmu,
- double mu, double muf, int nl, int nloops, double fdelm, bool prime) {
- if(nloops<0||nloops>4){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double (*as)(double),
-// double muf, int nl, int nloops, double fdelm, bool prime)
-// The function pointer passed should contain the adress of a function computing alpha_s in dependence of mu
-double CRunDec::mRS2mSI(double mRS, std::pair* mq, double (*as)(double),
- double muf, int nl, int nloops, double fdelm, bool prime) {
- if(as == NULL) {
- cout << "Pointer to as == NULL! Aborting..." << endl;
- RETURN
- }
- double mMS1 = 0;
- double mMS = mRS;
- double acc = 10e-6;
- while(fabs(mMS1 - mMS) > acc) {
- mMS1 = mMS;
- mMS = mRS2mMS(mRS, mq, as(mMS1), mMS1, muf, nl, nloops, fdelm, prime);
- }
- return mMS;
-}
-
-// Function: double CRunDec::mMS2mRGImod(double mMS, double asmu, int nl)
-// See 'mMS2mRGI' but 'Alphas/Pi -> AlphaS*2*Beta0/Pi'
-double CRunDec::mMS2mRGImod(double mMS, double asmu, int nl){
- if(nl<0||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fSetcx((2*bet0*asmu)/Pi, nl);
- return (double) mMS/cAsmu;
-}
-
-// Coefficients of eq.(22) of [RunDec]
-double CRunDec::fas5to6os(double A, double mass, double mu, double nlq, double nl){
- double lmM = log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=A*(-lmM)/6.;
- sum[2]=A*A*(-7./24. - (19.*lmM)/24. + (lmM*lmM)/36.);
- sum[3]=A*A*A*(-58933./124416. - (8521.*lmM)/1728. - (131.*lmM*lmM)/576. -
- (lmM*lmM*lmM)/216. +
- nlq*(2479./31104. + (409.*lmM)/1728. + Zeta2/9.) - (2.*Zeta2)/3. -
- (2.*Zeta2*log(2))/9. -
- (80507.*Zeta3)/27648.);
- sum[4]=
- A*A*A*A*(-141841753./24494400. -
- (7693.*lmM*lmM)/1152. - (8371.*lmM*lmM*lmM)/10368. + lmM*lmM*lmM*lmM/1296. -
- (644201.*Pi*Pi)/116640. - (71102219.*Pi*Pi*Pi*Pi)/195955200. - (49.*Zeta2)/18. -
- (49*log(2)*Zeta2)/54 + (49*Zeta3)/216 + (587*Pi*Pi*log(2))/486 +
- (9318467*Pi*Pi*Pi*Pi*log(2))/32659200 - (2913037*Pi*Pi*log(2)*log(2))/1306368 +
- (340853*Pi*Pi*log(2)*log(2)*log(2))/816480 + (3179149*log(2)*log(2)*log(2)*log(2))/1306368 -
- (340853*log(2)*log(2)*log(2)*log(2)*log(2))/1360800 + (3179149*A4)/54432 +
- (340853*A5)/11340 + lmM*(-19696909./746496. - (29*Zeta2)/9 -
- (29*log(2)*Zeta2)/27 + (59*Zeta3)/108 - (2529743*Zeta3)/165888) +
- nlq*nlq*(-140825./1492992. - (493*lmM*lmM)/20736 - (13*Pi*Pi)/972 +
- lmM*(-1679./186624. - Zeta2/27) - (19*Zeta3)/1728) -
- (2428169183*Zeta3)/87091200 + (1439*Pi*Pi*Zeta3)/1296 +
- nlq*(1773073./746496. + (6661*lmM*lmM)/10368 + (107*lmM*lmM*lmM)/1728 +
- (967*Pi*Pi)/1944 - (697709*Pi*Pi*Pi*Pi)/14929920 + (49*Zeta2)/108 +
- (11*Pi*Pi*log(2))/243 - (1709*Pi*Pi*log(2)*log(2))/124416 +
- (173*log(2)*log(2)*log(2)*log(2))/124416 + (173*A4)/5184 +
- (4756441*Zeta3)/995328 + lmM*(1110443./373248. + (41*Zeta2)/54 +
- (2*log(2)*Zeta2)/27 + (7*Zeta3)/27 + (110779*Zeta3)/82944) +
- (115*Zeta5)/576) - (40596749*Zeta5)/1451520);
-
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecAsDownOS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecAsDownOS(double als, double massth, double muth, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl <<" LOOPS"<fas5to6os(als/Pi, massth, muth, Nf,nl));
- return als*erg;
-}
-
-// Coefficients of eq.(25) of [RunDec]
-double CRunDec::fas6to5os(double A, double mass, double mu, double nlq, double nl){
- double lmM=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=A*(lmM)/6.;
- sum[2]=A*A*(7./24. + (19.*lmM)/24. + (lmM*lmM)/36.);
- sum[3]=A*A*A*(58933./124416. + (2.*Zeta2)/3. + (2.*Zeta2*log(2))/9. +
- (80507.*Zeta3)/27648. + (8941.*lmM)/1728. +
- (511.*lmM*lmM)/576. + (lmM*lmM*lmM)/216. +
- nlq*(-2479./31104. - Zeta2/9. - (409.*lmM)/1728));
- sum[4]=A*A*A*A*
- (2313952./382725. + (47039.*lmM*lmM)/3456. + (14149.*lmM*lmM*lmM)/10368. +
- lmM*lmM*lmM*lmM/1296. + (644201.*Pi*Pi)/116640. + (71102219.*Pi*Pi*Pi*Pi)/195955200. +
- (49.*Zeta2)/18. + (49.*log(2)*Zeta2)/54. - (49.*Zeta3)/216. - (587.*Pi*Pi*log(2))/486. -
- (9318467.*Pi*Pi*Pi*Pi*log(2))/32659200. + (2913037.*Pi*Pi*log(2)*log(2))/1306368. -
- (340853.*Pi*Pi*log(2)*log(2)*log(2))/816480. - (3179149.*log(2)*log(2)*log(2)*log(2))/1306368. +
- (340853.*log(2)*log(2)*log(2)*log(2)*log(2))/1360800. - (3179149.*A4)/54432. -
- (340853.*A5)/11340. + nlq*nlq*(140825./1492992. +
- (493.*lmM*lmM)/20736. + (13.*Pi*Pi)/972. + lmM*(1679./186624. + Zeta2/27.) +
- (19.*Zeta3)/1728.) + (2428169183.*Zeta3)/87091200. -
- (1439.*Pi*Pi*Zeta3)/1296. + lmM*(21084715./746496. + (35.*Zeta2)/9. +
- (35.*log(2)*Zeta2)/27. - (65.*Zeta3)/108. + (3022001.*Zeta3)/165888.) +
- nlq*(-1773073./746496. - (9115.*lmM*lmM)/10368. - (107.*lmM*lmM*lmM)/1728. -
- (967.*Pi*Pi)/1944. + (697709.*Pi*Pi*Pi*Pi)/14929920. - (49.*Zeta2)/108. -
- (11.*Pi*Pi*log(2))/243. + (1709.*Pi*Pi*log(2)*log(2))/124416. -
- (173.*log(2)*log(2)*log(2)*log(2))/124416. - (173.*A4)/5184. +
- lmM*(-1140191./373248. - (47.*Zeta2)/54. - (2.*log(2)*Zeta2)/27. - (7.*Zeta3)/27. -
- (110779.*Zeta3)/82944.) - (4756441.*Zeta3)/995328. -
- (115.*Zeta5)/576.) + (40596749.*Zeta5)/1451520.);
-
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecAsUpOS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecAsUpOS(double als, double massth, double muth, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl <<" LOOPS"<fas6to5os(als/Pi, massth, muth, Nf,nl));
- return als*erg;
-}
-
-
-double CRunDec::fas6to5ms(double A, double mass, double mu, double nlq, double nl){
- double lmm=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]= 1.;
- sum[1]=A*(lmm)/6.;
- sum[2]=A*A*(-11./72. + (11.*lmm)/24. + lmm*lmm/36.);
- sum[3]=A*A*A*(-564731./124416. + (2645.*lmm)/1728. + (167.*lmm*lmm)/576. + lmm*lmm*lmm/216. +
- (2633./31104. - (67.*lmm)/576. + lmm*lmm/36.)*nlq + (82043.*Zeta3)/27648.);
- sum[4]=A*A*A*A*(-1165152397./24494400. - (3031309*A4)/54432 - (340853*A5)/11340
- +(1837*lmm*lmm)/1152 + (2909*lmm*lmm*lmm)/10368 + lmm*lmm*lmm*lmm/1296
- -(3031309*log(2)*log(2)*log(2)*log(2))/1306368 + (340853*log(2)*log(2)*log(2)*log(2)*log(2))/1360800
- +(3031309*log(2)*log(2)*Pi*Pi)/1306368 - (340853*log(2)*log(2)*log(2)*Pi*Pi)/816480
- +(76940219*Pi*Pi*Pi*Pi)/195955200 - (9318467*log(2)*Pi*Pi*Pi*Pi)/32659200 + nlq*nlq*(271883./4478976.
- -(6865*lmm)/186624 + (77*lmm*lmm)/20736 - lmm*lmm*lmm/324 - (167*Zeta3)/5184)
- +(2362581983*Zeta3)/87091200 + lmm*(-11093717./746496. + (3022001*Zeta3)/165888) + nlq*(4770941./2239488.
- -(685*A4)/5184 + (277*lmm*lmm)/10368 + (271*lmm*lmm*lmm)/5184 - (685*log(2)*log(2)*log(2)*log(2))/124416
- +(685*log(2)*log(2)*Pi*Pi)/124416 + (541549*Pi*Pi*Pi*Pi)/14929920 + lmm*(141937./373248.
- -(110779*Zeta3)/82944) - (3645913*Zeta3)/995328 - (115*Zeta5)/576) + (12057583*Zeta5)/483840);
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecAsUpMS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecAsUpMS(double als, double massth, double muth, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl <<" LOOPS"<fas6to5ms(als/Pi, massth, muth, Nf,nl));
- return als*erg;
-}
-
-double CRunDec::fas5to6ms(double A, double mass, double mu, double nlq, double nl){
- double lmm=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]= 1.;
- sum[1]=A*(-lmm)/6.;
- sum[2]=A*A*(11./72. - (11*lmm)/24 + lmm*lmm/36.);
- sum[3]=A*A*A*(564731./124416. - (955.*lmm)/576. + (53.*lmm*lmm)/576. - lmm*lmm*lmm/216.
- +(-2633./31104. + (67.*lmm)/576. - lmm*lmm/36.)*nlq - (82043.*Zeta3)/27648.);
- sum[4]=A*A*A*A*(291716893./6123600. + (3031309*A4)/54432 + (340853*A5)/11340
- +(2177*lmm*lmm)/3456 - (1883*lmm*lmm*lmm)/10368 + lmm*lmm*lmm*lmm/1296
- +(3031309*log(2)*log(2)*log(2)*log(2))/1306368
- -(340853*log(2)*log(2)*log(2)*log(2)*log(2))/1360800
- -(3031309*log(2)*log(2)*Pi*Pi)/1306368 + (340853*log(2)*log(2)*log(2)*Pi*Pi)/816480
- -(76940219*Pi*Pi*Pi*Pi)/195955200 + (9318467*log(2)*Pi*Pi*Pi*Pi)/32659200
- +lmm*(7391699./746496. - (2529743*Zeta3)/165888) + nlq*nlq*(-271883./4478976.
- +(6865*lmm)/186624 - (77*lmm*lmm)/20736 + lmm*lmm*lmm/324 + (167*Zeta3)/5184)
- -(2362581983*Zeta3)/87091200 + nlq*(-4770941./2239488. + (685*A4)/5184
- -(1483*lmm*lmm)/10368 - (127*lmm*lmm*lmm)/5184 + (685*log(2)*log(2)*log(2)*log(2))/124416
- -(685*log(2)*log(2)*Pi*Pi)/124416 - (541549*Pi*Pi*Pi*Pi)/14929920 + (3645913*Zeta3)/995328
- +lmm*(-110341./373248. + (110779*Zeta3)/82944) + (115*Zeta5)/576) - (12057583*Zeta5)/483840);
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecAsDownMS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecAsDownMS(double als, double massth, double muth, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl <<" LOOPS"<fas5to6ms(als/Pi, massth, muth, Nf,nl));
- return als*erg;
-}
-
-
-double CRunDec::fas6to5si(double A, double mass, double mu, double nlq, double nl){
- double lmmu=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0] = 1.0;
- sum[1] = (A*lmmu)/6.;
- sum[2] = A*A*(-11./72. + (19.*lmmu)/24. + lmmu*lmmu/36.);
- sum[3] = A*A*A*(-564731./124416. + (2191.*lmmu)/576. + (511.*lmmu*lmmu)/576. + lmmu*lmmu*lmmu/216. +
- (2633./31104. - (281.*lmmu)/1728.)*nlq + (82043.*Zeta3)/27648.);
- sum[4] = A*A*A*A*(-1165152397./24494400. - (3031309*A4)/54432 - (340853*A5)/11340 + (33887*lmmu*lmmu)/3456
- +(14149*lmmu*lmmu*lmmu)/10368 + lmmu*lmmu*lmmu*lmmu/1296
- -(3031309*log(2)*log(2)*log(2)*log(2))/1306368 + (340853*log(2)*log(2)*log(2)*log(2)*log(2))/1360800
- +(3031309*log(2)*log(2)*Pi*Pi)/1306368 - (340853*log(2)*log(2)*log(2)*Pi*Pi)/816480
- +(76940219*Pi*Pi*Pi*Pi)/195955200 - (9318467*log(2)*Pi*Pi*Pi*Pi)/32659200 + nlq*nlq*(271883./4478976.
- -(8545*lmmu)/186624 + (79*lmmu*lmmu)/6912 - (167*Zeta3)/5184) + (2362581983*Zeta3)/87091200
- +lmmu*(-1531493./746496. + (2975921*Zeta3)/165888) + nlq*(4770941./2239488. - (685*A4)/5184
- -(515*lmmu*lmmu)/1152 - (107*lmmu*lmmu*lmmu)/1728 - (685*log(2)*log(2)*log(2)*log(2))/124416
- +(685*log(2)*log(2)*Pi*Pi)/124416 + (541549*Pi*Pi*Pi*Pi)/14929920 + lmmu*(-158687./373248.
- -(133819*Zeta3)/82944) - (3645913*Zeta3)/995328 - (115*Zeta5)/576) + (12057583*Zeta5)/483840);
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecAsUpSI(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecAsUpSI(double als, double massth, double muth, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl <<" LOOPS"<fas6to5si(als/Pi, massth, muth, Nf,nl));
- return als*erg;
-}
-
-double CRunDec::fas5to6si(double A, double mass, double mu, double nlq, double nl){
- double lmmu=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0] = 1;
- sum[1] = -(A*lmmu)/6.;
- sum[2] = A*A*(11./72. - (19.*lmmu)/24. + lmmu*lmmu/36.);
- sum[3] = A*A*A*(564731./124416. - (6793.*lmmu)/1728. - (131.*lmmu*lmmu)/576. - lmmu*lmmu*lmmu/216. +
- (-2633./31104. + (281.*lmmu)/1728.)*nlq - (82043.*Zeta3)/27648.);
- sum[4] = A*A*A*A*(291716893./6123600. + (3031309*A4)/54432 + (340853*A5)/11340 - (14023*lmmu*lmmu)/3456
- -(8371*lmmu*lmmu*lmmu)/10368 + lmmu*lmmu*lmmu*lmmu/1296
- +(3031309*log(2)*log(2)*log(2)*log(2))/1306368 - (340853*log(2)*log(2)*log(2)*log(2)*log(2))/1360800
- -(3031309*log(2)*log(2)*Pi*Pi)/1306368 + (340853*log(2)*log(2)*log(2)*Pi*Pi)/816480
- -(76940219*Pi*Pi*Pi*Pi)/195955200 + (9318467*log(2)*Pi*Pi*Pi*Pi)/32659200 + lmmu*(-2398621./746496.
- -(2483663*Zeta3)/165888) + nlq*nlq*(-271883./4478976. + (8545*lmmu)/186624 - (79*lmmu*lmmu)/6912
- +(167*Zeta3)/5184) - (2362581983*Zeta3)/87091200 + nlq*(-4770941./2239488. + (685*A4)/5184
- +(983*lmmu*lmmu)/3456 + (107*lmmu*lmmu*lmmu)/1728 + (685*log(2)*log(2)*log(2)*log(2))/124416
- -(685*log(2)*log(2)*Pi*Pi)/124416 - (541549*Pi*Pi*Pi*Pi)/14929920 + (3645913*Zeta3)/995328
- +lmmu*(190283./373248. + (133819*Zeta3)/82944) + (115*Zeta5)/576) - (12057583*Zeta5)/483840);
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecAsDownSI(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecAsDownSI(double als, double massth, double muth, int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl <<" LOOPS"<fas5to6si(als/Pi, massth, muth, Nf,nl));
- return als*erg;
-}
-
-
-// Coefficients of eq.(33) of [RunDec]
-double CRunDec::fmq6to5os(double A, double mass, double mu, double nlq, double nl){
- double lmM=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=0.;
- sum[2]=-A*A*(89./432. - (5.*lmM)/36. + lmM*lmM/12.);
- sum[3]=-A*A*A*(1871./2916. - (4.*A4)/9. - log(2)*log(2)*log(2)*log(2)/54. + (299.*lmM)/2592.
- + (299.*lmM*lmM)/432. + (35.*lmM*lmM*lmM)/216. + (1327.*nlq)/11664.
- - (53.*lmM*nlq)/432. - (lmM*lmM*lmM*nlq)/108. + (log(2)*log(2)*Zeta2)/9. - (407.*Zeta3)/864.
- - (5.*lmM*Zeta3)/6. - (2.*nlq*Zeta3)/27. + (103.*Zeta4)/72.);
- sum[4]=A*A*A*A*(-122722873709./3292047360. - (6941*lmM*lmM*lmM)/2592 -
- (1147*lmM*lmM*lmM*lmM)/3456 - (16187201*Pi*Pi*Pi*Pi)/52254720 + (787*Pi*Pi*Pi*Pi*Pi*Pi)/81648 +
- (5*Zeta2)/9 + (5*log(2)*Zeta2)/27 - (5*Zeta3)/108 + (145*Pi*Pi*Pi*Pi*log(2))/1944 -
- (1924649*Pi*Pi*log(2)*log(2))/4354560 + (59*Pi*Pi*log(2)*log(2)*log(2))/972 +
- (1924649*log(2)*log(2)*log(2)*log(2))/4354560 - (59*log(2)*log(2)*log(2)*log(2)*log(2))/1620 - (4*log(2)*log(2)*log(2)*log(2)*log(2)*log(2))/81 +
- (log(2)*log(2)*log(2)*log(2)*log(2)*log(16))/81 + (1924649*A4)/181440 +
- (118*A5)/27 + nlq*nlq*(-17671./124416. - (31*lmM*lmM)/1296 -
- lmM*lmM*lmM*lmM/864 + (7*Pi*Pi*Pi*Pi)/8640 + lmM*(3401./46656. - (7*Zeta3)/108) +
- (5*Zeta3)/864) + (443509931*Zeta3)/40642560 -
- (1061*Zeta3*Zeta3)/576 + lmM*lmM*(-34297./3456. + (175*Zeta3)/48) +
- lmM*(99319./41472. - (481*Pi*Pi*Pi*Pi)/2880 - (2*Zeta2)/3 - (2*log(2)*Zeta2)/9 +
- Zeta3/18 - (11*Pi*Pi*log(2)*log(2))/72 + (11*log(2)*log(2)*log(2)*log(2))/72 +
- (11*A4)/3 + (47317*Zeta3)/3072 - (575*Zeta5)/72) +
- nlq*(2403419./746496. + (10237*lmM*lmM)/10368 + (47*lmM*lmM*lmM)/288 +
- (17*lmM*lmM*lmM*lmM)/432 + (245*Pi*Pi*Pi*Pi)/62208 - (5*Zeta2)/54 -
- (49*Pi*Pi*Pi*Pi*log(2))/6480 + (49*Pi*Pi*log(2)*log(2))/2592 -
- (Pi*Pi*log(2)*log(2)*log(2))/162 - (49*log(2)*log(2)*log(2)*log(2))/2592 + log(2)*log(2)*log(2)*log(2)*log(2)/270 -
- (49*A4)/108 - (4*A5)/9 +
- lmM*(-26443./93312. + (163*Pi*Pi*Pi*Pi)/12960 + Zeta2/9 + (Pi*Pi*log(2)*log(2))/108 -
- log(2)*log(2)*log(2)*log(2)/108 - (2*A4)/9 - (599*Zeta3)/1728) +
- (1075*Zeta3)/1728 - (497*Zeta5)/288) + (59015*Zeta5)/1728);
-
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-double CRunDec::fmq6to5ms(double A, double mass, double mu, double nlq, double nl){
- double lmm=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=0.;
- sum[2]=A*A*(-89./432. + (5*lmm)/36 - lmm*lmm/12);
- sum[3]=A*A*A*(-2951./2916. + (16*A4+2./3.*log(2)*log(2)*log(2)*log(2)-2./3.*Pi*Pi*log(2)*log(2)-13./180.*Pi*Pi*Pi*Pi)/36
- -(155*lmm*lmm)/432 - (35*lmm*lmm*lmm)/216 - Pi*Pi*Pi*Pi/72
- +nl*(-1327./11664. + (53*lmm)/432 + lmm*lmm*lmm/108 + (2*Zeta3)/27)
- +lmm*(133./2592. + (5*Zeta3)/6) + (407*Zeta3)/864);
- sum[4]=A*A*A*A*(-131621265869./3292047360. + (1924649*A4)/181440 + (118*A5)/27
- -(3161*lmm*lmm*lmm)/2592 - (1147*lmm*lmm*lmm*lmm)/3456
- +(1924649*log(2)*log(2)*log(2)*log(2))/4354560 - (59*log(2)*log(2)*log(2)*log(2)*log(2))/1620
- +(log(16)*log(2)*log(2)*log(2)*log(2)*log(2))/81 - (4*log(2)*log(2)*log(2)*log(2)*log(2)*log(2))/81
- -(1924649*log(2)*log(2)*Pi*Pi)/4354560 + (59*log(2)*log(2)*log(2)*Pi*Pi)/972
- -(16187201*Pi*Pi*Pi*Pi)/52254720 + (145*log(2)*Pi*Pi*Pi*Pi)/1944 + (787*Pi*Pi*Pi*Pi*Pi*Pi)/81648
- +nlq*nlq*(-17671./124416. - (31*lmm*lmm)/1296 - lmm*lmm*lmm*lmm/864 + (7*Pi*Pi*Pi*Pi)/8640
- +lmm*(3401./46656. - (7*Zeta3)/108) + (5*Zeta3)/864) + (353193131*Zeta3)/40642560
- -(1061*Zeta3*Zeta3)/576 + lmm*lmm*(-16193./3456. + (175*Zeta3)/48) + lmm*(279367./41472.
- +(11*A4)/3 + (11*log(2)*log(2)*log(2)*log(2))/72 - (11*log(2)*log(2)*Pi*Pi)/72
- -(481*Pi*Pi*Pi*Pi)/2880 + (42197*Zeta3)/3072 - (575*Zeta5)/72) + nlq*(2261435./746496.
- -(49*A4)/108 - (4*A5)/9 + (8461*lmm*lmm)/10368 + (23*lmm*lmm*lmm)/288 + (17*lmm*lmm*lmm*lmm)/432
- -(49*log(2)*log(2)*log(2)*log(2))/2592 + log(2)*log(2)*log(2)*log(2)*log(2)/270
- +(49*log(2)*log(2)*Pi*Pi)/2592 - (log(2)*log(2)*log(2)*Pi*Pi)/162 + (245*Pi*Pi*Pi*Pi)/62208
- -(49*log(2)*Pi*Pi*Pi*Pi)/6480 + lmm*(-55315./93312. - (2*A4)/9 - log(2)*log(2)*log(2)*log(2)/108
- +(log(2)*log(2)*Pi*Pi)/108 + (163*Pi*Pi*Pi*Pi)/12960 - (599*Zeta3)/1728)
- +(1075*Zeta3)/1728 - (497*Zeta5)/288) + (59015*Zeta5)/1728);
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-double CRunDec::fmq6to5si(double A, double mass, double mu, double nlq, double nl){
- double lmmu=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=0.;
- sum[2]=A*A*(-89./432. + (5*lmmu)/36 - lmmu*lmmu/12);
- sum[3]=A*A*A*(-2951./2916. - (299*lmmu*lmmu)/432 - (35*lmmu*lmmu*lmmu)/216 - Pi*Pi*Pi*Pi/72 + (16*A4
- +(2*log(2)*log(2)*log(2)*log(2))/3 - (2*log(2)*log(2)*Pi*Pi)/3 - (13*Pi*Pi*Pi*Pi)/180)/36
- +nlq*(-1327./11664. + (53*lmmu)/432 + lmmu*lmmu*lmmu/108 + (2*Zeta3)/27) + lmmu*(853./2592.
- +(5*Zeta3)/6) + (407*Zeta3)/864);
- sum[4]=A*A*A*A*(-131621265869./3292047360. + (1924649*A4)/181440 + (118*A5)/27 - (6941*lmmu*lmmu*lmmu)/2592
- -(1147*lmmu*lmmu*lmmu*lmmu)/3456 + (1924649*log(2)*log(2)*log(2)*log(2))/4354560
- -(59*log(2)*log(2)*log(2)*log(2)*log(2))/1620 + (log(16)*log(2)*log(2)*log(2)*log(2)*log(2))/81
- -(4*log(2)*log(2)*log(2)*log(2)*log(2)*log(2))/81 - (1924649*log(2)*log(2)*Pi*Pi)/4354560
- +(59*log(2)*log(2)*log(2)*Pi*Pi)/972 - (16187201*Pi*Pi*Pi*Pi)/52254720
- +(145*log(2)*Pi*Pi*Pi*Pi)/1944 + (787*Pi*Pi*Pi*Pi*Pi*Pi)/81648 + nlq*nlq*(-17671./124416.
- -(31*lmmu*lmmu)/1296 - lmmu*lmmu*lmmu*lmmu/864 + (7*Pi*Pi*Pi*Pi)/8640 + lmmu*(3401./46656.
- -(7*Zeta3)/108) + (5*Zeta3)/864) + (353193131*Zeta3)/40642560 - (1061*Zeta3*Zeta3)/576
- +lmmu*lmmu*(-8531./1152. + (175*Zeta3)/48) + lmmu*(330503./41472. + (11*A4)/3
- +(11*log(2)*log(2)*log(2)*log(2))/72 - (11*log(2)*log(2)*Pi*Pi)/72 - (481*Pi*Pi*Pi*Pi)/2880
- +(47317*Zeta3)/3072 - (575*Zeta5)/72) + nlq*(2261435./746496. - (49*A4)/108 - (4*A5)/9
- +(8701*lmmu*lmmu)/10368 + (47*lmmu*lmmu*lmmu)/288 + (17*lmmu*lmmu*lmmu*lmmu)/432
- -(49*log(2)*log(2)*log(2)*log(2))/2592 + log(2)*log(2)*log(2)*log(2)*log(2)/270
- +(49*log(2)*log(2)*Pi*Pi)/2592 - (log(2)*log(2)*log(2)*Pi*Pi)/162 + (245*Pi*Pi*Pi*Pi)/62208
- -(49*log(2)*Pi*Pi*Pi*Pi)/6480 + lmmu*(-36019./93312. - (2*A4)/9 - log(2)*log(2)*log(2)*log(2)/108
- +(log(2)*log(2)*Pi*Pi)/108 + (163*Pi*Pi*Pi*Pi)/12960 - (599*Zeta3)/1728) + (1075*Zeta3)/1728
- -(497*Zeta5)/288) + (59015*Zeta5)/1728);
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecMqUpOS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecMqUpOS(double mq, double als, double massth, double muth,
- int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fmq6to5os(als/Pi, massth, muth, Nf,nl));
- return mq*erg;
-}
-
-// Function double CRunDec::DecMqUpMS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecMqUpMS(double mq, double als, double massth, double muth,
- int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fmq6to5ms(als/Pi, massth, muth, Nf,nl));
- return mq*erg;
-}
-
-// Function double CRunDec::DecMqUpSI(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecMqUpSI(double mq, double als, double massth, double muth,
- int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fmq6to5si(als/Pi, massth, muth, Nf,nl));
- return mq*erg;
-}
-
-// Coefficients of eq.(30) of [RunDec]
-double CRunDec::fmq5to6os(double A, double mass, double mu, double nlq, double nl){
- double lmM=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=0.;
- sum[2]=A*A*(89./432. - (5.*lmM)/36. + lmM*lmM/12.);
- sum[3]=A*A*A*(1871./2916. - (4.*A4)/9. - ln2*ln2*ln2*ln2/54. + (121.*lmM)/2592.
- +(319.*lmM*lmM)/432. + (29.*lmM*lmM*lmM)/216. + (1327.*nlq)/11664.
- -(53.*lmM*nlq)/432. - (lmM*lmM*lmM*nlq)/108. + (ln2*ln2*Zeta2)/9.
- -(407.*Zeta3)/864. - (5.*lmM*Zeta3)/6. - (2.*nlq*Zeta3)/27. + (103.*Zeta4)/72.);
- sum[4]=A*A*A*A*(122466970229./3292047360. - (1924649.*A4)/181440. - (118.*A5)/27.
- -(1924649.*ln2*ln2*ln2*ln2)/4354560. + (59.*ln2*ln2*ln2*ln2*ln2)/1620.
- -(1126487.*lmM)/373248. - (31.*A4*lmM)/9.-(31.*ln2*ln2*ln2*ln2*lmM)/216.
- +(104803.*lmM*lmM)/10368. + (1403.*lmM*lmM*lmM)/648. + (305.*lmM*lmM*lmM*lmM)/1152.
- -(2403419.*nlq)/746496. + (49.*A4*nlq)/108. + (4.*A5*nlq)/9. + (49.*ln2*ln2*ln2*ln2*nlq)/2592.
- -(ln2*ln2*ln2*ln2*ln2*nlq)/270. + (7045.*lmM*nlq)/31104. + (2.*A4*lmM*nlq)/9.
- +1./108.*ln2*ln2*ln2*ln2*lmM*nlq - (9601.*lmM*lmM*nlq)/10368. - (47.*lmM*lmM*lmM*nlq)/288.
- -(5.*lmM*lmM*lmM*lmM*nlq)/144. + (17671.*nlq*nlq)/124416. - (3401.*lmM*nlq*nlq)/46656.
- +(31.*lmM*lmM*nlq*nlq)/1296. + (lmM*lmM*lmM*lmM*nlq*nlq)/864. - (5.*Zeta2)/9.
- -(5.*ln2*Zeta2)/27. + (1924649.*ln2*ln2*Zeta2)/725760. - (59.*ln2*ln2*ln2*Zeta2)/162.
- +(2.*lmM*Zeta2)/3. + (2.*ln2*lmM*Zeta2)/9. + 31./36.*ln2*ln2*lmM*Zeta2 + (5.*nlq*Zeta2)/54.
- -49./432.*ln2*ln2*nlq*Zeta2 + 1./27.*ln2*ln2*ln2*nlq*Zeta2 - (lmM*nlq*Zeta2)/9.
- - 1./18.*ln2*ln2*lmM*nlq*Zeta2 - (441628331.*Zeta3)/40642560. - (420877.*lmM*Zeta3)/27648.
- -(155.*lmM*lmM*Zeta3)/48. - (1075.*nlq*Zeta3)/1728. + (221.*lmM*nlq*Zeta3)/576.
- -(5.*nlq*nlq*Zeta3)/864. + 7./108.*lmM*nlq*nlq*Zeta3 + (1061.*Zeta3*Zeta3)/576.
- +(16187201.*Zeta4)/580608. - (725.*ln2*Zeta4)/108. + (4123.*lmM*Zeta4)/288.
- -(1225.*nlq*Zeta4)/3456. + (49.*ln2*nlq*Zeta4)/72. - (163.*lmM*nlq*Zeta4)/144.
- -(7.*nlq*nlq*Zeta4)/96. - (59015.*Zeta5)/1728. + (575.*lmM*Zeta5)/72.
- +(497.*nlq*Zeta5)/288. - (3935.*Zeta6)/432.);
-
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-double CRunDec::fmq5to6ms(double A, double mass, double mu, double nlq, double nl){
- double lmm=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=0.;
- sum[2]=A*A*(89./432. - (5.*lmm)/36. + lmm*lmm/12.);
- sum[3]=A*A*A*(2951./2916. - (16*A4+2./3.*log(2)*log(2)*log(2)*log(2)-2./3.*Pi*Pi*log(2)*log(2)-13./180.*Pi*Pi*Pi*Pi)/36
- +(175*lmm*lmm)/432 + (29*lmm*lmm*lmm)/216 + Pi*Pi*Pi*Pi/72
- +lmm*(-311./2592. - (5*Zeta3)/6)
- +nlq*(1327./11664. - (53*lmm)/432 - lmm*lmm*lmm/108 - (2*Zeta3)/27) - (407*Zeta3)/864);
- sum[4]=A*A*A*A*(131968227029./3292047360. - (1924649*A4)/181440 - (118*A5)/27
- +(301*lmm*lmm*lmm)/324 + (305*lmm*lmm*lmm*lmm)/1152 - (1924649*log(2)*log(2)*log(2)*log(2))/4354560
- +(59*log(2)*log(2)*log(2)*log(2)*log(2))/1620 - (log(16)*log(2)*log(2)*log(2)*log(2)*log(2))/81
- +(4*log(2)*log(2)*log(2)*log(2)*log(2)*log(2))/81 + (1924649*log(2)*log(2)*Pi*Pi)/4354560
- -(59*log(2)*log(2)*log(2)*Pi*Pi)/972 + (16187201*Pi*Pi*Pi*Pi)/52254720
- -(145*log(2)*Pi*Pi*Pi*Pi)/1944 - (787*Pi*Pi*Pi*Pi*Pi*Pi)/81648 + lmm*lmm*(51163./10368.
- -(155*Zeta3)/48) + nlq*nlq*(17671./124416. + (31*lmm*lmm)/1296 + lmm*lmm*lmm*lmm/864
- -(7*Pi*Pi*Pi*Pi)/8640 + lmm*(-3401./46656. + (7*Zeta3)/108) - (5*Zeta3)/864)
- -(353193131*Zeta3)/40642560 + (1061*Zeta3*Zeta3)/576 - (59015*Zeta5)/1728
- +nlq*(-2261435./746496. + (49*A4)/108 + (4*A5)/9 - (7825*lmm*lmm)/10368 - (23*lmm*lmm*lmm)/288
- -(5*lmm*lmm*lmm*lmm)/144 + (49*log(2)*log(2)*log(2)*log(2))/2592
- -log(2)*log(2)*log(2)*log(2)*log(2)/270 - (49*log(2)*log(2)*Pi*Pi)/2592
- +(log(2)*log(2)*log(2)*Pi*Pi)/162 - (245*Pi*Pi*Pi*Pi)/62208 + (49*log(2)*Pi*Pi*Pi*Pi)/6480
- +lmm*(16669./31104. + (2*A4)/9 + log(2)*log(2)*log(2)*log(2)/108 - (log(2)*log(2)*Pi*Pi)/108
- -(163*Pi*Pi*Pi*Pi)/12960 + (221*Zeta3)/576) - (1075*Zeta3)/1728 + (497*Zeta5)/288)
- +lmm*(-2810855./373248. - (31*A4)/9 - (31*log(2)*log(2)*log(2)*log(2))/216
- +(31*log(2)*log(2)*Pi*Pi)/216 + (4123*Pi*Pi*Pi*Pi)/25920 - (373261*Zeta3)/27648 + (575*Zeta5)/72));
-
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-double CRunDec::fmq5to6si(double A, double mass, double mu, double nlq, double nl){
- double lmmu=log((mu*mu)/(mass*mass));
- if(nl==1){
- return 1;
- }
- double sum[5];
- sum[0]=1.;
- sum[1]=0.;
- sum[2]=A*A*(89./432. - (5*lmmu)/36 + lmmu*lmmu/12);
- sum[3]=A*A*A*(2951./2916. + (319*lmmu*lmmu)/432 + (29*lmmu*lmmu*lmmu)/216 + Pi*Pi*Pi*Pi/72 + (-16*A4
- -(2*log(2)*log(2)*log(2)*log(2))/3 + (2*log(2)*log(2)*Pi*Pi)/3 + (13*Pi*Pi*Pi*Pi)/180)/36
- +lmmu*(-1031./2592. - (5*Zeta3)/6) + nlq*(1327./11664. - (53*lmmu)/432 - lmmu*lmmu*lmmu/108
- -(2*Zeta3)/27) - (407*Zeta3)/864);
- sum[4]=A*A*A*A*(131968227029./3292047360. - (1924649*A4)/181440 - (118*A5)/27 + (1403*lmmu*lmmu*lmmu)/648
- +(305*lmmu*lmmu*lmmu*lmmu)/1152 - (1924649*log(2)*log(2)*log(2)*log(2))/4354560
- +(59*log(2)*log(2)*log(2)*log(2)*log(2))/1620 - (log(16)*log(2)*log(2)*log(2)*log(2)*log(2))/81
- +(4*log(2)*log(2)*log(2)*log(2)*log(2)*log(2))/81 + (1924649*log(2)*log(2)*Pi*Pi)/4354560
- -(59*log(2)*log(2)*log(2)*Pi*Pi)/972 + (16187201*Pi*Pi*Pi*Pi)/52254720 - (145*log(2)*Pi*Pi*Pi*Pi)/1944
- -(787*Pi*Pi*Pi*Pi*Pi*Pi)/81648 + lmmu*lmmu*(81763./10368. - (155*Zeta3)/48) + nlq*nlq*(17671./124416.
- +(31*lmmu*lmmu)/1296 + lmmu*lmmu*lmmu*lmmu/864 - (7*Pi*Pi*Pi*Pi)/8640 + lmmu*(-3401./46656.
- +(7*Zeta3)/108) - (5*Zeta3)/864) - (353193131*Zeta3)/40642560 + (1061*Zeta3*Zeta3)/576
- -(59015*Zeta5)/1728 + nlq*(-2261435./746496. + (49*A4)/108 + (4*A5)/9 - (8065*lmmu*lmmu)/10368
- -(47*lmmu*lmmu*lmmu)/288 - (5*lmmu*lmmu*lmmu*lmmu)/144 + (49*log(2)*log(2)*log(2)*log(2))/2592
- -log(2)*log(2)*log(2)*log(2)*log(2)/270 - (49*log(2)*log(2)*Pi*Pi)/2592 + (log(2)*log(2)*log(2)*Pi*Pi)/162
- -(245*Pi*Pi*Pi*Pi)/62208 + (49*log(2)*Pi*Pi*Pi*Pi)/6480 + lmmu*(10237./31104. + (2*A4)/9
- +log(2)*log(2)*log(2)*log(2)/108 - (log(2)*log(2)*Pi*Pi)/108 - (163*Pi*Pi*Pi*Pi)/12960 + (221*Zeta3)/576)
- -(1075*Zeta3)/1728 + (497*Zeta5)/288) + lmmu*(-3322343./373248. - (31*A4)/9
- -(31*log(2)*log(2)*log(2)*log(2))/216 + (31*log(2)*log(2)*Pi*Pi)/216 + (4123*Pi*Pi*Pi*Pi)/25920
- -(419341*Zeta3)/27648 + (575*Zeta5)/72));
-
- double erg=0.0;
- for(int i=1;i<=nl;i++){
- erg+=sum[i-1];
- }
- return erg;
-}
-
-// Function double CRunDec::DecMqDownOS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecMqDownOS(double mq, double als, double massth, double muth,
- int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fmq5to6os(als/Pi, massth, muth, Nf,nl));
- return mq*erg;
-}
-
-// Function double CRunDec::DecMqDownMS(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecMqDownMS(double mq, double als, double massth, double muth,
- int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fmq5to6ms(als/Pi, massth, muth, Nf,nl));
- return mq*erg;
-}
-
-// Function double CRunDec::DecMqDownSI(double als, double massth, double muth,
-// int nl)
-double CRunDec::DecMqDownSI(double mq, double als, double massth, double muth,
- int nl){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<fmq5to6si(als/Pi, massth, muth, Nf,nl));
- return mq*erg;
-}
-
-// Function double CRunDec::DecLambdaDown(double lam, double massth, int nl,
-// int nloops)
-double CRunDec::DecLambdaDown(double lam, double massth, int nl,
- int nloops){
- if(nl<1||nl>5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<5){
- cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<decpar[i+1].nf){
- help=decpar[i].nf;
- decpar[i].nf=decpar[i+1].nf;
- decpar[i+1].nf=help;
- help2=decpar[i].Mth;
- decpar[i].Mth=decpar[i+1].Mth;
- decpar[i+1].Mth=help2;
- help2=decpar[i].muth;
- decpar[i].muth=decpar[i+1].muth;
- decpar[i+1].muth=help2;
- }
- }
- zaehler-=1;
- }
- for(int i=3-n+2;i<=3;i++){
- if(decpar[i].nf-decpar[i-1].nf!=1){
- cout<<"WARNING: THERE IS A GAP IN NUMBER OF FLAVOURS. EXIT."<decpar[i+1].nf){
- help=decpar[i].nf;
- decpar[i].nf=decpar[i+1].nf;
- decpar[i+1].nf=help;
- help2=decpar[i].Mth;
- decpar[i].Mth=decpar[i+1].Mth;
- decpar[i+1].Mth=help2;
- help2=decpar[i].muth;
- decpar[i].muth=decpar[i+1].muth;
- decpar[i+1].muth=help2;
- }
- }
- zaehler-=1;
- }
- for(int i=3-n+2;i<=3;i++){
- if(decpar[i].nf-decpar[i-1].nf!=1){
- cout<<"WARNING: THERE IS A GAP IN NUMBER OF FLAVOURS. EXIT.";
- RETURN
- }
- }
-
- double erg1,erg2,erg3, erg4, erg5, erg6;
- int i;
- for(i=3-n+1;i<=3;i++){
- erg1= AlphasExact(asini, muini, decpar[i].muth,decpar[i].nf-1,nl);
- erg3= mMS2mMS(mqini,asini,erg1,decpar[i].nf-1,nl);
- erg2= DecAsUpOS(erg1,decpar[i].Mth,decpar[i].muth,decpar[i].nf-1,nl);
- erg4=DecMqUpOS(erg3,erg1,decpar[i].Mth,decpar[i].muth,decpar[i].nf-1,nl);
- asini=erg2;
- mqini=erg4;
- muini=decpar[i].muth;
- }
- erg5= AlphasExact(asini,muini,mu2,decpar[i-1].nf,nl);
- erg6= mMS2mMS(mqini,asini,erg5,decpar[i-1].nf,nl);
- for(int j=0;j<=3;j++){
- decpar[j].nf=0;
- decpar[j].Mth=0.;
- decpar[j].muth=0.;
- }
- return erg6;
-}
-
-// Function double CRunDec::mH2mL(double mqh, double ash, double mu1,
-// TriplenfMmu decpar[], double mu2, int nl)
-double CRunDec::mH2mL(double mqh, double ash, double mu1, TriplenfMmu decpar[],
- double mu2, int nl){
- int n=0;
- int help;
- double help2;
- double asini=ash;
- double muini=mu1;
- double mqini=mqh;
- for(int i=0; i<4; i++){
- if(decpar[i].nf!=0){
- n+=1;
- }
- }
- int zaehler=4;
- while (zaehler!=0){
- for(int i=0; iLamExpl(AlphaS, Mu, nl));
-}
-
-double CRunDec::LamImpl(double AlphaS, double Mu,int nf,int nl){
- SetConstants(nf);
- return (this->LamImpl(AlphaS, Mu, nl));
-}
-
-double CRunDec::AlphasLam(double Lambda, double Mu,int nf, int nl){
- SetConstants(nf);
- return (this->AlphasLam(Lambda, Mu, nl));
-}
-
-double CRunDec::AlphasExact(double AlphaS0, double Mu0, double MuEnd,
- int nf,int nl){
- SetConstants(nf);
- return (this->AlphasExact(AlphaS0,Mu0,MuEnd,nl));
-}
-
-double CRunDec::mMS2mMS(double mmu0, double AlphaSEnd, double AlphaS0,int nf,
- int nl){
- SetConstants(nf);
- return (this->mMS2mMS(mmu0, AlphaSEnd, AlphaS0, nl));
-}
-
-AsmMS CRunDec::AsmMSrunexact(double mMu, double AlphaS0, double Mu0,
- double MuEnd,int nf, int nl){
- SetConstants(nf);
- return (this->AsmMSrunexact(mMu, AlphaS0, Mu0, MuEnd, nl));
-}
-
-double CRunDec::mMS2mOS(double mMS, std::pair* mq, double asmu, double mu,int nf,
- int nl, double fdelm){
- SetConstants(nf);
- return (this->mMS2mOS(mMS, mq, asmu, mu, nl, fdelm));
-}
-
-double CRunDec::mOS2mMS(double mOS, std::pair* mq, double asmu, double Mu,int nf,
- int nl, double fdelm){
- SetConstants(nf);
- return (this->mOS2mMS(mOS, mq, asmu, Mu, nl, fdelm));
-}
-
-double CRunDec::mMS2mSI(double mMS, double asmu, double mu,int nf, int nl){
- SetConstants(nf);
- return (this->mMS2mSI(mMS, asmu, mu, nl));
-}
-
-double CRunDec::mRI2mMS(double mRI, double asmu,int nf, int nl){
- SetConstants(nf);
- return (this->mRI2mMS(mRI, asmu, nl));
-}
-
-double CRunDec::mMS2mRGI(double mMS, double asmu,int nf, int nl){
- SetConstants(nf);
- return (this->mMS2mRGI(mMS, asmu, nl));
-}
-
-double CRunDec::mRGI2mMS(double mRGI, double asmu,int nf, int nl){
- SetConstants(nf);
- return (this->mRGI2mMS(mRGI, asmu, nl));
-}
-
-double CRunDec::mOS2mSI(double mOS, std::pair* mq, double asM,int nf, int nl, double fdelm){
- SetConstants(nf);
- return (this->mOS2mSI(mOS, mq, asM, nl, fdelm));
-}
-
-
-double CRunDec::mOS2mMSrun(double mOS, std::pair* mq, double asmu, double mu,
- int nf, int nl){
- SetConstants(nf);
- return (this->mOS2mMSrun(mOS, mq, asmu, mu, nl));
-}
-
-double CRunDec::mMS2mOSrun(double mMS, std::pair* mq, double asmu, double mu,
- int nf, int nl){
- SetConstants(nf);
- return (this->mMS2mOSrun(mMS, mq, asmu, mu, nl));
-}
-
-double CRunDec::mMS2mRI(double mMS, double asmu,int nf, int nl){
- SetConstants(nf);
- return (this->mMS2mRI(mMS, asmu, nl));
-}
-
-double CRunDec::mOS2mMSit(double mOS, std::pair* mq, double asmu, double mu,
- int nf,int nl){
- SetConstants(nf);
- return (this->mOS2mMSit(mOS, mq, asmu, mu, nl));
-}
-
-double CRunDec::mMS2mRGImod(double mMS, double asmu,int nf, int nl){
- SetConstants(nf);
- return (this->mMS2mRGImod(mMS, asmu, nl));
-}
-
-
-double CRunDec::mOS2mRS(double mOS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops) {
- return (this->mOS2mRS(mOS, mq, asmu, mu, nuf, nl, nloops, false));
-}
-
-double CRunDec::mOS2mRSp(double mOS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops) {
- return (this->mOS2mRS(mOS, mq, asmu, mu, nuf, nl, nloops, true));
-}
-
-double CRunDec::mMS2mRS(double mMS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops, double fdelm) {
- return (this->mMS2mRS(mMS, mq, asmu, mu, nuf, nl, nloops, fdelm, false));
-}
-
-double CRunDec::mMS2mRSp(double mMS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops, double fdelm) {
- return (this->mMS2mRS(mMS, mq, asmu, mu, nuf, nl, nloops, fdelm, true));
-}
-
-double CRunDec::mRS2mMS(double mRS, std::pair* mq, double asmu, double mu, double muf, int nl, int nloops, double fdelm) {
- return (this->mRS2mMS(mRS, mq, asmu, mu, muf, nl, nloops, fdelm, false));
-}
-
-double CRunDec::mRSp2mMS(double mRS, std::pair* mq, double asmu, double mu, double muf, int nl, int nloops, double fdelm) {
- return (this->mRS2mMS(mRS, mq, asmu, mu, muf, nl, nloops, fdelm, true));
-}
-
-double CRunDec::mRS2mSI(double mRS, std::pair* mq, double (*as)(double), double muf, int nl, int nloops, double fdelm) {
- return (this->mRS2mSI(mRS, mq, as, muf, nl, nloops, fdelm, false));
-}
-
-double CRunDec::mRSp2mSI(double mRS, std::pair* mq, double (*as)(double), double muf, int nl, int nloops, double fdelm) {
- return (this->mRS2mSI(mRS, mq, as, muf, nl, nloops, fdelm, true));
-}
-
-double CRunDec::DecAsDownOS(double als, double massth, double muth,int nf,
- int nl){
- SetConstants(nf);
- return (this->DecAsDownOS(als, massth, muth, nl));
-}
-
-double CRunDec::DecAsUpOS(double als, double massth, double muth, int nf,
- int nl){
- SetConstants(nf);
- return (this->DecAsUpOS(als, massth, muth, nl));
-}
-
-double CRunDec::DecAsUpMS(double als, double massth, double muth, int nf,
- int nl){
- SetConstants(nf);
- return (this->DecAsUpMS(als, massth, muth, nl));
-}
-
-double CRunDec::DecAsDownMS(double als, double massth, double muth,int nf,
- int nl){
- SetConstants(nf);
- return (this->DecAsDownMS(als, massth, muth, nl));
-}
-
-double CRunDec::DecAsUpSI(double als, double massth, double muth, int nf,
- int nl){
- SetConstants(nf);
- return (this->DecAsUpSI(als, massth, muth, nl));
-}
-
-double CRunDec::DecAsDownSI(double als, double massth, double muth, int nf,
- int nl){
- SetConstants(nf);
- return (this->DecAsDownSI(als, massth, muth, nl));
-}
-
-double CRunDec::DecMqUpOS(double mq, double als, double massth, double muth,
- int nf, int nl){
- SetConstants(nf);
- return (this->DecMqUpOS(mq, als, massth, muth, nl));
-}
-
-double CRunDec::DecMqDownOS(double mq, double als, double massth,
- double muth,int nf, int nl){
- SetConstants(nf);
- return (this->DecMqDownOS(mq, als, massth, muth, nl));
-}
-
-double CRunDec::DecMqUpMS(double mq, double als, double massth, double muth,
- int nf, int nl){
- SetConstants(nf);
- return (this->DecMqUpMS(mq, als, massth, muth, nl));
-}
-
-double CRunDec::DecMqDownMS(double mq, double als, double massth,
- double muth,int nf, int nl){
- SetConstants(nf);
- return (this->DecMqDownMS(mq, als, massth, muth, nl));
-}
-
-double CRunDec::DecMqUpSI(double mq, double als, double massth, double muth,
- int nf, int nl){
- SetConstants(nf);
- return (this->DecMqUpSI(mq, als, massth, muth, nl));
-}
-
-double CRunDec::DecMqDownSI(double mq, double als, double massth,
- double muth,int nf, int nl){
- SetConstants(nf);
- return (this->DecMqDownSI(mq, als, massth, muth, nl));
-}
diff --git a/CRunDec3/CRunDec.3.0.h b/CRunDec3/CRunDec.3.0.h
deleted file mode 100644
index 95269b5..0000000
--- a/CRunDec3/CRunDec.3.0.h
+++ /dev/null
@@ -1,1017 +0,0 @@
-/*
- CRunDec.h
-
- Header file for CRunDec.cpp
-
- Author: Barbara Schmidt (Jan 2012)
- Florian Herren and Matthias Steinhauser (Jan 2016)
-*/
-
-/*
- Minor update: (Sep 2016)
- Remove "static" from definition and initialization of constants for Runge-Kutta procedure.
- The C++11 standard supports the new version, so a C++11 compliant compiler and maybe
- the switch "-std=c++11" should be used.
-*/
-
-/*
-License:
-
-CRunDec is free software: you can redistribute it and/or modify
-it under the terms of the GNU General Public License as published by
-the Free Software Foundation, either version 3 of the License, or
-(at your option) any later version.
-
-This program is distributed in the hope that it will be useful,
-but WITHOUT ANY WARRANTY; without even the implied warranty of
-MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-GNU General Public License for more details.
-
-You should have received a copy of the GNU General Public License
-along with this program. If not, see .
-*/
-
-
-#ifndef CRUNDEC_H_
-#define CRUNDEC_H_
-
-#include
-
-// Default return statement:
-#define RETURN return(0);
-// The following might be useful for Windows:
-//#define RETURN system("PAUSE"); exit(1);
-
-
-// Numerical values for input parameters:
-static struct RunDec_values {
- double asMz = 0.1181;
- double asMtau = 0.332;
- double Mz = 91.1876;
- double Mh = 125.09;
- double muc = 1.279;
- double mc3 = 0.986;
- double mub = 4.163;
- double Mtau = 1.77686;
- double Mc = 1.5;
- double Mb = 4.8;
- double Mt = 173.21;
-} NumDef;
-
-
-
-
-// Struct for triple {nf, Mth, Muth}:
-
-//! Structure containing: the number of light flavours (nf), the mass of the heavy quark (Mth) and the decoupling scale (muth)
-/*!
- This structure is used to pass information on decoupling thresholds to AlH2AlL, AlL2AlH, mH2mL and mL2mH.
-*/
-struct TriplenfMmu{
- int nf;
- double Mth;
- double muth;
-};
-
-
-//! Structure containing: \f$\alpha_s^{(n_f)}\f$ (Asexact) and \f$m_{MS}^{(n_f)}\f$ (mMSexact)
-/*!
- This structure is used to return the results of AsmMsexact.
-*/
-struct AsmMS{
- double Asexact;
- double mMSexact;
-};
-
-// class declaration of CRunDec:
-
-//! Main class, contains all functions
-/*!
- All RunDec functions are accessed via this class.
- For many of the functions versions with and without the parameters nl/nf exist.
- The ones with the parameters nl/nf call the ones without after calling SetNf(nl)/SetNf(nf). It is
- recommended not to use the ones without nl/nf since some of them change the values of nl/nf during their execution.
-
- Detailed descriptions are only provided for the the functions
- with parameters nl/nf.
-*/
-class CRunDec
-{
-private:
- // Aux. constants for implicit Runge-Kutta-Procedure:
- double a2=0.2, a3=0.3, a4=0.6, a5=1., a6=0.875;
-
- double b21=0.2, b31=3./40., b32=9./40., b41=0.3, b42=-0.9,
- b43=6./5.;
- double b51=-11./54., b52=2.5, b53=-70./27., b54=35./27.;
- double b61=1631./55296., b62=175./512., b63=575./13824.;
- double b64=44275./110592., b65=253./4096.;
-
- double c1=37./378., c2=0., c3=250./621., c4=125./594., c5=0.;
- double c6= 512./1771.;
-
- double dc1=37./378.-2825./27648., dc2=0.-0.,
- dc3=250./621.-18575./48384.;
- double dc4=125./594.-13525./55296., dc5=0.-277./14336.,
- dc6=512./1771.-0.25;
-
- // Coefficients for diff. equations:
- double Beta[5], B[5], Betap[5], Bp[5], Gamma[5], C[5], Nf;
-
- // Define constants (if not already done with constructor):
- void SetConstants(int n);
-
- // R.h.s. of diff. equations:
- friend double fSetdydx(CRunDec S, double A, int nl);
-
- friend double fSetdydxa1(CRunDec S, double x, double A);
- friend double fSetdydxM1(CRunDec S, double A, double M);
-
- friend double fSetdydxa2(CRunDec S, double x, double A);
- friend double fSetdydxM2(CRunDec S, double A, double M);
-
- friend double fSetdydxa3(CRunDec S, double x, double A);
- friend double fSetdydxM3(CRunDec S, double A, double M);
-
- friend double fSetdydxa4(CRunDec S, double x, double A);
- friend double fSetdydxM4(CRunDec S, double A, double M);
-
- friend double fSetdydxa5(CRunDec S, double x, double A);
- friend double fSetdydxM5(CRunDec S, double A, double M);
-
- // Additional aux. functions:
- int Abbruch(void);
- double fSetAsL(double Lambda, double Mu, int nl, double AlphaS);
- double fSetcx(double x, int nl);
- double fOsFromMs1(double mu, double M);
- double fOsFromMs2(double mu, double M, double nl);
- double fOsFromMs3(double mu, double M, double nl);
- double fOsFromMs4(double mu, double M, double nl, double err);
- double fMsFromOs1(double mu, double M);
- double fMsFromOs2(double mu, double M, double nl);
- double fMsFromOs3(double mu, double M, double nl);
- double fMsFromOs4(double mu, double M, double nl, double err);
- double fZmM(double n);
- double fZmInvM(double n);
- double deltamOS2mMS(double mOS, std::pair* mq,
- double asmu, double mu, int nlq, int nloops);
- double deltamMS2mOS(double mMS, std::pair* mq,
- double asmu, double mu, int nlq, int nloops);
- double fMsFromRi1(void);
- double fMsFromRi2(void);
- double fMsFromRi3(void);
- double fMumFromOs1(void);
- double fMumFromOs2(void);
- double fMumFromOs3(void);
- double fMumFromOs4(double err);
- double fRiFromMs(double alpha, double nl);
- double fMsFromRi(double alpha, double nl);
- double fHelpmOS2mMSit(double mMS,double mOS, std::pair* mq,
- double asmu, double mu, int nl);
- double fas5to6os(double alpha, double mass, double mu, double nlq, double nl);
- double fas6to5os(double alpha, double mass, double mu, double nlq, double nl);
- double fas5to6ms(double alpha, double mass, double mu, double nlq, double nl);
- double fas6to5ms(double alpha, double mass, double mu, double nlq, double nl);
- double fas6to5si(double alpha, double mass, double mu, double nlq, double nl);
- double fas5to6si(double alpha, double mass, double mu, double nlq, double nl);
- double fmq5to6os(double A, double mass, double mu, double nlq, double nl);
- double fmq6to5os(double A, double mass, double mu, double nlq, double nl);
- double fmq5to6ms(double A, double mass, double mu, double nlq, double nl);
- double fmq6to5ms(double A, double mass, double mu, double nlq, double nl);
- double fmq5to6si(double A, double mass, double mu, double nlq, double nl);
- double fmq6to5si(double A, double mass, double mu, double nlq, double nl);
-
- double PSdelta(double asmu, double muf, double mu, int nl, int nloops);
- double E1p(double mOS, double asmu, double mu, int nl, int nloops);
- double exOS2RS(double api, double mmu, double nnuf, int nnl, int nloops);
- double exOS2RSp(double api, double mmu, double nnuf, int nnl, int nloops);
- double mMS2mOSmod(double mMS, std::pair* mq,
- double asmu, double mu, int nf, int nloops, double err);
-
- double fRungeKuttaImpl(double &x, double y, double &htry, int nl,
- double (*f)(CRunDec, double, int));
- double fRKSchritt(double x, double y, double h, double &yerr,
- double (*f)(CRunDec, double, double));
- double PolyLog(unsigned int n, double x);
-
-
-public:
- // constructor:
- CRunDec();
- CRunDec(int);
-
- // Arrays and structs to store data:
- std::pair mq[4];
- TriplenfMmu nfMmu[4];
- AsmMS AM;
-
-
- //! GetNf returns the number of light flavours currently in use
- /*!
- \return \f$n_f\f$
- */
- int GetNf();
-
- //! SetNf sets the number of light flavours
- /*!
- \param nf \f$n_f\f$
- */
- void SetNf(int nf);
-
- // Functions for the running of alpha_s and m_q:
-
- //! LamExpl calculates \f$\Lambda^{(n_f)}\f$ from \f$\alpha_s^{(n_f)}\f$ explicitly solving for \f$\Lambda^{(n_f)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$\Lambda^{(n_f)}\f$
- */
- double LamExpl(double asmu, double mu, int nf, int nloops);
-
- //! LamImpl calculates \f$\Lambda^{(n_f)}\f$ from \f$\alpha_s^{(n_f)}\f$ implicitly solving for \f$\Lambda^{(n_f)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$\Lambda^{(n_f)}\f$
- */
- double LamImpl(double asmu, double mu,int nf,int nloops);
-
- //! AlphasLam calculates \f$\alpha_s^{(n_f)}\f$ from \f$\Lambda^{(n_f)}\f$
- /*!
- \param Lambda \f$\Lambda^{(n_f)}\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_f)}(\mu)\f$
- */
- double AlphasLam(double Lambda, double mu,int nf, int nloops);
-
- //! AlphasExact calculates \f$\alpha_s^{(n_f)}(\mu_1)\f$ from \f$\alpha_s^{(n_f)}(\mu_0)\f$
- /*!
- \param asmu0 \f$\alpha_s^{(n_f)}(\mu_0)\f$
- \param mu0 \f$\mu_0\f$
- \param mu1 \f$\mu_1\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_f)}(\mu_1)\f$
- */
- double AlphasExact(double asmu0, double mu0, double mu1, int nf,int nloops);
-
- //! mMS2mMS calculates \f$m_{MS}^{(n_f}(\mu_1)\f$ from \f$m_{MS}^{(n_f)}(\mu_0)\f$
- /*!
- \param mu0 \f$m_{MS}^{(n_f)}(\mu_0)\f$
- \param asmu0 \f$\alpha_s^{(n_f)}(\mu_0)\f$
- \param asmu1 \f$\alpha_s^{(n_f)}(\mu_1)\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_f)}(\mu_1)\f$
- */
- double mMS2mMS(double mu0, double asmu1, double asmu0,int nf, int nloops);
-
- //! AsmMSrunexact solves simultaneously the differential equations for \f$\alpha_s^{(n_f)}\f$ and \f$m_{MS}^{(n_f)}\f$
-
- //! Returns a structure of type AsmMS.
- /*!
- \param mmu \f$m_{MS}^{(n_f)}(\mu_0)\f$
- \param asmu0 \f$\alpha_s^{(n_f)}(\mu_0)\f$
- \param mu0 \f$\mu_0\f$
- \param mu1 \f$\mu_1\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$(\alpha_s^{(n_f)}(\mu_1),m_{MS}^{(n_f)}(\mu_1))\f$
- */
- AsmMS AsmMSrunexact(double mmu, double asmu0, double mu0, double mu1,
- int nf, int nloops);
-
- // Decoupling relations:
- //! DecLambdaUp calculates \f$\Lambda^{(n_l+1)}\f$ from \f$\Lambda^{(n_l)}\f$
- /*!
- \param lam \f$\Lambda^{(n_l)}\f$
- \param massth quark mass at which the matching is performed (\f$m_{MS}^{(n_l+1)}(m_{MS})\f$)
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\Lambda^{(n_l+1)}\f$
- */
- double DecLambdaUp(double lam, double massth, int nl, int nloops);
-
- //! DecLambdaDown calculates \f$\Lambda^{(n_l)}\f$ from \f$\Lambda^{(n_l+1)}\f$
- /*!
- \param lam \f$\Lambda^{(n_l+1)}\f$
- \param massth quark mass at which the matching is performed (\f$m_{MS}^{(n_l+1)}(m_{MS})\f$)
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\Lambda^{(n_l)}\f$
- */
- double DecLambdaDown(double lam, double massth, int nl, int nloops);
-
-
- //! DecAsDownOS calculates \f$\alpha_s^{(n_l)}\f$ from \f$\alpha_s^{(n_l+1)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_l+1)}(\mu)\f$
- \param massth mass of the heavy quark (\f$M_{OS}\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_l)}(\mu)\f$
- */
- double DecAsDownOS(double asmu, double massth, double muth, int nl, int nloops);
-
- //! DecAsUpOS calculates \f$\alpha_s^{(n_l+1)}\f$ from \f$\alpha_s^{(n_l)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param massth mass of the heavy quark (\f$M_{OS}\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_l+1)}(\mu)\f$
- */
- double DecAsUpOS(double asmu, double massth, double muth, int nl, int nloops);
-
- //! DecAsDownMS calculates \f$\alpha_s^{(n_l)}\f$ from \f$\alpha_s^{(n_l+1)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_l+1)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(\mu)\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_l)}(\mu)\f$
- */
- double DecAsDownMS(double asmu, double massth, double muth, int nl, int nloops);
-
- //! DecAsUpMS calculates \f$\alpha_s^{(n_l+1)}\f$ from \f$\alpha_s^{(n_l)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(\mu)\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_l+1)}(\mu)\f$
- */
- double DecAsUpMS(double asmu, double massth, double muth, int nl, int nloops);
-
- //! DecAsDownSI calculates \f$\alpha_s^{(n_l)}\f$ from \f$\alpha_s^{(n_l+1)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_l+1)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(m_{MS})\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_l)}(\mu)\f$
- */
- double DecAsDownSI(double asmu, double massth, double muth, int nl, int nloops);
-
- //! DecAsUpSI calculates \f$\alpha_s^{(n_l+1)}\f$ from \f$\alpha_s^{(n_l)}\f$
- /*!
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(m_{MS})\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$\alpha_s^{(n_l+1)}(\mu)\f$
- */
- double DecAsUpSI(double asmu, double massth, double muth, int nl, int nloops);
-
- //! DecMqUpOS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m_{MS}^{(n_l)}(\mu)\f$
- /*!
- \param mq \f$m_{MS}^{(n_l)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param massth mass of the heavy quark (\f$M_{OS}\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double DecMqUpOS(double mq, double asmu, double massth, double muth, int nl,
- int nloops);
-
- //! DecMqDownOS calculates \f$m_{MS}^{(n_l)}(\mu)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
- /*!
- \param mq \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_l+1)}(\mu)\f$
- \param massth mass of the heavy quark (\f$M_{OS}\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_l)}(\mu)\f$
- */
- double DecMqDownOS(double mq, double asmu, double massth, double muth, int nl,
- int nloops);
-
- //! DecMqUpMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m_{MS}^{(n_l)}(\mu)\f$
- /*!
- \param mq \f$m_{MS}^{(n_l)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(\mu)\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double DecMqUpMS(double mq, double asmu, double massth, double muth, int nl,
- int nloops);
-
- //! DecMqDownMS calculates \f$m_{MS}^{(n_l)}(\mu)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
- /*!
- \param mq \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_l+1)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(\mu)\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_l)}(\mu)\f$
- */
- double DecMqDownMS(double mq, double asmu, double massth, double muth, int nl,
- int nloops);
-
- //! DecMqUpSI calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m_{MS}^{(n_l)}(\mu)\f$
- /*!
- \param mq \f$m_{MS}^{(n_l)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(m_{MS})\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double DecMqUpSI(double mq, double asmu, double massth, double muth, int nl,
- int nloops);
-
- //! DecMqDownSI calculates \f$m_{MS}^{(n_l)}(\mu)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
- /*!
- \param mq \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_l+1)}(\mu)\f$
- \param massth mass of the heavy quark (\f$m_{MS}^{(n_l+1)}(m_{MS})\f$)
- \param muth scale \f$\mu\f$ at which the matching is performed
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_l)}(\mu)\f$
- */
- double DecMqDownSI(double mq, double asmu, double massth, double muth, int nl,
- int nloops);
-
- // Running and decoupling:
-
- //! AlL2AlH calculates \f$\alpha_s(\mu_2)\f$ from \f$\alpha_s(\mu_1)\f$ decoupling at intermediate scales, running from low to high
- /*!
- \param asl \f$\alpha_s(\mu_1)\f$
- \param mu1 \f$\mu_1\f$
- \param decpar arrays of triples indicating the number of flavours, the OS-mass and the scale at which decoupling is performed
- \param mu2 \f$\mu_2\f$
- \param nloops number of loops
- \return \f$\alpha_s(\mu_2)\f$
- */
- double AlL2AlH(double asl, double mu1, TriplenfMmu decpar[], double mu2, int nloops);
-
- //! AlH2AlL calculates \f$\alpha_s(\mu_2)\f$ from \f$\alpha_s(\mu_1)\f$ decoupling at intermediate scales, running from high to low
- /*!
- \param ash \f$\alpha_s(\mu_1)\f$
- \param mu1 \f$\mu_1\f$
- \param decpar arrays of triples indicating the number of flavours, the OS-mass and the scale at which decoupling is performed
- \param mu2 \f$\mu_2\f$
- \param nloops number of loops
- \return \f$\alpha_s(\mu_2)\f$
- */
- double AlH2AlL(double ash, double mu1, TriplenfMmu decpar[], double mu2, int nloops);
-
- //! mL2mH calculates \f$m_{MS}(\mu_2)\f$ from \f$m_{MS}(\mu_1)\f$ decoupling at intermediate scales, running from low to high
- /*!
- \param mql \f$m_{MS}(\mu_1)\f$
- \param asl \f$\alpha_s(\mu_1)\f$
- \param mu1 \f$\mu_1\f$
- \param decpar arrays of triples indicating the number of flavours, the OS-mass and the scale at which decoupling is performed
- \param mu2 \f$\mu_2\f$
- \param nloops number of loops
- \return \f$m_{MS}(\mu_2)\f$
- */
- double mL2mH(double mql, double asl, double mu1, TriplenfMmu decpar[], double mu2,
- int nloops);
-
- //! mH2mL calculates \f$m_{MS}(\mu_2)\f$ from \f$m_{MS}(\mu_1)\f$ decoupling at intermediate scales, running from high to low
- /*!
- \param mqh \f$m_{MS}(\mu_1)\f$
- \param ash \f$\alpha_s(\mu_1)\f$
- \param mu1 \f$\mu_1\f$
- \param decpar arrays of triples indicating the number of flavours, the OS-mass and the scale at which decoupling is performed
- \param mu2 \f$\mu_2\f$
- \param nloops number of loops
- \return \f$m_{MS}(\mu_2)\f$
- */
- double mH2mL(double mqh, double ash, double mu1, TriplenfMmu decpar[], double mu2,
- int nloops);
-
- // Mass relations
-
- //! mMS2mOS calculates \f$M_{OS}\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
- /*!
- \param mMS \f$m_{MS}^{(n_f)}(\mu)\f$
- \param mq pointer to pairs of light quark masses
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$M_{OS}\f$
- */
-
- double mMS2mOS(double mMS, std::pair* mq,
- double asmu, double mu,int nf, int nloops, double fdelm=1.0);
-
- //! mOS2mMS calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$M_{OS}\f$
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_f)}(\mu)\f$
- */
- double mOS2mMS(double mOS, std::pair* mq,
- double asmu, double mu,int nf, int nloops, double fdelm=1.0);
-
- //! mMS2mSI calculates \f$m_{MS}^{(n_f)}(m_{MS})\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
- /*!
- \param mMS \f$m_{MS}^{(n_f)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_f)}(m_{MS})\f$
- */
- double mMS2mSI(double mMS, double asmu, double mu, int nf, int nloops);
-
- //! mRI2mMS calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$m^{RI}(\mu)\f$
- /*!
- \param mRI \f$m^{RI}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_f)}(\mu)\f$
- */
- double mRI2mMS(double mRI, double asmu, int nf, int nloops);
-
- //! mMS2mRGI calculates \f$m^{RGI}\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
-
- //! for the difference between mMS2mRGImod and mMS2mRGI see arXiv:1201:6149
- /*!
- \param mMS \f$m_{MS}^{(n_f)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m^{RGI}\f$
- */
- double mMS2mRGI(double mMS, double asmu, int nf, int nloops);
-
- //! mRGI2mMS calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$m^{RGI}\f$
- /*!
- \param mRGI \f$m^{RGI}\f$
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_f)}(\mu)\f$
- */
- double mRGI2mMS(double mRGI, double asmu, int nf, int nloops);
-
- //! mOS2mSI calculates \f$m_{MS}^{(n_f)}(m_{MS})\f$ from \f$M_{OS}\f$
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses
- \param asM \f$\alpha_s^{(n_f)}(M_{OS})\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_f)}(m_{MS})\f$
- */
- double mOS2mSI(double mOS, std::pair* mq,
- double asM, int nf, int nloops, double fdelm=1.0);
-
- //! mOS2mMSrun calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$M_{OS}\f$
-
- //! Calculates \f$m_{MS}^{(n_f)}(m_{MS})\f$ in an intermediate step
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_f)}(\mu)\f$
- */
- double mOS2mMSrun(double mOS, std::pair* mq,
- double asmu, double mu,int nf, int nloops);
-
- //! mMS2mOSrun calculates \f$M_{OS}\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
-
- //! Calculates \f$m_{MS}^{(n_f)}(m_{MS})\f$ in an intermediate step
- /*!
- \param mMS \f$m_{MS}^{(n_f)}(\mu)\f$
- \param mq pointer to pairs of light quark masses
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$M_{OS}\f$
- */
- double mMS2mOSrun(double mMS, std::pair* mq,
- double asmu, double mu,int nf, int nloops);
-
- //! mMS2mRI calculates \f$m^{RI}(\mu)\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
- /*!
- \param mMS \f$m_{MS}^{(n_f)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m^{RI}(\mu)\f$
- */
- double mMS2mRI(double mMS, double asmu,int nf, int nloops);
-
- //! mOS2mMSit calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$M_{OS}\f$
-
- //! DEPRECATED, will be removed in future versions
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m_{MS}^{(n_f)}(\mu)\f$
- */
- double mOS2mMSit(double mOS, std::pair* mq,
- double asmu, double mu, int nf, int nloops);
-
- //! mMS2mRGImod calculates \f$m^{RGI}\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
-
- //! for the difference between mMS2mRGImod and mMS2mRGI see arXiv:1201:6149
- /*!
- \param mMS \f$m_{MS}^{(n_f)}(\mu)\f$
- \param asmu \f$\alpha_s^{(n_f)}(\mu)\f$
- \param nf \f$n_f\f$
- \param nloops number of loops
- \return \f$m^{RGI}\f$
- */
- double mMS2mRGImod(double mMS, double asmu, int nf, int nloops);
-
- //! mOS2mPS calculates \f$m^{PS}(\mu_f)\f$ from \f$M_{OS}\f$
- /*!
- \param mOS \f$M_{OS}\f$
-
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param muf \f$\mu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m^{PS}(\mu_f)\f$
- */
- double mOS2mPS(double mOS, std::pair* mq,
- double asmu, double mu, double muf, int nl, int nloops);
-
- //! mMS2mPS calculates \f$m^{PS}(\mu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
- /*!
- \param mMS \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param muf \f$\mu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m^{PS}(\mu_f)\f$
- */
- double mMS2mPS(double mMS, std::pair* mq,
- double asmu, double mu, double muf, int nl, int nloops, double fdelm=1.0);
-
- //! mPS2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{PS}(\mu_f)\f$
- /*!
- \param mPS \f$m^{PS}(\mu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param muf \f$\mu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double mPS2mMS(double mPS, std::pair* mq,
- double asmu, double mu, double muf, int nl, int nloops, double fdelm=1.0);
-
- //! mPS2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{PS}(\mu_f)\f$
- /*!
- \param mPS \f$m^{PS}(\mu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param as pointer to a function which calculates \f$\alpha_s^{(n_l)}(\mu)\f$ and takes \f$\mu\f$ as argument
- \param muf \f$\mu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
- */
- double mPS2mSI(double mPS, std::pair* mq,
- double (*as)(double), double muf, int nl, int nloops, double fdelm=1.0);
-
- //! mOS2m1S calculates \f$m^{1S}\f$ from \f$M_{OS}\f$
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m^{1S}\f$
- */
- double mOS2m1S(double mOS, std::pair* mq,
- double asmu, double mu, int nl, int nloops);
-
- //! mMS2m1S calculates \f$m^{1S}\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
- /*!
- \param mMS \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m^{1S}\f$
- */
- double mMS2m1S(double mMS, std::pair* mq,
- double asmu, double mu, int nl, int nloops, double fdelm=1.0);
-
- //! m1S2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{1S}\f$
- /*!
- \param m1S \f$m^{PS}\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double m1S2mMS(double m1S, std::pair* mq,
- double asmu, double mu, int nl, int nloops, double fdelm=1.0);
-
- //! m1S2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{1S}\f$
- /*!
- \param m1S \f$m^{1S}\f$
- \param mq pointer to pairs of light quark masses, not active
- \param as pointer to a function which calculates \f$\alpha_s^{(n_l)}(\mu)\f$ and takes \f$\mu\f$ as argument
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
- */
- double m1S2mSI(double m1S, std::pair* mq,
- double (*as)(double), int nl, int nloops, double fdelm=1.0);
-
- //! mOS2mRS calculates \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$ from \f$M_{OS}\f$
-
- //! This is a helper function, use the version without the parameter prime or mOS2mRSp
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ should be calculated
- \return \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
- */
- double mOS2mRS(double mOS, std::pair* mq, double asmu,
- double mu, double nuf, int nl, int nloops, bool prime);
-
- //! mMS2mRS calculates \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
-
- //! This is a helper function, use the version without the parameter prime or mMS2mRSp
- /*!
- \param mMS \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ should be calculated
- \return \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
- */
- double mMS2mRS(double mMS, std::pair* mq, double asmu,
- double mu, double nuf, int nl, int nloops, double fdelm, bool prime);
-
- //! mRS2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
-
- //! This is a helper function, use the version without the parameter prime or mRSp2mMS
- /*!
- \param mRS \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ is given
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double mRS2mMS(double mRS, std::pair* mq, double asmu,
- double mu, double nuf, int nl, int nloops, double fdelm, bool prime);
-
- //! mRS2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
-
- //! This is a helper function, use the version without the parameter prime or mRSp2mSI
- /*!
- \param mRS \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param as pointer to a function which calculates \f$\alpha_s^{(n_l)}(\mu)\f$ and takes \f$\mu\f$ as argument
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ is given
- \return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
- */
- double mRS2mSI(double mRS, std::pair* mq, double (*as)(double),
- double nuf, int nl, int nloops, double fdelm, bool prime);
-
-
- //! mOS2mRS calculates \f$m^{RS}(\nu_f)\f$ from \f$M_{OS}\f$
-
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m^{RS}(\nu_f)\f$
- */
- double mOS2mRS(double mOS, std::pair* mq,
- double asmu, double mu, double nuf, int nl, int nloops);
-
- //! mMS2mRS calculates \f$m^{RS}(\nu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
-
- /*!
- \param mMS \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m^{RS}(\nu_f)\f$
- */
- double mMS2mRS(double mMS, std::pair* mq,
- double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
-
- //! mRS2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{RS}(\nu_f)\f$
-
- /*!
- \param mRS \f$m^{RS}(\nu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double mRS2mMS(double mRS, std::pair* mq,
- double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
-
- //! mRS2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{RS}(\nu_f)\f$
-
- /*!
- \param mRS \f$m^{RS}(\nu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param as pointer to a function which calculates \f$\alpha_s^{(n_l)}(\mu)\f$ and takes \f$\mu\f$ as argument
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
- */
- double mRS2mSI(double mRS, std::pair* mq,
- double (*as)(double), double nuf, int nl, int nloops, double fdelm=1.0);
-
- //! mOS2mRSp calculates \f$m^{RS'}(\nu_f)\f$ from \f$M_{OS}\f$
-
- /*!
- \param mOS \f$M_{OS}\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \return \f$m^{RS'}(\nu_f)\f$
- */
- double mOS2mRSp(double mOS, std::pair* mq,
- double asmu, double mu, double nuf, int nl, int nloops);
-
- //! mMS2mRSp calculates \f$m^{RS'}(\nu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
-
- /*!
- \param mMS \f$m_{MS}^{(n_l+1)}(\mu)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m^{RS'}(\nu_f)\f$
- */
- double mMS2mRSp(double mMS, std::pair* mq,
- double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
-
- //! mRSp2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{RS'}(\nu_f)\f$
-
- /*!
- \param mRS \f$m^{RS'}(\nu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param asmu \f$\alpha_s^{(n_l)}(\mu)\f$
- \param mu \f$\mu\f$
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(\mu)\f$
- */
- double mRSp2mMS(double mRS, std::pair* mq,
- double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
-
- //! mRSp2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{RS'}(\nu_f)\f$
-
- /*!
- \param mRS \f$m^{RS'}(\nu_f)\f$
- \param mq pointer to pairs of light quark masses, not active
- \param as pointer to a function which calculates \f$\alpha_s^{(n_l)}(\mu)\f$ and takes \f$\mu\f$ as argument
- \param nuf \f$\nu_f\f$
- \param nl \f$n_l\f$
- \param nloops number of loops
- \param fdelm factor multiplying the non-logarithmic part of the 4-loop term
- \return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
- */
- double mRSp2mSI(double mRS, std::pair* mq,
- double (*as)(double), double nuf, int nl, int nloops, double fdelm=1.0);
-
-
- // Overload functions:
- double LamExpl(double asmu, double mu, int nloops);
- double LamImpl(double asmu, double mu,int nloops);
- double AlphasLam(double Lambda, double mu, int nloops);
- double AlphasExact(double asmu0, double mu0, double mu1, int nloops);
- double mMS2mMS(double mmu0, double asmu0, double asmu1, int nloops);
- AsmMS AsmMSrunexact(double mmu, double asmu0, double mu0, double mu1,
- int nloops);
-
- double DecAsDownOS(double asmu, double massth, double muth, int nloops);
- double DecAsUpOS(double asmu, double massth, double muth, int nloops);
- double DecAsDownMS(double asmu, double massth, double muth, int nloops);
- double DecAsUpMS(double asmu, double massth, double muth, int nloops);
- double DecAsDownSI(double asmu, double massth, double muth, int nloops);
- double DecAsUpSI(double asmu, double massth, double muth, int nloops);
-
- double DecMqUpOS(double mq, double asmu, double massth, double muth, int nloops);
- double DecMqDownOS(double mq, double asmu, double massth, double muth, int nloops);
- double DecMqUpMS(double mq, double asmu, double massth, double muth, int nloops);
- double DecMqDownMS(double mq, double asmu, double massth, double muth, int nloops);
- double DecMqUpSI(double mq, double asmu, double massth, double muth, int nloops);
- double DecMqDownSI(double mq, double asmu, double massth, double muth, int nloops);
-
-
- double mMS2mOS(double mMS, std::pair* mq,
- double asmu, double mu, int nloops, double fdelm=1.0);
- double mOS2mMS(double mOS, std::pair* mq,
- double asmu, double mu, int nloops, double fdelm=1.0);
-
- double mMS2mSI(double mMS, double asmu, double mu, int nloops);
- double mRI2mMS(double mRI, double asmu, int nloops);
- double mMS2mRGI(double mMS, double asmu, int nloops);
- double mRGI2mMS(double mRGI, double asmu, int nloops);
- double mOS2mSI(double mOS, std::pair* mq,
- double asM, int nloops, double fdelm=1.0);
- double mOS2mMSrun(double mOS, std::pair* mq,
- double asmu, double mu, int nloops);
- double mMS2mOSrun(double mMS, std::pair* mq,
- double asmu, double mu, int nloops);
- double mMS2mRI(double mMS, double asmu, int nloops);
- double mOS2mMSit(double mOS, std::pair* mq,
- double asmu, double mu, int nloops);
- double mMS2mRGImod(double mMS, double asmu, int nloops);
-
-
-
-};
-
-#endif
diff --git a/CRunDec3/CRunDec.3.1.cpp b/CRunDec3/CRunDec.cpp
similarity index 97%
rename from CRunDec3/CRunDec.3.1.cpp
rename to CRunDec3/CRunDec.cpp
index 6e98b17..44e9cf3 100644
--- a/CRunDec3/CRunDec.3.1.cpp
+++ b/CRunDec3/CRunDec.cpp
@@ -1,10 +1,19 @@
+/*
+ CRunDec.cpp, v3.1.1
+
+ Authors: Florian Herren and Matthias Steinhauser (May 2021)
+
+ Changes since v3.1:
+ September 2025: Bug fix in mMS2mOS and mOS2mMS: Missing nl^3*(Log[mu^2/m^2])^4 terms restored.
+ Thanks to Long Chen for pointing out that the terms are missing.
+*/
+
/*
CRunDec.cpp, v3.1
Authors: Florian Herren and Matthias Steinhauser (May 2021)
Changes since v3.0:
- Jul 2023: Scheme for the kinetic mass cahnged to integers (A=0, B=1, C=3, D=3)
May 2021: fix array lengths in AlphasLam and fSetAsL. Thanks to Florian Bernlochner for pointing out this bug.
Jul 2020: implemented mMS2mKIN and mKIN2mMS (arXiv:2005:06478)
Coefficients ctil[3][{3,4,5}] updated. Thanks to Christopher Lepenik for providing the new values.
@@ -123,7 +132,7 @@ SOFTWARE.
#include
#include
-#include "CRunDec.3.1.h"
+#include "CRunDec.h"
// Some constants:
#define cf 4./3.
@@ -132,7 +141,7 @@ SOFTWARE.
#define B4 -1.762800087073770864061897634679818807215137274389016762629478603776
#define A4 0.5174790616738993863307581618988629456223774751413792582443193479770
#define A5 0.5084005792422687074591088492585899413195411256648216487244977963526
-#define Pi 3.141592653589793238462643
+#define Pi M_PI
#define Zeta2 (Pi*Pi)/6.
#define Zeta3 1.20205690315959428539973816151144999076498629234049888179227155534
#define Zeta4 (Pi*Pi*Pi*Pi)/90.
@@ -833,7 +842,7 @@ double CRunDec::fMsFromOs4(double mu, double M, double nl, double err){
erg = - 3654.15040757339*err - 1524.2292266911543*lmM - 288.778291935394*lmM*lmM - 32.54735725308642*lmM*lmM*lmM -
1.85546875*lmM*lmM*lmM*lmM +
nl*nl*nl*(0. + 0.678141025604516*err + 0.3205521521864135*lmM + 0.0800290327210927*lmM*lmM +
- 0.010030864197530864*lmM*lmM*lmM) +
+ 0.010030864197530864*lmM*lmM*lmM + 0.001157407407407407*lmM*lmM*lmM*lmM) +
nl*nl*(0. - 43.48241924867489*err - 19.82672048099557*lmM - 4.482957520194182*lmM*lmM -
0.5270061728395061*lmM*lmM*lmM - 0.04108796296296297*lmM*lmM*lmM*lmM) +
nl*(0. + 756.9421565599532*err + 330.1770776731065*lmM + 67.99849534415492*lmM*lmM +
@@ -841,9 +850,9 @@ double CRunDec::fMsFromOs4(double mu, double M, double nl, double err){
return erg;
}
-// Function: double CRunDec::mOS2mMS(double mOS, RunDecPair mq[], double asmu,
+// Function: double CRunDec::mOS2mMS(double mOS, std::pair* mq, double asmu,
// double Mu, int nl, double fdelm)
-double CRunDec::mOS2mMS(double mOS, RunDecPair mq[], double asmu, double Mu,int nl, double fdelm){
+double CRunDec::mOS2mMS(double mOS, std::pair* mq, double asmu, double Mu,int nl, double fdelm){
if(nl<0||nl>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<* mq, double asmu, double mu, int nlq, int nloops){
double erg=0.0;
if(!mq)
return 0.0;
@@ -1057,7 +1066,7 @@ double CRunDec::deltamOS2mMS(double mOS, RunDecPair mq[], double asmu, double mu
}
// Compute 2 Loop contribution to mMS2mOS at 3 loops for MS light quarks
-double CRunDec::deltamMS2mOS(double mMS, RunDecPair mq[], double asmu, double mu, int nlq, int nloops){
+double CRunDec::deltamMS2mOS(double mMS, std::pair* mq, double asmu, double mu, int nlq, int nloops){
double erg=0.0;
double lmu = log(mu*mu/(mMS*mMS));
if(!mq)
@@ -1112,9 +1121,9 @@ double CRunDec::fZmInvM(double nl){
return erg;
}
-// Function: double CRunDec::mMS2mOS(double mMS, RunDecPair mq[], double asmu,
+// Function: double CRunDec::mMS2mOS(double mMS, std::pair* mq, double asmu,
// double mu,int nl)
-double CRunDec::mMS2mOS(double mMS, RunDecPair mq[], double asmu, double mu, int nl, double fdelm){
+double CRunDec::mMS2mOS(double mMS, std::pair* mq, double asmu, double mu, int nl, double fdelm){
if(nl<0||nl>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<* mq, double asmu,
// double mu, int nf, int nloop, double fdelm), needed for several other mass relations
-double CRunDec::mMS2mOSmod(double mMS, RunDecPair mq[], double asmu, double mu,int nf, int nloop, double fdelm){
+double CRunDec::mMS2mOSmod(double mMS, std::pair* mq, double asmu, double mu,int nf, int nloop, double fdelm){
if(nloop<0||nloop>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloop <<" LOOPS"<* mq, double asM,
// int nl)
-double CRunDec::mOS2mSI(double mOS, RunDecPair mq[], double asM, int nl, double fdelm){
+double CRunDec::mOS2mSI(double mOS, std::pair* mq, double asM, int nl, double fdelm){
if(nl<0||nl>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<* mq, double asmu,
// double mu, int nl)
-double CRunDec::mOS2mMSrun(double mOS, RunDecPair mq[], double asmu, double mu,
+double CRunDec::mOS2mMSrun(double mOS, std::pair* mq, double asmu, double mu,
int nl){
double asM=0.0;
asM= this-> AlphasExact(asmu, mu, mOS, nl);
@@ -1276,9 +1285,9 @@ double CRunDec::mOS2mMSrun(double mOS, RunDecPair mq[], double asmu, double mu,
return newM;
}
-// Function: double CRunDec::mMS2mOSrun(double mMS, RunDecPair mq[], double asmu,
+// Function: double CRunDec::mMS2mOSrun(double mMS, std::pair* mq, double asmu,
// double mu, int nl)
-double CRunDec::mMS2mOSrun(double mMS, RunDecPair mq[], double asmu, double mu,
+double CRunDec::mMS2mOSrun(double mMS, std::pair* mq, double asmu, double mu,
int nl){
double mNeu = mMS2mSI(mMS, asmu, mu, nl);
double asmNeu = AlphasExact(asmu, mu, mNeu, nl);
@@ -1321,7 +1330,7 @@ double CRunDec::mMS2mRI(double mMS, double asmu, int nl){
}
// Coefficients needed for the transformation of mOS to mMSit
-double CRunDec::fHelpmOS2mMSit(double mMS, double mOS, RunDecPair mq[], double asmu,
+double CRunDec::fHelpmOS2mMSit(double mMS, double mOS, std::pair* mq, double asmu,
double mu, int nl){
double sum[4];
double deltalight = deltamMS2mOS(mMS,mq,asmu/Pi,mu,Nf-1,nl);
@@ -1343,9 +1352,9 @@ double CRunDec::fHelpmOS2mMSit(double mMS, double mOS, RunDecPair mq[], double a
return erg;
}
-// Function: double CRunDec::mOS2mMSit(double mOS, RunDecPair mq[], double asmu,
+// Function: double CRunDec::mOS2mMSit(double mOS, std::pair* mq, double asmu,
// double mu, int nl)
-double CRunDec::mOS2mMSit(double mOS, RunDecPair mq[], double asmu, double mu,
+double CRunDec::mOS2mMSit(double mOS, std::pair* mq, double asmu, double mu,
int nl){
if(nl<0||nl>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nl<<" LOOPS"<* mq, double asmu, double mu, double muf, int nl, int nloops)
+double CRunDec::mOS2mPS(double mOS, std::pair* mq, double asmu, double mu, double muf, int nl, int nloops) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
// double mu, double muf, int nl, int nloops, double fdelm)
-double CRunDec::mMS2mPS(double mMS, RunDecPair mq[], double asmu,
+double CRunDec::mMS2mPS(double mMS, std::pair* mq, double asmu,
double mu, double muf, int nl, int nloops, double fdelm) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
// double mu, double muf, int nl, int nloops, double fdelm)
-double CRunDec::mPS2mMS(double mPS, RunDecPair mq[], double asmu,
+double CRunDec::mPS2mMS(double mPS, std::pair* mq, double asmu,
double mu, double muf, int nl, int nloops, double fdelm) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
// double (*as)(double), double muf, int nl, int nloops, double fdelm)
// The function pointer passed should contain the adress of a function computing alpha_s in dependence of mu
-double CRunDec::mPS2mSI(double mPS, RunDecPair mq[],
+double CRunDec::mPS2mSI(double mPS, std::pair* mq,
double (*as)(double), double muf, int nl, int nloops, double fdelm) {
if(as == NULL) {
cout << "Pointer to as == NULL! Aborting..." << endl;
@@ -1544,9 +1553,9 @@ double CRunDec::E1p(double mOS, double asmu, double mu, int nl, int nloops) {
return -E*(4.*asmu*asmu*mOS)/9.;
}
-// Function: double CRunDec::mOS2m1S(double mOS, RunDecPair mq[],
+// Function: double CRunDec::mOS2m1S(double mOS, std::pair* mq,
// double asmu, double mu, int nl, int nloops)
-double CRunDec::mOS2m1S(double mOS, RunDecPair mq[], double asmu, double mu, int nl, int nloops) {
+double CRunDec::mOS2m1S(double mOS, std::pair* mq, double asmu, double mu, int nl, int nloops) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
// double asmu, double mu, int nl, int nloops, double fdelm)
-double CRunDec::mMS2m1S(double mMS, RunDecPair mq[],
+double CRunDec::mMS2m1S(double mMS, std::pair* mq,
double asmu, double mu, int nl, int nloops, double fdelm) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
// double asmu, double mu, int nl, int nloops, double fdelm)
-double CRunDec::m1S2mMS(double m1S, RunDecPair mq[],
+double CRunDec::m1S2mMS(double m1S, std::pair* mq,
double asmu, double mu, int nl, int nloops, double fdelm) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq,
// double (*as)(double), int nl, int nloops, double fdelm)
// The function pointer passed should contain the adress of a function computing alpha_s in dependence of mu
-double CRunDec::m1S2mSI(double m1S, RunDecPair mq[],
+double CRunDec::m1S2mSI(double m1S, std::pair* mq,
double (*as)(double), int nl, int nloops, double fdelm) {
if(as == NULL) {
cout << "Pointer to as == NULL! Aborting..." << endl;
@@ -1911,9 +1920,9 @@ double CRunDec::exOS2RSp(double api, double mmu, double nnuf, int nnl, int nloop
return res;
}
-// Function: CRunDec::mOS2mRS(double mOS, RunDecPair mq[], double asmu,
+// Function: CRunDec::mOS2mRS(double mOS, std::pair* mq, double asmu,
// double mu, double nuf, int nl, int nloops, bool prime)
-double CRunDec::mOS2mRS(double mOS, RunDecPair mq[], double asmu,
+double CRunDec::mOS2mRS(double mOS, std::pair* mq, double asmu,
double mu, double nuf, int nl, int nloops, bool prime) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
// double mu, double nuf, int nl, int nloops, double fdelm, bool prime)
-double CRunDec::mMS2mRS(double mMS, RunDecPair mq[], double asmu,
+double CRunDec::mMS2mRS(double mMS, std::pair* mq, double asmu,
double mu, double nuf, int nl, int nloops, double fdelm, bool prime) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double asmu,
// double mu, double muf, int nl, int nloops, double fdelm, bool prime)
-double CRunDec::mRS2mMS(double mRS, RunDecPair mq[], double asmu,
+double CRunDec::mRS2mMS(double mRS, std::pair* mq, double asmu,
double mu, double muf, int nl, int nloops, double fdelm, bool prime) {
if(nloops<0||nloops>4){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops<<" LOOPS"<* mq, double (*as)(double),
// double muf, int nl, int nloops, double fdelm, bool prime)
// The function pointer passed should contain the adress of a function computing alpha_s in dependence of mu
-double CRunDec::mRS2mSI(double mRS, RunDecPair mq[], double (*as)(double),
+double CRunDec::mRS2mSI(double mRS, std::pair* mq, double (*as)(double),
double muf, int nl, int nloops, double fdelm, bool prime) {
if(as == NULL) {
cout << "Pointer to as == NULL! Aborting..." << endl;
@@ -3162,22 +3171,27 @@ double CRunDec::mkin2mMSD(double mkin, double apinlmus, double mus, double mufac
return mMS;
}
-// Function: double CRunDec::mMS2mKIN(double mMS, RunDecPair mq[],
-// double asmus, double mus, double muf, int nl, int nloops, int deccase)
-double CRunDec::mMS2mKIN(double mMS, RunDecPair mq[],
- double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, int deccase) {
+// Function: double CRunDec::mMS2mKIN(double mMS, std::pair* mq,
+// double asmus, double mus, double muf, int nl, int nloops, std::string deccase)
+double CRunDec::mMS2mKIN(double mMS, std::pair* mq,
+ double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, std::string deccase) {
if(nloops<0||nloops>3){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops <<" LOOPS"< 3){
+ int decc = 0;
+ if(deccase == "A") decc = 0;
+ else if(deccase == "B") decc = 1;
+ else if(deccase == "C") decc = 2;
+ else if(deccase == "D") decc = 3;
+ else if(deccase == "") decc = 4;
+ else {
cout<<"DECCASE "<< deccase <<" NOT IMPLEMENTED"<* mq, double asmus, double mus, double muf, int nloops, std::string deccase) {
+ if(deccase == "A") return mMS2mKIN(mMS, mq, asmus, mus, muf, 3, 3, nloops, deccase);
+ else if(deccase == "B") return mMS2mKIN(mMS, mq, asmus, mus, muf, 3, 3, nloops, deccase);
+ else if(deccase == "C") return mMS2mKIN(mMS, mq, asmus, mus, muf, 3, 4, nloops, deccase);
+ else if(deccase == "D") return mMS2mKIN(mMS, mq, asmus, mus, muf, 3, 3, nloops, deccase);
else {
cout<<"DECCASE "<< deccase <<" NOT IMPLEMENTED"<* mq,
// double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, std::string deccase)
-double CRunDec::mKIN2mMS(double mKIN, RunDecPair mq[],
- double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, int deccase) {
+double CRunDec::mKIN2mMS(double mKIN, std::pair* mq,
+ double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, std::string deccase) {
if(nloops<0||nloops>3){
cout<<"PROCEDURE IS NOT IMPLEMENTED FOR "<< nloops <<" LOOPS"< 3) {
+ int decc = 0;
+ if(deccase == "A") decc = 0;
+ else if(deccase == "B") decc = 1;
+ else if(deccase == "C") decc = 2;
+ else if(deccase == "D") decc = 3;
+ else if(deccase == "") decc = 4;
+ else {
cout<<"DECCASE "<< deccase <<" NOT IMPLEMENTED"<* mq, double asmus, double mus, double muf, int nloops, std::string deccase) {
+ if(deccase == "A") return mKIN2mMS(mKIN, mq, asmus, mus, muf, 3, 3, nloops, deccase);
+ else if(deccase == "B") return mKIN2mMS(mKIN, mq, asmus, mus, muf, 3, 3, nloops, deccase);
+ else if(deccase == "C") return mKIN2mMS(mKIN, mq, asmus, mus, muf, 3, 4, nloops, deccase);
+ else if(deccase == "D") return mKIN2mMS(mKIN, mq, asmus, mus, muf, 3, 3, nloops, deccase);
else {
cout<<"DECCASE "<< deccase <<" NOT IMPLEMENTED"<AsmMSrunexact(mMu, AlphaS0, Mu0, MuEnd, nl));
}
-double CRunDec::mMS2mOS(double mMS, RunDecPair mq[], double asmu, double mu,int nf,
+double CRunDec::mMS2mOS(double mMS, std::pair* mq, double asmu, double mu,int nf,
int nl, double fdelm){
SetConstants(nf);
return (this->mMS2mOS(mMS, mq, asmu, mu, nl, fdelm));
}
-double CRunDec::mOS2mMS(double mOS, RunDecPair mq[], double asmu, double Mu,int nf,
+double CRunDec::mOS2mMS(double mOS, std::pair* mq, double asmu, double Mu,int nf,
int nl, double fdelm){
SetConstants(nf);
return (this->mOS2mMS(mOS, mq, asmu, Mu, nl, fdelm));
@@ -4206,19 +4226,19 @@ double CRunDec::mRGI2mMS(double mRGI, double asmu,int nf, int nl){
return (this->mRGI2mMS(mRGI, asmu, nl));
}
-double CRunDec::mOS2mSI(double mOS, RunDecPair mq[], double asM,int nf, int nl, double fdelm){
+double CRunDec::mOS2mSI(double mOS, std::pair* mq, double asM,int nf, int nl, double fdelm){
SetConstants(nf);
return (this->mOS2mSI(mOS, mq, asM, nl, fdelm));
}
-double CRunDec::mOS2mMSrun(double mOS, RunDecPair mq[], double asmu, double mu,
+double CRunDec::mOS2mMSrun(double mOS, std::pair* mq, double asmu, double mu,
int nf, int nl){
SetConstants(nf);
return (this->mOS2mMSrun(mOS, mq, asmu, mu, nl));
}
-double CRunDec::mMS2mOSrun(double mMS, RunDecPair mq[], double asmu, double mu,
+double CRunDec::mMS2mOSrun(double mMS, std::pair* mq, double asmu, double mu,
int nf, int nl){
SetConstants(nf);
return (this->mMS2mOSrun(mMS, mq, asmu, mu, nl));
@@ -4229,7 +4249,7 @@ double CRunDec::mMS2mRI(double mMS, double asmu,int nf, int nl){
return (this->mMS2mRI(mMS, asmu, nl));
}
-double CRunDec::mOS2mMSit(double mOS, RunDecPair mq[], double asmu, double mu,
+double CRunDec::mOS2mMSit(double mOS, std::pair* mq, double asmu, double mu,
int nf,int nl){
SetConstants(nf);
return (this->mOS2mMSit(mOS, mq, asmu, mu, nl));
@@ -4241,35 +4261,35 @@ double CRunDec::mMS2mRGImod(double mMS, double asmu,int nf, int nl){
}
-double CRunDec::mOS2mRS(double mOS, RunDecPair mq[], double asmu, double mu, double nuf, int nl, int nloops) {
+double CRunDec::mOS2mRS(double mOS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops) {
return (this->mOS2mRS(mOS, mq, asmu, mu, nuf, nl, nloops, false));
}
-double CRunDec::mOS2mRSp(double mOS, RunDecPair mq[], double asmu, double mu, double nuf, int nl, int nloops) {
+double CRunDec::mOS2mRSp(double mOS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops) {
return (this->mOS2mRS(mOS, mq, asmu, mu, nuf, nl, nloops, true));
}
-double CRunDec::mMS2mRS(double mMS, RunDecPair mq[], double asmu, double mu, double nuf, int nl, int nloops, double fdelm) {
+double CRunDec::mMS2mRS(double mMS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops, double fdelm) {
return (this->mMS2mRS(mMS, mq, asmu, mu, nuf, nl, nloops, fdelm, false));
}
-double CRunDec::mMS2mRSp(double mMS, RunDecPair mq[], double asmu, double mu, double nuf, int nl, int nloops, double fdelm) {
+double CRunDec::mMS2mRSp(double mMS, std::pair* mq, double asmu, double mu, double nuf, int nl, int nloops, double fdelm) {
return (this->mMS2mRS(mMS, mq, asmu, mu, nuf, nl, nloops, fdelm, true));
}
-double CRunDec::mRS2mMS(double mRS, RunDecPair mq[], double asmu, double mu, double muf, int nl, int nloops, double fdelm) {
+double CRunDec::mRS2mMS(double mRS, std::pair* mq, double asmu, double mu, double muf, int nl, int nloops, double fdelm) {
return (this->mRS2mMS(mRS, mq, asmu, mu, muf, nl, nloops, fdelm, false));
}
-double CRunDec::mRSp2mMS(double mRS, RunDecPair mq[], double asmu, double mu, double muf, int nl, int nloops, double fdelm) {
+double CRunDec::mRSp2mMS(double mRS, std::pair* mq, double asmu, double mu, double muf, int nl, int nloops, double fdelm) {
return (this->mRS2mMS(mRS, mq, asmu, mu, muf, nl, nloops, fdelm, true));
}
-double CRunDec::mRS2mSI(double mRS, RunDecPair mq[], double (*as)(double), double muf, int nl, int nloops, double fdelm) {
+double CRunDec::mRS2mSI(double mRS, std::pair* mq, double (*as)(double), double muf, int nl, int nloops, double fdelm) {
return (this->mRS2mSI(mRS, mq, as, muf, nl, nloops, fdelm, false));
}
-double CRunDec::mRSp2mSI(double mRS, RunDecPair mq[], double (*as)(double), double muf, int nl, int nloops, double fdelm) {
+double CRunDec::mRSp2mSI(double mRS, std::pair* mq, double (*as)(double), double muf, int nl, int nloops, double fdelm) {
return (this->mRS2mSI(mRS, mq, as, muf, nl, nloops, fdelm, true));
}
diff --git a/CRunDec3/CRunDec.3.1.h b/CRunDec3/CRunDec.h
similarity index 92%
rename from CRunDec3/CRunDec.3.1.h
rename to CRunDec3/CRunDec.h
index eb07c7f..7672d73 100644
--- a/CRunDec3/CRunDec.3.1.h
+++ b/CRunDec3/CRunDec.h
@@ -80,18 +80,6 @@ struct TriplenfMmu{
double muth;
};
-//RunDecPair mq[]
-//! Structure containing: the number the mass of the heavy quark (first) and its scale (second)
-/*!
- This structure is used to pass information on decoupling thresholds.
- This substitutes the original argument pair* mq and facilitate the interface to python.
-*/
-struct RunDecPair{
- double first;
- double second;
-};
-
-
//! Structure containing: \f$\alpha_s^{(n_f)}\f$ (Asexact) and \f$m_{MS}^{(n_f)}\f$ (mMSexact)
/*!
@@ -118,20 +106,20 @@ class CRunDec
{
private:
// Aux. constants for implicit Runge-Kutta-Procedure:
- double a2=0.2, a3=0.3, a4=0.6, a5=1., a6=0.875;
+ const double a2=0.2, a3=0.3, a4=0.6, a5=1., a6=0.875;
- double b21=0.2, b31=3./40., b32=9./40., b41=0.3, b42=-0.9,
+ const double b21=0.2, b31=3./40., b32=9./40., b41=0.3, b42=-0.9,
b43=6./5.;
- double b51=-11./54., b52=2.5, b53=-70./27., b54=35./27.;
- double b61=1631./55296., b62=175./512., b63=575./13824.;
- double b64=44275./110592., b65=253./4096.;
+ const double b51=-11./54., b52=2.5, b53=-70./27., b54=35./27.;
+ const double b61=1631./55296., b62=175./512., b63=575./13824.;
+ const double b64=44275./110592., b65=253./4096.;
- double c1=37./378., c2=0., c3=250./621., c4=125./594., c5=0.;
- double c6= 512./1771.;
+ const double c1=37./378., c2=0., c3=250./621., c4=125./594., c5=0.;
+ const double c6= 512./1771.;
- double dc1=37./378.-2825./27648., dc2=0.-0.,
+ const double dc1=37./378.-2825./27648., dc2=0.-0.,
dc3=250./621.-18575./48384.;
- double dc4=125./594.-13525./55296., dc5=0.-277./14336.,
+ const double dc4=125./594.-13525./55296., dc5=0.-277./14336.,
dc6=512./1771.-0.25;
// Coefficients for diff. equations:
@@ -172,9 +160,9 @@ class CRunDec
double fMsFromOs4(double mu, double M, double nl, double err);
double fZmM(double n);
double fZmInvM(double n);
- double deltamOS2mMS(double mOS, RunDecPair mq[],
+ double deltamOS2mMS(double mOS, std::pair* mq,
double asmu, double mu, int nlq, int nloops);
- double deltamMS2mOS(double mMS, RunDecPair mq[],
+ double deltamMS2mOS(double mMS, std::pair* mq,
double asmu, double mu, int nlq, int nloops);
double fMsFromRi1(void);
double fMsFromRi2(void);
@@ -185,7 +173,7 @@ class CRunDec
double fMumFromOs4(double err);
double fRiFromMs(double alpha, double nl);
double fMsFromRi(double alpha, double nl);
- double fHelpmOS2mMSit(double mMS,double mOS, RunDecPair mq[],
+ double fHelpmOS2mMSit(double mMS,double mOS, std::pair* mq,
double asmu, double mu, int nl);
double fas5to6os(double alpha, double mass, double mu, double nlq, double nl);
double fas6to5os(double alpha, double mass, double mu, double nlq, double nl);
@@ -204,7 +192,7 @@ class CRunDec
double E1p(double mOS, double asmu, double mu, int nl, int nloops);
double exOS2RS(double api, double mmu, double nnuf, int nnl, int nloops);
double exOS2RSp(double api, double mmu, double nnuf, int nnl, int nloops);
- double mMS2mOSmod(double mMS, RunDecPair mq[],
+ double mMS2mOSmod(double mMS, std::pair* mq,
double asmu, double mu, int nf, int nloops, double err);
double mkin2mMSA(double mkin, double apinlmus, double mus, double mufac, int NLMSOS, int NLOSKIN, double mcMSmusmcin, double musmcin, int nloops);
double mkin2mMSB(double mkin, double apinlmus, double mus, double mufac, int NLMSOS, int NLOSKIN, double mcMSmusmcin, double musmcin, int nloops);
@@ -228,8 +216,7 @@ class CRunDec
CRunDec(int);
// Arrays and structs to store data:
- //std::pair mq[4];
- RunDecPair mq[4];
+ std::pair mq[4];
TriplenfMmu nfMmu[4];
AsmMS AM;
@@ -492,7 +479,7 @@ class CRunDec
\param nloops number of loops
\return \f$\alpha_s(\mu_2)\f$
*/
- double AlL2AlH(double asl, double mu1, TriplenfMmu *decpar, double mu2, int nloops);
+ double AlL2AlH(double asl, double mu1, TriplenfMmu decpar[], double mu2, int nloops);
//! AlH2AlL calculates \f$\alpha_s(\mu_2)\f$ from \f$\alpha_s(\mu_1)\f$ decoupling at intermediate scales, running from high to low
/*!
@@ -545,7 +532,7 @@ class CRunDec
\return \f$M_{OS}\f$
*/
- double mMS2mOS(double mMS, RunDecPair mq[],
+ double mMS2mOS(double mMS, std::pair* mq,
double asmu, double mu,int nf, int nloops, double fdelm=1.0);
//! mOS2mMS calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$M_{OS}\f$
@@ -559,7 +546,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_f)}(\mu)\f$
*/
- double mOS2mMS(double mOS, RunDecPair mq[],
+ double mOS2mMS(double mOS, std::pair* mq,
double asmu, double mu,int nf, int nloops, double fdelm=1.0);
//! mMS2mSI calculates \f$m_{MS}^{(n_f)}(m_{MS})\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
@@ -615,7 +602,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_f)}(m_{MS})\f$
*/
- double mOS2mSI(double mOS, RunDecPair mq[],
+ double mOS2mSI(double mOS, std::pair* mq,
double asM, int nf, int nloops, double fdelm=1.0);
//! mOS2mMSrun calculates \f$m_{MS}^{(n_f)}(\mu)\f$ from \f$M_{OS}\f$
@@ -630,7 +617,7 @@ class CRunDec
\param nloops number of loops
\return \f$m_{MS}^{(n_f)}(\mu)\f$
*/
- double mOS2mMSrun(double mOS, RunDecPair mq[],
+ double mOS2mMSrun(double mOS, std::pair* mq,
double asmu, double mu,int nf, int nloops);
//! mMS2mOSrun calculates \f$M_{OS}\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
@@ -645,7 +632,7 @@ class CRunDec
\param nloops number of loops
\return \f$M_{OS}\f$
*/
- double mMS2mOSrun(double mMS, RunDecPair mq[],
+ double mMS2mOSrun(double mMS, std::pair* mq,
double asmu, double mu,int nf, int nloops);
//! mMS2mRI calculates \f$m^{RI}(\mu)\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
@@ -670,7 +657,7 @@ class CRunDec
\param nloops number of loops
\return \f$m_{MS}^{(n_f)}(\mu)\f$
*/
- double mOS2mMSit(double mOS, RunDecPair mq[],
+ double mOS2mMSit(double mOS, std::pair* mq,
double asmu, double mu, int nf, int nloops);
//! mMS2mRGImod calculates \f$m^{RGI}\f$ from \f$m_{MS}^{(n_f)}(\mu)\f$
@@ -697,7 +684,7 @@ class CRunDec
\param nloops number of loops
\return \f$m^{PS}(\mu_f)\f$
*/
- double mOS2mPS(double mOS, RunDecPair mq[],
+ double mOS2mPS(double mOS, std::pair* mq,
double asmu, double mu, double muf, int nl, int nloops);
//! mMS2mPS calculates \f$m^{PS}(\mu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
@@ -712,7 +699,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m^{PS}(\mu_f)\f$
*/
- double mMS2mPS(double mMS, RunDecPair mq[],
+ double mMS2mPS(double mMS, std::pair* mq,
double asmu, double mu, double muf, int nl, int nloops, double fdelm=1.0);
//! mPS2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{PS}(\mu_f)\f$
@@ -727,7 +714,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(\mu)\f$
*/
- double mPS2mMS(double mPS, RunDecPair mq[],
+ double mPS2mMS(double mPS, std::pair* mq,
double asmu, double mu, double muf, int nl, int nloops, double fdelm=1.0);
//! mPS2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{PS}(\mu_f)\f$
@@ -741,7 +728,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
*/
- double mPS2mSI(double mPS, RunDecPair mq[],
+ double mPS2mSI(double mPS, std::pair* mq,
double (*as)(double), double muf, int nl, int nloops, double fdelm=1.0);
//! mOS2m1S calculates \f$m^{1S}\f$ from \f$M_{OS}\f$
@@ -754,7 +741,7 @@ class CRunDec
\param nloops number of loops
\return \f$m^{1S}\f$
*/
- double mOS2m1S(double mOS, RunDecPair mq[],
+ double mOS2m1S(double mOS, std::pair* mq,
double asmu, double mu, int nl, int nloops);
//! mMS2m1S calculates \f$m^{1S}\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
@@ -768,7 +755,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m^{1S}\f$
*/
- double mMS2m1S(double mMS, RunDecPair mq[],
+ double mMS2m1S(double mMS, std::pair* mq,
double asmu, double mu, int nl, int nloops, double fdelm=1.0);
//! m1S2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{1S}\f$
@@ -782,7 +769,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(\mu)\f$
*/
- double m1S2mMS(double m1S, RunDecPair mq[],
+ double m1S2mMS(double m1S, std::pair* mq,
double asmu, double mu, int nl, int nloops, double fdelm=1.0);
//! m1S2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{1S}\f$
@@ -795,7 +782,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
*/
- double m1S2mSI(double m1S, RunDecPair mq[],
+ double m1S2mSI(double m1S, std::pair* mq,
double (*as)(double), int nl, int nloops, double fdelm=1.0);
//! mOS2mRS calculates \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$ from \f$M_{OS}\f$
@@ -812,7 +799,7 @@ class CRunDec
\param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ should be calculated
\return \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
*/
- double mOS2mRS(double mOS, RunDecPair mq[], double asmu,
+ double mOS2mRS(double mOS, std::pair* mq, double asmu,
double mu, double nuf, int nl, int nloops, bool prime);
//! mMS2mRS calculates \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
@@ -830,7 +817,7 @@ class CRunDec
\param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ should be calculated
\return \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
*/
- double mMS2mRS(double mMS, RunDecPair mq[], double asmu,
+ double mMS2mRS(double mMS, std::pair* mq, double asmu,
double mu, double nuf, int nl, int nloops, double fdelm, bool prime);
//! mRS2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
@@ -848,7 +835,7 @@ class CRunDec
\param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ is given
\return \f$m_{MS}^{(n_l+1)}(\mu)\f$
*/
- double mRS2mMS(double mRS, RunDecPair mq[], double asmu,
+ double mRS2mMS(double mRS, std::pair* mq, double asmu,
double mu, double nuf, int nl, int nloops, double fdelm, bool prime);
//! mRS2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{RS}(\nu_f)\f$ or \f$m^{RS'}(\nu_f)\f$
@@ -865,7 +852,7 @@ class CRunDec
\param prime selects if \f$m^{RS}\f$ or \f$m^{RS'}\f$ is given
\return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
*/
- double mRS2mSI(double mRS, RunDecPair mq[], double (*as)(double),
+ double mRS2mSI(double mRS, std::pair* mq, double (*as)(double),
double nuf, int nl, int nloops, double fdelm, bool prime);
@@ -881,7 +868,7 @@ class CRunDec
\param nloops number of loops
\return \f$m^{RS}(\nu_f)\f$
*/
- double mOS2mRS(double mOS, RunDecPair mq[],
+ double mOS2mRS(double mOS, std::pair* mq,
double asmu, double mu, double nuf, int nl, int nloops);
//! mMS2mRS calculates \f$m^{RS}(\nu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
@@ -897,7 +884,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m^{RS}(\nu_f)\f$
*/
- double mMS2mRS(double mMS, RunDecPair mq[],
+ double mMS2mRS(double mMS, std::pair* mq,
double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
//! mRS2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{RS}(\nu_f)\f$
@@ -913,7 +900,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(\mu)\f$
*/
- double mRS2mMS(double mRS, RunDecPair mq[],
+ double mRS2mMS(double mRS, std::pair* mq,
double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
//! mRS2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{RS}(\nu_f)\f$
@@ -928,7 +915,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
*/
- double mRS2mSI(double mRS, RunDecPair mq[],
+ double mRS2mSI(double mRS, std::pair* mq,
double (*as)(double), double nuf, int nl, int nloops, double fdelm=1.0);
//! mOS2mRSp calculates \f$m^{RS'}(\nu_f)\f$ from \f$M_{OS}\f$
@@ -943,7 +930,7 @@ class CRunDec
\param nloops number of loops
\return \f$m^{RS'}(\nu_f)\f$
*/
- double mOS2mRSp(double mOS, RunDecPair mq[],
+ double mOS2mRSp(double mOS, std::pair* mq,
double asmu, double mu, double nuf, int nl, int nloops);
//! mMS2mRSp calculates \f$m^{RS'}(\nu_f)\f$ from \f$m_{MS}^{(n_l+1)}(\mu)\f$
@@ -959,7 +946,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m^{RS'}(\nu_f)\f$
*/
- double mMS2mRSp(double mMS, RunDecPair mq[],
+ double mMS2mRSp(double mMS, std::pair* mq,
double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
//! mRSp2mMS calculates \f$m_{MS}^{(n_l+1)}(\mu)\f$ from \f$m^{RS'}(\nu_f)\f$
@@ -975,7 +962,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(\mu)\f$
*/
- double mRSp2mMS(double mRS, RunDecPair mq[],
+ double mRSp2mMS(double mRS, std::pair* mq,
double asmu, double mu, double nuf, int nl, int nloops, double fdelm=1.0);
//! mRSp2mSI calculates \f$m_{MS}^{(n_l+1)}(m_{MS})\f$ from \f$m^{RS'}(\nu_f)\f$
@@ -990,7 +977,7 @@ class CRunDec
\param fdelm factor multiplying the non-logarithmic part of the 4-loop term
\return \f$m_{MS}^{(n_l+1)}(m_{MS})\f$
*/
- double mRSp2mSI(double mRS, RunDecPair mq[],
+ double mRSp2mSI(double mRS, std::pair* mq,
double (*as)(double), double nuf, int nl, int nloops, double fdelm=1.0);
@@ -1009,8 +996,8 @@ class CRunDec
\param deccase Specify if there is a light but massive quark
\return \f$m^{\mathrm{kin}}\f$
*/
- double mMS2mKIN(double mMS, RunDecPair mq[],
- double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, int deccase);
+ double mMS2mKIN(double mMS, std::pair* mq,
+ double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, std::string deccase);
//! mKIN2mMS calculates \f$m_{MS}^{(n_f)}(\mu_s)\f$ from \f$m^{\mathrm{kin}}\f$
@@ -1026,8 +1013,8 @@ class CRunDec
\param deccase Specify if there is a light but massive quark
\return \f$m_{MS}^{(n_f)}(\mu_s)\f$
*/
- double mKIN2mMS(double mKIN, RunDecPair mq[],
- double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, int deccase);
+ double mKIN2mMS(double mKIN, std::pair* mq,
+ double asmus, double mus, double muf, int nlmsos, int nloskin, int nloops, std::string deccase);
// Overload functions:
double LamExpl(double asmu, double mu, int nloops);
@@ -1053,28 +1040,28 @@ class CRunDec
double DecMqDownSI(double mq, double asmu, double massth, double muth, int nloops);
- double mMS2mOS(double mMS, RunDecPair mq[],
+ double mMS2mOS(double mMS, std::pair* mq,
double asmu, double mu, int nloops, double fdelm=1.0);
- double mOS2mMS(double mOS, RunDecPair mq[],
+ double mOS2mMS(double mOS, std::pair* mq,
double asmu, double mu, int nloops, double fdelm=1.0);
double mMS2mSI(double mMS, double asmu, double mu, int nloops);
double mRI2mMS(double mRI, double asmu, int nloops);
double mMS2mRGI(double mMS, double asmu, int nloops);
double mRGI2mMS(double mRGI, double asmu, int nloops);
- double mOS2mSI(double mOS, RunDecPair mq[],
+ double mOS2mSI(double mOS, std::pair* mq,
double asM, int nloops, double fdelm=1.0);
- double mOS2mMSrun(double mOS, RunDecPair mq[],
+ double mOS2mMSrun(double mOS, std::pair* mq,
double asmu, double mu, int nloops);
- double mMS2mOSrun(double mMS, RunDecPair mq[],
+ double mMS2mOSrun(double mMS, std::pair* mq,
double asmu, double mu, int nloops);
double mMS2mRI(double mMS, double asmu, int nloops);
- double mOS2mMSit(double mOS, RunDecPair mq[],
+ double mOS2mMSit(double mOS, std::pair* mq,
double asmu, double mu, int nloops);
double mMS2mRGImod(double mMS, double asmu, int nloops);
- double mMS2mKIN(double mMS, RunDecPair mq[], double asmus, double mus, double muf, int nloops, int deccase);
- double mKIN2mMS(double mKIN, RunDecPair mq[], double asmus, double mus, double muf, int nloops, int deccase);
+ double mMS2mKIN(double mMS, std::pair* mq, double asmus, double mus, double muf, int nloops, std::string deccase);
+ double mKIN2mMS(double mKIN, std::pair