singlekaon_xsec Class Reference
#include <singlekaon_xsec.h>
Public Member Functions | |
| void | init (double Etot, int type, int reac) |
| double | diffxsec (double Tlep, double Tkaon, double theta, double phikq) |
| double | Amatrix_NN (double theta, double phikq) |
| double | Amatrix_NP (double theta, double phikq) |
| double | Amatrix_PP (double theta, double phikq) |
Public Attributes | |
| double | threshold |
Private Attributes | |
| double | pi |
| double | amLam |
| double | am |
| double | amEta |
| double | aml |
| double | amSig |
| double | amk |
| double | ampi |
| double | Vus |
| int | ilep |
| int | ik |
| double | Enu |
| double | Ekaon |
| double | pkvec |
| double | GeVtocm |
| double | fpi |
| double | d |
| double | f |
| double | g |
| double | amup |
| double | amun |
| double | Fm1 |
| double | Fm2 |
| double | Elep |
| double | alepvec |
| double | aqvec |
| double | angkq |
| double | aq0 |
Detailed Description
Definition at line 6 of file singlekaon_xsec.h.
Member Function Documentation
◆ Amatrix_NN()
| double singlekaon_xsec::Amatrix_NN | ( | double | theta, |
| double | phikq ) |
Definition at line 129 of file singlekaon_xsec.cxx.
129 {
130
131 double sol = 0.;
132
133 double akk1=0., zdotq=0., qdotpk=0., akcrosk1=0., qcrospk=0.;
134 double zdotpk=0., azpk=0., aqkaon=0., akpk=0., apkk1=0.;
135 double C1=0., C2=0., C3=0., C4=0., /*C5=0., C6=0.,*/ C7=0., C8=0., C9=0.;
136 double aq2=0., gform=0., con=0., t1=0., t2=0., t3=0., /*t4=0.,*/ t5=0., t6=0.;
137
146 akpk=azpk + aqkaon/2.0;
147 apkk1=azpk - aqkaon/2.0;
148
153 //C5=d+3.*f;
154 //C6=1.0/(am*am+amk*amk-2.0*am*Ekaon-amLam*amLam);
159 gform=1.0/pow(1.0-aq2/(1.0*1.0),4);
160
162
163 t1=1.0;
164 t2=1.0;
165 t3=1.0;
166 //t4=1.0;
167 t5=1.0;
168 t6=1.0;
169
170 sol = con*gform*(4.*am*(-2.*(akk1*akk1)*
177 1.*apkk1)*C1*(C4*C4)*t3) -
179 (-2.*t2 +
180 C1*(C4*C4)*
189 Ekaon*(3.*(t2*t2) -
191 t2*t3 +
192 4.*(C1*C1)*(C4*C4*C4*C4)*
198 *(t3*t3))) -
199 8.*C1*C4*Fm1*t3*
203 (2.*t1 +
204 C1*
206 am*
208 2.*akpk*Fm1 +
209 (-2.*apkk1 +
213 (-2. +
216 amSig*
217 (-4.*Ekaon +
218 2.*(akpk - 1.*apkk1)*Fm1 +
220 ))*t3) +
221 2.*Ekaon*
223 (-1.*akpk + apkk1 -
225 Fm1*t3 -
227 amSig*
228 (t1 +
232 (C2*C7*t5 + C8*C9*t6))) +
237 C1*(C4*C4)*t3*
238 (2.*C1*(1. + (C4*C4))*
247 (t2 +
248 2.*C1*(C4*C4)*
250 ))) +
255 C4*
259 (-1.*C3*t2 -
260 2.*C1*C4*
261 (2. +
262 C4*
266 C1*C4*t3*
267 (-2.*C1*C4*
271 2.*am*
273 Fm1 -
284 C3*C4*
288 2.*C1*(C4*C4)*
293 t3) -
294 1.*
297 *(C2*C7*t5 + C8*C9*t6))) +
298 apkk1*(-2.*((-1.+C4)*(-1.+C4))*(t1*t1) +
299 2.*C1*C4*
302 2.*am*
306 Ekaon) -
309 2.*am*
311 Fm1)*t1*t3 +
312 C1*C4*t3*
313 (2.*C1*C4*
318 C4 -
319 1.*
323 2.*am*
325 2.*
327 C4 +
331 (-1. +
336 *(C2*C7*t5 + C8*C9*t6))) +
337 2.*(apkk1*apkk1)*C1*C4*t3*
339 (C2*C7*t5 + C8*C9*t6) +
342 2.*am*
343 (t1 - 1.*C4*t1 + C2*C7*Ekaon*t5 +
344 C8*C9*Ekaon*t6)))) +
345 2.*(akpk*akpk)*C1*C4*t3*
347 (C3*C4*(2.*t1 + t2 -
348 4.*C1*(C4*C4)*
350 + Fm1*
351 (4.*C1*C4*
353 t3 + C2*C7*t5 + C8*C9*t6)) -
355 (C2*C7*t5 + C8*C9*t6) +
358 2.*am*
359 ((1. + C4)*t1 -
360 1.*Ekaon*(C2*C7*t5 + C8*C9*t6))))) +
364 (C3*C4*(2.*t1 + t2 -
365 4.*C1*(C4*C4)*
367 - 1.*Fm1*(C2*C7*t5 + C8*C9*t6)) +
373 (2.*((-1.+C4)*(-1.+C4))*Ekaon +
375 (t3*t3) +
377 (2.*C1*((-1.+C4)*(-1.+C4))*C4*Ekaon*t3 +
380 Fm1*
381 (-2.*(-1. + C4)*t1 +
383 1.*C3*
384 (2.*Enu*
385 (t1*t2 +
386 C1*(C4*C4)*
389 (2.*t1 + t2)*t3 +
390 2.*(C1*C1)*(C4*C4*C4*C4)*
394 2.*am*
399 2.*apkk1*
400 (2.*t1 + t2 -
401 4.*C1*(C4*C4)*
403 t3))) -
404 1.*(C2*C7*t5 + C8*C9*t6)*
405 (-3.*C4*t1 +
407 + 2.*amSig*C1*C4*t3*
408 (2.*(-1.*((-1.+C4)*(-1.+C4)) +
411 (((-1.+C4)*(-1.+C4)) +
414 (C2*C7*t5 + C8*C9*t6)) +
415 2.*am*C1*C4*t3*
416 (-2.*
417 (((-1.+C4)*(-1.+C4)) +
420 Fm1*t3 -
424 (((-1.+C4)*(-1.+C4)) +
425 (4.*amSig*C4 -
428 (apkk1*Fm1 +
430 (C2*C7*t5 + C8*C9*t6)) +
432 (-3.*(C2*C7*t5 + C8*C9*t6) +
433 2.*Fm1*
434 ((2. + C4)*t1 +
435 Enu*(C2*C7*t5 + C8*C9*t6)))) -
436 1.*Elep*(2.*((1.+C4)*(1.+C4))*(t1*t1) +
439 (t3*t3) +
442 Ekaon*(t3*t3) +
447 (t2 +
449 t3)) -
450 2.*t1*
451 (2.*C1*C4*
456 (apkk1 - 1.*apkk1*C4 +
460 t3 +
461 C3*
462 (t2 -
464 t3))) +
465 (C2*C7*t5 + C8*C9*t6)*
467 (amSig*C1*C4*t3 + C2*C7*t5 + C8*C9*t6)
468 + 2.*apkk1*
469 (amSig*C1*C4*(1. + C4)*t3 + C2*C7*t5 +
470 C8*C9*t6)) +
473 (2.*C3*C4*
475 Fm1*(C2*C7*t5 + C8*C9*t6)) -
476 2.*
478 amSig*(C2*C7*t5 + C8*C9*t6) -
479 1.*apkk1*Fm1*(C2*C7*t5 + C8*C9*t6)))
480 + am*C1*C4*t3*
483 4.*C3*C4*Ekaon*
485 amSig*C3*
488 (C2*C7*t5 + C8*C9*t6)) -
489 2.*
491 t3 -
493 (C2*C7*t5 + C8*C9*t6) +
495 (4.*C1*C4*
496 (amSig*(1. + (C4*C4)) +
497 ((1.+C4)*(1.+C4))*Ekaon)*t3 +
498 C2*C7*t5 + C8*C9*t6)))) +
499 2.*apkk1*
501 (t3*t3) -
503 (t3*t3) +
504 (C2*C7*t5 + C8*C9*t6)*
505 (C4*(2.*t1 +
508 (C2*C7*t5 + C8*C9*t6)) +
509 amSig*C1*C4*t3*
513 C4*
516 1.*Enu*(C2*C7*t5 + C8*C9*t6))) +
517 am*C1*C4*t3*
519 C4*
522 1.*Enu*(C2*C7*t5 + C8*C9*t6)) +
523 Enu*
524 (C2*C7*t5 + C8*C9*t6 +
525 2.*Fm1*
526 (t1 -
527 1.*Ekaon*(C2*C7*t5 + C8*C9*t6))))))
530 (C4*C4)*
531 (-1. +
537 ) - 1.*Ekaon*
540 (-2.*C3*(C4*C4)*Ekaon -
544 4.*Fm1*
547 2.*Ekaon*
551 (3.*Ekaon*(t2*t2) -
552 4.*C1*(C4*C4)*
555 4.*(C1*C1)*(C4*C4*C4*C4)*
560 (t3*t3)) +
564 4.*(akpk*akpk)*C1*C4*Ekaon*
566 4.*(apkk1*apkk1)*C1*C4*Ekaon*
570 (C4*C4)*
573 Enu) +
576 4.*amSig*
577 (1. +
578 (C4*C4)*
579 (-1. +
589 2.*C1*
590 (((1.+C4)*(1.+C4))*Elep -
591 1.*((-1.+C4)*(-1.+C4))*Enu -
592 1.*amSig*
593 (-1. + (C4*C4) +
595 Fm1))*t3 +
600 ) - 1.*Ekaon*
601 (-4.*Enu*
602 ((1. - 3.*C4)*Fm1*t1 +
605 (4. +
606 (C4*C4)*
607 (-4. +
611 Fm1*
618 4.*Elep*
621 Fm1*
622 (t1 + 3.*C4*t1 -
624 2.*apkk1*
625 (-2.*(1. + C4)*Fm1*t1 +
627 C1*C4*
631 2.*C1*((1.+C4)*(1.+C4))*Ekaon*t3 +
636 2.*akpk*(C1*C4*
639 2.*Fm1*
640 (t1 - 1.*C4*t1 +
642 C4*(2.*C1*((-1.+C4)*(-1.+C4))*Ekaon*t3 +
646 C1*(C4*C4)*
647 (4.*(-1. + apkk1*C3)*Ekaon +
650 )))) +
656 (t2*t2) +
658 (akpk - 1.*apkk1)*C1*(C4*C4)*t2*t3 -
659 16.*amSig*((akpk - 1.*apkk1)*(akpk - 1.*apkk1))*(C1*C1)*(C4*C4*C4*C4)*
661 8.*C3*(-1.*Ekaon*t1*t2 +
662 C1*(C4*C4)*
664 akpk*amSig*(2.*t1 + t2))*t3 -
665 1.*(C1*C1)*(C4*C4*C4*C4)*
669 (t3*t3)) -
670 1.*(C2*C7*t5 + C8*C9*t6)*
672 - 4.*amSig*C1*C4*t3*
675 (C2*C7*t5 + C8*C9*t6)) +
677 (C4*Fm1*
678 (6.*t1 +
679 amSig*C1*
682 Fm1)*t3) -
684 (C2*C7*t5 + C8*C9*t6))) -
685 1.*Elep*(8.*((1.+C4)*(1.+C4))*(t1*t1) +
687 2.*apkk1)*(C3*C3)*(t2*t2) -
688 16.*C1*C4*
694 + 2.*Fm1*(C2*C7*t5 + C8*C9*t6)) -
700 2.*amSig*C1*C4*(1. + C4)*t3*
701 (C2*C7*t5 + C8*C9*t6) +
702 ((C2*C7*t5 + C8*C9*t6)*(C2*C7*t5 + C8*C9*t6)) +
704 (2.*apkk1*(C3*C3)*C4*
706 Fm1*
708 1.*C2*C7*t5 - 1.*C8*C9*t6))) +
713 (C3*C3)*(t2 + 2.*apkk1*C1*(C4*C4)*t3)) +
715 t3 - 1.*Fm1*(C2*C7*t5 + C8*C9*t6)))) -
716 4.*(8.*(akpk*akpk)*amSig*(C1*C1)*(C4*C4)*
719 (((1.+C4)*(1.+C4))*Ekaon -
723 ((C2*C7*t5 + C8*C9*t6)*
726 2.*apkk1*C1*C4*Fm1*t3*
728 1.*Ekaon*(C2*C7*t5 + C8*C9*t6))) +
729 2.*amSig*C1*C4*t3*
732 apkk1*
733 (-2.*
735 t1 +
737 (((1.+C4)*(1.+C4)) -
739 (1. + C4)*Ekaon*(C2*C7*t5 + C8*C9*t6)))
740 - 2.*akpk*C1*C4*t3*
742 (((-1.+C4)*(-1.+C4))*Ekaon +
745 (2.*(1. + C4)*
747 Ekaon*(C2*C7*t5 + C8*C9*t6)) +
748 amSig*
749 (2.*(((-1.+C4)*(-1.+C4)) -
752 (((-1.+C4)*(-1.+C4)) +
755 (C2*C7*t5 + C8*C9*t6)))) +
756 Enu*(8.*((-1.+C4)*(-1.+C4))*(t1*t1) -
761 + 8.*t1*
762 (2.*C1*C4*
768 ((C2*C7*t5 + C8*C9*t6)*(C2*C7*t5 + C8*C9*t6))) +
773 (C3*C3)*(t2 + 2.*apkk1*C1*(C4*C4)*t3)) +
774 2.*(C1*((-1.+C4)*(-1.+C4))*C4*t3 -
775 1.*Fm1*(C2*C7*t5 + C8*C9*t6))) -
780 (-8.*C1*C3*(C4*C4)*t3*
782 (C3*C3)*
783 (3.*(t2*t2) + 8.*apkk1*C1*(C4*C4)*t2*t3 +
785 + 4.*C1*C4*Fm1*t3*
787 C8*C9*t6)) +
789 (t3*t3) -
790 2.*amSig*C1*(-1. + C4)*C4*t3*
791 (C2*C7*t5 + C8*C9*t6) +
792 (C2*C7*t5 + C8*C9*t6)*
793 (4.*apkk1*C1*C4*Fm1*t3 + C2*C7*t5 +
794 C8*C9*t6)))) +
796 (t2*t2) -
797 2.*(-2.*(-1. + (C4*C4))*(t1*t1) -
798 2.*C1*C4*
799 (2.*akpk*
804 2.*apkk1*
805 (((1.+C4)*(1.+C4)) -
806 1.*
807 (amSig*(1. + C4) +
811 - 1.*Ekaon*
812 (4.*((1.+C4)*(1.+C4))*Elep +
813 4.*((-1.+C4)*(-1.+C4))*Enu +
818 t1*t3 +
819 C1*C4*t3*
821 (1. +
822 (C4*C4)*
823 (-1. +
825 (-2. +
829 )*t3 +
830 2.*(akpk*akpk)*C4*
835 (-1.*apkk1*t2 +
836 2.*Ekaon*
838 1.*Enu*t2 +
840 2.*apkk1*Ekaon*
842 t3 -
846 (C2*C7*t5 + C8*C9*t6)) -
848 (2.*(apkk1*apkk1)*(C3*C3)*C4*t2 +
849 apkk1*
850 (2.*(C3*C3)*C4*
855 1.*Fm1*(C2*C7*t5 + C8*C9*t6)) +
856 Ekaon*
861 t3 +
862 (1. +
864 (C2*C7*t5 + C8*C9*t6))) +
865 akpk*
868 2.*
870 Fm1*t3 +
874 (((-1.+C4)*(-1.+C4)) +
875 2.*
878 )*t3 -
879 1.*
880 (4.*apkk1*Fm1 +
881 Ekaon*
884 ) +
886 (C3*
888 + (C3*C3)*C4*
890 (t2 +
891 8.*C1*(C4*C4)*
892 (apkk1 -
894 *t3) +
895 2.*
896 (2.*apkk1*t2 +
901 Fm1*
903 C2*C7*t5 + C8*C9*t6))) +
906 (t2 +
907 C1*(C4*C4)*
911 t3) +
912 2.*
913 (-4.*C1*C4*Ekaon*
914 (((1.+C4)*(1.+C4))*Elep -
915 1.*((-1.+C4)*(-1.+C4))*Enu)*t3 +
917 2.*apkk1*C1*C4*
918 (((1.+C4)*(1.+C4)) +
922 (-1. + (C4*C4) +
925 2.*amSig*C1*C4*
926 (-1.*((-1.+C4)*(-1.+C4))*Enu +
927 4.*apkk1*C4*Fm1 -
928 1.*Ekaon*
929 (-1. + (C4*C4) +
930 2.*
932 Fm1) +
933 Elep*
935 t3 +
939 (C2*C7*t5 + C8*C9*t6)) +
941 (2.*C3*C4*
942 (-1.*t2 +
943 2.*C1*(C4*C4)*
947 (C3*C3)*C4*
948 (4.*(apkk1*apkk1)*C1*(C4*C4)*t3 +
949 apkk1*
950 (t2 -
954 (-2.*Ekaon*t2 +
957 2.*Fm1*
958 (C1*C4*
960 amSig*
961 (6. -
963 + 2.*Ekaon*
964 (-3. +
966 + C2*C7*t5 + C8*C9*t6)))))))));
967
968 if (sol <= 0.0)
969 std::cout << "Check Matrix and def. of Scalar Products" << std::endl;
970
971 return sol;
972}
References alepvec, am, amEta, amk, aml, ampi, amSig, angkq, aq0, aqvec, d, Ekaon, Elep, Enu, f, Fm1, fpi, g, pkvec, and Vus.
Referenced by diffxsec().
◆ Amatrix_NP()
| double singlekaon_xsec::Amatrix_NP | ( | double | theta, |
| double | phikq ) |
Definition at line 975 of file singlekaon_xsec.cxx.
975 {
976
977 double sol = 0.;
978
979 double akk1=0., zdotq=0., qdotpk=0., akcrosk1=0., qcrospk=0.;
980 double zdotpk=0., azpk=0., aqkaon=0., akpk=0., apkk1=0.;
981 double C1=0., C2=0., C3=0., C4=0., C5=0., C6=0., C7=0./*, C8=0., C9=0.*/;
982 double aq2=0., gform=0., con=0., t1=0., t2=0., t3=0., t4=0., t5=0./*, t6=0.*/;
983
992 akpk=azpk + aqkaon/2.0;
993 apkk1=azpk - aqkaon/2.0;
1001 //C8=am/(aml*aml-2.0*akk1+amk*amk-2.*aqkaon-amEta*amEta);
1002 //C9=d - 3.*f;
1004 gform=1.0/pow(1.0-aq2/(1.1*1.1),4);
1005
1007
1008 t1=1.0;
1009 t2=1.0; // !Full Term
1010 t3=1.0;
1011 t4=1.0;
1012 t5=1.0;
1013 //t6=1.0;
1014
1015 sol = con*gform*(0.4444444444444444*am*
1033 72.*akpk*amSig*C1*C4*t1*t3 -
1034 72.*amSig*apkk1*C1*C4*t1*t3 -
1035 72.*akpk*amSig*C1*C2*C4*t1*t3 +
1036 72.*amSig*apkk1*C1*C2*C4*t1*t3 +
1037 72.*akpk*amSig*C1*(C4*C4)*t1*t3 -
1038 72.*amSig*apkk1*C1*(C4*C4)*t1*t3 +
1039 72.*akpk*amSig*C1*C2*(C4*C4)*t1*t3 +
1040 72.*amSig*apkk1*C1*C2*(C4*C4)*t1*t3 -
1073 (C4*C4)*t2*t3 -
1082 Elep*t2*t3 -
1095 (C4*C4*C4*C4)*(t3*t3) -
1098 (C4*C4*C4*C4)*(t3*t3) +
1101 (C4*C4*C4*C4)*(t3*t3) +
1103 (C4*C4*C4*C4)*(t3*t3) -
1106 (C3*C3)*(C4*C4*C4*C4)*(t3*t3) +
1116 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
1119 Ekaon*(t3*t3) +
1122 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
1124 Ekaon*(t3*t3) -
1126 Ekaon*(t3*t3) -
1128 Ekaon*(t3*t3) -
1138 (C3*C3)*(C4*C4*C4*C4)*Elep*(t3*t3) +
1141 (C4*C4*C4*C4)*Elep*(t3*t3) +
1143 (C4*C4*C4*C4)*Elep*(t3*t3) +
1153 (C4*C4*C4*C4)*Enu*(t3*t3) +
1156 (C4*C4*C4*C4)*Enu*(t3*t3) -
1159 (C4*C4*C4*C4)*Enu*(t3*t3) -
1161 (C4*C4*C4*C4)*Enu*(t3*t3) +
1168 Fm1*(t3*t3) +
1174 Fm1*(t3*t3) -
1210 72.*amLam*apkk1*C5*C6*t1*t4 -
1211 72.*akpk*amLam*C2*C5*C6*t1*t4 +
1212 72.*amLam*apkk1*C2*C5*C6*t1*t4 -
1213 24.*akpk*amLam*(C5*C5)*C6*t1*t4 +
1214 24.*amLam*apkk1*(C5*C5)*C6*t1*t4 -
1215 24.*akpk*amLam*C2*(C5*C5)*C6*t1*t4 -
1216 24.*amLam*apkk1*C2*(C5*C5)*C6*t1*t4 +
1279 C6*t3*t4 +
1282 (C5*C5)*C6*t3*t4 +
1285 (C5*C5*C5)*C6*t3*t4 +
1289 C6*t3*t4 +
1291 (C5*C5*C5)*C6*t3*t4 -
1294 (C5*C5*C5)*C6*t3*t4 +
1296 (C5*C5*C5)*C6*t3*t4 -
1299 (C5*C5*C5)*C6*t3*t4 -
1301 (C5*C5*C5)*C6*t3*t4 +
1303 C6*t3*t4 -
1305 (C5*C5*C5)*C6*t3*t4 +
1307 C6*t3*t4 +
1309 (C5*C5*C5)*C6*t3*t4 +
1311 (C5*C5*C5)*C6*t3*t4 +
1317 Ekaon*t3*t4 +
1320 Ekaon*t3*t4 +
1323 Ekaon*t3*t4 +
1325 Ekaon*t3*t4 +
1327 C6*Ekaon*t3*t4 +
1330 Ekaon*t3*t4 +
1332 Ekaon*t3*t4 -
1335 (C5*C5*C5)*C6*Ekaon*t3*t4 -
1337 Ekaon*t3*t4 -
1339 C6*Ekaon*t3*t4 -
1341 (C5*C5*C5)*C6*Ekaon*t3*t4 +
1343 C6*Ekaon*t3*t4 +
1345 (C5*C5*C5)*C6*Ekaon*t3*t4 -
1355 C6*Elep*t3*t4 +
1358 (C4*C4)*(C5*C5*C5)*C6*Elep*t3*t4 -
1360 C6*Elep*t3*t4 +
1362 Elep*t3*t4 -
1364 (C5*C5*C5)*C6*Elep*t3*t4 +
1375 (C4*C4)*(C5*C5)*C6*Enu*t3*t4 +
1378 Enu*t3*t4 +
1381 Enu*t3*t4 -
1383 (C5*C5*C5)*C6*Enu*t3*t4 +
1385 C6*Enu*t3*t4 -
1388 C3*(C4*C4)*(C5*C5*C5)*C6*Enu*t3*t4 +
1393 Fm1*t3*t4 -
1400 C6*Fm1*t3*t4 +
1403 Fm1*t3*t4 +
1406 Fm1*t3*t4 -
1409 Fm1*t3*t4 -
1418 Fm1*t3*t4 -
1421 Fm1*t3*t4 +
1426 Fm2*t3*t4 -
1433 C6*Fm2*t3*t4 +
1436 Fm2*t3*t4 +
1439 Fm2*t3*t4 -
1442 Fm2*t3*t4 +
1451 Fm2*t3*t4 -
1454 Fm2*t3*t4 -
1460 Fm2*t3*t4 +
1463 Fm2*t3*t4 -
1465 Fm2*t3*t4 -
1467 Fm2*t3*t4 -
1470 Fm2*t3*t4 +
1481 Fm2*t3*t4 -
1506 (C6*C6)*(t4*t4) -
1517 Ekaon*(t4*t4) -
1522 Ekaon*(t4*t4) +
1526 Ekaon*(t4*t4) -
1540 (C6*C6)*Elep*(t4*t4) +
1553 Enu*(t4*t4) +
1562 Fm2*(t4*t4) -
1572 Fm2*(t4*t4) -
1603 (1. +
1604 (C4*C4)*
1611 6.*C1*C4*C5*C6*
1616 C4*
1622 (C5*C5)*(C6*C6)*
1629 1.*Ekaon*
1631 ((-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)*(-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)) +
1633 (-9.*(C1*C1)*(C4*C4)*
1634 (2.*C3*(C4*C4)*Ekaon +
1638 (t3*t3) +
1639 6.*C1*C4*C5*C6*
1640 (C3*C4*C5*
1641 (Ekaon*
1643 2.*(akpk - 1.*apkk1)*C5*
1646 (C5*C5)*(C6*C6)*
1647 (-2.*(C5*C5*C5)*Ekaon -
1651 ) +
1653 (apkk1*
1656 akpk*
1659 - 2.*Ekaon*
1660 (3.*C1*C4*
1663 C5*C6*
1664 (-1.*(-3. + C5)*Elep +
1680 12.*apkk1*C4*Fm1 +
1685 4.*amSig*
1686 (1. +
1687 (C4*C4)*
1694 - 1.*Ekaon*
1695 (4. +
1696 (C4*C4)*
1705 2.*Ekaon*
1706 (2.*(akpk*akpk)*
1708 2.*(apkk1*apkk1)*
1710 apkk1*
1711 (2. + 4.*C4 +
1712 (C4*C4)*
1718 1.*akpk*
1719 (2. - 4.*C4 +
1720 (C4*C4)*
1723 (-1. + apkk1*C3 +
1727 2.*Ekaon*
1730 2.*
1731 (((1.+C4)*(1.+C4))*Elep -
1732 1.*((-1.+C4)*(-1.+C4))*Enu) +
1733 amSig*
1734 (-2. +
1743 ))*(t3*t3) +
1744 C5*C6*t4*
1746 Fm2*t1 +
1747 24.*C2*
1748 ((akpk - 1.*apkk1)*Fm2 +
1750 t1 +
1752 (4.*Ekaon*
1753 (akpk +
1756 t4 +
1757 2.*(C5*C5*C5)*C6*
1760 + 2.*Ekaon*
1761 (-1.*akpk + apkk1 +
1765 (-2.*
1766 (3.*akpk - 3.*apkk1 +
1771 1.*Enu)) +
1772 2.*Ekaon*
1773 (2.*(akpk*akpk) +
1774 2.*apkk1*
1777 (apkk1 +
1778 2.*Ekaon*
1780 ) -
1781 1.*akpk*
1783 4.*
1785 )))*t4 -
1786 4.*(C5*C5)*
1789 + C6*
1795 ))*t4 +
1796 3.*akpk*
1799 + C5*
1802 (Enu*
1805 3.*amLam*
1808 Elep*
1812 (4.*amLam*
1816 3.*
1821 Elep*
1824 1.*Ekaon*
1825 (24.*t1 +
1826 C6*
1829 2.*apkk1*
1837 6.*C1*C4*t3*
1841 12.*C2*
1846 12.*akpk*C5*C6*Ekaon*t4 +
1848 12.*apkk1*C5*C6*Ekaon*t4 +
1849 4.*akpk*(C5*C5)*C6*Ekaon*t4 +
1850 4.*apkk1*(C5*C5)*C6*Ekaon*t4 +
1899 t4 +
1901 t4 -
1906 t4 +
1911 t4 -
1936 t4 -
1946 t4 +
1948 Enu*t4 +
1950 Enu*t4 -
1960 6.*akpk*
1967 t4 -
1968 2.*C5*C6*
1970 2.*C5*Enu +
1973 t4 +
1976 1.*(C5*C5*C5)*C6*
1978 1.*C3*(C5*C5)*C6*
1982 t4)) +
1983 Ekaon*
1988 12.*apkk1*C5*C6*t4 +
1990 4.*apkk1*(C5*C5)*C6*t4 +
1998 t4 +
2000 Enu*t4 -
2004 2.*akpk*
2006 2.*C5*(3. + C5)*C6*t4 +
2008 (3.*(C3*C3)*t2 - 2.*(C5*C5*C5)*C6*t4 +
2009 2.*C3*(C5*C5)*(-1. + 2.*apkk1*C5)*C6*
2010 t4)))))) -
2069 2.*am*(36.*(-1. + (C2*C2))*(t1*t1) +
2096 (t3*t3) -
2099 (t3*t3) +
2101 (t3*t3) +
2103 (t3*t3) -
2106 (t3*t3) -
2108 (t3*t3) +
2110 (t3*t3) -
2112 (t3*t3) +
2114 (t3*t3) -
2116 (t3*t3) -
2118 (t3*t3) +
2121 (t3*t3) +
2123 (t3*t3) -
2125 (t3*t3) -
2128 (t3*t3) -
2130 (t3*t3) -
2132 (t3*t3) -
2135 (t3*t3) +
2146 (t3*t3) +
2148 (t3*t3) +
2168 t3*t4 +
2177 t3*t4 -
2183 t3*t4 -
2197 Fm1*t3*t4 -
2200 Fm1*t3*t4 -
2203 t3*t4 -
2206 t3*t4 +
2208 t3*t4 +
2217 t3*t4 +
2220 Fm1*t3*t4 +
2223 Fm1*t3*t4 +
2225 t3*t4 +
2227 t3*t4 +
2230 Fm1*t3*t4 +
2232 t3*t4 +
2234 t3*t4 -
2236 t3*t4 +
2238 t3*t4 +
2240 t3*t4 +
2242 t3*t4 -
2245 t3*t4 -
2248 Fm1*t3*t4 +
2250 t3*t4 +
2252 t3*t4 +
2254 t3*t4 -
2256 t3*t4 +
2259 Fm2*t3*t4 +
2262 Fm2*t3*t4 -
2266 t3*t4 +
2268 t3*t4 -
2277 t3*t4 +
2280 Fm2*t3*t4 -
2283 Fm2*t3*t4 -
2285 t3*t4 -
2287 t3*t4 +
2290 Fm2*t3*t4 +
2293 Fm2*t3*t4 -
2295 t3*t4 -
2297 t3*t4 +
2299 t3*t4 -
2301 t3*t4 -
2304 t3*t4 -
2307 Fm2*t3*t4 -
2309 t3*t4 -
2311 t3*t4 +
2313 t3*t4 -
2315 t3*t4 -
2318 t3*t4 +
2320 Fm2*t3*t4 +
2322 Fm2*t3*t4 +
2324 Fm2*t3*t4 +
2331 t3*t4 +
2333 t3*t4 +
2335 t3*t4 +
2340 Fm2*t3*t4 +
2369 (t4*t4) +
2371 (t4*t4) +
2373 (t4*t4) +
2376 (t4*t4) +
2379 (t4*t4) -
2381 (t4*t4) +
2383 (t4*t4) +
2385 (t4*t4) +
2387 (t4*t4) +
2390 (t4*t4) +
2392 (t4*t4) +
2394 (t4*t4) -
2397 (t4*t4) -
2399 (t4*t4) +
2401 (t4*t4) -
2405 (t4*t4) +
2407 (t4*t4) +
2412 (t4*t4) +
2414 (t4*t4) +
2416 (t4*t4) +
2418 (t4*t4) +
2421 (t4*t4) +
2466 2.*akpk*
2467 (3.*C1*C4*
2472 )*t3 +
2473 C5*C6*
2477 C2*
2480 2.*apkk1*
2481 (3.*C1*C4*
2484 C2*
2487 1.*C5*C6*
2491 )*t4)) +
2497 (t3*t3) + 3.*apkk1*C3*(C5*C5*C5)*C6*t2*t4 +
2505 ((-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)*(-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)) +
2506 akpk*(3.*C1*C3*(C4*C4)*t3 - 1.*(C5*C5*C5)*C6*t4)*
2510 (9.*(C1*C1)*(C3*C3)*(C4*C4*C4*C4)*
2514 (t3*t3) +
2515 (C5*C5*C5)*C6*t4*
2516 (-6.*C3*t2 -
2517 1.*(C5*C5*C5)*C6*
2521 t4) +
2522 3.*C1*C3*(C4*C4)*t3*
2523 (6.*C3*t2 +
2524 (C5*C5*C5)*C6*
2525 (2.*apkk1 +
2528 amLam*
2530 )*t4))) +
2531 36.*akpk*C1*C2*C4*C7*Ekaon*t3*t5 -
2533 36.*apkk1*C1*C2*C4*C7*Ekaon*t3*t5 -
2534 36.*akpk*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 -
2535 36.*apkk1*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 +
2540 144.*akpk*apkk1*C1*C2*C4*C7*Fm1*t3*t5 -
2547 36.*akpk*C2*C5*C6*C7*Ekaon*t4*t5 -
2549 36.*apkk1*C2*C5*C6*C7*Ekaon*t4*t5 +
2550 12.*akpk*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
2551 12.*apkk1*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
2556 48.*akpk*apkk1*C2*C5*C6*C7*Fm2*t4*t5 -
2564 (t3*t3) +
2567 (t3*t3) +
2571 (t3*t3) -
2573 (t3*t3) +
2575 (t3*t3) -
2577 (t3*t3) -
2579 (t3*t3) +
2581 (t3*t3) +
2583 (t3*t3) -
2585 (t3*t3) -
2587 (t3*t3) +
2588 9.*(C3*C3)*
2590 2.*(akpk*akpk)*C1*(C4*C4)*t2*t3 -
2592 4.*akpk*apkk1*C1*(C4*C4)*t2*t3 +
2594 2.*(apkk1*apkk1)*C1*(C4*C4)*t2*t3 +
2600 (t3*t3) -
2603 (t3*t3) -
2605 (t3*t3) +
2608 (t3*t3) +
2612 (C4*C4*C4*C4)*(t3*t3)) -
2613 2.*Elep*
2614 (amSig*(akpk - 1.*apkk1)*C1*(C4*C4)*t3*
2616 Enu*
2619 (C1*C1)*(C4*C4*C4*C4)*(t3*t3))) -
2620 2.*Ekaon*
2621 (-1.*amSig*(akpk - 1.*apkk1)*C1*(C4*C4)*
2622 t2*t3 +
2623 Elep*
2624 ((t2*t2) +
2626 (C4*C4)*t2*t3 +
2628 (akpk - 1.*apkk1)*(C1*C1)*(C4*C4*C4*C4)*(t3*t3)
2629 ) -
2630 1.*Enu*
2631 ((t2*t2) +
2633 (C4*C4)*t2*t3 +
2635 (akpk - 1.*apkk1)*(C1*C1)*(C4*C4*C4*C4)*(t3*t3)
2636 ))) +
2640 t4 -
2642 t4 -
2644 t4 -
2646 t4 -
2648 t3*t4 -
2650 t4 -
2652 t4 -
2654 t4 -
2666 t4 +
2669 t3*t4 -
2671 t4 +
2673 t3*t4 +
2675 t3*t4 -
2677 t4 -
2679 Fm2*t3*t4 +
2681 t3*t4 +
2683 t3*t4 +
2685 t3*t4 -
2687 t4 -
2689 t4 -
2691 t3*t4 -
2702 (t4*t4) +
2704 (t4*t4) -
2706 (t4*t4) -
2708 (t4*t4) +
2710 (t4*t4) +
2712 (t4*t4) +
2715 (t4*t4) +
2718 (t4*t4) -
2720 (t4*t4) +
2723 (t4*t4) +
2725 (t4*t4) -
2727 (t4*t4) -
2729 (t4*t4) -
2734 (t4*t4) -
2741 (t4*t4) -
2744 (t4*t4) -
2746 (t4*t4) -
2748 (t4*t4) +
2750 (t4*t4) +
2752 (t4*t4) -
2754 (t4*t4) -
2756 (t4*t4) +
2757 3.*C3*
2758 (2.*(akpk*akpk)*(C5*C5*C5)*C6*
2761 (3.*C1*(C4*C4)*t3 +
2762 (C5*C5)*C6*
2764 2.*
2766 (t3*t3) +
2767 (C5*C5)*C6*
2768 ((apkk1*apkk1)*C5 -
2770 apkk1*C5*
2773 t4 +
2774 C1*(C4*C4)*Ekaon*t3*
2778 t4 +
2780 t4)) -
2782 (6.*(C1*C1)*(C4*C4*C4*C4)*
2786 (C5*C5)*C6*
2787 (-2. + apkk1*C5 +
2790 + C1*(C4*C4)*t3*
2791 (6.*t2 +
2792 (C5*C5)*C6*
2793 (-4.*(apkk1*apkk1)*C5 +
2798 amLam*
2800 2.*amSig*
2803 akpk*
2805 (C5*C5)*C6*
2811 t2*t4 -
2812 2.*C1*(C4*C4)*t3*
2813 (3.*t2 -
2814 2.*(C5*C5)*C6*
2818 (1. - 2.*apkk1*C5 +
2822 *t4))) +
2823 18.*C1*C2*C4*C7*Ekaon*t3*t5 -
2824 18.*akpk*C1*C2*C4*C7*Fm1*t3*t5 -
2826 18.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
2830 18.*C2*C5*C6*C7*Ekaon*t4*t5 -
2831 6.*akpk*C2*C5*C6*C7*Fm2*t4*t5 -
2833 6.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
2838 (9.*(C3*C3)*(Enu*(t2*t2) +
2839 2.*akpk*am*C1*(C4*C4)*t2*t3 +
2841 2.*akpk*amSig*C1*(C4*C4)*t2*t3 +
2843 2.*am*apkk1*C1*(C4*C4)*t2*t3 -
2844 2.*amSig*apkk1*C1*(C4*C4)*t2*t3 +
2861 (4.*t2 +
2862 C1*(C4*C4)*
2865 Elep*(-1.*(t2*t2) -
2870 Ekaon*
2871 (3.*(t2*t2) +
2872 4.*C1*(C4*C4)*
2874 (C1*C1)*(C4*C4*C4*C4)*
2881 - 6.*C3*(C5*C5*C5)*C6*
2900 2.*Ekaon*
2901 (apkk1 +
2902 Ekaon*
2908 amLam*
2909 (t2 +
2911 )) +
2917 (2.*t2 -
2918 1.*C1*(C4*C4)*
2919 (2.*apkk1 +
2922 amLam*
2924 )*t3)))*t4 +
2925 (C5*C5*C5*C5*C5*C5)*(C6*C6)*
2934 2.*am*
2943 + 2.*Ekaon*
2944 (apkk1 +
2945 2.*Ekaon*
2947 ))*(t4*t4)) -
2949 (3.*C1*C4*
2951 + C5*C6*
2954 (-9.*(C1*C1)*(C4*C4)*Fm1*
2959 amSig*
2960 (-4.*Ekaon +
2961 2.*(akpk - 1.*apkk1)*Fm1 +
2963 + am*
2965 2.*(akpk - 1.*apkk1)*Fm1 +
2966 4.*Ekaon*
2968 2.*amSig*
2971 C5*C6*Fm2*t4*
2972 (-1.*C5*C6*
2978 amLam*
2979 (2.*
2981 1.*apkk1*Fm2) +
2983 ) +
2984 am*
2986 2.*amLam*
2989 2.*
2990 ((akpk - 1.*apkk1)*Fm2 +
2991 2.*Ekaon*
2993 - 12.*C2*
2995 - 6.*C1*C4*t3*
2996 (C5*C6*
2998 1.*amSig*
3002 (-1.*Fm2 +
3003 Fm1*
3007 1.*amLam*
3010 amSig*
3011 (Fm2 +
3013 + am*
3014 (amSig*
3015 (-1.*Fm2 +
3016 Fm1*
3019 amLam*
3020 (-1.*Fm2 +
3021 Fm1*
3024 2.*
3026 Ekaon*
3027 (Fm2 +
3029 ))))*t4 +
3030 6.*C2*Fm1*
3032 + 2.*Ekaon*
3035 (6.*C2*t1 +
3036 (akpk - 1.*apkk1)*
3041 1.*apkk1*
3043 2.*Ekaon*
3047 ) -
3048 6.*am*
3050 C5*C6*Fm2*t4*
3052 2.*C2*C7*Ekaon*t5) +
3053 C1*C4*t3*
3057 + 2.*(Enu*(-36.*apkk1*(t1*t1) - 72.*apkk1*C2*(t1*t1) -
3104 (t3*t3) -
3108 (t3*t3) -
3112 (t3*t3) -
3119 (t3*t3) -
3121 (t3*t3) -
3125 (t3*t3) -
3127 (t3*t3) -
3129 (t3*t3) -
3131 (t3*t3) +
3170 t3*t4 +
3174 C6*t3*t4 -
3176 t3*t4 -
3181 t3*t4 +
3184 C6*t3*t4 +
3187 t3*t4 -
3189 t3*t4 -
3191 t3*t4 -
3197 Ekaon*t3*t4 +
3200 Ekaon*t3*t4 +
3202 t3*t4 +
3206 Ekaon*t3*t4 -
3209 Ekaon*t3*t4 +
3212 Ekaon*t3*t4 -
3214 Ekaon*t3*t4 -
3218 t3*t4 +
3224 t3*t4 -
3227 t3*t4 -
3230 t3*t4 +
3232 t3*t4 +
3234 Elep*t3*t4 -
3241 Fm1*t3*t4 -
3244 Fm1*t3*t4 -
3246 t3*t4 -
3248 t3*t4 +
3250 t3*t4 -
3252 Fm1*t3*t4 +
3254 Fm1*t3*t4 -
3256 Fm1*t3*t4 +
3258 Fm1*t3*t4 -
3260 Fm1*t3*t4 +
3262 Fm1*t3*t4 +
3264 Fm1*t3*t4 -
3266 Fm1*t3*t4 +
3272 t3*t4 +
3274 t3*t4 +
3276 Fm1*t3*t4 -
3278 Fm1*t3*t4 +
3280 Fm1*t3*t4 -
3282 Fm1*t3*t4 -
3289 Fm2*t3*t4 +
3291 t3*t4 -
3293 t3*t4 +
3295 t3*t4 -
3297 t3*t4 -
3299 Fm2*t3*t4 +
3301 Fm2*t3*t4 +
3303 t3*t4 -
3305 t3*t4 -
3307 Fm2*t3*t4 +
3309 Fm2*t3*t4 -
3311 Fm2*t3*t4 +
3313 Fm2*t3*t4 +
3319 t3*t4 -
3321 t3*t4 +
3323 Fm2*t3*t4 -
3325 Fm2*t3*t4 -
3327 Fm2*t3*t4 +
3329 Fm2*t3*t4 -
3348 (t4*t4) +
3352 (t4*t4) -
3356 (t4*t4) -
3363 (t4*t4) -
3365 (t4*t4) -
3369 (t4*t4) -
3371 (t4*t4) -
3373 (t4*t4) -
3375 (t4*t4) +
3380 36.*am*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 +
3381 36.*amSig*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 +
3382 36.*am*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 +
3383 36.*amSig*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 +
3388 36.*am*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 +
3389 36.*amLam*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 -
3390 12.*am*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 -
3391 12.*amLam*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 +
3399 t4 +
3401 t4 +
3403 t4 +
3405 t4 -
3413 3.*C3*
3414 (6.*t1*
3415 (t2 +
3418 (3.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4) -
3421 (3.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4) +
3423 (-3.*(C1*C1)*(C4*C4*C4*C4)*
3426 1.*(C5*C5)*C6*t2*t4 +
3427 C1*(C4*C4)*t3*
3428 (3.*t2 +
3429 (C5*C5)*C6*
3438 (9.*(C1*C1)*(C4*C4)*
3439 (-1. + 2.*C4 +
3440 (C4*C4)*
3441 (-1. + apkk1*C3 +
3444 (C5*C5)*C6*t4*
3445 (3.*C3*t2 +
3446 C6*
3447 (-9. - 6.*C5 - 1.*(C5*C5) +
3448 (C5*C5*C5)*
3451 3.*C1*C4*t3*
3452 (C5*C6*
3453 (6. + 2.*C5 +
3454 C4*
3455 (-6. - 2.*C5 +
3456 (C5*C5)*
3462 C3*C4*
3463 (3.*t2 +
3464 (C5*C5)*C6*
3466 ))) -
3467 6.*t1*
3468 (-1.*apkk1*(2. + C2)*
3472 (3.*C1*C4*
3473 (2.*(-1. + C4) +
3474 C2*
3477 3.*C3*
3478 (t2 -
3479 1.*C1*(C4*C4)*
3481 ) +
3482 C5*C6*
3483 (-6. - 2.*C5 +
3484 (C5*C5)*
3486 C2*
3489 1.*am*
3490 (3.*C1*C4*
3491 (apkk1*(2. + C2)*Fm1 +
3492 2.*Ekaon*
3493 (-2. + 2.*C2*(-1. + C4) + 2.*C4 +
3495 6.*C3*
3497 1.*Enu*t2 +
3498 amSig*apkk1*C1*(C4*C4)*t3 +
3499 apkk1*C1*(C4*C4)*Ekaon*t3) +
3500 C5*C6*
3501 (-2.*
3502 (6. + 2.*C5 - 1.*apkk1*(C5*C5) +
3503 2.*C2*(3. + C5))*Ekaon +
3504 apkk1*(2. + C2)*Fm2 +
3506 *t4) -
3508 (apkk1*
3509 (-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)
3510 + Ekaon*
3511 (3.*C1*C4*
3513 2.*C2*Fm1)*t3 +
3514 C5*C6*
3516 2.*C2*Fm2)*t4))) -
3518 (-1.*apkk1*
3519 (3.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4)*
3520 (3.*C1*C3*(C4*C4)*t3 - 1.*(C5*C5*C5)*C6*t4)
3521 + 2.*Ekaon*
3522 (9.*(C1*C1)*(C4*C4)*
3525 (t3*t3) -
3526 3.*C1*C4*C5*C6*
3529 2.*C4*
3532 (C5*C5)*(C6*C6)*
3535 *(t4*t4))) +
3537 (-18.*(C1*C1)*(C4*C4)*Ekaon*
3538 (-1. + (-1. + apkk1*C3)*(C4*C4) -
3541 2.*(C5*C5)*C6*t4*
3543 t2 +
3544 C6*Ekaon*
3545 (9. + (C5*C5) - 1.*apkk1*(C5*C5*C5) +
3548 3.*C1*C4*t3*
3549 (C5*C6*
3550 (2.*
3551 (2.*(3. + C5) +
3552 C4*(-6. - 2.*C5 + apkk1*(C5*C5)))*
3553 Ekaon +
3554 amLam*
3555 (apkk1*C4*(C5*C5) +
3556 2.*Ekaon*
3558 amSig*
3559 (-1.*apkk1*C4*(C5*C5) +
3560 2.*Ekaon*
3562 C3*C4*
3564 6.*Enu*t2 -
3565 1.*amLam*apkk1*(C5*C5)*C6*t4 +
3566 amSig*apkk1*(C5*C5)*C6*t4 +
3567 2.*apkk1*(C5*C5)*C6*Ekaon*t4))) -
3568 9.*apkk1*C2*C7*
3570 am*apkk1*
3572 (t3*t3) -
3576 t4 +
3578 ((3. + 2.*C5)*Fm1 +
3579 (-1. + 2.*C4)*Fm2)*t3*t4 +
3581 9.*C2*C7*
3589 (t3*t3) +
3594 t4 +
3596 C6*t3*t4 +
3598 t3*t4 +
3602 t4 -
3604 t4 -
3608 t4 -
3615 (9.*(C1*C1)*(C4*C4)*
3618 3.*C1*C4*C5*C6*
3623 (C5*C5)*(C6*C6)*
3627 (-9.*(C1*C1)*(C4*C4)*
3628 (-1. +
3629 (C4*C4)*
3630 (-1. + 3.*apkk1*C3 -
3635 C4*
3638 (t3*t3) +
3639 3.*C1*C4*C5*C6*
3640 (C4*
3641 (-6. +
3642 (C5*C5)*
3643 (3.*apkk1 -
3646 C5*
3647 (-2. + 3.*apkk1*C3 -
3655 2.*
3657 1.*C5*
3664 (C5*C5)*(C6*C6)*
3665 (9. + (C5*C5) +
3666 (C5*C5*C5)*
3667 (-3.*apkk1 +
3671 C5*
3674 (t4*t4)) +
3675 18.*amSig*C1*C2*C4*C7*t3*t5 -
3676 36.*amSig*C1*C2*(C4*C4)*C7*t3*t5 +
3677 54.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
3678 18.*amLam*C2*C5*C6*C7*t4*t5 +
3679 12.*amLam*C2*(C5*C5)*C6*C7*t4*t5 +
3680 18.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
3681 am*
3682 (18.*(C1*C1)*(C4*C4)*
3683 (2.*
3684 (1. - 2.*C4 +
3685 (1. - 1.*apkk1*C3)*(C4*C4))*Ekaon +
3686 amSig*
3687 (1. + (1. - 2.*apkk1*C3)*(C4*C4) +
3689 ) +
3693 3.*C1*C4*t3*
3694 (6.*amSig*C5*C6*t4 -
3695 6.*amSig*C4*C5*C6*t4 +
3696 2.*amSig*(C5*C5)*C6*t4 -
3697 2.*amSig*C4*(C5*C5)*C6*t4 +
3698 amSig*apkk1*C4*(C5*C5*C5)*C6*t4 +
3699 24.*C5*C6*Ekaon*t4 -
3700 24.*C4*C5*C6*Ekaon*t4 +
3701 8.*(C5*C5)*C6*Ekaon*t4 -
3702 8.*C4*(C5*C5)*C6*Ekaon*t4 +
3703 4.*apkk1*C4*(C5*C5*C5)*C6*Ekaon*t4 -
3708 amLam*C5*C6*
3709 (C4*
3710 (-6. - 2.*C5 +
3711 (C5*C5)*
3715 2.*
3719 ))*t4 +
3720 C3*C4*
3722 amLam*apkk1*(C5*C5)*C6*t4 +
3723 3.*amSig*apkk1*(C5*C5)*C6*t4 -
3726 4.*Ekaon*
3727 (3.*t2 + apkk1*(C5*C5)*C6*t4)) +
3728 6.*C2*C7*t5 - 12.*C2*C4*C7*t5) +
3729 2.*C5*C6*t4*
3730 (-3.*C3*C5*
3732 18.*C5*C6*Ekaon*t4 +
3733 12.*(C5*C5)*C6*Ekaon*t4 +
3734 2.*(C5*C5*C5)*C6*Ekaon*t4 -
3735 2.*apkk1*(C5*C5*C5*C5)*C6*Ekaon*t4 +
3739 amLam*C5*C6*
3740 (9. + (C5*C5) - 2.*apkk1*(C5*C5*C5) +
3743 9.*C2*C7*t5 + 6.*C2*C5*C7*t5)))))
3744 + 2.*(akpk*akpk)*
3747 (t3*t3) -
3748 3.*C1*C4*C5*C6*
3751 + 3.*Elep*
3754 amSig*C1*C4*t3*
3756 Fm1*
3757 (-6.*(-1. + C2)*t1 -
3759 6.*C2*(-1. + C4)*C7*t5) +
3760 amLam*C5*C6*t4*
3762 Fm2*
3763 (-2.*(-1. + C2)*t1 +
3765 2.*C2*(3. + C5)*C7*t5)) +
3768 3.*amSig*C1*C4*t3*
3771 C3*C4*
3772 (-6.*t1 - 3.*t2 -
3779 6.*C2*C7*Fm1*t5) +
3780 amLam*C5*C6*t4*
3783 C5*
3786 t4) +
3788 ) + am*
3792 3.*C1*C4*t3*
3795 C3*C4*
3796 (-6.*t1 - 3.*t2 +
3798 6.*C2*C7*Fm1*t5) +
3799 C5*C6*t4*
3800 (6.*(C5*C5)*t1 + 3.*C3*C5*t2 -
3803 6.*C2*C7*Fm2*t5)) +
3804 6.*Elep*
3805 (-1.*C5*C6*t4*
3806 (C2*(3. + C5)*C7*t5 +
3807 Fm2*
3808 ((-1. + C2)*t1 - 2.*C2*C7*Ekaon*t5))
3809 + C1*C4*t3*
3813 t4 - 3.*C2*C7*t5 + 3.*C2*C4*C7*t5 +
3814 Fm1*
3815 (-3.*(-1. + C2)*t1 +
3816 Ekaon*
3817 (-1.*amLam*(-3. + C5)*C5*C6*t4 +
3818 amSig*(-3. + C5)*C5*C6*t4 +
3819 6.*C2*C7*t5)))))) +
3825 72.*am*apkk1*C1*C2*C4*t1*t3 +
3826 72.*amSig*apkk1*C1*C2*C4*t1*t3 -
3829 72.*am*apkk1*C1*(C4*C4)*t1*t3 -
3830 72.*amSig*apkk1*C1*(C4*C4)*t1*t3 -
3878 (t3*t3) +
3880 (t3*t3) +
3882 (t3*t3) +
3895 (t3*t3) -
3897 (t3*t3) +
3899 Ekaon*(t3*t3) +
3901 (t3*t3) -
3904 (t3*t3) +
3906 (t3*t3) -
3908 (t3*t3) +
3910 (t3*t3) +
3912 Enu*(t3*t3) -
3914 (t3*t3) +
3916 (t3*t3) -
3919 (t3*t3) -
3921 (t3*t3) -
3924 (t3*t3) -
3926 (t3*t3) -
3928 (t3*t3) +
3930 (t3*t3) -
3932 (t3*t3) +
3934 (t3*t3) +
3936 (t3*t3) +
3938 (t3*t3) -
3940 (t3*t3) -
3942 (t3*t3) +
3944 (t3*t3) +
3946 (t3*t3) -
3950 (t3*t3) -
3956 (t3*t3) -
3966 72.*am*apkk1*C2*C5*C6*t1*t4 +
3967 72.*amLam*apkk1*C2*C5*C6*t1*t4 +
3970 24.*am*apkk1*(C5*C5)*C6*t1*t4 +
3971 24.*amLam*apkk1*(C5*C5)*C6*t1*t4 +
4016 t3*t4 -
4022 C6*t3*t4 +
4025 t3*t4 +
4028 C6*t3*t4 -
4031 (C5*C5)*C6*t3*t4 -
4033 t3*t4 -
4035 C6*t3*t4 -
4037 (C5*C5)*C6*t3*t4 -
4039 (C5*C5)*C6*t3*t4 -
4043 C6*t3*t4 -
4045 t3*t4 -
4047 C6*t3*t4 -
4049 C6*t3*t4 +
4054 t3*t4 -
4057 t3*t4 +
4060 t3*t4 -
4063 C6*Ekaon*t3*t4 -
4066 Ekaon*t3*t4 +
4069 C6*Ekaon*t3*t4 +
4071 Ekaon*t3*t4 -
4073 (C5*C5)*C6*Ekaon*t3*t4 +
4076 Ekaon*t3*t4 +
4078 Ekaon*t3*t4 -
4080 C6*Ekaon*t3*t4 -
4082 Enu*t3*t4 +
4085 (C5*C5)*C6*Enu*t3*t4 -
4087 C6*Enu*t3*t4 -
4089 (C5*C5)*C6*Enu*t3*t4 -
4091 t3*t4 +
4094 (C5*C5*C5)*C6*Enu*t3*t4 -
4096 Enu*t3*t4 -
4098 C6*Enu*t3*t4 +
4100 Enu*t3*t4 -
4106 t3*t4 -
4108 Enu*t3*t4 -
4114 Enu*t3*t4 +
4118 Fm1*t3*t4 +
4121 C6*Fm1*t3*t4 +
4125 C6*Fm1*t3*t4 +
4127 t3*t4 +
4129 Fm1*t3*t4 -
4132 Fm1*t3*t4 +
4134 t3*t4 -
4136 t3*t4 -
4138 t3*t4 +
4140 Fm1*t3*t4 -
4145 Fm1*t3*t4 +
4147 Fm1*t3*t4 +
4151 t3*t4 -
4153 Fm1*t3*t4 -
4155 Fm1*t3*t4 +
4158 Fm1*t3*t4 -
4160 t3*t4 +
4162 t3*t4 -
4164 t3*t4 +
4166 t3*t4 +
4168 Fm1*t3*t4 -
4170 Fm1*t3*t4 +
4174 Fm2*t3*t4 -
4177 C6*Fm2*t3*t4 -
4181 C6*Fm2*t3*t4 -
4183 t3*t4 -
4185 Fm2*t3*t4 -
4189 Fm2*t3*t4 -
4191 t3*t4 +
4193 t3*t4 +
4195 Fm2*t3*t4 +
4200 Fm2*t3*t4 -
4202 Fm2*t3*t4 -
4206 t3*t4 +
4208 Fm2*t3*t4 +
4210 Fm2*t3*t4 -
4213 Fm2*t3*t4 -
4215 t3*t4 +
4217 t3*t4 +
4219 t3*t4 -
4221 t3*t4 +
4223 Fm2*t3*t4 -
4225 Fm2*t3*t4 -
4227 Fm2*t3*t4 +
4236 Fm2*t3*t4 -
4242 Fm2*t3*t4 -
4257 (t4*t4) +
4271 (t4*t4) +
4273 (t4*t4) +
4275 (t4*t4) -
4278 (t4*t4) +
4280 (t4*t4) -
4283 (t4*t4) +
4285 (t4*t4) -
4287 (t4*t4) +
4289 (t4*t4) -
4292 (t4*t4) +
4294 (t4*t4) +
4297 (t4*t4) +
4299 (t4*t4) -
4302 (t4*t4) -
4304 (t4*t4) +
4306 (t4*t4) -
4309 (t4*t4) +
4311 (t4*t4) +
4313 (t4*t4) +
4315 (t4*t4) -
4317 (t4*t4) -
4319 (t4*t4) +
4321 (t4*t4) -
4327 (t4*t4) -
4335 (amLam*C5*C6*
4337 Fm2*
4338 (2.*(1. + C2)*t1 +
4340 amSig*C1*C4*t3*
4342 Fm1*
4343 (6.*(1. + C2)*t1 -
4345 2.*am*
4346 ((1. + C2)*C5*C6*Fm2*t1*t4 +
4347 C1*C4*t3*
4350 Fm1*
4351 (3.*(1. + C2)*t1 -
4353 C6*Ekaon*t4)))) +
4355 144.*apkk1*(C2*C2)*C7*t1*t5 -
4372 t3*t5 +
4376 t3*t5 -
4393 t4*t5 +
4398 144.*apkk1*(C2*C2)*(C7*C7)*Ekaon*(t5*t5) +
4400 144.*apkk1*(C2*C2)*(C7*C7)*Enu*(t5*t5) -
4401 1.*Elep*(36.*((-1. + C2)*(-1. + C2))*(t1*t1) -
4411 (t3*t3) +
4421 t4 +
4423 t4 -
4425 (C5*C5)*C6*t3*t4 +
4428 (C5*C5*C5)*C6*t3*t4 +
4430 t4 -
4432 t4 +
4434 t4 -
4436 t4 +
4438 t3*t4 -
4440 t3*t4 +
4442 t3*t4 -
4444 t3*t4 -
4446 t4 +
4448 t4 -
4450 t4 +
4452 t4 +
4454 t3*t4 -
4456 t3*t4 +
4458 t3*t4 -
4469 (9.*(C1*C1)*(C4*C4)*
4470 (-1. + 2.*C4 +
4472 3.*C1*C4*C5*
4473 (2.*(3. + C5) +
4474 C4*
4477 (C5*C5)*
4478 (-9. - 6.*C5 - 1.*(C5*C5) +
4480 6.*t1*
4484 3.*C3*
4485 (2.*t2 +
4490 1.*am*
4491 (3.*C1*(2. + C2)*C4*Fm1*t3 +
4492 C5*C6*
4494 t4)) +
4495 2.*
4497 (3.*C1*C4*(1. + C4)*t3 -
4498 1.*(-3. + C5)*C5*C6*t4) -
4499 1.*apkk1*(1. + C2)*
4503 (3.*C1*C4*
4506 1.*C5*C6*
4509 + am*
4510 (-3.*C1*C4*
4511 (2.*(-1. + C2)*(1. + C4)*Ekaon +
4512 apkk1*(1. + C2)*Fm1)*t3 +
4513 3.*amSig*C1*C4*
4516 1.*C5*C6*
4517 (-2.*(-1. + C2)*(-3. + C5)*Ekaon +
4518 apkk1*(1. + C2)*Fm2 +
4519 amLam*
4522 ) -
4525 (3.*(C1*C1)*(C4*C4)*
4526 (3. - 2.*C4 + 3.*(C4*C4))*(t3*t3) +
4527 2.*C1*C4*C5*
4528 (9. + C5 - 3.*C4*(1. + C5))*C6*t3*t4
4529 + (C5*C5)*(9. + 2.*C5 + (C5*C5))*(C6*C6)*
4530 (t4*t4)) -
4531 4.*Ekaon*
4532 (9.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
4533 (t3*t3) -
4534 3.*C1*C4*C5*C6*
4535 (2.*(1. + C4)*(-3. + C5)*Ekaon +
4536 amSig*
4540 amLam*
4544 ((-3. + C5)*(-3. + C5))*(C5*C5)*(C6*C6)*Ekaon*(t4*t4)
4545 ) +
4547 (-9.*(C1*C1)*C3*(C4*C4*C4*C4)*
4549 (C5*C5)*C6*t4*
4550 (6.*C3*t2 +
4551 (C5*C5*C5)*C6*
4553 3.*C1*(C4*C4)*t3*
4554 ((C5*C5*C5)*C6*
4556 2.*Ekaon*
4558 + C3*
4559 (6.*t2 +
4560 (C5*C5)*C6*
4562 2.*Ekaon*
4564 )) -
4567 36.*amSig*apkk1*C1*C2*C4*C7*t3*t5 +
4568 36.*amSig*apkk1*C1*C2*(C4*C4)*C7*t3*t5 +
4573 36.*amLam*apkk1*C2*C5*C6*C7*t4*t5 -
4574 12.*amLam*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
4578 144.*apkk1*(C2*C2)*(C7*C7)*(t5*t5) -
4579 1.*am*
4582 (t3*t3) -
4584 (t3*t3) -
4585 3.*amSig*C1*C4*t3*
4586 (C5*C6*
4588 2.*Ekaon*
4591 C3*C4*
4592 (6.*t2 -
4595 (12.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4))
4596 ) +
4597 3.*C1*C4*t3*
4604 C3*C4*
4606 4.*Ekaon*
4608 6.*C2*C7*t5 +
4610 2.*C5*C6*t4*
4611 (-3.*amLam*C3*C5*t2 -
4612 6.*C3*C5*Ekaon*t2 +
4617 )) +
4618 2.*
4621 (t3*t3) +
4623 (t3*t3) +
4625 t4 -
4627 t3*t4 +
4629 t3*t4 -
4631 Ekaon*t3*t4 +
4633 Fm1*t3*t4 -
4635 Fm1*t3*t4 +
4637 Fm1*t3*t4 -
4639 Fm1*t3*t4 -
4641 t3*t4 +
4643 t3*t4 +
4645 Fm2*t3*t4 -
4647 Fm2*t3*t4 +
4649 (t4*t4) +
4651 (t4*t4) +
4653 18.*apkk1*C1*C2*C4*C7*t3*t5 -
4654 18.*apkk1*C1*C2*(C4*C4)*C7*t3*t5 +
4656 t5 - 18.*apkk1*C2*C5*C6*C7*t4*t5 +
4657 6.*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
4659 t5 +
4661 (18.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
4662 (t3*t3) +
4663 3.*amSig*C1*C4*t3*
4664 (6.*C1*C4*(1. + (C4*C4))*t3 +
4665 C5*C6*
4669 2.*C5*C6*t4*
4670 (amLam*C5*(9. + (C5*C5))*C6*t4 +
4671 9.*C5*C6*Ekaon*t4 -
4672 6.*(C5*C5)*C6*Ekaon*t4 +
4673 (C5*C5*C5)*C6*Ekaon*t4 + 9.*C2*C7*t5) -
4674 3.*C1*C4*t3*
4675 (amLam*C5*C6*
4679 2.*
4680 (-6.*(1. + C4)*C5*C6*Ekaon*t4 +
4681 2.*(1. + C4)*(C5*C5)*C6*Ekaon*t4 -
4682 3.*C2*C7*t5))))))))));
4683
4684 if (sol <= 0.0)
4685 std::cout << "Check Matrix and def. of Scalar Products" << std::endl;
4686
4687 return sol;
4688}
References alepvec, am, amk, aml, amLam, ampi, amSig, angkq, aq0, aqvec, d, Ekaon, Elep, Enu, f, Fm1, Fm2, fpi, g, pkvec, and Vus.
Referenced by diffxsec().
◆ Amatrix_PP()
| double singlekaon_xsec::Amatrix_PP | ( | double | theta, |
| double | phikq ) |
Definition at line 4691 of file singlekaon_xsec.cxx.
4691 {
4692 double sol = 0.;
4693
4694 double akk1=0., zdotq=0., qdotpk=0., akcrosk1=0., qcrospk=0.;
4695 double zdotpk=0., azpk=0., aqkaon=0., akpk=0., apkk1=0.;
4696 double C1=0., C2=0., C3=0., C4=0., C5=0., C6=0., C7=0., C8=0., C9=0.;
4697 double aq2=0., gform=0., con=0., t1=0., t2=0., t3=0., t4=0., t5=0., t6=0.;
4698
4707 akpk=azpk + aqkaon/2.0;
4708 apkk1=azpk - aqkaon/2.0;
4719 gform=1.0/pow(1.0-aq2/(1.1*1.1),4);
4720
4722
4723 t1=1.0;
4724 t2=1.0; // !Full Term
4725 t3=1.0;
4726 t4=1.0;
4727 t5=1.0;
4728 t6=1.0;
4729
4730 sol = con*gform*(0.4444444444444444*am*
4748 144.*akpk*amSig*C1*C4*t1*t3 +
4749 144.*amSig*apkk1*C1*C4*t1*t3 -
4750 144.*akpk*amSig*C1*(C4*C4)*t1*t3 +
4751 144.*amSig*apkk1*C1*(C4*C4)*t1*t3 +
4787 (C4*C4)*t2*t3 -
4797 Elep*t2*t3 +
4810 (C4*C4*C4*C4)*(t3*t3) -
4813 (C4*C4*C4*C4)*(t3*t3) +
4816 (C4*C4*C4*C4)*(t3*t3) +
4818 (C4*C4*C4*C4)*(t3*t3) -
4821 (C3*C3)*(C4*C4*C4*C4)*(t3*t3) +
4831 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
4834 Ekaon*(t3*t3) +
4837 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
4839 Ekaon*(t3*t3) -
4841 Ekaon*(t3*t3) -
4843 Ekaon*(t3*t3) -
4853 (C3*C3)*(C4*C4*C4*C4)*Elep*(t3*t3) +
4856 (C4*C4*C4*C4)*Elep*(t3*t3) +
4858 (C4*C4*C4*C4)*Elep*(t3*t3) +
4868 (C4*C4*C4*C4)*Enu*(t3*t3) +
4871 (C4*C4*C4*C4)*Enu*(t3*t3) -
4874 (C4*C4*C4*C4)*Enu*(t3*t3) -
4876 (C4*C4*C4*C4)*Enu*(t3*t3) +
4883 Fm1*(t3*t3) +
4889 Fm1*(t3*t3) -
4925 144.*akpk*amLam*C5*C6*t1*t4 -
4926 144.*amLam*apkk1*C5*C6*t1*t4 -
4927 48.*akpk*amLam*(C5*C5)*C6*t1*t4 +
4928 48.*amLam*apkk1*(C5*C5)*C6*t1*t4 +
4965 (C5*C5)*C6*t2*t4 +
4974 Elep*t2*t4 +
4997 C6*t3*t4 -
5000 (C5*C5)*C6*t3*t4 -
5003 (C4*C4)*(C5*C5)*C6*t3*t4 +
5005 C6*t3*t4 +
5008 (C3*C3)*(C4*C4)*(C5*C5)*C6*t3*t4 +
5010 C6*t3*t4 +
5013 (C3*C3)*(C4*C4)*(C5*C5)*C6*t3*t4 -
5015 (C5*C5)*C6*t3*t4 +
5017 (C4*C4)*(C5*C5)*C6*t3*t4 -
5019 (C5*C5)*C6*t3*t4 -
5021 (C5*C5)*C6*t3*t4 -
5023 (C5*C5)*C6*t3*t4 -
5029 Ekaon*t3*t4 -
5032 Ekaon*t3*t4 -
5035 Ekaon*t3*t4 -
5037 Ekaon*t3*t4 -
5039 C6*Ekaon*t3*t4 +
5042 (C4*C4)*(C5*C5)*C6*Ekaon*t3*t4 +
5044 C6*Ekaon*t3*t4 +
5046 (C5*C5)*C6*Ekaon*t3*t4 +
5048 (C5*C5)*C6*Ekaon*t3*t4 -
5050 (C5*C5)*C6*Ekaon*t3*t4 -
5052 (C5*C5)*C6*Ekaon*t3*t4 +
5062 (C5*C5)*C6*Elep*t3*t4 -
5065 (C3*C3)*(C4*C4)*(C5*C5)*C6*Elep*t3*t4 +
5067 (C5*C5)*C6*Elep*t3*t4 -
5069 Elep*t3*t4 +
5071 (C4*C4)*(C5*C5)*C6*Elep*t3*t4 -
5082 (C4*C4)*(C5*C5)*C6*Enu*t3*t4 -
5084 Enu*t3*t4 +
5087 (C3*C3)*(C4*C4)*(C5*C5)*C6*Enu*t3*t4 -
5089 (C5*C5)*C6*Enu*t3*t4 +
5091 Enu*t3*t4 -
5093 (C4*C4)*(C5*C5)*C6*Enu*t3*t4 -
5098 Fm1*t3*t4 +
5105 C6*Fm1*t3*t4 -
5108 Fm1*t3*t4 -
5111 Fm1*t3*t4 +
5114 Fm1*t3*t4 +
5123 Fm1*t3*t4 +
5126 Fm1*t3*t4 -
5131 Fm2*t3*t4 +
5138 C6*Fm2*t3*t4 -
5141 Fm2*t3*t4 -
5144 Fm2*t3*t4 +
5147 Fm2*t3*t4 -
5156 Fm2*t3*t4 +
5159 Fm2*t3*t4 +
5165 Fm2*t3*t4 -
5168 Fm2*t3*t4 +
5170 Fm2*t3*t4 +
5172 Fm2*t3*t4 +
5175 Fm2*t3*t4 -
5186 Fm2*t3*t4 +
5208 (C5*C5*C5*C5)*(C6*C6)*(t4*t4) +
5211 (C6*C6)*(t4*t4) +
5213 (C6*C6)*(t4*t4) -
5216 (C5*C5*C5*C5)*(C6*C6)*(t4*t4) +
5226 (C6*C6)*Ekaon*(t4*t4) +
5229 Ekaon*(t4*t4) +
5232 (C6*C6)*Ekaon*(t4*t4) +
5234 Ekaon*(t4*t4) -
5237 (C6*C6)*Ekaon*(t4*t4) -
5247 (C5*C5*C5*C5)*(C6*C6)*Elep*(t4*t4) +
5250 (C6*C6)*Elep*(t4*t4) +
5252 (C6*C6)*Elep*(t4*t4) +
5262 Enu*(t4*t4) +
5265 Enu*(t4*t4) -
5268 (C6*C6)*Enu*(t4*t4) -
5270 (C6*C6)*Enu*(t4*t4) +
5277 Fm2*(t4*t4) -
5287 Fm2*(t4*t4) -
5318 (1. +
5319 (C4*C4)*
5326 6.*C1*C4*C5*C6*
5331 C4*
5332 (C5*
5337 (C5*C5)*(C6*C6)*
5338 (9. +
5339 (C5*C5)*
5347 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
5349 (9.*(C1*C1)*(C4*C4)*
5350 (2.*C3*(C4*C4)*Ekaon +
5354 (t3*t3) +
5355 6.*C1*C4*C5*C6*
5356 (2.*C3*C4*C5*Ekaon +
5360 t4 +
5361 (C5*C5)*(C6*C6)*
5362 (2.*C3*(C5*C5)*Ekaon +
5366 (t4*t4)) -
5368 (apkk1*
5371 akpk*
5374 2.*Ekaon*
5375 (3.*C1*C4*
5378 C5*C6*
5394 12.*apkk1*C4*Fm1 +
5399 4.*amSig*
5400 (1. +
5401 (C4*C4)*
5408 - 1.*Ekaon*
5409 (4. +
5410 (C4*C4)*
5419 2.*Ekaon*
5420 (2.*(akpk*akpk)*
5422 2.*(apkk1*apkk1)*
5424 apkk1*
5425 (2. + 4.*C4 +
5426 (C4*C4)*
5432 1.*akpk*
5433 (2. - 4.*C4 +
5434 (C4*C4)*
5437 (-1. + apkk1*C3 +
5441 2.*Ekaon*
5444 2.*
5445 (((1.+C4)*(1.+C4))*Elep -
5446 1.*((-1.+C4)*(-1.+C4))*Enu) +
5447 amSig*
5448 (-2. +
5457 ))*(t3*t3) +
5458 6.*C1*C4*t3*
5491 (amLam*
5492 (6. +
5496 amSig*
5497 (6. +
5501 2.*
5502 (-1.*Enu*
5505 Elep*
5508 Ekaon*
5511 t1 + 12.*akpk*C5*C6*t4 +
5513 12.*apkk1*C5*C6*t4 +
5514 4.*akpk*(C5*C5)*C6*t4 +
5515 4.*apkk1*(C5*C5)*C6*t4 -
5526 4.*Elep*
5528 Fm1*
5529 (6.*t1 -
5538 t4 +
5540 t4 +
5547 Elep*t4 -
5555 t4 -
5557 Enu*t4 +
5571 6.*akpk*
5578 t4 -
5579 2.*C5*C6*
5581 2.*C5*Enu +
5584 t4 +
5586 (-1.*(C5*C5)*C6*
5589 (3.*t2 +
5590 (-1.*akpk + apkk1)*(C5*C5)*C6*t4)))
5591 + Ekaon*
5596 12.*apkk1*C5*C6*t4 +
5598 4.*apkk1*(C5*C5)*C6*t4 +
5601 t4 -
5605 t4 -
5614 2.*akpk*
5616 2.*C5*(3. + C5)*C6*t4 +
5618 (3.*C3*t2 - 2.*(C5*C5)*C6*t4 +
5619 2.*apkk1*C3*(C5*C5)*C6*t4))))) +
5620 C5*C6*t4*
5622 Fm2 +
5623 Ekaon*
5626 t1 +
5627 (C5*C5*C5)*C6*
5632 Ekaon*
5636 4.*amLam*
5644 2.*Ekaon*
5647 apkk1*
5651 2.*Ekaon*
5652 (amLam*
5655 2.*Enu +
5659 1.*akpk*
5662 (-1. + apkk1*C3 +
5664 4.*(C5*C5)*C6*
5671 )*t4 +
5686 12.*akpk*
5690 (2.*C6*
5694 amLam*
5699 t4)) -
5700 1.*Ekaon*
5714 6.*akpk*
5715 (6.*C6*t4 +
5717 ))))) -
5867 (t3*t3) -
5870 (t3*t3) +
5872 (t3*t3) +
5874 (t3*t3) -
5877 (t3*t3) -
5879 (t3*t3) +
5881 (t3*t3) -
5883 (t3*t3) +
5885 (t3*t3) -
5887 (t3*t3) -
5889 (t3*t3) +
5892 (t3*t3) +
5894 (t3*t3) -
5896 (t3*t3) -
5899 (t3*t3) -
5901 (t3*t3) -
5903 (t3*t3) -
5906 (t3*t3) +
5917 (t3*t3) +
5919 (t3*t3) +
5939 t3*t4 -
5948 t3*t4 +
5954 t3*t4 +
5968 Fm1*t3*t4 +
5971 Fm1*t3*t4 +
5974 t3*t4 +
5977 t3*t4 -
5979 t3*t4 -
5988 t3*t4 -
5991 Fm1*t3*t4 -
5994 Fm1*t3*t4 -
5996 t3*t4 -
5998 t3*t4 -
6001 Fm1*t3*t4 -
6003 t3*t4 -
6005 t3*t4 +
6007 t3*t4 -
6009 t3*t4 -
6011 t3*t4 -
6013 t3*t4 +
6016 t3*t4 +
6019 Fm1*t3*t4 -
6021 t3*t4 -
6023 t3*t4 -
6025 t3*t4 +
6027 t3*t4 -
6030 Fm2*t3*t4 -
6033 Fm2*t3*t4 +
6037 t3*t4 -
6039 t3*t4 +
6048 t3*t4 -
6051 Fm2*t3*t4 +
6054 Fm2*t3*t4 +
6056 t3*t4 +
6058 t3*t4 -
6061 Fm2*t3*t4 -
6064 Fm2*t3*t4 +
6066 t3*t4 +
6068 t3*t4 -
6070 t3*t4 +
6072 t3*t4 +
6075 t3*t4 +
6078 Fm2*t3*t4 +
6080 t3*t4 +
6082 t3*t4 -
6084 t3*t4 +
6086 t3*t4 +
6089 t3*t4 -
6091 Fm2*t3*t4 -
6093 Fm2*t3*t4 -
6095 Fm2*t3*t4 -
6102 t3*t4 -
6104 t3*t4 -
6106 t3*t4 -
6111 Fm2*t3*t4 +
6140 (t4*t4) +
6142 (t4*t4) +
6144 (t4*t4) +
6147 (t4*t4) +
6150 (t4*t4) -
6152 (t4*t4) +
6154 (t4*t4) +
6156 (t4*t4) +
6158 (t4*t4) +
6161 (t4*t4) +
6163 (t4*t4) +
6165 (t4*t4) -
6168 (t4*t4) -
6170 (t4*t4) +
6172 (t4*t4) -
6176 (t4*t4) +
6178 (t4*t4) +
6183 (t4*t4) +
6185 (t4*t4) +
6187 (t4*t4) +
6189 (t4*t4) +
6192 (t4*t4) +
6237 2.*akpk*
6238 (3.*C1*C4*
6243 t3 +
6244 C5*C6*
6249 - 2.*apkk1*
6250 (3.*C1*C4*
6254 t3 +
6255 C5*C6*
6260 (18.*apkk1*C1*(C4*C4)*t2*t3 +
6265 6.*apkk1*(C5*C5)*C6*t2*t4 +
6270 t3*t4 -
6273 akpk*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
6274 (-6.*t2 +
6277 (9.*(C1*C1)*(C4*C4*C4*C4)*
6281 (t3*t3) +
6282 (C5*C5)*C6*t4*
6283 (-12.*t2 -
6284 1.*(C5*C5)*C6*
6288 t4) -
6289 3.*C1*(C4*C4)*t3*
6290 (12.*t2 +
6291 (C5*C5)*C6*
6292 (2.*apkk1 +
6295 amLam*
6297 )*t4))) -
6298 18.*akpk*C1*C2*C4*C7*Ekaon*t3*t5 +
6300 18.*apkk1*C1*C2*C4*C7*Ekaon*t3*t5 +
6301 18.*akpk*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 +
6302 18.*apkk1*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 -
6307 72.*akpk*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
6314 18.*akpk*C2*C5*C6*C7*Ekaon*t4*t5 -
6316 18.*apkk1*C2*C5*C6*C7*Ekaon*t4*t5 +
6317 6.*akpk*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
6318 6.*apkk1*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
6323 24.*akpk*apkk1*C2*C5*C6*C7*Fm2*t4*t5 -
6330 18.*akpk*C1*C4*C8*C9*Ekaon*t3*t6 -
6332 18.*apkk1*C1*C4*C8*C9*Ekaon*t3*t6 -
6333 18.*akpk*C1*(C4*C4)*C8*C9*Ekaon*t3*t6 -
6334 18.*apkk1*C1*(C4*C4)*C8*C9*Ekaon*t3*t6 +
6339 72.*akpk*apkk1*C1*C4*C8*C9*Fm1*t3*t6 -
6346 18.*akpk*C5*C6*C8*C9*Ekaon*t4*t6 +
6348 18.*apkk1*C5*C6*C8*C9*Ekaon*t4*t6 -
6349 6.*akpk*(C5*C5)*C6*C8*C9*Ekaon*t4*t6 -
6350 6.*apkk1*(C5*C5)*C6*C8*C9*Ekaon*t4*t6 -
6355 24.*akpk*apkk1*C5*C6*C8*C9*Fm2*t4*t6 +
6363 (t3*t3) +
6366 (t3*t3) +
6370 (t3*t3) -
6372 (t3*t3) +
6374 (t3*t3) -
6376 (t3*t3) -
6378 (t3*t3) +
6380 (t3*t3) +
6382 (t3*t3) -
6384 (t3*t3) -
6386 (t3*t3) +
6398 t4 -
6401 t3*t4 +
6403 t4 -
6405 t3*t4 -
6407 t3*t4 +
6409 t4 +
6411 Fm2*t3*t4 -
6413 t3*t4 -
6415 t3*t4 -
6417 t3*t4 +
6419 t4 +
6421 t4 +
6423 t3*t4 -
6430 (t4*t4) -
6433 (t4*t4) -
6435 (t4*t4) -
6437 (t4*t4) +
6439 (t4*t4) +
6441 (t4*t4) -
6443 (t4*t4) -
6445 (t4*t4) -
6446 2.*C3*
6448 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
6449 2.*akpk*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
6450 (-3.*t2 +
6452 ) +
6453 Ekaon*
6455 (t3*t3) +
6456 (C5*C5)*C6*t4*
6459 6.*C1*(C4*C4)*t3*
6463 (9.*(C1*C1)*(C4*C4*C4*C4)*
6467 (C5*C5)*C6*t4*
6468 (-6.*t2 +
6469 (C5*C5)*C6*
6473 3.*C1*(C4*C4)*t3*
6474 (6.*t2 +
6475 (C5*C5)*C6*
6479 *t4))) +
6480 (C3*C3)*
6482 36.*(akpk*akpk)*C1*(C4*C4)*t2*t3 +
6484 72.*akpk*apkk1*C1*(C4*C4)*t2*t3 -
6486 36.*(apkk1*apkk1)*C1*(C4*C4)*t2*t3 -
6492 (t3*t3) -
6495 (t3*t3) -
6497 (t3*t3) +
6500 (t3*t3) + 12.*(akpk*akpk)*(C5*C5)*C6*t2*t4 +
6502 24.*akpk*apkk1*(C5*C5)*C6*t2*t4 -
6504 12.*(apkk1*apkk1)*(C5*C5)*C6*t2*t4 -
6509 t4 +
6511 t3*t4 -
6513 t3*t4 +
6515 Enu*t3*t4 +
6517 Enu*t3*t4 -
6519 C6*Enu*t3*t4 -
6521 C6*Enu*t3*t4 +
6523 t4 -
6528 (t4*t4) -
6531 (t4*t4) -
6533 (t4*t4) +
6536 (t4*t4) +
6538 (36.*(t2*t2) -
6540 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6543 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6544 - 1.*
6545 ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))) -
6546 2.*Elep*
6547 (2.*(akpk - 1.*apkk1)*
6548 (3.*amSig*C1*(C4*C4)*t3 +
6549 amLam*(C5*C5)*C6*t4)*
6550 (-3.*t2 +
6552 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
6553 Enu*
6554 (36.*(t2*t2) -
6556 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6559 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6560 - 1.*
6563 2.*Ekaon*
6564 (6.*(akpk - 1.*apkk1)*t2*
6565 (3.*amSig*C1*(C4*C4)*t3 +
6566 amLam*(C5*C5)*C6*t4) +
6567 Enu*
6568 (-36.*(t2*t2) +
6570 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) -
6571 1.*(akpk - 1.*apkk1)*
6573 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6574 + ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))) +
6575 Elep*
6576 (36.*(t2*t2) -
6578 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6579 (akpk - 1.*apkk1)*
6581 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6582 + ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))))) -
6583 9.*C1*C2*C4*C7*Ekaon*t3*t5 +
6584 9.*akpk*C1*C2*C4*C7*Fm1*t3*t5 +
6586 9.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 -
6590 9.*C2*C5*C6*C7*Ekaon*t4*t5 -
6591 3.*akpk*C2*C5*C6*C7*Fm2*t4*t5 -
6593 3.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
6597 9.*C1*C4*C8*C9*Ekaon*t3*t6 -
6598 9.*akpk*C1*C4*C8*C9*Fm1*t3*t6 -
6600 9.*apkk1*C1*C4*C8*C9*Fm1*t3*t6 +
6604 9.*C5*C6*C8*C9*Ekaon*t4*t6 +
6605 3.*akpk*C5*C6*C8*C9*Fm2*t4*t6 +
6607 3.*apkk1*C5*C6*C8*C9*Fm2*t4*t6 -
6614 36.*akpk*amSig*C1*(C4*C4)*t2*t3 -
6616 36.*am*apkk1*C1*(C4*C4)*t2*t3 +
6617 36.*amSig*apkk1*C1*(C4*C4)*t2*t3 -
6632 12.*akpk*am*(C5*C5)*C6*t2*t4 -
6634 12.*akpk*amLam*(C5*C5)*C6*t2*t4 -
6636 12.*am*apkk1*(C5*C5)*C6*t2*t4 +
6637 12.*amLam*apkk1*(C5*C5)*C6*t2*t4 -
6644 C6*t3*t4 +
6648 C6*t3*t4 +
6651 t3*t4 -
6655 t3*t4 -
6660 Enu*t3*t4 -
6662 t3*t4 +
6677 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
6679 (9.*(C1*C1)*(C4*C4*C4*C4)*
6682 6.*C1*(C4*C4)*t3*
6683 (-12.*t2 +
6684 (C5*C5)*C6*
6687 )*t4) +
6688 (C5*C5)*C6*t4*
6689 (-24.*t2 +
6690 (C5*C5)*C6*
6693 + Ekaon*
6694 (108.*(t2*t2) -
6706 t4 -
6708 t4 +
6710 t4 +
6713 t4 +
6715 t4 -
6726 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) -
6727 2.*akpk*
6729 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
6730 1.*
6731 ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))) +
6732 Elep*(-36.*(t2*t2) +
6734 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6737 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
6739 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
6740 2.*am*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
6741 (3.*amSig*C1*(C4*C4)*t3 +
6742 amLam*(C5*C5)*C6*t4) +
6743 ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4))))) +
6745 (3.*C1*C4*
6747 + C5*C6*
6749 t4) +
6755 amSig*
6756 (-4.*Ekaon +
6757 2.*(akpk - 1.*apkk1)*Fm1 +
6759 + am*
6761 2.*(akpk - 1.*apkk1)*Fm1 +
6762 4.*Ekaon*
6764 2.*amSig*
6767 6.*C1*C4*t3*
6783 (-1.*Fm2 +
6784 Fm1*
6792 am*
6799 amSig*C5*C6*
6800 (-1.*Fm2 +
6801 Fm1*
6804 amLam*C5*C6*
6805 (-1.*Fm2 +
6806 Fm1*
6810 + C5*C6*Fm2*t4*
6823 am*
6825 2.*amLam*C5*C6*
6828 2.*
6829 (C5*C6*
6830 ((akpk - 1.*apkk1)*Fm2 +
6831 2.*Ekaon*
6833 3.*C2*C7*t5 + 3.*C8*C9*t6)))) -
6834 2.*Ekaon*
6837 (-12.*f*t1 +
6838 (akpk - 1.*apkk1)*
6842 (akpk*
6844 apkk1*
6846 2.*Ekaon*
6850 t4)) -
6851 6.*am*
6853 C5*C6*Fm2*t4*
6854 (-2.*f*t1 +
6855 Ekaon*
6856 (2.*amLam*C5*C6*t4 + C2*C7*t5 -
6857 1.*C8*C9*t6)) +
6858 C1*C4*t3*
6860 1.*Ekaon*
6863 3.*Fm1*(C2*C7*t5 - 1.*C8*C9*t6)))))
6864 )) + 2.*(Enu*
6912 (t3*t3) -
6916 (t3*t3) -
6920 (t3*t3) -
6927 (t3*t3) -
6929 (t3*t3) -
6933 (t3*t3) -
6935 (t3*t3) -
6937 (t3*t3) -
6939 (t3*t3) +
6978 t3*t4 -
6982 C6*t3*t4 +
6984 t3*t4 +
6989 t3*t4 -
6992 C6*t3*t4 -
6995 t3*t4 +
6997 t3*t4 +
6999 t3*t4 +
7005 Ekaon*t3*t4 -
7008 Ekaon*t3*t4 -
7010 t3*t4 -
7014 Ekaon*t3*t4 +
7017 Ekaon*t3*t4 -
7020 Ekaon*t3*t4 +
7022 Ekaon*t3*t4 +
7026 t3*t4 -
7032 t3*t4 +
7035 t3*t4 +
7038 t3*t4 -
7040 t3*t4 -
7042 Elep*t3*t4 +
7049 Fm1*t3*t4 +
7052 Fm1*t3*t4 +
7054 t3*t4 +
7056 t3*t4 -
7058 t3*t4 +
7060 Fm1*t3*t4 -
7062 Fm1*t3*t4 +
7064 Fm1*t3*t4 -
7066 Fm1*t3*t4 +
7068 Fm1*t3*t4 -
7070 Fm1*t3*t4 -
7072 Fm1*t3*t4 +
7074 Fm1*t3*t4 -
7080 t3*t4 -
7082 t3*t4 -
7084 Fm1*t3*t4 +
7086 Fm1*t3*t4 -
7088 Fm1*t3*t4 +
7090 Fm1*t3*t4 +
7097 Fm2*t3*t4 -
7099 t3*t4 +
7101 t3*t4 -
7103 t3*t4 +
7105 t3*t4 +
7107 Fm2*t3*t4 -
7109 Fm2*t3*t4 -
7111 t3*t4 +
7113 t3*t4 +
7115 Fm2*t3*t4 -
7117 Fm2*t3*t4 +
7119 Fm2*t3*t4 -
7121 Fm2*t3*t4 -
7127 t3*t4 +
7129 t3*t4 -
7131 Fm2*t3*t4 +
7133 Fm2*t3*t4 +
7135 Fm2*t3*t4 -
7137 Fm2*t3*t4 +
7156 (t4*t4) +
7160 (t4*t4) -
7164 (t4*t4) -
7171 (t4*t4) -
7173 (t4*t4) -
7177 (t4*t4) -
7179 (t4*t4) -
7181 (t4*t4) -
7183 (t4*t4) +
7188 18.*am*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 -
7189 18.*amSig*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 -
7190 18.*am*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 -
7191 18.*amSig*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 -
7196 18.*am*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 +
7197 18.*amLam*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 -
7198 6.*am*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 -
7199 6.*amLam*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 +
7204 18.*am*(apkk1*apkk1)*C1*C4*C8*C9*t3*t6 +
7205 18.*amSig*(apkk1*apkk1)*C1*C4*C8*C9*t3*t6 +
7206 18.*am*(apkk1*apkk1)*C1*(C4*C4)*C8*C9*t3*t6 +
7207 18.*amSig*(apkk1*apkk1)*C1*(C4*C4)*C8*C9*t3*t6 +
7212 18.*am*(apkk1*apkk1)*C5*C6*C8*C9*t4*t6 -
7213 18.*amLam*(apkk1*apkk1)*C5*C6*C8*C9*t4*t6 +
7214 6.*am*(apkk1*apkk1)*(C5*C5)*C6*C8*C9*t4*t6 +
7215 6.*amLam*(apkk1*apkk1)*(C5*C5)*C6*C8*C9*t4*t6 -
7220 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
7222 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
7223 (3.*t2 +
7225 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
7226 12.*t1*
7227 (6.*t2 -
7229 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
7231 (9.*(C1*C1)*(C4*C4*C4*C4)*
7234 (C5*C5)*C6*t4*
7235 (6.*t2 +
7236 (C5*C5)*C6*
7239 6.*C1*(C4*C4)*t3*
7240 (3.*t2 +
7241 (C5*C5)*C6*
7244 ) -
7247 (C2*C7*t5 - 1.*C8*C9*t6)) +
7253 (t3*t3) -
7255 (t3*t3) -
7260 (t3*t3) -
7262 (t3*t3) +
7264 (t3*t3) +
7267 (t3*t3) -
7270 (t3*t3) -
7272 (t3*t3) +
7277 t3*t4 -
7279 C6*Ekaon*t3*t4 +
7284 t4 +
7288 Elep*t3*t4 +
7290 Enu*t3*t4 +
7292 t3*t4 -
7294 t3*t4 +
7296 Fm1*t3*t4 -
7298 Fm1*t3*t4 -
7301 t3*t4 -
7303 t4 -
7305 Fm1*t3*t4 -
7307 Fm1*t3*t4 -
7309 t4 +
7311 t4 +
7313 Fm2*t3*t4 -
7315 Fm2*t3*t4 -
7318 t3*t4 +
7320 t4 +
7322 Fm2*t3*t4 +
7324 Fm2*t3*t4 +
7326 t3*t4 +
7328 (t4*t4) -
7330 (t4*t4) -
7335 (t4*t4) -
7337 (t4*t4) +
7339 (t4*t4) +
7342 (t4*t4) +
7345 (t4*t4) -
7347 (t4*t4) +
7349 (9.*(C1*C1)*(C4*C4)*
7350 (-1. + 2.*C4 +
7351 (C4*C4)*
7352 (-1. + apkk1*C3 +
7355 3.*C1*C4*t3*
7356 (C5*C6*
7357 (6. + 2.*C5 - 2.*C4*(3. + C5) -
7362 2.*C3*C4*
7363 (3.*t2 +
7364 (C5*C5)*C6*
7366 )) +
7367 (C5*C5)*C6*t4*
7368 (-1.*C6*
7369 (9. + 6.*C5 + (C5*C5) +
7371 C3*
7372 (6.*t2 +
7373 (C5*C5)*C6*
7375 ))) -
7376 12.*t1*
7377 (apkk1*(2. + f)*
7381 (apkk1*C3*
7382 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
7383 Ekaon*
7384 (-3.*C1*C4*
7387 C5*C6*
7391 (-3.*C1*C4*
7392 (-2. + 2.*C4*(1. + f) +
7394 - 1.*C5*C6*
7395 (6. + 2.*C5*(1. + f) -
7397 *t4 +
7398 C3*
7399 (6.*t2 +
7401 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))) -
7402 1.*am*
7403 (-3.*C1*C4*
7405 2.*Ekaon*
7408 1.*C5*C6*
7410 2.*Ekaon*
7413 2.*C3*
7415 3.*amSig*apkk1*C1*(C4*C4)*t3 +
7416 amLam*apkk1*(C5*C5)*C6*t4 +
7417 Ekaon*
7418 (6.*t2 + 3.*apkk1*C1*(C4*C4)*t3 +
7419 apkk1*(C5*C5)*C6*t4)))) -
7420 9.*am*apkk1*C1*C2*C4*C7*t3*t5 -
7421 9.*amSig*apkk1*C1*C2*C4*C7*t3*t5 +
7423 9.*am*apkk1*C2*C5*C6*C7*t4*t5 +
7424 9.*amLam*apkk1*C2*C5*C6*C7*t4*t5 -
7426 9.*am*apkk1*C1*C4*C8*C9*t3*t6 +
7427 9.*amSig*apkk1*C1*C4*C8*C9*t3*t6 -
7429 9.*am*apkk1*C5*C6*C8*C9*t4*t6 -
7430 9.*amLam*apkk1*C5*C6*C8*C9*t4*t6 +
7434 (-1. + 2.*C4 + (-1. + apkk1*C3)*(C4*C4))*
7435 (t3*t3) -
7438 t4 +
7445 (t4*t4) +
7447 (9.*(C1*C1)*(C4*C4)*
7450 6.*C1*C4*C5*C6*
7453 t4 +
7454 (C5*C5)*(C6*C6)*
7458 (-9.*(C1*C1)*(C4*C4)*
7459 (-1. +
7460 (C4*C4)*
7461 (-1. + 3.*apkk1*C3 -
7466 C4*
7469 (t3*t3) -
7470 6.*C1*C4*C5*C6*
7477 C4*
7478 (-3. +
7479 C5*
7480 (-1. + 3.*apkk1*C3 -
7488 (C5*C5)*(C6*C6)*
7489 (9. +
7490 (C5*C5)*
7491 (1. - 3.*apkk1*C3 +
7496 C5*
7499 (t4*t4)) -
7500 27.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
7501 9.*amLam*C2*C5*C6*C7*t4*t5 +
7502 6.*amLam*C2*(C5*C5)*C6*C7*t4*t5 +
7503 9.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
7504 27.*apkk1*C1*C4*C8*C9*Fm1*t3*t6 -
7505 9.*amLam*C5*C6*C8*C9*t4*t6 -
7506 6.*amLam*(C5*C5)*C6*C8*C9*t4*t6 -
7507 9.*apkk1*C5*C6*C8*C9*Fm2*t4*t6 -
7508 3.*amSig*C1*C4*t3*
7509 (2.*amLam*C5*
7510 (3. + C5 +
7511 C4*(-3. - 1.*C5 + apkk1*C3*C5))*C6*
7512 t4 +
7513 apkk1*C5*C6*
7515 2.*C4*Fm2)*t4 -
7516 3.*(-1. + 2.*C4)*
7517 (C2*C7*t5 - 1.*C8*C9*t6)) +
7518 am*
7519 (18.*(C1*C1)*(C4*C4)*
7520 (2.*
7521 (1. - 2.*C4 +
7522 (1. - 1.*apkk1*C3)*(C4*C4))*Ekaon +
7523 amSig*
7524 (1. + (1. - 2.*apkk1*C3)*(C4*C4) +
7526 ) +
7530 3.*C1*C4*t3*
7531 (6.*amSig*C5*C6*t4 -
7532 6.*amSig*C4*C5*C6*t4 +
7533 2.*amSig*(C5*C5)*C6*t4 -
7534 2.*amSig*C4*(C5*C5)*C6*t4 +
7535 24.*C5*C6*Ekaon*t4 -
7536 24.*C4*C5*C6*Ekaon*t4 +
7537 8.*(C5*C5)*C6*Ekaon*t4 -
7538 8.*C4*(C5*C5)*C6*Ekaon*t4 -
7543 2.*amLam*C5*C6*
7545 C5*
7548 C4*
7551 4.*C3*C4*
7553 3.*Enu*t2 +
7554 amLam*apkk1*(C5*C5)*C6*t4 +
7555 amSig*apkk1*(C5*C5)*C6*t4 +
7556 2.*apkk1*(C5*C5)*C6*Ekaon*t4 -
7559 3.*C2*C7*t5 - 6.*C2*C4*C7*t5 -
7560 3.*C8*C9*t6 + 6.*C4*C8*C9*t6) +
7561 C5*C6*t4*
7562 (36.*C5*C6*Ekaon*t4 +
7563 24.*(C5*C5)*C6*Ekaon*t4 +
7564 4.*(C5*C5*C5)*C6*Ekaon*t4 +
7565 2.*amLam*C5*C6*
7566 (9. + 6.*C5 + (C5*C5) +
7570 2.*C3*C5*
7572 2.*amLam*apkk1*(C5*C5)*C6*t4 -
7575 2.*Ekaon*
7576 (6.*t2 + apkk1*(C5*C5)*C6*t4)) +
7577 9.*C2*C7*t5 + 6.*C2*C5*C7*t5 -
7578 9.*C8*C9*t6 - 6.*C5*C8*C9*t6)))))
7579 + 2.*(akpk*akpk)*
7582 (t3*t3) +
7584 t4 +
7586 + 3.*Elep*
7588 (C2*C7*t5 - 1.*C8*C9*t6) +
7589 amSig*C1*C4*t3*
7591 Fm1*
7592 (12.*(-1. + f)*t1 +
7594 3.*(-1. + C4)*(C2*C7*t5 - 1.*C8*C9*t6))
7595 - 1.*amLam*C5*C6*t4*
7597 Fm2*
7598 (4.*(-1. + f)*t1 +
7600 (3. + C5)*(C2*C7*t5 - 1.*C8*C9*t6))) -
7604 3.*amSig*C1*C4*t3*
7605 (2.*C3*C4*
7606 (-6.*t1 - 3.*t2 +
7610 Fm1*
7614 + 3.*C2*C7*t5 - 3.*C8*C9*t6)) +
7615 amLam*C5*C6*t4*
7616 (2.*C3*C5*
7617 (-6.*t1 - 3.*t2 +
7621 Fm2*
7625 - 3.*C2*C7*t5 + 3.*C8*C9*t6))) +
7629 C5*C6*t4*
7630 (C3*
7631 (6.*C5*(2.*t1 + t2) -
7633 Fm2*
7635 3.*C2*C7*t5 - 3.*C8*C9*t6)) +
7636 3.*C1*C4*t3*
7637 (2.*C3*C4*
7638 (6.*t1 + 3.*t2 -
7640 Fm1*
7642 3.*C2*C7*t5 + 3.*C8*C9*t6))) +
7643 3.*Elep*
7644 (C5*C6*t4*
7645 (-1.*(3. + C5)*
7646 (C2*C7*t5 - 1.*C8*C9*t6) +
7647 Fm2*
7648 (4.*t1 - 4.*f*t1 +
7649 2.*C2*C7*Ekaon*t5 -
7650 2.*C8*C9*Ekaon*t6)) +
7651 C1*C4*t3*
7655 t4 + 3.*C2*C7*t5 - 3.*C2*C4*C7*t5 -
7656 3.*C8*C9*t6 + 3.*C4*C8*C9*t6 +
7657 2.*Fm1*
7658 (6.*(-1. + f)*t1 +
7659 Ekaon*
7660 (amLam*(-3. + C5)*C5*C6*t4 -
7661 1.*amSig*(-3. + C5)*C5*C6*t4 -
7662 3.*C2*C7*t5 + 3.*C8*C9*t6)))))) -
7664 (amLam*C5*C6*
7666 1.*Fm2*
7667 (4.*(1. + f)*t1 -
7669 amSig*C1*C4*t3*
7671 Fm1*
7672 (12.*(1. + f)*t1 -
7674 2.*am*
7676 C1*C4*t3*
7679 1.*Fm1*
7680 (6.*(1. + f)*t1 -
7682 C6*Ekaon*t4)))) +
7684 (t3*t3) +
7686 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
7687 3.*amSig*C1*C4*t3*
7688 (2.*C3*C4*
7689 (-6.*t1 - 3.*t2 +
7694 1.*amLam*C5*C6*t4*
7695 (2.*C3*C5*
7696 (-6.*t1 - 3.*t2 +
7700 3.*Fm2*(C2*C7*t5 - 1.*C8*C9*t6)) +
7702 C5*C6*t4*
7703 (C3*
7704 (6.*C5*(2.*t1 + t2) -
7707 3.*C1*C4*t3*
7708 (2.*C3*C4*
7709 (6.*t1 + 3.*t2 -
7711 3.*Fm1*(C2*C7*t5 - 1.*C8*C9*t6)))) +
7713 72.*amSig*C1*(C4*C4)*t1*t3 +
7730 (t3*t3) +
7732 (t3*t3) +
7734 (t3*t3) -
7737 (t3*t3) - 72.*amLam*C5*C6*t1*t4 -
7738 24.*amLam*(C5*C5)*C6*t1*t4 -
7755 t4 -
7757 t4 +
7763 t3*t4 -
7765 t4 +
7767 t3*t4 +
7773 t3*t4 -
7775 t4 -
7777 t3*t4 -
7779 t3*t4 +
7782 Fm2*t3*t4 +
7794 (t4*t4) -
7796 (t4*t4) -
7799 (t4*t4) -
7805 (9.*(C1*C1)*(C4*C4)*
7806 (2.*((-1.+C4)*(-1.+C4))*Ekaon +
7809 *(t3*t3) +
7810 C5*C6*t4*
7811 (2.*C5*((3. + C5)*(3. + C5))*C6*Ekaon*t4 +
7814 3.*Fm2*
7815 (8.*(1. + f)*t1 -
7817 + 3.*C1*C4*t3*
7818 (C5*C6*
7820 2.*Ekaon*
7821 (-6. - 2.*C5 + 2.*C4*(3. + C5) +
7823 Fm1*
7824 (24.*(1. + f)*t1 -
7825 1.*C5*C6*
7828 2.*Ekaon*
7831 2.*C3*
7832 ((3.*am*C1*(C4*C4)*t3 +
7833 3.*amSig*C1*(C4*C4)*t3 +
7836 2.*apkk1*
7837 (6.*t1 + 3.*t2 -
7844 Enu*
7846 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
7848 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
7849 (-6.*t2 +
7851 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
7854 (C4*C4*C4*C4)*(t3*t3) +
7855 (C5*C5)*C6*t4*
7856 (-6.*t2 +
7858 - 6.*C1*(C4*C4)*t3*
7859 (3.*t2 +
7861 C6*t4)) -
7862 2.*am*
7864 (t3*t3) +
7866 (t3*t3) +
7867 6.*C1*(C4*C4)*t3*
7868 (-3.*Ekaon*t2 +
7871 3.*amSig*C1*(C4*C4)*t3*
7875 (C5*C5)*C6*t4*
7880 12.*t1*
7881 (6.*t2 +
7883 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) -
7885 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
7886 am*
7887 (3.*amSig*C1*(C4*C4)*t3 +
7888 6.*C1*(C4*C4)*Ekaon*t3 +
7890 - 108.*C2*C7*f*t1*t5 -
7901 18.*(C2*C2)*(C7*C7)*Ekaon*(t5*t5) -
7902 18.*(C2*C2)*(C7*C7)*Enu*(t5*t5) +
7903 108.*C8*C9*f*t1*t6 +
7914 36.*C2*C7*C8*C9*Ekaon*t5*t6 +
7915 36.*C2*C7*C8*C9*Enu*t5*t6 +
7916 18.*(C8*C8)*(C9*C9)*Ekaon*(t6*t6) -
7917 18.*(C8*C8)*(C9*C9)*Enu*(t6*t6) +
7918 am*(18.*(C1*C1)*(C4*C4)*
7924 C4*
7927 (t3*t3) +
7928 2.*C5*C6*t4*
7929 (-12.*
7938 t4 +
7941 C5*
7944 - 9.*C2*C7*Ekaon*t5 +
7945 3.*C2*C5*C7*Ekaon*t5 -
7946 3.*apkk1*C2*C7*Fm2*t5 +
7948 9.*C8*C9*Ekaon*t6 -
7949 3.*C5*C8*C9*Ekaon*t6 +
7950 3.*apkk1*C8*C9*Fm2*t6 -
7952 3.*C1*C4*t3*
7953 (24.*
7960 4.*C5*
7967 4.*C4*
7970 2.*
7974 (2.*Ekaon*
7976 amSig*
7979 + 3.*
7982 (C2*C7*t5 - 1.*C8*C9*t6))))) +
7993 (t3*t3) +
8003 t4 -
8005 t4 +
8007 (C5*C5)*C6*t3*t4 -
8009 t4 +
8011 t4 -
8013 t4 +
8015 t4 -
8017 t3*t4 +
8019 t3*t4 -
8021 t3*t4 +
8023 t3*t4 +
8025 t4 -
8027 t4 +
8029 t4 -
8031 t4 -
8033 t3*t4 +
8035 t3*t4 -
8037 t3*t4 +
8047 (t4*t4) -
8049 (9.*(C1*C1)*(C4*C4)*
8050 (-1. + 2.*C4 +
8052 6.*C1*C4*C5*
8053 (3. + C5 +
8055 t3*t4 +
8056 (C5*C5)*
8057 (-9. - 6.*C5 +
8061 (3.*(C1*C1)*(C4*C4)*
8062 (3. - 2.*C4 + 3.*(C4*C4))*(t3*t3) +
8063 2.*C1*C4*C5*
8064 (-9. - 1.*C5 + 3.*C4*(1. + C5))*C6*
8065 t3*t4 +
8066 (C5*C5)*(9. + 2.*C5 + (C5*C5))*(C6*C6)*
8067 (t4*t4)) -
8068 4.*Ekaon*
8069 (9.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
8070 (t3*t3) +
8071 3.*C1*C4*C5*C6*
8072 (2.*(1. + C4)*(-3. + C5)*Ekaon +
8073 amSig*
8077 amLam*
8081 ((-3. + C5)*(-3. + C5))*(C5*C5)*(C6*C6)*Ekaon*(t4*t4)
8082 ) -
8084 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
8085 (-12.*t2 +
8087 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))) -
8088 12.*t1*
8089 (-2.*
8091 (3.*C1*C4*(1. + C4)*t3 +
8092 (-3. + C5)*C5*C6*t4) -
8093 1.*apkk1*(1. + f)*
8097 (3.*C1*C4*
8100 C5*C6*
8103 am*
8104 (-3.*C1*C4*
8107 3.*amSig*C1*C4*
8110 C5*C6*
8113 amLam*
8117 ((2. + f)*
8122 C3*
8123 (12.*t2 -
8125 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))))
8128 18.*amSig*apkk1*C1*C2*C4*C7*t3*t5 -
8129 18.*amSig*apkk1*C1*C2*(C4*C4)*C7*t3*t5 -
8134 18.*amLam*apkk1*C2*C5*C6*C7*t4*t5 -
8135 6.*amLam*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
8139 36.*apkk1*(C2*C2)*(C7*C7)*(t5*t5) -
8142 18.*amSig*apkk1*C1*C4*C8*C9*t3*t6 +
8143 18.*amSig*apkk1*C1*(C4*C4)*C8*C9*t3*t6 +
8148 18.*amLam*apkk1*C5*C6*C8*C9*t4*t6 +
8149 6.*amLam*apkk1*(C5*C5)*C6*C8*C9*t4*t6 -
8153 72.*apkk1*C2*C7*C8*C9*t5*t6 +
8155 36.*apkk1*(C8*C8)*(C9*C9)*(t6*t6) -
8156 1.*am*
8159 (t3*t3) +
8161 (t3*t3) +
8162 6.*amSig*C1*C4*t3*
8163 (C5*C6*Ekaon*
8165 2.*C4*Fm2)*t4 +
8166 2.*C3*C4*
8167 (-3.*t2 +
8170 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))) +
8171 3.*C1*C4*t3*
8174 2.*C4*Fm2)*t4 +
8175 4.*C3*C4*
8176 (-6.*Ekaon*t2 +
8180 (C2*C7*t5 - 1.*C8*C9*t6)) +
8181 C5*C6*t4*
8183 4.*C3*C5*Ekaon*
8185 4.*amLam*C3*C5*
8188 (C2*C7*t5 - 1.*C8*C9*t6))) +
8189 2.*
8192 (t3*t3) +
8194 (t3*t3) -
8196 t4 +
8198 t3*t4 -
8200 t3*t4 +
8202 Ekaon*t3*t4 -
8204 Fm1*t3*t4 +
8206 Fm1*t3*t4 -
8208 Fm1*t3*t4 +
8210 Fm1*t3*t4 +
8212 t3*t4 -
8214 t3*t4 -
8216 Fm2*t3*t4 +
8218 Fm2*t3*t4 +
8220 (t4*t4) +
8222 (t4*t4) +
8224 9.*apkk1*C1*C2*C4*C7*t3*t5 +
8225 9.*apkk1*C1*C2*(C4*C4)*C7*t3*t5 -
8227 t5 - 9.*apkk1*C2*C5*C6*C7*t4*t5 +
8228 3.*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
8230 t5 - 9.*apkk1*C1*C4*C8*C9*t3*t6 -
8231 9.*apkk1*C1*(C4*C4)*C8*C9*t3*t6 +
8233 t6 + 9.*apkk1*C5*C6*C8*C9*t4*t6 -
8234 3.*apkk1*(C5*C5)*C6*C8*C9*t4*t6 -
8236 t6 +
8238 (18.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
8239 (t3*t3) +
8240 3.*amSig*C1*C4*t3*
8241 (6.*C1*C4*(1. + (C4*C4))*t3 +
8242 C5*C6*
8246 + C5*C6*t4*
8247 (2.*amLam*C5*(9. + (C5*C5))*C6*t4 +
8248 18.*C5*C6*Ekaon*t4 -
8249 12.*(C5*C5)*C6*Ekaon*t4 +
8250 2.*(C5*C5*C5)*C6*Ekaon*t4 +
8251 9.*C2*C7*t5 - 9.*C8*C9*t6) +
8252 3.*C1*C4*t3*
8253 (-12.*(1. + C4)*C5*C6*Ekaon*t4 +
8254 4.*(1. + C4)*(C5*C5)*C6*Ekaon*t4 +
8255 amLam*C5*C6*
8259 3.*C2*C7*t5 + 3.*C8*C9*t6))))) -
8260 2.*apkk1*
8269 t4 +
8271 t4 +
8273 t4 +
8280 (9.*(C1*C1)*(C4*C4)*(1. + (C4*C4))*(t3*t3) +
8281 6.*C1*C4*C5*(-3. + C4*C5)*C6*t3*t4 +
8282 (C5*C5)*(9. + (C5*C5))*(C6*C6)*(t4*t4)) +
8283 72.*C2*C7*f*t1*t5 +
8291 18.*(C2*C2)*(C7*C7)*Ekaon*(t5*t5) +
8292 18.*(C2*C2)*(C7*C7)*Enu*(t5*t5) -
8293 72.*C8*C9*f*t1*t6 -
8301 36.*C2*C7*C8*C9*Ekaon*t5*t6 -
8302 36.*C2*C7*C8*C9*Enu*t5*t6 -
8303 18.*(C8*C8)*(C9*C9)*Ekaon*(t6*t6) +
8304 18.*(C8*C8)*(C9*C9)*Enu*(t6*t6) -
8305 3.*amSig*C1*C4*t3*
8316 C4*
8317 (-24.*t1 +
8320 (6.*C1*t3 +
8323 - 1.*am*
8325 (t3*t3) -
8326 3.*C1*C4*t3*
8336 3.*C2*C7*Enu*t5 +
8338 3.*C8*C9*Enu*t6 -
8340 C4*
8342 Enu*
8345 3.*C2*C7*t5 - 3.*C8*C9*t6))) +
8346 C5*C6*t4*
8348 3.*C5*
8351 3.*
8353 Enu*
8354 (-3.*C2*C7*t5 + 3.*C8*C9*t6 +
8355 2.*Fm2*
8356 (2.*t1 + C2*C7*Ekaon*t5 -
8357 1.*C8*C9*Ekaon*t6))))))))));
8358
8359 if (sol <= 0.0)
8360 std::cout << "Check Matrix and def. of Scalar Products" << std::endl;
8361
8362 return sol;
8363}
References alepvec, am, amEta, amk, aml, amLam, ampi, amSig, angkq, aq0, aqvec, d, Ekaon, Elep, Enu, f, Fm1, Fm2, fpi, g, pkvec, and Vus.
Referenced by diffxsec().
◆ diffxsec()
| double singlekaon_xsec::diffxsec | ( | double | Tlep, |
| double | Tkaon, | ||
| double | theta, | ||
| double | phikq ) |
Definition at line 90 of file singlekaon_xsec.cxx.
90 {
91
93 if (Tkaon > tkmax) return 0.;
94
97
98 double a1=0., check=0., amat2=0.;
99
106
107 double ddifflep;
108 if (fabs(check) <= 1.0) {
109 angkq = check;
111 amat2 = Amatrix_NN(theta, phikq);
113 amat2 = Amatrix_NP(theta, phikq);
115 amat2 = Amatrix_PP(theta, phikq);
116 else
118
120
121 } else {
122 ddifflep = 0.0;
123 }
124
126}
double Amatrix_NN(double theta, double phikq)
Definition singlekaon_xsec.cxx:129
double Amatrix_PP(double theta, double phikq)
Definition singlekaon_xsec.cxx:4691
double Amatrix_NP(double theta, double phikq)
Definition singlekaon_xsec.cxx:975
References alepvec, am, Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), amk, aml, angkq, aq0, aqvec, Ekaon, Elep, Enu, GeVtocm, ik, pi, and pkvec.
◆ init()
| void singlekaon_xsec::init | ( | double | Etot, |
| int | type, | ||
| int | reac ) |
Definition at line 3 of file singlekaon_xsec.cxx.
3 {
4 Enu = Etot;
5 ilep = type;
6 ik = reac;
7
8 bool verbose = false;
9
10 // Fundamental physics parameters, see PDG booklet
11 pi = std::acos(-1.0);
13 //am=0.939565346; // Nucleon mass (was always set to neutron?)
18
19 // For help on these numbers, see [arXiv:1004:5484]
29
30 // Set the lepton mass
32 aml = 0.510998928*1.e-3;
34 aml = 0.1056583715;
36 aml = 1.77682;
37 } else {
39 return;
40 }
41
42 if (verbose) {
43 std::cout << "Cross-section calculation initialised." << std::endl;
45 }
46
47 std::string lepstr;
51 else lepstr="ERROR";
52
53 // Set reaction parameters
55 if (verbose) {
56 std::cout << " - Reaction:\t nu + n ---> " << lepstr << " + n + K" << std::endl;
57 }
58 amSig = 1.197449;
59 amk = 0.493677;
60 ampi = 0.1349766;
62 if (verbose) {
63 std::cout << " - Reaction:\t nu + n ---> " << lepstr << " + p + K0" << std::endl;
64 }
65 amSig = 1.192642;
66 amk = 0.497614;
67 ampi = 0.13957018;
69 if (verbose) {
70 std::cout << " - Reaction:\t nu + p ---> " << lepstr << " + p + K" << std::endl;
71 }
72 amSig = 1.192642;
73 amk = 0.493677;
74 ampi = 0.134766;
75 } else {
77 return;
78 }
79
81 if (verbose) {
82 std::cout << std::fixed;
84 std::cout << std::scientific;
85 }
86
87 return;
88}
References am, amEta, amk, aml, amLam, ampi, amSig, amun, amup, d, Enu, f, Fm1, Fm2, fpi, g, GeVtocm, ik, ilep, pi, threshold, and Vus.
Referenced by main().
Member Data Documentation
◆ alepvec
|
private |
Definition at line 21 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ am
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), diffxsec(), and init().
◆ amEta
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_PP(), and init().
◆ amk
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), diffxsec(), and init().
◆ aml
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), diffxsec(), and init().
◆ amLam
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NP(), Amatrix_PP(), and init().
◆ ampi
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ amSig
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ amun
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by init().
◆ amup
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by init().
◆ angkq
|
private |
Definition at line 21 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ aq0
|
private |
Definition at line 21 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ aqvec
|
private |
Definition at line 21 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ d
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ Ekaon
|
private |
Definition at line 15 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ Elep
|
private |
Definition at line 21 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ Enu
|
private |
Definition at line 15 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), diffxsec(), and init().
◆ f
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ Fm1
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ Fm2
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by Amatrix_NP(), Amatrix_PP(), and init().
◆ fpi
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ g
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
◆ GeVtocm
|
private |
Definition at line 18 of file singlekaon_xsec.h.
Referenced by diffxsec(), and init().
◆ ik
|
private |
Definition at line 12 of file singlekaon_xsec.h.
Referenced by diffxsec(), and init().
◆ ilep
|
private |
Definition at line 12 of file singlekaon_xsec.h.
Referenced by init().
◆ pi
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by diffxsec(), and init().
◆ pkvec
|
private |
Definition at line 15 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and diffxsec().
◆ threshold
| double singlekaon_xsec::threshold |
Definition at line 26 of file singlekaon_xsec.h.
Referenced by init().
◆ Vus
|
private |
Definition at line 9 of file singlekaon_xsec.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and init().
The documentation for this class was generated from the following files:
Generated on for GENIEGenerator by
1.14.0