int hbBlocksMDL(char*fname,int * nHch)
{
char buff[100];
sprintf(buff,"cat UMSSM_spectr.dat UMSSM_decay.dat > %s",fname);
system(buff);
FILE * f=fopen(fname,"a");
int pdg,i;
char *h[3]={"h1","h2","h3"};
if(!f) return 0;
fprintf(f,"Block HiggsBoundsInputHiggsCouplingsBosons\n");
fprintf(f,"# Effective coupling normalised to SM one and squared\n");
fprintf(f,"# For (*) normalized on Sin(2*W)\n");
for(i=0;i<3;i++)
{ pdg=pNum(h[i]);
fprintf(f," %12.4E 3 %d 24 24 # %s-W-W \n", SQR(slhaVal("REDCOUP",0.,2,1+i,4)),pdg,h[i]);
fprintf(f," %12.4E 3 %d 23 23 # %s-Z-Z \n", SQR(slhaVal("REDCOUP",0.,2,1+i,5)),pdg,h[i]);
fprintf(f," %12.4E 3 %d 21 21 # %s-gluon-gluon\n",SQR(slhaVal("REDCOUP",0.,2,1+i,6)),pdg,h[i]);
fprintf(f," %12.4E 3 %d 22 22 # %s-gamma-gamma\n",SQR(slhaVal("REDCOUP",0.,2,1+i,7)),pdg,h[i]);
}
{ pdg=pNum("ha");
fprintf(f," %12.4E 3 %d 24 24 # %s-W-W \n", SQR(slhaVal("REDCOUP",0.,2,4,4)),pdg,"ha");
fprintf(f," %12.4E 3 %d 23 23 # %s-Z-Z \n", SQR(slhaVal("REDCOUP",0.,2,4,5)),pdg,"ha");
fprintf(f," %12.4E 3 %d 21 21 # %s-gluon-gluon\n",SQR(slhaVal("REDCOUP",0.,2,4,6)),pdg,"ha");
fprintf(f," %12.4E 3 %d 22 22 # %s-gamma-gamma\n",SQR(slhaVal("REDCOUP",0.,2,4,7)),pdg,"ha");
}
fprintf(f,"Block HiggsBoundsInputHiggsCouplingsFermions\n");
fprintf(f,"# Effective coupling normalised to SM one and squared\n");
for(i=0;i<3;i++)
{ pdg=pNum(h[i]);
fprintf(f," %12.4E 0. 3 %d 5 5 # %s-b-b \n",SQR(slhaVal("REDCOUP",0.,2,1+i,3)),pdg,h[i]);
fprintf(f," %12.4E 0. 3 %d 6 6 # %s-top-top\n",SQR(slhaVal("REDCOUP",0.,2,1+i,1)),pdg,h[i]);
fprintf(f," %12.4E 0. 3 %d 15 15 # %s-tau-tau\n",SQR(slhaVal("REDCOUP",0.,2,1+i,2)),pdg,h[i]);
}
{ pdg=pNum("ha");
fprintf(f," 0. %12.4E 3 %d 5 5 # %s-b-b \n",SQR(slhaVal("REDCOUP",0.,2,4,3)),pdg,"ha");
fprintf(f," 0. %12.4E 3 %d 6 6 # %s-top-top\n",SQR(slhaVal("REDCOUP",0.,2,4,1)),pdg,"ha");
fprintf(f," 0. %12.4E 3 %d 15 15 # %s-tau-tau\n",SQR(slhaVal("REDCOUP",0.,2,4,2)),pdg,"ha");
}
fclose(f);
if(nHch) *nHch=1;
return 4;
}
static int channels(aChannel* Chann, int *i, double *w, int*pdg, char* txt, int len)
{ int k,j;
if(!Chann || *i<1) return 0;
for(k=0;k<*i;k++) if(Chann[k].weight==0) return 0;
k--;
*w=Chann[k].weight;
txt[0]=0;
for(j=0;j<2;j++) sprintf(txt+strlen(txt),"%s ",Chann[k].prtcl[j]);
strcpy(txt+strlen(txt),"-> ");
for(j=2;j<4;j++) sprintf(txt+strlen(txt),"%s ",Chann[k].prtcl[j]);
for(j=0;j<4;j++) pdg[j]=pNum(Chann[k].prtcl[j]);
if(Chann[k].prtcl[4])
{ sprintf(txt+strlen(txt),"%s ",Chann[k].prtcl[j]);
pdg[j]=pNum(Chann[k].prtcl[j]);
} else pdg[j]=0;
for(j=strlen(txt);j<len;j++) txt[j]=' ';
return 1;
}
void smodels(double Pcm, int nf,double csMinFb, char*fileName,int wrt)
{
int SMP[16]={1,2,3,4,5,6, 11,12,13,14,15,16, 21,22,23,24};
int i,j;
FILE*f=fopen(fileName,"w");
int np=0;
char**plist=NULL;
int smH=-1;
char* gluname=NULL;
char* phname=NULL;
char* bname=NULL;
char* Bname=NULL;
char* lname=NULL;
char* Lname=NULL;
// find SM Higgs
for(i=0;i<nModelParticles;i++)
{
if(ModelPrtcls[i].NPDG== 21) gluname=ModelPrtcls[i].name;
if(ModelPrtcls[i].NPDG== 22) phname=ModelPrtcls[i].name;
if(ModelPrtcls[i].NPDG== 5) { bname=ModelPrtcls[i].name; Bname=ModelPrtcls[i].aname;}
if(ModelPrtcls[i].NPDG== -5) { bname=ModelPrtcls[i].aname; Bname=ModelPrtcls[i].name; }
if(ModelPrtcls[i].NPDG== 15) { lname=ModelPrtcls[i].name; Lname=ModelPrtcls[i].aname;}
if(ModelPrtcls[i].NPDG==-15) { Lname=ModelPrtcls[i].aname; lname=ModelPrtcls[i].name; }
}
//printf("gluname %s bname %s lname %s\n", gluname,bname,lname);
if(gluname && bname && lname)
for(smH=0;smH<nModelParticles;smH++) if( ModelPrtcls[smH].spin2==0 && ModelPrtcls[smH].cdim==1
&& ModelPrtcls[smH].name[0]!='~' && strcmp(ModelPrtcls[smH].name,ModelPrtcls[smH].aname)==0 )
{ double w,ggBr,bbBr,llBr, hMass=pMass(ModelPrtcls[smH].name);
txtList L;
double ggBrSM=0.073, bbBrSM=0.60,llBrSM=0.063,wSM=4.24E-3;
double prec=0.9;
char chan[50];
if(hMass<123 || hMass>128) continue;
w=pWidth(ModelPrtcls[smH].name,&L);
sprintf(chan,"%s,%s",gluname,gluname);
ggBr=findBr(L, chan);
sprintf(chan,"%s,%s",lname,Lname);
llBr=findBr(L, chan);
sprintf(chan,"%s,%s",bname,Bname);
bbBr=findBr(L, chan);
if(ggBr==0) { bbBr*=w/(w+0.073*0.00424); llBr*=w/(w+ggBrSM*wSM);}
if( bbBrSM*prec< bbBr && bbBr<bbBrSM*(2-prec) && llBrSM*prec< llBr && llBr<llBrSM*(2-prec)) break;
}
if(smH<nModelParticles) printf("SM HIGGS=%s\n",ModelPrtcls[smH].name);
else printf("NO SM-like HIGGS in the model\n");
fprintf(f,"BLOCK MASS\n");
for(i=0;i<nModelParticles;i++) if(pMass(ModelPrtcls[i].name) <Pcm)
{
for(j=0;j<16;j++) if(abs(ModelPrtcls[i].NPDG)==SMP[j]) break;
if(j==16 )
{
np++;
plist=realloc(plist,np*sizeof(char*));
plist[np-1]=ModelPrtcls[i].name;
if(strcmp(ModelPrtcls[i].name,ModelPrtcls[i].aname))
{ np++;
plist=realloc(plist,np*sizeof(char*));
plist[np-1]=ModelPrtcls[i].aname;
}
fprintf(f," %d %E # %s \n",ModelPrtcls[i].NPDG,findValW(ModelPrtcls[i].mass),ModelPrtcls[i].name);
}
}
fprintf(f,"\n");
for(i=0;i<nModelParticles;i++)
{ for(j=0;j<16;j++) if(ModelPrtcls[i].NPDG==SMP[j]) break;
if(j==16) slhaDecayPrint(ModelPrtcls[i].name,1,f);
}
for(i=0;i<np;i++) for(j=i;j<np;j++) if(pMass(plist[i])+pMass(plist[j])<Pcm)
if(plist[i][0]=='~' && plist[j][0]=='~')
{ int q31,q32,q3,c1,c2;
qNumbers(plist[i], NULL, &q31,&c1);
qNumbers(plist[j], NULL, &q32,&c2);
q3=q31+q32;
if(q3<0) { q3*=-1; if(abs(c1)==3) c1*=-1; if(abs(c2)==3) c2*=-1;}
if(c1>c2){ int c=c1; c1=c2;c2=c;}
if ( (c2==1 || (c1==1 && c2==8) || (c1==-3 && c2==3) || (c1==8 && c2==8) )
&& (q3!=0 && q3 !=3) ) continue;
if ( ((c1==-3 && c2== 3)||(c1== 1 && c2== 1)||
(c1== 8 && c2== 8)||(c1== 1 && c2== 8)) && (q3!=0 && q3!=3) ) continue;
if ( ((c1== 3 && c2== 8)||(c1== 1 && c2== 3)) && (q3!=2) ) continue;
if ( ((c1==-3 && c2== 8)||(c1==-3 && c2== 1)) && (q3!=1) ) continue;
if ( (c1== 3 && c2== 3) && (q3!=4 && q3!=1) ) continue;
if ( (c1==-3 && c2==-3) && (q3!=2) ) continue;
{ double dcs;
double Qf=0.5*(pMass(plist[i])+pMass(plist[j]));
dcs=hCollider(Pcm,1,nf,Qf,Qf,plist[i],plist[j],0,wrt);
if(dcs>csMinFb*0.001)
{
fprintf(f,"XSECTION %E 2212 2212 2 %d %d\n",2*Pcm, pNum(plist[i]),pNum(plist[j]));
/*pb*/ fprintf(f,"0 0 0 0 0 0 %E micrOMEGAs 3.6\n\n", dcs);
}
}
}
fclose(f);
free(plist);
f=fopen("particles.py","w");
fprintf(f,"#!/usr/bin/env python\n");
fprintf(f,"rOdd ={\n");
for(np=0,i=0;i<nModelParticles;i++) if(ModelPrtcls[i].name[0]=='~' && pMass(ModelPrtcls[i].name) <Pcm )
{
if(np) fprintf(f,",\n");
fprintf(f, " %d : \"%s\",\n", ModelPrtcls[i].NPDG,ModelPrtcls[i].name);
fprintf(f, " %d : \"%s\"" , -ModelPrtcls[i].NPDG,ModelPrtcls[i].aname);
np++;
}
fprintf(f,"\n}\n");
fprintf(f,"rEven ={\n");
for(np=0,i=0;i<nModelParticles;i++) if(ModelPrtcls[i].name[0]!='~' && pMass(ModelPrtcls[i].name) <Pcm )
{
for(j=0;j<16;j++) if(abs(ModelPrtcls[i].NPDG)==SMP[j]) break;
if(j==16 )
{ if(np) fprintf(f,",\n");
if(ModelPrtcls[i].NPDG==smH)
{
fprintf(f, " %d : \"higgs\",\n", ModelPrtcls[i].NPDG);
fprintf(f, " %d : \"higgs\"\n", -ModelPrtcls[i].NPDG);
} else
{ char * n=ModelPrtcls[i].name;
char * an=ModelPrtcls[i].aname;
if(strcmp( n,"higgs")==0) n="!higgs";
if(strcmp(an,"higgs")==0) an="!higgs";
fprintf(f, " %d : \"%s\",\n", ModelPrtcls[i].NPDG,n);
fprintf(f, " %d : \"%s\"" , -ModelPrtcls[i].NPDG,an);
}
np++;
}
}
fprintf(f,",\n"
" 23 : \"Z\",\n"
" -23 : \"Z\",\n"
" 22 : \"photon\",\n"
" -22 : \"photon\",\n"
" 24 : \"W+\",\n"
" -24 : \"W-\",\n"
" 16 : \"nu\",\n"
" -16 : \"nu\",\n"
" 15 : \"ta-\",\n"
" -15 : \"ta+\",\n"
" 14 : \"nu\",\n"
" -14 : \"nu\",\n"
" 13 : \"mu-\",\n"
" -13 : \"mu+\",\n"
" 12 : \"nu\",\n"
" -12 : \"nu\",\n"
" 11 : \"e-\",\n"
" -11 : \"e+\",\n"
" 5 : \"b\",\n"
" -5 : \"b\",\n"
" 6 : \"t+\",\n"
" -6 : \"t-\",\n"
" 1 : \"jet\",\n"
" 2 : \"jet\",\n"
" 3 : \"jet\",\n"
" 4 : \"jet\",\n"
" 21 : \"jet\",\n"
" -21 : \"jet\",\n"
" -1 : \"jet\",\n"
" -2 : \"jet\",\n"
" -3 : \"jet\",\n"
" -4 : \"jet\"" );
fprintf(f,"\n}\n");
fprintf(f,
"\nptcDic = {\"e\" : [\"e+\", \"e-\"],\n"
" \"mu\" : [\"mu+\", \"mu-\"],\n"
" \"ta\" : [\"ta+\", \"ta-\"],\n"
" \"l+\" : [\"e+\", \"mu+\"],\n"
" \"l-\" : [\"e-\", \"mu-\"],\n"
" \"l\" : [\"e-\", \"mu-\", \"e+\", \"mu+\"],\n"
" \"W\" : [\"W+\", \"W-\"],\n"
" \"t\" : [\"t+\", \"t-\"],\n"
" \"L+\" : [\"e+\", \"mu+\", \"ta+\"],\n"
" \"L-\" : [\"e-\", \"mu-\", \"ta-\"],\n"
" \"L\" : [\"e+\", \"mu+\", \"ta+\", \"e-\", \"mu-\", \"ta-\"]}\n"
);
fprintf(f,"qNumbers ={\n");
for(np=0,i=0;i<nModelParticles;i++) if(pMass(ModelPrtcls[i].name) <Pcm )
{
for(j=0;j<16;j++) if(abs(ModelPrtcls[i].NPDG)==SMP[j]) break;
if(j==16 )
{ if(np) fprintf(f,",\n");
fprintf(f, " %d : [%d,%d,%d]", ModelPrtcls[i].NPDG, ModelPrtcls[i].spin2, ModelPrtcls[i].q3, ModelPrtcls[i].cdim);
np++;
}
}
fprintf(f,"\n}\n");
fclose(f);
}
int HBblocks(char * fname)
{ FILE * f=fopen(fname,"a");
double tb,sb,cb,alpha,sa,ca,ta,samb,camb,dMb,MbHl,MbSM,MbH,MbH3,Q;
if(!f) return 1;
Q=findValW("Q");
if(slhaDecayExists(pNum("h")) <0) slhaDecayPrint("h", f);
if(slhaDecayExists(pNum("H")) <0) slhaDecayPrint("H", f);
if(slhaDecayExists(pNum("H3"))<0) slhaDecayPrint("H3",f);
if(slhaDecayExists(pNum("t")) <0) slhaDecayPrint("t", f);
if(slhaDecayExists(pNum("H+"))<0) slhaDecayPrint("H+",f);
tb=findValW("tB");
sb=tb/sqrt(1+tb*tb);
cb=1/sqrt(1+tb*tb);
alpha=findValW("alpha");
sa=sin(alpha);
ca=cos(alpha);
ta=sa/ca;
samb=sa*cb-ca*sb;
camb=ca*cb+sa*sb;
dMb=findValW("dMb");
MbSM=findValW("Mb");
MbH= MbSM/(1+dMb)*(1+dMb*ta/tb);
MbH3=MbSM/(1+dMb)*(1-dMb/tb/tb);
MbHl=MbSM/(1+dMb)*(1-dMb/ta/tb);
fprintf(f,"Block HiggsBoundsInputHiggsCouplingsBosons\n");
fprintf(f,"# Effective coupling normalised to SM one and squared\n");
fprintf(f,"# For (*) normalized on Sin(2*W)\n");
fprintf(f," %12.4E 3 25 24 24 # higgs-W-W \n", SQR(samb) );
fprintf(f," %12.4E 3 25 23 23 # higgs-Z-Z \n", SQR(samb) );
fprintf(f," %12.4E 3 25 25 23 # higgs-higgs-Z \n", 0. );
{ assignVal("Q",pMass("h"));
calcMainFunc();
fprintf(f," %12.4E 3 25 21 21 # higgs-gluon-gluon\n", SQR(findValW("LGGh")/findValW("LGGSM")) );
fprintf(f," %12.4E 3 25 22 22 # higgs-gamma-gamma\n", SQR(findValW("LAAh")/findValW("LAASM")) );
}
fprintf(f," %12.4E 3 35 24 24 # higgs-W-W \n", SQR(camb) );
fprintf(f," %12.4E 3 35 23 23 # higgs-Z-Z \n", SQR(camb) );
fprintf(f," %12.4E 3 35 25 23 # higgs-higgs-Z \n", 0. );
fprintf(f," %12.4E 3 35 35 23 # higgs-higgs-Z \n", 0. );
{ assignVal("Q",pMass("H"));
calcMainFunc();
fprintf(f," %12.4E 3 35 21 21 # higgs-gluon-gluon\n",SQR(findValW("LGGH")/findValW("LGGSM")) );
fprintf(f," %12.4E 3 35 22 22 # higgs-gamma-gamma\n",SQR(findValW("LAAH")/findValW("LAASM")) );
}
fprintf(f," %12.4E 3 36 24 24 # higgs-W-W \n", 0. );
fprintf(f," %12.4E 3 36 23 23 # higgs-Z-Z \n", 0. );
{ assignVal("Q",pMass("H3"));
calcMainFunc();
fprintf(f," %12.4E 3 36 21 21 # higgs-gluon-gluon\n",SQR(findValW("LGGH3")/2/findValW("LGGSM")) );
fprintf(f," %12.4E 3 36 22 22 # higgs-gamma-gamma\n",SQR(findValW("LAAH3")/2/findValW("LAASM")) );
}
fprintf(f," %12.4E 3 36 25 23 #*higgs-higgs-Z \n", SQR(camb) );
fprintf(f," %12.4E 3 36 35 23 #*higgs-higgs-Z \n", SQR(samb) );
fprintf(f," %12.4E 3 36 36 23 #* higgs-higgs-Z \n", 0. );
fprintf(f,"Block HiggsBoundsInputHiggsCouplingsFermions\n");
fprintf(f,"# Effective coupling normalised to SM one and squared\n");
fprintf(f," %12.4E %12.4E 3 25 5 5 # higgs-b-b \n" ,SQR((sa/cb)*(MbHl/MbSM)),0.);
fprintf(f," %12.4E %12.4E 3 25 6 6 # higgs-top-top \n",SQR(ca/sb) ,0.);
fprintf(f," %12.4E %12.4E 3 25 15 15 # higgs-tau-tau \n",SQR(sa/cb) ,0.);
fprintf(f," %12.4E %12.4E 3 35 5 5 # higgs-b-b \n" ,SQR((ca/cb)*(MbH/MbSM)) ,0.);
fprintf(f," %12.4E %12.4E 3 35 6 6 # higgs-top-top \n",SQR(sa/sb) ,0.);
fprintf(f," %12.4E %12.4E 3 35 15 15 # higgs-tau-tau \n",SQR(ca/cb) ,0.);
fprintf(f," %12.4E %12.4E 3 36 5 5 # higgs-b-b \n" ,0.,SQR(tb*(MbH3/MbSM)));
fprintf(f," %12.4E %12.4E 3 36 6 6 # higgs-top-top \n",0.,SQR(1/tb) );
fprintf(f," %12.4E %12.4E 3 36 15 15 # higgs-tau-tau \n",0.,SQR(tb) );
assignValW("Q",Q);
calcMainFunc();
fclose(f);
return 0;
}
