int main(int n, char **args) { double alphaSMZ,mbp,mbmb,mtp; double Q; int nf; double McMc=1.4; sscanf(args[1],"%lf",&alphaSMZ); sscanf(args[2],"%lf",&McMc); sscanf(args[3],"%lf",&mbmb); sscanf(args[4],"%lf",&mtp); sscanf(args[5],"%lf",&Q); initQCD(alphaSMZ,McMc,mbmb,mtp); printf("qMass : %E %E %E %E\n", qMass[3], qMass[4], qMass[5], qMass[6]); printf("lambda: %E %E %E %E\n", lambda[3],lambda[4],lambda[5],lambda[6]); printf("qMin=%E\n",qMin); printf("MbPole=%f\n", MbPole); printf("MtMt=%f\n", qMass[6]); printf("qmass[6]=%f\n",qMass[6]); printf("alphaS(%f)=%f\n",Q,alphaQCD(Q)); printf("MbRun=%f MbEff=%f \n",MbRun(Q),MbEff(Q) ); printf("MtRun=%f MtEff=%f \n",MtRun(Q),MtEff(Q) ); }
double initQCD5(double MZalphaS, double McMc, double MbMb, double MtP) { double lmbd= initQCD(MZalphaS,McMc, MbMb,MtP); lambda[6]=lambda[5]; nfMax=5; return lmbd; }
double alphaQCD(double Q) { if(notInitialized) initQCD(0.1172,1.2,4.23,171.4); if(Q<qMin) return 1; return alpha3(Q,lambda[NF(Q)],NF(Q)); }
double initqcd_(double * alfsMZ, double * McMc, double * MbMb, double * Mtp) { return initQCD(*alfsMZ,*McMc,*MbMb,*Mtp); }
double alphaQCD(double Q) { if(notInitialized) initQCD(0.1184,1.2,4.23,173.07); if(Q<qLim) Q=qLim; if(Q<qMin) return 1; return alpha3(Q,lambda[NF(Q)],NF(Q)); }