static int input(int y0, char*hlp, char*directive, char*text, int cur, int lim)
{
if(errTxt)
{
goto_xy(1,y0+1);
scrcolor(FGmain,BGmain);
clr_eol();
print("%s","Error: ");
scrcolor(Red,BGmain);
print("%s",errTxt);
scrcolor(FGmain,BGmain);
if(blind) { printf("%s\n",errTxt); sortie(110);}
be_be();
errTxt=NULL;
}
for(;;)
{ int rc;
goto_xy(1,y0);
scrcolor(FGmain,BGmain);
print("%s",directive);
rc=str_redact(text,cur, lim);
switch(rc)
{ case KB_UP:
case KB_DOWN:
case KB_PAGED:
case KB_PAGEU: prtcllist(rc); continue;
case KB_F1: show_help(hlp); continue;
}
return rc;
}
}
static void f8_key_prog(int x)
{
static char FUNC[75]="2*2 % Press ESC to finish, F1 for help, ENTER to calculate";
int npos=1;
void * pscr;
get_text(1,20,80,21,&pscr);
scrcolor(Red,White);
goto_xy(1,20); print("CALC : ");
goto_xy(1,21); print("result=");
scrcolor(Black,White);
for(;;)
{ double res;
int err;
int key;
goto_xy(9,20); key=str_redact(FUNC,npos,70);
if(key==KB_ESC) break;
else if(key==KB_F1) show_help("n_calc");
else if(key==KB_ENTER)
{ goto_xy(8,21);
err=calcExpression(FUNC,rd_num,&res);
goto_xy(9,21);
if(err) {print("Erorr: %s",errmesstxt(err)); npos=rderrpos;}
else print("%E ",res);
}
}
put_text(&pscr);
}
static int correctVar(char* txt,int x,int y, int type, void * var,int clear)
{ int npos,key,err;
int xx;
void * pscr;
char * buff;
int * I;
long * L;
double * D;
int maxLen;
get_text(x,y,maxCol(),y,&pscr);
scrcolor(White,Black);
goto_xy(x,y);
print(txt);
xx=where_x();
if(type<0) {buff=var; maxLen=-type;} else
{ buff=malloc(100);
switch (type)
{
case 'L': L=var; sprintf(buff,"%ld",*L); break;
case 'D': D=var; sprintf(buff,"%lg",*D); break;
case 'I': I=var; sprintf(buff,"%d" ,*I); break;
}
maxLen=MAX(15,strlen(buff));
}
npos = 1;
do
{
goto_xy(xx,y);
key = str_redact(buff,npos,maxLen);
if (key == KB_ESC)
{
put_text(&pscr);
if(type>=0) free(buff);
return 0;
}
err=0;
if (key == KB_ENTER)
{
/* trim(buff); */
if(type<0) err=1; else
{
switch (type)
{
case 'L': err=sscanf(buff,"%ld",L); break;
case 'D': err=sscanf(buff,"%lf",D); break;
case 'I': err=sscanf(buff,"%d" ,I); break;
}
}
if (err<=0) messanykey(10,10," incorrect number");
}
} while (err<=0);
if(clear)put_text(&pscr); else free(pscr);
if(type>=0) free(buff);
return 1;
}
int viewresults(void)
{
int k,kmenu;
void * pscr = NULL;
shortstr newname;
int dirStat=checkDir("results");
if(dirStat==0){messanykey(10,15,"directory RESULTS is empty"); return 1;}
kmenu = 1;
label_1:
menu1(10,10,"","\010"
" View "
" Delete "
" Rename ","s_res",&pscr,&kmenu);
switch (kmenu)
{
case 0:
return 0;
case 1:
viewDir("results");
break;
case 2:
if(dirStat==2)
{ char mess[]="Can not clean dir 'results' because it contains the LOCK file";
if(blind) { printf("%s\n",mess); sortie(102);}
else { messanykey(3,13,mess); break;}
}
if ( mess_y_n( 6,13," Delete files ") )
{ struct dirent **namelist;
int n,i;
n = scandir("./results", &namelist, NULL, NULL);
for(i=0; i<n;i++)
{
char buff[100];
if(strcmp(namelist[i]->d_name,"aux") && strcmp(namelist[i]->d_name,"..") && strcmp(namelist[i]->d_name,".") )
{
sprintf(buff,"rm -rf results/%s",namelist[i]->d_name);
if(unlink(buff+7)) system(buff);
}
}
free(namelist);
}
put_text(&pscr);
return 1;
case 3:
strcpy(newname," ");
while(1)
{ void * pscr3;
get_text(1,maxRow(),maxCol(),maxRow(),&pscr3);
goto_xy(1,maxRow());
print("Enter new name: ");
k = str_redact(newname,1,30);
if (k == KB_ESC)
{ goto_xy(1,24);
clr_eol();
goto label_1;
}
if (k == KB_ENTER)
{
trim(newname);
if(rename("results",newname)==0)
{ char command[200];
mkdir("results",0755);
sprintf(command," cp -rp %s/aux results",newname);
system(command);
put_text(&pscr);
put_text(&pscr3);
return 1;
}
else messanykey(10,15," Can't rename the directory");
}
put_text(&pscr3);
}
}
goto label_1;
}
int qcdmen_(void)
{
void * pscr=NULL;
int mode;
int returnCode=0;
initStrFun(0);
L10:{ char strmen[]="\030"
" parton dist. alpha OFF "
" alpha(MZ)= ZZZZ "
" nf = NF "
" order= NNLO "
" mb(mb)= MbMb "
" Mtop(pole)= Mtp "
" Alpha(Q) plot "
" Qren = RRR "
" Qpdf1= FF1 "
" Qpdf2= FF2 "
" Qshow= FFS ";
if(alphaPDF)
{ int k=0;
//printf("alphaPDF=%d sf_alpha[0]=%p sf_alpha[1]=%p\n", alphaPDF,sf_alpha[0],sf_alpha[1]);
switch(alphaPDF)
{ case 1: if(sf_alpha[0]) k=1; else if(sf_alpha[1]) k=2; break;
case 2: if(sf_alpha[1]) k=2; else if(sf_alpha[0]) k=1; break;
}
if(k) improveStr(strmen,"OFF","pdf%d",k);
}
improveStr(strmen,"ZZZZ","%.4f", alphaMZ);
improveStr(strmen,"NF","%d",alphaNF);
if(alphaOrder==1) improveStr(strmen,"NNLO","%-.4s","LO");
else if(alphaOrder==2) improveStr(strmen,"NNLO","%-.4s","NLO");
else alphaOrder=3;
improveStr(strmen,"MbMb","%.3f", MbMb);
improveStr(strmen,"Mtp","%.2f", Mtp);
improveStr(strmen,"RRR","%-.16s", Rscale_str);
improveStr(strmen,"FF1","%-.16s", F1scale_str);
improveStr(strmen,"FF2","%-.16s", F2scale_str);
improveStr(strmen,"FFS","%-.16s", Sscale_str);
menu1(54,8,"QCD alpha",strmen,"n_alpha",&pscr,&mode);
}
switch (mode)
{ case 0: if(returnCode) init_alpha(); return returnCode;
case 1:
{ char alphaMen[100]="\006"
" OFF "
" pdf1 "
" pdf2 ";
int k=0;
menu1(54,12,"alpha",alphaMen,"",NULL,&k);
if(k)alphaPDF=k-1;
if(alphaPDF && !sf_alpha[alphaPDF-1]) messanykey(20,20,"WARNING! This pdf does not define alphaQCD ");
}
break;
case 2:
{ double alphaMZ_old=alphaMZ;
if(correctDouble(3,15,"Enter new value ",&alphaMZ,1)) returnCode=1;
if(alphaMZ>0 && alphaMZ<0.3) returnCode=1; else
{ alphaMZ=alphaMZ_old;
messanykey(5,15,"Your input is out of alphaMZ range");
}
}
break;
case 3:
{ int NF_old=alphaNF;
if(correctInt(3,15,"Enter new value ",&alphaNF,1))
{
if(alphaNF<=6 && alphaNF>=3) returnCode=1;
else { messanykey(5,15,"NF out of range"); alphaNF=NF_old;}
}
}
break;
case 4:
{ char lomen[]="\010"
" LO "
" NLO "
" NNLO ";
void *pscrlo=NULL;
int k=0;
menu1(52,12,"",lomen,"",&pscrlo,&k);
if(k) { alphaOrder=k; returnCode=1; put_text(&pscrlo);}
}
break;
case 5: correctDouble(3,15,"Enter new value ",&MbMb,1); break;
case 6: correctDouble(3,15,"Enter new value ",&Mtp,1); break;
case 7:
{ void * screen;
int i;
static double qMin=1, qMax=1000;
static int nPoints=100;
if(returnCode) init_alpha();
get_text(1,1,maxCol(),maxRow(),&screen);
if(correctDouble(40 ,15 ,"Q_min=",&qMin,0)&& qMin>=0.5
&& correctDouble(40 ,16 ,"Q_max=",&qMax,0)&& qMax>qMin
&& correctInt(33,17,"number of points=" ,&nPoints,0)
&& nPoints>3&& nPoints<=150)
{ double *f[3]={NULL,NULL,NULL};
double *ff[3]={NULL,NULL,NULL};
char buff[3][100];
char* Y[3]={buff[0],buff[1],buff[2]};
switch(alphaOrder)
{ case 1: sprintf(Y[0],"MSbar LO"); break;
case 2: sprintf(Y[0],"MSbar NLO"); break;
case 3: sprintf(Y[0],"MSbar NNLO"); break;
default:sprintf(Y[0],"MSbar");
}
int N,k;
f[0]=(double*) malloc(nPoints*sizeof(double));
int xLog= (qMin>0 && qMax/qMin >10)? 1 : 0;
for(i=0;i<nPoints;i++)
{ double z=(i+0.5)/(double)(nPoints),q;
if(xLog) q=pow(qMin,1-z)*pow(qMax,z); else q=qMin*(1-z)+qMax*z;
f[0][i]=alpha_0(q);
// printf("i=%d %E\n",i,f[0][i]);
}
N=1;
for(k=0;k<2;k++) if(sf_alpha[k])
{
char buff[300];
strFunName(k+1,Y[N]);
char *p=strstr(Y[N],"(proton");
if(p) p[0]=0; else
{ p=strstr(Y[N],"(anti-proton");
if(p) p[0]=0;
}
f[N]=(double*) malloc(nPoints*sizeof(double));
for(i=0;i<nPoints;i++)
{ double z=(i+0.5)/(double)(nPoints),q;
if(xLog) q=pow(qMin,1-z)*pow(qMax,z);
else q=qMin*(1-z)+qMax*z;
f[N][i]=(*sf_alpha[k])(q);
}
N++;
}
// printf("N=%d Y[0]=%s\n Y[1]=%s\n Y[2]=%s\n", N,Y[0],Y[1],Y[2]);
// plot_Nar(NULL, "Alpha(Q)", log10(qMin), log10(qMax),"log10(Q/GeV)", nPoints, N, f,ff,Y);
plot_Nar(NULL, "Alpha(Q)", qMin,qMax,"Q/GeV", nPoints, xLog, N, f,ff,Y);
for(k=0;k<N;k++) free(f[k]);
} else messanykey(40,18,
" Correct input is \n"
" 0.5<= Q_min <Q_max\n"
" number of points <=150 and >=4");
put_text(&screen);
} break;
case 8:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Renorm. scale: ");
if(str_redact(Rscale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str, mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
case 9:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Fct1.scale: ");
if(str_redact(F1scale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str, mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
case 10:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Fct1.scale: ");
if(str_redact(F2scale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str, mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
case 11:
{ int npos=1,rc;
do
{
char mess[200];
goto_xy(2,12); print("Shworing scale: ");
if(str_redact(Sscale_str,npos,60)==KB_ENTER) returnCode=1;
goto_xy(2,12); clr_eol();
rc=initScales(Rscale_str,F1scale_str,F2scale_str,Sscale_str,mess);
if(rc) messanykey(10,10,mess);
} while(rc);
}
break;
}
goto L10;
}
int viewresults(void)
{
int k,kmenu;
void * pscr = NULL;
shortstr newname;
int dirStat=checkDir("results");
if(dirStat==0){messanykey(10,15,"directory RESULTS is empty"); return 1;}
kmenu = 1;
label_1:
menu1(10,10,"","\010"
" View "
" Delete "
" Rename ","s_res",&pscr,&kmenu);
switch (kmenu)
{
case 0:
return 0;
case 1:
viewDir("results");
break;
case 2:
if(dirStat==2)
{ char mess[]="Can not clean dir 'results' because it contains the LOCK file";
if(blind) { printf("%s\n",mess); sortie(102);}
else { messanykey(3,13,mess); break;}
}
if ( mess_y_n( 6,13," Delete files ") ) system("rm -r results; mkdir results");
put_text(&pscr);
return 1;
case 3:
strcpy(newname," ");
while(1)
{ void * pscr3;
get_text(1,maxRow(),maxCol(),maxRow(),&pscr3);
goto_xy(1,maxRow());
print("Enter new name: ");
k = str_redact(newname,1,30);
if (k == KB_ESC)
{ goto_xy(1,24);
clr_eol();
goto label_1;
}
if (k == KB_ENTER)
{
trim(newname);
if(rename("results",newname)==0)
{
mkdir("results",-1);
put_text(&pscr);
put_text(&pscr3);
return 1;
}
else messanykey(10,15," Can't rename the directory");
}
put_text(&pscr3);
}
}
goto label_1;
}
