//-------------------------------------------
//   calculate QCD coupling constant alphaS()
//--------------------------------------------
double alphaS(double Q) {
return 12.*3.141593/(33-2*5)/log(Q*Q/0.1/0.1);
}