34 : sym(symmetrize)
35{
36
37
38 mSigma = 0.55;
39 cmSigma = 1.20;
40 gSigma = 9.4;
41
42 mOmega = 0.783;
43 cmOmega = 0.808;
44 gOmega = 10.95;
45
46 mPion = 0.138;
47 cmPion = 0.51;
48 gPion = 7.27;
49
50 mNucleon = 0.938;
51
52
53
54 m42 = 4. * mNucleon * mNucleon;
55 mPion2 = mPion * mPion;
56 cmPion2 = cmPion * cmPion;
57 dPion1 = cmPion2-mPion2;
58 dPion2 = dPion1 * dPion1;
59 cm6gp = 1.5 * (cmPion2*cmPion2*cmPion2) * (gPion*gPion*gPion*gPion) * m42 * m42 / dPion2;
60
61 cPion_3 = -(cm6gp/3.);
62 cPion_2 = -(cm6gp * mPion2/dPion1);
63 cPion_1 = -(cm6gp * mPion2 * (2. * cmPion2 + mPion2) / dPion2);
64 cPion_m = -(cm6gp * cmPion2 * mPion2 / dPion2);
65 cPion_L = -(cm6gp * 2. * cmPion2 * mPion2 * (cmPion2 + mPion2) / dPion2 / dPion1);
66 cPion_0 = -(cPion_3 + cPion_2 + cPion_1 + cPion_m);
67
68
69
71
72 mSigma2 = mSigma * mSigma;
73 cmSigma2 = cmSigma * cmSigma;
74 cmSigma4 = cmSigma2 * cmSigma2;
75 cmSigma6 = cmSigma2 * cmSigma4;
76 dSigma1 = m42 - cmSigma2;
77 dSigma2 = m42 - mSigma2;
78 dSigma3 = cmSigma2 - mSigma2;
79
80 G4double dSigma1Sq = dSigma1 * dSigma1;
81 G4double dSigma2Sq = dSigma2 * dSigma2;
82 G4double dSigma3Sq = dSigma3 * dSigma3;
83
84 cm2gs = 0.5 * cmSigma2 * gSigmaSq*gSigmaSq / dSigma3Sq;
85
86
87 cSigma_3 = -(cm2gs * dSigma1Sq / 3.);
88 cSigma_2 = -(cm2gs * cmSigma2 * dSigma1 * dSigma2 / dSigma3);
89 cSigma_1 = -(cm2gs * cmSigma4 * (2. * dSigma1 + dSigma2) * dSigma2 / dSigma3Sq);
90 cSigma_m = -(cm2gs * cmSigma6 * dSigma2Sq / mSigma2 / dSigma3Sq);
91 cSigma_L = -(cm2gs * cmSigma6 * dSigma2 * (dSigma1 + dSigma2) * 2. / (dSigma3 * dSigma3Sq));
92 cSigma_0 = -(cSigma_3 + cSigma_2 + cSigma_1 + cSigma_m);
93
94
95
97
98 mOmega2 = mOmega * mOmega;
99 cmOmega2 = cmOmega * cmOmega;
100 cmOmega4 = cmOmega2 * cmOmega2;
101 cmOmega6 = cmOmega2 * cmOmega4;
102 dOmega1 = m42 - cmOmega2;
103 dOmega2 = m42 - mOmega2;
104 dOmega3 = cmOmega2 - mOmega2;
105 sOmega1 = cmOmega2 + mOmega2;
106
107 G4double dOmega3Sq = dOmega3 * dOmega3;
108
109 cm2go = 0.5 * cmOmega2 * gOmegaSq * gOmegaSq / dOmega3Sq;
110
111 cOmega_3 = cm2go / 3.;
112 cOmega_2 = -(cm2go * cmOmega2 / dOmega3);
113 cOmega_1 = cm2go * cmOmega4 / dOmega3Sq;
114 cOmega_m = cm2go * cmOmega6 / (dOmega3Sq * mOmega2);
115 cOmega_L = -(cm2go * cmOmega6 * 4. / (dOmega3 * dOmega3Sq));
116
117
118
119 G4double fac1Tmp = (gSigma * gOmega * cmSigma2 * cmOmega2);
120 fac1 = -(fac1Tmp * fac1Tmp * m42);
121 dMix1 = cmOmega2 - cmSigma2;
122 dMix2 = cmOmega2 - mSigma2;
123 dMix3 = cmSigma2 - mOmega2;
124
128
129 cMix_o1 = fac1 / (cmOmega2 * dMix1Sq * dMix2 * dOmega3);
130 cMix_s1 = fac1 / (cmSigma2 * dMix1Sq * dMix3 * dSigma3);
131 cMix_Omega = fac1 / (dOmega3Sq * dMix3Sq * (mOmega2 - mSigma2));
132 cMix_sm = fac1 / (dSigma3Sq * dMix2Sq * (mSigma2 - mOmega2));
133 fac2 = (-fac1) / (dMix1*dMix1Sq * dOmega3Sq * dMix2Sq);
134 fac3 = (-fac1) / (dMix1*dMix1Sq * dSigma3Sq * dMix3Sq);
135
136 cMix_oLc = fac2 * (3. * cmOmega2*cmOmega4 - cmOmega4 * cmSigma2
137 - 2. * cmOmega4 * mOmega2 - 2. * cmOmega4 * mSigma2
138 + cmOmega2 * mOmega2 * mSigma2 + cmSigma2 * mOmega2 * mSigma2
139 - 4. * cmOmega4 * m42 + 2. * cmOmega2 * cmSigma2 * m42
140 + 3. * cmOmega2 * mOmega2 * m42 - cmSigma2 * mOmega2 * m42
141 + 3. * cmOmega2 * mSigma2 * m42 - cmSigma2 * mSigma2 * m42
142 - 2. * mOmega2 * mSigma2 * m42);
143
144 cMix_oLs = fac2 * (8. * cmOmega4 - 4. * cmOmega2 * cmSigma2
145 - 6. * cmOmega2 * mOmega2 + 2. * cmSigma2 * mOmega2
146 - 6. * cmOmega2 * mSigma2 + 2. * cmSigma2 * mSigma2
147 + 4. * mOmega2 * mSigma2);
148
149 cMix_sLc = fac3 * (cmOmega2 * cmSigma4 - 3. * cmSigma6
150 + 2. * cmSigma4 * mOmega2 + 2. * cmSigma4 * mSigma2
151 - cmOmega2 * mOmega2 * mSigma2 - cmSigma2 * mOmega2 * mSigma2
152 - 2. * cmOmega2 * cmSigma2 * m42 + 4. * cmSigma4 * m42
153 + cmOmega2 * mOmega2 * m42 - 3. * cmSigma2 * mOmega2 * m42
154 + cmOmega2 * mSigma2 * m42 - 3. * cmSigma2 * mSigma2 * m42
155 + 2. * mOmega2 * mSigma2 * m42);
156
157 cMix_sLs = fac3 * (4. * cmOmega2 * cmSigma2 - 8. * cmSigma4
158 - 2. * cmOmega2 * mOmega2 + 6. * cmSigma2 * mOmega2
159 - 2. * cmOmega2 * mSigma2 + 6. * cmSigma2 * mSigma2
160 - 4. * mOmega2 * mSigma2);
161}