double McEff(double Q) { return McRun(Q)*sqrt(1+DeltaQCD(Q)); }
double mcrun_(double *q){ return McRun(*q);}