60 {
61
63
65 EvtVector4R pv, ps, ppr;
66
72
73
74 EvtHelSys angles( ppr, pv );
75 double theta = angles.getHelAng( 1 );
76 double phi = angles.getHelAng( 2 );
77 double gamma = 0;
78
81 double F10 = F01;
82 double F11 = F00;
83
84 vertex( 0, 0, 0,
Djmn( 2, 0, 0, phi, theta, gamma ) * F00 );
85
86 vertex( 0, 0, 1,
Djmn( 2, 0, 1, phi, theta, gamma ) * F01 );
87 vertex( 0, 1, 0,
Djmn( 2, 0, -1, phi, theta, gamma ) * F10 );
88 vertex( 0, 1, 1,
Djmn( 2, 0, 0, phi, theta, gamma ) * F11 );
89
90 vertex( 1, 0, 0,
Djmn( 2, 2, 0, phi, theta, gamma ) * F00 );
91
92 vertex( 1, 0, 1,
Djmn( 2, 2, 1, phi, theta, gamma ) * F01 );
93 vertex( 1, 1, 0,
Djmn( 2, 2, -1, phi, theta, gamma ) * F10 );
94 vertex( 1, 1, 1,
Djmn( 2, 2, 0, phi, theta, gamma ) * F11 );
95
96 vertex( 2, 0, 0,
Djmn( 2, -2, 0, phi, theta, gamma ) * F00 );
97
98 vertex( 2, 0, 1,
Djmn( 2, -2, 1, phi, theta, gamma ) * F01 );
99 vertex( 2, 1, 0,
Djmn( 2, -2, -1, phi, theta, gamma ) * F10 );
100 vertex( 2, 1, 1,
Djmn( 2, -2, 0, phi, theta, gamma ) * F11 );
101
102 vertex( 3, 0, 0,
Djmn( 2, 1, 0, phi, theta, gamma ) * F00 );
103
104 vertex( 3, 0, 1,
Djmn( 2, 1, 1, phi, theta, gamma ) * F01 );
105 vertex( 3, 1, 0,
Djmn( 2, 1, -1, phi, theta, gamma ) * F10 );
106 vertex( 3, 1, 1,
Djmn( 2, 1, 0, phi, theta, gamma ) * F11 );
107
108 vertex( 4, 0, 0,
Djmn( 2, -1, 0, phi, theta, gamma ) * F00 );
109
110 vertex( 4, 0, 1,
Djmn( 2, -1, 1, phi, theta, gamma ) * F01 );
111 vertex( 4, 1, 0,
Djmn( 2, -1, -1, phi, theta, gamma ) * F10 );
112 vertex( 4, 1, 1,
Djmn( 2, -1, 0, phi, theta, gamma ) * F11 );
113
114 return;
115}
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)