10 #include <gsl/gsl_integration.h>
11 #include <gsl/gsl_sf_dilog.h>
18 myGTHDM(GeneralTHDM_i),
65 double eta = M_PI/2.0-(beta-alpha1);
72 double yl_h = 1.0 + eta*
sl.
real();
73 double yl_H = -
sl.
real() + eta;
76 double rmu_hSM, rmu_h, rmu_H, rmu_A, rmu_Hp;
77 double part_hSM, part_h, part_H, part_A, part_Hp;
81 throw std::runtime_error(
"The implemented approximation for g-2_\\mu only works for Higgs masses above 1 GeV.");
86 throw std::runtime_error(
"(g-2) is only aviable in the A2HDM.");
91 throw std::runtime_error(
"(g-2) is only aviable in the CP-conserving limit.");
101 part_hSM=rmu_hSM*(-7.0/6.0-
log(rmu_hSM));
102 part_h=yl_h*yl_h*rmu_h*(-7.0/6.0-
log(rmu_h));
103 part_H=yl_H*yl_H*rmu_H*(-7.0/6.0-
log(rmu_H));
104 part_A=yl_A*yl_A*rmu_A*(11.0/6.0+
log(rmu_A));
105 part_Hp=-yl_A*yl_A*rmu_Hp*1.0/6.0;
118 __f_params ¶ms= *reinterpret_cast<__f_params *>(p);
119 double integ = (2.0*x*(1.0-x)-1.0) *
log(x*(1.0-x)/params.
a) / (x*(1.0-x)-params.
a);
124 __f_params ¶ms= *reinterpret_cast<__f_params *>(p);
125 double integ =
log(x*(1.0-x)/params.
a) / (x*(1.0-x)-params.
a);
136 gsl_integration_glfixed_table * w
137 = gsl_integration_glfixed_table_alloc(100);
141 F.params = reinterpret_cast<void *>(¶ms);
143 result = gsl_integration_glfixed (&F, 0, 1, w);
145 gsl_integration_glfixed_table_free (w);
156 gsl_integration_glfixed_table * w
157 = gsl_integration_glfixed_table_alloc(100);
161 F.params = reinterpret_cast<void *>(¶ms);
163 result = gsl_integration_glfixed (&F, 0, 1, w);
165 gsl_integration_glfixed_table_free (w);
212 ml=0.0; mm=0.0; mh=0.0;
220 else if(m3<=m2 && m2<=m1)
227 else if(m2<=m1 && m1<=m3)
234 else if(m2<=m3 && m3<=m1)
241 else if(m3<=m1 && m1<=m2)
249 else if(m1<=m2 && m2<=m3)
259 - 2.0*ml*ml*mm*mm - 2.0*mm*mm*mh*mh - 2.0*mh*mh*ml*ml);
260 gslpp::complex ap = 1.0/(2.0*mh*mh)*( mh*mh + ml*ml - mm*mm - lambda);
261 gslpp::complex am = 1.0/(2.0*mh*mh)*( mh*mh - ml*ml + mm*mm - lambda);
275 double aFNphoton, aFNZ;
295 double SW4 = sW2*sW2;
317 double eta = M_PI/2.0-(beta-alpha1);
322 double mTAU2 = mTAU*mTAU;
326 double yt_h = 1.0 + eta*
su.
real();
327 double yb_h = 1.0 + eta*
sd.
real();
328 double ytau_h = 1.0 + eta*
sl.
real();
330 double yt_H = -
su.
real() + eta;
331 double yb_H = -
sd.
real() + eta;
332 double ytau_H = -
sl.
real() + eta;
335 double yb_A = -
sd.
real();
336 double ytau_A = -
sl.
real();
338 double rtau_hSM, rtau_h, rtau_H, rtau_A, rt_hSM, rt_h, rt_H, rt_A, rb_hSM, rb_h, rb_H, rb_A;
339 double stau_hSM, stau_h, stau_H, stau_A, st_hSM, st_h, st_H, st_A, sb_hSM, sb_h, sb_H, sb_A;
341 double part_hSM_photon, part_h_photon, part_H_photon, part_A_photon;
342 double part_hSM_Z, part_h_Z, part_H_Z, part_A_Z;
347 throw std::runtime_error(
"The implemented approximation for g-2_\\mu only works for Higgs masses above 1 GeV.");
352 rtau_hSM=mTAU*mTAU/
mH1_2;
353 rtau_h=mTAU*mTAU/
mH1_2;
354 rtau_H=mTAU*mTAU/
mH2_2;
355 rtau_A=mTAU*mTAU/
mH3_2;
365 stau_hSM=mTAU*mTAU/(
mH1_2 - MZ*MZ);
366 stau_h=mTAU*mTAU/(
mH1_2-MZ*MZ);
367 stau_H=mTAU*mTAU/(
mH2_2-MZ*MZ);
368 stau_A=mTAU*mTAU/(
mH3_2-MZ*MZ);
369 st_hSM=mt*mt/(15647.5081-MZ*MZ);
370 st_h=mt*mt/(
mH1_2-MZ*MZ);
371 st_H=mt*mt/(
mH2_2-MZ*MZ);
372 st_A=mt*mt/(
mH3_2-MZ*MZ);
373 sb_hSM=mb*mb/(
mH1_2 -MZ*MZ);
374 sb_h=mb*mb/(
mH1_2-MZ*MZ);
375 sb_H=mb*mb/(
mH2_2-MZ*MZ);
376 sb_A=mb*mb/(
mH3_2-MZ*MZ);
379 part_h_photon=ytau_h*(ytau_h*rtau_h*
gscalar2(rtau_h)+(4.0/3.0)*yt_h*rt_h*
gscalar2(rt_h)+(1.0/3.0)*yb_h*rb_h*
gscalar2(rb_h));
380 part_H_photon=ytau_H*(ytau_H*rtau_H*
gscalar2(rtau_H)+(4.0/3.0)*yt_H*rt_H*
gscalar2(rt_H)+(1.0/3.0)*yb_H*rb_H*
gscalar2(rb_H));
383 part_hSM_Z = stau_hSM*(3.0/4.0 - cW2)*(3.0/4.0 - cW2)*(
gscalar2(rtau_hSM)-
gscalar2(mTAU*mTAU/(MZ*MZ)))/(cW2*(1-cW2)) + (-2.0)*st_hSM*(3.0/4.0 - cW2)*(-5.0/12.0 + 2.0/3.0*cW2)*(
gscalar2(rt_hSM)-
gscalar2(mt*mt/(MZ*MZ)))/(cW2*(1-cW2))+ (1.0)*sb_hSM*(3.0/4.0 - cW2)*(1.0/12.0 + 1.0/3.0*cW2)*(
gscalar2(rb_hSM)-
gscalar2(mb*mb/(MZ*MZ)))/(cW2*(1-cW2));
384 part_h_Z = ytau_h*(ytau_h*stau_h*(3.0/4.0 - cW2)*(3.0/4.0 - cW2)*(
gscalar2(rtau_h)-
gscalar2(mTAU*mTAU/(MZ*MZ)))/(cW2*(1-cW2)) + (-2.0)*yt_h*st_h*(3.0/4.0 - cW2)*(-5.0/12.0 + 2.0/3.0*cW2)*(
gscalar2(rt_h)-
gscalar2(mt*mt/(MZ*MZ)))/(cW2*(1-cW2))+ (1.0)*yb_h*sb_h*(3.0/4.0 - cW2)*(1.0/12.0 + 1.0/3.0*cW2)*(
gscalar2(rb_h)-
gscalar2(mb*mb/(MZ*MZ)))/(cW2*(1-cW2)));
385 part_H_Z = ytau_H*(ytau_H*stau_h*(3.0/4.0 - cW2)*(3.0/4.0 - cW2)*(
gscalar2(rtau_H)-
gscalar2(mTAU*mTAU/(MZ*MZ)))/(cW2*(1-cW2)) + (-2.0)*yt_H*st_H*(3.0/4.0 - cW2)*(-5.0/12.0 + 2.0/3.0*cW2)*(
gscalar2(rt_H)-
gscalar2(mt*mt/(MZ*MZ)))/(cW2*(1-cW2))+ (1.0)*yb_H*sb_H*(3.0/4.0 - cW2)*(1.0/12.0 + 1.0/3.0*cW2)*(
gscalar2(rb_H)-
gscalar2(mb*mb/(MZ*MZ)))/(cW2*(1-cW2)));
389 aFNphoton =
GF*
mMU*
mMU/(4.0*pi*pi*pi*
sqrt(2.0)) * aem * (-part_hSM_photon+part_h_photon+part_H_photon+part_A_photon);
390 aFNZ =
GF*
mMU*
mMU/(4.0*pi*pi*pi*
sqrt(2.0)) * aem * (-part_hSM_Z+part_h_Z+part_H_Z+part_A_Z);
393 aFC =
sl*((aem*aem*mTAU*mTAU*
mMU*
mMU*
395 mTAU*mTAU/MW2 + (-1.0/2.0 + mTAU*mTAU/mHp2)*
396 log (mTAU*mTAU/mHp2) - (-1.0/2.0 + mTAU*mTAU/MW2)*
397 log (mTAU*mTAU/MW2) + (mTAU*
398 mTAU*(-1.0 + mTAU*mTAU/mHp2)*(-pi*pi/6.0 +
401 mTAU*(-1.0 + mTAU*mTAU/MW2)*(-pi*pi/6.0 +
403 MW2*(mHp2 - MW2)*pi*pi*SW4)- (3.0*aem*aem*
mMU*
mMU*mt2*
404 su*(mb2/mHp2 - mt2/mHp2 - mb2/MW2 +
405 mt2/MW2 + (13.0/12.0 + mb2/mHp2)*
406 log (mb2/mHp2) + (13.0/12.0 - mt2/mHp2)*
log (mt2/mHp2) +
407 (-
log (mb2/mHp2)*
log (mb2/mHp2) +
log (mt2/mHp2)*
log (mt2/mHp2))/
408 3.0 - (13.0/12.0 + mb2/MW2)*
409 log (mb2/MW2) - (13.0/12.0 - mt2/MW2)*
410 log (mt2/MW2) + (
log (mb2/MW2)*
log (mb2/MW2) -
411 log (mt2/MW2)*
log (mt2/MW2))/
412 3.0 - (4.0*(-1.0 - mb2/mHp2 + mt2/mHp2)*
414 1.0))/(3.0*(1.0 + (-(mb2/mHp2) + mt2/mHp2)*(-(mb2/mHp2) +
417 mt2/mHp2))) + ((-((mb2*(5.0/3.0 + mb2/mHp2))/
418 mHp2) + (mt2*(-8.0/3.0 + mt2/mHp2))/mHp2)*
420 1.0))/(1.0 + (-(mb2/mHp2) + mt2/mHp2)*(-(mb2/mHp2) +
422 2.0*(mb2/mHp2 + mt2/mHp2))+ ((5.0*mb2)/(3.0*mHp2) - (8.0*
423 mt2)/(3.0*mHp2) + (-(mb2/mHp2) +
424 mt2/mHp2)*(-(mb2/mHp2) +
425 mt2/mHp2))*(-(
log (mb2/mt2)*
log (mt2/mHp2))/2.0 +
426 PolyLog.
Li2 (1.0 - mb2/mt2) - ((-(mb2/mHp2) + mt2/mHp2)*
428 1.0))/(1.0 + (-(mb2/mHp2) + mt2/mHp2)*(-(mb2/mHp2) +
430 2.0*(mb2/mHp2 + mt2/mHp2))) + (4.0*(-1.0 - mb2/MW2 +
433 1.0))/(3.0*(1.0 + (-(mb2/MW2) + mt2/MW2)*(-(mb2/MW2) +
436 mt2/MW2))) - ((-((mb2*(5.0/3.0 + mb2/MW2))/
437 MW2) + (mt2*(-8.0/3.0 + mt2/MW2))/MW2)*
439 1.0))/(1.0 + (-(mb2/MW2) + mt2/MW2)*(-(mb2/MW2) +
441 2.0*(mb2/MW2 + mt2/MW2))- ((-(mb2/MW2) +
442 mt2/MW2)*(-(mb2/MW2) + mt2/MW2) + (5.0*mb2)/(3.0*
443 MW2) - (8.0*mt2)/(3.0*
444 MW2))*(-(
log (mb2/mt2)*
log (mt2/MW2))/2.0 +
445 PolyLog.
Li2 (1.0 - mb2/mt2) - ((-(mb2/MW2) + mt2/MW2)*
447 1.0))/(1.0 + (-(mb2/MW2) + mt2/MW2)*(-(mb2/MW2) +
448 mt2/MW2) - 2.0*(mb2/MW2 + mt2/MW2)))))/(32.0*
449 MW2*(mHp2 - MW2)*pi*pi*SW4) + (3.0*aem*aem*mb2*
mMU*
mMU*
450 sd*(mb2/mHp2 - mt2/mHp2 - mb2/MW2 +
451 mt2/MW2 + (1.0/12.0 + mb2/mHp2)*
452 log (mb2/mHp2) + (1.0/12.0 - mt2/mHp2)*
log (mt2/mHp2) - (1.0/12.0 + mb2/MW2)*
453 log (mb2/MW2) - (1.0/12.0 - mt2/MW2)*
455 MW2) + ((-((mb2*(-1.0/3.0 + mb2/mHp2))/
456 mHp2) + (mt2*(-2.0/3.0 + mt2/mHp2))/mHp2)*
458 1.0))/(1.0 + (-(mb2/mHp2) + mt2/mHp2)*(-(mb2/mHp2) +
460 2.0*(mb2/mHp2 + mt2/mHp2)) + (-mb2/(3.0*mHp2) - (2.0*
461 mt2)/(3.0*mHp2) + (-(mb2/mHp2) +
462 mt2/mHp2)*(-(mb2/mHp2) +
464 mt2/mHp2))*(-(
log (mb2/mt2)*
log (mt2/mHp2))/2.0 +
465 PolyLog.
Li2 (1.0 - mb2/mt2) - ((-(mb2/mHp2) + mt2/mHp2)*
467 1.0))/(1.0 + (-(mb2/mHp2) + mt2/mHp2)*(-(mb2/mHp2) +
470 mt2/mHp2))) - ((-((mb2*(-1.0/3.0 + mb2/MW2))/
471 MW2) + (mt2*(-2.0/3.0 + mt2/MW2))/MW2)*
473 1.0))/(1.0 + (-(mb2/MW2) + mt2/MW2)*(-(mb2/MW2) +
475 2.0*(mb2/MW2 + mt2/MW2)) - ((-(mb2/MW2) +
476 mt2/MW2)*(-(mb2/MW2) + mt2/MW2) -
477 mb2/(3.0*MW2) - (2.0*mt2)/(3.0*
478 MW2))*(-(
log (mb2/mt2)*
log (mt2/MW2))/2.0 +
479 PolyLog.
Li2 (1.0 - mb2/mt2) - ((-(mb2/MW2) + mt2/MW2)*
481 1.0))/(1.0 + (-(mb2/MW2) + mt2/MW2)*(-(mb2/MW2) +
482 mt2/MW2) - 2.0*(mb2/MW2 + mt2/MW2)))))/(32.0*
483 MW2*(mHp2 - MW2)*pi*pi*SW4)) ;
506 double mHp=
sqrt(mHp2);
523 double Lambda5 = (2.0*M2)/(vev*vev);
532 double eta = M_PI/2.0-(beta-alpha1);
534 double CW =
sqrt(CW2);
535 double CW4 = CW2*CW2;
536 double CW6 = CW4*CW2;
537 double CW8 = CW4*CW4;
538 double CW10 = CW8*CW2;
539 double CW12 = CW10*CW2;
540 double CW14 = CW10*CW4;
543 double MZ4 = MZ2*MZ2;
544 double MZ6 = MZ4*MZ2;
545 double MZ8 = MZ4*MZ4;
546 double MZ10 = MZ8*MZ2;
547 double MZ12 = MZ8*MZ4;
548 double MZ14 = MZ10*MZ4;
549 double MZ16 = MZ8*MZ8;
552 double MH4 = MH2*MH2;
553 double MH6 = MH4*MH2;
554 double MH8 = MH4*MH4;
555 double MH10 = MH8*MH2;
556 double MH12 = MH10*MH2;
557 double MH14 = MH12*MH2;
558 double MH16 = MH10*MH4*MH2;
563 double mH3_8 = mH3_6*
mH3_2;
564 double mH3_10 = mH3_8*
mH3_2;
566 double mHp4 = mHp2*mHp2;
567 double mHp3 = mHp2*
sqrt(mHp2);
568 double mHp6 = mHp4*mHp2;
569 double mHp8 = mHp6*mHp2;
570 double mHp10 = mHp8*mHp2;
575 double mH2_6 = mH2_4*
mH2_2;
576 double mH2_8 = mH2_4*mH2_4;
577 double mH2_10 = mH2_8*
mH2_2;
578 double mH2_12 = mH2_10*
mH2_2;
585 aEWadd = (eta*aem*aem*
mMU*
mMU*(-24.0*(1.0 - 4.0*CW2)*(1.0 - 4.0*CW2)*(1.0- 3.0*CW2 + 2.0*CW4)*MH16*
586 pi*pi - 24.0*(1.0- 4.0*CW2)*(1.0- 4.0*CW2)*(1.0- 3.0*CW2 + 2.0*CW4)*MH14*MZ2*
587 (-5.0 - 8.0*CW2 +
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi +
588 256.0*CW4*MZ16*(5.0 - 12.0*CW2 + 8.0*CW4 +
590 CW6*(-2.0 + 9.0*
negsquareroot((MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4))))*
591 (pi - 4.0*CW2*pi)*(pi - 4.0*CW2*pi) - 24.0*(1.0 - 4.0*CW2)*(1.0 - 4.0*CW2)*(1.0 - 3.0*CW2 + 2*CW4)*MH12*
592 MZ4*(6 + (4.0 + 16.0*CW4 - 3.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4) -
593 8.0*CW2*(-5.0 +
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4)))*pi*pi) -
594 MH10*MZ6*(-24.0 + 68.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4)*pi*pi +
595 12288.0*CW12*(-5.0 +
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi +
597 (MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) -
598 6144.0*CW10*(6.0 + (-16.0 +
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) -
601 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) +
604 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
607 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi)) -
608 2*MH4*MH2*MZ10*(-128.0*CW12*(1824 + (57.0*
negsquareroot(1.0 - 4.0*CW2) -
610 10.0*(24.0 + (1.0 - 2*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) +
611 8.0*CW2*(-12.0 + (-43.0+ 16.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) +
614 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
617 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
620 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) +
623 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi)) +
624 MH8*MZ8*(-20.0*(-30.0 + (3.0 +
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) +
625 12288.0*CW12*(-6.0 + (-4.0 + 3.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*
626 pi*pi) - 128.0*CW10*(-624 + (-768.0 + 3.0*
negsquareroot(1.0 - 4.0*CW2) +
630 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
633 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) +
636 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
639 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi)) -
641 64.0*CW4*(-72.0+ (1.0 + 6.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) -
642 16.0*CW2*(-120.0 + (19.0 + 10.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) +
643 128.0*CW8*(-108.0 + (1.0 + 6.0*
negsquareroot((MH2*(MH2 - 4.0*CW2*MZ2))/
644 (CW4*MZ4)))*pi*pi) - CW6*(-17664.0 +
646 907.0*
negsquareroot((MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi)) -
647 4.0*MH4*MZ12*(-20.0*pi*pi + 9216.0*CW14*(24.0 +
negsquareroot(1.0 - 4.0*CW2)*pi*pi) +
648 8.0*CW2*(-120 + (21.0 + 10.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) -
649 16.0*CW4*(-474.0 + (5.0 + 47.0*
negsquareroot((MH2*(MH2 - 4.0*MZ2))/MZ4))*pi*pi) +
652 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
653 16.0*CW12*(18432.0 + (-416.0 + 1404.0*
negsquareroot(1.0 - 4.0*CW2) -
655 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) +
656 8.0*CW10*(5784.0 + (-1616.0 + 2190.0*
negsquareroot(1.0 - 4.0*CW2) -
658 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi) -
661 (MH2*(MH2 - 4.0*CW2*MZ2))/(CW4*MZ4)))*pi*pi)))*
sl)/
662 (4608.0*CW4*(1.0 - 4.0*CW2)*(1.0 - 4.0*CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MH4*MZ8*(MH2 - MZ2)*
663 (MH2 - 4.0*CW2*MZ2)*(MH2 - 4.0*CW2*MZ2)*pi*pi) +
664 (eta*aem*aem*
mMU*
mMU*((71.0 - 278.0*CW2 + 120.0*CW4)*MH4 +
665 6.0*CW2*(9.0 - 88.0*CW2 + 112.0*CW4)*MH2*MZ2 -
666 128.0*CW4*(10.0 - 49.0*CW2 + 36.0*CW4)*MZ4)*
sl*
log(CW2))/
667 (768.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + 4.0*CW2)*MH2*MZ2*(-MH2 + 4.0*CW2*MZ2)*
668 pi*pi) - (eta*aem*aem*(-3.0 + 4.0*CW2)*
mMU*
mMU*((1.0 - 5.0*CW2 + 10.0*CW4)*MH2 +
669 (-7.0 + 61.0*CW2 - 162.0*CW4 + 96.0*CW6)*MZ2)*
sl*
log(CW2)*
log(CW2))/
670 (256.0*CW2*
negpow((1.0 - 4.0*CW2), 3.0/2.0)*(-1.0 + CW2)*(-1.0 + CW2)*MZ2*(-MH2 + MZ2)*pi*pi) +
671 (eta*aem*aem*(-3.0 + 4.0*CW2)*
mMU*
mMU*((1.0 - 5.0*CW2 + 10.0*CW4)*MH2 +
672 (-7.0 + 61.0*CW2 - 162.0*CW4 + 96.0*CW6)*MZ2)*
sl*
log((1.0 -
negsquareroot(1.0 - 4.0*CW2))/2.0)*
log((1.0 -
negsquareroot(1.0 - 4.0*CW2))/2.0))/(128.0*CW2*
negpow((1.0 - 4.0*CW2),3.0/2.0)*(-1.0 + CW2)*(-1.0 + CW2)*MZ2*(-MH2 + MZ2)*pi*pi)
674 (6.0*(1.0 - 3.0*CW2 + 2.0*CW4)*MH4*MH2 - 12.0*(1.0 - 3.0*CW2 + 2.0*CW4)*MH4*MZ2 +
675 (5.0 - 12.0*CW2 + 8.0*CW4)*MH2*MZ4 + 2.0*(5.0 - 12.0*CW2 + 8.0*CW4)*MZ6)*
678 (192.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*MZ6*pi*pi) +
680 (-16.0*MH10 + (16.0 + 99.0*CW2)*MH8*MZ2 - 3.0*CW2*(15.0 + 38.0*CW2)*MH4*MH2*
681 MZ4 + 2.0*CW4*(-231.0 + 113.0*CW2)*MH4*MZ6 - 2.0*CW6*(-907.0 + 768.0*CW2)*
682 MH2*MZ8 + 1152.0*CW8*(-1.0 + CW2)*MZ10)*
sl*
686 MH2/(CW2*MZ2))/2.0))/(384.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*MH4*MZ2*
687 (MH2 - MZ2)*(MH2 - 4.0*CW2*MZ2)*(MH2 - 4.0*CW2*MZ2)*pi*pi) +
688 (eta*aem*aem*
mMU*
mMU*(-48.0*(1.0 - 3.0*CW2 + 2*CW4)*MH8 +
689 (10.0 + 73.0*CW2 - 512.0*CW4 + 384.0*CW6)*MH4*MH2*MZ2 +
690 (-60.0 + 151.0*CW2 + 258.0*CW4 - 256.0*CW6)*MH4*MZ4 +
691 2*(40 + 24.0*CW2 - 149.0*CW4 + 256.0*CW6)*MH2*MZ6 +
692 64*CW2*(-5.0 + 12.0*CW2 - 28.0*CW4 + 18.0*CW6)*MZ8)*
sl*
693 log(MH2/MZ2))/(768.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*MZ4*(MH2 - MZ2)*
694 (MH2 - 4.0*CW2*MZ2)*pi*pi) - (eta*aem*aem*(-1.0 + 2.0*CW2)*MH4*
mMU*
mMU*
695 (MH2 - 4.0*MZ2)*
sl*
log(MH2/MZ2)*
log(MH2/MZ2))/(64*CW4*(-1.0 + CW2)*MZ8*
696 pi*pi) - (eta*aem*aem*(-3.0 + 4.0*CW2)*
mMU*
mMU*((1.0 - 5.0*CW2 + 10.0*CW4)*MH2 +
697 (-7.0 + 61.0*CW2 - 162.0*CW4 + 96.0*CW6)*MZ2)*
sl*
700 (-1.0 + CW2)*MZ2*(-MH2 + MZ2)*pi*pi) -
701 (eta*aem*aem*
mMU*
mMU*(12.0*(1.0 - 3.0*CW2 + 2.0*CW4)*MH10 -
702 48.0*(1.0 - 3.0*CW2 + 2*CW4)*MH8*MZ2 + (17.0 - 48.0*CW2 + 32.0*CW4)*MH4*MH2*
703 MZ4 - 3.0*(5.0 - 12.0*CW2 + 8.0*CW4)*MH2*MZ8 - 4.0*(5.0 - 12.0*CW2 + 8.0*CW4)*
704 MZ10)*
sl*
PolyLog.
Li2(1.0 - MH2/MZ2))/(192.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH4*
706 (6.0*(1.0 - 3.0*CW2 + 2.0*CW4)*MH4*MH2 - 12.0*(1.0 - 3.0*CW2 + 2.0*CW4)*MH4*MZ2 +
707 (5.0 - 12.0*CW2 + 8.0*CW4)*MH2*MZ4 + 2.0*(5.0 - 12.0*CW2 + 8.0*CW4)*MZ6)*
709 (192.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*MZ6*pi*pi) +
711 (6.0*(1.0 - 3.0*CW2 + 2*CW4)*MH4*MH2 - 12.0*(1.0 - 3.0*CW2 + 2.0*CW4)*MH4*MZ2 +
712 (5.0 - 12.0*CW2 + 8.0*CW4)*MH2*MZ4 + 2.0*(5.0 - 12.0*CW2 + 8.0*CW4)*MZ6)*
714 (192.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*MZ6*pi*pi) +
715 (eta*aem*aem*
mMU*
mMU*(-16.0*MH4*MH2 + 3.0*CW2*MH4*MZ2 + 12.0*CW4*MH2*MZ4 +
717 (384.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH4*MZ4*pi*pi) +
719 (16.0*MH10 - (16.0 + 99.0*CW2)*MH8*MZ2 + 3.0*CW2*(15.0 + 38.0*CW2)*MH4*MH2*
720 MZ4 + 2*CW4*(231.0 - 113.0*CW2)*MH4*MZ6 + 2.0*CW6*(-907.0 + 768.0*CW2)*
721 MH2*MZ8 - 1152.0*CW8*(-1.0 + CW2)*MZ10)*
sl*
723 MH2/(CW2*MZ2))/2.0))/(384.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*MH4*MZ2*
724 (MH2 - MZ2)*(MH2 - 4.0*CW2*MZ2)*(MH2 - 4.0*CW2*MZ2)*pi*pi) +
726 (16.0*MH10 - (16.0 + 99.0*CW2)*MH8*MZ2 + 3.0*CW2*(15.0 + 38.0*CW2)*MH4*MH2*
727 MZ4 + 2.0*CW4*(231.0 - 113.0*CW2)*MH4*MZ6 + 2.0*CW6*(-907.0 + 768.0*CW2)*
728 MH2*MZ8 - 1152.0*CW8*(-1.0+ CW2)*MZ10)*
sl*
730 MH2/(CW2*MZ2))/2.0))/(384.0*CW2*(-1.0+ CW2)*(-1.0+ CW2)*MH4*MZ2*
731 (MH2 - MZ2)*(MH2 - 4.0*CW2*MZ2)*(MH2 - 4.0*CW2*MZ2)*pi*pi);
737 gslpp::complex aNonYuk12 = (aem*aem*
mMU*
mMU*(15.0*(47.0 - 78.0*CW2 + 132.0*CW4 - 256.0*CW6 + 128.0*CW8)*mH3_4 -
738 80.0*CW2*(1.0 + 7.0*CW2 - 16.0*CW4 + 8.0*CW6)*
mH3_2*MZ2 + 24.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MZ4))/
739 (8640.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
mH3_2*MZ4*pi*pi);
741 gslpp::complex aNonYuklim = (aem*aem*(3.0 - 6.0*CW2 + 4.0*CW4)*
mH2_2*
mMU*
mMU)/(192.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MZ4*pi*pi);
747 (1152.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*(-2.0*
mH3_2 + MZ2 -
negsquareroot(-4.0*
mH3_2*MZ2 + MZ4))*pi*pi)
751 (1152.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*(-2.0*
mH3_2 + MZ2 -
negsquareroot(-4.0*
mH3_2*MZ2 + MZ4))*pi*pi)
752 -((aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
mMU*
mMU*(-4.0*
mH3_2 + CW2*MZ2)*
756 (1152.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*(-2.0*
mH3_2 + CW2*MZ2 -
761 (1152.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*(2.0*
mH3_2 - CW2*MZ2 +
770 aNonYuk = aNonYuk12+aNonYuklim+aNonYuk6;
775 aNonYuk12 = (aem*aem*
mMU*
776 mMU*(-15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*
mH3_2*
778 15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*mH2_4*
780 mHp2*((50.0 - 64.0*CW2 + 8.0*CW4)*mHp2 +
781 CW2*(13.0 - 43.0*CW2 + 50.0*CW4 - 20.0*CW6)*MZ2) +
782 2.0*(15.0*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp6 +
783 30.0*CW2*(-17.0 + 41.0*CW2 - 80.0*CW4 + 184.0*CW6 -
784 192.0*CW8 + 64.0*CW10)*mHp4*MZ2 -
785 80.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 - 8.0*CW2 +
787 24.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6) +
789 mHp2*(-50.0*mHp2 + CW2*(64.0*mHp2 - 13.0*MZ2) +
790 4.0*CW8*(4.0*
mH2_2 + 5.0*MZ2) -
791 2.0*CW6*(16.0*
mH2_2 + 25.0*MZ2) +
792 CW4*(10.0*
mH2_2 - 8.0*mHp2 + 43.0*MZ2))))/(17280.0*
793 CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*MZ6*pi*pi) +
796 2.0*CW2*mHp2*MZ2 +
mH3_2*(-4.0*mHp2 + CW2*MZ2) +
797 mH2_2*(-4.0*mHp2 + CW2*MZ2)))/(128.0*
798 CW6*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*pi*pi) - (aem*
801 log(
mH3_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
802 pi*pi) - (aem*aem*(
mH2_2 - mHp2)*
mMU*
803 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
804 log(
mH2_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
807 aNonYuk3 = (aem*aem*
mMU*
808 mMU*((-7.0 + 14.0*CW2 - 4.0*CW4 + 5.0*CW6 - 16.0*CW8 + 8.0*CW10)*
810 mH3_8*((7.0 - 14.0*CW2 + 4.0*CW4 - 20.0*CW6 + 64.0*CW8 -
812 mH2_2 + (28.0 - 56.0*CW2 + 16.0*CW4 - 5.0*CW6 +
813 16.0*CW8 - 8.0*CW10)*mHp2 -
814 3.0*CW2*(-19.0 + 26.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*MZ2) +
815 mH3_6*(-14.0*mHp2*(2.0*
mH2_2 + 3.0*mHp2) +
816 CW2*(84.0*mHp4 - 93.0*mHp2*MZ2 +
817 mH2_2*(56.0*mHp2 - 57*MZ2)) +
818 8.0*CW10*(6.0*mH2_4 + 3.0*MZ2*(mHp2 + MZ2) +
819 mH2_2*(4.0*mHp2 + 3.0*MZ2)) -
820 4.0*CW8*(24.0*mH2_4 + 4.0*
mH2_2*(4.0*mHp2 + 3.0*MZ2) +
821 3.0*MZ2*(4.0*mHp2 + 5.0*MZ2)) +
822 CW6*(30.0*mH2_4 +
mH2_2*(20.0*mHp2 + 3.0*MZ2) +
823 3.0*MZ2*(mHp2 + 13.0*MZ2)) +
824 CW4*(
mH2_2*(-16.0*mHp2 + 78.0*MZ2) -
825 3.0*(8.0*mHp4 - 38*mHp2*MZ2 + MZ4))) +
827 4.0*CW10*(8.0*mH2_6 + 6.0*
mH2_2*mHp2*MZ2 +
828 6.0*mHp2*MZ4 - MZ6 + 6.0*mH2_4*(2.0*mHp2 - MZ2)) -
829 CW6*(20.0*mH2_6 + 12*mHp4*MZ2 + 69.0*mHp2*MZ4 +
830 2.0*MZ6 + 3.0*
mH2_2*MZ2*(mHp2 - 10.0*MZ2) +
831 15.0*mH2_4*(2.0*mHp2 - MZ2)) +
832 CW2*(-56.0*mHp6 + 15.0*mHp4*MZ2 +
833 mH2_2*(-84.0*mHp4 + 93.0*mHp2*MZ2)) +
834 2.0*CW8*(32.0*mH2_6 + 36*mHp2*MZ4 - MZ6 +
835 24.0*mH2_4*(2.0*mHp2 - MZ2) +
836 6.0*
mH2_2*(4.0*mHp2*MZ2 - MZ4)) +
837 CW4*(3.0*
mH2_2*(8.0*mHp4 - 38*mHp2*MZ2 - 3.0*MZ4) +
838 2.0*mHp2*(8.0*mHp4 + 3.0*mHp2*MZ2 + 6.0*MZ4))) +
840 2.0*CW8*(8.0*mH2_8 + 6.0*mHp4*MZ4 -
mH2_2*MZ6 -
841 mHp2*MZ6 + 8*mH2_6*(4.0*mHp2 - 3.0*MZ2) +
842 24.0*mH2_4*MZ2*(-mHp2 + MZ2)) -
843 CW4*mHp4*(4.0*mHp4 + 42*mHp2*MZ2 + 21.0*MZ4 +
844 2.0*
mH2_2*(8.0*mHp2 + 3.0*MZ2)) +
845 CW2*(7.0*mHp6*(2.0*mHp2 + 3.0*MZ2) +
846 mH2_2*(56.0*mHp6 - 15.0*mHp4*MZ2)) +
847 4.0*CW10*(2.0*mH2_8 -
mH2_2*MZ6 - mHp2*MZ6 +
848 mH2_6*(8.0*mHp2 - 6.0*MZ2) +
849 mH2_4*(-6.0*mHp2*MZ2 + 6.0*MZ4)) +
850 CW6*(5.0*mH2_8 + 5.0*mH2_6*(4.0*mHp2 - 3.0*MZ2) +
851 15.0*mH2_4*MZ2*(-mHp2 + MZ2) +
852 2.0*mHp2*MZ2*(6.0*mHp4 + 21.0*mHp2*MZ2 + MZ4) +
853 2.0*
mH2_2*(6.0*mHp4*MZ2 + MZ6))) +
854 mH2_2*mHp2*(7.0*mHp6 - 7.0*CW2*(2.0*mHp6 + 3.0*mHp4*MZ2) +
855 CW4*(4.0*mHp6 + 42*mHp4*MZ2 + 21.0*mHp2*MZ4) -
856 4.0*CW10*(2.0*mH2_6 - 6.0*mH2_4*MZ2 + 6.0*
mH2_2*MZ4 -
858 CW8*(8.0*mH2_6 - 24.0*mH2_4*MZ2 + 24.0*
mH2_2*MZ4 +
859 6.0*mHp2*MZ4 - MZ6) -
860 CW6*(5.0*mH2_6 - 15.0*mH2_4*MZ2 + 15.0*
mH2_2*MZ4 +
861 2.0*(6.0*mHp4*MZ2 + 21.0*mHp2*MZ4 + MZ6))))*
863 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 -
866 aNonYuk6 = (aem*aem*(
mH3_2 - mHp2)*
869 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
871 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
873 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
874 aem*(1.0 - 2.0*CW2)*(1.0 - 2.0*CW2)*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
878 mHp2)))/(1152.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
879 MZ6*(-2.0*mHp2 + MZ2 -
negsquareroot (-4.0*mHp2*MZ2 + MZ4))*pi*pi) - (aem*
880 aem*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
892 MZ2)))/(2.0*mH3*mH2)))/(1152.0*
893 CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
897 pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
mMU*
899 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
901 mH3_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
903 2.0*
mH3_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
906 2.0*
mH3_2*(mHp2 + CW2*MZ2))))*
909 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
910 CW2*MZ2)))/(2.0*mH3*mHp)))/(1152.0*
911 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
912 MZ8*(
mH3_2 + mHp2 - CW2*MZ2 +
914 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
915 CW2*MZ2)))*pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
917 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
919 mH2_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
921 2.0*
mH2_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
924 2.0*
mH2_2*(mHp2 + CW2*MZ2))))*
927 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
928 CW2*MZ2)))/(2.0*mH2*mHp)))/(1152.0*
929 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
930 MZ8*(
mH2_2 + mHp2 - CW2*MZ2 +
932 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
933 CW2*MZ2)))*pi*pi) - (aem*aem*(
mH3_2 - mHp2)*
mMU*
935 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) - (aem*
937 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
939 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
943 CW2*mHp*MZ2)*(mHp3 - CW2*mHp*MZ2) +
944 mH3_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
TF (mH3, mHp, CW*MZ))/(64.0*
945 CW8*(-1.0 + CW2)*(-1.0 + CW2)*
946 MZ8*(
mH3_2*
mH3_2 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
947 2.0*
mH3_2*(mHp2 + CW2*MZ2))*pi*pi) + (aem*aem*(
mH2_2 - mHp2)*
949 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
950 CW2*mHp*MZ2) +
mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
951 TF(mH2, mHp, CW*MZ))/(64.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*
952 MZ8*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
953 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi);
957 ((7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mH3_8 +
958 mH3_6*(-4.0*(7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*
mH2_2 -
959 3.0*CW2*(7.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*MZ2) +
960 3.0*mH3_4*(2.0*(7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mH2_4 +
961 CW2*(-5.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*
mH2_2*MZ2 + CW4*(7.0 - 12.0*CW2 + 8.0*CW4)*
962 MZ4) +
mH2_2*((7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mH2_6 -
963 3.0*CW2*(19.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*mH2_4*MZ2 +
964 3.0*CW4*(13.0 - 12.0*CW2 + 8.0*CW4)*
mH2_2*MZ4 - 2.0*CW6*(1.0+ 2.0*CW2)*MZ6) +
965 mH3_2*(-4.0*(7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mH2_6 +
966 3.0*CW2*(31.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*mH2_4*MZ2 - 36.0*CW4*
mH2_2*MZ4 +
967 2.0*CW6*(1.0 + 2.0*CW2)*MZ6))*
log(
mH2_2/MZ2)))/(2304.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*
973 aNonYuk = aNonYuk12+aNonYuk3+ aNonYuklim +aNonYuk6;
976 else if (
mH3_2 == mHp2){
978 aNonYuk12 = (aem*aem*
mMU*
979 mMU*(-15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*
mH3_2*
981 15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*mH2_4*
983 mHp2*((50.0 - 64.0*CW2 + 8.0*CW4)*mHp2 +
984 CW2*(13.0 - 43.0*CW2 + 50.0*CW4 - 20.0*CW6)*MZ2) +
985 2.0*(15.0*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp6 +
986 30.0*CW2*(-17.0 + 41.0*CW2 - 80.0*CW4 + 184.0*CW6 -
987 192.0*CW8 + 64.0*CW10)*mHp4*MZ2 -
988 80.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 - 8.0*CW2 +
990 24.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6) +
992 mHp2*(-50.0*mHp2 + CW2*(64.0*mHp2 - 13.0*MZ2) +
993 4.0*CW8*(4.0*
mH2_2 + 5.0*MZ2) -
994 2.0*CW6*(16.0*
mH2_2 + 25.0*MZ2) +
995 CW4*(10.0*
mH2_2 - 8.0*mHp2 + 43.0*MZ2))))/(17280.0*
996 CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*MZ6*pi*pi) +
999 2.0*CW2*mHp2*MZ2 +
mH3_2*(-4.0*mHp2 + CW2*MZ2) +
1000 mH2_2*(-4.0*mHp2 + CW2*MZ2)))/(128.0*
1001 CW6*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*pi*pi) - (aem*
1004 log(
mH3_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
1005 pi*pi) - (aem*aem*(
mH2_2 - mHp2)*
mMU*
1006 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1007 log(
mH2_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
1011 aNonYuk4 = (aem*aem*
mMU*
1015 mHp2)*(14.0*mH2_4 + 7.0*mHp2*(2.0*mHp2 + 3.0*MZ2) +
1016 mH2_2*(-28.0*mHp2 + 57*MZ2)) -
1017 4.0*CW10*(
mH2_2 - mHp2)*(2.0*mH3_8 -
1018 2.0*mH3_6*(4.0*
mH2_2 + 3.0*MZ2) +
1020 mH3_2*(8.0*mH2_6 - 6.0*mH2_4*MZ2 + MZ6) +
1021 mH2_2*(2.0*mH2_6 - 6.0*mH2_4*MZ2 + 6.0*
mH2_2*MZ4 +
1023 CW4*(
mH3_2*(-4.0*mH2_8 + 2.0*mH2_6*(8.0*mHp2 - 39*MZ2) +
1024 2.0*
mH2_2*(8.0*mHp6 + 3.0*mHp4*MZ2) +
1025 mH2_4*(-24.0*mHp4 + 114*mHp2*MZ2 + 9*MZ4) -
1026 mHp4*(4.0*mHp4 + 42*mHp2*MZ2 + 21.0*MZ4)) +
1027 mH2_2*(4.0*mH2_8 + 4.0*mHp8 + 42*mHp6*MZ2 +
1028 21.0*mHp4*MZ4 + mH2_6*(-16.0*mHp2 + 78.0*MZ2) +
1029 3.0*mH2_4*(8.0*mHp4 - 38*mHp2*MZ2 + MZ4) -
1030 2.0*
mH2_2*(8.0*mHp6 + 3.0*mHp4*MZ2 + 6.0*mHp2*MZ4))) +
1031 2.0*CW8*(8.0*mH3_8*(
mH2_2 - mHp2) -
1032 8*mH3_6*(
mH2_2 - mHp2)*(4.0*
mH2_2 + 3.0*MZ2) +
1035 mH3_2*(-32.0*mH2_8 - 6.0*mHp4*MZ4 -
mH2_2*MZ6 +
1036 mHp2*MZ6 + 6.0*mH2_4*MZ2*(-4.0*mHp2 + MZ2) +
1037 8*mH2_6*(4.0*mHp2 + 3.0*MZ2)) +
1038 mH2_2*(8.0*mH2_8 + 6.0*mHp4*MZ4 - mHp2*MZ6 -
1039 8*mH2_6*(mHp2 + 3.0*MZ2) +
1040 6.0*mH2_4*(4.0*mHp2*MZ2 + 5.0*MZ4) +
1041 mH2_2*(-36*mHp2*MZ4 + MZ6))) +
1042 CW6*(-5.0*mH3_8*(
mH2_2 - mHp2) +
1043 5.0*mH3_6*(
mH2_2 - mHp2)*(4.0*
mH2_2 + 3.0*MZ2) -
1046 mH3_2*(20.0*mH2_8 + 3.0*mH2_4*MZ2*(mHp2 - 10.0*MZ2) -
1047 mH2_6*(20.0*mHp2 + 3.0*MZ2) +
1048 2.0*mHp2*MZ2*(6.0*mHp4 + 21.0*mHp2*MZ2 + MZ4) -
1049 2.0*
mH2_2*(6.0*mHp4*MZ2 + MZ6)) +
1050 mH2_2*(-5.0*mH2_8 + mH2_6*(5.0*mHp2 + 3.0*MZ2) -
1051 3.0*mH2_4*MZ2*(mHp2 + 13.0*MZ2) -
1052 2.0*mHp2*MZ2*(6.0*mHp4 + 21.0*mHp2*MZ2 + MZ4) +
1053 mH2_2*(12*mHp4*MZ2 + 69.0*mHp2*MZ4 + 2.0*MZ6))))*
1055 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 -
1058 aNonYuk6 = (aem*aem*(
mH3_2 - mHp2)*
1061 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1063 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1065 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1066 aem*(1.0 - 2.0*CW2)*(1.0 - 2.0*CW2)*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
1070 mHp2)))/(1152.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1071 MZ6*(-2.0*mHp2 + MZ2 -
negsquareroot (-4.0*mHp2*MZ2 + MZ4))*pi*pi) - (aem*
1072 aem*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
1084 MZ2)))/(2.0*mH3*mH2)))/(1152.0*
1085 CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1089 pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
mMU*
1091 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1093 mH3_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
1095 2.0*
mH3_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
1098 2.0*
mH3_2*(mHp2 + CW2*MZ2))))*
1101 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1102 CW2*MZ2)))/(2.0*mH3*mHp)))/(1152.0*
1103 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1104 MZ8*(
mH3_2 + mHp2 - CW2*MZ2 +
1106 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1107 CW2*MZ2)))*pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
1109 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1111 mH2_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
1113 2.0*
mH2_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
1116 2.0*
mH2_2*(mHp2 + CW2*MZ2))))*
1119 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1120 CW2*MZ2)))/(2.0*mH2*mHp)))/(1152.0*
1121 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1122 MZ8*(
mH2_2 + mHp2 - CW2*MZ2 +
1124 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1125 CW2*MZ2)))*pi*pi) - (aem*aem*(
mH3_2 - mHp2)*
mMU*
1127 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) - (aem*
1129 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1131 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1135 CW2*mHp*MZ2)*(mHp3 - CW2*mHp*MZ2) +
1136 mH3_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
TF (mH3, mHp, CW*MZ))/(64.0*
1137 CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1138 MZ8*(
mH3_2*
mH3_2 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1139 2.0*
mH3_2*(mHp2 + CW2*MZ2))*pi*pi) + (aem*aem*(
mH2_2 - mHp2)*
1141 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
1142 CW2*mHp*MZ2) +
mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
1143 TF(mH2, mHp, CW*MZ))/(64.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1144 MZ8*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1145 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi);
1147 aNonYuklim = -((aem*aem*
mMU*
mMU*(12*CW4*mHp2*(-
mH2_2 + mHp2)*MZ4 +
1148 ((7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mH2_8 +
1149 mH2_6*(-4.0*(7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mHp2 -
1150 3.0*CW2*(7.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*MZ2) +
1151 3.0*mH2_4*(2.0*(7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mHp4 +
1152 CW2*(-5.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*mHp2*MZ2 + CW4*(7.0 - 12.0*CW2 + 8.0*CW4)*
1153 MZ4) + mHp2*((7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mHp6 -
1154 3.0*CW2*(19.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*mHp4*MZ2 +
1155 3.0*CW4*(13.0 - 12.0*CW2 + 8.0*CW4)*mHp2*MZ4 - 2.0*CW6*(1.0 + 2.0*CW2)*MZ6) +
1156 mH2_2*(-4.0*(7.0 - 7.0*CW2 - 3.0*CW4 - 8.0*CW6 + 8.0*CW8)*mHp6 +
1157 3.0*CW2*(31.0 - 7.0*CW2 - 8.0*CW4 + 8.0*CW6)*mHp4*MZ2 - 36.0*CW4*mHp2*MZ4 +
1158 2.0*CW6*(1.0 + 2.0*CW2)*MZ6))*
log(mHp2/MZ2)))/(2304.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1159 (
mH2_2 - mHp2)*MZ8*pi*pi));
1163 aNonYuk = aNonYuk12+aNonYuk4+aNonYuklim+aNonYuk6;
1169 aNonYuk12 = (aem*aem*
mMU*
1170 mMU*(-15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*
mH3_2*
1172 15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*mH2_4*
1174 mHp2*((50.0 - 64.0*CW2 + 8.0*CW4)*mHp2 +
1175 CW2*(13.0 - 43.0*CW2 + 50.0*CW4 - 20.0*CW6)*MZ2) +
1176 2.0*(15.0*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp6 +
1177 30.0*CW2*(-17.0 + 41.0*CW2 - 80.0*CW4 + 184.0*CW6 -
1178 192.0*CW8 + 64.0*CW10)*mHp4*MZ2 -
1179 80.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 - 8.0*CW2 +
1181 24.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6) +
1183 mHp2*(-50.0*mHp2 + CW2*(64.0*mHp2 - 13.0*MZ2) +
1184 4.0*CW8*(4.0*
mH2_2 + 5.0*MZ2) -
1185 2.0*CW6*(16.0*
mH2_2 + 25.0*MZ2) +
1186 CW4*(10.0*
mH2_2 - 8.0*mHp2 + 43.0*MZ2))))/(17280.0*
1187 CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*MZ6*pi*pi) +
1190 2.0*CW2*mHp2*MZ2 +
mH3_2*(-4.0*mHp2 + CW2*MZ2) +
1191 mH2_2*(-4.0*mHp2 + CW2*MZ2)))/(128.0*
1192 CW6*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*pi*pi) - (aem*
1195 log(
mH3_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
1196 pi*pi) - (aem*aem*(
mH2_2 - mHp2)*
mMU*
1197 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1198 log(
mH2_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
1201 aNonYuk5 = (aem*aem*
mMU*
1202 mMU*(-((7.0 - 14.0*CW2 + 4.0*CW4)*
1203 mH3_8*(
mH2_2 - mHp2)) + (7.0 - 14.0*CW2 + 4.0*CW4)*
1204 mH3_6*(
mH2_2 - mHp2)*(4.0*mHp2 + 3.0*CW2*MZ2) -
1206 mH3_2*(
mH2_2 - mHp2)*(2.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4 +
1207 CW2*(-5.0 - 2.0*CW2 + 4.0*CW4)*mHp2*MZ2 +
1208 CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*MZ4) +
1209 mHp2*((7.0 - 14.0*CW2 + 4.0*CW4)*
1210 mH2_8 - (7.0 - 14.0*CW2 + 4.0*CW4)*
1211 mH2_6*(4.0*mHp2 + 3.0*CW2*MZ2) +
1212 3.0*mH2_4*(2.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4 +
1213 CW2*(-5.0 - 2.0*CW2 + 4.0*CW4)*mHp2*MZ2 +
1214 CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*MZ4) +
1215 mH2_2*(-5.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 +
1216 6.0*CW2*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp4*MZ2 -
1217 3.0*CW4*(29 - 50.0*CW2 + 12*CW4)*mHp2*MZ4 +
1218 4.0*CW6*(1.0 + CW2 - 2.0*CW4)*MZ6) +
1219 2.0*mHp2*((7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 -
1220 3.0*CW2*(19.0 - 26.0*CW2 + 4.0*CW4)*mHp4*MZ2 +
1221 3.0*CW4*(13.0 - 20.0*CW2 + 4.0*CW4)*mHp2*MZ4 +
1222 2.0*CW6*(-1.0 - CW2 + 2.0*CW4)*MZ6)) +
1223 mH3_2*((-7.0 + 14.0*CW2 - 4.0*CW4)*
1224 mH2_8 + (7.0 - 14.0*CW2 + 4.0*CW4)*
1225 mH2_6*(4.0*mHp2 + 3.0*CW2*MZ2) -
1226 3.0*mH2_4*(2.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4 +
1227 CW2*(-5.0 - 2.0*CW2 + 4.0*CW4)*mHp2*MZ2 +
1228 CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*MZ4) +
1229 mHp2*(-5.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 +
1230 6.0*CW2*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp4*MZ2 -
1231 3.0*CW4*(29 - 50.0*CW2 + 12*CW4)*mHp2*MZ4 +
1232 4.0*CW6*(1.0 + CW2 - 2.0*CW4)*MZ6) +
1233 2.0*
mH2_2*(4.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 -
1234 3.0*CW2*(31 - 38*CW2 + 4.0*CW4)*mHp4*MZ2 +
1235 6.0*CW4*(8 - 15.0*CW2 + 4.0*CW4)*mHp2*MZ4 +
1236 2.0*CW6*(-1.0 - CW2 + 2.0*CW4)*MZ6)))*
1237 log(mHp2/MZ2))/(2304.0*
1238 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 -
1239 mHp2)*(-
mH2_2 + mHp2)*MZ8*pi*pi);
1241 aNonYuk6 = (aem*aem*(
mH3_2 - mHp2)*
1244 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1246 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1248 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1249 aem*(1.0 - 2.0*CW2)*(1.0 - 2.0*CW2)*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
1253 mHp2)))/(1152.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1254 MZ6*(-2.0*mHp2 + MZ2 -
negsquareroot (-4.0*mHp2*MZ2 + MZ4))*pi*pi) - (aem*
1255 aem*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
1267 MZ2)))/(2.0*mH3*mH2)))/(1152.0*
1268 CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1272 pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
mMU*
1274 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1276 mH3_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
1278 2.0*
mH3_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
1281 2.0*
mH3_2*(mHp2 + CW2*MZ2))))*
1284 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1285 CW2*MZ2)))/(2.0*mH3*mHp)))/(1152.0*
1286 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1287 MZ8*(
mH3_2 + mHp2 - CW2*MZ2 +
1289 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1290 CW2*MZ2)))*pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
1292 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1294 mH2_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
1296 2.0*
mH2_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
1299 2.0*
mH2_2*(mHp2 + CW2*MZ2))))*
1302 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1303 CW2*MZ2)))/(2.0*mH2*mHp)))/(1152.0*
1304 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1305 MZ8*(
mH2_2 + mHp2 - CW2*MZ2 +
1307 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1308 CW2*MZ2)))*pi*pi) - (aem*aem*(
mH3_2 - mHp2)*
mMU*
1310 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) - (aem*
1312 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1314 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1318 CW2*mHp*MZ2)*(mHp3 - CW2*mHp*MZ2) +
1319 mH3_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
TF (mH3, mHp, CW*MZ))/(64.0*
1320 CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1321 MZ8*(
mH3_2*
mH3_2 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1322 2.0*
mH3_2*(mHp2 + CW2*MZ2))*pi*pi) + (aem*aem*(
mH2_2 - mHp2)*
1324 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
1325 CW2*mHp*MZ2) +
mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
1326 TF(mH2, mHp, CW*MZ))/(64.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1327 MZ8*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1328 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi);
1331 aNonYuklim = (aem*aem*
mMU*
mMU*(6.0*CW4*(-1.0 + 7.0*CW2 - 10.0*CW4 + 4.0*CW6)*
mH3_2*(
mH3_2 - mHp2)*MZ4 +
1332 ((-7.0 + 14.0*CW2 - 4.0*CW4)*mH3_8 + mH3_6*(4.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp2 +
1333 3.0*CW2*(19.0 - 26*CW2 + 4*CW4)*MZ2) - 3*mH3_4*(2*(7 - 14*CW2 + 4*CW4)*mHp4 +
1334 CW2*(31.0 - 38.0*CW2 + 4.0*CW4)*mHp2*MZ2 + CW4*(1.0 - 8.0*CW2 + 4.0*CW4)*MZ4) +
1335 mHp2*((-7.0 + 14.0*CW2 - 4.0*CW4)*mHp6 + 3.0*CW2*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4*MZ2 -
1336 3*CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp2*MZ4 + 2.0*CW6*(1.0 + CW2 - 2.0*CW4)*MZ6) +
1337 mH3_2*(4.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 + 3.0*CW2*(5.0 + 2*CW2 - 4.0*CW4)*mHp4*MZ2 +
1338 6.0*CW4*(2.0 - 9.0*CW2 + 4.0*CW4)*mHp2*MZ4 + 2.0*CW6*(-1.0 - CW2 + 2.0*CW4)*MZ6))*
1339 log(
mH3_2/MZ2)))/(1152.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 - mHp2)*MZ8*pi*pi);
1343 aNonYuk = aNonYuk12+ aNonYuklim +aNonYuk5+aNonYuk6;
1349 aNonYuk = (aem*aem*
mMU*
1350 mMU*(-15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*
mH3_2*
1352 15.0*(-25.0 + 32.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*mH2_4*
1354 mHp2*((50.0 - 64.0*CW2 + 8.0*CW4)*mHp2 +
1355 CW2*(13.0 - 43.0*CW2 + 50.0*CW4 - 20.0*CW6)*MZ2) +
1356 2.0*(15.0*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp6 +
1357 30.0*CW2*(-17.0 + 41.0*CW2 - 80.0*CW4 + 184.0*CW6 -
1358 192.0*CW8 + 64.0*CW10)*mHp4*MZ2 -
1359 80.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 - 8.0*CW2 +
1361 24.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*MZ6) +
1363 mHp2*(-50.0*mHp2 + CW2*(64.0*mHp2 - 13.0*MZ2) +
1364 4.0*CW8*(4.0*
mH2_2 + 5.0*MZ2) -
1365 2.0*CW6*(16.0*
mH2_2 + 25.0*MZ2) +
1366 CW4*(10.0*
mH2_2 - 8.0*mHp2 + 43.0*MZ2))))/(17280.0*
1367 CW6*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*MZ6*pi*pi) +
1370 2.0*CW2*mHp2*MZ2 +
mH3_2*(-4.0*mHp2 + CW2*MZ2) +
1371 mH2_2*(-4.0*mHp2 + CW2*MZ2)))/(128.0*
1372 CW6*(-1.0 + CW2)*(-1.0 + CW2)*MZ6*pi*pi) - (aem*
1375 log(
mH3_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
1376 pi*pi) - (aem*aem*(
mH2_2 - mHp2)*
mMU*
1377 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1378 log(
mH2_2/mHp2))/(128.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*
1379 pi*pi)) + (aem*aem*
mMU*
1380 mMU*((-7.0 + 14.0*CW2 - 4.0*CW4 + 5.0*CW6 - 16.0*CW8 + 8.0*CW10)*
1382 mH3_8*((7.0 - 14.0*CW2 + 4.0*CW4 - 20.0*CW6 + 64.0*CW8 -
1384 mH2_2 + (28.0 - 56.0*CW2 + 16.0*CW4 - 5.0*CW6 +
1385 16.0*CW8 - 8.0*CW10)*mHp2 -
1386 3.0*CW2*(-19.0 + 26.0*CW2 + CW4 - 16.0*CW6 + 8.0*CW8)*MZ2) +
1387 mH3_6*(-14.0*mHp2*(2.0*
mH2_2 + 3.0*mHp2) +
1388 CW2*(84.0*mHp4 - 93.0*mHp2*MZ2 +
1389 mH2_2*(56.0*mHp2 - 57*MZ2)) +
1390 8.0*CW10*(6.0*mH2_4 + 3.0*MZ2*(mHp2 + MZ2) +
1391 mH2_2*(4.0*mHp2 + 3.0*MZ2)) -
1392 4.0*CW8*(24.0*mH2_4 + 4.0*
mH2_2*(4.0*mHp2 + 3.0*MZ2) +
1393 3.0*MZ2*(4.0*mHp2 + 5.0*MZ2)) +
1394 CW6*(30.0*mH2_4 +
mH2_2*(20.0*mHp2 + 3.0*MZ2) +
1395 3.0*MZ2*(mHp2 + 13.0*MZ2)) +
1396 CW4*(
mH2_2*(-16.0*mHp2 + 78.0*MZ2) -
1397 3.0*(8.0*mHp4 - 38*mHp2*MZ2 + MZ4))) +
1399 4.0*CW10*(8.0*mH2_6 + 6.0*
mH2_2*mHp2*MZ2 +
1400 6.0*mHp2*MZ4 - MZ6 + 6.0*mH2_4*(2.0*mHp2 - MZ2)) -
1401 CW6*(20.0*mH2_6 + 12*mHp4*MZ2 + 69.0*mHp2*MZ4 +
1402 2.0*MZ6 + 3.0*
mH2_2*MZ2*(mHp2 - 10.0*MZ2) +
1403 15.0*mH2_4*(2.0*mHp2 - MZ2)) +
1404 CW2*(-56.0*mHp6 + 15.0*mHp4*MZ2 +
1405 mH2_2*(-84.0*mHp4 + 93.0*mHp2*MZ2)) +
1406 2.0*CW8*(32.0*mH2_6 + 36*mHp2*MZ4 - MZ6 +
1407 24.0*mH2_4*(2.0*mHp2 - MZ2) +
1408 6.0*
mH2_2*(4.0*mHp2*MZ2 - MZ4)) +
1409 CW4*(3.0*
mH2_2*(8.0*mHp4 - 38*mHp2*MZ2 - 3.0*MZ4) +
1410 2.0*mHp2*(8.0*mHp4 + 3.0*mHp2*MZ2 + 6.0*MZ4))) +
1412 2.0*CW8*(8.0*mH2_8 + 6.0*mHp4*MZ4 -
mH2_2*MZ6 -
1413 mHp2*MZ6 + 8*mH2_6*(4.0*mHp2 - 3.0*MZ2) +
1414 24.0*mH2_4*MZ2*(-mHp2 + MZ2)) -
1415 CW4*mHp4*(4.0*mHp4 + 42*mHp2*MZ2 + 21.0*MZ4 +
1416 2.0*
mH2_2*(8.0*mHp2 + 3.0*MZ2)) +
1417 CW2*(7.0*mHp6*(2.0*mHp2 + 3.0*MZ2) +
1418 mH2_2*(56.0*mHp6 - 15.0*mHp4*MZ2)) +
1419 4.0*CW10*(2.0*mH2_8 -
mH2_2*MZ6 - mHp2*MZ6 +
1420 mH2_6*(8.0*mHp2 - 6.0*MZ2) +
1421 mH2_4*(-6.0*mHp2*MZ2 + 6.0*MZ4)) +
1422 CW6*(5.0*mH2_8 + 5.0*mH2_6*(4.0*mHp2 - 3.0*MZ2) +
1423 15.0*mH2_4*MZ2*(-mHp2 + MZ2) +
1424 2.0*mHp2*MZ2*(6.0*mHp4 + 21.0*mHp2*MZ2 + MZ4) +
1425 2.0*
mH2_2*(6.0*mHp4*MZ2 + MZ6))) +
1426 mH2_2*mHp2*(7.0*mHp6 - 7.0*CW2*(2.0*mHp6 + 3.0*mHp4*MZ2) +
1427 CW4*(4.0*mHp6 + 42*mHp4*MZ2 + 21.0*mHp2*MZ4) -
1428 4.0*CW10*(2.0*mH2_6 - 6.0*mH2_4*MZ2 + 6.0*
mH2_2*MZ4 -
1430 CW8*(8.0*mH2_6 - 24.0*mH2_4*MZ2 + 24.0*
mH2_2*MZ4 +
1431 6.0*mHp2*MZ4 - MZ6) -
1432 CW6*(5.0*mH2_6 - 15.0*mH2_4*MZ2 + 15.0*
mH2_2*MZ4 +
1433 2.0*(6.0*mHp4*MZ2 + 21.0*mHp2*MZ4 + MZ6))))*
1435 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 -
1440 mHp2)*(14.0*mH2_4 + 7.0*mHp2*(2.0*mHp2 + 3.0*MZ2) +
1441 mH2_2*(-28.0*mHp2 + 57*MZ2)) -
1442 4.0*CW10*(
mH2_2 - mHp2)*(2.0*mH3_8 -
1443 2.0*mH3_6*(4.0*
mH2_2 + 3.0*MZ2) +
1445 mH3_2*(8.0*mH2_6 - 6.0*mH2_4*MZ2 + MZ6) +
1446 mH2_2*(2.0*mH2_6 - 6.0*mH2_4*MZ2 + 6.0*
mH2_2*MZ4 +
1448 CW4*(
mH3_2*(-4.0*mH2_8 + 2.0*mH2_6*(8.0*mHp2 - 39*MZ2) +
1449 2.0*
mH2_2*(8.0*mHp6 + 3.0*mHp4*MZ2) +
1450 mH2_4*(-24.0*mHp4 + 114*mHp2*MZ2 + 9*MZ4) -
1451 mHp4*(4.0*mHp4 + 42*mHp2*MZ2 + 21.0*MZ4)) +
1452 mH2_2*(4.0*mH2_8 + 4.0*mHp8 + 42*mHp6*MZ2 +
1453 21.0*mHp4*MZ4 + mH2_6*(-16.0*mHp2 + 78.0*MZ2) +
1454 3.0*mH2_4*(8.0*mHp4 - 38*mHp2*MZ2 + MZ4) -
1455 2.0*
mH2_2*(8.0*mHp6 + 3.0*mHp4*MZ2 + 6.0*mHp2*MZ4))) +
1456 2.0*CW8*(8.0*mH3_8*(
mH2_2 - mHp2) -
1457 8*mH3_6*(
mH2_2 - mHp2)*(4.0*
mH2_2 + 3.0*MZ2) +
1460 mH3_2*(-32.0*mH2_8 - 6.0*mHp4*MZ4 -
mH2_2*MZ6 +
1461 mHp2*MZ6 + 6.0*mH2_4*MZ2*(-4.0*mHp2 + MZ2) +
1462 8*mH2_6*(4.0*mHp2 + 3.0*MZ2)) +
1463 mH2_2*(8.0*mH2_8 + 6.0*mHp4*MZ4 - mHp2*MZ6 -
1464 8*mH2_6*(mHp2 + 3.0*MZ2) +
1465 6.0*mH2_4*(4.0*mHp2*MZ2 + 5.0*MZ4) +
1466 mH2_2*(-36*mHp2*MZ4 + MZ6))) +
1467 CW6*(-5.0*mH3_8*(
mH2_2 - mHp2) +
1468 5.0*mH3_6*(
mH2_2 - mHp2)*(4.0*
mH2_2 + 3.0*MZ2) -
1471 mH3_2*(20.0*mH2_8 + 3.0*mH2_4*MZ2*(mHp2 - 10.0*MZ2) -
1472 mH2_6*(20.0*mHp2 + 3.0*MZ2) +
1473 2.0*mHp2*MZ2*(6.0*mHp4 + 21.0*mHp2*MZ2 + MZ4) -
1474 2.0*
mH2_2*(6.0*mHp4*MZ2 + MZ6)) +
1475 mH2_2*(-5.0*mH2_8 + mH2_6*(5.0*mHp2 + 3.0*MZ2) -
1476 3.0*mH2_4*MZ2*(mHp2 + 13.0*MZ2) -
1477 2.0*mHp2*MZ2*(6.0*mHp4 + 21.0*mHp2*MZ2 + MZ4) +
1478 mH2_2*(12*mHp4*MZ2 + 69.0*mHp2*MZ4 + 2.0*MZ6))))*
1480 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 -
1482 mMU*(-((7.0 - 14.0*CW2 + 4.0*CW4)*
1483 mH3_8*(
mH2_2 - mHp2)) + (7.0 - 14.0*CW2 + 4.0*CW4)*
1484 mH3_6*(
mH2_2 - mHp2)*(4.0*mHp2 + 3.0*CW2*MZ2) -
1486 mH3_2*(
mH2_2 - mHp2)*(2.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4 +
1487 CW2*(-5.0 - 2.0*CW2 + 4.0*CW4)*mHp2*MZ2 +
1488 CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*MZ4) +
1489 mHp2*((7.0 - 14.0*CW2 + 4.0*CW4)*
1490 mH2_8 - (7.0 - 14.0*CW2 + 4.0*CW4)*
1491 mH2_6*(4.0*mHp2 + 3.0*CW2*MZ2) +
1492 3.0*mH2_4*(2.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4 +
1493 CW2*(-5.0 - 2.0*CW2 + 4.0*CW4)*mHp2*MZ2 +
1494 CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*MZ4) +
1495 mH2_2*(-5.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 +
1496 6.0*CW2*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp4*MZ2 -
1497 3.0*CW4*(29 - 50.0*CW2 + 12*CW4)*mHp2*MZ4 +
1498 4.0*CW6*(1.0 + CW2 - 2.0*CW4)*MZ6) +
1499 2.0*mHp2*((7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 -
1500 3.0*CW2*(19.0 - 26.0*CW2 + 4.0*CW4)*mHp4*MZ2 +
1501 3.0*CW4*(13.0 - 20.0*CW2 + 4.0*CW4)*mHp2*MZ4 +
1502 2.0*CW6*(-1.0 - CW2 + 2.0*CW4)*MZ6)) +
1503 mH3_2*((-7.0 + 14.0*CW2 - 4.0*CW4)*
1504 mH2_8 + (7.0 - 14.0*CW2 + 4.0*CW4)*
1505 mH2_6*(4.0*mHp2 + 3.0*CW2*MZ2) -
1506 3.0*mH2_4*(2.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp4 +
1507 CW2*(-5.0 - 2.0*CW2 + 4.0*CW4)*mHp2*MZ2 +
1508 CW4*(7.0 - 14.0*CW2 + 4.0*CW4)*MZ4) +
1509 mHp2*(-5.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 +
1510 6.0*CW2*(25.0 - 32.0*CW2 + 4.0*CW4)*mHp4*MZ2 -
1511 3.0*CW4*(29 - 50.0*CW2 + 12*CW4)*mHp2*MZ4 +
1512 4.0*CW6*(1.0 + CW2 - 2.0*CW4)*MZ6) +
1513 2.0*
mH2_2*(4.0*(7.0 - 14.0*CW2 + 4.0*CW4)*mHp6 -
1514 3.0*CW2*(31 - 38*CW2 + 4.0*CW4)*mHp4*MZ2 +
1515 6.0*CW4*(8 - 15.0*CW2 + 4.0*CW4)*mHp2*MZ4 +
1516 2.0*CW6*(-1.0 - CW2 + 2.0*CW4)*MZ6)))*
1517 log(mHp2/MZ2))/(2304.0*
1518 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*(
mH3_2 -
1519 mHp2)*(-
mH2_2 + mHp2)*MZ8*pi*pi) + (aem*aem*(
mH3_2 - mHp2)*
1522 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1524 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1526 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1527 aem*(1.0 - 2.0*CW2)*(1.0 - 2.0*CW2)*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
1531 mHp2)))/(1152.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1532 MZ6*(-2.0*mHp2 + MZ2 -
negsquareroot (-4.0*mHp2*MZ2 + MZ4))*pi*pi) - (aem*
1533 aem*(5.0 - 16.0*CW2 + 8.0*CW4)*
mMU*
1545 MZ2)))/(2.0*mH3*mH2)))/(1152.0*
1546 CW2*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1550 pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
mMU*
1552 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1554 mH3_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
1556 2.0*
mH3_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
1559 2.0*
mH3_2*(mHp2 + CW2*MZ2))))*
1562 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1563 CW2*MZ2)))/(2.0*mH3*mHp)))/(1152.0*
1564 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1565 MZ8*(
mH3_2 + mHp2 - CW2*MZ2 +
1567 mHp2 + (
mH3_2 + mHp2 - CW2*MZ2)*(
mH3_2 + mHp2 -
1568 CW2*MZ2)))*pi*pi) + (aem*aem*(7.0 - 14.0*CW2 + 4.0*CW4)*
1570 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1572 mH2_2*(-2.0*mHp2 - 2.0*CW2*MZ2 +
1574 2.0*
mH2_2*(mHp2 + CW2*MZ2))) + (mHp2 - CW2*MZ2)*(mHp2 -
1577 2.0*
mH2_2*(mHp2 + CW2*MZ2))))*
1580 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1581 CW2*MZ2)))/(2.0*mH2*mHp)))/(1152.0*
1582 CW8*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + CW2)*
1583 MZ8*(
mH2_2 + mHp2 - CW2*MZ2 +
1585 mHp2 + (
mH2_2 + mHp2 - CW2*MZ2)*(
mH2_2 + mHp2 -
1586 CW2*MZ2)))*pi*pi) - (aem*aem*(
mH3_2 - mHp2)*
mMU*
1588 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) - (aem*
1590 mMU*(mH2_4 - 2.0*
mH2_2*mHp2 + mHp4 - CW2*mHp2*MZ2)*
1592 CW8*(-1.0 + CW2)*(-1.0 + CW2)*MZ8*pi*pi) + (aem*
1596 CW2*mHp*MZ2)*(mHp3 - CW2*mHp*MZ2) +
1597 mH3_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
TF (mH3, mHp, CW*MZ))/(64.0*
1598 CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1599 MZ8*(
mH3_2*
mH3_2 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1600 2.0*
mH3_2*(mHp2 + CW2*MZ2))*pi*pi) + (aem*aem*(
mH2_2 - mHp2)*
1602 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
1603 CW2*mHp*MZ2) +
mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
1604 TF(mH2, mHp, CW*MZ))/(64.0*CW8*(-1.0 + CW2)*(-1.0 + CW2)*
1605 MZ8*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1606 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi);
1886 pi) + (aem*aem*(MH2 + 2.0*mHp2)*
mMU*
1887 mMU*((-3.0 + 2.0*CW2)*MH2 - 8.0*CW2*(-1.0 + CW2)*MZ2)*
1888 log (MH2/MZ2))/(128.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*
1889 MZ2*(MH2 - MZ2)*pi*pi) - (aem*aem*(MH2 + 2.0*mHp2)*
mMU*
mMU*
1890 log (mHp2/MZ2))/(16.0*CW2*(-1.0 + CW2)*MH2*MZ2*pi*pi) + (aem*aem*
1891 mHp2*(MH2 + 2.0*mHp2)*
mMU*
1892 mMU*((-3.0 + 2.0*CW2)*MH2 - 8.0*CW2*(-1.0 + CW2)*MZ2)*
1893 TF (MH, mHp, mHp))/(64.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
1894 MH4*(MH2 - 4.0*mHp2)*MZ2*(MH2 - MZ2)*pi*pi) - (aem*
1895 aem*(3.0 - 10.0*CW2 + 8.0*CW4)*mHp2*(MH2 + 2.0*mHp2)*
mMU*
mMU*
1896 TF (mHp, mHp, MZ))/(64.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
1897 MZ4*(-MH2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi);
1902 8.0*CW2*(-1.0 + CW2)*mHp2*MZ2)*(-1.0 + tanb*tanb))/(64.0*
1903 CW4*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*MZ4*pi*pi*tanb) - (aem*aem*
1907 CW2*(-1.0 + CW2)*(-1.0 + CW2)*
1908 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*
1909 pi*tanb) + (aem*aem*mH2_4*
mMU*
1910 mMU*(-7.0*mHp6 + 8.0*CW2*mHp4*MZ2 - 4.0*CW4*mHp2*MZ4 +
1911 2.0*
mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*(-1.0 +
1913 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
1914 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*
1917 CW2*(-1.0 + CW2)*(-1.0 + CW2)*
1918 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*
1919 pi*tanb) - (aem*aem*mH2_4*
mMU*
1920 mMU*(-7.0*mHp6 + 8.0*CW2*mHp4*MZ2 - 4.0*CW4*mHp2*MZ4 +
1921 2.0*
mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*(-1.0 +
1923 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
1924 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*
1925 pi*tanb) + (aem*aem*
mMU*
1926 mMU*(2.0*CW2*mH2_4*MZ2 +
1927 8.0*CW2*(-1.0 + CW2)*mHp2*MZ2*(-mHp2 + CW2*MZ2) +
1928 mH2_2*(-2.0*mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
1929 log (mHp2/MZ2))/(128.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*
1930 MZ4*(-mHp2 + CW2*MZ2)*pi*pi*tanb) + (aem*aem*
mH2_2*
mMU*
1931 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
1932 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
1934 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ6*pi*pi*tanb) - (aem*aem*
1936 log (mHp2/MZ2))/(256.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*
1937 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*
1940 2.0*mH2_4*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1941 mH2_2*(mHp6 - 8.0*CW2*mHp4*MZ2 + 4.0*CW4*mHp2*MZ4))*(-1.0 +
1942 tanb*tanb)*
log (mHp2/MZ2)*
log (mHp2/MZ2))/(512.0*
1943 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
1944 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*
1948 2.0*mH2_4*MZ2 + (-3.0 + 2.0*CW2)*
mH2_2*mHp2*MZ2 -
1949 8.0*CW2*(-1.0 + CW2)*mHp2*MZ4)*(-1.0 + tanb*tanb))/(128.0*
1950 CW4*(-1.0 + CW2)*(-1.0 + CW2)*mHp2*MZ4*(-
mH2_2 + MZ2)*pi*
1951 pi*tanb) - (3.0*aem*aem*mH2_4*mHp2*
mMU*
1952 mMU*(-1.0 + tanb*tanb)*
log(mHp2/MZ2))/(256.0*
1953 CW2*(-1.0 + CW2)*(-1.0 + CW2)*
1954 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
1956 mMU*(2.0*
mH2_2 - 2.0*mHp2 + CW2*MZ2)*(-1.0 +
1957 tanb*tanb))/(128.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
1958 MZ4*(-mHp2 + CW2*MZ2)*pi*pi*tanb) - (aem*aem*
mH2_2*
mMU*
1959 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
1960 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
1961 log (
mH2_2/mHp2))/(128.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*
1962 mHp4*MZ6*pi*pi*tanb) + (aem*aem*mH2_4*(
mH2_2 - 2.0*mHp2)*
1964 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ2*(-mHp2 + CW2*MZ2)*pi*
1965 pi*tanb) - (aem*aem*mH2_4*(
mH2_2 - 2.0*mHp2)*
mMU*
1966 mMU*(-1.0 + tanb*tanb)*
log (mHp2/MZ2))/(128.0*
1967 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ2*(-mHp2 + CW2*MZ2)*pi*
1969 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
1970 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
1972 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ6*pi*pi*tanb)+ (aem*aem*
1974 2.0*CW4*mHp6*MZ2*(
mH2_2 + 4.0*(mHp2 + MZ2)) +
1975 CW2*(-mH2_10 + 2.0*mH2_4*mHp4*MZ2 + 8.0*mHp8*MZ2 +
1976 mH2_8*(4.0*mHp2 + MZ2) -
1977 2.0*mH2_6*(mHp4 + 2.0*mHp2*MZ2) +
1978 mH2_2*(2.0*mHp8 + 3.0*mHp6*MZ2)))*(-1.0 + tanb*tanb)*
1979 TF (mH2, mHp, mHp))/(64.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
mH2_2*
1980 mHp4*(
mH2_2 - 4.0*mHp2)*MZ2*(
mH2_2 - MZ2)*(-mHp2 + CW2*MZ2)*
1983 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
1984 CW2*mHp*MZ2) +
mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*(-1.0 +
1985 tanb*tanb)*
TF(mH2, mHp, CW*MZ))/(64.0*
1986 CW6*(-1.0 + CW2)*(-1.0 + CW2)*
1988 CW2*MZ2)*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
1989 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi*tanb)+ (aem*
1990 aem*(3.0 - 10.0*CW2 + 8.0*CW4)*
mH2_2*mHp2*
mMU*
1991 mMU*(-1.0 + tanb*tanb)*
TF (mHp, mHp, MZ))/(64.0*
1992 CW4*(-1.0 + CW2)*(-1.0 + CW2)*
1993 MZ4*(-
mH2_2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi*tanb));
1997 mMU*((-3.0 + 2.0*CW2)*MH2 - 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 +
1998 tanb*tanb)*
log(MH2/MZ2))/(128.0*
1999 CW4*(-1.0 + CW2)*(-1.0 + CW2)*MZ2*(-MH2 + MZ2)*pi*pi*
2000 tanb) + (aem*aem*
mMU*
2001 mMU*((3.0 - 2.0*CW2)*
mH2_2 + 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 +
2003 CW4*(-1.0 + CW2)*(-1.0 + CW2)*MZ2*(-
mH2_2 + MZ2)*pi*pi*
2004 tanb) + (aem*aem*mHp2*
mMU*
2005 mMU*((3.0 - 2.0*CW2)*MH2 + 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 +
2006 tanb*tanb)*
TF (MH, mHp, mHp))/(64.0*
2007 CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*(MH2 - 4.0*mHp2)*
2008 MZ2*(MH2 - MZ2)*pi*pi*tanb) + (aem*aem*mHp2*
mMU*
2010 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 + tanb*tanb)*
2011 TF (mH2, mHp, mHp))/(64.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
2013 tanb) + (aem*aem*(3.0 - 10.0*CW2 + 8.0*CW4)*(MH2 -
mH2_2)*
2014 mHp2*
mMU*
mMU*(-1.0 + tanb*tanb)*
TF (mHp, mHp, MZ))/(64.0*
2015 CW4*(-1.0 + CW2)*(-1.0 + CW2)*
2016 MZ2*(-MH2 + MZ2)*(-
mH2_2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi*
2017 tanb)+
sl*((aem*aem*
mMU*
2018 mMU*(144.0*MH6*mH2_4*MZ4 - 144.0*CW2*MH4*mH2_4*MZ6 -
2019 1152.0*CW2*(-1.0 + CW2)*mH2_4*mHp4*MZ6 +
2020 MH2*(3.0*(40.0 - 41.0*CW2 + 64.0*CW4)*mH2_4*mHp2*MZ6 +
2022 MZ4*((1.0 - 3.0*CW2 + 2.0*CW4)*mHp2 + CW2*MZ2) +
2023 24.0*(1.0 - 3.0*CW2 + 2.0*CW4)*mH2_10*mHp2*pi*pi +
2024 16.0*(5.0 - 12.0*CW2 + 8.0*CW4 + 2.0*CW6)*mHp2*MZ10*pi*
2027 2.0*(1.0 - 3.0*CW2 + 2.0*CW4)*mHp2*pi*pi) +
2029 MZ6*(288*CW6*MZ2 + 5.0*MZ2*(-8.0 + pi*pi) -
2030 12.0*CW2*(8.0*mHp2 + MZ2*(-8.0 + pi*pi)) + 2.0*CW4*(48.0*mHp2 +
2031 MZ2*(-184.0 + 5.0*pi*pi))))))/(4608.0*
2032 CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*mH2_4*mHp2*MZ8*pi*
2034 mMU*(MH4 + mH2_4 - 3.0*CW2*mHp2*MZ2 + 3.0*CW4*MZ4 +
2035 MH2*(
mH2_2 - mHp2 - 3.0*CW2*MZ2) -
2036 mH2_2*(mHp2 + 3.0*CW2*MZ2))*
log (CW2)*
log (CW2))/(128.0*
2037 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ2*(-mHp2 + CW2*MZ2)*pi*
2041 CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2042 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2043 pi*pi) + (aem*aem*
mMU*
2044 mMU*(-4.0*MH6*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) +
2046 MZ2*(-3.0*mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2047 CW2*MZ2)*(-mHp2 + CW2*MZ2) +
2048 MH4*(7.0*mHp6 + 2.0*CW2*mHp4*MZ2 - 10.0*CW4*mHp2*MZ4 +
2051 2.0*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2)*(4.0*mHp4 +
2052 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4)))*
2054 log (MH2/mHp2))/(512.0*
2055 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2056 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2057 pi*pi) + (aem*aem*
mMU*
2058 mMU*(4.0*mH2_6*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2060 MZ2*(-3.0*mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2061 CW2*MZ2)*(-mHp2 + CW2*MZ2) -
2062 mH2_4*(9.0*mHp6 + 2.0*CW2*mHp4*MZ2 - 10.0*CW4*mHp2*MZ4 +
2065 2.0*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2)*(4.0*mHp4 +
2066 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4)))*
2069 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2070 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2074 CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2075 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2076 pi*pi) - (aem*aem*
mMU*
2077 mMU*(-4.0*MH6*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) +
2079 MZ2*(-3.0*mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2080 CW2*MZ2)*(-mHp2 + CW2*MZ2) +
2081 MH4*(7.0*mHp6 + 2.0*CW2*mHp4*MZ2 - 10.0*CW4*mHp2*MZ4 +
2084 2.0*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2)*(4.0*mHp4 +
2085 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4)))*
2087 log (MH2/MZ2))/(512.0*
2088 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2089 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2090 pi*pi) + (aem*aem*
mMU*
2092 MZ8*(-3.0*mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2093 CW2*MZ2)*(-mHp2 + CW2*MZ2) -
2094 4.0*mH2_6*(mHp2 - CW2*MZ2)*(mHp2 -
2095 CW2*MZ2)*((2.0 - 6.0*CW2 + 4.0*CW4)*mHp6 -
2096 2.0*CW2*(1.0 - 3.0*CW2 + 2.0*CW4)*mHp4*MZ2 + CW2*MZ6) +
2098 MZ2*(32*(1.0 - 3.0*CW2 + 2.0*CW4)*mHp10 -
2099 96.0*CW2*(1.0 - 3.0*CW2 + 2.0*CW4)*mHp8*MZ2 +
2100 3.0*CW2*(3.0 + 32*CW2 - 96.0*CW4 + 64.0*CW6)*mHp6*
2101 MZ4 + 2.0*CW4*(1.0 - 16.0*CW2 + 48.0*CW4 - 32*CW6)*
2102 mHp4*MZ6 - 10.0*CW6*mHp2*MZ8 + 6.0*CW8*MZ10) -
2105 2.0*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2)*(4.0*mHp4 +
2106 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4)))*
2109 CW4*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2110 MZ8*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2111 pi*pi) + (aem*aem*(MH2 -
mH2_2)*
mMU*
2113 8.0*CW2*(-1.0 + CW2)*mHp4*MZ2*(-mHp2 + CW2*MZ2) +
2115 2.0*CW2*mHp2*MZ2 + 2.0*CW4*MZ4))*
log (mHp2/MZ2))/(64.0*
2116 CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH2*
mH2_2*mHp2*
2117 MZ4*(-mHp2 + CW2*MZ2)*pi*pi) - (aem*aem*
mMU*
2118 mMU*(2.0*mHp8 + MH6*(mHp2 + CW2*MZ2) -
2119 MH2*(3.0*mHp6 + 4.0*CW2*mHp4*MZ2))*
log (MH2/mHp2)*
2120 log (mHp2/MZ2))/(128.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2121 MZ6*pi*pi)+ (aem*aem*
mMU*
2122 mMU*(2.0*mHp8 + mH2_6*(mHp2 + CW2*MZ2) -
2124 log (mHp2/MZ2))/(128.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2127 CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2128 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2129 pi*pi) - (aem*aem*(MH2 -
mH2_2)*
mMU*
2131 2.0*CW2*
mH2_2*mHp4*MZ2 + 22.0*CW2*mHp6*MZ2 -
2132 10.0*CW4*
mH2_2*mHp2*MZ4 - 26*CW4*mHp4*MZ4 +
2133 6.0*CW6*
mH2_2*MZ6 + 18.0*CW6*mHp2*MZ6 - 6.0*CW8*MZ8 +
2134 MH2*(7.0*mHp6 + 2.0*CW2*mHp4*MZ2 - 10.0*CW4*mHp2*MZ4 +
2135 6.0*CW6*MZ6))*
log(mHp2/MZ2)*
log(mHp2/MZ2))/(512.0*
2136 CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2137 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2141 mH2_4*MZ2*((5.0 + 26*CW2 - 16.0*CW4)*mHp2 -
2143 3.0*mH2_6*((3.0 - 12.0*CW2 + 8.0*CW4)*
2144 mHp2 + (2.0 + CW2)*MZ2) +
2146 MZ2*((9.0 - 6.0*CW2)*mHp2 + (15 - 59*CW2 + 62*CW4)*
2149 MZ4*(-5.0*MZ2 + 18.0*CW6*MZ2 +
2150 CW4*(6.0*mHp2 - 28*MZ2) -
2151 6.0*CW2*(mHp2 - 2.0*MZ2))))/(192*
2152 CW4*(-1.0 + CW2)*(-1.0 + CW2)*
mH2_2*mHp2*
2153 MZ4*(
mH2_2 - MZ2)*pi*pi) + (aem*aem*
2155 log (mHp2/MZ2))/(256.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2156 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2159 MZ2)*(-(aem*aem*
mMU*
2160 mMU*(-8.0*MH10- 128.0*CW4*(-1.0 + CW2)*mHp4*MZ6 +
2162 MZ4*(8.0*(-5.0 + 4.0*CW2)*mHp2 +
2163 CW2*(33 - 32*CW2)*MZ2) +
2164 MH6*MZ2*((10.0 + 11.0*CW2)*mHp2 -
2165 4.0*CW2*(9.0 + 4.0*CW2)*MZ2) -
2166 4.0*MH8*(mHp2 - (2.0 + 9.0*CW2)*MZ2) +
2167 MH4*MZ2*((12.0 - 8.0*CW2)*mHp4 +
2168 5.0*CW2*(-11.0 + 6.0*CW2)*mHp2*MZ2 +
2169 16.0*CW4*MZ4)))/(256.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
2170 MH2*mHp2*MZ4*(MH2 - MZ2)*(MH2 - 4.0*CW2*MZ2)*pi*
2171 pi) + (aem*aem*MH2*(5.0*MH2 +
mH2_2)*mHp2*
mMU*
mMU*
2172 log (mHp2/MZ2))/(256.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2173 MZ2*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2176 mMU*(24.0*(1.0 - 4.0*CW2)*MH6*
mH2_2*mHp2 +
2177 12.0*CW2*(-1.0 + 4.0*CW2)*MH4*
mH2_2*
2178 MZ2*(8.0*mHp2 + CW2*MZ2) + 2.0*CW2*
2179 MZ2*(-48.0*(-1.0 + 4.0*CW2)*mH2_6*mHp2 +
2180 24.0*CW4*(-1.0 + 4.0*CW2)*mH2_4*MZ4 +
2181 CW2*(133 - 520*CW2 + 240*CW4)*
mH2_2*mHp2*
2182 MZ2*(-mHp2 + CW2*MZ2) +
2183 64.0*CW4*(10.0 - 49*CW2 + 36*CW4)*mHp2*
2184 MZ4*(-mHp2 + CW2*MZ2)) +
2185 MH2*(24.0*(-1.0 + 4.0*CW2)*mH2_6*mHp2 +
2186 12.0*CW4*(1.0 - 4.0*CW2)*mH2_4*MZ4 -
2187 32*CW4*(10.0 - 49*CW2 + 36*CW4)*mHp2*
2188 MZ4*(-mHp2 + CW2*MZ2) +
2190 MZ2*((71.0 - 278.0*CW2 + 120*CW4)*mHp4 +
2191 CW2*(-71.0 + 278.0*CW2 - 120*CW4)*mHp2*MZ2 +
2192 48.0*CW4*(1.0 - 4.0*CW2)*MZ4))))/(768*
2193 CW4*(-1.0 + CW2)*(-1.0 + CW2)*(-1.0 + 4.0*CW2)*
mH2_2*
2194 mHp2*MZ4*(-mHp2 + CW2*MZ2)*(-MH2 + 4.0*CW2*MZ2)*pi*
2196 mMU*(2.0*mHp8 + MH6*(mHp2 + CW2*MZ2) -
2197 MH2*(3.0*mHp6 + 4.0*CW2*mHp4*MZ2))*
2198 log (MH2/mHp2))/(128.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*
2199 mHp4*MZ6*pi*pi) - (aem*aem*
mMU*
2200 mMU*(2.0*mHp8 + mH2_6*(mHp2 + CW2*MZ2) -
2201 mH2_2*(3.0*mHp6 + 4.0*CW2*mHp4*MZ2))*
2202 log (
mH2_2/mHp2))/(128.0*CW6*(-1.0 + CW2)*(-1.0 + CW2)*
2203 mHp4*MZ6*pi*pi) + (aem*aem*
mMU*
2204 mMU*(-2.0*MH6 + 3.0*CW2*mHp4*MZ2 - 4.0*CW4*mHp2*MZ4 +
2205 CW6*MZ6 + MH4*(mHp2 + 3.0*CW2*MZ2) +
2206 MH2*(4.0*mHp4 + 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4))*
2207 log (MH2/MZ2))/(128.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2208 MZ2*(-mHp2 + CW2*MZ2)*pi*pi) - (aem*aem*
mMU*
2209 mMU*(-2.0*mH2_6 + 3.0*CW2*mHp4*MZ2 - 4.0*CW4*mHp2*MZ4 +
2210 CW6*MZ6 + mH2_4*(mHp2 + 3.0*CW2*MZ2) +
2211 mH2_2*(4.0*mHp4 + 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4))*
2212 log (
mH2_2/MZ2))/(128.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2213 mHp4*MZ2*(-mHp2 + CW2*MZ2)*pi*pi) + (aem*aem*
mMU*
2214 mMU*(MH4*(mHp2 + 3.0*CW2*MZ2) +
2215 MH2*(-4.0*mHp4 + 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4) -
2216 mH2_2*(-4.0*mHp4 + 3.0*CW2*mHp2*MZ2 - 3.0*CW4*MZ4 +
2217 mH2_2*(mHp2 + 3.0*CW2*MZ2)))*
2218 log (mHp2/MZ2))/(128.0*CW2*(-1.0 + CW2)*(-1.0 + CW2)*
2219 mHp4*MZ2*(-mHp2 + CW2*MZ2)*pi*pi)) + (aem*aem*
mMU*
2220 mMU*(2.0*mHp8 + MH6*(mHp2 + CW2*MZ2) -
2221 MH2*(3.0*mHp6 + 4.0*CW2*mHp4*MZ2))*
2223 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ6*pi*pi) - (aem*aem*
2224 mMU*
mMU*(2.0*mHp8 + mH2_6*(mHp2 + CW2*MZ2) -
2225 mH2_2*(3.0*mHp6 + 4.0*CW2*mHp4*MZ2))*
2227 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ6*pi*pi)+ (aem*aem*
2228 mMU*
mMU*(2.0*MH6*(mHp2 + CW2*MZ2) -
2229 CW4*MZ4*(mHp4 - 6.0*CW2*mHp2*MZ2 + 2.0*CW4*MZ4) -
2230 2.0*MH4*(mHp4 + 4.0*CW2*mHp2*MZ2 + 3.0*CW4*MZ4) +
2231 MH2*(4.0*CW2*mHp4*MZ2 + 6.0*CW6*MZ6))*
2233 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*MZ6*pi*pi)+ (aem*aem*
2234 mMU*
mMU*(12.0*(1.0 - 3.0*CW2 + 2.0*CW4)*mH2_10 -
2235 48.0*(1.0 - 3.0*CW2 + 2.0*CW4)*mH2_8*
2236 MZ2 + (17.0 - 48.0*CW2 + 32.0*CW4)*mH2_6*MZ4 -
2237 3.0*(5.0 - 12.0*CW2 + 8.0*CW4)*
mH2_2*MZ8 -
2238 4.0*(5.0 - 12.0*CW2 + 8.0*CW4)*MZ10)*
2240 CW4*(-1.0 + CW2)*(-1.0 + CW2)*mH2_4*MZ8*pi*pi) + (aem*aem*
2241 mMU*
mMU*(-6.0*CW6*
mH2_2*mHp4*MZ6 - 8.0*CW8*mHp4*MZ8 -
2242 3.0*mH2_10*(mHp2 + CW2*MZ2) +
2243 3.0*mH2_8*(mHp4 + 4.0*CW2*mHp2*MZ2 + 3.0*CW4*MZ4) +
2244 mH2_6*(2.0*CW2*mHp4*MZ2 - 9.0*CW6*MZ6) +
2245 3.0*mH2_4*(-3.0*CW6*mHp2*MZ6 + CW8*MZ8))*
2247 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mH2_4*mHp4*MZ6*pi*pi) - (aem*
2248 aem*(-3.0 + 4.0*CW2)*
mMU*
2249 mMU*((1.0 - 5.0*CW2 + 10.0*CW4)*
2250 mH2_2 + (-7.0 + 61*CW2 - 162*CW4 + 96.0*CW6)*MZ2)*
2251 TF (CW, CW, 1.0))/(128.0*
2252 CW2*(1.0 - 4.0*CW2)*(1.0 - 4.0*CW2)*(-1.0 + CW2)*(-1.0 +
2253 CW2)*MZ2*(-
mH2_2 + MZ2)*pi*pi) + (aem*
2254 aem*(MH2 + 2.0*mHp2)*
mMU*
2255 mMU*(3.0*MH2*mHp8 - 8.0*CW6*mHp6*MZ4 +
2256 2.0*CW4*mHp6*MZ2*(MH2 + 4.0*(mHp2 + MZ2)) +
2257 CW2*(MH10 - 2.0*MH4*mHp4*MZ2 - 8.0*mHp8*MZ2 -
2258 MH8*(4.0*mHp2 + MZ2) +
2259 2.0*MH6*(mHp4 + 2.0*mHp2*MZ2) -
2260 MH2*(2.0*mHp8 + 3.0*mHp6*MZ2)))*
2261 TF (MH, mHp, mHp))/(64.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*MH4*
2262 mHp4*(MH2 - 4.0*mHp2)*MZ2*(MH2 - MZ2)*(-mHp2 + CW2*MZ2)*pi*
2264 mMU*(MH8*(-mHp4 + CW4*MZ4) +
2265 MH6*(mHp6 + CW2*mHp4*MZ2 - 2.0*CW4*mHp2*MZ4 -
2266 4.0*CW6*MZ6) + (mHp2 - CW2*MZ2)*(mHp2 -
2267 CW2*MZ2)*(mHp2 - CW2*MZ2)*(2.0*mHp6 +
2268 2.0*CW2*mHp4*MZ2 + 2.0*CW4*mHp2*MZ4 - CW6*MZ6) +
2269 MH4*(3.0*mHp8 + 3.0*CW2*mHp6*MZ2 + CW4*mHp4*MZ4 -
2270 CW6*mHp2*MZ6 + 6.0*CW8*MZ8) +
2271 MH2*(-5.0*mHp10 + CW4*mHp6*MZ4 + 8.0*CW8*mHp2*MZ8 -
2272 4.0*CW10*MZ10))*
TF (MH, mHp, CW*MZ))/(64.0*
2273 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2275 CW2*MZ2)*(MH4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2276 2.0*MH2*(mHp2 + CW2*MZ2))*pi*pi) + (aem*aem*
mMU*
2277 mMU*(-2.0*MH10 + 2.0*CW8*MZ8*(-mHp2 + CW2*MZ2) +
2278 2.0*MH8*(mHp2 + 9.0*CW2*MZ2) -
2279 4.0*MH6*(4.0*CW2*mHp2*MZ2 + 11.0*CW4*MZ4) +
2280 MH4*(37*CW4*mHp2*MZ4 + 3.0*CW6*MZ6) +
2281 2.0*MH2*(-9.0*CW6*mHp2*MZ6 + 25*CW8*MZ8))*
2282 TF (MH, CW*MZ, CW*MZ))/(128.0*
2283 CW6*(-1.0 + CW2)*(-1.0 + CW2)*MH2*
2284 MZ6*(MH2 - 4.0*CW2*MZ2)*(MH2 - 4.0*CW2*MZ2)*(-mHp2 +
2285 CW2*MZ2)*pi*pi) + (aem*aem*(
mH2_2 + 2.0*mHp2)*
mMU*
2286 mMU*(-3.0*
mH2_2*mHp8 + 8.0*CW6*mHp6*MZ4 -
2287 2.0*CW4*mHp6*MZ2*(
mH2_2 + 4.0*(mHp2 + MZ2)) +
2289 CW2*(-mH2_10 + 2.0*mH2_4*mHp4*MZ2 + 8.0*mHp8*MZ2 +
2290 mH2_8*(4.0*mHp2 + MZ2) -
2291 2.0*mH2_6*(mHp4 + 2.0*mHp2*MZ2) +
2292 mH2_2*(2.0*mHp8 + 3.0*mHp6*MZ2)))*
2293 TF (mH2, mHp, mHp))/(64.0*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
2294 mH2_4*mHp4*(
mH2_2 - 4.0*mHp2)*
2295 MZ2*(
mH2_2 - MZ2)*(-mHp2 + CW2*MZ2)*pi*pi) - (aem*aem*
mMU*
2296 mMU*(mH2_8*(-mHp4 + CW4*MZ4) +
2297 mH2_6*(mHp6 + CW2*mHp4*MZ2 - 2.0*CW4*mHp2*MZ4 -
2298 4.0*CW6*MZ6) + (mHp2 - CW2*MZ2)*(mHp2 -
2299 CW2*MZ2)*(mHp2 - CW2*MZ2)*(2.0*mHp6 +
2300 2.0*CW2*mHp4*MZ2 + 2.0*CW4*mHp2*MZ4 - CW6*MZ6) +
2301 mH2_4*(3.0*mHp8 + 3.0*CW2*mHp6*MZ2 + CW4*mHp4*MZ4 -
2302 CW6*mHp2*MZ6 + 6.0*CW8*MZ8) +
2303 mH2_2*(-5.0*mHp10 + CW4*mHp6*MZ4 + 8.0*CW8*mHp2*MZ8 -
2304 4.0*CW10*MZ10))*
TF (mH2, mHp, CW*MZ))/(64.0*
2305 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mHp4*
2307 CW2*MZ2)*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2308 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi) + (aem*aem*
mMU*
2309 mMU*(6.0*(1.0 - 3.0*CW2 + 2.0*CW4)*mH2_6 -
2310 12.0*(1.0 - 3.0*CW2 + 2.0*CW4)*mH2_4*
2311 MZ2 + (5.0 - 12.0*CW2 + 8.0*CW4)*
mH2_2*MZ4 +
2312 2.0*(5.0 - 12.0*CW2 + 8.0*CW4)*MZ6)*
2313 TF (mH2, MZ, MZ))/(192*CW4*(-1.0 + CW2)*(-1.0 + CW2)*
mH2_2*
2314 MZ8*pi*pi) + (aem*aem*
mMU*
2316 2.0*CW6*(-95 + 78*CW2)*
mH2_2*MZ8*(-mHp2 + CW2*MZ2) +
2317 144*CW8*(-1.0 + CW2)*MZ10*(-mHp2 + CW2*MZ2) +
2319 MZ4*((-2.0 + 5.0*CW2)*mHp2 + CW2*(-7.0 + CW2)*MZ2) -
2320 3.0*mH2_10*(mHp2 + (1.0 + 5.0*CW2)*MZ2) +
2322 MZ6*((17.0 + 10.0*CW2)*mHp2 - CW2*(29 + 10.0*CW2)*MZ2) +
2323 mH2_8*MZ2*((3.0 + 4.0*CW2)*mHp2 +
2324 CW2*(15 + 14*CW2)*MZ2))*
TF (mH2, CW*MZ, CW*MZ))/(192*
2325 CW6*(-1.0 + CW2)*(-1.0 + CW2)*mH2_4*
2326 MZ6*(
mH2_2 - MZ2)*(
mH2_2 - 4.0*CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2327 pi*pi) + (aem*aem*(3.0 - 10.0*CW2 + 8.0*CW4)*(MH2 -
mH2_2)*
2328 mHp2*
mMU*
mMU*(2.0*mHp2 + MZ2)*
TF(mHp, mHp, MZ))/(64.0*
2329 CW4*(-1.0 + CW2)*(-1.0 + CW2)*
2330 MZ4*(MH2 - MZ2)*(-
mH2_2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi)));
2336 mMU*((-3.0 + 2.0*CW2)*MH2 - 8.0*CW2*(-1.0 + CW2)*MZ2)*
2337 log (MH2/MZ2))/(128.0*CW2*(-1.0 + CW2)*MH2*(MH2 - MZ2)*pi*pi*
2338 pi) - (aem*
mMU*
mMU*
log(mHp2/MZ2))/(16.0*MH2*pi*pi*pi) + (aem*
2339 mHp2*
mMU*
mMU*((-3.0 + 2.0*CW2)*MH2 - 8.0*CW2*(-1.0 + CW2)*MZ2)*
2340 TF (MH, mHp, mHp))/(64.0*CW2*(-1.0 + CW2)*
2341 MH4*(MH2 - 4.0*mHp2)*(MH2 - MZ2)*pi*pi*
2342 pi) - (aem*(3.0 - 10.0*CW2 + 8.0*CW4)*mHp2*
mMU*
mMU*
2343 TF (mHp, mHp, MZ))/(64.0*CW2*(-1.0 + CW2)*
2344 MZ2*(-MH2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi*pi) +
sl*((aem*
mMU*
2346 8.0*CW2*(-1.0 + CW2)*mHp2*MZ2)*(-1.0 + tanb*tanb))/(128.0*
2347 CW2*(-1.0 + CW2)*
mH2_2*mHp2*MZ2*pi*pi*pi*tanb) - (aem*
2350 CW2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2351 CW2*MZ2)*pi*pi*pi*tanb) + (aem*
mH2_2*
mMU*
2352 mMU*(-7.0*mHp6 + 8.0*CW2*mHp4*MZ2 - 4.0*CW4*mHp2*MZ4 +
2353 2.0*
mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*(-1.0 +
2356 log(
mH2_2/mHp2))/(1024*(-1.0 + CW2)*
2357 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2358 pi*pi*pi*tanb) + (aem*
mH3_2*mHp2*
mMU*
2359 mMU*(-1.0 + tanb*tanb)*
2361 CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*pi*pi*
2363 mMU*(-7.0*mHp6 + 8.0*CW2*mHp4*MZ2 - 4.0*CW4*mHp2*MZ4 +
2364 2.0*
mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*(-1.0 +
2367 log (
mH2_2/MZ2))/(1024*(-1.0 + CW2)*
2368 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2369 pi*pi*pi*tanb) + (aem*
mMU*
2370 mMU*(2.0*CW2*mH2_4*MZ2 +
2371 8.0*CW2*(-1.0 + CW2)*mHp2*MZ2*(-mHp2 + CW2*MZ2) +
2372 mH2_2*(-2.0*mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
2373 log (mHp2/MZ2))/(256.0*CW2*(-1.0 + CW2)*
mH2_2*mHp2*
2374 MZ2*(-mHp2 + CW2*MZ2)*pi*pi*pi*tanb) + (aem*
mMU*
2375 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
2376 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
2377 log (
mH2_2/mHp2)*
log (mHp2/MZ2))/(256.0*CW4*(-1.0 + CW2)*
2378 mHp4*MZ4*pi*pi*pi*tanb) - (aem*
mH3_2*mHp2*
mMU*
2380 log (mHp2/MZ2))/(512.0*(-1.0 + CW2)*(-mHp2 +
2381 CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*pi*pi*
2384 2.0*mH2_4*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2385 mH2_2*(mHp6 - 8.0*CW2*mHp4*MZ2 +
2386 4.0*CW4*mHp2*MZ4))*(-1.0 + tanb*tanb)*
2388 log (mHp2/MZ2))/(1024*(-1.0 + CW2)*
2389 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*
2394 2.0*mH2_4*MZ2 + (3.0 - 2.0*CW2)*
mH2_2*mHp2*MZ2 +
2395 8.0*CW2*(-1.0 + CW2)*mHp2*MZ4)*(-1.0 +
2396 tanb*tanb))/(256.0*CW2*(-1.0 + CW2)*
mH2_2*mHp2*
2397 MZ2*(
mH2_2 - MZ2)*pi*pi*pi*tanb) - (3.0*aem*
mH2_2*mHp2*
2398 mMU*
mMU*(-1.0 + tanb*tanb)*
2399 log (mHp2/MZ2))/(512.0*(-1.0 + CW2)*(-mHp2 +
2400 CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*pi*pi*
2403 mMU*(2.0*
mH2_2 - 2.0*mHp2 + CW2*MZ2)*(-1.0 +
2404 tanb*tanb))/(256.0*CW2*(-1.0 + CW2)*
2405 MZ2*(-mHp2 + CW2*MZ2)*pi*pi*pi*tanb)
2407 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
2408 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
2409 log (
mH2_2/mHp2))/(256.0*CW4*(-1.0 + CW2)*mHp4*MZ4*pi*pi*
2411 mMU*(-1.0 + tanb*tanb)*
2412 log (
mH2_2/MZ2))/(256.0*(-1.0 + CW2)*
2413 mHp4*(-mHp2 + CW2*MZ2)*pi*pi*pi*tanb) - (aem*
2415 log (mHp2/MZ2))/(256.0*(-1.0 + CW2)*mHp4*(-mHp2 + CW2*MZ2)*
2416 pi*pi*pi*tanb) ) -(aem*
mMU*
2417 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
2418 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
2420 pi*pi*pi*tanb) + (aem*
mMU*
2421 mMU*(-3.0*
mH2_2*mHp8 + 8.0*CW6*mHp6*MZ4 -
2422 2.0*CW4*mHp6*MZ2*(
mH2_2 + 4.0*(mHp2 + MZ2)) +
2423 CW2*(-mH2_10 + 2.0*mH2_4*mHp4*MZ2 + 8.0*mHp8*MZ2 +
2424 mH2_8*(4.0*mHp2 + MZ2) -
2425 2.0*mH2_6*(mHp4 + 2.0*mHp2*MZ2) +
2426 mH2_2*(2.0*mHp8 + 3.0*mHp6*MZ2)))*(-1.0 + tanb*tanb)*
2427 TF (mH2, mHp, mHp))/(128.0*CW2*(-1.0 + CW2)*mH2_4*
2428 mHp4*(
mH2_2 - 4.0*mHp2)*(
mH2_2 - MZ2)*(-mHp2 + CW2*MZ2)*
2429 pi*pi*pi*tanb) - (aem*
mMU*
2431 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
2433 mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*(-1.0 + tanb*tanb)*
2434 TF (mH2, mHp, CW*MZ))/(128.0*CW4*(-1.0 + CW2)*
2436 CW2*MZ2)*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2437 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi*pi*
2438 tanb)+(aem*(3.0 - 10.0*CW2 + 8.0*CW4)*mHp2*
mMU*
2439 mMU*(-1.0 + tanb*tanb)*
TF (mHp, mHp, MZ))/(128.0*
2440 CW2*(-1.0 + CW2)*MZ2*(-
mH2_2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*
2441 pi*pi*tanb)) + eta*(-(aem*(MH2 -
mH2_2)*
mMU*
mMU*(-1.0 + tanb*tanb))/(16.0*MH2*
2443 mMU*((3.0 - 2.0*CW2)*MH2 + 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 +
2444 tanb*tanb)*
log (MH2/MZ2))/(256.0*CW2*(-1.0 + CW2)*
2445 MH2*(MH2 - MZ2)*pi*pi*pi*tanb) + (aem*
mMU*
2447 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 + tanb*tanb)*
2448 log (
mH2_2/MZ2))/(256.0*CW2*(-1.0 + CW2)*
2451 pi*pi*pi*tanb) + (aem*mHp2*
mMU*
2452 mMU*((3.0 - 2.0*CW2)*MH2 + 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 +
2453 tanb*tanb)*
TF (MH, mHp, mHp))/(128.0*CW2*(-1.0 + CW2)*
2454 MH4*(MH2 - 4.0*mHp2)*(MH2 - MZ2)*pi*pi*pi*tanb)+ (aem*mHp2*
2456 8.0*CW2*(-1.0 + CW2)*MZ2)*(-1.0 + tanb*tanb)*
2457 TF (mH2, mHp, mHp))/(128.0*CW2*(-1.0 + CW2)*
2458 mH2_4*(
mH2_2 - 4.0*mHp2)*(
mH2_2 - MZ2)*pi*pi*pi*
2459 tanb) + (aem*(3.0 - 10.0*CW2 + 8.0*CW4)*(MH2 -
mH2_2)*mHp2*
2460 mMU*
mMU*(-1.0 + tanb*tanb)*
TF (mHp, mHp, MZ))/(128.0*
2462 MZ2*(-MH2 + MZ2)*(-
mH2_2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi*pi*
2464 mMU*(MH2*
mH2_2 + 8.0*CW2*(-1.0 + CW2)*mHp2*MZ2))/(64.0*
2465 CW2*(-1.0 + CW2)*MH2*
mH2_2*mHp2*MZ2*pi*pi*pi) - (aem*
2467 mMU*(-3.0*mHp6 + 4.0*CW2*mHp4*MZ2 - 2.0*CW4*mHp2*MZ4 +
2468 MH2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*
2470 log (MH2/mHp2))/(256.0*(-1.0 + CW2)*
2471 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2472 CW2*MZ2)*pi*pi*pi) + (aem*
mH2_2*
mMU*
2473 mMU*(-3.0*mHp6 + 4.0*CW2*mHp4*MZ2 - 2.0*CW4*mHp2*MZ4 +
2474 mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*
2476 log (
mH2_2/mHp2))/(256.0*(-1.0 + CW2)*
2477 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2478 CW2*MZ2)*pi*pi*pi) + (aem*MH2*
mMU*
2479 mMU*(-3.0*mHp6 + 4.0*CW2*mHp4*MZ2 - 2.0*CW4*mHp2*MZ4 +
2480 MH2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*
2482 log (MH2/MZ2))/(256.0*(-1.0 + CW2)*
2483 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2484 CW2*MZ2)*pi*pi*pi) - (aem*
mH2_2*
mMU*
2485 mMU*(-3.0*mHp6 + 4.0*CW2*mHp4*MZ2 - 2.0*CW4*mHp2*MZ4 +
2486 mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2))*
2488 log (
mH2_2/MZ2))/(256.0*(-1.0 + CW2)*
2490 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2491 CW2*MZ2)*pi*pi*pi) - (aem*(MH2 -
mH2_2)*
mMU*
2493 4.0*(-1.0 + CW2)*mHp2*(mHp2 - CW2*MZ2))*
2494 log (mHp2/MZ2))/(64.0*(-1.0 + CW2)*MH2*
mH2_2*
2495 mHp2*(-mHp2 + CW2*MZ2)*pi*pi*pi) - (aem*
mMU*
2496 mMU*(mHp6 + MH4*(mHp2 + CW2*MZ2) -
2497 2.0*MH2*(mHp4 + CW2*mHp2*MZ2))*
log (MH2/mHp2)*
2498 log (mHp2/MZ2))/(128.0*CW4*(-1.0 + CW2)*mHp4*MZ4*pi*pi*
2500 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
2502 log (mHp2/MZ2))/(128.0*CW4*(-1.0 + CW2)*mHp4*MZ4*pi*pi*
2504 mMU*(-(MH4*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2)) +
2505 MH2*(mHp6 - 4.0*CW2*mHp4*MZ2 + 2.0*CW4*mHp2*MZ4) +
2506 mH2_2*(
mH2_2*(mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2507 mHp2*(mHp4 - 4.0*CW2*mHp2*MZ2 + 2.0*CW4*MZ4)))*
2509 log(mHp2/MZ2))/(256.0*(-1.0 + CW2)*
2510 mHp4*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 +
2511 CW2*MZ2)*pi*pi*pi) +
2515 2.0*MH4*MZ2 + (-3.0 + 2.0*CW2)*MH2*mHp2*MZ2 -
2516 8.0*CW2*(-1.0 + CW2)*mHp2*MZ4))/(128.0*
2517 CW2*(-1.0 + CW2)*MH2*mHp2*MZ2*(MH2 - MZ2)*pi*pi*
2518 pi) + (aem*MH2*mHp2*
mMU*
mMU*
2519 log (mHp2/MZ2))/(128.0*(-1.0 + CW2)*(-mHp2 +
2520 CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*pi*
2524 2.0*mH2_4*MZ2 + (3.0 - 2.0*CW2)*
mH2_2*mHp2*MZ2 +
2525 8.0*CW2*(-1.0 + CW2)*mHp2*MZ4))/(128.0*
2526 CW2*(-1.0 + CW2)*
mH2_2*mHp2*MZ2*(
mH2_2 - MZ2)*pi*pi*
2528 log (mHp2/MZ2))/(128.0*(-1.0 + CW2)*(-mHp2 +
2529 CW2*MZ2)*(-mHp2 + CW2*MZ2)*(-mHp2 + CW2*MZ2)*pi*pi*
2531 CW2*(-1.0 + CW2)*MZ2*(-mHp2 + CW2*MZ2)*pi*pi*
2533 mMU*(mHp6 + MH4*(mHp2 + CW2*MZ2) -
2534 2.0*MH2*(mHp4 + CW2*mHp2*MZ2))*
2535 log (MH2/mHp2))/(128.0*CW4*(-1.0 + CW2)*mHp4*MZ4*pi*
2537 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
2538 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*
2539 log (
mH2_2/mHp2))/(128.0*CW4*(-1.0 + CW2)*mHp4*MZ4*
2540 pi*pi*pi) - (aem*MH2*(MH2 - 2.0*mHp2)*
mMU*
mMU*
2541 log (MH2/MZ2))/(128.0*(-1.0 + CW2)*
2542 mHp4*(-mHp2 + CW2*MZ2)*pi*pi*pi) + (aem*
2544 log (
mH2_2/MZ2))/(128.0*(-1.0 + CW2)*
2545 mHp4*(-mHp2 + CW2*MZ2)*pi*pi*
2546 pi) + (aem*(MH2 -
mH2_2)*(MH2 +
mH2_2 - 2.0*mHp2)*
2547 mMU*
mMU*
log (mHp2/MZ2))/(128.0*(-1.0 + CW2)*
2548 mHp4*(-mHp2 + CW2*MZ2)*pi*pi*pi)) + (aem*
mMU*
2549 mMU*(mHp6 + MH4*(mHp2 + CW2*MZ2) -
2550 2.0*MH2*(mHp4 + CW2*mHp2*MZ2))*
2551 PolyLog.
Li2 (1.0 - MH2/mHp2))/(64.0*CW4*(-1.0 + CW2)*
2552 mHp4*MZ4*pi*pi*pi) - (aem*
mMU*
2553 mMU*(mHp6 + mH2_4*(mHp2 + CW2*MZ2) -
2554 2.0*
mH2_2*(mHp4 + CW2*mHp2*MZ2))*
2556 mHp4*MZ4*pi*pi*pi) + (aem*
mMU*
2557 mMU*(3.0*MH2*mHp8 - 8.0*CW6*mHp6*MZ4 +
2558 2.0*CW4*mHp6*MZ2*(MH2 + 4.0*(mHp2 + MZ2)) +
2559 CW2*(MH10 - 2.0*MH4*mHp4*MZ2 - 8.0*mHp8*MZ2 -
2560 MH8*(4.0*mHp2 + MZ2) +
2561 2.0*MH6*(mHp4 + 2.0*mHp2*MZ2) -
2562 MH2*(2.0*mHp8 + 3.0*mHp6*MZ2)))*
2563 TF (MH, mHp, mHp))/(64.0*CW2*(-1.0 + CW2)*MH4*
2564 mHp4*(MH2 - 4.0*mHp2)*(MH2 - MZ2)*(-mHp2 + CW2*MZ2)*pi*
2567 MH4*(3.0*mHp2 + CW2*MZ2) + (mHp3 - CW2*mHp*MZ2)*(mHp3 -
2568 CW2*mHp*MZ2) + MH2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
2569 TF (MH, mHp, CW*MZ))/(64.0*CW4*(-1.0 + CW2)*
2571 CW2*MZ2)*(MH4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2572 2.0*MH2*(mHp2 + CW2*MZ2))*pi*pi*pi) + (aem*
mMU*
2573 mMU*(-3.0*
mH2_2*mHp8 + 8.0*CW6*mHp6*MZ4 -
2574 2.0*CW4*mHp6*MZ2*(
mH2_2 + 4.0*(mHp2 + MZ2)) +
2575 CW2*(-mH2_10 + 2.0*mH2_4*mHp4*MZ2 + 8.0*mHp8*MZ2 +
2576 mH2_8*(4.0*mHp2 + MZ2) -
2577 2.0*mH2_6*(mHp4 + 2.0*mHp2*MZ2) +
2578 mH2_2*(2.0*mHp8 + 3.0*mHp6*MZ2)))*
2579 TF (mH2, mHp, mHp))/(64.0*CW2*(-1.0 + CW2)*mH2_4*
2580 mHp4*(
mH2_2 - 4.0*mHp2)*(
mH2_2 - MZ2)*(-mHp2 +
2581 CW2*MZ2)*pi*pi*pi) - (aem*
mMU*
2583 mH2_4*(3.0*mHp2 + CW2*MZ2) + (mHp3 -
2584 CW2*mHp*MZ2)*(mHp3 - CW2*mHp*MZ2) +
2585 mH2_2*(-3.0*mHp4 + CW2*mHp2*MZ2))*
2586 TF (mH2, mHp, CW*MZ))/(64.0*CW4*(-1.0 + CW2)*
2588 CW2*MZ2)*(mH2_4 + (mHp2 - CW2*MZ2)*(mHp2 - CW2*MZ2) -
2589 2.0*
mH2_2*(mHp2 + CW2*MZ2))*pi*pi*
2590 pi) + (aem*(3.0 - 10.0*CW2 + 8.0*CW4)*(MH2 -
mH2_2)*
2591 mHp2*
mMU*
mMU*
TF (mHp, mHp, MZ))/(64.0*CW2*(-1.0 + CW2)*
2592 MZ2*(MH2 - MZ2)*(-
mH2_2 + MZ2)*(-4.0*mHp2 + MZ2)*pi*pi*
2596 aYuk = aYuk2+aYuk3+aYuk4+aYuk5;
2620 double gminus2muNLO = gminus2muNLOFvalue + gminus2muNLOBvalue;
2641 vmcgminus2mu = StandardModelMatching::CMgminus2mu();
2662 std::stringstream out;
2664 throw std::runtime_error(
"GeneralTHDMMatching::CMgminus2mu(): order " + out.str() +
" not implemented.\nOnly leading order (LO) or next-to-leading order (NLO) are allowed.");
2668 return(vmcgminus2mu);
2675 double xt = x_t(Mut);
2681 double xHW=mHp2/(MW*MW);
2686 double SWH=xtH*((2.0*xHW-8.0)*
log(xtH)/((1.0-xHW)*(1.0-xtH)*(1.0-xtH))+6.0*xHW*
log(xt)/((1.0-xHW)*(1.0-xt)*(1.0-xt))-(8.0-2.0*xt)/((1.0-xt)*(1.0-xtH)))*
su*
su;
2687 double SHH=xtH*((1.0+xtH)/((1.0-xtH)*(1.0-xtH))+2.0*xtH*
log(xtH)/((1.0-xtH)*(1.0-xtH)*(1.0-xtH)))*
su*
su*
su*
su;
2688 double C1bsSRR = 4.0*mb*mb*xt*xtH*
sd*
su/(mHp2*(xtH-1.0)*(xtH-1.0)*(xtH-1.0))
2689 * (
sd*
su*(2.0*(xtH-1.0)-(xtH+1.0)*
log(xtH))
2690 + (2.0*xt*xt*(xtH-1.0)*(xtH-1.0)*(xtH-1.0)*
log(xt)/(xt-1.0)
2691 +2.0*xt*(xtH-1.0)*((xt-xtH)*(xtH-1.0)+(xtH-xt*xtH)*
log(xtH)))/((xt-1.0)*(xt-xtH)*xtH));
2703 std::stringstream out;
2705 throw std::runtime_error(
"THDMMatching::CMdbs2(): order " + out.str() +
"not implemented");
2719 std::stringstream out;
2721 throw std::runtime_error(
"THDMMatching::CMdbs2(): order " + out.str() +
"not implemented");
2769 double C10 =
su.
abs2()*xt*xt/8*(1/(xHp-xt) + xHp/((xHp-xt)*(xHp-xt))*(
log(xt)-
log(xHp)));
2776 - xHp/(xHp-1)*
log(xHp))
2801 + ((xt*xHp)/(xHp - xt))*(
log(xHp) -
log(xt))) +
2802 su.
abs2()*(1/(6*(xHp - xt)))*((xHp*xHp - 8*xHp*xt - 17*xt*xt)/ 6
2803 + ((xt*xt*(3*xHp + xt))/(xHp - xt))*(
log(xHp) -
log(xt))))
2804 + sW2*(xt/(6*(xHp - xt)*(xHp - xt)))*(
sd*
su.
conjugate()*((5*xt - 3*xHp)/2
2805 + ((xHp*(2*xHp - 3*xt))/(xHp - xt))*(
log(xHp) -
log(xt))) -
2806 su.
abs2()*(1/(6*(xHp - xt)))*(((4*xHp*xHp*xHp - 12*xHp*xHp*xHp*xt + 9*xHp*xt*xt + 3*xt*xt*xt)/(xHp - xt))* (
log(xHp) -
log(xt))
2807 - (17*xHp*xHp - 64*xHp*xt + 71.0*xt*xt)/6));
2812 (xHp*(
log(xHp) -
log(xt)) + xt - xHp);
2825 double f1 = (1.0/(2.0*(xHp - xt)))*(-xHp + xt + xHp*
log(xHp) -
2834 double f2 = (1.0/(2.0*(xHp - xt)))*(xt - ((xHp*xt)/(xHp - xt))*(
log(xHp) -
log(xt)));
2842 double f3= (1.0/(2.0*(xHp - xt)))*
2843 (xHp - (xHp*xHp*
log(xHp))/(xHp - xt) + (xt*(2*xHp - xt)*
log(xt))/
2853 double f4 = (1.0/(4*(xHp - xt)*(xHp - xt)))*((xt*(3.0*xHp - xt))/2.0 -
2854 ((xHp*xHp*xt)/(xHp - xt))*(
log(xHp) -
log(xt)));
2862 double f5 = (1.0/(4*(xHp - xt)*(xHp - xt)))*((xt*(xHp - 3.0*xt))/2.0 -
2863 ((xHp*xt*(xHp - 2.0*xt))/(xHp - xt))*(
log(xHp) -
log(xt)));
2871 double f6 = (1.0/(2.0*(xHp - xt)))*
2872 ((xt*(xt*xt - 3.0*xHp*xt + 9*xHp - 5*xt - 2.0))/(4*(xt - 1.0)*(xt - 1.0)) +
2873 ((xHp*(xHp*xt - 3.0*xHp + 2.0*xt))/(2.0*(xHp - 1.0)*(xHp - xt)))*
2874 log(xHp) + ((xHp*xHp*(-2.0*xt*xt*xt + 6*xt*xt - 9*xt + 2.0) +
2875 3.0*xHp*xt*xt*(xt*xt - 2.0*xt + 3.0) - xt*xt*(2.0*xt*xt*xt - 3.0*xt*xt +
2876 3.0*xt + 1.0))/(2.0*(xt - 1.0)*(xt - 1.0)*(xt - 1.0)*(xHp - xt)))*
log(xt));
2885 double f7= (1.0/(2.0*(xHp - xt)))*
2886 (((xt*xt + xt - 8)*(xHp - xt))/(4*(xt - 1.0)*(xt - 1.0)) -
2887 ((xHp*(xHp + 2.0))/(2.0*(xHp - 1.0)))*
log(xHp) +
2888 ((xHp*(xt*xt*xt - 3.0*xt*xt + 3.0*xt + 2.0) + 3.0*(xt - 2.0)*xt*xt)/
2889 (2.0*(xt - 1.0)*(xt - 1.0)*(xt - 1.0)))*
log(xt));
2896 double f8 = (1.0/(4*(xHp - xt)))*((xt*
log(xt))/(xt - 1.0) -
2897 (xHp*
log(xHp))/(xHp - 1.0));
2906 double f9 = (1.0/(8*(xHp - xt)))*(xHp/(xHp - 1.0) +
2907 (xt*xt*
log(xt))/((xt - 1.0)*(xHp - xt)) -
2908 ((xHp*(xHp*xt + xHp - 2.0*xt))/((xHp - 1.0)*(xHp - 1.0)*(xHp - xt)))*
2915 double f10 = (1.0/(8*(xHp - xt)))*
2916 ((xHp - xt)/((xHp - 1.0)*(xt - 1.0)) + ((xt*(xt - 2.0))/(xt - 1.0)*(xt - 1.0))*
2917 log(xt) - ((xHp*(xHp - 2.0))/(xHp - 1.0)*(xHp - 1.0))*
log(xHp));
2927 +
su.
abs2()*(xt*xt/(2.0*(xHp-xt)*(xHp-xt))*(
log(xHp)-
log(xt)) + (xHp-3.0*xt)/(4.0*(xHp-xt))));
2935 +
su.
abs2()*(xt/(2.0*(xHp - xt)*(xHp - xt)))*((xHp + xt)/2.0 - ((xHp*xt)/(xHp - xt))*(
log(xHp) -
log(xt)));
2968 gslpp::complex GeneralTHDMMatching::CphiU(
double xHp,
double xt,
double vev,
double xphi,
double mu,
double Ri1,
double Ri2,
double Ri3,
double mHi_2,
double lambda3,
double Relambda7,
double Imlambda7,
gslpp::complex su,
gslpp::complex sd){
2971 (Ri2 + i*Ri3)*mu + (vev*vev/mHi_2)*
lambdaHHphi(
lambda3, Relambda7, Imlambda7, Ri1, Ri2, Ri3)*
g0(xHp, xt,
su,
sd) + Ri1*((1/(2.0*xphi))*
g1a(xHp, xt,
su,
sd)) + Ri2*(1/(2.0*xphi)*
g2a(xHp, xt,
su,
sd)) + i*Ri3*(1/(2.0*xphi))*
g3a(xHp, xt,
su,
sd));
3003 double xt = (Mt_muw*Mt_muw)/(MW*MW);
3004 double xHp = (mHp2)/(MW*MW);
3008 double mH1_2=mHl*mHl;
3046 double xphi1 =
mH1_2/(MW*MW);
3047 double xphi2 =
mH2_2/(MW*MW);
3048 double xphi3 =
mH3_2/(MW*MW);
3062 gslpp::complex CSphi1U = yl1.
real()*
CphiU(xHp, xt, vev, xphi1, mu, R11, R12, R13,
mH1_2,
lambda3, Relambda7, Imlambda7,
su,
sd);
3063 gslpp::complex CSphi2U = yl2.
real()*
CphiU(xHp, xt, vev, xphi2, mu, R21, R22, R23,
mH2_2,
lambda3, Relambda7, Imlambda7,
su,
sd);
3064 gslpp::complex CSphi3U = yl3.
real()*
CphiU(xHp, xt, vev, xphi3, mu, R31, R32, R33,
mH3_2,
lambda3, Relambda7, Imlambda7,
su,
sd);
3066 gslpp::complex CPphi1U = i*yl1.
imag()*
CphiU(xHp, xt, vev, xphi1, mu, R11, R12, R13,
mH1_2,
lambda3, Relambda7, Imlambda7,
su,
sd);
3067 gslpp::complex CPphi2U = i*yl2.
imag()*
CphiU(xHp, xt, vev, xphi2, mu, R21, R22, R23,
mH2_2,
lambda3, Relambda7, Imlambda7,
su,
sd);
3068 gslpp::complex CPphi3U = i*yl3.
imag()*
CphiU(xHp, xt, vev, xphi3, mu, R31, R32, R33,
mH3_2,
lambda3, Relambda7, Imlambda7,
su,
sd);
3073 gslpp::complex CPphiU = CPboxU + CPZU + CPphi1U + CPphi2U + CPphi3U;
3082 std::stringstream out;
3084 throw std::runtime_error(
"GeneralTHDMMatching::CMBMll(): scheme " + out.str() +
"not implemented");
3097 std::stringstream out;
3099 throw std::runtime_error(
"GeneralTHDMMatching::CMBMll(): order " + out.str() +
"not implemeted");
3101 vmcBMll.push_back(
mcBMll);
3115 vmcbtaunu = StandardModelMatching::CMbtaunu(meson_i);
3124 else throw std::runtime_error(
"GeneralTHDMMatching::CMbtaunu(): meson not implemeted");
3127 std::stringstream out;
3129 throw std::runtime_error(
"THDMMatching::CMbtaunu(): order " + out.str() +
"not implemented");
3143 throw std::runtime_error(
"bsgamma is only available in the ATHDM at the moment.");
3162 for (
int j=0; j<8; j++){
3166 for (
int j=0; j<8; j++){
3170 for (
int j=0; j<8; j++){
3175 std::stringstream out;
3177 throw std::runtime_error(
"THDMMatching::CMbsg(): order " + out.str() +
"not implemented");
3180 vmcbsg.push_back(
mcbsg);
3203 std::stringstream out;
3205 throw std::runtime_error(
"order" + out.str() +
"not implemeted");
3211 double x = mt*mt/mhp/mhp;
3219 double xm3 = xm2*xm;
3220 double xm4 = xm3*xm;
3221 double xm6 = xm4*xm2;
3222 double xm8 = xm4*xm4;
3223 double xo = 1. - 1./x;
3225 double xo4 = xo*xo2*xo;
3226 double xo6 = xo4*xo2;
3227 double xo8 = xo4*xo4;
3230 double Lx3 = Lx2*Lx;
3231 double Li2 = gsl_sf_dilog(1.-1./x);
3233 double abssu =
su.
abs();
3234 double abssu2 = abssu*abssu;
3237 double lstmu = 2. *
log(mu/mt);
3239 double n70ct = - ( (7. - 5.*x - 8.*x2)/(36.*xm3) + (2.*x - 3.*x2)*Lx/(6.*xm4) ) * x/2.;
3240 double n70fr = ( (3.*x - 5.*x2)/(6.*xm2) + (2.*x - 3.*x2)*Lx/(3.*xm3) ) / 2.;
3242 double n80ct = - ( (2. + 5.*x - x2)/(12.*xm3) + (x*Lx)/(2.*xm4) ) * x/2.;
3243 double n80fr = ( (3.*x - x2)/(2.*xm2) + (x*Lx)/xm3 ) / 2.;
3245 double n41ct = (-16.*x + 29.*x2 - 7.*x3)/(36.*xm3) + (-2.*x + 3.*x2)*Lx/(6.*xm4);
3247 double n71ct = (797.*x - 5436.*x2 + 7569.*x3 - 1202.*x4)/(486.*xm4) +
3248 (36.*x2 - 74.*x3 + 16.*x4)*
Li2/(9.*xm4) +
3249 ((7.*x - 463.*x2 + 807.*x3 - 63.*x4)*Lx)/(81.*xm3*xm2);
3250 double cd71ct = (-31.*x - 18.*x2 + 135.*x3 - 14.*x4)/(27.*xm4) + (-28.*x2 + 46.*x3 + 6.*x4)*Lx/(9.*xm3*xm2);
3251 double n71fr = (28.*x - 52.*x2 + 8.*x3)/(3.*xm3) + (-48.*x + 112.*x2 - 32.*x3)*
Li2/(9.*xm3) +
3252 (66.*x - 128.*x2 + 14.*x3)*Lx/(9.*xm4);
3253 double cd71fr = (42.*x - 94.*x2 + 16.*x3)/(9.*xm3) + (32.*x - 56.*x2 - 12.*x3)*Lx/(9.*xm4);
3255 double n81ct = (1130.*x - 18153.*x2 + 7650.*x3 - 4451.*x4)/(1296.*xm4) +
3256 (30.*x2 - 17.*x3 + 13.*x4)*
Li2/(6.*xm4) +
3257 (-2.*x - 2155.*x2 + 321.*x3 - 468.*x4)*Lx/(216.*xm3*xm2);
3258 double cd81ct = (-38.*x - 261.*x2 + 18.*x3 - 7.*x4)/(36.*xm4) + (-31.*x2 - 17.*x3)*Lx/(6.*xm3*xm2);
3259 double n81fr = (143.*x - 44.*x2 + 29.*x3)/(8.*xm3) + (-36.*x + 25.*x2 - 17.*x3)*
Li2/(6.*xm3) +
3260 (165.*x - 7.*x2 + 34.*x3)*Lx/(12.*xm4);
3261 double cd81fr = (81.*x - 16.*x2 + 7.*x3)/(6.*xm3) + (19.*x + 17.*x2)*Lx/(3.*xm4);
3263 double n32ct = (10.*x4 + 30.*x2 - 20.*x)/(27.*xm4) *
Li2 +
3264 (30.*x3 - 66.*x2 - 56.*x)/(81.*xm4) * Lx + (6.*x3 - 187.*x2 + 213.*x)/(-81.*xm3);
3265 double cd32ct = (-30.*x2 + 20.*x)/(27.*xm4)*Lx + (-35.*x3 + 145.*x2 - 80.*x)/(-81.*xm3);
3267 double n42ct = (515.*x4 - 906.*x3 + 99.*x2 + 182.*x)/(54.*xm4) *
Li2 +
3268 (1030.*x4 - 2763.*x3 - 15.*x2 + 980.*x)/(-108.*xm3*xm2)*Lx +
3269 (-29467.*x4 + 68142.*x3 - 6717.*x2 - 18134.*x)/(1944.*xm4);
3270 double cd42ct = (-375.*x3 - 95.*x2 + 182.*x)/(-54.*xm3*xm2)*Lx +
3271 (133.*x4 - 108.*x3 + 4023.*x2 - 2320.*x)/(324.*xm4);
3273 double cd72ct = -(x * (67930.*x4 - 470095.*x3 + 1358478.*x2 - 700243.*x + 54970.))/(-2187.*xm3*xm2) +
3274 (x * (10422.*x4 - 84390.*x3 + 322801.*x2 - 146588.*x + 1435.))/(729.*xm6)*Lx +
3275 (2.*x2 * (260.*x3 - 1515.*x2 + 3757.*x - 1446.))/(-27. * xm3*xm2) *
Li2;
3276 double ce72ct = (x * (-518.*x4 + 3665.*x3 - 17397.*x2 + 3767.*x + 1843.))/(-162.*xm3*xm2) +
3277 (x2 * (-63.*x3 + 532.*x2 + 2089.*x - 1118.))/(27.*xm6)*Lx;
3278 double cd72fr = -(x * (3790.*x3 - 22511.*x2 + 53614.*x - 21069.))/(81.*xm4) -
3279 (2.*x * (-1266.*x3 + 7642.*x2 - 21467.*x + 8179.))/(-81.*xm3*xm2)*Lx +
3280 (8.*x * (139.*x3 - 612.*x2 + 1103.*x - 342.))/(27.*xm4) *
Li2;
3281 double ce72fr = -(x * (284.*x3 - 1435.*x2 + 4304.*x - 1425.))/(27.*xm4) -
3282 (2.*x * (63.*x3 - 397.*x2 - 970.*x + 440.))/(-27.*xm3*xm2)*Lx;
3284 double cd82ct = -(x * (522347.*x4 - 2423255.*x3 + 2706021.*x2 - 5930609.*x + 148856))/(-11664.*xm3*xm2) +
3285 (x * (51948.*x4 - 233781.*x3 + 48634.*x2 - 698693.*x + 2452.))/(1944.*xm6)*Lx +
3286 (x2 * (481.*x3 - 1950.*x2 + 1523.*x - 2550.))/(-18.*xm3*xm2) *
Li2;
3287 double ce82ct = (x * (-259.*x4 + 1117.*x3 + 2925.*x2 + 28411.*x + 2366.))/(-216.*xm3*xm2) -
3288 (x2 * (139.*x2 + 2938.*x + 2683.))/(36.*xm6)*Lx;
3289 double cd82fr = -(x * (1463.*x3 - 5794.*x2 + 5543.*x - 15036.))/(27.*xm4) -
3290 (x * (-1887.*x3 + 7115.*x2 + 2519.*x + 19901.))/(-54.*xm3*xm2)*Lx -
3291 (x * (-629.*x3 + 2178.*x2 - 1729.*x + 2196.))/(18.*xm4) *
Li2;
3292 double ce82fr = -(x * (259.*x3 - 947.*x2 - 251.*x - 5973.))/(36.*xm4)-
3293 (x * (139.*x2 + 2134.*x + 1183.))/(-18.*xm3*xm2)*Lx;
3297 n72ct = -274.2/x5 - 72.13*Lx/x5 + 24.41/x4 - 168.3*Lx/x4 + 79.15/x3 -
3298 103.8*Lx/x3 + 47.09/x2 - 38.12*Lx/x2 + 15.35/x - 8.753*Lx/x + 3.970;
3300 n72ct = 1.283 + 0.7158 * xo + 0.4119 * xo2 + 0.2629 * xo*xo2 + 0.1825 * xo4 +
3301 0.1347 * xo*xo4 + 0.1040 * xo6 + 0.0831 * xo*xo6 + 0.06804 * xo8 +
3302 0.05688 * xo*xo8 + 0.04833 * xo2*xo8 + 0.04163 * xo*xo2*xo8 + 0.03625 * xo4*xo8 +
3303 0.03188 * xo*xo4*xo8 + 0.02827 * xo6*xo8 + 0.02525 * xo*xo6*xo8 + 0.02269 * xo8*xo8;
3304 else if (mhp < 520.)
3305 n72ct = 1.283 - 0.7158 * xm - 0.3039 * xm2 - 0.1549 * xm3 - 0.08625 * xm4 -
3306 0.05020 * xm3*xm2 - 0.02970 * xm6 - 0.01740 * xm3*xm4 - 0.009752 * xm8 -
3307 0.004877 * xm3*xm6 - 0.001721 * xm2*xm8 + 0.0003378 * xm3*xm8 + 0.001679 * xm4*xm8 +
3308 0.002542 * xm*xm4*xm8 + 0.003083 * xm6*xm8 + 0.003404 * xm*xm6*xm8 + 0.003574 * xm8*xm8;
3309 else n72ct = -823.0*x5 + 42.30*x5*Lx3 - 412.4*x5*Lx2 - 3362*x5*Lx -
3310 1492*x4 - 23.26*x4*Lx3 - 541.4*x4*Lx2 - 2540*x4*Lx -
3311 1158*x3 - 34.50*x3*Lx3 - 348.2*x3*Lx2 - 1292*x3*Lx -
3312 480.9*x2 - 20.73*x2*Lx3 - 112.4*x2*Lx2 - 396.1*x2*Lx -
3313 8.278*x + 0.9225*x*Lx2 + 4.317*x*Lx;
3317 n72fr = -( 194.3/x5 + 101.1*Lx/x5 - 24.97/x4 + 168.4*Lx/x4 - 78.90/x3 +
3318 106.2*Lx/x3 - 49.32/x2 + 38.43*Lx/x2 - 12.91/x + 9.757*Lx/x + 8.088 );
3320 n72fr = -( 12.82 - 1.663 * xo - 0.8852 * xo2 - 0.4827 * xo*xo2 - 0.2976 * xo4 -
3321 0.2021 * xo*xo4 - 0.1470 * xo6 - 0.1125 * xo*xo6 - 0.08931 * xo8 -
3322 0.07291 * xo*xo8 - 0.06083 * xo2*xo8 - 0.05164 * xo*xo2*xo8 - 0.04446 * xo4*xo8 -
3323 0.03873 * xo*xo4*xo8 - 0.03407 * xo6*xo8 - 0.03023 * xo*xo6*xo8 - 0.02702 * xo8*xo8 );
3324 else if (mhp < 400.)
3325 n72fr = -( 12.82 + 1.663 * xm + 0.7780 * xm2 + 0.3755 * xm3 + 0.1581 * xm4 +
3326 0.03021 * xm3*xm2 - 0.04868 * xm6 - 0.09864 * xm3*xm4 - 0.1306 * xm8 -
3327 0.1510 * xm3*xm6 - 0.1637 * xm2*xm8 - 0.1712 * xm3*xm8 - 0.1751 * xm4*xm8 -
3328 0.1766 * xm*xm4*xm8 - 0.1763 * xm6*xm8 - 0.1748 * xm*xm6*xm8 - 0.1724 * xm8*xm8 );
3329 else n72fr = -( 2828 * x5 - 66.63 * x5*Lx3 + 469.4 * x5*Lx2 + 1986 * x5*Lx +
3330 1480 * x4 + 36.08 * x4*Lx3 + 323.2 * x4*Lx2 + 169.9 * x4*Lx +
3331 166.7 * x3 + 19.73 * x3*Lx3 - 46.61 * x3*Lx2 - 826.2 * x3*Lx -
3332 524.1 * x2 - 8.889 * x2*Lx3 - 195.7 * x2*Lx2 - 870.3 * x2*Lx -
3333 572.2 * x - 20.94 * x*Lx3 - 123.5 * x*Lx2 - 453.5 * x*Lx );
3337 n82ct = 826.2/x5 - 300.7*Lx/x5 + 96.35/x4 + 91.89*Lx/x4 - 66.39/x3 +
3338 78.58*Lx/x3 - 39.76/x2 + 20.02*Lx/x2 - 5.214/x + 2.278;
3340 n82ct = 1.188 + 0.4078 * xo + 0.2002 * xo2 + 0.1190 * xo*xo2 + 0.07861 * xo4 +
3341 0.05531 * xo*xo4 + 0.04061 * xo6 + 0.03075 * xo*xo6 + 0.02386 * xo8 +
3342 0.01888 * xo*xo8 + 0.01520 * xo2*xo8 + 0.01241 * xo*xo2*xo8 + 0.01026 * xo4*xo8 +
3343 0.008575 * xo*xo4*xo8 + 0.007238 * xo6*xo8 + 0.006164 * xo*xo6*xo8 + 0.005290 * xo8*xo8;
3344 else if (mhp < 600.)
3345 n82ct = 1.188 - 0.4078 * xm - 0.2076 * xm2 - 0.1265 * xm3 - 0.08570 * xm4 -
3346 0.06204 * xm3*xm2 - 0.04689 * xm6 - 0.03652 * xm3*xm4 - 0.02907 * xm8 -
3347 0.02354 * xm3*xm6 - 0.01933 * xm2*xm8 - 0.01605 * xm3*xm8 - 0.01345 * xm4*xm8 -
3348 0.01137 * xm*xm4*xm8 - 0.009678 * xm6*xm8 - 0.008293 * xm*xm6*xm8 - 0.007148 * xm8*xm8;
3349 else n82ct = -19606 * x5 - 226.7 * x5*Lx3 - 5251 * x5*Lx2 - 26090 * x5*Lx -
3350 9016 * x4 - 143.4 * x4*Lx3 - 2244 * x4*Lx2 - 10102 * x4*Lx -
3351 3357 * x3 - 66.32 * x3*Lx3 - 779.6 * x3*Lx2 - 3077 * x3*Lx -
3352 805.5 * x2 - 22.98 * x2*Lx3 - 169.1 * x2*Lx2 - 602.7 * x2*Lx +
3353 0.7437 * x + 0.6908 * x*Lx2 + 3.238 * x*Lx;
3357 n82fr = -( -1003/x5 + 476.9*Lx/x5 - 205.7/x4 - 71.62*Lx/x4 + 62.26/x3 -
3358 110.7*Lx/x3 + 63.74/x2 - 35.42*Lx/x2 + 10.89/x - 3.174 );
3360 n82fr = -( -0.6110 - 1.095 * xo - 0.4463 * xo2 - 0.2568 * xo*xo2 - 0.1698 * xo4 -
3361 0.1197 * xo*xo4 - 0.08761 * xo6 - 0.06595 * xo*xo6 - 0.05079 * xo8 -
3362 0.03987 * xo*xo8 - 0.03182 * xo2*xo8 - 0.02577 * xo*xo2*xo8 - 0.02114 * xo4*xo8 -
3363 0.01754 * xo*xo4*xo8 - 0.01471 * xo6*xo8 - 0.01244 * xo*xo6*xo8 - 0.01062 * xo8*xo8 );
3364 else if (mhp < 520.)
3365 n82fr = -( -0.6110 + 1.095 * xm + 0.6492 * xm2 + 0.4596 * xm3 + 0.3569 * xm4 +
3366 0.2910 * xm3*xm2 + 0.2438 * xm6 + 0.2075 * xm3*xm4 + 0.1785 * xm8 +
3367 0.1546 * xm3*xm6 + 0.1347 * xm2*xm8 + 0.1177 * xm3*xm8 + 0.1032 * xm4*xm8 +
3368 0.09073 * xm*xm4*xm8 + 0.07987 * xm6*xm8 + 0.07040 * xm*xm6*xm8 + 0.06210 * xm8*xm8 );
3369 else n82fr = -( -15961 * x5 + 1003 * x5*Lx3 - 2627 * x5*Lx2 - 29962 * x5*Lx -
3370 11683 * x4 + 54.66 * x4*Lx3 - 2777 * x4*Lx2 - 17770 * x4*Lx -
3371 6481 * x3 - 40.68 * x3*Lx3 - 1439 * x3*Lx2 - 7906 * x3*Lx -
3372 2943 * x2 - 31.83 * x2*Lx3 - 612.6 * x2*Lx2 - 2770 * x2*Lx -
3373 929.8 * x - 19.80 * x*Lx3 - 174.7 * x*Lx2 - 658.4 * x*Lx );
3382 CWbsgArrayNNLO[6] = -( n72ct + cd72ct * lstmu + ce72ct * lstmu * lstmu) * abssu2
3383 +( n72fr + cd72fr * lstmu + ce72fr * lstmu * lstmu) * susd
3385 CWbsgArrayNNLO[7] = -( n82ct + cd82ct * lstmu + ce82ct * lstmu * lstmu) * abssu2
3386 +( n82fr + cd82fr * lstmu + ce82fr * lstmu * lstmu) * susd
3397 std::stringstream out;
3399 throw std::runtime_error(
"order" + out.str() +
"not implemeted");
3425 std::stringstream out;
3427 throw std::runtime_error(
"order" + out.str() +
"not implemeted");