63 {
64
66
70
77
78
84
89
93
95 M2*=M2;
98
99 double norm=1.0/sqrt(2*M2+4*
m2-4*
m2*
m2/M2);
102
103
109
110 EvtComplex A000,A001,A010,A011,A100,A101,A110,A111,A200,A201,A210,A211;
111 A000 = norm_e*(eps0*L11);
112 A001 = norm_e*(eps0*L12);
113 A010 = norm_e*(eps0*L21);
114 A011 = norm_e*(eps0*L22);
115
116 A100 = norm_e*(eps1*L11);
117 A101 = norm_e*(eps1*L12);
118 A110 = norm_e*(eps1*L21);
119 A111 = norm_e*(eps1*L22);
120
121 A200 = norm_e*(eps2*L11);
122 A201 = norm_e*(eps2*L12);
123 A210 = norm_e*(eps2*L21);
124 A211 = norm_e*(eps2*L22);
125
126 vertex(0,0,0,norm_edm*
cont(p4eps0,l11) + A000);
127 vertex(0,0,1,norm_edm*
cont(p4eps0,l12) + A001);
128 vertex(0,1,0,norm_edm*
cont(p4eps0,l21) + A010);
129 vertex(0,1,1,norm_edm*
cont(p4eps0,l22) + A011);
130
131 vertex(1,0,0,norm_edm*
cont(p4eps1,l11) + A100);
132 vertex(1,0,1,norm_edm*
cont(p4eps1,l12) + A101);
133 vertex(1,1,0,norm_edm*
cont(p4eps1,l21) + A110);
134 vertex(1,1,1,norm_edm*
cont(p4eps1,l22) + A111);
135
136 vertex(2,0,0,norm_edm*
cont(p4eps2,l11) + A200);
137 vertex(2,0,1,norm_edm*
cont(p4eps2,l12) + A201);
138 vertex(2,1,0,norm_edm*
cont(p4eps2,l21) + A210);
139 vertex(2,1,1,norm_edm*
cont(p4eps2,l22) + A211);
140
141 return;
142
143}
double sin(const BesAngle a)
double cos(const BesAngle a)
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
EvtTensor4C EvtLeptonTg5Current(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
EvtVector4C EvtLeptonVCurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
EvtComplex cont(const EvtTensor4C &t1, const EvtTensor4C &t2)
void vertex(const EvtComplex &)
virtual EvtDiracSpinor spParent(int) const
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
virtual EvtVector4C eps(int i) const
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)