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;