static int mappingMenu(void)
{
int r=0;
int mode=1;
void * pscr=NULL;
char menutxt[]="\030"
" Kinematics "
" Regularization ";
for(;;)
{
menu1(54,6,"",menutxt,"n_map_*",&pscr, &mode);
switch(mode)
{ case 0: return r;
case 1: r=r|entkin_(); break;
case 2: do r=r|(2*edittable(1,4,®Tab,1,"n_reg",0));while(fillRegArray());
break;
}
}
}
int monte_carlo_menu(void)
{
static int r=0;
int mode=1;
void * pscr=NULL;
void (*quit)(int)=f3_key[7];
char menutxt[]="\030"
" Subprocess "
" IN state "
" Model parameters "
" Constraints "
" QCD alpha & scales "
" Breit-Wigner "
" Aliases "
" Cuts "
" Phase space mapping "
" Monte Carlo simulation "
" Easy ";
if(nout_int!=2 ) menutxt[menutxt[0]*10+1]=0;
if(nin_int==1) improveStr(menutxt,"Easy", "Total width");
else improveStr(menutxt,"Easy", "1D integration");
wrtprc_();
for(;;)
{
infor();
f3_key[7]=quit;
menu1(54,4,"",menutxt,"n_mc_*",&pscr, &mode);
if(mode==1||mode==2||mode==3||mode==5||mode==7)f3_key[7]=NULL;
switch (mode)
{
case 0: return 0;
case 1: r=r|3*sub_men__(); break;
case 2: r=r|in_setting(); break;
case 3: r=r|change_parameter(54,7,0); break;
case 4: { int modeC=1;
for(;;)
{ char menuC[]="\030"
" All Constraints "
" Masses,Widths,Branching";
void * pscrC=NULL;
menu1(54,6,"",menuC,"n_constr_*",&pscrC, &modeC);
switch(modeC)
{ case 0: break;
case 1: show_depend(54,7); break;
case 2: show_spectrum(54,9); break;
}
if(!modeC) break;
} break;
}
case 5: r=r|qcdmen_(); break;
case 6: r=r|w_men__(); break;
case 7: do r=r|(3*edittable(1,4,&compTab,1,"n_comp",0)); while (fillCompositeArray());
break;
case 8: do r=r|(3*edittable(1,4,&cutTab,1,"n_cut",0)); while (fillCutArray());
break;
case 9: r=r|mappingMenu(); break;
case 10:
if(nout_int==1 && !sf_num[0] && !sf_num[1] )
{ if(blind) return 1;
messanykey(15,15,"Phase space integration for 2->1 processes\n needs distribution functions.");
break;
}
if(checkEnergy())
{
if(blind==1)
{ char fname[50];
int i,j;
sprintf(fname,"events_%d.txt", nSess);
FILE * f=fopen(fname,"w");
fprintf(f,"#%s\n", VERSION_);
fprintf(f,"#Type %d -> %d\n", nin_int,nout_int);
fprintf(f,"#Initial_state ");
if(nin_int==1) fprintf(f," P1=0\n");
else
{ fprintf(f," P1_3=0 P2_3=0\n");
wrt_sf__(f);
}
fprintf(f,"#PROCESS ");
for(i=1;i<=nin_int+nout_int; i++)
{ int pcode;
char * pname=pinf_int(Nsub,i,NULL,&pcode);
fprintf(f," %d(%s)", pcode, pname);
if(i==nin_int) fprintf(f," ->");
}
fprintf(f,"\n");
fprintf(f,"#MASSES ");
for(i=0;i<nin_int+nout_int;i++)
{ REAL m;
pinf_int(Nsub,i+1,&m,NULL);
fprintf(f," %.10E", (double)m);
}
fprintf(f,"\n");
fprintf(f,"#Cross_section(Width) %E\n",0.);
fprintf(f,"#Number_of_events %10d\n",0);
fprintf(f,"#Sum_of_weights %12.4E %12.4E \n",0.,0.);
fprintf(f,"#Events ");
if(nin_int==2) fprintf(f," P1_3 [Gev] P2_3 [Gev] ");
for(i=1;i<=nout_int; i++) for(j=1;j<=3;j++)
fprintf(f," P%d_%d [Gev] ",i+nin_int,j);
integral.old=1;
fclose(f);
return 1;
}
messanykey(15,15,"Energy is too small!");
break;
}
if(fillCutArray())
{ if(blind) return 2;
messanykey(15,15,"Can not evaluate cut limits");
break;
}
case 11:
if(mode==11)
{ void (*f10_tmp)(int);
w_sess__(NULL);
f10_tmp=f3_key[7];
f3_key[7]=f10_key_prog_for22;
if(nin_int==1) decay12(); else
{ REAL m1,m2, Pcm;
pinf_int(Nsub,1,&m1,NULL);
pinf_int(Nsub,2,&m2,NULL);
if(sf_num[0] && sf_mass[0]>m1) m1= sf_mass[0];
if(sf_num[1] && sf_mass[1]>m2) m2= sf_mass[1];
incomkin(m1,m2,inP1,inP2,NULL,&Pcm,NULL);
if(sf_num[0]||sf_num[1]||nCuts)
messanykey(10,10,"Structure functions and cuts are ignored\n");
cs_numcalc(Pcm);
}
f3_key[7]= f10_tmp;
r_sess__(NULL);
break;
} else if(fillRegArray())
{
if(blind) return 3;
messanykey(15,15,
"Can not evaluate regularization paremeters");
break;
}
if(mode==10) runVegas();
r=0;
break;
}
//printf("r=%d\n",r);
if(r) clearEventMax();
if(r&2) clearGrid();
if(r&1)newSession();
}
}
int main(int argc,char** argv)
{
/*===================================
0 - First start.
1 - Model menu.
2 - Enter process menu.
3 - Feynman diagrams menu; squaring.
4 - Squared diagram menu; symbolic calculation.
5 - Write results; new process.
10 -Restart symbolic calculations
==========================================================*/
/* 0-Start; 1-Restart; 2-Heap Error,3-Edit Model,4-UserBreak */
void *pscr1=NULL,*pscr2=NULL,*pscr3=NULL,*pscr4=NULL,*pscr5=NULL;
int k1=1,k2=1,k3=1,k4=1,k5=1;
int n;
int pid=0;
char LOCKtxt[]="Directory 'results/' contains the .lock file created by n_calchep.\n"
"To continue you may a)Close the n_calchep session, or b)Rename 'results/',\n";
blind=0;
strcpy(pathtocomphep,argv[0]);
for(n=strlen(pathtocomphep)-1; n && pathtocomphep[n]!=f_slash; n--);
pathtocomphep[n-3]=0;
for ( n=1;n<argc;n++)
{ if (strcmp(argv[n],"-blind")==0 && n<argc-1 )
{ blind=1;
inkeyString=argv[++n];
}
if (strcmp(argv[n],"+blind")==0 ) blind=2;
}
if(!writeLockFile(".lock"))
{ fprintf(stderr,"locked by other s_calchep. See .lock file\n");
exit(100);
}
strcpy(pathtouser,"");
sprintf(pathtohelp,"%shelp%c",pathtocomphep,f_slash);
outputDir="results/";
{ char * icon=(char *) malloc(strlen(pathtocomphep)+10);
int err;
sprintf(icon,"%sicon",pathtocomphep);
err=start1("CalcHEP/symb",icon,"calchep.ini",&xw_error);
if(err && blind==0)
{ printf("Error:You have launched interactive session for version compiled without X11 library.\n");
printf("Options: a) Use blind session; b) Recompile CalcHEP with X11\n");
exit(66);
}
free(icon);
}
fillModelMenu();
f3_key[0]=f3_key_prog; f3_mess[0]="Model";
f3_key[1]=f4_key_prog; f3_mess[1]="Diagrams";
f3_key[2]=f5_key_prog; f3_mess[2]="Switches";
f3_key[3]=f6_key_prog; f3_mess[3]="Results";
f3_mess[4]="Del";
f3_mess[5]="UnDel";
f3_key[6]=f9_key_prog; f3_mess[6]="Ref";
f3_key[7]=f10_key_prog; f3_mess[7]="Quit";
restoreent(&menulevel);
if(!blind && menulevel<2) cheplabel();
switch (menulevel)
{
case 10:
case 6:
case 5:
case 4:
case 3: readModelFiles("./models",n_model);
modelinfo();
loadModel(0,forceUG);
processinfo();
diagramsinfo();
case 2: k1=n_model;
case 1: break;
}
switch (menulevel)
{
case 2: menuhelp();
goto label_20;
case 3: goto label_31;
case 4: goto label_40;
case 5: goto label_50;
case 6:
case 10: goto restart2;
}
label_10: /* Menu2(ModelMenu): */
f3_key[0]=NULL; /*models*/
f3_key[1]=NULL; /*diagrams*/
menulevel = 1;
forceUG=0;
menuhelp();
for(;;)
{
showheap();
k1=n_model;
menu1(56,4,"",modelmenu,"s_1",&pscr1,&k1);
n_model=k1;
if(n_model == 0)
{
if( mess_y_n(56,4,"Quit session")) { saveent(menulevel); goto exi; }
}
else if(n_model == maxmodel+1)
{
clrbox(1,4,55,18);
makenewmodel();
menuhelp();
}
else if (n_model > 0)
{
put_text(&pscr1);
goto label_20;
}
}
label_20: /* Menu3:Enter Process */
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel = 2;
saveent(menulevel);
readModelFiles("./models",n_model);
modelinfo();
k2 = 1;
do
{ char strmen[]="\026"
" Enter Process "
" Force Unit.Gauge OFF "
" Edit model "
" Delete model ";
if(forceUG)improveStr(strmen,"OFF","ON");
menu1(56,4,"",strmen,"s_2_*",&pscr2,&k2);
switch (k2)
{
case 0: goto_xy(1,1); clr_eol(); goto label_10;
case 2: forceUG=!forceUG; modelinfo(); break;
case 3: editModel(1); break;
case 4:
if ( deletemodel(n_model) )
{
goto_xy(1,1);
clr_eol();
n_model=1;
fillModelMenu();
goto label_10;
}
else readModelFiles("./models",n_model);
}
} while (k2 != 1);
loadModel(0,forceUG);
label_21:
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel=2;
errorcode=enter(); /* Enter a process */
newCodes=0;
showheap();
if (errorcode) /* 'Esc' pressed */ { menuhelp(); goto label_20;}
errorcode=construct(); /* unSquared diagrams */
if (errorcode)
{ if(blind)
{ printf("Processes of this type are absent\n"); sortie(111); }
else
{ messanykey(5,22,"Processes of this type are absent");
clrbox(1,19,80,24);
goto label_21;
}
}
else if(!blind)
{ int dirStat=checkDir("results");
if(dirStat!=0)
{ messanykey( 10,10,"There are files in directory 'results/'.\n"
"To continue you has to clean or rename this directory.");
viewresults();
if(checkDir("results")!=0) goto label_21;
}
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
goto label_31;
}
label_30: /* Menu4: Squaring,...*/
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
label_31:
f3_key[0]=f3_key_prog;
f3_key[1]=NULL;
menulevel=3;
k3 = 1;
do
{
menu1(56,4,"","\026"
" View diagrams "
/* " Amplitude calculation" */
" Squaring technique "
" Write down processes ","s_squa_*",&pscr3,&k3);
switch (k3)
{
case 0: clrbox(1,2,55,11); menuhelp(); goto label_20;
case 1: viewfeyndiag(1); break;
case 3: { FILE*f=fopen("results/list_prc.txt","w");
int k,ndel,ncalc,nrest;
char process[100];
long recpos;
menup=fopen(MENUP_NAME,"r");
for(k=1;;k++)
{ int err=rd_menu(1,k,process,&ndel,&ncalc,&nrest,&recpos);
if(!err) break;
trim(process);
fprintf(f,"%s\n",process);
}
fclose(f);
fclose(menup);
messanykey(20,14,"See file 'results/list_prc.txt'");
} break;
/* case 2: messanykey(10,10,"Not implemented yet"); Amplitudes(); */
}
} while (k3 != 2);
if (!squaring()) goto label_30; /* process is absent */
clear_tmp();
saveent(menulevel);
restoreent(&menulevel);
label_40: /* Menu5: View squared diagrams..... */
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
sq_diagramsinfo(); /* ???????? */
k4=1;
saveent(menulevel);
pscr4=NULL;
for(;;)
{ int res;
menu1(56,4,"","\026"
" View squared diagrams"
" Symbolic calculations"
" Make&Launch n_calchep"
" Make n_calchep "
" REDUCE program "
,"s_calc_*",&pscr4,&k4);
res=checkDir("results");
if(res==1)
{
int n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id) unLockFile(n_calchep_id); else res=2;
}
switch (k4)
{ case 0:
if (mess_y_n(50,3,"Return to previous menu?"))goto label_30;
break;
case 1:
viewsqdiagr(); break;
case 2: /* Compute all diagrams */
restart2:
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
calcallproc();
sq_diagramsinfo();
if(!continuetest()) break;
showheap();
put_text(&pscr4);
break;
case 3:
{ static char keystr[12]="]{{[{";
if(res==2) messanykey(3,10,LOCKtxt); else inkeyString=keystr;
}
break;
case 4:
if(res==2) messanykey(3,10,LOCKtxt); else
{ char command[100];
FILE *f=fopen("results/EXTLIB","w");
fprintf(f,"EXTLIB=\"%s\"\nexport EXTLIB\n",EXTLIB);
fclose(f);
saveent(menulevel);
{ char command[100];
sprintf(command,"cd results; $CALCHEP/bin/make_VandP ../models %d",n_model);
system(command);
}
finish();
sortie(22);
}
case 5: mk_reduceprograms(); break;
}
if(k4==2 && continuetest()) { put_text(&pscr4); break; }
}
label_50:
k5=1;
pscr5=NULL;
menulevel=5;
saveent(menulevel);
for(;;)
{ int n_calchep_id;
menu1(56,4,"","\026"
" C code "
" C-compiler "
" Edit Linker "
" REDUCE code "
" MATHEMATICA code "
" FORM code "
" Enter new process " ,"s_out_*",&pscr5,&k5);
if((k5==1||k5==2)&&pid)
{ int epid=waitpid(pid,NULL,WNOHANG);
if(epid) pid=0; else
{
messanykey(3,10,"This option is frozen while n_calchep runs or is compiled");
continue;
}
}
switch (k5)
{ case 0: goto label_40; break;
case 1:
c_prog(); newCodes=0; saveent(menulevel);
break;
case 2:
if(newCodes) { c_prog(); newCodes=0; saveent(menulevel);}
pid=fork();
if(pid==0)
{
char * command;
n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id)
{
command=(char*)malloc(strlen(pathtocomphep)+strlen(EXTLIB)+100);
sprintf(command,"EXTLIB=\"%s\"\n"
"export EXTLIB\n"
"cd results\n"
"xterm -e %s/make__n_calchep %d",
EXTLIB,pathtocomphep,n_model);
system(command);
sprintf(command,"cd results; ./n_calchep");
unLockFile(n_calchep_id);
system(command);
free(command);
}
exit(0);
}
break;
case 3: if(edittable(1,1,&modelTab[4],1," ",0))
{ char fName[STRSIZ];
readEXTLIB();
sprintf(fName,"%smodels%c%s%d.mdl",pathtouser,f_slash,mdFls[4],n_model);
writetable( &modelTab[4],fName);
}
break;
case 4: makeReduceOutput(); break;
case 5: makeMathOutput(); break;
case 6: makeFormOutput(); break;
case 7:
put_text(&pscr5);
clrbox(1,2,55,11);
menuhelp();
goto label_20;
}
}
exi:
finish();
sortie(0);
return 0;
}
void editModel( int edit)
{ int n,i,j;
void * pscr = NULL;
int edited=0,renamed=0;
int ok=1;
char tabMenu[STRSIZ], tabName[80];
char menuName[30];
char * tabhelp[5]={ "s_mdl_1", "s_mdl_2", "s_mdl_3", "s_mdl_4","s_mdl_5"};
n = 1;
cont:
if(edit) f3_key[7]=NULL;
do
{
if(edit) strcpy(menuName,"Edit model"); else strcpy(menuName,"View model");
strcpy(tabMenu,"\017");
for (i=0;i<5;i++)
{ strcpy(tabName,modelTab[i].headln);
trim(tabName);
for(j=strlen(tabName);j<14;j++) tabName[j]=' ';
tabName[14]=0;
strcat(tabMenu," "); strcat(tabMenu,tabName);
}
if(edit) strcat(tabMenu," RENAME "
" CHECK MODEL ");
if(edit) menu1(56,7,"",tabMenu,"s_mdl_e",&pscr,&n);
else menu1(17,15,"",tabMenu,"s_mdl_v",&pscr,&n);
switch(n)
{ case 0:
if( edited || renamed )
{
if(mess_y_n(15,19," Save corrections ?") )
{
if (ok||loadModel(1,0) ) writeModelFiles(n_model); else goto cont;
if(renamed)
{ int i,size=modelmenu[0];
for(i=0;modelTab[0].mdlName[i]&&i<size-1;i++)
modelmenu[(n_model-1)*size +2+i]=modelTab[0].mdlName[i];
for(;i<size-1;i++) modelmenu[(n_model-1)*size +2+i]=' ';
modelinfo();
}
}else readModelFiles("./models",n_model);
}
f3_key[7]=f10_key_prog;
return;
case 1: case 2: case 3: case 4: case 5:
edited=edittable(1,1,&modelTab[n-1],1,tabhelp[n-1],!edit) || edited;
if(edited) ok=0;
break;
case 6:
if(correctStr (10,17,"New name:",modelTab[0].mdlName,29,1))
{ trim(modelTab[0].mdlName);
for(i=1;i<5;i++)strcpy(modelTab[i].mdlName,modelTab[0].mdlName);
renamed=1;
}
break;
case 7:
if (!loadModel(1,0) ) goto cont;
else {ok=1; messanykey(10,15,"The model is Ok");}
break;
}
}while(n);
}
int main(int argc,char** argv)
{
/*===================================
0 - First start.
1 - Model menu.
2 - Enter process menu.
3 - Feynman diagrams menu; squaring.
4 - Squared diagram menu; symbolic calculation.
5 - Write results; new process.
10 -Restart symbolic calculations
==========================================================*/
/* 0-Start; 1-Restart; 2-Heap Error,3-Edit Model,4-UserBreak */
void *pscr1=NULL,*pscr2=NULL,*pscr3=NULL,*pscr4=NULL,*pscr5=NULL;
int k1=1,k2=1,k3=1,k4=1,k5=1;
int n;
int pid=0;
char LOCKtxt[]="Directory 'results/' contains the .lock file created by n_calchep.\n"
"To continue you may a)Close the n_calchep session, or b)Rename 'results/',\n";
blind=0;
strcpy(pathtocalchep,argv[0]);
for(n=strlen(pathtocalchep)-1; n && pathtocalchep[n]!=f_slash; n--);
pathtocalchep[n-3]=0;
for ( n=1;n<argc;n++)
{ if (strcmp(argv[n],"-blind")==0 && n<argc-1 )
{ blind=1;
inkeyString=argv[++n];
}
if (strcmp(argv[n],"+blind")==0 ) blind=2;
if (strcmp(argv[n],"--version")==0 ) { printf("%s\n", VERSION_); exit(0);}
}
if(!writeLockFile(".lock"))
{ fprintf(stderr,"locked by other s_calchep. See .lock file\n");
exit(100);
}
strcpy(pathtouser,"");
sprintf(pathtohelp,"%shelp%c",pathtocalchep,f_slash);
outputDir="results/";
{ char * icon=(char *) malloc(strlen(pathtocalchep)+20);
char title[30];
int err;
sprintf(icon,"%s/include/icon",pathtocalchep);
sprintf(title,"CalcHEP_%s/symb", VERSION);
err=start1(title,icon,"calchep.ini",&xw_error);
if(err && blind==0)
{ printf("Error:You have launched interactive session for version compiled without X11 library.\n");
printf(" Presumably X11 development package is not installed in your computer.\n");
printf(" In this case directory /usr/include/X11/ is empty.\n");
printf("Options: a) Use blind session; b) Update Linux and recompile CalcHEP \n");
printf("Name of needed package\n");
printf(" libX11-devel Fedora/Scientific Linux/CYGWIN/Darwin(MAC)\n");
printf(" libX11-dev Ubuntu/Debian\n");
printf(" xorg-x11-devel SUSE\n");
exit(66);
}
free(icon);
}
fillModelMenu();
f3_key[0]=f3_key_prog; f3_mess[0]="Model";
f3_key[1]=f4_key_prog; f3_mess[1]="Diagrams";
f3_key[2]=f5_key_prog; f3_mess[2]="Switches";
f3_key[3]=f6_key_prog; f3_mess[3]="Results";
f3_mess[4]="Del";
f3_mess[5]="UnDel";
f3_key[6]=f9_key_prog; f3_mess[6]="Ref";
f3_key[7]=f10_key_prog; f3_mess[7]="Quit";
restoreent(&menulevel);
if(!blind && menulevel<2) cheplabel();
switch (menulevel)
{
case 10:
case 6:
case 5:
case 4:
case 3: readModelFiles("./models",n_model);
modelinfo();
loadModel(0,forceUG);
processinfo();
diagramsinfo();
k1=n_model;
break;
case 2: k1=n_model;
readModelFiles("./models",n_model);
break;
}
switch (menulevel)
{
case 2: menuhelp();
goto label_20;
case 3: goto label_31;
case 4: goto label_40;
case 5: goto label_50;
case 6:
case 10: goto restart2;
}
label_10: /* Menu2(ModelMenu): */
f3_key[0]=NULL; /*models*/
f3_key[1]=NULL; /*diagrams*/
menulevel = 1;
forceUG=0;
menuhelp();
for(;;)
{
showheap();
k1=n_model;
menu1(56,4,"",modelmenu,"s_1",&pscr1,&k1);
if(k1 == 0)
{
if( mess_y_n(56,4,"Quit session")) {n_model=0; saveent(menulevel); goto exi; }
}
else if(k1 == maxmodel+1)
{
clrbox(1,4,55,18);
makenewmodel();
menuhelp();
}
else if (k1 > 0)
{ int err=0;
put_text(&pscr1);
if(k1!=n_model || ldModelStatus==0) { err=readModelFiles("./models",k1);}
n_model=k1;
if(err){ if(blind) sortie(133); else goto label_10;} else goto label_20;
}
}
label_20: /* Menu3:Enter Process */
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel = 2;
saveent(menulevel);
modelinfo();
k2 = 1;
do
{ char strmen[]="\026"
" Enter Process "
" Force Unit.Gauge= OFF"
" Edit model "
" Numerical Evaluation "
"======================"
" Delete model ";
if(forceUG)improveStr(strmen,"OFF","ON");
menu1(56,4,"",strmen,"s_2_*",&pscr2,&k2);
switch (k2)
{
case 0: goto_xy(1,1); clr_eol(); goto label_10;
case 2: forceUG=!forceUG; modelinfo(); break;
case 3: editModel(1); break;
case 4: numcheck();
case 5: break;
case 6:
if(deletemodel(n_model))
{
goto_xy(1,1);
clr_eol();
n_model=1;
ldModelStatus=0;
fillModelMenu();
goto label_10;
}
}
} while (k2 != 1);
if(!loadModel(0,forceUG)) goto label_20;
label_21:
f3_key[0]=NULL;
f3_key[1]=NULL;
menulevel=2;
checkAuxDir(n_model);
errorcode=enter(); /* Enter a process */
newCodes=0;
showheap();
if (errorcode) /* 'Esc' pressed */ { menuhelp(); goto label_20;}
errorcode=construct(); /* unSquared diagrams */
if (errorcode)
{ if(blind)
{ printf("Processes of this type are absent\n"); sortie(111); }
else
{ messanykey(5,22,"Processes of this type are absent");
clrbox(1,19,80,24);
goto label_21;
}
}
else if(!blind)
{ int dirStat=checkDir("results");
if(dirStat!=0)
{ messanykey( 10,10,"There are files in directory 'results/'.\n"
"To continue you has to clean or rename this directory.");
viewresults();
if(checkDir("results")!=0) goto label_21;
}
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
goto label_31;
}
label_30: /* Menu4: Squaring,...*/
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
label_31:
f3_key[0]=f3_key_prog;
f3_key[1]=NULL;
menulevel=3;
k3 = 1;
do
{
menu1(56,4,"","\026"
" View diagrams "
/* " Amplitude calculation" */
" Square diagrams "
" Write down processes ","s_squa_*",&pscr3,&k3);
switch (k3)
{
case 0: clrbox(1,2,55,11); menuhelp(); goto label_20;
case 1: viewfeyndiag(1); break;
case 3: { FILE*f=fopen("results/list_prc.txt","w");
int k,ndel,ncalc,nrest;
char process[100];
long recpos;
menup=fopen(MENUP_NAME,"r");
for(k=1;;k++)
{ int err=rd_menu(1,k,process,&ndel,&ncalc,&nrest,&recpos);
if(!err) break;
trim(process);
fprintf(f,"%s\n",process);
}
fclose(f);
fclose(menup);
messanykey(20,14,"See file 'results/list_prc.txt'");
} break;
/* case 2: messanykey(10,10,"Not implemented yet"); Amplitudes(); */
}
} while (k3 != 2);
if (!squaring()) goto label_30; /* process is absent */
clear_tmp();
saveent(menulevel);
restoreent(&menulevel);
label_40: /* Menu5: View squared diagrams..... */
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
clr_scr(FGmain,BGmain);
modelinfo();
processinfo();
diagramsinfo();
sq_diagramsinfo(); /* ???????? */
k4=1;
saveent(menulevel);
pscr4=NULL;
for(;;)
{ int res;
menu1(56,4,"","\026"
" View squared diagrams"
" Symbolic calculations"
" Make&Launch n_calchep"
" Make n_calchep "
" REDUCE program "
,"s_calc_*",&pscr4,&k4);
res=checkDir("results");
if(res==1)
{
int n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id) unLockFile(n_calchep_id); else res=2;
}
switch (k4)
{ case 0:
if (mess_y_n(50,3,"Return to previous menu?"))goto label_30;
break;
case 1:
viewsqdiagr(); break;
case 2: /* Compute all diagrams */
restart2:
f3_key[0]=f3_key_prog;
f3_key[1]=f4_key_prog;
menulevel=4;
if(!nPROCSS ) calcallproc(); else
{ int *pids=malloc(sizeof(int)*nPROCSS);
int *pipes=malloc(2*sizeof(int)*nPROCSS);
int **qd=malloc(sizeof(int*)*nPROCSS);
int totD=sqDiagList(qd, nPROCSS);
int totC,nProc;
int k;
fflush(NULL);
for(k=0;k<nPROCSS;k++)
{ int* kpipe=pipes+2*k;
pipe(kpipe);
pids[k]=fork();
if(pids[k]==0)
{
close(kpipe[0]);
calcWithFork(k,qd[k],kpipe[1]);
exit(0);
}
}
for(k=0;k<nPROCSS;k++) close(pipes[2*k+1]);
infoLine(0.);
for(nProc=nPROCSS,totC=0;nProc;)
{ int one;
int err;
nProc=0;
for(k=0;k<nPROCSS;k++) if(pids[k])
{ if(waitpid(pids[k],NULL,WNOHANG)==0)
{ nProc++;
if(read(pipes[2*k],&one,sizeof(int)))totC+=one;
} else
{ for(;read(pipes[2*k],&one,sizeof(int))>0;) totC+=one;
pids[k]=0;
}
}
if(infoLine((double)totC/(double)totD)) for(k=0;k<nPROCSS;k++) if(pids[k]) kill(pids[k],SIGUSR1);
}
infoLine(2);
for(k=0;k<nPROCSS;k++)
{ char ctlgName[100];
char command[200];
sprintf(ctlgName,"%s_%d",CATALOG_NAME,k);
if(access(ctlgName,R_OK)==0)
{ sprintf(command," cat %s >> %s", ctlgName,CATALOG_NAME);
system(command);
unlink(ctlgName);
}
}
if(totC) newCodes=1;
updateMenuQ();
for(k=0;k<nPROCSS;k++) close(pipes[2*k]);
free(pids);
free(pipes);
for(k=0;k<nPROCSS;k++) free(qd[k]);
free(qd);
}
sq_diagramsinfo();
if(!continuetest()) break;
showheap();
put_text(&pscr4);
break;
case 3:
{ static char keystr[12]="]{{[{";
if(res==2) messanykey(3,10,LOCKtxt); else inkeyString=keystr;
}
break;
case 4:
if(res==2) messanykey(3,10,LOCKtxt); else
{ char command[100];
saveent(menulevel);
chdir("results");
makeVandP(0,"../models",n_model,2,pathtocalchep);
finish();
sortie(22);
}
case 5: mk_reduceprograms(); break;
}
if(k4==2 && continuetest()) { put_text(&pscr4); break; }
}
label_50:
k5=1;
pscr5=NULL;
menulevel=5;
saveent(menulevel);
sq_diagramsinfo();
for(;;)
{ int n_calchep_id;
menu1(56,4,"","\026"
" C code "
" C-compiler "
" Edit Linker "
" REDUCE code "
" MATHEMATICA code "
" FORM code "
" Enter new process " ,"s_out_*",&pscr5,&k5);
if((k5==1||k5==2)&&pid)
{ int epid=waitpid(pid,NULL,WNOHANG);
if(epid) pid=0; else
{
messanykey(3,10,"This option is frozen while n_calchep runs or is compiled");
continue;
}
}
switch (k5)
{ case 0: goto label_40; break;
case 1:
c_prog(); newCodes=0; saveent(menulevel);
break;
case 2:
if(newCodes) { c_prog(); newCodes=0; saveent(menulevel);}
n_calchep_id=setLockFile("results/.lock");
if(n_calchep_id)
{
if(nPROCSS)
{ chdir("results");
makeVandP(0,"../models",n_model,2,pathtocalchep);
pCompile();
if(access("./n_calchep",X_OK)==0)
{
fflush(NULL);
pid=fork();
if(pid==0)
{
system("./n_calchep");
exit(0);
}
} else messanykey(15,15,"n_calchep is not generated");
chdir("../");
} else
{ fflush(NULL);
pid=fork();
if(pid==0)
{ if(chdir("results")==0)
{ char*command=malloc(100+strlen(pathtocalchep));
makeVandP(0,"../models",n_model,2,pathtocalchep);
sprintf(command,"xterm -e %s/sbin/make__n_calchep %d", pathtocalchep,n_model);
system(command);
free(command);
system("./n_calchep");
}
exit(0);
}
}
unLockFile(n_calchep_id);
}
break;
case 3: if(edittable(1,1,&modelTab[4],1," ",0))
{ char fName[STRSIZ];
sprintf(fName,"%smodels%c%s%d.mdl",pathtouser,f_slash,mdFls[4],n_model);
writetable( &modelTab[4],fName);
}
break;
case 4: makeReduceOutput(); break;
case 5: makeMathOutput(); break;
case 6: makeFormOutput(); break;
case 7:
put_text(&pscr5);
clrbox(1,2,55,11);
menuhelp();
goto label_20;
}
}
exi:
finish();
sortie(0);
return 0;
}
int main(int argc,char** argv)
{
int n;
char icon[200];
char title[30];
blind=0;
for( n=1;n<argc;n++)
{
if(strcmp(argv[n],"--version")==0) { printf("%s\n",VERSION_); exit(0); }
if(strcmp(argv[n],"-blind")==0&& n<argc-1 )
{ blind=1;
inkeyString=argv[++n];
} else if (strcmp(argv[n],"+blind")==0 ) blind=2;
}
if(!writeLockFile(".lock"))
{ fprintf(stderr,"locked by other n_calchep. See .lock\n");
exit(100);
}
setenv("CALCHEP",rootDir,0);
sprintf(pathtocalchep,"%s%c",rootDir,f_slash);
sprintf(pathtohelp,"%shelp%c",pathtocalchep,f_slash);
sprintf(icon,"%s/include/icon",pathtocalchep);
sprintf(title,"CalcHEP_%s/num", VERSION);
f3_key[2]=f5_key_prog; f3_mess[2]="Options";
f3_key[3]=f6_key_prog; f3_mess[3]="Results";
f3_key[5]=f8_key_prog; f3_mess[5]="Calc";
f3_key[6]=f9_key_prog; f3_mess[6]="Ref";
f3_key[7]=f10_key_prog; f3_mess[7]="Quit";
{ int size=100;
for(;;)
{ compDir=realloc(compDir,size+20);
if(getcwd(compDir,size)) break; else size*=2;
}
strcat(compDir,"/aux");
libDir=malloc(strlen(compDir)+20);
sprintf(libDir,"%s/so_generated",compDir);
modelNum=1;
calchepDir=getenv("CALCHEP");
if(!calchepDir) calchepDir=interface_ext.CALCHEP;
ForceUG=interface_ext.forceUG;
}
/* ** initialization of the session */
link_process(PtrInterface_ext);
start1(title,icon,"calchep.ini;../calchep.ini",&xw_error);
nPROCSS=sysconf(_SC_NPROCESSORS_ONLN);
if(r_sess__(NULL)==-1)
{
char buff[200];
int pdg[11]={21,1,-1,2,-2,3,-3,4,-4,5,-5};
int k,ok=0;
strcpy(buff,"{p*\\09");
for(k=0;k<11;k++)
{ char*ch=pdg2name(pdg[k]);
if(ch)
{ if(ok) strcat(buff,",");
strcat(buff,ch);
ok=1;
}
}
strcat(buff,"}");
if(ok)
{ int blind0=blind;
char*inkeyString0=inkeyString;
inkeyString=buff;
blind=1;
edittable(1,4,&compTab,1,"n_comp",0);
blind=blind0;
inkeyString=inkeyString0;
fillCompositeArray();
strcpy(buff,"{%\\09T(p*)\\09 50{%\\09J(p*,p*)\\09 0.5}");
inkeyString=buff;
blind=1;
edittable(1,4,&cutTab,1,"n_cut",0);
blind=blind0;
inkeyString=inkeyString0;
}
}
{ char *ch=getenv("nParProc");
if(ch) sscanf(ch,"%d",&nPROCSS);
}
goto_xy(10,10); print("Calculation of constraints. Please, be patient.");
escpressed();
n_comphep();
finish();
sortie(0);
}
void editHist(void) {do edittable(1,4,&histTab,1,"n_distrib",0); while(correctHistList());}
int add_hist(FILE *f, char *procname)
{
table histTab2=histTab;
histRec * histPtr2=histPtr;
histRec *r,*r2;
char procname2[100];
int i,err=0;
histPtr=NULL;
histTab.strings=NULL;
linelist qLl;
histRec*histPtr_;
err=rdr_hist2(f,procname2);
if(err) return err;
if(!histPtr2){ strcpy(procname,procname2); return 0; }
if(!histPtr ){ histTab=histTab2; histPtr=histPtr2; return 0;}
if(strcmp(histTab.format,histTab2.format))
{ int l[6],l2[6];
char *ch,*ch2;
for(i=0,ch=histTab.format,ch2=histTab2.format; i<6; i++,ch++,ch2++)
{ ch =strchr(ch, '|'); l[i]=ch -histTab.format;
ch2=strchr(ch2,'|'); l2[i]=ch2-histTab2.format;
}
for(i=5;i;i--) {l[i]-=l[i-1]+1; l2[i]-=l2[i-1]+1;}
{ int blind_=blind;
char*inkeyString_=inkeyString;
char buff1[20]="",buff2[20];
blind=1;
for(i=0;i<6;i++,strcat(buff1,"\\09"))
{
if(l[i]<l2[i])
{ sprintf(buff2,"%s\\13%d{}",buff1,l2[i]);
inkeyString=buff2;
edittable(1,4,&histTab,1,"n_distrib",0);
}else if(l[i]>l2[i])
{ sprintf(buff2,"%s\\13%d{}",buff1,l[i]);
inkeyString=buff2;
edittable(1,4,&histTab2,1,"n_distrib",0);
}
}
blind=blind_;
inkeyString=inkeyString_;
correctHistList();
{
table histTab3=histTab;
histRec * histPtr3=histPtr;
histTab=histTab2;
histPtr=histPtr2;
correctHistList();
histTab2=histTab;
histPtr2=histPtr;
histTab=histTab3;
histPtr=histPtr3;
}
}
}
for(r2=histPtr2;r2;r2=r2->next) for(r=histPtr;r;r=r->next)
if(r->mother && strcmp2(r->mother->line,r2->mother->line)==0)
{ double n1=r->nPoints;
double n2=r2->nPoints;
double n=n1+n2;
int Ntot;
if(n1)
{ if(strcmp(r->key[1],"0")==0) Ntot=300; else Ntot=900;
if(n2) for(i=0;i<Ntot;i++)
{
r2->ff[i]=(n/n1)*(n/n1)*(r->ff[i] - r->f[i]*r->f[i]/n1 )
+(n/n2)*(n/n2)*(r2->ff[i]- r2->f[i]*r2->f[i]/n2);
r2->f[i]=(n/n1)*r->f[i]+(n/n2)*r2->f[i];
r2->ff[i]+=(r2->f[i])*(r2->f[i])/n;
}
else for(i=0;i<Ntot;i++)
{
r2->f[i]= r->f[i];
r2->ff[i]=r2->ff[i];
}
}
r2->nPoints=n;
if(r->mother->next) r->mother->next->pred = r->mother->pred;
if(r->mother->pred) r->mother->pred->next = r->mother->next;
else histTab.strings=r->mother->next;
free(r->mother);
r->mother=NULL;
break;
}
correctHistList();
qLl=histTab.strings;
if(qLl)
{
for( ;qLl->next; qLl=qLl->next) continue;
qLl->next=histTab2.strings;
histTab2.strings=histTab.strings;
for(histPtr_=histPtr; histPtr_->next; histPtr_=histPtr_->next) continue;
histPtr_->next=histPtr2;
histPtr2=histPtr;
}
histTab=histTab2;
histPtr=histPtr2;
return err;
}
int monte_carlo_menu(void)
{
static int r=0;
int mode=1;
void * pscr=NULL;
void * pscr_mem=NULL;
void (*quit)(int)=f3_key[7];
char menutxt[]="\030"
" Subprocess "
" IN state "
" Model parameters "
" Constraints "
" QCD coupling "
" Breit-Wigner "
" Aliases "
" Cuts "
" Phase space mapping "
" Monte Carlo simulation "
" Easy ";
if(nout_int!=2 ) menutxt[menutxt[0]*10+1]=0;
if(nin_int==1) improveStr(menutxt,"Easy", "Total width");
else improveStr(menutxt,"Easy", "1D intergration");
get_text(1,10,80,24,&pscr_mem);
wrtprc_();
for(;;)
{
infor();
f3_key[7]=quit;
menu1(54,4,"",menutxt,"n_mc_*",&pscr, &mode);
if(mode==1||mode==2||mode==3||mode==5||mode==7)f3_key[7]=NULL;
switch (mode)
{
case 0: put_text(&pscr_mem); return 0;
case 1: r=r|3*sub_men__(); break;
case 2: r=r|in_setting(); break;
case 3: r=r|change_parameter(54,7,0); break;
case 4: { int modeC=1;
for(;;)
{ char menuC[]="\030"
" All Constraints "
" Masses,Widths,Branching";
void * pscrC=NULL;
menu1(54,6,"",menuC,"n_constr_*",&pscrC, &modeC);
switch(modeC)
{ case 0: put_text(&pscr_mem); break;
case 1: show_depend(54,7); break;
case 2: show_spectrum(54,9); break;
}
if(!modeC) break;
} break;
}
case 5: r=r|qcdmen_(); break;
case 6: r=r|w_men__(); break;
case 7: do r=r|(3*edittable(1,4,&compTab,1,"n_comp",0)); while (fillCompositeArray());
break;
case 8: do r=r|(3*edittable(1,4,&cutTab,1,"n_cut",0)); while (fillCutArray());
break;
case 9: r=r|mappingMenu(); break;
case 10:
if(nout_int==1 && !sf_num[0] && !sf_num[1] )
{ if(blind) return 1;
messanykey(15,15,"Phase space integration for 2->1 processes\n needs distribution functions.");
break;
}
if(checkEnergy())
{ if(blind) return 1;
messanykey(15,15,"Energy is too small!");
break;
}
if(fillCutArray())
{ if(blind) return 2;
messanykey(15,15,"Can not evaluate cuts limlts");
break;
}
case 11:
if(mode==11)
{ void (*f10_tmp)(int);
w_sess__(NULL);
f10_tmp=f3_key[7];
f3_key[7]=f10_key_prog_for22;
if(nin_int==1) decay12(); else
{ REAL m1,m2, Pcm;
pinf_int(Nsub,1,&m1,NULL);
pinf_int(Nsub,2,&m2,NULL);
incomkin(m1,m2,inP1,inP2,NULL,&Pcm,NULL);
if(sf_num[0]||sf_num[1]||nCuts)
messanykey(10,10,"Structure functions and cuts are ignored\n");
cs_numcalc(Pcm);
}
f3_key[7]= f10_tmp;
r_sess__(NULL);
break;
} else if(fillRegArray())
{
if(blind) return 3;
messanykey(15,15,
"Can not evaluate regularization paremeters");
break;
}
if(mode==10) runVegas();
r=0;
break;
}
if(r) clearEventMax();
if(r&2) clearGrid();
if(r&1)newSession();
}
}
