genie::AlamSimoAtharVacasSKPXSec2014 Class Reference
Differential cross section model for single kaon production. More...
#include <AlamSimoAtharVacasSKPXSec2014.h>
Inheritance diagram for genie::AlamSimoAtharVacasSKPXSec2014:
Collaboration diagram for genie::AlamSimoAtharVacasSKPXSec2014:
Public Member Functions | |
| AlamSimoAtharVacasSKPXSec2014 () | |
| AlamSimoAtharVacasSKPXSec2014 (string config) | |
| virtual | ~AlamSimoAtharVacasSKPXSec2014 () |
| double | XSec (const Interaction *i, KinePhaseSpace_t k) const |
| Compute the cross section for the input interaction. | |
| double | Integral (const Interaction *i) const |
| bool | ValidProcess (const Interaction *i) const |
| Can this cross section algorithm handle the input process? | |
| void | Configure (const Registry &config) |
| void | Configure (string param_set) |
| Public Member Functions inherited from genie::XSecAlgorithmI | |
| virtual | ~XSecAlgorithmI () |
| virtual bool | ValidKinematics (const Interaction *i) const |
| Is the input kinematical point a physically allowed one? | |
| Public Member Functions inherited from genie::Algorithm | |
| virtual | ~Algorithm () |
| virtual void | FindConfig (void) |
| virtual const Registry & | GetConfig (void) const |
| Registry * | GetOwnedConfig (void) |
| virtual const AlgId & | Id (void) const |
| Get algorithm ID. | |
| virtual AlgStatus_t | GetStatus (void) const |
| Get algorithm status. | |
| virtual bool | AllowReconfig (void) const |
| virtual AlgCmp_t | Compare (const Algorithm *alg) const |
| Compare with input algorithm. | |
| virtual void | SetId (const AlgId &id) |
| Set algorithm ID. | |
| virtual void | SetId (string name, string config) |
| const Algorithm * | SubAlg (const RgKey ®istry_key) const |
| void | AdoptConfig (void) |
| void | AdoptSubstructure (void) |
| virtual void | Print (ostream &stream) const |
| Print algorithm info. | |
Private Member Functions | |
| void | LoadConfig (void) |
| double | Amatrix_NN (double theta, double phikq) const |
| double | Amatrix_NP (double theta, double phikq) const |
| double | Amatrix_PP (double theta, double phikq) const |
Private Attributes | |
| const XSecIntegratorI * | fXSecIntegrator |
| cross section integrator | |
| double | amLam |
| double | amEta |
| double | Vus |
| double | fpi |
| double | d |
| double | f |
| double | g |
| double | amup |
| double | amun |
| double | Fm1 |
| double | Fm2 |
| int | leptonPDG |
| int | reactionType |
| double | aml |
| double | amSig |
| double | amk |
| double | ampi |
| double | am |
| double | Enu |
| double | Ekaon |
| double | pkvec |
| double | Elep |
| double | alepvec |
| double | aqvec |
| double | angkq |
| double | aq0 |
Additional Inherited Members | |
| Static Public Member Functions inherited from genie::Algorithm | |
| static string | BuildParamVectKey (const std::string &comm_name, unsigned int i) |
| static string | BuildParamVectSizeKey (const std::string &comm_name) |
| static string | BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j) |
| static string | BuildParamMatRowSizeKey (const std::string &comm_name) |
| static string | BuildParamMatColSizeKey (const std::string &comm_name) |
| Protected Member Functions inherited from genie::XSecAlgorithmI | |
| XSecAlgorithmI () | |
| XSecAlgorithmI (string name) | |
| XSecAlgorithmI (string name, string config) | |
| Protected Member Functions inherited from genie::Algorithm | |
| Algorithm () | |
| Algorithm (string name) | |
| Algorithm (string name, string config) | |
| void | Initialize (void) |
| void | DeleteConfig (void) |
| void | DeleteSubstructure (void) |
| Registry * | ExtractLocalConfig (const Registry &in) const |
| Registry * | ExtractLowerConfig (const Registry &in, const string &alg_key) const |
| Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key. | |
| template<class T> | |
| bool | GetParam (const RgKey &name, T &p, bool is_top_call=true) const |
| template<class T> | |
| bool | GetParamDef (const RgKey &name, T &p, const T &def) const |
| template<class T> | |
| int | GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const |
| Handle to load vectors of parameters. | |
| int | GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const |
| template<class T> | |
| int | GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const |
| Handle to load matrix of parameters. | |
| template<class T> | |
| int | GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const |
| int | GetParamMatKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const |
| int | AddTopRegistry (Registry *rp, bool owns=true) |
| add registry with top priority, also update ownership | |
| int | AddLowRegistry (Registry *rp, bool owns=true) |
| add registry with lowest priority, also update ownership | |
| int | MergeTopRegistry (const Registry &r) |
| int | AddTopRegisties (const vector< Registry * > &rs, bool owns=false) |
| Add registries with top priority, also udated Ownerships. | |
| Protected Attributes inherited from genie::Algorithm | |
| bool | fAllowReconfig |
| bool | fOwnsSubstruc |
| true if it owns its substructure (sub-algs,...) | |
| AlgId | fID |
| algorithm name and configuration set | |
| vector< Registry * > | fConfVect |
| vector< bool > | fOwnerships |
| ownership for every registry in fConfVect | |
| AlgStatus_t | fStatus |
| algorithm execution status | |
| AlgMap * | fOwnedSubAlgMp |
| local pool for owned sub-algs (taken out of the factory pool) | |
Detailed Description
Differential cross section model for single kaon production.
- References:\n Weak Kaon Production off the Nucleon
- M. Rafi Alam, I. Ruiz Simo, M. Sajjad Athar and M.J. Vicente Vacas Physical Review D82 (2010) 033001; arXiv:1004.5484 [hep-ph]
- Created:\n 2014-02-14
- License:\n Copyright (c) 2003-2025, The GENIE Collaboration
- For the full text of the license visit http://copyright.genie-mc.org
Definition at line 31 of file AlamSimoAtharVacasSKPXSec2014.h.
Constructor & Destructor Documentation
◆ AlamSimoAtharVacasSKPXSec2014() [1/2]
| AlamSimoAtharVacasSKPXSec2014::AlamSimoAtharVacasSKPXSec2014 | ( | ) |
Definition at line 32 of file AlamSimoAtharVacasSKPXSec2014.cxx.
32 :
34{
35
36}
References genie::XSecAlgorithmI::XSecAlgorithmI().
◆ AlamSimoAtharVacasSKPXSec2014() [2/2]
| AlamSimoAtharVacasSKPXSec2014::AlamSimoAtharVacasSKPXSec2014 | ( | string | config | ) |
Definition at line 38 of file AlamSimoAtharVacasSKPXSec2014.cxx.
References genie::XSecAlgorithmI::XSecAlgorithmI().
◆ ~AlamSimoAtharVacasSKPXSec2014()
|
virtual |
Definition at line 44 of file AlamSimoAtharVacasSKPXSec2014.cxx.
45{
46
47}
Member Function Documentation
◆ Amatrix_NN()
|
private |
Definition at line 345 of file AlamSimoAtharVacasSKPXSec2014.cxx.
346{
347
348 double sol = 0.;
349
350 double akk1=0., zdotq=0., qdotpk=0., akcrosk1=0., qcrospk=0.;
351 double zdotpk=0., azpk=0., aqkaon=0., akpk=0., apkk1=0.;
352 double C1=0., C2=0., C3=0., C4=0., /*C5=0., C6=0.,*/ C7=0., C8=0., C9=0.;
353 double aq2=0., gform=0., con=0., t1=0., t2=0., t3=0., /*t4=0.,*/ t5=0., t6=0.;
354
363 akpk=azpk + aqkaon/2.0;
364 apkk1=azpk - aqkaon/2.0;
365
370 //C5=d+3.*f;
371 //C6=1.0/(am*am+amk*amk-2.0*am*Ekaon-amLam*amLam);
376 gform=1.0/pow(1.0-aq2/(1.0*1.0),4);
377
379
380 t1=1.0;
381 t2=1.0;
382 t3=1.0;
383 //t4=1.0;
384 t5=1.0;
385 t6=1.0;
386
387 sol = con*gform*(4.*am*(-2.*(akk1*akk1)*
394 1.*apkk1)*C1*(C4*C4)*t3) -
396 (-2.*t2 +
397 C1*(C4*C4)*
406 Ekaon*(3.*(t2*t2) -
408 t2*t3 +
409 4.*(C1*C1)*(C4*C4*C4*C4)*
415 *(t3*t3))) -
416 8.*C1*C4*Fm1*t3*
420 (2.*t1 +
421 C1*
423 am*
425 2.*akpk*Fm1 +
426 (-2.*apkk1 +
430 (-2. +
433 amSig*
434 (-4.*Ekaon +
435 2.*(akpk - 1.*apkk1)*Fm1 +
437 ))*t3) +
438 2.*Ekaon*
440 (-1.*akpk + apkk1 -
442 Fm1*t3 -
444 amSig*
445 (t1 +
449 (C2*C7*t5 + C8*C9*t6))) +
454 C1*(C4*C4)*t3*
455 (2.*C1*(1. + (C4*C4))*
464 (t2 +
465 2.*C1*(C4*C4)*
467 ))) +
472 C4*
476 (-1.*C3*t2 -
477 2.*C1*C4*
478 (2. +
479 C4*
483 C1*C4*t3*
484 (-2.*C1*C4*
488 2.*am*
490 Fm1 -
501 C3*C4*
505 2.*C1*(C4*C4)*
510 t3) -
511 1.*
514 *(C2*C7*t5 + C8*C9*t6))) +
515 apkk1*(-2.*((-1.+C4)*(-1.+C4))*(t1*t1) +
516 2.*C1*C4*
519 2.*am*
523 Ekaon) -
526 2.*am*
528 Fm1)*t1*t3 +
529 C1*C4*t3*
530 (2.*C1*C4*
535 C4 -
536 1.*
540 2.*am*
542 2.*
544 C4 +
548 (-1. +
553 *(C2*C7*t5 + C8*C9*t6))) +
554 2.*(apkk1*apkk1)*C1*C4*t3*
556 (C2*C7*t5 + C8*C9*t6) +
559 2.*am*
560 (t1 - 1.*C4*t1 + C2*C7*Ekaon*t5 +
561 C8*C9*Ekaon*t6)))) +
562 2.*(akpk*akpk)*C1*C4*t3*
564 (C3*C4*(2.*t1 + t2 -
565 4.*C1*(C4*C4)*
567 + Fm1*
568 (4.*C1*C4*
570 t3 + C2*C7*t5 + C8*C9*t6)) -
572 (C2*C7*t5 + C8*C9*t6) +
575 2.*am*
576 ((1. + C4)*t1 -
577 1.*Ekaon*(C2*C7*t5 + C8*C9*t6))))) +
581 (C3*C4*(2.*t1 + t2 -
582 4.*C1*(C4*C4)*
584 - 1.*Fm1*(C2*C7*t5 + C8*C9*t6)) +
590 (2.*((-1.+C4)*(-1.+C4))*Ekaon +
592 (t3*t3) +
594 (2.*C1*((-1.+C4)*(-1.+C4))*C4*Ekaon*t3 +
597 Fm1*
598 (-2.*(-1. + C4)*t1 +
600 1.*C3*
601 (2.*Enu*
602 (t1*t2 +
603 C1*(C4*C4)*
606 (2.*t1 + t2)*t3 +
607 2.*(C1*C1)*(C4*C4*C4*C4)*
611 2.*am*
616 2.*apkk1*
617 (2.*t1 + t2 -
618 4.*C1*(C4*C4)*
620 t3))) -
621 1.*(C2*C7*t5 + C8*C9*t6)*
622 (-3.*C4*t1 +
624 + 2.*amSig*C1*C4*t3*
625 (2.*(-1.*((-1.+C4)*(-1.+C4)) +
628 (((-1.+C4)*(-1.+C4)) +
631 (C2*C7*t5 + C8*C9*t6)) +
632 2.*am*C1*C4*t3*
633 (-2.*
634 (((-1.+C4)*(-1.+C4)) +
637 Fm1*t3 -
641 (((-1.+C4)*(-1.+C4)) +
642 (4.*amSig*C4 -
645 (apkk1*Fm1 +
647 (C2*C7*t5 + C8*C9*t6)) +
649 (-3.*(C2*C7*t5 + C8*C9*t6) +
650 2.*Fm1*
651 ((2. + C4)*t1 +
652 Enu*(C2*C7*t5 + C8*C9*t6)))) -
653 1.*Elep*(2.*((1.+C4)*(1.+C4))*(t1*t1) +
656 (t3*t3) +
659 Ekaon*(t3*t3) +
664 (t2 +
666 t3)) -
667 2.*t1*
668 (2.*C1*C4*
673 (apkk1 - 1.*apkk1*C4 +
677 t3 +
678 C3*
679 (t2 -
681 t3))) +
682 (C2*C7*t5 + C8*C9*t6)*
684 (amSig*C1*C4*t3 + C2*C7*t5 + C8*C9*t6)
685 + 2.*apkk1*
686 (amSig*C1*C4*(1. + C4)*t3 + C2*C7*t5 +
687 C8*C9*t6)) +
690 (2.*C3*C4*
692 Fm1*(C2*C7*t5 + C8*C9*t6)) -
693 2.*
695 amSig*(C2*C7*t5 + C8*C9*t6) -
696 1.*apkk1*Fm1*(C2*C7*t5 + C8*C9*t6)))
697 + am*C1*C4*t3*
700 4.*C3*C4*Ekaon*
702 amSig*C3*
705 (C2*C7*t5 + C8*C9*t6)) -
706 2.*
708 t3 -
710 (C2*C7*t5 + C8*C9*t6) +
712 (4.*C1*C4*
713 (amSig*(1. + (C4*C4)) +
714 ((1.+C4)*(1.+C4))*Ekaon)*t3 +
715 C2*C7*t5 + C8*C9*t6)))) +
716 2.*apkk1*
718 (t3*t3) -
720 (t3*t3) +
721 (C2*C7*t5 + C8*C9*t6)*
722 (C4*(2.*t1 +
725 (C2*C7*t5 + C8*C9*t6)) +
726 amSig*C1*C4*t3*
730 C4*
733 1.*Enu*(C2*C7*t5 + C8*C9*t6))) +
734 am*C1*C4*t3*
736 C4*
739 1.*Enu*(C2*C7*t5 + C8*C9*t6)) +
740 Enu*
741 (C2*C7*t5 + C8*C9*t6 +
742 2.*Fm1*
743 (t1 -
744 1.*Ekaon*(C2*C7*t5 + C8*C9*t6))))))
747 (C4*C4)*
748 (-1. +
754 ) - 1.*Ekaon*
757 (-2.*C3*(C4*C4)*Ekaon -
761 4.*Fm1*
764 2.*Ekaon*
768 (3.*Ekaon*(t2*t2) -
769 4.*C1*(C4*C4)*
772 4.*(C1*C1)*(C4*C4*C4*C4)*
777 (t3*t3)) +
781 4.*(akpk*akpk)*C1*C4*Ekaon*
783 4.*(apkk1*apkk1)*C1*C4*Ekaon*
787 (C4*C4)*
790 Enu) +
793 4.*amSig*
794 (1. +
795 (C4*C4)*
796 (-1. +
806 2.*C1*
807 (((1.+C4)*(1.+C4))*Elep -
808 1.*((-1.+C4)*(-1.+C4))*Enu -
809 1.*amSig*
810 (-1. + (C4*C4) +
812 Fm1))*t3 +
817 ) - 1.*Ekaon*
818 (-4.*Enu*
819 ((1. - 3.*C4)*Fm1*t1 +
822 (4. +
823 (C4*C4)*
824 (-4. +
828 Fm1*
835 4.*Elep*
838 Fm1*
839 (t1 + 3.*C4*t1 -
841 2.*apkk1*
842 (-2.*(1. + C4)*Fm1*t1 +
844 C1*C4*
848 2.*C1*((1.+C4)*(1.+C4))*Ekaon*t3 +
853 2.*akpk*(C1*C4*
856 2.*Fm1*
857 (t1 - 1.*C4*t1 +
859 C4*(2.*C1*((-1.+C4)*(-1.+C4))*Ekaon*t3 +
863 C1*(C4*C4)*
864 (4.*(-1. + apkk1*C3)*Ekaon +
867 )))) +
873 (t2*t2) +
875 (akpk - 1.*apkk1)*C1*(C4*C4)*t2*t3 -
876 16.*amSig*((akpk - 1.*apkk1)*(akpk - 1.*apkk1))*(C1*C1)*(C4*C4*C4*C4)*
878 8.*C3*(-1.*Ekaon*t1*t2 +
879 C1*(C4*C4)*
881 akpk*amSig*(2.*t1 + t2))*t3 -
882 1.*(C1*C1)*(C4*C4*C4*C4)*
886 (t3*t3)) -
887 1.*(C2*C7*t5 + C8*C9*t6)*
889 - 4.*amSig*C1*C4*t3*
892 (C2*C7*t5 + C8*C9*t6)) +
894 (C4*Fm1*
895 (6.*t1 +
896 amSig*C1*
899 Fm1)*t3) -
901 (C2*C7*t5 + C8*C9*t6))) -
902 1.*Elep*(8.*((1.+C4)*(1.+C4))*(t1*t1) +
904 2.*apkk1)*(C3*C3)*(t2*t2) -
905 16.*C1*C4*
911 + 2.*Fm1*(C2*C7*t5 + C8*C9*t6)) -
917 2.*amSig*C1*C4*(1. + C4)*t3*
918 (C2*C7*t5 + C8*C9*t6) +
919 ((C2*C7*t5 + C8*C9*t6)*(C2*C7*t5 + C8*C9*t6)) +
921 (2.*apkk1*(C3*C3)*C4*
923 Fm1*
925 1.*C2*C7*t5 - 1.*C8*C9*t6))) +
930 (C3*C3)*(t2 + 2.*apkk1*C1*(C4*C4)*t3)) +
932 t3 - 1.*Fm1*(C2*C7*t5 + C8*C9*t6)))) -
933 4.*(8.*(akpk*akpk)*amSig*(C1*C1)*(C4*C4)*
936 (((1.+C4)*(1.+C4))*Ekaon -
940 ((C2*C7*t5 + C8*C9*t6)*
943 2.*apkk1*C1*C4*Fm1*t3*
945 1.*Ekaon*(C2*C7*t5 + C8*C9*t6))) +
946 2.*amSig*C1*C4*t3*
949 apkk1*
950 (-2.*
952 t1 +
954 (((1.+C4)*(1.+C4)) -
956 (1. + C4)*Ekaon*(C2*C7*t5 + C8*C9*t6)))
957 - 2.*akpk*C1*C4*t3*
959 (((-1.+C4)*(-1.+C4))*Ekaon +
962 (2.*(1. + C4)*
964 Ekaon*(C2*C7*t5 + C8*C9*t6)) +
965 amSig*
966 (2.*(((-1.+C4)*(-1.+C4)) -
969 (((-1.+C4)*(-1.+C4)) +
972 (C2*C7*t5 + C8*C9*t6)))) +
973 Enu*(8.*((-1.+C4)*(-1.+C4))*(t1*t1) -
978 + 8.*t1*
979 (2.*C1*C4*
985 ((C2*C7*t5 + C8*C9*t6)*(C2*C7*t5 + C8*C9*t6))) +
990 (C3*C3)*(t2 + 2.*apkk1*C1*(C4*C4)*t3)) +
991 2.*(C1*((-1.+C4)*(-1.+C4))*C4*t3 -
992 1.*Fm1*(C2*C7*t5 + C8*C9*t6))) -
997 (-8.*C1*C3*(C4*C4)*t3*
999 (C3*C3)*
1000 (3.*(t2*t2) + 8.*apkk1*C1*(C4*C4)*t2*t3 +
1002 + 4.*C1*C4*Fm1*t3*
1004 C8*C9*t6)) +
1006 (t3*t3) -
1007 2.*amSig*C1*(-1. + C4)*C4*t3*
1008 (C2*C7*t5 + C8*C9*t6) +
1009 (C2*C7*t5 + C8*C9*t6)*
1010 (4.*apkk1*C1*C4*Fm1*t3 + C2*C7*t5 +
1011 C8*C9*t6)))) +
1013 (t2*t2) -
1014 2.*(-2.*(-1. + (C4*C4))*(t1*t1) -
1015 2.*C1*C4*
1016 (2.*akpk*
1021 2.*apkk1*
1022 (((1.+C4)*(1.+C4)) -
1023 1.*
1024 (amSig*(1. + C4) +
1028 - 1.*Ekaon*
1029 (4.*((1.+C4)*(1.+C4))*Elep +
1030 4.*((-1.+C4)*(-1.+C4))*Enu +
1035 t1*t3 +
1036 C1*C4*t3*
1038 (1. +
1039 (C4*C4)*
1040 (-1. +
1042 (-2. +
1046 )*t3 +
1047 2.*(akpk*akpk)*C4*
1052 (-1.*apkk1*t2 +
1053 2.*Ekaon*
1055 1.*Enu*t2 +
1057 2.*apkk1*Ekaon*
1059 t3 -
1063 (C2*C7*t5 + C8*C9*t6)) -
1065 (2.*(apkk1*apkk1)*(C3*C3)*C4*t2 +
1066 apkk1*
1067 (2.*(C3*C3)*C4*
1072 1.*Fm1*(C2*C7*t5 + C8*C9*t6)) +
1073 Ekaon*
1078 t3 +
1079 (1. +
1081 (C2*C7*t5 + C8*C9*t6))) +
1082 akpk*
1085 2.*
1087 Fm1*t3 +
1091 (((-1.+C4)*(-1.+C4)) +
1092 2.*
1095 )*t3 -
1096 1.*
1097 (4.*apkk1*Fm1 +
1098 Ekaon*
1101 ) +
1103 (C3*
1105 + (C3*C3)*C4*
1107 (t2 +
1108 8.*C1*(C4*C4)*
1109 (apkk1 -
1111 *t3) +
1112 2.*
1113 (2.*apkk1*t2 +
1118 Fm1*
1120 C2*C7*t5 + C8*C9*t6))) +
1123 (t2 +
1124 C1*(C4*C4)*
1128 t3) +
1129 2.*
1130 (-4.*C1*C4*Ekaon*
1131 (((1.+C4)*(1.+C4))*Elep -
1132 1.*((-1.+C4)*(-1.+C4))*Enu)*t3 +
1134 2.*apkk1*C1*C4*
1135 (((1.+C4)*(1.+C4)) +
1139 (-1. + (C4*C4) +
1142 2.*amSig*C1*C4*
1143 (-1.*((-1.+C4)*(-1.+C4))*Enu +
1144 4.*apkk1*C4*Fm1 -
1145 1.*Ekaon*
1146 (-1. + (C4*C4) +
1147 2.*
1149 Fm1) +
1150 Elep*
1152 t3 +
1156 (C2*C7*t5 + C8*C9*t6)) +
1158 (2.*C3*C4*
1159 (-1.*t2 +
1160 2.*C1*(C4*C4)*
1164 (C3*C3)*C4*
1165 (4.*(apkk1*apkk1)*C1*(C4*C4)*t3 +
1166 apkk1*
1167 (t2 -
1171 (-2.*Ekaon*t2 +
1174 2.*Fm1*
1175 (C1*C4*
1177 amSig*
1178 (6. -
1180 + 2.*Ekaon*
1181 (-3. +
1183 + C2*C7*t5 + C8*C9*t6)))))))));
1184
1185 if (sol <= 0.0) {
1187 << "Matrix element is 0 or negative! This is bad and really should throw an exception!";
1188 }
1189
1190 return sol;
1191}
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE,...
Definition Messenger.h:96
double pkvec
Definition AlamSimoAtharVacasSKPXSec2014.h:72
double Ekaon
Definition AlamSimoAtharVacasSKPXSec2014.h:72
double alepvec
Definition AlamSimoAtharVacasSKPXSec2014.h:77
double amSig
Definition AlamSimoAtharVacasSKPXSec2014.h:71
double amEta
Definition AlamSimoAtharVacasSKPXSec2014.h:63
double angkq
Definition AlamSimoAtharVacasSKPXSec2014.h:77
double aqvec
Definition AlamSimoAtharVacasSKPXSec2014.h:77
References alepvec, am, amEta, amk, aml, ampi, amSig, angkq, aq0, aqvec, d, Ekaon, Elep, Enu, f, Fm1, fpi, g, LOG, pERROR, pkvec, and Vus.
Referenced by XSec().
◆ Amatrix_NP()
|
private |
Definition at line 1194 of file AlamSimoAtharVacasSKPXSec2014.cxx.
1195{
1196
1197 double sol = 0.;
1198
1199 double akk1=0., zdotq=0., qdotpk=0., akcrosk1=0., qcrospk=0.;
1200 double zdotpk=0., azpk=0., aqkaon=0., akpk=0., apkk1=0.;
1201 double C1=0., C2=0., C3=0., C4=0., C5=0., C6=0., C7=0./*, C8=0., C9=0.*/;
1202 double aq2=0., gform=0., con=0., t1=0., t2=0., t3=0., t4=0., t5=0./*, t6=0.*/;
1203
1212 akpk=azpk + aqkaon/2.0;
1213 apkk1=azpk - aqkaon/2.0;
1221 //C8=am/(aml*aml-2.0*akk1+amk*amk-2.*aqkaon-amEta*amEta);
1222 //C9=d - 3.*f;
1224 gform=1.0/pow(1.0-aq2/(1.1*1.1),4);
1225
1227
1228 t1=1.0;
1229 t2=1.0; // !Full Term
1230 t3=1.0;
1231 t4=1.0;
1232 t5=1.0;
1233 //t6=1.0;
1234
1235 sol = con*gform*(0.4444444444444444*am*
1253 72.*akpk*amSig*C1*C4*t1*t3 -
1254 72.*amSig*apkk1*C1*C4*t1*t3 -
1255 72.*akpk*amSig*C1*C2*C4*t1*t3 +
1256 72.*amSig*apkk1*C1*C2*C4*t1*t3 +
1257 72.*akpk*amSig*C1*(C4*C4)*t1*t3 -
1258 72.*amSig*apkk1*C1*(C4*C4)*t1*t3 +
1259 72.*akpk*amSig*C1*C2*(C4*C4)*t1*t3 +
1260 72.*amSig*apkk1*C1*C2*(C4*C4)*t1*t3 -
1293 (C4*C4)*t2*t3 -
1302 Elep*t2*t3 -
1315 (C4*C4*C4*C4)*(t3*t3) -
1318 (C4*C4*C4*C4)*(t3*t3) +
1321 (C4*C4*C4*C4)*(t3*t3) +
1323 (C4*C4*C4*C4)*(t3*t3) -
1326 (C3*C3)*(C4*C4*C4*C4)*(t3*t3) +
1336 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
1339 Ekaon*(t3*t3) +
1342 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
1344 Ekaon*(t3*t3) -
1346 Ekaon*(t3*t3) -
1348 Ekaon*(t3*t3) -
1358 (C3*C3)*(C4*C4*C4*C4)*Elep*(t3*t3) +
1361 (C4*C4*C4*C4)*Elep*(t3*t3) +
1363 (C4*C4*C4*C4)*Elep*(t3*t3) +
1373 (C4*C4*C4*C4)*Enu*(t3*t3) +
1376 (C4*C4*C4*C4)*Enu*(t3*t3) -
1379 (C4*C4*C4*C4)*Enu*(t3*t3) -
1381 (C4*C4*C4*C4)*Enu*(t3*t3) +
1388 Fm1*(t3*t3) +
1394 Fm1*(t3*t3) -
1430 72.*amLam*apkk1*C5*C6*t1*t4 -
1431 72.*akpk*amLam*C2*C5*C6*t1*t4 +
1432 72.*amLam*apkk1*C2*C5*C6*t1*t4 -
1433 24.*akpk*amLam*(C5*C5)*C6*t1*t4 +
1434 24.*amLam*apkk1*(C5*C5)*C6*t1*t4 -
1435 24.*akpk*amLam*C2*(C5*C5)*C6*t1*t4 -
1436 24.*amLam*apkk1*C2*(C5*C5)*C6*t1*t4 +
1499 C6*t3*t4 +
1502 (C5*C5)*C6*t3*t4 +
1505 (C5*C5*C5)*C6*t3*t4 +
1509 C6*t3*t4 +
1511 (C5*C5*C5)*C6*t3*t4 -
1514 (C5*C5*C5)*C6*t3*t4 +
1516 (C5*C5*C5)*C6*t3*t4 -
1519 (C5*C5*C5)*C6*t3*t4 -
1521 (C5*C5*C5)*C6*t3*t4 +
1523 C6*t3*t4 -
1525 (C5*C5*C5)*C6*t3*t4 +
1527 C6*t3*t4 +
1529 (C5*C5*C5)*C6*t3*t4 +
1531 (C5*C5*C5)*C6*t3*t4 +
1537 Ekaon*t3*t4 +
1540 Ekaon*t3*t4 +
1543 Ekaon*t3*t4 +
1545 Ekaon*t3*t4 +
1547 C6*Ekaon*t3*t4 +
1550 Ekaon*t3*t4 +
1552 Ekaon*t3*t4 -
1555 (C5*C5*C5)*C6*Ekaon*t3*t4 -
1557 Ekaon*t3*t4 -
1559 C6*Ekaon*t3*t4 -
1561 (C5*C5*C5)*C6*Ekaon*t3*t4 +
1563 C6*Ekaon*t3*t4 +
1565 (C5*C5*C5)*C6*Ekaon*t3*t4 -
1575 C6*Elep*t3*t4 +
1578 (C4*C4)*(C5*C5*C5)*C6*Elep*t3*t4 -
1580 C6*Elep*t3*t4 +
1582 Elep*t3*t4 -
1584 (C5*C5*C5)*C6*Elep*t3*t4 +
1595 (C4*C4)*(C5*C5)*C6*Enu*t3*t4 +
1598 Enu*t3*t4 +
1601 Enu*t3*t4 -
1603 (C5*C5*C5)*C6*Enu*t3*t4 +
1605 C6*Enu*t3*t4 -
1608 C3*(C4*C4)*(C5*C5*C5)*C6*Enu*t3*t4 +
1613 Fm1*t3*t4 -
1620 C6*Fm1*t3*t4 +
1623 Fm1*t3*t4 +
1626 Fm1*t3*t4 -
1629 Fm1*t3*t4 -
1638 Fm1*t3*t4 -
1641 Fm1*t3*t4 +
1646 Fm2*t3*t4 -
1653 C6*Fm2*t3*t4 +
1656 Fm2*t3*t4 +
1659 Fm2*t3*t4 -
1662 Fm2*t3*t4 +
1671 Fm2*t3*t4 -
1674 Fm2*t3*t4 -
1680 Fm2*t3*t4 +
1683 Fm2*t3*t4 -
1685 Fm2*t3*t4 -
1687 Fm2*t3*t4 -
1690 Fm2*t3*t4 +
1701 Fm2*t3*t4 -
1726 (C6*C6)*(t4*t4) -
1737 Ekaon*(t4*t4) -
1742 Ekaon*(t4*t4) +
1746 Ekaon*(t4*t4) -
1760 (C6*C6)*Elep*(t4*t4) +
1773 Enu*(t4*t4) +
1782 Fm2*(t4*t4) -
1792 Fm2*(t4*t4) -
1823 (1. +
1824 (C4*C4)*
1831 6.*C1*C4*C5*C6*
1836 C4*
1842 (C5*C5)*(C6*C6)*
1849 1.*Ekaon*
1851 ((-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)*(-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)) +
1853 (-9.*(C1*C1)*(C4*C4)*
1854 (2.*C3*(C4*C4)*Ekaon +
1858 (t3*t3) +
1859 6.*C1*C4*C5*C6*
1860 (C3*C4*C5*
1861 (Ekaon*
1863 2.*(akpk - 1.*apkk1)*C5*
1866 (C5*C5)*(C6*C6)*
1867 (-2.*(C5*C5*C5)*Ekaon -
1871 ) +
1873 (apkk1*
1876 akpk*
1879 - 2.*Ekaon*
1880 (3.*C1*C4*
1883 C5*C6*
1884 (-1.*(-3. + C5)*Elep +
1900 12.*apkk1*C4*Fm1 +
1905 4.*amSig*
1906 (1. +
1907 (C4*C4)*
1914 - 1.*Ekaon*
1915 (4. +
1916 (C4*C4)*
1925 2.*Ekaon*
1926 (2.*(akpk*akpk)*
1928 2.*(apkk1*apkk1)*
1930 apkk1*
1931 (2. + 4.*C4 +
1932 (C4*C4)*
1938 1.*akpk*
1939 (2. - 4.*C4 +
1940 (C4*C4)*
1943 (-1. + apkk1*C3 +
1947 2.*Ekaon*
1950 2.*
1951 (((1.+C4)*(1.+C4))*Elep -
1952 1.*((-1.+C4)*(-1.+C4))*Enu) +
1953 amSig*
1954 (-2. +
1963 ))*(t3*t3) +
1964 C5*C6*t4*
1966 Fm2*t1 +
1967 24.*C2*
1968 ((akpk - 1.*apkk1)*Fm2 +
1970 t1 +
1972 (4.*Ekaon*
1973 (akpk +
1976 t4 +
1977 2.*(C5*C5*C5)*C6*
1980 + 2.*Ekaon*
1981 (-1.*akpk + apkk1 +
1985 (-2.*
1986 (3.*akpk - 3.*apkk1 +
1991 1.*Enu)) +
1992 2.*Ekaon*
1993 (2.*(akpk*akpk) +
1994 2.*apkk1*
1997 (apkk1 +
1998 2.*Ekaon*
2000 ) -
2001 1.*akpk*
2003 4.*
2005 )))*t4 -
2006 4.*(C5*C5)*
2009 + C6*
2015 ))*t4 +
2016 3.*akpk*
2019 + C5*
2022 (Enu*
2025 3.*amLam*
2028 Elep*
2032 (4.*amLam*
2036 3.*
2041 Elep*
2044 1.*Ekaon*
2045 (24.*t1 +
2046 C6*
2049 2.*apkk1*
2057 6.*C1*C4*t3*
2061 12.*C2*
2066 12.*akpk*C5*C6*Ekaon*t4 +
2068 12.*apkk1*C5*C6*Ekaon*t4 +
2069 4.*akpk*(C5*C5)*C6*Ekaon*t4 +
2070 4.*apkk1*(C5*C5)*C6*Ekaon*t4 +
2119 t4 +
2121 t4 -
2126 t4 +
2131 t4 -
2156 t4 -
2166 t4 +
2168 Enu*t4 +
2170 Enu*t4 -
2180 6.*akpk*
2187 t4 -
2188 2.*C5*C6*
2190 2.*C5*Enu +
2193 t4 +
2196 1.*(C5*C5*C5)*C6*
2198 1.*C3*(C5*C5)*C6*
2202 t4)) +
2203 Ekaon*
2208 12.*apkk1*C5*C6*t4 +
2210 4.*apkk1*(C5*C5)*C6*t4 +
2218 t4 +
2220 Enu*t4 -
2224 2.*akpk*
2226 2.*C5*(3. + C5)*C6*t4 +
2228 (3.*(C3*C3)*t2 - 2.*(C5*C5*C5)*C6*t4 +
2229 2.*C3*(C5*C5)*(-1. + 2.*apkk1*C5)*C6*
2230 t4)))))) -
2289 2.*am*(36.*(-1. + (C2*C2))*(t1*t1) +
2316 (t3*t3) -
2319 (t3*t3) +
2321 (t3*t3) +
2323 (t3*t3) -
2326 (t3*t3) -
2328 (t3*t3) +
2330 (t3*t3) -
2332 (t3*t3) +
2334 (t3*t3) -
2336 (t3*t3) -
2338 (t3*t3) +
2341 (t3*t3) +
2343 (t3*t3) -
2345 (t3*t3) -
2348 (t3*t3) -
2350 (t3*t3) -
2352 (t3*t3) -
2355 (t3*t3) +
2366 (t3*t3) +
2368 (t3*t3) +
2388 t3*t4 +
2397 t3*t4 -
2403 t3*t4 -
2417 Fm1*t3*t4 -
2420 Fm1*t3*t4 -
2423 t3*t4 -
2426 t3*t4 +
2428 t3*t4 +
2437 t3*t4 +
2440 Fm1*t3*t4 +
2443 Fm1*t3*t4 +
2445 t3*t4 +
2447 t3*t4 +
2450 Fm1*t3*t4 +
2452 t3*t4 +
2454 t3*t4 -
2456 t3*t4 +
2458 t3*t4 +
2460 t3*t4 +
2462 t3*t4 -
2465 t3*t4 -
2468 Fm1*t3*t4 +
2470 t3*t4 +
2472 t3*t4 +
2474 t3*t4 -
2476 t3*t4 +
2479 Fm2*t3*t4 +
2482 Fm2*t3*t4 -
2486 t3*t4 +
2488 t3*t4 -
2497 t3*t4 +
2500 Fm2*t3*t4 -
2503 Fm2*t3*t4 -
2505 t3*t4 -
2507 t3*t4 +
2510 Fm2*t3*t4 +
2513 Fm2*t3*t4 -
2515 t3*t4 -
2517 t3*t4 +
2519 t3*t4 -
2521 t3*t4 -
2524 t3*t4 -
2527 Fm2*t3*t4 -
2529 t3*t4 -
2531 t3*t4 +
2533 t3*t4 -
2535 t3*t4 -
2538 t3*t4 +
2540 Fm2*t3*t4 +
2542 Fm2*t3*t4 +
2544 Fm2*t3*t4 +
2551 t3*t4 +
2553 t3*t4 +
2555 t3*t4 +
2560 Fm2*t3*t4 +
2589 (t4*t4) +
2591 (t4*t4) +
2593 (t4*t4) +
2596 (t4*t4) +
2599 (t4*t4) -
2601 (t4*t4) +
2603 (t4*t4) +
2605 (t4*t4) +
2607 (t4*t4) +
2610 (t4*t4) +
2612 (t4*t4) +
2614 (t4*t4) -
2617 (t4*t4) -
2619 (t4*t4) +
2621 (t4*t4) -
2625 (t4*t4) +
2627 (t4*t4) +
2632 (t4*t4) +
2634 (t4*t4) +
2636 (t4*t4) +
2638 (t4*t4) +
2641 (t4*t4) +
2686 2.*akpk*
2687 (3.*C1*C4*
2692 )*t3 +
2693 C5*C6*
2697 C2*
2700 2.*apkk1*
2701 (3.*C1*C4*
2704 C2*
2707 1.*C5*C6*
2711 )*t4)) +
2717 (t3*t3) + 3.*apkk1*C3*(C5*C5*C5)*C6*t2*t4 +
2725 ((-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)*(-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)) +
2726 akpk*(3.*C1*C3*(C4*C4)*t3 - 1.*(C5*C5*C5)*C6*t4)*
2730 (9.*(C1*C1)*(C3*C3)*(C4*C4*C4*C4)*
2734 (t3*t3) +
2735 (C5*C5*C5)*C6*t4*
2736 (-6.*C3*t2 -
2737 1.*(C5*C5*C5)*C6*
2741 t4) +
2742 3.*C1*C3*(C4*C4)*t3*
2743 (6.*C3*t2 +
2744 (C5*C5*C5)*C6*
2745 (2.*apkk1 +
2748 amLam*
2750 )*t4))) +
2751 36.*akpk*C1*C2*C4*C7*Ekaon*t3*t5 -
2753 36.*apkk1*C1*C2*C4*C7*Ekaon*t3*t5 -
2754 36.*akpk*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 -
2755 36.*apkk1*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 +
2760 144.*akpk*apkk1*C1*C2*C4*C7*Fm1*t3*t5 -
2767 36.*akpk*C2*C5*C6*C7*Ekaon*t4*t5 -
2769 36.*apkk1*C2*C5*C6*C7*Ekaon*t4*t5 +
2770 12.*akpk*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
2771 12.*apkk1*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
2776 48.*akpk*apkk1*C2*C5*C6*C7*Fm2*t4*t5 -
2784 (t3*t3) +
2787 (t3*t3) +
2791 (t3*t3) -
2793 (t3*t3) +
2795 (t3*t3) -
2797 (t3*t3) -
2799 (t3*t3) +
2801 (t3*t3) +
2803 (t3*t3) -
2805 (t3*t3) -
2807 (t3*t3) +
2808 9.*(C3*C3)*
2810 2.*(akpk*akpk)*C1*(C4*C4)*t2*t3 -
2812 4.*akpk*apkk1*C1*(C4*C4)*t2*t3 +
2814 2.*(apkk1*apkk1)*C1*(C4*C4)*t2*t3 +
2820 (t3*t3) -
2823 (t3*t3) -
2825 (t3*t3) +
2828 (t3*t3) +
2832 (C4*C4*C4*C4)*(t3*t3)) -
2833 2.*Elep*
2834 (amSig*(akpk - 1.*apkk1)*C1*(C4*C4)*t3*
2836 Enu*
2839 (C1*C1)*(C4*C4*C4*C4)*(t3*t3))) -
2840 2.*Ekaon*
2841 (-1.*amSig*(akpk - 1.*apkk1)*C1*(C4*C4)*
2842 t2*t3 +
2843 Elep*
2844 ((t2*t2) +
2846 (C4*C4)*t2*t3 +
2848 (akpk - 1.*apkk1)*(C1*C1)*(C4*C4*C4*C4)*(t3*t3)
2849 ) -
2850 1.*Enu*
2851 ((t2*t2) +
2853 (C4*C4)*t2*t3 +
2855 (akpk - 1.*apkk1)*(C1*C1)*(C4*C4*C4*C4)*(t3*t3)
2856 ))) +
2860 t4 -
2862 t4 -
2864 t4 -
2866 t4 -
2868 t3*t4 -
2870 t4 -
2872 t4 -
2874 t4 -
2886 t4 +
2889 t3*t4 -
2891 t4 +
2893 t3*t4 +
2895 t3*t4 -
2897 t4 -
2899 Fm2*t3*t4 +
2901 t3*t4 +
2903 t3*t4 +
2905 t3*t4 -
2907 t4 -
2909 t4 -
2911 t3*t4 -
2922 (t4*t4) +
2924 (t4*t4) -
2926 (t4*t4) -
2928 (t4*t4) +
2930 (t4*t4) +
2932 (t4*t4) +
2935 (t4*t4) +
2938 (t4*t4) -
2940 (t4*t4) +
2943 (t4*t4) +
2945 (t4*t4) -
2947 (t4*t4) -
2949 (t4*t4) -
2954 (t4*t4) -
2961 (t4*t4) -
2964 (t4*t4) -
2966 (t4*t4) -
2968 (t4*t4) +
2970 (t4*t4) +
2972 (t4*t4) -
2974 (t4*t4) -
2976 (t4*t4) +
2977 3.*C3*
2978 (2.*(akpk*akpk)*(C5*C5*C5)*C6*
2981 (3.*C1*(C4*C4)*t3 +
2982 (C5*C5)*C6*
2984 2.*
2986 (t3*t3) +
2987 (C5*C5)*C6*
2988 ((apkk1*apkk1)*C5 -
2990 apkk1*C5*
2993 t4 +
2994 C1*(C4*C4)*Ekaon*t3*
2998 t4 +
3000 t4)) -
3002 (6.*(C1*C1)*(C4*C4*C4*C4)*
3006 (C5*C5)*C6*
3007 (-2. + apkk1*C5 +
3010 + C1*(C4*C4)*t3*
3011 (6.*t2 +
3012 (C5*C5)*C6*
3013 (-4.*(apkk1*apkk1)*C5 +
3018 amLam*
3020 2.*amSig*
3023 akpk*
3025 (C5*C5)*C6*
3031 t2*t4 -
3032 2.*C1*(C4*C4)*t3*
3033 (3.*t2 -
3034 2.*(C5*C5)*C6*
3038 (1. - 2.*apkk1*C5 +
3042 *t4))) +
3043 18.*C1*C2*C4*C7*Ekaon*t3*t5 -
3044 18.*akpk*C1*C2*C4*C7*Fm1*t3*t5 -
3046 18.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
3050 18.*C2*C5*C6*C7*Ekaon*t4*t5 -
3051 6.*akpk*C2*C5*C6*C7*Fm2*t4*t5 -
3053 6.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
3058 (9.*(C3*C3)*(Enu*(t2*t2) +
3059 2.*akpk*am*C1*(C4*C4)*t2*t3 +
3061 2.*akpk*amSig*C1*(C4*C4)*t2*t3 +
3063 2.*am*apkk1*C1*(C4*C4)*t2*t3 -
3064 2.*amSig*apkk1*C1*(C4*C4)*t2*t3 +
3081 (4.*t2 +
3082 C1*(C4*C4)*
3085 Elep*(-1.*(t2*t2) -
3090 Ekaon*
3091 (3.*(t2*t2) +
3092 4.*C1*(C4*C4)*
3094 (C1*C1)*(C4*C4*C4*C4)*
3101 - 6.*C3*(C5*C5*C5)*C6*
3120 2.*Ekaon*
3121 (apkk1 +
3122 Ekaon*
3128 amLam*
3129 (t2 +
3131 )) +
3137 (2.*t2 -
3138 1.*C1*(C4*C4)*
3139 (2.*apkk1 +
3142 amLam*
3144 )*t3)))*t4 +
3145 (C5*C5*C5*C5*C5*C5)*(C6*C6)*
3154 2.*am*
3163 + 2.*Ekaon*
3164 (apkk1 +
3165 2.*Ekaon*
3167 ))*(t4*t4)) -
3169 (3.*C1*C4*
3171 + C5*C6*
3174 (-9.*(C1*C1)*(C4*C4)*Fm1*
3179 amSig*
3180 (-4.*Ekaon +
3181 2.*(akpk - 1.*apkk1)*Fm1 +
3183 + am*
3185 2.*(akpk - 1.*apkk1)*Fm1 +
3186 4.*Ekaon*
3188 2.*amSig*
3191 C5*C6*Fm2*t4*
3192 (-1.*C5*C6*
3198 amLam*
3199 (2.*
3201 1.*apkk1*Fm2) +
3203 ) +
3204 am*
3206 2.*amLam*
3209 2.*
3210 ((akpk - 1.*apkk1)*Fm2 +
3211 2.*Ekaon*
3213 - 12.*C2*
3215 - 6.*C1*C4*t3*
3216 (C5*C6*
3218 1.*amSig*
3222 (-1.*Fm2 +
3223 Fm1*
3227 1.*amLam*
3230 amSig*
3231 (Fm2 +
3233 + am*
3234 (amSig*
3235 (-1.*Fm2 +
3236 Fm1*
3239 amLam*
3240 (-1.*Fm2 +
3241 Fm1*
3244 2.*
3246 Ekaon*
3247 (Fm2 +
3249 ))))*t4 +
3250 6.*C2*Fm1*
3252 + 2.*Ekaon*
3255 (6.*C2*t1 +
3256 (akpk - 1.*apkk1)*
3261 1.*apkk1*
3263 2.*Ekaon*
3267 ) -
3268 6.*am*
3270 C5*C6*Fm2*t4*
3272 2.*C2*C7*Ekaon*t5) +
3273 C1*C4*t3*
3277 + 2.*(Enu*(-36.*apkk1*(t1*t1) - 72.*apkk1*C2*(t1*t1) -
3324 (t3*t3) -
3328 (t3*t3) -
3332 (t3*t3) -
3339 (t3*t3) -
3341 (t3*t3) -
3345 (t3*t3) -
3347 (t3*t3) -
3349 (t3*t3) -
3351 (t3*t3) +
3390 t3*t4 +
3394 C6*t3*t4 -
3396 t3*t4 -
3401 t3*t4 +
3404 C6*t3*t4 +
3407 t3*t4 -
3409 t3*t4 -
3411 t3*t4 -
3417 Ekaon*t3*t4 +
3420 Ekaon*t3*t4 +
3422 t3*t4 +
3426 Ekaon*t3*t4 -
3429 Ekaon*t3*t4 +
3432 Ekaon*t3*t4 -
3434 Ekaon*t3*t4 -
3438 t3*t4 +
3444 t3*t4 -
3447 t3*t4 -
3450 t3*t4 +
3452 t3*t4 +
3454 Elep*t3*t4 -
3461 Fm1*t3*t4 -
3464 Fm1*t3*t4 -
3466 t3*t4 -
3468 t3*t4 +
3470 t3*t4 -
3472 Fm1*t3*t4 +
3474 Fm1*t3*t4 -
3476 Fm1*t3*t4 +
3478 Fm1*t3*t4 -
3480 Fm1*t3*t4 +
3482 Fm1*t3*t4 +
3484 Fm1*t3*t4 -
3486 Fm1*t3*t4 +
3492 t3*t4 +
3494 t3*t4 +
3496 Fm1*t3*t4 -
3498 Fm1*t3*t4 +
3500 Fm1*t3*t4 -
3502 Fm1*t3*t4 -
3509 Fm2*t3*t4 +
3511 t3*t4 -
3513 t3*t4 +
3515 t3*t4 -
3517 t3*t4 -
3519 Fm2*t3*t4 +
3521 Fm2*t3*t4 +
3523 t3*t4 -
3525 t3*t4 -
3527 Fm2*t3*t4 +
3529 Fm2*t3*t4 -
3531 Fm2*t3*t4 +
3533 Fm2*t3*t4 +
3539 t3*t4 -
3541 t3*t4 +
3543 Fm2*t3*t4 -
3545 Fm2*t3*t4 -
3547 Fm2*t3*t4 +
3549 Fm2*t3*t4 -
3568 (t4*t4) +
3572 (t4*t4) -
3576 (t4*t4) -
3583 (t4*t4) -
3585 (t4*t4) -
3589 (t4*t4) -
3591 (t4*t4) -
3593 (t4*t4) -
3595 (t4*t4) +
3600 36.*am*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 +
3601 36.*amSig*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 +
3602 36.*am*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 +
3603 36.*amSig*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 +
3608 36.*am*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 +
3609 36.*amLam*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 -
3610 12.*am*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 -
3611 12.*amLam*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 +
3619 t4 +
3621 t4 +
3623 t4 +
3625 t4 -
3633 3.*C3*
3634 (6.*t1*
3635 (t2 +
3638 (3.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4) -
3641 (3.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4) +
3643 (-3.*(C1*C1)*(C4*C4*C4*C4)*
3646 1.*(C5*C5)*C6*t2*t4 +
3647 C1*(C4*C4)*t3*
3648 (3.*t2 +
3649 (C5*C5)*C6*
3658 (9.*(C1*C1)*(C4*C4)*
3659 (-1. + 2.*C4 +
3660 (C4*C4)*
3661 (-1. + apkk1*C3 +
3664 (C5*C5)*C6*t4*
3665 (3.*C3*t2 +
3666 C6*
3667 (-9. - 6.*C5 - 1.*(C5*C5) +
3668 (C5*C5*C5)*
3671 3.*C1*C4*t3*
3672 (C5*C6*
3673 (6. + 2.*C5 +
3674 C4*
3675 (-6. - 2.*C5 +
3676 (C5*C5)*
3682 C3*C4*
3683 (3.*t2 +
3684 (C5*C5)*C6*
3686 ))) -
3687 6.*t1*
3688 (-1.*apkk1*(2. + C2)*
3692 (3.*C1*C4*
3693 (2.*(-1. + C4) +
3694 C2*
3697 3.*C3*
3698 (t2 -
3699 1.*C1*(C4*C4)*
3701 ) +
3702 C5*C6*
3703 (-6. - 2.*C5 +
3704 (C5*C5)*
3706 C2*
3709 1.*am*
3710 (3.*C1*C4*
3711 (apkk1*(2. + C2)*Fm1 +
3712 2.*Ekaon*
3713 (-2. + 2.*C2*(-1. + C4) + 2.*C4 +
3715 6.*C3*
3717 1.*Enu*t2 +
3718 amSig*apkk1*C1*(C4*C4)*t3 +
3719 apkk1*C1*(C4*C4)*Ekaon*t3) +
3720 C5*C6*
3721 (-2.*
3722 (6. + 2.*C5 - 1.*apkk1*(C5*C5) +
3723 2.*C2*(3. + C5))*Ekaon +
3724 apkk1*(2. + C2)*Fm2 +
3726 *t4) -
3728 (apkk1*
3729 (-3.*C1*C3*(C4*C4)*t3 + (C5*C5*C5)*C6*t4)
3730 + Ekaon*
3731 (3.*C1*C4*
3733 2.*C2*Fm1)*t3 +
3734 C5*C6*
3736 2.*C2*Fm2)*t4))) -
3738 (-1.*apkk1*
3739 (3.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4)*
3740 (3.*C1*C3*(C4*C4)*t3 - 1.*(C5*C5*C5)*C6*t4)
3741 + 2.*Ekaon*
3742 (9.*(C1*C1)*(C4*C4)*
3745 (t3*t3) -
3746 3.*C1*C4*C5*C6*
3749 2.*C4*
3752 (C5*C5)*(C6*C6)*
3755 *(t4*t4))) +
3757 (-18.*(C1*C1)*(C4*C4)*Ekaon*
3758 (-1. + (-1. + apkk1*C3)*(C4*C4) -
3761 2.*(C5*C5)*C6*t4*
3763 t2 +
3764 C6*Ekaon*
3765 (9. + (C5*C5) - 1.*apkk1*(C5*C5*C5) +
3768 3.*C1*C4*t3*
3769 (C5*C6*
3770 (2.*
3771 (2.*(3. + C5) +
3772 C4*(-6. - 2.*C5 + apkk1*(C5*C5)))*
3773 Ekaon +
3774 amLam*
3775 (apkk1*C4*(C5*C5) +
3776 2.*Ekaon*
3778 amSig*
3779 (-1.*apkk1*C4*(C5*C5) +
3780 2.*Ekaon*
3782 C3*C4*
3784 6.*Enu*t2 -
3785 1.*amLam*apkk1*(C5*C5)*C6*t4 +
3786 amSig*apkk1*(C5*C5)*C6*t4 +
3787 2.*apkk1*(C5*C5)*C6*Ekaon*t4))) -
3788 9.*apkk1*C2*C7*
3790 am*apkk1*
3792 (t3*t3) -
3796 t4 +
3798 ((3. + 2.*C5)*Fm1 +
3799 (-1. + 2.*C4)*Fm2)*t3*t4 +
3801 9.*C2*C7*
3809 (t3*t3) +
3814 t4 +
3816 C6*t3*t4 +
3818 t3*t4 +
3822 t4 -
3824 t4 -
3828 t4 -
3835 (9.*(C1*C1)*(C4*C4)*
3838 3.*C1*C4*C5*C6*
3843 (C5*C5)*(C6*C6)*
3847 (-9.*(C1*C1)*(C4*C4)*
3848 (-1. +
3849 (C4*C4)*
3850 (-1. + 3.*apkk1*C3 -
3855 C4*
3858 (t3*t3) +
3859 3.*C1*C4*C5*C6*
3860 (C4*
3861 (-6. +
3862 (C5*C5)*
3863 (3.*apkk1 -
3866 C5*
3867 (-2. + 3.*apkk1*C3 -
3875 2.*
3877 1.*C5*
3884 (C5*C5)*(C6*C6)*
3885 (9. + (C5*C5) +
3886 (C5*C5*C5)*
3887 (-3.*apkk1 +
3891 C5*
3894 (t4*t4)) +
3895 18.*amSig*C1*C2*C4*C7*t3*t5 -
3896 36.*amSig*C1*C2*(C4*C4)*C7*t3*t5 +
3897 54.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
3898 18.*amLam*C2*C5*C6*C7*t4*t5 +
3899 12.*amLam*C2*(C5*C5)*C6*C7*t4*t5 +
3900 18.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
3901 am*
3902 (18.*(C1*C1)*(C4*C4)*
3903 (2.*
3904 (1. - 2.*C4 +
3905 (1. - 1.*apkk1*C3)*(C4*C4))*Ekaon +
3906 amSig*
3907 (1. + (1. - 2.*apkk1*C3)*(C4*C4) +
3909 ) +
3913 3.*C1*C4*t3*
3914 (6.*amSig*C5*C6*t4 -
3915 6.*amSig*C4*C5*C6*t4 +
3916 2.*amSig*(C5*C5)*C6*t4 -
3917 2.*amSig*C4*(C5*C5)*C6*t4 +
3918 amSig*apkk1*C4*(C5*C5*C5)*C6*t4 +
3919 24.*C5*C6*Ekaon*t4 -
3920 24.*C4*C5*C6*Ekaon*t4 +
3921 8.*(C5*C5)*C6*Ekaon*t4 -
3922 8.*C4*(C5*C5)*C6*Ekaon*t4 +
3923 4.*apkk1*C4*(C5*C5*C5)*C6*Ekaon*t4 -
3928 amLam*C5*C6*
3929 (C4*
3930 (-6. - 2.*C5 +
3931 (C5*C5)*
3935 2.*
3939 ))*t4 +
3940 C3*C4*
3942 amLam*apkk1*(C5*C5)*C6*t4 +
3943 3.*amSig*apkk1*(C5*C5)*C6*t4 -
3946 4.*Ekaon*
3947 (3.*t2 + apkk1*(C5*C5)*C6*t4)) +
3948 6.*C2*C7*t5 - 12.*C2*C4*C7*t5) +
3949 2.*C5*C6*t4*
3950 (-3.*C3*C5*
3952 18.*C5*C6*Ekaon*t4 +
3953 12.*(C5*C5)*C6*Ekaon*t4 +
3954 2.*(C5*C5*C5)*C6*Ekaon*t4 -
3955 2.*apkk1*(C5*C5*C5*C5)*C6*Ekaon*t4 +
3959 amLam*C5*C6*
3960 (9. + (C5*C5) - 2.*apkk1*(C5*C5*C5) +
3963 9.*C2*C7*t5 + 6.*C2*C5*C7*t5)))))
3964 + 2.*(akpk*akpk)*
3967 (t3*t3) -
3968 3.*C1*C4*C5*C6*
3971 + 3.*Elep*
3974 amSig*C1*C4*t3*
3976 Fm1*
3977 (-6.*(-1. + C2)*t1 -
3979 6.*C2*(-1. + C4)*C7*t5) +
3980 amLam*C5*C6*t4*
3982 Fm2*
3983 (-2.*(-1. + C2)*t1 +
3985 2.*C2*(3. + C5)*C7*t5)) +
3988 3.*amSig*C1*C4*t3*
3991 C3*C4*
3992 (-6.*t1 - 3.*t2 -
3999 6.*C2*C7*Fm1*t5) +
4000 amLam*C5*C6*t4*
4003 C5*
4006 t4) +
4008 ) + am*
4012 3.*C1*C4*t3*
4015 C3*C4*
4016 (-6.*t1 - 3.*t2 +
4018 6.*C2*C7*Fm1*t5) +
4019 C5*C6*t4*
4020 (6.*(C5*C5)*t1 + 3.*C3*C5*t2 -
4023 6.*C2*C7*Fm2*t5)) +
4024 6.*Elep*
4025 (-1.*C5*C6*t4*
4026 (C2*(3. + C5)*C7*t5 +
4027 Fm2*
4028 ((-1. + C2)*t1 - 2.*C2*C7*Ekaon*t5))
4029 + C1*C4*t3*
4033 t4 - 3.*C2*C7*t5 + 3.*C2*C4*C7*t5 +
4034 Fm1*
4035 (-3.*(-1. + C2)*t1 +
4036 Ekaon*
4037 (-1.*amLam*(-3. + C5)*C5*C6*t4 +
4038 amSig*(-3. + C5)*C5*C6*t4 +
4039 6.*C2*C7*t5)))))) +
4045 72.*am*apkk1*C1*C2*C4*t1*t3 +
4046 72.*amSig*apkk1*C1*C2*C4*t1*t3 -
4049 72.*am*apkk1*C1*(C4*C4)*t1*t3 -
4050 72.*amSig*apkk1*C1*(C4*C4)*t1*t3 -
4098 (t3*t3) +
4100 (t3*t3) +
4102 (t3*t3) +
4115 (t3*t3) -
4117 (t3*t3) +
4119 Ekaon*(t3*t3) +
4121 (t3*t3) -
4124 (t3*t3) +
4126 (t3*t3) -
4128 (t3*t3) +
4130 (t3*t3) +
4132 Enu*(t3*t3) -
4134 (t3*t3) +
4136 (t3*t3) -
4139 (t3*t3) -
4141 (t3*t3) -
4144 (t3*t3) -
4146 (t3*t3) -
4148 (t3*t3) +
4150 (t3*t3) -
4152 (t3*t3) +
4154 (t3*t3) +
4156 (t3*t3) +
4158 (t3*t3) -
4160 (t3*t3) -
4162 (t3*t3) +
4164 (t3*t3) +
4166 (t3*t3) -
4170 (t3*t3) -
4176 (t3*t3) -
4186 72.*am*apkk1*C2*C5*C6*t1*t4 +
4187 72.*amLam*apkk1*C2*C5*C6*t1*t4 +
4190 24.*am*apkk1*(C5*C5)*C6*t1*t4 +
4191 24.*amLam*apkk1*(C5*C5)*C6*t1*t4 +
4236 t3*t4 -
4242 C6*t3*t4 +
4245 t3*t4 +
4248 C6*t3*t4 -
4251 (C5*C5)*C6*t3*t4 -
4253 t3*t4 -
4255 C6*t3*t4 -
4257 (C5*C5)*C6*t3*t4 -
4259 (C5*C5)*C6*t3*t4 -
4263 C6*t3*t4 -
4265 t3*t4 -
4267 C6*t3*t4 -
4269 C6*t3*t4 +
4274 t3*t4 -
4277 t3*t4 +
4280 t3*t4 -
4283 C6*Ekaon*t3*t4 -
4286 Ekaon*t3*t4 +
4289 C6*Ekaon*t3*t4 +
4291 Ekaon*t3*t4 -
4293 (C5*C5)*C6*Ekaon*t3*t4 +
4296 Ekaon*t3*t4 +
4298 Ekaon*t3*t4 -
4300 C6*Ekaon*t3*t4 -
4302 Enu*t3*t4 +
4305 (C5*C5)*C6*Enu*t3*t4 -
4307 C6*Enu*t3*t4 -
4309 (C5*C5)*C6*Enu*t3*t4 -
4311 t3*t4 +
4314 (C5*C5*C5)*C6*Enu*t3*t4 -
4316 Enu*t3*t4 -
4318 C6*Enu*t3*t4 +
4320 Enu*t3*t4 -
4326 t3*t4 -
4328 Enu*t3*t4 -
4334 Enu*t3*t4 +
4338 Fm1*t3*t4 +
4341 C6*Fm1*t3*t4 +
4345 C6*Fm1*t3*t4 +
4347 t3*t4 +
4349 Fm1*t3*t4 -
4352 Fm1*t3*t4 +
4354 t3*t4 -
4356 t3*t4 -
4358 t3*t4 +
4360 Fm1*t3*t4 -
4365 Fm1*t3*t4 +
4367 Fm1*t3*t4 +
4371 t3*t4 -
4373 Fm1*t3*t4 -
4375 Fm1*t3*t4 +
4378 Fm1*t3*t4 -
4380 t3*t4 +
4382 t3*t4 -
4384 t3*t4 +
4386 t3*t4 +
4388 Fm1*t3*t4 -
4390 Fm1*t3*t4 +
4394 Fm2*t3*t4 -
4397 C6*Fm2*t3*t4 -
4401 C6*Fm2*t3*t4 -
4403 t3*t4 -
4405 Fm2*t3*t4 -
4409 Fm2*t3*t4 -
4411 t3*t4 +
4413 t3*t4 +
4415 Fm2*t3*t4 +
4420 Fm2*t3*t4 -
4422 Fm2*t3*t4 -
4426 t3*t4 +
4428 Fm2*t3*t4 +
4430 Fm2*t3*t4 -
4433 Fm2*t3*t4 -
4435 t3*t4 +
4437 t3*t4 +
4439 t3*t4 -
4441 t3*t4 +
4443 Fm2*t3*t4 -
4445 Fm2*t3*t4 -
4447 Fm2*t3*t4 +
4456 Fm2*t3*t4 -
4462 Fm2*t3*t4 -
4477 (t4*t4) +
4491 (t4*t4) +
4493 (t4*t4) +
4495 (t4*t4) -
4498 (t4*t4) +
4500 (t4*t4) -
4503 (t4*t4) +
4505 (t4*t4) -
4507 (t4*t4) +
4509 (t4*t4) -
4512 (t4*t4) +
4514 (t4*t4) +
4517 (t4*t4) +
4519 (t4*t4) -
4522 (t4*t4) -
4524 (t4*t4) +
4526 (t4*t4) -
4529 (t4*t4) +
4531 (t4*t4) +
4533 (t4*t4) +
4535 (t4*t4) -
4537 (t4*t4) -
4539 (t4*t4) +
4541 (t4*t4) -
4547 (t4*t4) -
4555 (amLam*C5*C6*
4557 Fm2*
4558 (2.*(1. + C2)*t1 +
4560 amSig*C1*C4*t3*
4562 Fm1*
4563 (6.*(1. + C2)*t1 -
4565 2.*am*
4566 ((1. + C2)*C5*C6*Fm2*t1*t4 +
4567 C1*C4*t3*
4570 Fm1*
4571 (3.*(1. + C2)*t1 -
4573 C6*Ekaon*t4)))) +
4575 144.*apkk1*(C2*C2)*C7*t1*t5 -
4592 t3*t5 +
4596 t3*t5 -
4613 t4*t5 +
4618 144.*apkk1*(C2*C2)*(C7*C7)*Ekaon*(t5*t5) +
4620 144.*apkk1*(C2*C2)*(C7*C7)*Enu*(t5*t5) -
4621 1.*Elep*(36.*((-1. + C2)*(-1. + C2))*(t1*t1) -
4631 (t3*t3) +
4641 t4 +
4643 t4 -
4645 (C5*C5)*C6*t3*t4 +
4648 (C5*C5*C5)*C6*t3*t4 +
4650 t4 -
4652 t4 +
4654 t4 -
4656 t4 +
4658 t3*t4 -
4660 t3*t4 +
4662 t3*t4 -
4664 t3*t4 -
4666 t4 +
4668 t4 -
4670 t4 +
4672 t4 +
4674 t3*t4 -
4676 t3*t4 +
4678 t3*t4 -
4689 (9.*(C1*C1)*(C4*C4)*
4690 (-1. + 2.*C4 +
4692 3.*C1*C4*C5*
4693 (2.*(3. + C5) +
4694 C4*
4697 (C5*C5)*
4698 (-9. - 6.*C5 - 1.*(C5*C5) +
4700 6.*t1*
4704 3.*C3*
4705 (2.*t2 +
4710 1.*am*
4711 (3.*C1*(2. + C2)*C4*Fm1*t3 +
4712 C5*C6*
4714 t4)) +
4715 2.*
4717 (3.*C1*C4*(1. + C4)*t3 -
4718 1.*(-3. + C5)*C5*C6*t4) -
4719 1.*apkk1*(1. + C2)*
4723 (3.*C1*C4*
4726 1.*C5*C6*
4729 + am*
4730 (-3.*C1*C4*
4731 (2.*(-1. + C2)*(1. + C4)*Ekaon +
4732 apkk1*(1. + C2)*Fm1)*t3 +
4733 3.*amSig*C1*C4*
4736 1.*C5*C6*
4737 (-2.*(-1. + C2)*(-3. + C5)*Ekaon +
4738 apkk1*(1. + C2)*Fm2 +
4739 amLam*
4742 ) -
4745 (3.*(C1*C1)*(C4*C4)*
4746 (3. - 2.*C4 + 3.*(C4*C4))*(t3*t3) +
4747 2.*C1*C4*C5*
4748 (9. + C5 - 3.*C4*(1. + C5))*C6*t3*t4
4749 + (C5*C5)*(9. + 2.*C5 + (C5*C5))*(C6*C6)*
4750 (t4*t4)) -
4751 4.*Ekaon*
4752 (9.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
4753 (t3*t3) -
4754 3.*C1*C4*C5*C6*
4755 (2.*(1. + C4)*(-3. + C5)*Ekaon +
4756 amSig*
4760 amLam*
4764 ((-3. + C5)*(-3. + C5))*(C5*C5)*(C6*C6)*Ekaon*(t4*t4)
4765 ) +
4767 (-9.*(C1*C1)*C3*(C4*C4*C4*C4)*
4769 (C5*C5)*C6*t4*
4770 (6.*C3*t2 +
4771 (C5*C5*C5)*C6*
4773 3.*C1*(C4*C4)*t3*
4774 ((C5*C5*C5)*C6*
4776 2.*Ekaon*
4778 + C3*
4779 (6.*t2 +
4780 (C5*C5)*C6*
4782 2.*Ekaon*
4784 )) -
4787 36.*amSig*apkk1*C1*C2*C4*C7*t3*t5 +
4788 36.*amSig*apkk1*C1*C2*(C4*C4)*C7*t3*t5 +
4793 36.*amLam*apkk1*C2*C5*C6*C7*t4*t5 -
4794 12.*amLam*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
4798 144.*apkk1*(C2*C2)*(C7*C7)*(t5*t5) -
4799 1.*am*
4802 (t3*t3) -
4804 (t3*t3) -
4805 3.*amSig*C1*C4*t3*
4806 (C5*C6*
4808 2.*Ekaon*
4811 C3*C4*
4812 (6.*t2 -
4815 (12.*C1*(C4*C4)*t3 - 1.*(C5*C5)*C6*t4))
4816 ) +
4817 3.*C1*C4*t3*
4824 C3*C4*
4826 4.*Ekaon*
4828 6.*C2*C7*t5 +
4830 2.*C5*C6*t4*
4831 (-3.*amLam*C3*C5*t2 -
4832 6.*C3*C5*Ekaon*t2 +
4837 )) +
4838 2.*
4841 (t3*t3) +
4843 (t3*t3) +
4845 t4 -
4847 t3*t4 +
4849 t3*t4 -
4851 Ekaon*t3*t4 +
4853 Fm1*t3*t4 -
4855 Fm1*t3*t4 +
4857 Fm1*t3*t4 -
4859 Fm1*t3*t4 -
4861 t3*t4 +
4863 t3*t4 +
4865 Fm2*t3*t4 -
4867 Fm2*t3*t4 +
4869 (t4*t4) +
4871 (t4*t4) +
4873 18.*apkk1*C1*C2*C4*C7*t3*t5 -
4874 18.*apkk1*C1*C2*(C4*C4)*C7*t3*t5 +
4876 t5 - 18.*apkk1*C2*C5*C6*C7*t4*t5 +
4877 6.*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
4879 t5 +
4881 (18.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
4882 (t3*t3) +
4883 3.*amSig*C1*C4*t3*
4884 (6.*C1*C4*(1. + (C4*C4))*t3 +
4885 C5*C6*
4889 2.*C5*C6*t4*
4890 (amLam*C5*(9. + (C5*C5))*C6*t4 +
4891 9.*C5*C6*Ekaon*t4 -
4892 6.*(C5*C5)*C6*Ekaon*t4 +
4893 (C5*C5*C5)*C6*Ekaon*t4 + 9.*C2*C7*t5) -
4894 3.*C1*C4*t3*
4895 (amLam*C5*C6*
4899 2.*
4900 (-6.*(1. + C4)*C5*C6*Ekaon*t4 +
4901 2.*(1. + C4)*(C5*C5)*C6*Ekaon*t4 -
4902 3.*C2*C7*t5))))))))));
4903
4904 if (sol <= 0.0) {
4906 << "Matrix element is 0 or negative! This is bad and really should throw an exception!";
4907 }
4908
4909 return sol;
4910}
double amLam
Definition AlamSimoAtharVacasSKPXSec2014.h:63
References alepvec, am, amk, aml, amLam, ampi, amSig, angkq, aq0, aqvec, d, Ekaon, Elep, Enu, f, Fm1, Fm2, fpi, g, LOG, pERROR, pkvec, and Vus.
Referenced by XSec().
◆ Amatrix_PP()
|
private |
Definition at line 4914 of file AlamSimoAtharVacasSKPXSec2014.cxx.
4915{
4916 double sol = 0.;
4917
4918 double akk1=0., zdotq=0., qdotpk=0., akcrosk1=0., qcrospk=0.;
4919 double zdotpk=0., azpk=0., aqkaon=0., akpk=0., apkk1=0.;
4920 double C1=0., C2=0., C3=0., C4=0., C5=0., C6=0., C7=0., C8=0., C9=0.;
4921 double aq2=0., gform=0., con=0., t1=0., t2=0., t3=0., t4=0., t5=0., t6=0.;
4922
4931 akpk=azpk + aqkaon/2.0;
4932 apkk1=azpk - aqkaon/2.0;
4943 gform=1.0/pow(1.0-aq2/(1.1*1.1),4);
4944
4946
4947 t1=1.0;
4948 t2=1.0; // !Full Term
4949 t3=1.0;
4950 t4=1.0;
4951 t5=1.0;
4952 t6=1.0;
4953
4954 sol = con*gform*(0.4444444444444444*am*
4972 144.*akpk*amSig*C1*C4*t1*t3 +
4973 144.*amSig*apkk1*C1*C4*t1*t3 -
4974 144.*akpk*amSig*C1*(C4*C4)*t1*t3 +
4975 144.*amSig*apkk1*C1*(C4*C4)*t1*t3 +
5011 (C4*C4)*t2*t3 -
5021 Elep*t2*t3 +
5034 (C4*C4*C4*C4)*(t3*t3) -
5037 (C4*C4*C4*C4)*(t3*t3) +
5040 (C4*C4*C4*C4)*(t3*t3) +
5042 (C4*C4*C4*C4)*(t3*t3) -
5045 (C3*C3)*(C4*C4*C4*C4)*(t3*t3) +
5055 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
5058 Ekaon*(t3*t3) +
5061 (C4*C4*C4*C4)*Ekaon*(t3*t3) +
5063 Ekaon*(t3*t3) -
5065 Ekaon*(t3*t3) -
5067 Ekaon*(t3*t3) -
5077 (C3*C3)*(C4*C4*C4*C4)*Elep*(t3*t3) +
5080 (C4*C4*C4*C4)*Elep*(t3*t3) +
5082 (C4*C4*C4*C4)*Elep*(t3*t3) +
5092 (C4*C4*C4*C4)*Enu*(t3*t3) +
5095 (C4*C4*C4*C4)*Enu*(t3*t3) -
5098 (C4*C4*C4*C4)*Enu*(t3*t3) -
5100 (C4*C4*C4*C4)*Enu*(t3*t3) +
5107 Fm1*(t3*t3) +
5113 Fm1*(t3*t3) -
5149 144.*akpk*amLam*C5*C6*t1*t4 -
5150 144.*amLam*apkk1*C5*C6*t1*t4 -
5151 48.*akpk*amLam*(C5*C5)*C6*t1*t4 +
5152 48.*amLam*apkk1*(C5*C5)*C6*t1*t4 +
5189 (C5*C5)*C6*t2*t4 +
5198 Elep*t2*t4 +
5221 C6*t3*t4 -
5224 (C5*C5)*C6*t3*t4 -
5227 (C4*C4)*(C5*C5)*C6*t3*t4 +
5229 C6*t3*t4 +
5232 (C3*C3)*(C4*C4)*(C5*C5)*C6*t3*t4 +
5234 C6*t3*t4 +
5237 (C3*C3)*(C4*C4)*(C5*C5)*C6*t3*t4 -
5239 (C5*C5)*C6*t3*t4 +
5241 (C4*C4)*(C5*C5)*C6*t3*t4 -
5243 (C5*C5)*C6*t3*t4 -
5245 (C5*C5)*C6*t3*t4 -
5247 (C5*C5)*C6*t3*t4 -
5253 Ekaon*t3*t4 -
5256 Ekaon*t3*t4 -
5259 Ekaon*t3*t4 -
5261 Ekaon*t3*t4 -
5263 C6*Ekaon*t3*t4 +
5266 (C4*C4)*(C5*C5)*C6*Ekaon*t3*t4 +
5268 C6*Ekaon*t3*t4 +
5270 (C5*C5)*C6*Ekaon*t3*t4 +
5272 (C5*C5)*C6*Ekaon*t3*t4 -
5274 (C5*C5)*C6*Ekaon*t3*t4 -
5276 (C5*C5)*C6*Ekaon*t3*t4 +
5286 (C5*C5)*C6*Elep*t3*t4 -
5289 (C3*C3)*(C4*C4)*(C5*C5)*C6*Elep*t3*t4 +
5291 (C5*C5)*C6*Elep*t3*t4 -
5293 Elep*t3*t4 +
5295 (C4*C4)*(C5*C5)*C6*Elep*t3*t4 -
5306 (C4*C4)*(C5*C5)*C6*Enu*t3*t4 -
5308 Enu*t3*t4 +
5311 (C3*C3)*(C4*C4)*(C5*C5)*C6*Enu*t3*t4 -
5313 (C5*C5)*C6*Enu*t3*t4 +
5315 Enu*t3*t4 -
5317 (C4*C4)*(C5*C5)*C6*Enu*t3*t4 -
5322 Fm1*t3*t4 +
5329 C6*Fm1*t3*t4 -
5332 Fm1*t3*t4 -
5335 Fm1*t3*t4 +
5338 Fm1*t3*t4 +
5347 Fm1*t3*t4 +
5350 Fm1*t3*t4 -
5355 Fm2*t3*t4 +
5362 C6*Fm2*t3*t4 -
5365 Fm2*t3*t4 -
5368 Fm2*t3*t4 +
5371 Fm2*t3*t4 -
5380 Fm2*t3*t4 +
5383 Fm2*t3*t4 +
5389 Fm2*t3*t4 -
5392 Fm2*t3*t4 +
5394 Fm2*t3*t4 +
5396 Fm2*t3*t4 +
5399 Fm2*t3*t4 -
5410 Fm2*t3*t4 +
5432 (C5*C5*C5*C5)*(C6*C6)*(t4*t4) +
5435 (C6*C6)*(t4*t4) +
5437 (C6*C6)*(t4*t4) -
5440 (C5*C5*C5*C5)*(C6*C6)*(t4*t4) +
5450 (C6*C6)*Ekaon*(t4*t4) +
5453 Ekaon*(t4*t4) +
5456 (C6*C6)*Ekaon*(t4*t4) +
5458 Ekaon*(t4*t4) -
5461 (C6*C6)*Ekaon*(t4*t4) -
5471 (C5*C5*C5*C5)*(C6*C6)*Elep*(t4*t4) +
5474 (C6*C6)*Elep*(t4*t4) +
5476 (C6*C6)*Elep*(t4*t4) +
5486 Enu*(t4*t4) +
5489 Enu*(t4*t4) -
5492 (C6*C6)*Enu*(t4*t4) -
5494 (C6*C6)*Enu*(t4*t4) +
5501 Fm2*(t4*t4) -
5511 Fm2*(t4*t4) -
5542 (1. +
5543 (C4*C4)*
5550 6.*C1*C4*C5*C6*
5555 C4*
5556 (C5*
5561 (C5*C5)*(C6*C6)*
5562 (9. +
5563 (C5*C5)*
5571 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
5573 (9.*(C1*C1)*(C4*C4)*
5574 (2.*C3*(C4*C4)*Ekaon +
5578 (t3*t3) +
5579 6.*C1*C4*C5*C6*
5580 (2.*C3*C4*C5*Ekaon +
5584 t4 +
5585 (C5*C5)*(C6*C6)*
5586 (2.*C3*(C5*C5)*Ekaon +
5590 (t4*t4)) -
5592 (apkk1*
5595 akpk*
5598 2.*Ekaon*
5599 (3.*C1*C4*
5602 C5*C6*
5618 12.*apkk1*C4*Fm1 +
5623 4.*amSig*
5624 (1. +
5625 (C4*C4)*
5632 - 1.*Ekaon*
5633 (4. +
5634 (C4*C4)*
5643 2.*Ekaon*
5644 (2.*(akpk*akpk)*
5646 2.*(apkk1*apkk1)*
5648 apkk1*
5649 (2. + 4.*C4 +
5650 (C4*C4)*
5656 1.*akpk*
5657 (2. - 4.*C4 +
5658 (C4*C4)*
5661 (-1. + apkk1*C3 +
5665 2.*Ekaon*
5668 2.*
5669 (((1.+C4)*(1.+C4))*Elep -
5670 1.*((-1.+C4)*(-1.+C4))*Enu) +
5671 amSig*
5672 (-2. +
5681 ))*(t3*t3) +
5682 6.*C1*C4*t3*
5715 (amLam*
5716 (6. +
5720 amSig*
5721 (6. +
5725 2.*
5726 (-1.*Enu*
5729 Elep*
5732 Ekaon*
5735 t1 + 12.*akpk*C5*C6*t4 +
5737 12.*apkk1*C5*C6*t4 +
5738 4.*akpk*(C5*C5)*C6*t4 +
5739 4.*apkk1*(C5*C5)*C6*t4 -
5750 4.*Elep*
5752 Fm1*
5753 (6.*t1 -
5762 t4 +
5764 t4 +
5771 Elep*t4 -
5779 t4 -
5781 Enu*t4 +
5795 6.*akpk*
5802 t4 -
5803 2.*C5*C6*
5805 2.*C5*Enu +
5808 t4 +
5810 (-1.*(C5*C5)*C6*
5813 (3.*t2 +
5814 (-1.*akpk + apkk1)*(C5*C5)*C6*t4)))
5815 + Ekaon*
5820 12.*apkk1*C5*C6*t4 +
5822 4.*apkk1*(C5*C5)*C6*t4 +
5825 t4 -
5829 t4 -
5838 2.*akpk*
5840 2.*C5*(3. + C5)*C6*t4 +
5842 (3.*C3*t2 - 2.*(C5*C5)*C6*t4 +
5843 2.*apkk1*C3*(C5*C5)*C6*t4))))) +
5844 C5*C6*t4*
5846 Fm2 +
5847 Ekaon*
5850 t1 +
5851 (C5*C5*C5)*C6*
5856 Ekaon*
5860 4.*amLam*
5868 2.*Ekaon*
5871 apkk1*
5875 2.*Ekaon*
5876 (amLam*
5879 2.*Enu +
5883 1.*akpk*
5886 (-1. + apkk1*C3 +
5888 4.*(C5*C5)*C6*
5895 )*t4 +
5910 12.*akpk*
5914 (2.*C6*
5918 amLam*
5923 t4)) -
5924 1.*Ekaon*
5938 6.*akpk*
5939 (6.*C6*t4 +
5941 ))))) -
6091 (t3*t3) -
6094 (t3*t3) +
6096 (t3*t3) +
6098 (t3*t3) -
6101 (t3*t3) -
6103 (t3*t3) +
6105 (t3*t3) -
6107 (t3*t3) +
6109 (t3*t3) -
6111 (t3*t3) -
6113 (t3*t3) +
6116 (t3*t3) +
6118 (t3*t3) -
6120 (t3*t3) -
6123 (t3*t3) -
6125 (t3*t3) -
6127 (t3*t3) -
6130 (t3*t3) +
6141 (t3*t3) +
6143 (t3*t3) +
6163 t3*t4 -
6172 t3*t4 +
6178 t3*t4 +
6192 Fm1*t3*t4 +
6195 Fm1*t3*t4 +
6198 t3*t4 +
6201 t3*t4 -
6203 t3*t4 -
6212 t3*t4 -
6215 Fm1*t3*t4 -
6218 Fm1*t3*t4 -
6220 t3*t4 -
6222 t3*t4 -
6225 Fm1*t3*t4 -
6227 t3*t4 -
6229 t3*t4 +
6231 t3*t4 -
6233 t3*t4 -
6235 t3*t4 -
6237 t3*t4 +
6240 t3*t4 +
6243 Fm1*t3*t4 -
6245 t3*t4 -
6247 t3*t4 -
6249 t3*t4 +
6251 t3*t4 -
6254 Fm2*t3*t4 -
6257 Fm2*t3*t4 +
6261 t3*t4 -
6263 t3*t4 +
6272 t3*t4 -
6275 Fm2*t3*t4 +
6278 Fm2*t3*t4 +
6280 t3*t4 +
6282 t3*t4 -
6285 Fm2*t3*t4 -
6288 Fm2*t3*t4 +
6290 t3*t4 +
6292 t3*t4 -
6294 t3*t4 +
6296 t3*t4 +
6299 t3*t4 +
6302 Fm2*t3*t4 +
6304 t3*t4 +
6306 t3*t4 -
6308 t3*t4 +
6310 t3*t4 +
6313 t3*t4 -
6315 Fm2*t3*t4 -
6317 Fm2*t3*t4 -
6319 Fm2*t3*t4 -
6326 t3*t4 -
6328 t3*t4 -
6330 t3*t4 -
6335 Fm2*t3*t4 +
6364 (t4*t4) +
6366 (t4*t4) +
6368 (t4*t4) +
6371 (t4*t4) +
6374 (t4*t4) -
6376 (t4*t4) +
6378 (t4*t4) +
6380 (t4*t4) +
6382 (t4*t4) +
6385 (t4*t4) +
6387 (t4*t4) +
6389 (t4*t4) -
6392 (t4*t4) -
6394 (t4*t4) +
6396 (t4*t4) -
6400 (t4*t4) +
6402 (t4*t4) +
6407 (t4*t4) +
6409 (t4*t4) +
6411 (t4*t4) +
6413 (t4*t4) +
6416 (t4*t4) +
6461 2.*akpk*
6462 (3.*C1*C4*
6467 t3 +
6468 C5*C6*
6473 - 2.*apkk1*
6474 (3.*C1*C4*
6478 t3 +
6479 C5*C6*
6484 (18.*apkk1*C1*(C4*C4)*t2*t3 +
6489 6.*apkk1*(C5*C5)*C6*t2*t4 +
6494 t3*t4 -
6497 akpk*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
6498 (-6.*t2 +
6501 (9.*(C1*C1)*(C4*C4*C4*C4)*
6505 (t3*t3) +
6506 (C5*C5)*C6*t4*
6507 (-12.*t2 -
6508 1.*(C5*C5)*C6*
6512 t4) -
6513 3.*C1*(C4*C4)*t3*
6514 (12.*t2 +
6515 (C5*C5)*C6*
6516 (2.*apkk1 +
6519 amLam*
6521 )*t4))) -
6522 18.*akpk*C1*C2*C4*C7*Ekaon*t3*t5 +
6524 18.*apkk1*C1*C2*C4*C7*Ekaon*t3*t5 +
6525 18.*akpk*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 +
6526 18.*apkk1*C1*C2*(C4*C4)*C7*Ekaon*t3*t5 -
6531 72.*akpk*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
6538 18.*akpk*C2*C5*C6*C7*Ekaon*t4*t5 -
6540 18.*apkk1*C2*C5*C6*C7*Ekaon*t4*t5 +
6541 6.*akpk*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
6542 6.*apkk1*C2*(C5*C5)*C6*C7*Ekaon*t4*t5 +
6547 24.*akpk*apkk1*C2*C5*C6*C7*Fm2*t4*t5 -
6554 18.*akpk*C1*C4*C8*C9*Ekaon*t3*t6 -
6556 18.*apkk1*C1*C4*C8*C9*Ekaon*t3*t6 -
6557 18.*akpk*C1*(C4*C4)*C8*C9*Ekaon*t3*t6 -
6558 18.*apkk1*C1*(C4*C4)*C8*C9*Ekaon*t3*t6 +
6563 72.*akpk*apkk1*C1*C4*C8*C9*Fm1*t3*t6 -
6570 18.*akpk*C5*C6*C8*C9*Ekaon*t4*t6 +
6572 18.*apkk1*C5*C6*C8*C9*Ekaon*t4*t6 -
6573 6.*akpk*(C5*C5)*C6*C8*C9*Ekaon*t4*t6 -
6574 6.*apkk1*(C5*C5)*C6*C8*C9*Ekaon*t4*t6 -
6579 24.*akpk*apkk1*C5*C6*C8*C9*Fm2*t4*t6 +
6587 (t3*t3) +
6590 (t3*t3) +
6594 (t3*t3) -
6596 (t3*t3) +
6598 (t3*t3) -
6600 (t3*t3) -
6602 (t3*t3) +
6604 (t3*t3) +
6606 (t3*t3) -
6608 (t3*t3) -
6610 (t3*t3) +
6622 t4 -
6625 t3*t4 +
6627 t4 -
6629 t3*t4 -
6631 t3*t4 +
6633 t4 +
6635 Fm2*t3*t4 -
6637 t3*t4 -
6639 t3*t4 -
6641 t3*t4 +
6643 t4 +
6645 t4 +
6647 t3*t4 -
6654 (t4*t4) -
6657 (t4*t4) -
6659 (t4*t4) -
6661 (t4*t4) +
6663 (t4*t4) +
6665 (t4*t4) -
6667 (t4*t4) -
6669 (t4*t4) -
6670 2.*C3*
6672 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
6673 2.*akpk*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
6674 (-3.*t2 +
6676 ) +
6677 Ekaon*
6679 (t3*t3) +
6680 (C5*C5)*C6*t4*
6683 6.*C1*(C4*C4)*t3*
6687 (9.*(C1*C1)*(C4*C4*C4*C4)*
6691 (C5*C5)*C6*t4*
6692 (-6.*t2 +
6693 (C5*C5)*C6*
6697 3.*C1*(C4*C4)*t3*
6698 (6.*t2 +
6699 (C5*C5)*C6*
6703 *t4))) +
6704 (C3*C3)*
6706 36.*(akpk*akpk)*C1*(C4*C4)*t2*t3 +
6708 72.*akpk*apkk1*C1*(C4*C4)*t2*t3 -
6710 36.*(apkk1*apkk1)*C1*(C4*C4)*t2*t3 -
6716 (t3*t3) -
6719 (t3*t3) -
6721 (t3*t3) +
6724 (t3*t3) + 12.*(akpk*akpk)*(C5*C5)*C6*t2*t4 +
6726 24.*akpk*apkk1*(C5*C5)*C6*t2*t4 -
6728 12.*(apkk1*apkk1)*(C5*C5)*C6*t2*t4 -
6733 t4 +
6735 t3*t4 -
6737 t3*t4 +
6739 Enu*t3*t4 +
6741 Enu*t3*t4 -
6743 C6*Enu*t3*t4 -
6745 C6*Enu*t3*t4 +
6747 t4 -
6752 (t4*t4) -
6755 (t4*t4) -
6757 (t4*t4) +
6760 (t4*t4) +
6762 (36.*(t2*t2) -
6764 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6767 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6768 - 1.*
6769 ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))) -
6770 2.*Elep*
6771 (2.*(akpk - 1.*apkk1)*
6772 (3.*amSig*C1*(C4*C4)*t3 +
6773 amLam*(C5*C5)*C6*t4)*
6774 (-3.*t2 +
6776 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
6777 Enu*
6778 (36.*(t2*t2) -
6780 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6783 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6784 - 1.*
6787 2.*Ekaon*
6788 (6.*(akpk - 1.*apkk1)*t2*
6789 (3.*amSig*C1*(C4*C4)*t3 +
6790 amLam*(C5*C5)*C6*t4) +
6791 Enu*
6792 (-36.*(t2*t2) +
6794 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) -
6795 1.*(akpk - 1.*apkk1)*
6797 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6798 + ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))) +
6799 Elep*
6800 (36.*(t2*t2) -
6802 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6803 (akpk - 1.*apkk1)*
6805 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))
6806 + ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))))) -
6807 9.*C1*C2*C4*C7*Ekaon*t3*t5 +
6808 9.*akpk*C1*C2*C4*C7*Fm1*t3*t5 +
6810 9.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 -
6814 9.*C2*C5*C6*C7*Ekaon*t4*t5 -
6815 3.*akpk*C2*C5*C6*C7*Fm2*t4*t5 -
6817 3.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
6821 9.*C1*C4*C8*C9*Ekaon*t3*t6 -
6822 9.*akpk*C1*C4*C8*C9*Fm1*t3*t6 -
6824 9.*apkk1*C1*C4*C8*C9*Fm1*t3*t6 +
6828 9.*C5*C6*C8*C9*Ekaon*t4*t6 +
6829 3.*akpk*C5*C6*C8*C9*Fm2*t4*t6 +
6831 3.*apkk1*C5*C6*C8*C9*Fm2*t4*t6 -
6838 36.*akpk*amSig*C1*(C4*C4)*t2*t3 -
6840 36.*am*apkk1*C1*(C4*C4)*t2*t3 +
6841 36.*amSig*apkk1*C1*(C4*C4)*t2*t3 -
6856 12.*akpk*am*(C5*C5)*C6*t2*t4 -
6858 12.*akpk*amLam*(C5*C5)*C6*t2*t4 -
6860 12.*am*apkk1*(C5*C5)*C6*t2*t4 +
6861 12.*amLam*apkk1*(C5*C5)*C6*t2*t4 -
6868 C6*t3*t4 +
6872 C6*t3*t4 +
6875 t3*t4 -
6879 t3*t4 -
6884 Enu*t3*t4 -
6886 t3*t4 +
6901 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
6903 (9.*(C1*C1)*(C4*C4*C4*C4)*
6906 6.*C1*(C4*C4)*t3*
6907 (-12.*t2 +
6908 (C5*C5)*C6*
6911 )*t4) +
6912 (C5*C5)*C6*t4*
6913 (-24.*t2 +
6914 (C5*C5)*C6*
6917 + Ekaon*
6918 (108.*(t2*t2) -
6930 t4 -
6932 t4 +
6934 t4 +
6937 t4 +
6939 t4 -
6950 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) -
6951 2.*akpk*
6953 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
6954 1.*
6955 ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)))) +
6956 Elep*(-36.*(t2*t2) +
6958 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
6961 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
6963 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
6964 2.*am*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
6965 (3.*amSig*C1*(C4*C4)*t3 +
6966 amLam*(C5*C5)*C6*t4) +
6967 ((3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4)*(3.*amSig*C1*(C4*C4)*t3 + amLam*(C5*C5)*C6*t4))))) +
6969 (3.*C1*C4*
6971 + C5*C6*
6973 t4) +
6979 amSig*
6980 (-4.*Ekaon +
6981 2.*(akpk - 1.*apkk1)*Fm1 +
6983 + am*
6985 2.*(akpk - 1.*apkk1)*Fm1 +
6986 4.*Ekaon*
6988 2.*amSig*
6991 6.*C1*C4*t3*
7007 (-1.*Fm2 +
7008 Fm1*
7016 am*
7023 amSig*C5*C6*
7024 (-1.*Fm2 +
7025 Fm1*
7028 amLam*C5*C6*
7029 (-1.*Fm2 +
7030 Fm1*
7034 + C5*C6*Fm2*t4*
7047 am*
7049 2.*amLam*C5*C6*
7052 2.*
7053 (C5*C6*
7054 ((akpk - 1.*apkk1)*Fm2 +
7055 2.*Ekaon*
7057 3.*C2*C7*t5 + 3.*C8*C9*t6)))) -
7058 2.*Ekaon*
7061 (-12.*f*t1 +
7062 (akpk - 1.*apkk1)*
7066 (akpk*
7068 apkk1*
7070 2.*Ekaon*
7074 t4)) -
7075 6.*am*
7077 C5*C6*Fm2*t4*
7078 (-2.*f*t1 +
7079 Ekaon*
7080 (2.*amLam*C5*C6*t4 + C2*C7*t5 -
7081 1.*C8*C9*t6)) +
7082 C1*C4*t3*
7084 1.*Ekaon*
7087 3.*Fm1*(C2*C7*t5 - 1.*C8*C9*t6)))))
7088 )) + 2.*(Enu*
7136 (t3*t3) -
7140 (t3*t3) -
7144 (t3*t3) -
7151 (t3*t3) -
7153 (t3*t3) -
7157 (t3*t3) -
7159 (t3*t3) -
7161 (t3*t3) -
7163 (t3*t3) +
7202 t3*t4 -
7206 C6*t3*t4 +
7208 t3*t4 +
7213 t3*t4 -
7216 C6*t3*t4 -
7219 t3*t4 +
7221 t3*t4 +
7223 t3*t4 +
7229 Ekaon*t3*t4 -
7232 Ekaon*t3*t4 -
7234 t3*t4 -
7238 Ekaon*t3*t4 +
7241 Ekaon*t3*t4 -
7244 Ekaon*t3*t4 +
7246 Ekaon*t3*t4 +
7250 t3*t4 -
7256 t3*t4 +
7259 t3*t4 +
7262 t3*t4 -
7264 t3*t4 -
7266 Elep*t3*t4 +
7273 Fm1*t3*t4 +
7276 Fm1*t3*t4 +
7278 t3*t4 +
7280 t3*t4 -
7282 t3*t4 +
7284 Fm1*t3*t4 -
7286 Fm1*t3*t4 +
7288 Fm1*t3*t4 -
7290 Fm1*t3*t4 +
7292 Fm1*t3*t4 -
7294 Fm1*t3*t4 -
7296 Fm1*t3*t4 +
7298 Fm1*t3*t4 -
7304 t3*t4 -
7306 t3*t4 -
7308 Fm1*t3*t4 +
7310 Fm1*t3*t4 -
7312 Fm1*t3*t4 +
7314 Fm1*t3*t4 +
7321 Fm2*t3*t4 -
7323 t3*t4 +
7325 t3*t4 -
7327 t3*t4 +
7329 t3*t4 +
7331 Fm2*t3*t4 -
7333 Fm2*t3*t4 -
7335 t3*t4 +
7337 t3*t4 +
7339 Fm2*t3*t4 -
7341 Fm2*t3*t4 +
7343 Fm2*t3*t4 -
7345 Fm2*t3*t4 -
7351 t3*t4 +
7353 t3*t4 -
7355 Fm2*t3*t4 +
7357 Fm2*t3*t4 +
7359 Fm2*t3*t4 -
7361 Fm2*t3*t4 +
7380 (t4*t4) +
7384 (t4*t4) -
7388 (t4*t4) -
7395 (t4*t4) -
7397 (t4*t4) -
7401 (t4*t4) -
7403 (t4*t4) -
7405 (t4*t4) -
7407 (t4*t4) +
7412 18.*am*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 -
7413 18.*amSig*(apkk1*apkk1)*C1*C2*C4*C7*t3*t5 -
7414 18.*am*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 -
7415 18.*amSig*(apkk1*apkk1)*C1*C2*(C4*C4)*C7*t3*t5 -
7420 18.*am*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 +
7421 18.*amLam*(apkk1*apkk1)*C2*C5*C6*C7*t4*t5 -
7422 6.*am*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 -
7423 6.*amLam*(apkk1*apkk1)*C2*(C5*C5)*C6*C7*t4*t5 +
7428 18.*am*(apkk1*apkk1)*C1*C4*C8*C9*t3*t6 +
7429 18.*amSig*(apkk1*apkk1)*C1*C4*C8*C9*t3*t6 +
7430 18.*am*(apkk1*apkk1)*C1*(C4*C4)*C8*C9*t3*t6 +
7431 18.*amSig*(apkk1*apkk1)*C1*(C4*C4)*C8*C9*t3*t6 +
7436 18.*am*(apkk1*apkk1)*C5*C6*C8*C9*t4*t6 -
7437 18.*amLam*(apkk1*apkk1)*C5*C6*C8*C9*t4*t6 +
7438 6.*am*(apkk1*apkk1)*(C5*C5)*C6*C8*C9*t4*t6 +
7439 6.*amLam*(apkk1*apkk1)*(C5*C5)*C6*C8*C9*t4*t6 -
7444 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
7446 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
7447 (3.*t2 +
7449 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
7450 12.*t1*
7451 (6.*t2 -
7453 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
7455 (9.*(C1*C1)*(C4*C4*C4*C4)*
7458 (C5*C5)*C6*t4*
7459 (6.*t2 +
7460 (C5*C5)*C6*
7463 6.*C1*(C4*C4)*t3*
7464 (3.*t2 +
7465 (C5*C5)*C6*
7468 ) -
7471 (C2*C7*t5 - 1.*C8*C9*t6)) +
7477 (t3*t3) -
7479 (t3*t3) -
7484 (t3*t3) -
7486 (t3*t3) +
7488 (t3*t3) +
7491 (t3*t3) -
7494 (t3*t3) -
7496 (t3*t3) +
7501 t3*t4 -
7503 C6*Ekaon*t3*t4 +
7508 t4 +
7512 Elep*t3*t4 +
7514 Enu*t3*t4 +
7516 t3*t4 -
7518 t3*t4 +
7520 Fm1*t3*t4 -
7522 Fm1*t3*t4 -
7525 t3*t4 -
7527 t4 -
7529 Fm1*t3*t4 -
7531 Fm1*t3*t4 -
7533 t4 +
7535 t4 +
7537 Fm2*t3*t4 -
7539 Fm2*t3*t4 -
7542 t3*t4 +
7544 t4 +
7546 Fm2*t3*t4 +
7548 Fm2*t3*t4 +
7550 t3*t4 +
7552 (t4*t4) -
7554 (t4*t4) -
7559 (t4*t4) -
7561 (t4*t4) +
7563 (t4*t4) +
7566 (t4*t4) +
7569 (t4*t4) -
7571 (t4*t4) +
7573 (9.*(C1*C1)*(C4*C4)*
7574 (-1. + 2.*C4 +
7575 (C4*C4)*
7576 (-1. + apkk1*C3 +
7579 3.*C1*C4*t3*
7580 (C5*C6*
7581 (6. + 2.*C5 - 2.*C4*(3. + C5) -
7586 2.*C3*C4*
7587 (3.*t2 +
7588 (C5*C5)*C6*
7590 )) +
7591 (C5*C5)*C6*t4*
7592 (-1.*C6*
7593 (9. + 6.*C5 + (C5*C5) +
7595 C3*
7596 (6.*t2 +
7597 (C5*C5)*C6*
7599 ))) -
7600 12.*t1*
7601 (apkk1*(2. + f)*
7605 (apkk1*C3*
7606 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
7607 Ekaon*
7608 (-3.*C1*C4*
7611 C5*C6*
7615 (-3.*C1*C4*
7616 (-2. + 2.*C4*(1. + f) +
7618 - 1.*C5*C6*
7619 (6. + 2.*C5*(1. + f) -
7621 *t4 +
7622 C3*
7623 (6.*t2 +
7625 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))) -
7626 1.*am*
7627 (-3.*C1*C4*
7629 2.*Ekaon*
7632 1.*C5*C6*
7634 2.*Ekaon*
7637 2.*C3*
7639 3.*amSig*apkk1*C1*(C4*C4)*t3 +
7640 amLam*apkk1*(C5*C5)*C6*t4 +
7641 Ekaon*
7642 (6.*t2 + 3.*apkk1*C1*(C4*C4)*t3 +
7643 apkk1*(C5*C5)*C6*t4)))) -
7644 9.*am*apkk1*C1*C2*C4*C7*t3*t5 -
7645 9.*amSig*apkk1*C1*C2*C4*C7*t3*t5 +
7647 9.*am*apkk1*C2*C5*C6*C7*t4*t5 +
7648 9.*amLam*apkk1*C2*C5*C6*C7*t4*t5 -
7650 9.*am*apkk1*C1*C4*C8*C9*t3*t6 +
7651 9.*amSig*apkk1*C1*C4*C8*C9*t3*t6 -
7653 9.*am*apkk1*C5*C6*C8*C9*t4*t6 -
7654 9.*amLam*apkk1*C5*C6*C8*C9*t4*t6 +
7658 (-1. + 2.*C4 + (-1. + apkk1*C3)*(C4*C4))*
7659 (t3*t3) -
7662 t4 +
7669 (t4*t4) +
7671 (9.*(C1*C1)*(C4*C4)*
7674 6.*C1*C4*C5*C6*
7677 t4 +
7678 (C5*C5)*(C6*C6)*
7682 (-9.*(C1*C1)*(C4*C4)*
7683 (-1. +
7684 (C4*C4)*
7685 (-1. + 3.*apkk1*C3 -
7690 C4*
7693 (t3*t3) -
7694 6.*C1*C4*C5*C6*
7701 C4*
7702 (-3. +
7703 C5*
7704 (-1. + 3.*apkk1*C3 -
7712 (C5*C5)*(C6*C6)*
7713 (9. +
7714 (C5*C5)*
7715 (1. - 3.*apkk1*C3 +
7720 C5*
7723 (t4*t4)) -
7724 27.*apkk1*C1*C2*C4*C7*Fm1*t3*t5 +
7725 9.*amLam*C2*C5*C6*C7*t4*t5 +
7726 6.*amLam*C2*(C5*C5)*C6*C7*t4*t5 +
7727 9.*apkk1*C2*C5*C6*C7*Fm2*t4*t5 +
7728 27.*apkk1*C1*C4*C8*C9*Fm1*t3*t6 -
7729 9.*amLam*C5*C6*C8*C9*t4*t6 -
7730 6.*amLam*(C5*C5)*C6*C8*C9*t4*t6 -
7731 9.*apkk1*C5*C6*C8*C9*Fm2*t4*t6 -
7732 3.*amSig*C1*C4*t3*
7733 (2.*amLam*C5*
7734 (3. + C5 +
7735 C4*(-3. - 1.*C5 + apkk1*C3*C5))*C6*
7736 t4 +
7737 apkk1*C5*C6*
7739 2.*C4*Fm2)*t4 -
7740 3.*(-1. + 2.*C4)*
7741 (C2*C7*t5 - 1.*C8*C9*t6)) +
7742 am*
7743 (18.*(C1*C1)*(C4*C4)*
7744 (2.*
7745 (1. - 2.*C4 +
7746 (1. - 1.*apkk1*C3)*(C4*C4))*Ekaon +
7747 amSig*
7748 (1. + (1. - 2.*apkk1*C3)*(C4*C4) +
7750 ) +
7754 3.*C1*C4*t3*
7755 (6.*amSig*C5*C6*t4 -
7756 6.*amSig*C4*C5*C6*t4 +
7757 2.*amSig*(C5*C5)*C6*t4 -
7758 2.*amSig*C4*(C5*C5)*C6*t4 +
7759 24.*C5*C6*Ekaon*t4 -
7760 24.*C4*C5*C6*Ekaon*t4 +
7761 8.*(C5*C5)*C6*Ekaon*t4 -
7762 8.*C4*(C5*C5)*C6*Ekaon*t4 -
7767 2.*amLam*C5*C6*
7769 C5*
7772 C4*
7775 4.*C3*C4*
7777 3.*Enu*t2 +
7778 amLam*apkk1*(C5*C5)*C6*t4 +
7779 amSig*apkk1*(C5*C5)*C6*t4 +
7780 2.*apkk1*(C5*C5)*C6*Ekaon*t4 -
7783 3.*C2*C7*t5 - 6.*C2*C4*C7*t5 -
7784 3.*C8*C9*t6 + 6.*C4*C8*C9*t6) +
7785 C5*C6*t4*
7786 (36.*C5*C6*Ekaon*t4 +
7787 24.*(C5*C5)*C6*Ekaon*t4 +
7788 4.*(C5*C5*C5)*C6*Ekaon*t4 +
7789 2.*amLam*C5*C6*
7790 (9. + 6.*C5 + (C5*C5) +
7794 2.*C3*C5*
7796 2.*amLam*apkk1*(C5*C5)*C6*t4 -
7799 2.*Ekaon*
7800 (6.*t2 + apkk1*(C5*C5)*C6*t4)) +
7801 9.*C2*C7*t5 + 6.*C2*C5*C7*t5 -
7802 9.*C8*C9*t6 - 6.*C5*C8*C9*t6)))))
7803 + 2.*(akpk*akpk)*
7806 (t3*t3) +
7808 t4 +
7810 + 3.*Elep*
7812 (C2*C7*t5 - 1.*C8*C9*t6) +
7813 amSig*C1*C4*t3*
7815 Fm1*
7816 (12.*(-1. + f)*t1 +
7818 3.*(-1. + C4)*(C2*C7*t5 - 1.*C8*C9*t6))
7819 - 1.*amLam*C5*C6*t4*
7821 Fm2*
7822 (4.*(-1. + f)*t1 +
7824 (3. + C5)*(C2*C7*t5 - 1.*C8*C9*t6))) -
7828 3.*amSig*C1*C4*t3*
7829 (2.*C3*C4*
7830 (-6.*t1 - 3.*t2 +
7834 Fm1*
7838 + 3.*C2*C7*t5 - 3.*C8*C9*t6)) +
7839 amLam*C5*C6*t4*
7840 (2.*C3*C5*
7841 (-6.*t1 - 3.*t2 +
7845 Fm2*
7849 - 3.*C2*C7*t5 + 3.*C8*C9*t6))) +
7853 C5*C6*t4*
7854 (C3*
7855 (6.*C5*(2.*t1 + t2) -
7857 Fm2*
7859 3.*C2*C7*t5 - 3.*C8*C9*t6)) +
7860 3.*C1*C4*t3*
7861 (2.*C3*C4*
7862 (6.*t1 + 3.*t2 -
7864 Fm1*
7866 3.*C2*C7*t5 + 3.*C8*C9*t6))) +
7867 3.*Elep*
7868 (C5*C6*t4*
7869 (-1.*(3. + C5)*
7870 (C2*C7*t5 - 1.*C8*C9*t6) +
7871 Fm2*
7872 (4.*t1 - 4.*f*t1 +
7873 2.*C2*C7*Ekaon*t5 -
7874 2.*C8*C9*Ekaon*t6)) +
7875 C1*C4*t3*
7879 t4 + 3.*C2*C7*t5 - 3.*C2*C4*C7*t5 -
7880 3.*C8*C9*t6 + 3.*C4*C8*C9*t6 +
7881 2.*Fm1*
7882 (6.*(-1. + f)*t1 +
7883 Ekaon*
7884 (amLam*(-3. + C5)*C5*C6*t4 -
7885 1.*amSig*(-3. + C5)*C5*C6*t4 -
7886 3.*C2*C7*t5 + 3.*C8*C9*t6)))))) -
7888 (amLam*C5*C6*
7890 1.*Fm2*
7891 (4.*(1. + f)*t1 -
7893 amSig*C1*C4*t3*
7895 Fm1*
7896 (12.*(1. + f)*t1 -
7898 2.*am*
7900 C1*C4*t3*
7903 1.*Fm1*
7904 (6.*(1. + f)*t1 -
7906 C6*Ekaon*t4)))) +
7908 (t3*t3) +
7910 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
7911 3.*amSig*C1*C4*t3*
7912 (2.*C3*C4*
7913 (-6.*t1 - 3.*t2 +
7918 1.*amLam*C5*C6*t4*
7919 (2.*C3*C5*
7920 (-6.*t1 - 3.*t2 +
7924 3.*Fm2*(C2*C7*t5 - 1.*C8*C9*t6)) +
7926 C5*C6*t4*
7927 (C3*
7928 (6.*C5*(2.*t1 + t2) -
7931 3.*C1*C4*t3*
7932 (2.*C3*C4*
7933 (6.*t1 + 3.*t2 -
7935 3.*Fm1*(C2*C7*t5 - 1.*C8*C9*t6)))) +
7937 72.*amSig*C1*(C4*C4)*t1*t3 +
7954 (t3*t3) +
7956 (t3*t3) +
7958 (t3*t3) -
7961 (t3*t3) - 72.*amLam*C5*C6*t1*t4 -
7962 24.*amLam*(C5*C5)*C6*t1*t4 -
7979 t4 -
7981 t4 +
7987 t3*t4 -
7989 t4 +
7991 t3*t4 +
7997 t3*t4 -
7999 t4 -
8001 t3*t4 -
8003 t3*t4 +
8006 Fm2*t3*t4 +
8018 (t4*t4) -
8020 (t4*t4) -
8023 (t4*t4) -
8029 (9.*(C1*C1)*(C4*C4)*
8030 (2.*((-1.+C4)*(-1.+C4))*Ekaon +
8033 *(t3*t3) +
8034 C5*C6*t4*
8035 (2.*C5*((3. + C5)*(3. + C5))*C6*Ekaon*t4 +
8038 3.*Fm2*
8039 (8.*(1. + f)*t1 -
8041 + 3.*C1*C4*t3*
8042 (C5*C6*
8044 2.*Ekaon*
8045 (-6. - 2.*C5 + 2.*C4*(3. + C5) +
8047 Fm1*
8048 (24.*(1. + f)*t1 -
8049 1.*C5*C6*
8052 2.*Ekaon*
8055 2.*C3*
8056 ((3.*am*C1*(C4*C4)*t3 +
8057 3.*amSig*C1*(C4*C4)*t3 +
8060 2.*apkk1*
8061 (6.*t1 + 3.*t2 -
8068 Enu*
8070 ((3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*(3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) -
8072 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
8073 (-6.*t2 +
8075 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)) +
8078 (C4*C4*C4*C4)*(t3*t3) +
8079 (C5*C5)*C6*t4*
8080 (-6.*t2 +
8082 - 6.*C1*(C4*C4)*t3*
8083 (3.*t2 +
8085 C6*t4)) -
8086 2.*am*
8088 (t3*t3) +
8090 (t3*t3) +
8091 6.*C1*(C4*C4)*t3*
8092 (-3.*Ekaon*t2 +
8095 3.*amSig*C1*(C4*C4)*t3*
8099 (C5*C5)*C6*t4*
8104 12.*t1*
8105 (6.*t2 +
8107 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) -
8109 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4) +
8110 am*
8111 (3.*amSig*C1*(C4*C4)*t3 +
8112 6.*C1*(C4*C4)*Ekaon*t3 +
8114 - 108.*C2*C7*f*t1*t5 -
8125 18.*(C2*C2)*(C7*C7)*Ekaon*(t5*t5) -
8126 18.*(C2*C2)*(C7*C7)*Enu*(t5*t5) +
8127 108.*C8*C9*f*t1*t6 +
8138 36.*C2*C7*C8*C9*Ekaon*t5*t6 +
8139 36.*C2*C7*C8*C9*Enu*t5*t6 +
8140 18.*(C8*C8)*(C9*C9)*Ekaon*(t6*t6) -
8141 18.*(C8*C8)*(C9*C9)*Enu*(t6*t6) +
8142 am*(18.*(C1*C1)*(C4*C4)*
8148 C4*
8151 (t3*t3) +
8152 2.*C5*C6*t4*
8153 (-12.*
8162 t4 +
8165 C5*
8168 - 9.*C2*C7*Ekaon*t5 +
8169 3.*C2*C5*C7*Ekaon*t5 -
8170 3.*apkk1*C2*C7*Fm2*t5 +
8172 9.*C8*C9*Ekaon*t6 -
8173 3.*C5*C8*C9*Ekaon*t6 +
8174 3.*apkk1*C8*C9*Fm2*t6 -
8176 3.*C1*C4*t3*
8177 (24.*
8184 4.*C5*
8191 4.*C4*
8194 2.*
8198 (2.*Ekaon*
8200 amSig*
8203 + 3.*
8206 (C2*C7*t5 - 1.*C8*C9*t6))))) +
8217 (t3*t3) +
8227 t4 -
8229 t4 +
8231 (C5*C5)*C6*t3*t4 -
8233 t4 +
8235 t4 -
8237 t4 +
8239 t4 -
8241 t3*t4 +
8243 t3*t4 -
8245 t3*t4 +
8247 t3*t4 +
8249 t4 -
8251 t4 +
8253 t4 -
8255 t4 -
8257 t3*t4 +
8259 t3*t4 -
8261 t3*t4 +
8271 (t4*t4) -
8273 (9.*(C1*C1)*(C4*C4)*
8274 (-1. + 2.*C4 +
8276 6.*C1*C4*C5*
8277 (3. + C5 +
8279 t3*t4 +
8280 (C5*C5)*
8281 (-9. - 6.*C5 +
8285 (3.*(C1*C1)*(C4*C4)*
8286 (3. - 2.*C4 + 3.*(C4*C4))*(t3*t3) +
8287 2.*C1*C4*C5*
8288 (-9. - 1.*C5 + 3.*C4*(1. + C5))*C6*
8289 t3*t4 +
8290 (C5*C5)*(9. + 2.*C5 + (C5*C5))*(C6*C6)*
8291 (t4*t4)) -
8292 4.*Ekaon*
8293 (9.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
8294 (t3*t3) +
8295 3.*C1*C4*C5*C6*
8296 (2.*(1. + C4)*(-3. + C5)*Ekaon +
8297 amSig*
8301 amLam*
8305 ((-3. + C5)*(-3. + C5))*(C5*C5)*(C6*C6)*Ekaon*(t4*t4)
8306 ) -
8308 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4)*
8309 (-12.*t2 +
8311 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))) -
8312 12.*t1*
8313 (-2.*
8315 (3.*C1*C4*(1. + C4)*t3 +
8316 (-3. + C5)*C5*C6*t4) -
8317 1.*apkk1*(1. + f)*
8321 (3.*C1*C4*
8324 C5*C6*
8327 am*
8328 (-3.*C1*C4*
8331 3.*amSig*C1*C4*
8334 C5*C6*
8337 amLam*
8341 ((2. + f)*
8346 C3*
8347 (12.*t2 -
8349 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))))
8352 18.*amSig*apkk1*C1*C2*C4*C7*t3*t5 -
8353 18.*amSig*apkk1*C1*C2*(C4*C4)*C7*t3*t5 -
8358 18.*amLam*apkk1*C2*C5*C6*C7*t4*t5 -
8359 6.*amLam*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
8363 36.*apkk1*(C2*C2)*(C7*C7)*(t5*t5) -
8366 18.*amSig*apkk1*C1*C4*C8*C9*t3*t6 +
8367 18.*amSig*apkk1*C1*(C4*C4)*C8*C9*t3*t6 +
8372 18.*amLam*apkk1*C5*C6*C8*C9*t4*t6 +
8373 6.*amLam*apkk1*(C5*C5)*C6*C8*C9*t4*t6 -
8377 72.*apkk1*C2*C7*C8*C9*t5*t6 +
8379 36.*apkk1*(C8*C8)*(C9*C9)*(t6*t6) -
8380 1.*am*
8383 (t3*t3) +
8385 (t3*t3) +
8386 6.*amSig*C1*C4*t3*
8387 (C5*C6*Ekaon*
8389 2.*C4*Fm2)*t4 +
8390 2.*C3*C4*
8391 (-3.*t2 +
8394 (3.*C1*(C4*C4)*t3 + (C5*C5)*C6*t4))) +
8395 3.*C1*C4*t3*
8398 2.*C4*Fm2)*t4 +
8399 4.*C3*C4*
8400 (-6.*Ekaon*t2 +
8404 (C2*C7*t5 - 1.*C8*C9*t6)) +
8405 C5*C6*t4*
8407 4.*C3*C5*Ekaon*
8409 4.*amLam*C3*C5*
8412 (C2*C7*t5 - 1.*C8*C9*t6))) +
8413 2.*
8416 (t3*t3) +
8418 (t3*t3) -
8420 t4 +
8422 t3*t4 -
8424 t3*t4 +
8426 Ekaon*t3*t4 -
8428 Fm1*t3*t4 +
8430 Fm1*t3*t4 -
8432 Fm1*t3*t4 +
8434 Fm1*t3*t4 +
8436 t3*t4 -
8438 t3*t4 -
8440 Fm2*t3*t4 +
8442 Fm2*t3*t4 +
8444 (t4*t4) +
8446 (t4*t4) +
8448 9.*apkk1*C1*C2*C4*C7*t3*t5 +
8449 9.*apkk1*C1*C2*(C4*C4)*C7*t3*t5 -
8451 t5 - 9.*apkk1*C2*C5*C6*C7*t4*t5 +
8452 3.*apkk1*C2*(C5*C5)*C6*C7*t4*t5 +
8454 t5 - 9.*apkk1*C1*C4*C8*C9*t3*t6 -
8455 9.*apkk1*C1*(C4*C4)*C8*C9*t3*t6 +
8457 t6 + 9.*apkk1*C5*C6*C8*C9*t4*t6 -
8458 3.*apkk1*(C5*C5)*C6*C8*C9*t4*t6 -
8460 t6 +
8462 (18.*(C1*C1)*(C4*C4)*((1.+C4)*(1.+C4))*Ekaon*
8463 (t3*t3) +
8464 3.*amSig*C1*C4*t3*
8465 (6.*C1*C4*(1. + (C4*C4))*t3 +
8466 C5*C6*
8470 + C5*C6*t4*
8471 (2.*amLam*C5*(9. + (C5*C5))*C6*t4 +
8472 18.*C5*C6*Ekaon*t4 -
8473 12.*(C5*C5)*C6*Ekaon*t4 +
8474 2.*(C5*C5*C5)*C6*Ekaon*t4 +
8475 9.*C2*C7*t5 - 9.*C8*C9*t6) +
8476 3.*C1*C4*t3*
8477 (-12.*(1. + C4)*C5*C6*Ekaon*t4 +
8478 4.*(1. + C4)*(C5*C5)*C6*Ekaon*t4 +
8479 amLam*C5*C6*
8483 3.*C2*C7*t5 + 3.*C8*C9*t6))))) -
8484 2.*apkk1*
8493 t4 +
8495 t4 +
8497 t4 +
8504 (9.*(C1*C1)*(C4*C4)*(1. + (C4*C4))*(t3*t3) +
8505 6.*C1*C4*C5*(-3. + C4*C5)*C6*t3*t4 +
8506 (C5*C5)*(9. + (C5*C5))*(C6*C6)*(t4*t4)) +
8507 72.*C2*C7*f*t1*t5 +
8515 18.*(C2*C2)*(C7*C7)*Ekaon*(t5*t5) +
8516 18.*(C2*C2)*(C7*C7)*Enu*(t5*t5) -
8517 72.*C8*C9*f*t1*t6 -
8525 36.*C2*C7*C8*C9*Ekaon*t5*t6 -
8526 36.*C2*C7*C8*C9*Enu*t5*t6 -
8527 18.*(C8*C8)*(C9*C9)*Ekaon*(t6*t6) +
8528 18.*(C8*C8)*(C9*C9)*Enu*(t6*t6) -
8529 3.*amSig*C1*C4*t3*
8540 C4*
8541 (-24.*t1 +
8544 (6.*C1*t3 +
8547 - 1.*am*
8549 (t3*t3) -
8550 3.*C1*C4*t3*
8560 3.*C2*C7*Enu*t5 +
8562 3.*C8*C9*Enu*t6 -
8564 C4*
8566 Enu*
8569 3.*C2*C7*t5 - 3.*C8*C9*t6))) +
8570 C5*C6*t4*
8572 3.*C5*
8575 3.*
8577 Enu*
8578 (-3.*C2*C7*t5 + 3.*C8*C9*t6 +
8579 2.*Fm2*
8580 (2.*t1 + C2*C7*Ekaon*t5 -
8581 1.*C8*C9*Ekaon*t6))))))))));
8582
8583 if (sol <= 0.0) {
8585 << "Matrix element is 0 or negative! This is bad and really should throw an exception!";
8586 }
8587
8588 return sol;
8589}
References alepvec, am, amEta, amk, aml, amLam, ampi, amSig, angkq, aq0, aqvec, d, Ekaon, Elep, Enu, f, Fm1, Fm2, fpi, g, LOG, pERROR, pkvec, and Vus.
Referenced by XSec().
◆ Configure() [1/2]
|
virtual |
Configure the algorithm with an external registry The registry is merged with the top level registry if it is owned, Otherwise a copy of it is added with the highest priority
Reimplemented from genie::Algorithm.
Definition at line 182 of file AlamSimoAtharVacasSKPXSec2014.cxx.
183{
184 Algorithm::Configure(config);
185 this->LoadConfig();
186}
void LoadConfig(void)
Definition AlamSimoAtharVacasSKPXSec2014.cxx:194
virtual void Configure(const Registry &config)
Definition Algorithm.cxx:62
References genie::Algorithm::Configure(), and LoadConfig().
◆ Configure() [2/2]
|
virtual |
Configure the algorithm from the AlgoConfigPool based on param_set string given in input An algorithm contains a vector of registries coming from different xml configuration files, which are loaded according a very precise prioriy This methods will load a number registries in order of priority: 1) "Tunable" parameter set from CommonParametes. This is loaded with the highest prioriry and it is designed to be used for tuning procedure Usage not expected from the user. 2) For every string defined in "CommonParame" the corresponding parameter set will be loaded from CommonParameter.xml 3) parameter set specified by the config string and defined in the xml file of the algorithm 4) if config is not "Default" also the Default parameter set from the same xml file will be loaded Effectively this avoids the repetion of a parameter when it is not changed in the requested configuration
Reimplemented from genie::Algorithm.
Definition at line 188 of file AlamSimoAtharVacasSKPXSec2014.cxx.
References genie::Algorithm::Configure(), and LoadConfig().
◆ Integral()
|
virtual |
Integrate the model over the kinematic phase space available to the input interaction (kinematical cuts can be included)
Implements genie::XSecAlgorithmI.
Definition at line 166 of file AlamSimoAtharVacasSKPXSec2014.cxx.
167{
169 return xsec;
170}
const XSecIntegratorI * fXSecIntegrator
cross section integrator
Definition AlamSimoAtharVacasSKPXSec2014.h:52
References fXSecIntegrator.
◆ LoadConfig()
|
private |
Definition at line 194 of file AlamSimoAtharVacasSKPXSec2014.cxx.
195{
196 fXSecIntegrator =
198 assert(fXSecIntegrator);
199
203
205 // fpi is 0.0924 in Athar's code, use the same one that is already in UserPhysicsOptions
207
211 // we really want the anomolous moment, but the one in UserPhysicsOptions is the full moment, despite the name
213 amup -= 1 ;
215
218
219}
bool GetParam(const RgKey &name, T &p, bool is_top_call=true) const
const Algorithm * SubAlg(const RgKey ®istry_key) const
Definition Algorithm.cxx:345
TParticlePDG * Find(int pdgc, bool must_exist=true)
Definition PDGLibrary.cxx:86
static const double kNeutronMass
Definition Framework/Conventions/Constants.h:76
References am, amEta, amLam, amun, amup, d, f, genie::PDGLibrary::Find(), Fm1, Fm2, fpi, fXSecIntegrator, g, genie::Algorithm::GetParam(), genie::PDGLibrary::Instance(), genie::constants::kGF, genie::constants::kNeutronMass, genie::kPdgEta, genie::kPdgLambda, genie::Algorithm::SubAlg(), and Vus.
Referenced by Configure(), and Configure().
◆ ValidProcess()
|
virtual |
Can this cross section algorithm handle the input process?
Implements genie::XSecAlgorithmI.
Definition at line 172 of file AlamSimoAtharVacasSKPXSec2014.cxx.
173{
175
176 const ProcessInfo & proc_info = interaction->ProcInfo();
178
179 return true;
180}
References genie::ProcessInfo::IsSingleKaon(), genie::ProcessInfo::IsWeakCC(), genie::kISkipProcessChk, and genie::Interaction::ProcInfo().
Referenced by XSec().
◆ XSec()
|
virtual |
Compute the cross section for the input interaction.
Implements genie::XSecAlgorithmI.
Definition at line 49 of file AlamSimoAtharVacasSKPXSec2014.cxx.
51{
52 // Check whether interaction is valid
55
56 const InitialState & init_state = interaction -> InitState();
57 const Kinematics & kinematics = interaction -> Kine();
58 Target * target = interaction->InitStatePtr()->TgtPtr();
61
62 // Initialisation begins here
63 // --------------------------
64 // mutable variables used in the matrix element calculation
66 leptonPDG = interaction->FSPrimLeptonPdg();
67
68 // Fill reaction type (NN=1, NP=2, PP=3)
70 int numProtons = interaction->ExclTagPtr()->NProtons();
71 int numNeutrons = interaction->ExclTagPtr()->NNeutrons();
72 int kaonPDG = interaction->ExclTagPtr()->StrangeHadronPdg();
76 else {
78 << "Reaction not defined! This should NEVER happen!"
79 << "\n - numProtons = " << numProtons
80 << "\n - numNeutrons = " << numNeutrons
81 << "\n - kaonPDG = " << kaonPDG;
82 return 0.;
83 }
84
89
90 // Set lepton mass
92
93 double theta = TMath::ACos(costheta);
94
95 // Set reaction parameters, which are mutables used in the matrix element calculations
101 }
107 }
113 }
114 else {
115 return 0.;
116 }
117
118 // Cross section calculation begins here
119 // -------------------------------------
121 if (Tkaon > tkmax) return 0.;
122
125
126 double xsec = 0.;
127
132
133 // aqvec is the magnitude of the three-momentum transfer to hadron system
135
136 // this check is basically the longitudinal component of the kaon momentum (in the momentum transfer frame) divided by the total
137 // if it is larger than 1, the kinematics are non-physical so we should return zero
139
140 double amat2; // the matrix element
141 if (fabs(check) <= 1.0) { // so it has to be smaller than 1, but it could be negative if the kaon backscatters in com frame
142 angkq = check;
146 else return 0.;
148 }
149 else {
150 xsec = 0.;
151 }
152
153 // the matrix element calculation is for d4sigma/dtk dpl dcosthetal dphi_kq
154 // we have T_l instead of p_l so we multiply by dp/dT = E/p
156
157 // xsec is now the nucleon-level cross section
158 // There are no fancy nuclear effects for this model, so the nucleus cross section is just
159 // the nucleon XS times the number of nucleons of the appropriate isospin for the process selected
161 else xsec *= nTargetProtons; // PP
162
163 return xsec;
164}
double Amatrix_NN(double theta, double phikq) const
Definition AlamSimoAtharVacasSKPXSec2014.cxx:345
bool ValidProcess(const Interaction *i) const
Can this cross section algorithm handle the input process?
Definition AlamSimoAtharVacasSKPXSec2014.cxx:172
int leptonPDG
Definition AlamSimoAtharVacasSKPXSec2014.h:70
int reactionType
Definition AlamSimoAtharVacasSKPXSec2014.h:70
double Amatrix_NP(double theta, double phikq) const
Definition AlamSimoAtharVacasSKPXSec2014.cxx:1194
double Amatrix_PP(double theta, double phikq) const
Definition AlamSimoAtharVacasSKPXSec2014.cxx:4914
virtual bool ValidKinematics(const Interaction *i) const
Is the input kinematical point a physically allowed one?
Definition XSecAlgorithmI.cxx:40
static const double kPionMass
Definition Framework/Conventions/Constants.h:73
static const double kPi0Mass
Definition Framework/Conventions/Constants.h:74
static const double kProtonMass
Definition Framework/Conventions/Constants.h:75
References alepvec, am, Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), amk, aml, ampi, amSig, angkq, aq0, aqvec, Ekaon, Elep, Enu, genie::Interaction::ExclTagPtr(), genie::PDGLibrary::Find(), genie::Interaction::FSPrimLeptonPdg(), genie::Interaction::InitStatePtr(), genie::PDGLibrary::Instance(), genie::kKVctl, genie::kKVphikq, genie::kKVTk, genie::kKVTl, genie::constants::kNeutronMass, genie::kPdgK0, genie::kPdgKP, genie::kPdgSigma0, genie::kPdgSigmaM, genie::constants::kPi, genie::constants::kPi0Mass, genie::constants::kPionMass, genie::constants::kProtonMass, genie::kRfHitNucRest, leptonPDG, LOG, genie::Target::N(), genie::XclsTag::NNeutrons(), genie::XclsTag::NProtons(), pERROR, pkvec, genie::InitialState::ProbeE(), reactionType, genie::XclsTag::StrangeHadronPdg(), genie::InitialState::TgtPtr(), genie::XSecAlgorithmI::ValidKinematics(), ValidProcess(), and genie::Target::Z().
Member Data Documentation
◆ alepvec
|
private |
Definition at line 77 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ am
|
private |
Definition at line 71 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), LoadConfig(), and XSec().
◆ amEta
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_PP(), and LoadConfig().
◆ amk
|
private |
Definition at line 71 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ aml
|
mutableprivate |
Definition at line 71 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ amLam
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ ampi
|
private |
Definition at line 71 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ amSig
|
private |
Definition at line 71 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ amun
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by LoadConfig().
◆ amup
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by LoadConfig().
◆ angkq
|
private |
Definition at line 77 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ aq0
|
private |
Definition at line 77 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ aqvec
|
private |
Definition at line 77 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ d
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ Ekaon
|
private |
Definition at line 72 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ Elep
|
mutableprivate |
Definition at line 77 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ Enu
|
mutableprivate |
Definition at line 72 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ f
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ Fm1
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ Fm2
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ fpi
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ fXSecIntegrator
|
private |
cross section integrator
Definition at line 52 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Integral(), and LoadConfig().
◆ g
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and LoadConfig().
◆ leptonPDG
|
mutableprivate |
Definition at line 70 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by XSec().
◆ pkvec
|
private |
Definition at line 72 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and XSec().
◆ reactionType
|
private |
Definition at line 70 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by XSec().
◆ Vus
|
private |
Definition at line 63 of file AlamSimoAtharVacasSKPXSec2014.h.
Referenced by Amatrix_NN(), Amatrix_NP(), Amatrix_PP(), and LoadConfig().
The documentation for this class was generated from the following files:
Generated on for GENIEGenerator by
1.14.0