60 {
62 EvtVector4R p4[100];
63
65
66 int i, more;
67 int ndaugjs;
68 int kf[100];
69 EvtId id[100];
71
72 double px[100], py[100], pz[100], e[100];
73
74 if ( p->
getNDaug() != 0 ) {
return; }
75 do {
77
78 for ( i = 0; i < ndaugjs; i++ )
79 {
80
82
84
85
86
87
88
89 if ( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] >= e[i] * e[i] )
90 e[i] = sqrt( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] ) + 0.0000000000001;
91
92 p4[i].
set( e[i], px[i], py[i], pz[i] );
93 }
94
96
97 more = ( ( channel != -1 ) && ( channel != p->
getChannel() ) );
98
99 } while ( more );
100
102
103 for ( i = 0; i < ndaugjs; i++ ) p->
getDaug( i )->
init(
id[i], p4[i] );
104
105 return;
106}
************Class m_ypar INTEGER m_KeyWgt INTEGER m_KeyIHVP INTEGER m_KeyGPS INTEGER m_IsBeamPolarized INTEGER m_EvtGenInterface DOUBLE PRECISION m_Emin DOUBLE PRECISION m_sphot DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_q2 DOUBLE PRECISION m_PolBeam2 DOUBLE PRECISION m_xErrPb *COMMON c_KK2f $ !CMS energy average $ !Spin Polarization vector first beam $ !Spin Polarization vector second beam $ !Beam energy spread[GeV] $ !minimum hadronization energy[GeV] $ !input READ never touch them !$ !debug facility $ !maximum weight $ !inverse alfaQED $ !minimum real photon energy
static int inChannelList(EvtId parent, int ndaug, EvtId *daugs)
static EvtId evtIdFromStdHep(int stdhep)
static EvtSpinType::spintype getSpinType(EvtId i)
void makeDaughters(int ndaug, EvtId *id)
virtual void init(EvtId part_n, const EvtVector4R &p4)=0
EvtParticle * getDaug(int i)
static void pythiaInit(int f)
static void pythiacont(double *, int *, int *, double *, double *, double *, double *)
void set(int i, double d)