-
Notifications
You must be signed in to change notification settings - Fork 0
/
sm.c
121 lines (96 loc) · 5.1 KB
/
sm.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
#include "src/include.h"
/*--------------------------------------------------------*/
/* Calculation of the observables in the Standard Model */
/*--------------------------------------------------------*/
int main()
{
struct parameters param;
Init_param(¶m);
slha_adjust(¶m);
param.SM=1;
printf("\n");
printf("SuperIso v3.4 - F. Mahmoudi\n\n");
printf("Standard Model predictions\n\n");
printf("Observable\t\t\tValue\n\n");
double C0b[11],C0spec[11],C1b[11],C1spec[11],C0w[11],C1w[11],C2w[11],C2b[11],Cpb[11];
double complex CQpb[3],CQ0b[3],CQ1b[3];
CQ0b[1]=CQ0b[2]=CQ1b[1]=CQ1b[2]=CQpb[1]=CQpb[2]=0.;
double obs[Nobs_BKsll+1];
double mu_W=2.*param.mass_W;
double mu_b=param.mass_b_pole/2.;
CW_calculator(C0w,C1w,C2w,mu_W,¶m);
C_calculator_base1(C0w,C1w,C2w,mu_W,C0b,C1b,C2b,mu_b,¶m);
Cprime_calculator(Cpb,CQpb,mu_W,mu_b,¶m);
printf("BR(b->s gamma)\t\t\t%.3e\n",bsgamma(C0b,C1b,C2b,Cpb,mu_b,mu_W,¶m));
double lambda_h=0.5;
double mu_spec=sqrt(lambda_h*mu_b);
C_calculator_base2(C0w,C1w,mu_W,C0b,C1b,mu_b,¶m);
C_calculator_base2(C0w,C1w,mu_W,C0spec,C1spec,mu_spec,¶m);
Cprime_calculator(Cpb,CQpb,mu_W,mu_b,¶m);
printf("delta0(B->K* gamma)\t\t%.3e\n\n",delta0(C0b,C0spec,C1b,C1spec,Cpb,¶m,mu_b,mu_spec));
mu_b=param.mass_b_pole;
C_calculator_base1(C0w,C1w,C2w,mu_W,C0b,C1b,C2b,mu_b,¶m);
Cprime_calculator(Cpb,CQpb,mu_W,mu_b,¶m);
printf("BR(Bs->mu mu)\t\t\t%.3e\n",Bsmumu(C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("BR(Bs->mu mu)_untag\t\t%.3e\n",Bsmumu_untag(C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("BR(Bd->mu mu)\t\t\t%.3e\n\n",Bdmumu(C0b,C1b,C2b,CQ0b,CQ1b,¶m,mu_b));
printf("BR(B->K* mu mu)_low\t\t%.3e\n",BRBKstarmumu(1.,6.,obs,C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("AFB(B->K* mu mu)_low\t\t%.3e\n",obs[1]);
printf("FL(B->K* mu mu)_low\t\t%.3e\n",obs[2]);
printf("P1=AT1(B->K* mu mu)_low\t\t%.3e\n",obs[4]);
printf("AT2(B->K* mu mu)_low\t\t%.3e\n",obs[5]);
printf("AT3(B->K* mu mu)_low\t\t%.3e\n",obs[6]);
printf("AT4(B->K* mu mu)_low\t\t%.3e\n",obs[7]);
printf("AT5(B->K* mu mu)_low\t\t%.3e\n",obs[8]);
printf("P4=HT1(B->K* mu mu)_low\t\t%.3e\n",obs[9]);
printf("P5=HT2(B->K* mu mu)_low\t\t%.3e\n",obs[10]);
printf("HT3(B->K* mu mu)_low\t\t%.3e\n",obs[11]);
printf("P2(B->K* mu mu)_low\t\t%.3e\n",obs[14]);
printf("P3(B->K* mu mu)_low\t\t%.3e\n",obs[15]);
printf("P6(B->K* mu mu)_low\t\t%.3e\n",obs[16]);
printf("P4'(B->K* mu mu)_low\t\t%.3e\n",obs[17]);
printf("P5'(B->K* mu mu)_low\t\t%.3e\n",obs[18]);
printf("P6'(B->K* mu mu)_low\t\t%.3e\n",obs[19]);
printf("P8(B->K* mu mu)_low\t\t%.3e\n",obs[20]);
printf("P8'(B->K* mu mu)_low\t\t%.3e\n",obs[21]);
printf("AI(B->K* mu mu)_low\t\t%.3e\n\n",AI_BKstarmumu_lowq2(C0b,C1b,C2b,¶m,mu_b));
printf("BR(B->K* mu mu)_high\t\t%.3e\n",BRBKstarmumu_highq2(obs,C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("AFB(B->K* mu mu)_high\t\t%.3e\n",obs[1]);
printf("FL(B->K* mu mu)_high\t\t%.3e\n",obs[2]);
printf("AT1(B->K* mu mu)_high\t\t%.3e\n",obs[4]);
printf("P1=AT2(B->K* mu mu)_high\t%.3e\n",obs[5]);
printf("AT3(B->K* mu mu)_high\t\t%.3e\n",obs[6]);
printf("AT4(B->K* mu mu)_high\t\t%.3e\n",obs[7]);
printf("AT5(B->K* mu mu)_high\t\t%.3e\n",obs[8]);
printf("P4=HT1(B->K* mu mu)_high\t%.3e\n",obs[9]);
printf("P5=HT2(B->K* mu mu)_high\t%.3e\n",obs[10]);
printf("HT3(B->K* mu mu)_high\t\t%.3e\n",obs[11]);
printf("P2(B->K* mu mu)_high\t\t%.3e\n",obs[14]);
printf("P3(B->K* mu mu)_high\t\t%.3e\n",obs[15]);
printf("P6(B->K* mu mu)_high\t\t%.3e\n",obs[16]);
printf("P4'(B->K* mu mu)_high\t\t%.3e\n",obs[17]);
printf("P5'(B->K* mu mu)_high\t\t%.3e\n",obs[18]);
printf("P6'(B->K* mu mu)_high\t\t%.3e\n",obs[19]);
printf("P8(B->K* mu mu)_high\t\t%.3e\n",obs[20]);
printf("P8'(B->K* mu mu)_high\t\t%.3e\n",obs[21]);
printf("AI(B->K* mu mu)_high\t\t%.3e\n\n",AI_BKstarmumu_highq2(C0b,C1b,C2b,¶m,mu_b));
double smin=pow(2.*param.mass_mu,2.);
double smax=pow(param.m_Bd-param.m_Kstar,2.)*0.999;
BRBKstarmumu(smin,smax,obs,C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b);
printf("q0^2(AFB(B->K* mu mu))\t\t%.3e\n",obs[0]);
printf("q0^2(AI(B->K* mu mu))\t\t%.3e\n\n",AI_BKstarmumu_zero(C0b,C1b,C2b,¶m,mu_b));
printf("BR(B->Xs mu mu)_low\t\t%.3e\n",BRBXsmumu_lowq2(C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("BR(B->Xs mu mu)_high\t\t%.3e\n",BRBXsmumu_highq2(C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("q0^2(AFB(B->Xs mu mu)\t\t%.3e\n",A_BXsmumu_zero(C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("BR(B->Xs tau tau)_high\t\t%.3e\n\n",BRBXstautau_highq2(C0b,C1b,C2b,CQ0b,CQ1b,Cpb,CQpb,¶m,mu_b));
printf("BR(B->tau nu)\t\t\t%.3e\n",Btaunu(¶m));
printf("R(B->tau nu)\t\t\t%.3e\n",RBtaunu(¶m));
printf("BR(B->D tau nu)\t\t\t%.3e\n",BDtaunu(¶m));
printf("BR(B->D tau nu)/BR(B->D e nu)\t%.3e\n",BDtaunu_BDenu(¶m));
printf("BR(Ds->tau nu)\t\t\t%.3e\n",Dstaunu(¶m));
printf("BR(Ds->mu nu)\t\t\t%.3e\n",Dsmunu(¶m));
printf("BR(D->mu nu)\t\t\t%.3e\n",Dmunu(¶m));
printf("BR(K->mu nu)/BR(pi->mu nu)\t%.3e\n",Kmunu_pimunu(¶m));
printf("Rmu23(K->mu nu)\t\t\t%.3e\n\n",Rmu23(¶m));
return 1;
}