162 {
163
164
167
169 EvtId other_b;
170
172
178
179
180
181
182 EvtComplex amp;
183
185
186 EvtComplex A_a1p, Abar_a1p, A_a2p, Abar_a2p;
187 EvtComplex A_a1m, Abar_a1m, A_a2m, Abar_a2m;
188
191
194
196 Abar_a1m =
198
200 Abar_a2m =
202
203 EvtComplex a2p_amp =
EvtAmpA2( mom1, mom2, mom3, mom4 ) +
204 EvtAmpA2( mom1, mom4, mom3, mom2 ) +
206
207 EvtComplex a2m_amp =
EvtAmpA2( mom2, mom3, mom4, mom1 ) +
208 EvtAmpA2( mom2, mom1, mom4, mom3 ) +
210
211 EvtComplex a1p_amp =
EvtAmpA1( mom1, mom2, mom3, mom4 ) +
212 EvtAmpA1( mom1, mom4, mom3, mom2 ) +
214
215 EvtComplex a1m_amp =
EvtAmpA1( mom2, mom3, mom4, mom1 ) +
216 EvtAmpA1( mom2, mom1, mom4, mom3 ) +
218
219 A = A_a2p * a2p_amp + A_a1p * a1p_amp + A_a2m * a2m_amp + A_a1m * a1m_amp;
220 Abar = Abar_a2p * a2p_amp + Abar_a1p * a1p_amp + Abar_a2m * a2m_amp + Abar_a1m * a1m_amp;
221
222 if ( other_b == B0B )
223 {
227 }
228 if ( other_b == B0 )
229 {
233 }
234
236
237 return;
238}
EvtComplex EvtAmpA2(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
EvtComplex EvtAmpA1(const EvtVector4R &p4pi1, const EvtVector4R &p4pi2, const EvtVector4R &p4pi3, const EvtVector4R &p4pi4)
double sin(const BesAngle a)
double cos(const BesAngle a)
static void OtherB(EvtParticle *p, double &t, EvtId &otherb)
void vertex(const EvtComplex &)
static EvtId getId(const std::string &name)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)