58 double M_H_pow_2 =
pow(M_H,2.);
59 double M_W_pow_2 =
pow(M_W,2.);
60 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
66 double m_u[4],m_u_pow_2[4];
68 for (
int i = 0; i<3; ++i) {
69 for (
int j = 0; j<3; j++) {
78 for(
int i =0;i<3;++i) {
79 m_u_pow_2[i]=
pow(m_u[i],2.);
84 double D0h,D0h_c,D2h,D2h_c;
88 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
89 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
90 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
91 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
92 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
94 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
96 C1_chargedhiggs +=
pow(g_2,4.)*CKM_product*m_u_pow_2[i]*m_u_pow_2[j]*(2.*
pow(M_W,2.)*D0h*
pow(tbeta,-2.) - D2h_c*
pow(tbeta,-4.) - 2*D2h*
pow(tbeta,-2.))/(128.*
pow(
PI,2.)*
pow(M_W,4.));
103 double M_D[6],M_D_pow_2[6],dm[6];
106 double Mg_pow_2 =
pow(Mg,2);
116 for(
int i = 0; i<6; ++i) {
117 for(
int j = 0; j<6; ++j) {
121 for(
int i = 0; i<6; ++i) {
122 M_D_pow_2[i]=
pow(M_D[i],2);
131 for(
int i =0; i<6; ++i) {
132 for(
int j = 0; j<6; ++j) {
133 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
134 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
135 C1_gluino += -Z_D[2][i]*Z_D[2][j]*
pow(g_3,4.)*(D0g*
pow(Mg, 2) + 11.*D2g)*conj(Z_D[0][i])*conj(Z_D[0][j])/(144.*
pow(
PI, 2));
143 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
144 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
160 for(
int i=0; i<2; ++i){
161 M_ch_pow_2[i]=
pow(M_ch[i],2);
163 for(
int i=0; i<6; ++i){
164 M_U_pow_2[i]=
pow(M_U[i],2);
166 for(
int i=0; i <2; ++i) {
167 for(
int j=0; j<2; ++j) {
171 for(
int i=2; i<2; ++i) {
172 for(
int j=0; j<2; ++j) {
176 for(
int i=0; i<6; ++i) {
177 for(
int j=0; j<6; ++j) {
192 double common = sqrt(2.)/v2;
197 double otherc = sqrt(2.)/v1;
210 for(
int i = 0; i<6; ++i) {
211 for(
int j=0; j<6; ++j) {
212 for(
int a =0; a<2; ++a) {
213 for (
int b=0; b<2; ++b) {
214 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
215 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
216 for(
int k=0; k<3; ++k) {
219 C1_chargino += -D2ch*
pow(V_CKM[k][2], 2)*(-Q_e*Z_p[0][a]*conj(Z_U[k][i])*swi + Yu[k]*Z_p[1][a]*conj(Z_U[k+3][i]))*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*(-Q_e*Z_U[k][i]*conj(Z_p[0][b])*swi + Z_U[k+3][i]*conj(Yu[k])*conj(Z_p[1][b]))*(-Q_e*Z_U[k][j]*conj(Z_p[0][a])*swi + Z_U[k+3][j]*conj(Yu[k])*conj(Z_p[1][a]))*
pow(conj(V_CKM[k][0]), 2)/(32.0*
pow(
PI, 2));
231 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
236 double otherc = sqrt(2.)/v1;
243 M_ch0[0]=fabs(temp_ch0[0]);
244 M_ch0[1]=fabs(temp_ch0[1]);
245 M_ch0[2]=fabs(temp_ch0[2]);
246 M_ch0[3]=fabs(temp_ch0[3]);
248 for(
int i=0; i<4; ++i){
249 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
253 for(
int i=0; i<6; ++i) {
254 M_D_pow_2[i]=
pow(M_D[i],2);
259 for(
int i=0; i<4; ++i) {
260 for(
int j=0; j<4; ++j) {
266 for(
int i=0; i<4; ++i){
268 for(
int j=0; j<4; ++j) {
279 for(
int i=0; i<6; ++i) {
280 for(
int j=0; j<6; ++j) {
281 for(
int a=0; a<4; ++a) {
282 for (
int b=0; b<4; ++b) {
283 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
284 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
285 C1_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*
286 (-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*
287 (-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*
288 (-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(64.0*
pow(
PI, 2)) -
289 D2ne*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*
290 (-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*
291 (-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*
292 (-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
308 double M_g_pow_2 =
pow(M_g,2.);
316 for(
int i=0; i<6; ++i) {
317 for(
int j=0; j<6; ++j) {
318 for(
int a=0; a<4; ++a) {
319 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
320 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
322 C1_mixed += -D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[0][i]*Z_D[0][j]*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[2][i])/(2.0*cw*sw) +
323 conj(Yd[2])*conj(Z_D[5][i])*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[2][j])/(2.0*cw*sw) + conj(Yd[2])*conj(Z_D[5][j])*conj(Z_N[2][a])) +
324 Z_D[2][i]*Z_D[2][j]*
pow(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]), 2))/(96.0*
pow(
PI, 2)) -
325 D2mix*Z_D[2][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) +
326 conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0][i])/(8.0*
pow(
PI, 2));
336 return C1_chargedhiggs+C1_gluino+C1_chargino+C1_neutralino;
397 double M_H_pow_2 =
pow(M_H,2.);
398 double M_W_pow_2 =
pow(M_W,2.);
399 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
405 double m_u[4],m_u_pow_2[4];
407 for (
int i = 0; i<3; ++i) {
408 for (
int j = 0; j<3; j++) {
417 for(
int i =0;i<3;++i) {
418 m_u_pow_2[i]=
pow(m_u[i],2.);
424 double D0h,D0h_c,D2h,D2h_c;
428 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
429 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
430 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
431 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
432 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
434 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
436 Cp1_chargedhiggs += -
pow(g_2,4.)*
pow(m_b,2.)*
pow(m_q,2.)*CKM_product*(D2h_c*
pow(tbeta,4.)+2.*D2h*
pow(tbeta,2.))/(128.*
pow(
PI,2.)*
pow(M_W,4.));
442 double M_D[6],M_D_pow_2[6],dm[6];
445 double Mg_pow_2 =
pow(Mg,2);
455 for(
int i = 0; i<6; ++i) {
456 for(
int j = 0; j<6; ++j) {
460 for(
int i = 0; i<6; ++i) {
461 M_D_pow_2[i]=
pow(M_D[i],2);
470 for(
int i =0; i<6; ++i) {
471 for(
int j = 0; j<6; ++j) {
472 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
473 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
475 Cp1_gluino += -Z_D[5][i]*Z_D[5][j]*
pow(g_3,4.)*(D0g*
pow(Mg, 2.) + 11.*D2g)*conj(Z_D[0+3][i])*conj(Z_D[0+3][j])/(144.*
pow(
PI, 2.));
484 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
485 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
501 for(
int i=0; i<2; ++i){
502 M_ch_pow_2[i]=
pow(M_ch[i],2);
504 for(
int i=0; i<6; ++i){
505 M_U_pow_2[i]=
pow(M_U[i],2);
507 for(
int i=0; i <2; ++i) {
508 for(
int j=0; j<2; ++j) {
512 for(
int i=2; i<2; ++i) {
513 for(
int j=0; j<2; ++j) {
517 for(
int i=0; i<6; ++i) {
518 for(
int j=0; j<6; ++j) {
533 double common = sqrt(2.)/v2;
538 double otherc = sqrt(2.)/v1;
551 for(
int i = 0; i<6; ++i) {
552 for(
int j=0; j<6; ++j) {
553 for(
int a =0; a<2; ++a) {
554 for (
int b=0; b<2; ++b) {
555 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
556 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
557 for(
int k=0; k<3; ++k) {
559 Cp1_chargino += -D2ch*
pow(V_CKM[k][2], 2)*
pow(Yd[2], 2)*Z_m[1][a]*Z_m[1][b]*Z_U[k][i]*Z_U[k][j]*
pow(conj(V_CKM[k][0]), 2)*
pow(conj(Yd[0]), 2)*conj(Z_m[1][a])*conj(Z_m[1][b])*conj(Z_U[k][i])*conj(Z_U[k][j])/(32.0*
pow(
PI, 2));
572 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
577 double otherc = sqrt(2.)/v1;
584 M_ch0[0]=fabs(temp_ch0[0]);
585 M_ch0[1]=fabs(temp_ch0[1]);
586 M_ch0[2]=fabs(temp_ch0[2]);
587 M_ch0[3]=fabs(temp_ch0[3]);
589 for(
int i=0; i<4; ++i){
590 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
594 for(
int i=0; i<6; ++i) {
595 M_D_pow_2[i]=
pow(M_D[i],2);
600 for(
int i=0; i<4; ++i) {
601 for(
int j=0; j<4; ++j) {
607 for(
int i=0; i<4; ++i){
609 for(
int j=0; j<4; ++j) {
621 for(
int i=0; i<6; ++i) {
622 for(
int j=0; j<6; ++j) {
623 for(
int a=0; a<4; ++a) {
624 for (
int b=0; b<4; ++b) {
625 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
626 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
628 Cp1_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][b])/(3.0*cw) +
629 Yd[2]*Z_D[2][j]*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][i])/(3.0*cw) +
630 Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][i]))/(64.0*
pow(
PI, 2)) - D2ne*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][b])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
631 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
632 Yd[2]*Z_D[5][i]*Z_N[2][a])/(32.0*
pow(
PI, 2));
649 double M_g_pow_2 =
pow(M_g,2.);
657 for(
int i=0; i<6; ++i) {
658 for(
int j=0; j<6; ++j) {
659 for(
int a=0; a<4; ++a) {
660 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
661 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
663 Cp1_mixed += -D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[0+3][i]*Z_D[0+3][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][j])/(3.0*cw) +
664 Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][j])) + Z_D[5][i]*Z_D[5][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
665 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j])))/(96.0*
pow(
PI, 2)) - D2mix*Z_D[5][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
666 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*conj(Z_D[0+3][i])/(8.0*
pow(
PI, 2));
735 double M_H_pow_2 =
pow(M_H,2.);
736 double M_W_pow_2 =
pow(M_W,2.);
737 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
743 double m_u[4],m_u_pow_2[4];
745 for (
int i = 0; i<3; ++i) {
746 for (
int j = 0; j<3; j++) {
755 for(
int i =0;i<3;++i) {
756 m_u_pow_2[i]=
pow(m_u[i],2.);
762 double D0h,D0h_c,D2h,D2h_c;
766 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
767 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
768 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
769 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
770 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
772 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
774 C2_chargedhiggs += -
pow(g_2,4.)*
pow(m_q,2.)*CKM_product*
pow(m_u[i],2)*
pow(m_u[j],2)*(D0h_c - 2*D0h)/(128.*
pow(
PI,2.)*
pow(M_W,4.));
781 double M_D[6],M_D_pow_2[6],dm[6];
784 double Mg_pow_2 =
pow(Mg,2);
794 for(
int i = 0; i<6; ++i) {
795 for(
int j = 0; j<6; ++j) {
799 for(
int i = 0; i<6; ++i) {
800 M_D_pow_2[i]=
pow(M_D[i],2);
809 for(
int i =0; i<6; ++i) {
810 for(
int j = 0; j<6; ++j) {
811 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
812 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
813 C2_gluino += -17.*D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[2][j]*
pow(g_3, 4.)*conj(Z_D[0+3][i])*conj(Z_D[0+3][j])/(288.*
pow(
PI, 2));
822 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
823 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
839 for(
int i=0; i<2; ++i){
840 M_ch_pow_2[i]=
pow(M_ch[i],2);
842 for(
int i=0; i<6; ++i){
843 M_U_pow_2[i]=
pow(M_U[i],2);
845 for(
int i=0; i <2; ++i) {
846 for(
int j=0; j<2; ++j) {
850 for(
int i=2; i<2; ++i) {
851 for(
int j=0; j<2; ++j) {
855 for(
int i=0; i<6; ++i) {
856 for(
int j=0; j<6; ++j) {
871 double common = sqrt(2.)/v2;
876 double otherc = sqrt(2.)/v1;
890 for(
int i = 0; i<6; ++i) {
891 for(
int j=0; j<6; ++j) {
892 for(
int a =0; a<2; ++a) {
893 for (
int b=0; b<2; ++b) {
894 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
895 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
896 for(
int k=0; k<3; ++k) {
911 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
916 double otherc = sqrt(2.)/v1;
923 M_ch0[0]=fabs(temp_ch0[0]);
924 M_ch0[1]=fabs(temp_ch0[1]);
925 M_ch0[2]=fabs(temp_ch0[2]);
926 M_ch0[3]=fabs(temp_ch0[3]);
928 for(
int i=0; i<4; ++i){
929 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
933 for(
int i=0; i<6; ++i) {
934 M_D_pow_2[i]=
pow(M_D[i],2);
939 for(
int i=0; i<4; ++i) {
940 for(
int j=0; j<4; ++j) {
946 for(
int i=0; i<4; ++i){
948 for(
int j=0; j<4; ++j) {
960 for(
int i=0; i<6; ++i) {
961 for(
int j=0; j<6; ++j) {
962 for(
int a=0; a<4; ++a) {
963 for (
int b=0; b<4; ++b) {
964 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
965 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
967 C2_neutralino += D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][j])/(3.0*cw) +
968 Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/
969 (2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])/(32.0*
pow(
PI, 2));
985 double M_g_pow_2 =
pow(M_g,2.);
993 for(
int i=0; i<6; ++i) {
994 for(
int j=0; j<6; ++j) {
995 for(
int a=0; a<4; ++a) {
996 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
997 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
1000 C2_mixed += D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[2][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
1001 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*conj(Z_D[2][j]) + 3.0*Z_D[2][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 -
1002 Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[0+3][i]))/(48.0*
pow(
PI, 2)) + D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
1003 Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[0+3][i])*conj(Z_D[0+3][j])/(48.0*
pow(
PI, 2));
1070 double M_H_pow_2 =
pow(M_H,2.);
1071 double M_W_pow_2 =
pow(M_W,2.);
1072 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
1078 double m_u[4],m_u_pow_2[4];
1080 for (
int i = 0; i<3; ++i) {
1081 for (
int j = 0; j<3; j++) {
1090 for(
int i =0;i<3;++i) {
1091 m_u_pow_2[i]=
pow(m_u[i],2.);
1097 double D0h,D0h_c,D2h,D2h_c;
1101 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
1102 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
1103 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
1104 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
1105 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
1107 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
1109 Cp2_chargedhiggs += -
pow(g_2,4.)*
pow(m_b,2.)*CKM_product*
pow(m_u[i],2)*
pow(m_u[j],2)*(D0h_c-2.*D0h)/(128.*
pow(
PI,2.)*
pow(M_W,4.));
1115 double M_D[6],M_D_pow_2[6],dm[6];
1118 double Mg_pow_2 =
pow(Mg,2);
1128 for(
int i = 0; i<6; ++i) {
1129 for(
int j = 0; j<6; ++j) {
1133 for(
int i = 0; i<6; ++i) {
1134 M_D_pow_2[i]=
pow(M_D[i],2);
1143 for(
int i =0; i<6; ++i) {
1144 for(
int j = 0; j<6; ++j) {
1145 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
1146 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
1148 Cp2_gluino += -17.*D0g*
pow(Mg, 2.)*Z_D[5][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[0][i])*conj(Z_D[0][j])/(288.*
pow(
PI, 2.));
1156 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
1157 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
1173 for(
int i=0; i<2; ++i){
1174 M_ch_pow_2[i]=
pow(M_ch[i],2);
1176 for(
int i=0; i<6; ++i){
1177 M_U_pow_2[i]=
pow(M_U[i],2);
1179 for(
int i=0; i <2; ++i) {
1180 for(
int j=0; j<2; ++j) {
1184 for(
int i=2; i<2; ++i) {
1185 for(
int j=0; j<2; ++j) {
1189 for(
int i=0; i<6; ++i) {
1190 for(
int j=0; j<6; ++j) {
1205 double common = sqrt(2.)/v2;
1210 double otherc = sqrt(2.)/v1;
1224 for(
int i = 0; i<6; ++i) {
1225 for(
int j=0; j<6; ++j) {
1226 for(
int a =0; a<2; ++a) {
1227 for (
int b=0; b<2; ++b) {
1228 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
1229 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
1230 for(
int k=0; k<3; ++k) {
1246 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
1251 double otherc = sqrt(2.)/v1;
1258 M_ch0[0]=fabs(temp_ch0[0]);
1259 M_ch0[1]=fabs(temp_ch0[1]);
1260 M_ch0[2]=fabs(temp_ch0[2]);
1261 M_ch0[3]=fabs(temp_ch0[3]);
1263 for(
int i=0; i<4; ++i){
1264 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
1268 for(
int i=0; i<6; ++i) {
1269 M_D_pow_2[i]=
pow(M_D[i],2);
1274 for(
int i=0; i<4; ++i) {
1275 for(
int j=0; j<4; ++j) {
1281 for(
int i=0; i<4; ++i){
1282 if(temp_ch0[i]<0.) {
1283 for(
int j=0; j<4; ++j) {
1301 for(
int i=0; i<6; ++i) {
1302 for(
int j=0; j<6; ++j) {
1303 for(
int a=0; a<4; ++a) {
1304 for (
int b=0; b<4; ++b) {
1305 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
1306 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
1307 C1_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*
1308 (-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*
1309 (-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*
1310 (-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(64.0*
pow(
PI, 2)) -
1311 D2ne*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*
1312 (-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*
1313 (-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*
1314 (-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
1317 Cp2_neutralino += D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
1318 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 -
1319 conj(Z_N[1][b])*cw)*conj(Z_D[0][j])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][j])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
1335 double M_g_pow_2 =
pow(M_g,2.);
1343 for(
int i=0; i<6; ++i) {
1344 for(
int j=0; j<6; ++j) {
1345 for(
int a=0; a<4; ++a) {
1346 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
1347 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
1349 Cp2_mixed += D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[5][i]*
pow(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) +
1350 conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]), 2)*conj(Z_D[5][j]) + 3.0*Z_D[5][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
1351 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0][i]))/(48.0*
pow(
PI, 2)) +
1352 D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
1353 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*conj(Z_D[0][i])*conj(Z_D[0][j])/(48.0*
pow(
PI, 2));
1421 double M_H_pow_2 =
pow(M_H,2.);
1422 double M_W_pow_2 =
pow(M_W,2.);
1423 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
1429 double m_u[4],m_u_pow_2[4];
1431 for (
int i = 0; i<3; ++i) {
1432 for (
int j = 0; j<3; j++) {
1441 for(
int i =0;i<3;++i) {
1442 m_u_pow_2[i]=
pow(m_u[i],2.);
1449 double D0h,D0h_c,D2h,D2h_c;
1453 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
1454 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
1455 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
1456 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
1457 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
1459 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
1461 C3_chargedhiggs += 0.;
1468 double M_D[6],M_D_pow_2[6],dm[6];
1471 double Mg_pow_2 =
pow(Mg,2);
1481 for(
int i = 0; i<6; ++i) {
1482 for(
int j = 0; j<6; ++j) {
1486 for(
int i = 0; i<6; ++i) {
1487 M_D_pow_2[i]=
pow(M_D[i],2);
1496 for(
int i =0; i<6; ++i) {
1497 for(
int j = 0; j<6; ++j) {
1498 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
1499 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
1501 C3_gluino += -D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[2][j]*
pow(g_3,4.)*conj(Z_D[0+3][i])*conj(Z_D[0+3][j])/(96.*
pow(
PI, 2));
1510 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
1511 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
1527 for(
int i=0; i<2; ++i){
1528 M_ch_pow_2[i]=
pow(M_ch[i],2);
1530 for(
int i=0; i<6; ++i){
1531 M_U_pow_2[i]=
pow(M_U[i],2);
1533 for(
int i=0; i <2; ++i) {
1534 for(
int j=0; j<2; ++j) {
1538 for(
int i=2; i<2; ++i) {
1539 for(
int j=0; j<2; ++j) {
1543 for(
int i=0; i<6; ++i) {
1544 for(
int j=0; j<6; ++j) {
1559 double common = sqrt(2.)/v2;
1564 double otherc = sqrt(2.)/v1;
1577 for(
int i = 0; i<6; ++i) {
1578 for(
int j=0; j<6; ++j) {
1579 for(
int a =0; a<2; ++a) {
1580 for (
int b=0; b<2; ++b) {
1581 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
1582 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
1583 for(
int k=0; k<3; ++k) {
1585 C3_chargino += -D0ch*M_ch[a]*M_ch[b]*
pow(V_CKM[k][2], 2)*Z_m[1][a]*Z_m[1][b]*Z_U[k][i]*Z_U[k][j]*(-Q_e*Z_p[0][a]*conj(Z_U[k][i])*swi + Yu[k]*Z_p[1][a]*conj(Z_U[k+3][i]))*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*
pow(conj(V_CKM[k][0]), 2)*
pow(conj(Yd[0]), 2)/(32.0*
pow(
PI, 2));
1598 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
1603 double otherc = sqrt(2.)/v1;
1610 M_ch0[0]=fabs(temp_ch0[0]);
1611 M_ch0[1]=fabs(temp_ch0[1]);
1612 M_ch0[2]=fabs(temp_ch0[2]);
1613 M_ch0[3]=fabs(temp_ch0[3]);
1615 for(
int i=0; i<4; ++i){
1616 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
1620 for(
int i=0; i<6; ++i) {
1621 M_D_pow_2[i]=
pow(M_D[i],2);
1626 for(
int i=0; i<4; ++i) {
1627 for(
int j=0; j<4; ++j) {
1633 for(
int i=0; i<4; ++i){
1634 if(temp_ch0[i]<0.) {
1635 for(
int j=0; j<4; ++j) {
1647 for(
int i=0; i<6; ++i) {
1648 for(
int j=0; j<6; ++j) {
1649 for(
int a=0; a<4; ++a) {
1650 for (
int b=0; b<4; ++b) {
1651 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
1652 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
1654 C3_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*((-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) +
1655 Yd[2]*Z_D[5][j]*Z_N[2][b]) - (-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
1656 Yd[2]*Z_D[5][j]*Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
1657 Yd[2]*Z_D[5][i]*Z_N[2][a])/(32.0*
pow(
PI, 2));
1674 double M_g_pow_2 =
pow(M_g,2.);
1682 for(
int i=0; i<6; ++i) {
1683 for(
int j=0; j<6; ++j) {
1684 for(
int a=0; a<4; ++a) {
1685 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
1686 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
1688 C3_mixed += D0mix*M_ch0[a]*M_g*Z_D[2][i]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
1689 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*conj(Z_D[2][j])/(48.0*
pow(
PI, 2)) - D0mix*M_ch0[a]*M_g*Z_D[2][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
1690 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[0+3][i])/(48.0*
pow(
PI, 2)) +
1691 D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
1692 Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[0+3][i])*conj(Z_D[0+3][j])/(48.0*
pow(
PI, 2));
1759 double M_H_pow_2 =
pow(M_H,2.);
1760 double M_W_pow_2 =
pow(M_W,2.);
1761 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
1767 double m_u[4],m_u_pow_2[4];
1769 for (
int i = 0; i<3; ++i) {
1770 for (
int j = 0; j<3; j++) {
1779 for(
int i =0;i<3;++i) {
1780 m_u_pow_2[i]=
pow(m_u[i],2.);
1786 double D0h,D0h_c,D2h,D2h_c;
1790 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
1791 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
1792 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
1793 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
1794 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
1796 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
1798 Cp3_chargedhiggs += 0.;
1804 double M_D[6],M_D_pow_2[6],dm[6];
1807 double Mg_pow_2 =
pow(Mg,2);
1817 for(
int i = 0; i<6; ++i) {
1818 for(
int j = 0; j<6; ++j) {
1822 for(
int i = 0; i<6; ++i) {
1823 M_D_pow_2[i]=
pow(M_D[i],2);
1832 for(
int i =0; i<6; ++i) {
1833 for(
int j = 0; j<6; ++j) {
1834 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
1835 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
1837 Cp3_gluino += -D0g*
pow(Mg, 2)*Z_D[5][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[0][i])*conj(Z_D[0][j])/(96.*
pow(
PI, 2));
1845 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
1846 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
1862 for(
int i=0; i<2; ++i){
1863 M_ch_pow_2[i]=
pow(M_ch[i],2);
1865 for(
int i=0; i<6; ++i){
1866 M_U_pow_2[i]=
pow(M_U[i],2);
1868 for(
int i=0; i <2; ++i) {
1869 for(
int j=0; j<2; ++j) {
1873 for(
int i=2; i<2; ++i) {
1874 for(
int j=0; j<2; ++j) {
1878 for(
int i=0; i<6; ++i) {
1879 for(
int j=0; j<6; ++j) {
1894 double common = sqrt(2.)/v2;
1899 double otherc = sqrt(2.)/v1;
1911 for(
int i = 0; i<6; ++i) {
1912 for(
int j=0; j<6; ++j) {
1913 for(
int a =0; a<2; ++a) {
1914 for (
int b=0; b<2; ++b) {
1915 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
1916 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
1917 for(
int k=0; k<3; ++k) {
1919 Cp3_chargino += -D0ch*M_ch[a]*M_ch[b]*
pow(V_CKM[k][2], 2)*
pow(Yd[2], 2)*(-Q_e*Z_U[k][i]*conj(Z_p[0][b])*swi + Z_U[k+3][i]*conj(Yu[k])*conj(Z_p[1][b]))*(-Q_e*Z_U[k][j]*conj(Z_p[0][a])*swi + Z_U[k+3][j]*conj(Yu[k])*conj(Z_p[1][a]))*
pow(conj(V_CKM[k][0]), 2)*conj(Z_m[1][a])*conj(Z_m[1][b])*conj(Z_U[k][i])*conj(Z_U[k][j])/(32.0*
pow(
PI, 2));
1931 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
1936 double otherc = sqrt(2.)/v1;
1943 M_ch0[0]=fabs(temp_ch0[0]);
1944 M_ch0[1]=fabs(temp_ch0[1]);
1945 M_ch0[2]=fabs(temp_ch0[2]);
1946 M_ch0[3]=fabs(temp_ch0[3]);
1948 for(
int i=0; i<4; ++i){
1949 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
1953 for(
int i=0; i<6; ++i) {
1954 M_D_pow_2[i]=
pow(M_D[i],2);
1959 for(
int i=0; i<4; ++i) {
1960 for(
int j=0; j<4; ++j) {
1966 for(
int i=0; i<4; ++i){
1967 if(temp_ch0[i]<0.) {
1968 for(
int j=0; j<4; ++j) {
1979 for(
int i=0; i<6; ++i) {
1980 for(
int j=0; j<6; ++j) {
1981 for(
int a=0; a<4; ++a) {
1982 for (
int b=0; b<4; ++b) {
1983 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
1984 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
1986 Cp3_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][j])/(2.0*cw*sw) +
1987 conj(Yd[0])*conj(Z_D[0+3][j])*conj(Z_N[2][b])) + (-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][b])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
1988 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a])))*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
1989 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
2005 double M_g_pow_2 =
pow(M_g,2.);
2013 for(
int i=0; i<6; ++i) {
2014 for(
int j=0; j<6; ++j) {
2015 for(
int a=0; a<4; ++a) {
2016 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
2017 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
2019 Cp3_mixed += D0mix*M_ch0[a]*M_g*Z_D[5][i]*
pow(g_3, 2)*
pow(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) +
2020 conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]), 2)*conj(Z_D[5][j])/(48.0*
pow(
PI, 2)) - D0mix*M_ch0[a]*M_g*Z_D[5][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
2021 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) +
2022 conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0][i])/(48.0*
pow(
PI, 2)) + D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
2023 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*conj(Z_D[0][i])*conj(Z_D[0][j])/(48.0*
pow(
PI, 2));
2089 double M_H_pow_2 =
pow(M_H,2.);
2090 double M_W_pow_2 =
pow(M_W,2.);
2091 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
2097 double m_u[4],m_u_pow_2[4];
2099 for (
int i = 0; i<3; ++i) {
2100 for (
int j = 0; j<3; j++) {
2109 for(
int i =0;i<3;++i) {
2110 m_u_pow_2[i]=
pow(m_u[i],2.);
2116 double D0h,D0h_c,D2h,D2h_c;
2120 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
2121 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
2122 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
2123 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
2124 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
2126 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
2128 C4_chargedhiggs +=
pow(g_2,4.)*CKM_product*(m_b*m_q*D2h*
pow(tbeta,2.)/
pow(M_W,2.) - m_b*m_q*
pow(m_u[i],2)*
pow(m_u[j],2)*(D0h_c + D0h*(
pow(tbeta,2.) +
pow(tbeta,-2.)))/(4*
pow(M_W,4.)))/(16.*
pow(
PI,2.));
2135 double M_D[6],M_D_pow_2[6],dm[6];
2138 double Mg_pow_2 =
pow(Mg,2);
2148 for(
int i = 0; i<6; ++i) {
2149 for(
int j = 0; j<6; ++j) {
2153 for(
int i = 0; i<6; ++i) {
2154 M_D_pow_2[i]=
pow(M_D[i],2);
2163 for(
int i =0; i<6; ++i) {
2164 for(
int j = 0; j<6; ++j) {
2165 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
2166 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
2168 C4_gluino += -7.*D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[0][i])*conj(Z_D[0+3][j])/(48.*
pow(
PI, 2)) + D2g*
pow(g_3,4.)*(6.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[0][i])*conj(Z_D[0+3][j]) + 11.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[0][j])*conj(Z_D[0+3][i]))/(72.*
pow(
PI, 2));
2177 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
2178 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
2194 for(
int i=0; i<2; ++i){
2195 M_ch_pow_2[i]=
pow(M_ch[i],2);
2197 for(
int i=0; i<6; ++i){
2198 M_U_pow_2[i]=
pow(M_U[i],2);
2200 for(
int i=0; i <2; ++i) {
2201 for(
int j=0; j<2; ++j) {
2205 for(
int i=2; i<2; ++i) {
2206 for(
int j=0; j<2; ++j) {
2210 for(
int i=0; i<6; ++i) {
2211 for(
int j=0; j<6; ++j) {
2226 double common = sqrt(2.)/v2;
2231 double otherc = sqrt(2.)/v1;
2243 for(
int i = 0; i<6; ++i) {
2244 for(
int j=0; j<6; ++j) {
2245 for(
int a =0; a<2; ++a) {
2246 for (
int b=0; b<2; ++b) {
2247 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
2248 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
2249 for(
int k=0; k<3; ++k) {
2251 C4_chargino += D2ch*
pow(V_CKM[k][2], 2)*Yd[2]*Z_m[1][a]*Z_U[k][j]*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*(-Q_e*Z_U[k][i]*conj(Z_p[0][b])*swi + Z_U[k+3][i]*conj(Yu[k])*conj(Z_p[1][b]))*
pow(conj(V_CKM[k][0]), 2)*conj(Yd[0])*conj(Z_m[1][a])*conj(Z_U[k][i])/(8.0*
pow(
PI, 2));
2264 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
2269 double otherc = sqrt(2.)/v1;
2276 M_ch0[0]=fabs(temp_ch0[0]);
2277 M_ch0[1]=fabs(temp_ch0[1]);
2278 M_ch0[2]=fabs(temp_ch0[2]);
2279 M_ch0[3]=fabs(temp_ch0[3]);
2281 for(
int i=0; i<4; ++i){
2282 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
2286 for(
int i=0; i<6; ++i) {
2287 M_D_pow_2[i]=
pow(M_D[i],2);
2292 for(
int i=0; i<4; ++i) {
2293 for(
int j=0; j<4; ++j) {
2299 for(
int i=0; i<4; ++i){
2300 if(temp_ch0[i]<0.) {
2301 for(
int j=0; j<4; ++j) {
2313 for(
int i=0; i<6; ++i) {
2314 for(
int j=0; j<6; ++j) {
2315 for(
int a=0; a<4; ++a) {
2316 for (
int b=0; b<4; ++b) {
2317 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
2318 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
2320 C4_neutralino += D2ne*((-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) +
2321 Yd[2]*Z_D[5][j]*Z_N[2][b]) + (-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
2322 Yd[2]*Z_D[5][j]*Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) +
2323 conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(8.0*
pow(
PI, 2));
2340 double M_g_pow_2 =
pow(M_g,2.);
2348 for(
int i=0; i<6; ++i) {
2349 for(
int j=0; j<6; ++j) {
2350 for(
int a=0; a<4; ++a) {
2351 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
2352 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
2354 C4_mixed += D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2355 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0+3][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
2356 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[0][i]))/(16.0*
pow(
PI, 2)) -
2357 D2mix*
pow(g_3, 2)*(-3.0*Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2358 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a])) - 3.0*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][i])/(3.0*cw) +
2359 Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][i]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[0][j])*conj(Z_D[0+3][i]))/(24.0*
pow(
PI, 2)) -
2360 D2mix*
pow(g_3, 2)*(-Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2361 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a])) - (-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
2362 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[0][j])*conj(Z_D[0+3][i]))/(24.0*
pow(
PI, 2)) -
2363 D2mix*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
2364 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*conj(Z_D[0][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
2365 Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0+3][i]))/(24.0*
pow(
PI, 2));
2377 scalar_t delta_d[6][6],delta_d_LL[3][3],delta_d_LR[3][3],delta_d_RL[3][3],delta_d_RR[3][3];
2378 scalar_t delta_u[6][6],delta_u_LL[3][3],delta_u_LR[3][3],delta_u_RL[3][3],delta_u_RR[3][3];
2380 double m_av = (M_U[0]+M_U[1]+M_U[2]+M_U[3]+M_U[4]+M_U[5]+M_D[0]+M_D[1]+M_D[2]+M_D[3]+M_D[4]+M_D[5])/12.;
2382 getDelta(delta_d,Z_D,M_D,m_av,delta_d_LL,delta_d_LR,delta_d_RL,delta_d_RR);
2383 getDelta(delta_u,Z_U,M_U,m_av,delta_u_LL,delta_u_LR,delta_u_RL,delta_u_RR);
2389 double x_mu =
pow(abs(mu),2)/
pow(m_av,2);
2391 double x_g =
pow(M_g,2)/
pow(m_av,2);
2392 double alpha_s =
pow(g_3,2.)/(4.*
PI) ;
2393 double alpha_2 =
pow(g_2,2.)/(4.*
PI);
2394 double eps=2*alpha_s*mu*M_g*
f(x_g)/(3.*
PI*
pow(m_av,2));
2400 scalar_t a2 = -alpha_s*
pow(M_g,2)*delta_d_LL[2][0]*delta_d_RR[2][0]*
pow(
h1(x_g),2);
2401 scalar_t a3 = alpha_2*
pow(m_t,2)*A_t*M_g*
h1(x_g)*
h3(x_mu)*delta_d_RR[2][0]*V_tb*conj(V_tq)/(
pow(M_W,2));
2402 scalar_t a4 = alpha_2*M_2*M_g*delta_u_LL[2][0]*delta_d_RR[2][0]*
h1(x_g)*
h4(x_2,x_g);
2404 scalar_t C4_higgspenguin = a1*(a2+a3+a4);
2465 double M_H_pow_2 =
pow(M_H,2.);
2466 double M_W_pow_2 =
pow(M_W,2.);
2467 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
2473 double m_u[4],m_u_pow_2[4];
2475 for (
int i = 0; i<3; ++i) {
2476 for (
int j = 0; j<3; j++) {
2485 for(
int i =0;i<3;++i) {
2486 m_u_pow_2[i]=
pow(m_u[i],2.);
2493 double D0h,D0h_c,D2h,D2h_c;
2497 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
2498 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
2499 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
2500 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
2501 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
2503 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][0])*conj(V_CKM[j][0]);
2505 C5_chargedhiggs +=
pow(g_2,4.)*m_b*m_q*CKM_product*
pow(m_u[i],2.)*(D2h_c-2.*D2h)/(32.*
pow(
PI,2.)*
pow(M_W,4.));
2512 double M_D[6],M_D_pow_2[6],dm[6];
2515 double Mg_pow_2 =
pow(Mg,2);
2525 for(
int i = 0; i<6; ++i) {
2526 for(
int j = 0; j<6; ++j) {
2530 for(
int i = 0; i<6; ++i) {
2531 M_D_pow_2[i]=
pow(M_D[i],2);
2540 for(
int i =0; i<6; ++i) {
2541 for(
int j = 0; j<6; ++j) {
2542 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
2543 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
2545 C5_gluino += -D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[0][i])*conj(Z_D[0+3][j])/(144.*
pow(
PI, 2)) + 5.*D2g*
pow(g_3,4.)*(-2.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[0][i])*conj(Z_D[0+3][j]) + 3.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[0][j])*conj(Z_D[0+3][i]))/(72.*
pow(
PI, 2));
2554 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
2555 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
2571 for(
int i=0; i<2; ++i){
2572 M_ch_pow_2[i]=
pow(M_ch[i],2);
2574 for(
int i=0; i<6; ++i){
2575 M_U_pow_2[i]=
pow(M_U[i],2);
2577 for(
int i=0; i <2; ++i) {
2578 for(
int j=0; j<2; ++j) {
2582 for(
int i=2; i<2; ++i) {
2583 for(
int j=0; j<2; ++j) {
2587 for(
int i=0; i<6; ++i) {
2588 for(
int j=0; j<6; ++j) {
2603 double common = sqrt(2.)/v2;
2608 double otherc = sqrt(2.)/v1;
2622 for(
int i = 0; i<6; ++i) {
2623 for(
int j=0; j<6; ++j) {
2624 for(
int a =0; a<2; ++a) {
2625 for (
int b=0; b<2; ++b) {
2626 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
2627 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
2628 for(
int k=0; k<3; ++k) {
2630 C5_chargino += -D0ch*M_ch[a]*M_ch[b]*
pow(V_CKM[k][2], 2)*Yd[2]*Z_m[1][b]*Z_U[k][i]*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*(-Q_e*Z_U[k][j]*conj(Z_p[0][a])*swi + Z_U[k+3][j]*conj(Yu[k])*conj(Z_p[1][a]))*
pow(conj(V_CKM[k][0]), 2)*conj(Yd[0])*conj(Z_m[1][a])*conj(Z_U[k][i])/(16.0*
pow(
PI, 2));
2643 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
2648 double otherc = sqrt(2.)/v1;
2655 M_ch0[0]=fabs(temp_ch0[0]);
2656 M_ch0[1]=fabs(temp_ch0[1]);
2657 M_ch0[2]=fabs(temp_ch0[2]);
2658 M_ch0[3]=fabs(temp_ch0[3]);
2660 for(
int i=0; i<4; ++i){
2661 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
2665 for(
int i=0; i<6; ++i) {
2666 M_D_pow_2[i]=
pow(M_D[i],2);
2671 for(
int i=0; i<4; ++i) {
2672 for(
int j=0; j<4; ++j) {
2678 for(
int i=0; i<4; ++i){
2679 if(temp_ch0[i]<0.) {
2680 for(
int j=0; j<4; ++j) {
2692 for(
int i=0; i<6; ++i) {
2693 for(
int j=0; j<6; ++j) {
2694 for(
int a=0; a<4; ++a) {
2695 for (
int b=0; b<4; ++b) {
2696 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
2697 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
2699 C5_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][i])/(3.0*cw) +
2700 Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2701 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))/(16.0*
pow(
PI, 2)) - D2ne*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
2702 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][b]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0- Z_N[1][a]*cw)/(2.0*cw*sw) +
2703 Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][b]))/(8.0*
pow(
PI, 2));
2720 double M_g_pow_2 =
pow(M_g,2.);
2728 for(
int i=0; i<6; ++i) {
2729 for(
int j=0; j<6; ++j) {
2730 for(
int a=0; a<4; ++a) {
2731 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
2732 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
2734 C5_mixed += -D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2735 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0+3][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
2736 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[0][i]))/(48.0*
pow(
PI, 2)) +
2737 D2mix*
pow(g_3, 2)*(-Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][i]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2738 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a])) - (-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][i])/(3*cw) +
2739 Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][i]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[0][j])*conj(Z_D[0+3][i]))/(24.0*
pow(
PI, 2)) +
2740 D2mix*
pow(g_3, 2)*(-Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2741 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a])) - (-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
2742 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[0][j])*conj(Z_D[0+3][i]))/(8.0*
pow(
PI, 2)) +
2743 D2mix*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[0+3][j])/(3.0*cw) +
2744 Z_N[2][a]*conj(Yd[0])*conj(Z_D[0][j]))*conj(Z_D[0][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
2745 conj(Z_N[1][a])*cw)*conj(Z_D[0][i])/(2.0*cw*sw) + conj(Yd[0])*conj(Z_D[0+3][i])*conj(Z_N[2][a]))*conj(Z_D[0+3][i]))/(8.0*
pow(
PI, 2));
2813 double M_H_pow_2 =
pow(M_H,2.);
2814 double M_W_pow_2 =
pow(M_W,2.);
2815 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
2821 double m_u[4],m_u_pow_2[4];
2823 for (
int i = 0; i<3; ++i) {
2824 for (
int j = 0; j<3; j++) {
2833 for(
int i =0;i<3;++i) {
2834 m_u_pow_2[i]=
pow(m_u[i],2.);
2839 double D0h,D0h_c,D2h,D2h_c;
2843 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
2844 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
2845 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
2846 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
2847 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
2849 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
2851 C1_chargedhiggs +=
pow(g_2,4.)*CKM_product*m_u_pow_2[i]*m_u_pow_2[j]*(2.*
pow(M_W,2.)*D0h*
pow(tbeta,-2.) - D2h_c*
pow(tbeta,-4.) - 2*D2h*
pow(tbeta,-2.))/(128.*
pow(
PI,2.)*
pow(M_W,4.));
2857 double M_D[6],M_D_pow_2[6],dm[6];
2860 double Mg_pow_2 =
pow(Mg,2);
2870 for(
int i = 0; i<6; ++i) {
2871 for(
int j = 0; j<6; ++j) {
2875 for(
int i = 0; i<6; ++i) {
2876 M_D_pow_2[i]=
pow(M_D[i],2);
2885 for(
int i =0; i<6; ++i) {
2886 for(
int j = 0; j<6; ++j) {
2887 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
2888 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
2889 C1_gluino += -Z_D[2][i]*Z_D[2][j]*
pow(g_3,4.)*(D0g*
pow(Mg, 2) + 11.*D2g)*conj(Z_D[1][i])*conj(Z_D[1][j])/(144.*
pow(
PI, 2));
2898 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
2899 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
2915 for(
int i=0; i<2; ++i){
2916 M_ch_pow_2[i]=
pow(M_ch[i],2);
2918 for(
int i=0; i<6; ++i){
2919 M_U_pow_2[i]=
pow(M_U[i],2);
2921 for(
int i=0; i <2; ++i) {
2922 for(
int j=0; j<2; ++j) {
2926 for(
int i=2; i<2; ++i) {
2927 for(
int j=0; j<2; ++j) {
2931 for(
int i=0; i<6; ++i) {
2932 for(
int j=0; j<6; ++j) {
2947 double common = sqrt(2.)/v2;
2952 double otherc = sqrt(2.)/v1;
2966 for(
int i = 0; i<6; ++i) {
2967 for(
int j=0; j<6; ++j) {
2968 for(
int a =0; a<2; ++a) {
2969 for (
int b=0; b<2; ++b) {
2970 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
2971 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
2972 for(
int k=0; k<3; ++k) {
2975 C1_chargino += -D2ch*
pow(V_CKM[k][2], 2)*(-Q_e*Z_p[0][a]*conj(Z_U[k][i])*swi + Yu[k]*Z_p[1][a]*conj(Z_U[k+3][i]))*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*(-Q_e*Z_U[k][i]*conj(Z_p[0][b])*swi + Z_U[k+3][i]*conj(Yu[k])*conj(Z_p[1][b]))*(-Q_e*Z_U[k][j]*conj(Z_p[0][a])*swi + Z_U[k+3][j]*conj(Yu[k])*conj(Z_p[1][a]))*
pow(conj(V_CKM[k][1]), 2)/(32.0*
pow(
PI, 2));
2988 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
2993 double otherc = sqrt(2.)/v1;
3000 M_ch0[0]=fabs(temp_ch0[0]);
3001 M_ch0[1]=fabs(temp_ch0[1]);
3002 M_ch0[2]=fabs(temp_ch0[2]);
3003 M_ch0[3]=fabs(temp_ch0[3]);
3005 for(
int i=0; i<4; ++i){
3006 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
3010 for(
int i=0; i<6; ++i) {
3011 M_D_pow_2[i]=
pow(M_D[i],2);
3016 for(
int i=0; i<4; ++i) {
3017 for(
int j=0; j<4; ++j) {
3023 for(
int i=0; i<4; ++i){
3024 if(temp_ch0[i]<0.) {
3025 for(
int j=0; j<4; ++j) {
3037 for(
int i=0; i<6; ++i) {
3038 for(
int j=0; j<6; ++j) {
3039 for(
int a=0; a<4; ++a) {
3040 for (
int b=0; b<4; ++b) {
3041 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
3042 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
3043 C1_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*
3044 (-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*
3045 (-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*
3046 (-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][b]))/(64.0*
pow(
PI, 2)) -
3047 D2ne*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*
3048 (-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*
3049 (-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*
3050 (-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
3067 double M_g_pow_2 =
pow(M_g,2.);
3075 for(
int i=0; i<6; ++i) {
3076 for(
int j=0; j<6; ++j) {
3077 for(
int a=0; a<4; ++a) {
3078 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
3079 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
3081 C1_mixed += -D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[1][i]*Z_D[1][j]*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[2][i])/(2.0*cw*sw) +
3082 conj(Yd[2])*conj(Z_D[5][i])*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[2][j])/(2.0*cw*sw) + conj(Yd[2])*conj(Z_D[5][j])*conj(Z_N[2][a])) +
3083 Z_D[2][i]*Z_D[2][j]*
pow(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]), 2))/(96.0*
pow(
PI, 2)) -
3084 D2mix*Z_D[2][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) +
3085 conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1][i])/(8.0*
pow(
PI, 2));
3152 double M_H_pow_2 =
pow(M_H,2.);
3153 double M_W_pow_2 =
pow(M_W,2.);
3154 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
3160 double m_u[4],m_u_pow_2[4];
3162 for (
int i = 0; i<3; ++i) {
3163 for (
int j = 0; j<3; j++) {
3172 for(
int i =0;i<3;++i) {
3173 m_u_pow_2[i]=
pow(m_u[i],2.);
3180 double D0h,D0h_c,D2h,D2h_c;
3184 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
3185 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
3186 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
3187 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
3188 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
3190 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
3192 Cp1_chargedhiggs += -
pow(g_2,4.)*
pow(m_b,2.)*
pow(m_q,2.)*CKM_product*(D2h_c*
pow(tbeta,4.)+2.*D2h*
pow(tbeta,2.))/(128.*
pow(
PI,2.)*
pow(M_W,4.));
3199 double M_D[6],M_D_pow_2[6],dm[6];
3202 double Mg_pow_2 =
pow(Mg,2);
3212 for(
int i = 0; i<6; ++i) {
3213 for(
int j = 0; j<6; ++j) {
3217 for(
int i = 0; i<6; ++i) {
3218 M_D_pow_2[i]=
pow(M_D[i],2);
3227 for(
int i =0; i<6; ++i) {
3228 for(
int j = 0; j<6; ++j) {
3229 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
3230 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
3231 Cp1_gluino += -Z_D[5][i]*Z_D[5][j]*
pow(g_3,4.)*(D0g*
pow(Mg, 2.) + 11.*D2g)*conj(Z_D[1+3][i])*conj(Z_D[1+3][j])/(144.*
pow(
PI, 2.));
3239 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
3240 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
3256 for(
int i=0; i<2; ++i){
3257 M_ch_pow_2[i]=
pow(M_ch[i],2);
3259 for(
int i=0; i<6; ++i){
3260 M_U_pow_2[i]=
pow(M_U[i],2);
3262 for(
int i=0; i <2; ++i) {
3263 for(
int j=0; j<2; ++j) {
3267 for(
int i=2; i<2; ++i) {
3268 for(
int j=0; j<2; ++j) {
3272 for(
int i=0; i<6; ++i) {
3273 for(
int j=0; j<6; ++j) {
3288 double common = sqrt(2.)/v2;
3293 double otherc = sqrt(2.)/v1;
3307 for(
int i = 0; i<6; ++i) {
3308 for(
int j=0; j<6; ++j) {
3309 for(
int a =0; a<2; ++a) {
3310 for (
int b=0; b<2; ++b) {
3311 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
3312 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
3313 for(
int k=0; k<3; ++k) {
3315 Cp1_chargino += -D2ch*
pow(V_CKM[k][2], 2)*
pow(Yd[2], 2)*Z_m[1][a]*Z_m[1][b]*Z_U[k][i]*Z_U[k][j]*
pow(conj(V_CKM[k][1]), 2)*
pow(conj(Yd[1]), 2)*conj(Z_m[1][a])*conj(Z_m[1][b])*conj(Z_U[k][i])*conj(Z_U[k][j])/(32.0*
pow(
PI, 2));
3330 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
3335 double otherc = sqrt(2.)/v1;
3342 M_ch0[0]=fabs(temp_ch0[0]);
3343 M_ch0[1]=fabs(temp_ch0[1]);
3344 M_ch0[2]=fabs(temp_ch0[2]);
3345 M_ch0[3]=fabs(temp_ch0[3]);
3347 for(
int i=0; i<4; ++i){
3348 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
3352 for(
int i=0; i<6; ++i) {
3353 M_D_pow_2[i]=
pow(M_D[i],2);
3358 for(
int i=0; i<4; ++i) {
3359 for(
int j=0; j<4; ++j) {
3365 for(
int i=0; i<4; ++i){
3366 if(temp_ch0[i]<0.) {
3367 for(
int j=0; j<4; ++j) {
3379 for(
int i=0; i<6; ++i) {
3380 for(
int j=0; j<6; ++j) {
3381 for(
int a=0; a<4; ++a) {
3382 for (
int b=0; b<4; ++b) {
3383 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
3384 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
3386 Cp1_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][b])/(3.0*cw) +
3387 Yd[2]*Z_D[2][j]*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][i])/(3.0*cw) +
3388 Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][i]))/(64.0*
pow(
PI, 2)) - D2ne*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][b])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
3389 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
3390 Yd[2]*Z_D[5][i]*Z_N[2][a])/(32.0*
pow(
PI, 2));
3407 double M_g_pow_2 =
pow(M_g,2.);
3415 for(
int i=0; i<6; ++i) {
3416 for(
int j=0; j<6; ++j) {
3417 for(
int a=0; a<4; ++a) {
3418 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
3419 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
3421 Cp1_mixed += -D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[1+3][i]*Z_D[1+3][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][j])/(3.0*cw) +
3422 Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][j])) + Z_D[5][i]*Z_D[5][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
3423 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j])))/(96.0*
pow(
PI, 2)) - D2mix*Z_D[5][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
3424 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*conj(Z_D[1+3][i])/(8.0*
pow(
PI, 2));
3492 double M_H_pow_2 =
pow(M_H,2.);
3493 double M_W_pow_2 =
pow(M_W,2.);
3494 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
3500 double m_u[4],m_u_pow_2[4];
3502 for (
int i = 0; i<3; ++i) {
3503 for (
int j = 0; j<3; j++) {
3512 for(
int i =0;i<3;++i) {
3513 m_u_pow_2[i]=
pow(m_u[i],2.);
3519 double D0h,D0h_c,D2h,D2h_c;
3523 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
3524 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
3525 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
3526 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
3527 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
3529 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
3531 C2_chargedhiggs += -
pow(g_2,4.)*
pow(m_q,2.)*CKM_product*
pow(m_u[i],2)*
pow(m_u[j],2)*(D0h_c - 2*D0h)/(128.*
pow(
PI,2.)*
pow(M_W,4.));
3538 double M_D[6],M_D_pow_2[6],dm[6];
3541 double Mg_pow_2 =
pow(Mg,2);
3551 for(
int i = 0; i<6; ++i) {
3552 for(
int j = 0; j<6; ++j) {
3556 for(
int i = 0; i<6; ++i) {
3557 M_D_pow_2[i]=
pow(M_D[i],2);
3566 for(
int i =0; i<6; ++i) {
3567 for(
int j = 0; j<6; ++j) {
3568 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
3569 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
3570 C2_gluino += -17.*D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[2][j]*
pow(g_3, 4.)*conj(Z_D[1+3][i])*conj(Z_D[1+3][j])/(288.*
pow(
PI, 2));
3579 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
3580 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
3596 for(
int i=0; i<2; ++i){
3597 M_ch_pow_2[i]=
pow(M_ch[i],2);
3599 for(
int i=0; i<6; ++i){
3600 M_U_pow_2[i]=
pow(M_U[i],2);
3602 for(
int i=0; i <2; ++i) {
3603 for(
int j=0; j<2; ++j) {
3607 for(
int i=2; i<2; ++i) {
3608 for(
int j=0; j<2; ++j) {
3612 for(
int i=0; i<6; ++i) {
3613 for(
int j=0; j<6; ++j) {
3628 double common = sqrt(2.)/v2;
3633 double otherc = sqrt(2.)/v1;
3646 for(
int i = 0; i<6; ++i) {
3647 for(
int j=0; j<6; ++j) {
3648 for(
int a =0; a<2; ++a) {
3649 for (
int b=0; b<2; ++b) {
3650 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
3651 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
3652 for(
int k=0; k<3; ++k) {
3667 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
3672 double otherc = sqrt(2.)/v1;
3679 M_ch0[0]=fabs(temp_ch0[0]);
3680 M_ch0[1]=fabs(temp_ch0[1]);
3681 M_ch0[2]=fabs(temp_ch0[2]);
3682 M_ch0[3]=fabs(temp_ch0[3]);
3684 for(
int i=0; i<4; ++i){
3685 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
3689 for(
int i=0; i<6; ++i) {
3690 M_D_pow_2[i]=
pow(M_D[i],2);
3695 for(
int i=0; i<4; ++i) {
3696 for(
int j=0; j<4; ++j) {
3702 for(
int i=0; i<4; ++i){
3703 if(temp_ch0[i]<0.) {
3704 for(
int j=0; j<4; ++j) {
3716 for(
int i=0; i<6; ++i) {
3717 for(
int j=0; j<6; ++j) {
3718 for(
int a=0; a<4; ++a) {
3719 for (
int b=0; b<4; ++b) {
3720 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
3721 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
3723 C2_neutralino += D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][j])/(3.0*cw) +
3724 Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/
3725 (2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])/(32.0*
pow(
PI, 2));
3743 double M_g_pow_2 =
pow(M_g,2.);
3751 for(
int i=0; i<6; ++i) {
3752 for(
int j=0; j<6; ++j) {
3753 for(
int a=0; a<4; ++a) {
3754 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
3755 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
3757 C2_mixed += D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[2][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
3758 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*conj(Z_D[2][j]) + 3.0*Z_D[2][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 -
3759 Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[1+3][i]))/(48.0*
pow(
PI, 2)) + D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
3760 Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[1+3][i])*conj(Z_D[1+3][j])/(48.0*
pow(
PI, 2));
3829 double M_H_pow_2 =
pow(M_H,2.);
3830 double M_W_pow_2 =
pow(M_W,2.);
3831 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
3837 double m_u[4],m_u_pow_2[4];
3839 for (
int i = 0; i<3; ++i) {
3840 for (
int j = 0; j<3; j++) {
3849 for(
int i =0;i<3;++i) {
3850 m_u_pow_2[i]=
pow(m_u[i],2.);
3856 double D0h,D0h_c,D2h,D2h_c;
3860 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
3861 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
3862 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
3863 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
3864 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
3866 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
3868 Cp2_chargedhiggs += -
pow(g_2,4.)*
pow(m_b,2.)*CKM_product*
pow(m_u[i],2)*
pow(m_u[j],2)*(D0h_c-2.*D0h)/(128.*
pow(
PI,2.)*
pow(M_W,4.));
3875 double M_D[6],M_D_pow_2[6],dm[6];
3878 double Mg_pow_2 =
pow(Mg,2);
3888 for(
int i = 0; i<6; ++i) {
3889 for(
int j = 0; j<6; ++j) {
3893 for(
int i = 0; i<6; ++i) {
3894 M_D_pow_2[i]=
pow(M_D[i],2);
3903 for(
int i =0; i<6; ++i) {
3904 for(
int j = 0; j<6; ++j) {
3905 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
3906 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
3908 Cp2_gluino += -17.*D0g*
pow(Mg, 2.)*Z_D[5][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[1][i])*conj(Z_D[1][j])/(288.*
pow(
PI, 2.));
3917 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
3918 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
3934 for(
int i=0; i<2; ++i){
3935 M_ch_pow_2[i]=
pow(M_ch[i],2);
3937 for(
int i=0; i<6; ++i){
3938 M_U_pow_2[i]=
pow(M_U[i],2);
3940 for(
int i=0; i <2; ++i) {
3941 for(
int j=0; j<2; ++j) {
3945 for(
int i=2; i<2; ++i) {
3946 for(
int j=0; j<2; ++j) {
3950 for(
int i=0; i<6; ++i) {
3951 for(
int j=0; j<6; ++j) {
3966 double common = sqrt(2.)/v2;
3971 double otherc = sqrt(2.)/v1;
3983 for(
int i = 0; i<6; ++i) {
3984 for(
int j=0; j<6; ++j) {
3985 for(
int a =0; a<2; ++a) {
3986 for (
int b=0; b<2; ++b) {
3987 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
3988 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
3989 for(
int k=0; k<3; ++k) {
4004 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
4009 double otherc = sqrt(2.)/v1;
4016 M_ch0[0]=fabs(temp_ch0[0]);
4017 M_ch0[1]=fabs(temp_ch0[1]);
4018 M_ch0[2]=fabs(temp_ch0[2]);
4019 M_ch0[3]=fabs(temp_ch0[3]);
4021 for(
int i=0; i<4; ++i){
4022 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
4026 for(
int i=0; i<6; ++i) {
4027 M_D_pow_2[i]=
pow(M_D[i],2);
4032 for(
int i=0; i<4; ++i) {
4033 for(
int j=0; j<4; ++j) {
4039 for(
int i=0; i<4; ++i){
4040 if(temp_ch0[i]<0.) {
4041 for(
int j=0; j<4; ++j) {
4053 for(
int i=0; i<6; ++i) {
4054 for(
int j=0; j<6; ++j) {
4055 for(
int a=0; a<4; ++a) {
4056 for (
int b=0; b<4; ++b) {
4057 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
4058 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
4060 Cp2_neutralino += D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
4061 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 -
4062 conj(Z_N[1][b])*cw)*conj(Z_D[1][j])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][j])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
4080 double M_g_pow_2 =
pow(M_g,2.);
4088 for(
int i=0; i<6; ++i) {
4089 for(
int j=0; j<6; ++j) {
4090 for(
int a=0; a<4; ++a) {
4091 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
4092 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
4094 Cp2_mixed += D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[5][i]*
pow(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) +
4095 conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]), 2)*conj(Z_D[5][j]) + 3.0*Z_D[5][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
4096 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1][i]))/(48.0*
pow(
PI, 2)) +
4097 D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
4098 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*conj(Z_D[1][i])*conj(Z_D[1][j])/(48.0*
pow(
PI, 2));
4168 double M_H_pow_2 =
pow(M_H,2.);
4169 double M_W_pow_2 =
pow(M_W,2.);
4170 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
4176 double m_u[4],m_u_pow_2[4];
4178 for (
int i = 0; i<3; ++i) {
4179 for (
int j = 0; j<3; j++) {
4188 for(
int i =0;i<3;++i) {
4189 m_u_pow_2[i]=
pow(m_u[i],2.);
4196 double D0h,D0h_c,D2h,D2h_c;
4200 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
4201 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
4202 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
4203 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
4204 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
4206 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
4208 C3_chargedhiggs += 0.;
4215 double M_D[6],M_D_pow_2[6],dm[6];
4218 double Mg_pow_2 =
pow(Mg,2);
4228 for(
int i = 0; i<6; ++i) {
4229 for(
int j = 0; j<6; ++j) {
4233 for(
int i = 0; i<6; ++i) {
4234 M_D_pow_2[i]=
pow(M_D[i],2);
4243 for(
int i =0; i<6; ++i) {
4244 for(
int j = 0; j<6; ++j) {
4245 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
4246 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
4248 C3_gluino += -D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[2][j]*
pow(g_3,4.)*conj(Z_D[1+3][i])*conj(Z_D[1+3][j])/(96.*
pow(
PI, 2));
4257 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
4258 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
4274 for(
int i=0; i<2; ++i){
4275 M_ch_pow_2[i]=
pow(M_ch[i],2);
4277 for(
int i=0; i<6; ++i){
4278 M_U_pow_2[i]=
pow(M_U[i],2);
4280 for(
int i=0; i <2; ++i) {
4281 for(
int j=0; j<2; ++j) {
4285 for(
int i=2; i<2; ++i) {
4286 for(
int j=0; j<2; ++j) {
4290 for(
int i=0; i<6; ++i) {
4291 for(
int j=0; j<6; ++j) {
4306 double common = sqrt(2.)/v2;
4311 double otherc = sqrt(2.)/v1;
4323 for(
int i = 0; i<6; ++i) {
4324 for(
int j=0; j<6; ++j) {
4325 for(
int a =0; a<2; ++a) {
4326 for (
int b=0; b<2; ++b) {
4327 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
4328 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
4329 for(
int k=0; k<3; ++k) {
4332 C3_chargino += -D0ch*M_ch[a]*M_ch[b]*
pow(V_CKM[k][2], 2)*Z_m[1][a]*Z_m[1][b]*Z_U[k][i]*Z_U[k][j]*(-Q_e*Z_p[0][a]*conj(Z_U[k][i])*swi + Yu[k]*Z_p[1][a]*conj(Z_U[k+3][i]))*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*
pow(conj(V_CKM[k][1]), 2)*
pow(conj(Yd[1]), 2)/(32.0*
pow(
PI, 2));
4346 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
4351 double otherc = sqrt(2.)/v1;
4358 M_ch0[0]=fabs(temp_ch0[0]);
4359 M_ch0[1]=fabs(temp_ch0[1]);
4360 M_ch0[2]=fabs(temp_ch0[2]);
4361 M_ch0[3]=fabs(temp_ch0[3]);
4363 for(
int i=0; i<4; ++i){
4364 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
4368 for(
int i=0; i<6; ++i) {
4369 M_D_pow_2[i]=
pow(M_D[i],2);
4374 for(
int i=0; i<4; ++i) {
4375 for(
int j=0; j<4; ++j) {
4381 for(
int i=0; i<4; ++i){
4382 if(temp_ch0[i]<0.) {
4383 for(
int j=0; j<4; ++j) {
4395 for(
int i=0; i<6; ++i) {
4396 for(
int j=0; j<6; ++j) {
4397 for(
int a=0; a<4; ++a) {
4398 for (
int b=0; b<4; ++b) {
4399 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
4400 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
4402 C3_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*((-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) +
4403 Yd[2]*Z_D[5][j]*Z_N[2][b]) - (-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
4404 Yd[2]*Z_D[5][j]*Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
4405 Yd[2]*Z_D[5][i]*Z_N[2][a])/(32.0*
pow(
PI, 2));
4423 double M_g_pow_2 =
pow(M_g,2.);
4431 for(
int i=0; i<6; ++i) {
4432 for(
int j=0; j<6; ++j) {
4433 for(
int a=0; a<4; ++a) {
4434 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
4435 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
4437 C3_mixed += D0mix*M_ch0[a]*M_g*Z_D[2][i]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
4438 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*conj(Z_D[2][j])/(48.0*
pow(
PI, 2)) - D0mix*M_ch0[a]*M_g*Z_D[2][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
4439 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[1+3][i])/(48.0*
pow(
PI, 2)) +
4440 D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
4441 Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[1+3][i])*conj(Z_D[1+3][j])/(48.0*
pow(
PI, 2));
4508 double M_H_pow_2 =
pow(M_H,2.);
4509 double M_W_pow_2 =
pow(M_W,2.);
4510 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
4516 double m_u[4],m_u_pow_2[4];
4518 for (
int i = 0; i<3; ++i) {
4519 for (
int j = 0; j<3; j++) {
4528 for(
int i =0;i<3;++i) {
4529 m_u_pow_2[i]=
pow(m_u[i],2.);
4535 double D0h,D0h_c,D2h,D2h_c;
4539 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
4540 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
4541 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
4542 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
4543 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
4545 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
4547 Cp3_chargedhiggs += 0.;
4553 double M_D[6],M_D_pow_2[6],dm[6];
4556 double Mg_pow_2 =
pow(Mg,2);
4566 for(
int i = 0; i<6; ++i) {
4567 for(
int j = 0; j<6; ++j) {
4571 for(
int i = 0; i<6; ++i) {
4572 M_D_pow_2[i]=
pow(M_D[i],2);
4580 for(
int i =0; i<6; ++i) {
4581 for(
int j = 0; j<6; ++j) {
4582 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
4583 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
4585 Cp3_gluino += -D0g*
pow(Mg, 2)*Z_D[5][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[1][i])*conj(Z_D[1][j])/(96.*
pow(
PI, 2));
4593 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
4594 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
4610 for(
int i=0; i<2; ++i){
4611 M_ch_pow_2[i]=
pow(M_ch[i],2);
4613 for(
int i=0; i<6; ++i){
4614 M_U_pow_2[i]=
pow(M_U[i],2);
4616 for(
int i=0; i <2; ++i) {
4617 for(
int j=0; j<2; ++j) {
4621 for(
int i=2; i<2; ++i) {
4622 for(
int j=0; j<2; ++j) {
4626 for(
int i=0; i<6; ++i) {
4627 for(
int j=0; j<6; ++j) {
4642 double common = sqrt(2.)/v2;
4647 double otherc = sqrt(2.)/v1;
4659 for(
int i = 0; i<6; ++i) {
4660 for(
int j=0; j<6; ++j) {
4661 for(
int a =0; a<2; ++a) {
4662 for (
int b=0; b<2; ++b) {
4663 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
4664 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
4665 for(
int k=0; k<3; ++k) {
4667 Cp3_chargino += -D0ch*M_ch[a]*M_ch[b]*
pow(V_CKM[k][2], 2)*
pow(Yd[2], 2)*(-Q_e*Z_U[k][i]*conj(Z_p[0][b])*swi + Z_U[k+3][i]*conj(Yu[k])*conj(Z_p[1][b]))*(-Q_e*Z_U[k][j]*conj(Z_p[0][a])*swi + Z_U[k+3][j]*conj(Yu[k])*conj(Z_p[1][a]))*
pow(conj(V_CKM[k][1]), 2)*conj(Z_m[1][a])*conj(Z_m[1][b])*conj(Z_U[k][i])*conj(Z_U[k][j])/(32.0*
pow(
PI, 2));
4679 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
4684 double otherc = sqrt(2.)/v1;
4691 M_ch0[0]=fabs(temp_ch0[0]);
4692 M_ch0[1]=fabs(temp_ch0[1]);
4693 M_ch0[2]=fabs(temp_ch0[2]);
4694 M_ch0[3]=fabs(temp_ch0[3]);
4696 for(
int i=0; i<4; ++i){
4697 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
4701 for(
int i=0; i<6; ++i) {
4702 M_D_pow_2[i]=
pow(M_D[i],2);
4707 for(
int i=0; i<4; ++i) {
4708 for(
int j=0; j<4; ++j) {
4714 for(
int i=0; i<4; ++i){
4715 if(temp_ch0[i]<0.) {
4716 for(
int j=0; j<4; ++j) {
4727 for(
int i=0; i<6; ++i) {
4728 for(
int j=0; j<6; ++j) {
4729 for(
int a=0; a<4; ++a) {
4730 for (
int b=0; b<4; ++b) {
4731 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
4732 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
4734 Cp3_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][j])/(2.0*cw*sw) +
4735 conj(Yd[1])*conj(Z_D[1+3][j])*conj(Z_N[2][b])) + (-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][b])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][b]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
4736 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a])))*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
4737 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][b]))/(32.0*
pow(
PI, 2));
4753 double M_g_pow_2 =
pow(M_g,2.);
4761 for(
int i=0; i<6; ++i) {
4762 for(
int j=0; j<6; ++j) {
4763 for(
int a=0; a<4; ++a) {
4764 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
4765 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
4767 Cp3_mixed += D0mix*M_ch0[a]*M_g*Z_D[5][i]*
pow(g_3, 2)*
pow(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) +
4768 conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]), 2)*conj(Z_D[5][j])/(48.0*
pow(
PI, 2)) - D0mix*M_ch0[a]*M_g*Z_D[5][j]*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
4769 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) +
4770 conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1][i])/(48.0*
pow(
PI, 2)) + D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) +
4771 Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*conj(Z_D[1][i])*conj(Z_D[1][j])/(48.0*
pow(
PI, 2));
4838 double M_H_pow_2 =
pow(M_H,2.);
4839 double M_W_pow_2 =
pow(M_W,2.);
4840 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
4846 double m_u[4],m_u_pow_2[4];
4848 for (
int i = 0; i<3; ++i) {
4849 for (
int j = 0; j<3; j++) {
4858 for(
int i =0;i<3;++i) {
4859 m_u_pow_2[i]=
pow(m_u[i],2.);
4865 double D0h,D0h_c,D2h,D2h_c;
4869 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
4870 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
4871 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
4872 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
4873 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
4875 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
4877 C4_chargedhiggs +=
pow(g_2,4.)*CKM_product*(m_b*m_q*D2h*
pow(tbeta,2.)/
pow(M_W,2.) - m_b*m_q*
pow(m_u[i],2)*
pow(m_u[j],2)*(D0h_c + D0h*(
pow(tbeta,2.) +
pow(tbeta,-2.)))/(4*
pow(M_W,4.)))/(16.*
pow(
PI,2.));
4884 double M_D[6],M_D_pow_2[6],dm[6];
4887 double Mg_pow_2 =
pow(Mg,2);
4897 for(
int i = 0; i<6; ++i) {
4898 for(
int j = 0; j<6; ++j) {
4902 for(
int i = 0; i<6; ++i) {
4903 M_D_pow_2[i]=
pow(M_D[i],2);
4912 for(
int i =0; i<6; ++i) {
4913 for(
int j = 0; j<6; ++j) {
4914 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
4915 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
4917 C4_gluino += -7.*D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[1][i])*conj(Z_D[1+3][j])/(48.*
pow(
PI, 2)) + D2g*
pow(g_3,4.)*(6.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[1][i])*conj(Z_D[1+3][j]) + 11.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[1][j])*conj(Z_D[1+3][i]))/(72.*
pow(
PI, 2));
4926 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
4927 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
4943 for(
int i=0; i<2; ++i){
4944 M_ch_pow_2[i]=
pow(M_ch[i],2);
4946 for(
int i=0; i<6; ++i){
4947 M_U_pow_2[i]=
pow(M_U[i],2);
4949 for(
int i=0; i <2; ++i) {
4950 for(
int j=0; j<2; ++j) {
4954 for(
int i=2; i<2; ++i) {
4955 for(
int j=0; j<2; ++j) {
4959 for(
int i=0; i<6; ++i) {
4960 for(
int j=0; j<6; ++j) {
4975 double common = sqrt(2.)/v2;
4980 double otherc = sqrt(2.)/v1;
4992 for(
int i = 0; i<6; ++i) {
4993 for(
int j=0; j<6; ++j) {
4994 for(
int a =0; a<2; ++a) {
4995 for (
int b=0; b<2; ++b) {
4996 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
4997 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
4998 for(
int k=0; k<3; ++k) {
5000 C4_chargino += D2ch*
pow(V_CKM[k][2], 2)*Yd[2]*Z_m[1][a]*Z_U[k][j]*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*(-Q_e*Z_U[k][i]*conj(Z_p[0][b])*swi + Z_U[k+3][i]*conj(Yu[k])*conj(Z_p[1][b]))*
pow(conj(V_CKM[k][1]), 2)*conj(Yd[1])*conj(Z_m[1][a])*conj(Z_U[k][i])/(8.0*
pow(
PI, 2));
5013 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
5018 double otherc = sqrt(2.)/v1;
5025 M_ch0[0]=fabs(temp_ch0[0]);
5026 M_ch0[1]=fabs(temp_ch0[1]);
5027 M_ch0[2]=fabs(temp_ch0[2]);
5028 M_ch0[3]=fabs(temp_ch0[3]);
5030 for(
int i=0; i<4; ++i){
5031 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
5035 for(
int i=0; i<6; ++i) {
5036 M_D_pow_2[i]=
pow(M_D[i],2);
5041 for(
int i=0; i<4; ++i) {
5042 for(
int j=0; j<4; ++j) {
5048 for(
int i=0; i<4; ++i){
5049 if(temp_ch0[i]<0.) {
5050 for(
int j=0; j<4; ++j) {
5062 for(
int i=0; i<6; ++i) {
5063 for(
int j=0; j<6; ++j) {
5064 for(
int a=0; a<4; ++a) {
5065 for (
int b=0; b<4; ++b) {
5066 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
5067 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
5069 C4_neutralino += D2ne*((-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) +
5070 Yd[2]*Z_D[5][j]*Z_N[2][b]) + (-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
5071 Yd[2]*Z_D[5][j]*Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) +
5072 conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][b]))/(8.0*
pow(
PI, 2));
5090 double M_g_pow_2 =
pow(M_g,2.);
5098 for(
int i=0; i<6; ++i) {
5099 for(
int j=0; j<6; ++j) {
5100 for(
int a=0; a<4; ++a) {
5101 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
5102 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
5104 C4_mixed += D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5105 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1+3][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
5106 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[1][i]))/(16.0*
pow(
PI, 2)) -
5107 D2mix*
pow(g_3, 2)*(-3.0*Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5108 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a])) - 3.0*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][i])/(3.0*cw) +
5109 Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][i]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[1][j])*conj(Z_D[1+3][i]))/(24.0*
pow(
PI, 2)) -
5110 D2mix*
pow(g_3, 2)*(-Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5111 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a])) - (-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
5112 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[1][j])*conj(Z_D[1+3][i]))/(24.0*
pow(
PI, 2)) -
5113 D2mix*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
5114 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*conj(Z_D[1][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) +
5115 Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 - conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1+3][i]))/(24.0*
pow(
PI, 2));
5126 scalar_t delta_d[6][6],delta_d_LL[3][3],delta_d_LR[3][3],delta_d_RL[3][3],delta_d_RR[3][3];
5127 scalar_t delta_u[6][6],delta_u_LL[3][3],delta_u_LR[3][3],delta_u_RL[3][3],delta_u_RR[3][3];
5129 double m_av = (M_U[0]+M_U[1]+M_U[2]+M_U[3]+M_U[4]+M_U[5]+M_D[0]+M_D[1]+M_D[2]+M_D[3]+M_D[4]+M_D[5])/12.;
5131 getDelta(delta_d,Z_D,M_D,m_av,delta_d_LL,delta_d_LR,delta_d_RL,delta_d_RR);
5132 getDelta(delta_u,Z_U,M_U,m_av,delta_u_LL,delta_u_LR,delta_u_RL,delta_u_RR);
5138 double x_mu =
pow(abs(mu),2)/
pow(m_av,2);
5140 double x_g =
pow(M_g,2)/
pow(m_av,2);
5141 double alpha_s =
pow(g_3,2.)/(4.*
PI) ;
5142 double alpha_2 =
pow(g_2,2.)/(4.*
PI);
5143 double eps=2*alpha_s*mu*M_g*
f(x_g)/(3.*
PI*
pow(m_av,2));
5149 scalar_t a2 = -alpha_s*
pow(M_g,2)*delta_d_LL[2][1]*delta_d_RR[2][1]*
pow(
h1(x_g),2);
5150 scalar_t a3 = alpha_2*
pow(m_t,2)*A_t*M_g*
h1(x_g)*
h3(x_mu)*delta_d_RR[2][1]*V_tb*conj(V_tq)/(
pow(M_W,2));
5151 scalar_t a4 = alpha_2*M_2*M_g*delta_u_LL[2][1]*delta_d_RR[2][1]*
h1(x_g)*
h4(x_2,x_g);
5153 scalar_t C4_higgspenguin = a1*(a2+a3+a4);
5214 double M_H_pow_2 =
pow(M_H,2.);
5215 double M_W_pow_2 =
pow(M_W,2.);
5216 std::array<std::array<scalar_t, 3>, 3> V_CKM {};
5222 double m_u[4],m_u_pow_2[4];
5224 for (
int i = 0; i<3; ++i) {
5225 for (
int j = 0; j<3; j++) {
5234 for(
int i =0;i<3;++i) {
5235 m_u_pow_2[i]=
pow(m_u[i],2.);
5242 double D0h,D0h_c,D2h,D2h_c;
5246 for(
int i = 0; i<3; i++)
for(
int j=1; j<3; j++) {
5247 D0h =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
5248 D0h_c =
D0(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
5249 D2h =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_W_pow_2);
5250 D2h_c =
D2p(m_u_pow_2[i],m_u_pow_2[j],M_H_pow_2,M_H_pow_2);
5252 CKM_product = V_CKM[i][3]*V_CKM[j][3]*conj(V_CKM[i][1])*conj(V_CKM[j][1]);
5254 C5_chargedhiggs +=
pow(g_2,4.)*m_b*m_q*CKM_product*
pow(m_u[i],2.)*(D2h_c-2.*D2h)/(32.*
pow(
PI,2.)*
pow(M_W,4.));
5261 double M_D[6],M_D_pow_2[6],dm[6];
5264 double Mg_pow_2 =
pow(Mg,2);
5274 for(
int i = 0; i<6; ++i) {
5275 for(
int j = 0; j<6; ++j) {
5279 for(
int i = 0; i<6; ++i) {
5280 M_D_pow_2[i]=
pow(M_D[i],2);
5290 for(
int i =0; i<6; ++i) {
5291 for(
int j = 0; j<6; ++j) {
5292 D2g =
D2p(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
5293 D0g =
D0(M_D_pow_2[i],M_D_pow_2[j],Mg_pow_2,Mg_pow_2);
5295 C5_gluino += -D0g*
pow(Mg, 2)*Z_D[2][i]*Z_D[5][j]*
pow(g_3,4.)*conj(Z_D[1][i])*conj(Z_D[1+3][j])/(144.*
pow(
PI, 2)) + 5.*D2g*
pow(g_3,4.)*(-2.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[1][i])*conj(Z_D[1+3][j]) + 3.*Z_D[2][i]*Z_D[5][j]*conj(Z_D[1][j])*conj(Z_D[1+3][i]))/(72.*
pow(
PI, 2));
5304 double M_ch[2],M_ch_pow_2[2],M_U[6],M_U_pow_2[6];
5305 scalar_t Z_p[2][2],Z_m[2][2],Z_U[6][6];
5321 for(
int i=0; i<2; ++i){
5322 M_ch_pow_2[i]=
pow(M_ch[i],2);
5324 for(
int i=0; i<6; ++i){
5325 M_U_pow_2[i]=
pow(M_U[i],2);
5327 for(
int i=0; i <2; ++i) {
5328 for(
int j=0; j<2; ++j) {
5332 for(
int i=2; i<2; ++i) {
5333 for(
int j=0; j<2; ++j) {
5337 for(
int i=0; i<6; ++i) {
5338 for(
int j=0; j<6; ++j) {
5353 double common = sqrt(2.)/v2;
5358 double otherc = sqrt(2.)/v1;
5371 for(
int i = 0; i<6; ++i) {
5372 for(
int j=0; j<6; ++j) {
5373 for(
int a =0; a<2; ++a) {
5374 for (
int b=0; b<2; ++b) {
5375 D0ch =
D0(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
5376 D2ch =
D2p(M_U_pow_2[i],M_U_pow_2[j],M_ch_pow_2[a],M_ch_pow_2[b]);
5377 for(
int k=0; k<3; ++k) {
5379 C5_chargino += -D0ch*M_ch[a]*M_ch[b]*
pow(V_CKM[k][2], 2)*Yd[2]*Z_m[1][b]*Z_U[k][i]*(-Q_e*Z_p[0][b]*conj(Z_U[k][j])*swi + Yu[k]*Z_p[1][b]*conj(Z_U[k+3][j]))*(-Q_e*Z_U[k][j]*conj(Z_p[0][a])*swi + Z_U[k+3][j]*conj(Yu[k])*conj(Z_p[1][a]))*
pow(conj(V_CKM[k][1]), 2)*conj(Yd[1])*conj(Z_m[1][a])*conj(Z_U[k][i])/(16.0*
pow(
PI, 2));
5392 double M_ch0[4],M_ch0_pow_2[4],M_D[6],M_D_pow_2[6];
5397 double otherc = sqrt(2.)/v1;
5404 M_ch0[0]=fabs(temp_ch0[0]);
5405 M_ch0[1]=fabs(temp_ch0[1]);
5406 M_ch0[2]=fabs(temp_ch0[2]);
5407 M_ch0[3]=fabs(temp_ch0[3]);
5409 for(
int i=0; i<4; ++i){
5410 M_ch0_pow_2[i]=
pow(M_ch0[i],2);
5414 for(
int i=0; i<6; ++i) {
5415 M_D_pow_2[i]=
pow(M_D[i],2);
5420 for(
int i=0; i<4; ++i) {
5421 for(
int j=0; j<4; ++j) {
5427 for(
int i=0; i<4; ++i){
5428 if(temp_ch0[i]<0.) {
5429 for(
int j=0; j<4; ++j) {
5441 for(
int i=0; i<6; ++i) {
5442 for(
int j=0; j<6; ++j) {
5443 for(
int a=0; a<4; ++a) {
5444 for (
int b=0; b<4; ++b) {
5445 D0ne =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
5446 D2ne =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_ch0_pow_2[b]);
5448 C5_neutralino += -D0ne*M_ch0[a]*M_ch0[b]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][i])/(3.0*cw) +
5449 Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][b]*sw/3.0 - Z_N[1][b]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][b])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5450 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))/(16.0*
pow(
PI, 2)) - D2ne*(-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
5451 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][b]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][b]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0- Z_N[1][a]*cw)/(2.0*cw*sw) +
5452 Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][b])*sw/3.0 - conj(Z_N[1][b])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][b]))/(8.0*
pow(
PI, 2));
5468 double M_g_pow_2 =
pow(M_g,2.);
5476 for(
int i=0; i<6; ++i) {
5477 for(
int j=0; j<6; ++j) {
5478 for(
int a=0; a<4; ++a) {
5479 D0mix =
D0(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
5480 D2mix =
D2p(M_D_pow_2[i],M_D_pow_2[j],M_ch0_pow_2[a],M_g_pow_2);
5482 C5_mixed += -D0mix*M_ch0[a]*M_g*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5483 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1+3][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
5484 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[1][i]))/(48.0*
pow(
PI, 2)) +
5485 D2mix*
pow(g_3, 2)*(-Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][i])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][i]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5486 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a])) - (-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[5][i])/(3*cw) +
5487 Z_N[2][a]*conj(Yd[2])*conj(Z_D[2][i]))*(-sqrt(2)*Q_e*Z_D[2][j]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][j]*Z_N[2][a])*conj(Z_D[1][j])*conj(Z_D[1+3][i]))/(24.0*
pow(
PI, 2)) +
5488 D2mix*
pow(g_3, 2)*(-Z_D[2][j]*Z_D[5][i]*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) + Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5489 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a])) - (-sqrt(2)*Q_e*Z_D[5][j]*conj(Z_N[0][a])/(3.0*cw) +
5490 Yd[2]*Z_D[2][j]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*conj(Z_D[1][j])*conj(Z_D[1+3][i]))/(8.0*
pow(
PI, 2)) +
5491 D2mix*
pow(g_3, 2)*(Z_D[2][j]*(-sqrt(2)*Q_e*Z_D[5][i]*conj(Z_N[0][a])/(3.0*cw) + Yd[2]*Z_D[2][i]*conj(Z_N[2][a]))*(-sqrt(2)*Q_e*Z_N[0][a]*conj(Z_D[1+3][j])/(3.0*cw) +
5492 Z_N[2][a]*conj(Yd[1])*conj(Z_D[1][j]))*conj(Z_D[1][i]) + Z_D[5][j]*(-sqrt(2)*Q_e*Z_D[2][i]*(Z_N[0][a]*sw/3.0 - Z_N[1][a]*cw)/(2.0*cw*sw) + Yd[2]*Z_D[5][i]*Z_N[2][a])*(-sqrt(2)*Q_e*(conj(Z_N[0][a])*sw/3.0 -
5493 conj(Z_N[1][a])*cw)*conj(Z_D[1][i])/(2.0*cw*sw) + conj(Yd[1])*conj(Z_D[1+3][i])*conj(Z_N[2][a]))*conj(Z_D[1+3][i]))/(8.0*
pow(
PI, 2));