63 {
64
66
68 EvtVector4R pv, ps, ppr;
69
75
76
77 EvtHelSys angles( ppr, pv );
78 double theta = angles.getHelAng( 1 );
79 double phi = angles.getHelAng( 2 );
80 double gamma = 0.;
81
82 double A10 = 1.0, A11 = 2.08, A12 = 3.03;
83 double Af10 = A10, Af1f1 = A11, Af1f2 = A12;
84 int j1, j2, j3;
85 for ( j1 = 0; j1 <= 2; )
86 {
87 for ( j2 = 0; j2 <= 1; )
88 {
89 for ( j3 = 0; j3 <= 4; )
90 {
92 j3++;
93 }
94 j2++;
95 }
96 j1++;
97 }
98
99 vertex( 0, 0, 1,
Djmn( 1, 1, -1, phi, theta, gamma ) * A12 );
100 vertex( 0, 0, 3,
Djmn( 1, 1, 0, phi, theta, gamma ) * A11 );
101 vertex( 0, 0, 0,
Djmn( 1, 1, 1, phi, theta, gamma ) * A10 );
102 vertex( 0, 1, 0,
Djmn( 1, 1, -1, phi, theta, gamma ) * Af10 );
103 vertex( 0, 1, 4,
Djmn( 1, 1, 0, phi, theta, gamma ) * Af1f1 );
104 vertex( 0, 1, 2,
Djmn( 1, 1, 1, phi, theta, gamma ) * Af1f2 );
105
106 vertex( 1, 0, 1,
Djmn( 1, -1, -1, phi, theta, gamma ) * A12 );
107 vertex( 1, 0, 3,
Djmn( 1, -1, 0, phi, theta, gamma ) * A11 );
108 vertex( 1, 0, 0,
Djmn( 1, -1, 1, phi, theta, gamma ) * A10 );
109 vertex( 1, 1, 0,
Djmn( 1, -1, -1, phi, theta, gamma ) * Af10 );
110 vertex( 1, 1, 4,
Djmn( 1, -1, 0, phi, theta, gamma ) * Af1f1 );
111 vertex( 1, 1, 2,
Djmn( 1, -1, 1, phi, theta, gamma ) * Af1f2 );
112
113 vertex( 2, 0, 1,
Djmn( 1, 0, -1, phi, theta, gamma ) * A12 );
114 vertex( 2, 0, 3,
Djmn( 1, 0, 0, phi, theta, gamma ) * A11 );
115 vertex( 2, 0, 0,
Djmn( 1, 0, 1, phi, theta, gamma ) * A10 );
116 vertex( 2, 1, 0,
Djmn( 1, 0, -1, phi, theta, gamma ) * Af10 );
117 vertex( 2, 1, 4,
Djmn( 1, 0, 0, phi, theta, gamma ) * Af1f1 );
118 vertex( 2, 1, 2,
Djmn( 1, 0, 1, phi, theta, gamma ) * Af1f2 );
119
120 return;
121}
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
void vertex(const EvtComplex &)
const EvtVector4R & getP4() const
EvtParticle * getDaug(int i)
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)