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);
}
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 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();
}
}
