REAL S34_out(double GG, REAL * momenta,int * err)
{REAL ans=0;
REAL DP[10];
REAL* V=va_out;
REAL mass[3],width[3];
char * Qtxt[3];
REAL Q0[3],Q1[3],Q2[3];
width[2]=V[7]; mass[2]=V[29]; Qtxt[2]="\3\4";
width[1]=V[13]; mass[1]=V[12]; Qtxt[1]="\1\2";
*err=*err|prepDen(2,nin_out,BWrange_out*
BWrange_out,mass,width,Qtxt,momenta,Q0,Q1,Q2);
sprod_(5, momenta, DP);
{int i; for(i=0;i<1;i++)
{ REAL r=Farr[i](GG,DP,Q0,Q1,Q2);
if(r>Fmax) Fmax=r;
ans+=r;
}}
return ans;
}
REAL S167_out(double GG, REAL * momenta,int * err)
{REAL ans=0;
REAL DP[15];
REAL* V=va_out;
REAL mass[4],width[4];
char * Qtxt[4];
REAL Q0[4],Q1[4],Q2[4];
width[2]=V[7]; mass[2]=V[53]; Qtxt[2]="\3\4";
width[1]=V[12]; mass[1]=V[11]; Qtxt[1]="\1\2\5";
width[3]=0.; mass[3]=V[53]; Qtxt[3]="\2\5";
*err=*err|prepDen(3,nin_out,BWrange_out*
BWrange_out,mass,width,Qtxt,momenta,Q0,Q1,Q2);
sprod_(6, momenta, DP);
{int i; for(i=0;i<1;i++)
{ REAL r=Farr[i](GG,DP,Q0,Q1,Q2);
if(r>Fmax) Fmax=r;
ans+=r;
}}
return ans;
}
Real S1_ext(Real * momenta,int * err)
{ Real ans=0;
Real s0max=0;
int i;
Real mass[3],Q[3];
Real* V=va_ext;
double width[3];
for(i=0; i<nin_ext; i++) s0max+=momenta[4*i];
width[2]=0.;
mass[2]=V[10];
Q[2]=mass[2]*mass[2]-sqrMom(nin_ext,"\1\3",momenta);
width[1]=0.;
mass[1]=V[10];
Q[1]=mass[1]*mass[1]-sqrMom(nin_ext,"\1\4",momenta);
//
// JH - Lambda = 500GeV
// double Lambda2 = 4.e+6; // Lambda^2
double Lambda2 = 2.5e+5; // Lambda^2
double ss = 4.*momenta[0]*momenta[4];
// JH - power 2.5
// double FormFac = pow(1.0+ss/Lambda2,2.);
double FormFac = pow(1.0+ss/Lambda2,2.5);
//
DkappA = V[4]/FormFac;
LambA = V[5]/FormFac;
a0w = V[6]/FormFac;
aCw = V[7]/FormFac;
//
printf("DkappA=%g, LambdA=%g, a0w=%g, aCw=%g\ns=%f\n",DkappA,LambA,a0w,aCw,ss);
//
*err=*err|prepDen(2, BWrange_ext*BWrange_ext, s0max,mass,width, Q);
for(i=0; i<6; i++)
{ Real r=Farr[i]();
if(r>Fmax) Fmax=r;
ans+=r;
}
return ans;
}
REAL S7_out(double GG, REAL * momenta,int * err)
{REAL ans=0;
REAL DP[15];
REAL* V=va_out;
REAL mass[8],width[8];
char * Qtxt[8];
REAL Q0[8],Q1[8],Q2[8];
width[7]=0.; mass[7]=V[9]; Qtxt[7]="\2\3";
width[6]=0.; mass[6]=V[8]; Qtxt[6]="\2\3";
width[3]=V[10]; mass[3]=V[9]; Qtxt[3]="\1\2";
width[5]=0.; mass[5]=V[8]; Qtxt[5]="\1\2";
width[2]=V[6]; mass[2]=V[23]; Qtxt[2]="\4\5";
width[1]=V[10]; mass[1]=V[9]; Qtxt[1]="\1\2\3";
width[4]=0.; mass[4]=0; Qtxt[4]="\1\3";
*err=*err|prepDen(7,nin_out,BWrange_out*
BWrange_out,mass,width,Qtxt,momenta,Q0,Q1,Q2);
sprod_(6, momenta, DP);
{int i; for(i=0;i<21;i++)
{ REAL r=Farr[i](GG,DP,Q0,Q1,Q2);
if(r>Fmax) Fmax=r;
ans+=r;
}}
return ans;
}