59 {
60
61
64
66 EvtId other_b;
67
69
70 EvtParticle *pip, *pim, *pi0;
71
73
74
78
79 EvtVector4R p4[3];
80
83 int iset;
84
86
87 if ( first == 1 )
88 {
89 iset = 10000;
91 }
92 else { iset = 0; }
93
94 double p4piplus[4], p4piminus[4], p4gamm1[4], p4gamm2[4];
95
96 double realA, imgA, realbarA, imgbarA;
97
98 evt3pions_( &
alpha, &iset, p4piplus, p4piminus, p4gamm1, p4gamm2, &realA, &imgA, &realbarA,
99 &imgbarA );
100
101 p4[0].
set( p4piplus[3], p4piplus[0], p4piplus[1], p4piplus[2] );
102 p4[1].
set( p4piminus[3], p4piminus[0], p4piminus[1], p4piminus[2] );
103 p4[2].
set( p4gamm1[3] + p4gamm2[3], p4gamm1[0] + p4gamm2[0], p4gamm1[1] + p4gamm2[1],
104 p4gamm1[2] + p4gamm2[2] );
105
107 {
110 }
111 else
112 {
115 }
116
118
119 EvtComplex amp;
120
121 EvtComplex
A( realA, imgA );
122 EvtComplex Abar( realbarA, imgbarA );
123
124 if ( other_b == B0B )
125 {
128 }
129 if ( other_b == B0 )
130 {
133 }
134
136
137 return;
138}
void evt3pions_(double *, int *, double *, double *, double *, double *, double *, double *, double *, double *)
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)
void makeDaughters(int ndaug, EvtId *id)
virtual void init(EvtId part_n, const EvtVector4R &p4)=0
EvtParticle * getDaug(int i)
void set(int i, double d)