74{
77
79
80 int i,more;
81 int ndaugjs;
82 int kf[100];
85
86 double px[100],py[100],pz[100],e[100];
87
89 do{
91
92 for(i=0;i<ndaugjs;i++)
93 {
94
96
98
99
100
101
102
103 if (px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i]>=e[i]*e[i])
104 e[i]=sqrt(px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i])+0.0000000000001;
105
106 p4[i].set(e[i],px[i],py[i],pz[i]);
107
108 }
109
111
112 more=((channel!=-1)&&(channel!=p->
getChannel()));
113
114 }while(more);
115
117
118 for(i=0;i<ndaugjs;i++)
120
121 return ;
122}
************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 *)
double double double * p4