76 {
77
80
81 if (first){
83
84 float val=0.6;
87 }
88 char vak[20];
90 std::string temp(vak);
91 lugive_(temp.c_str(),strlen(temp.c_str()));
92 }
94
96
97 int flavor;
98
99 int i,more;
100 int ndaugjs;
101 int kf[100];
104
106
107 double px[100],py[100],pz[100],e[100];
108
109 if ( p->
getNDaug() != 0 ) {
return;}
110 do{
111
113
114 for(i=0;i<ndaugjs;i++){
115
117
119
120
121
122
123
124 if (px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i]>=e[i]*e[i]){
125
126 e[i]=sqrt(px[i]*px[i]+py[i]*py[i]+pz[i]*pz[i])+0.0000000000001;
127
128 }
129
130 p4[i].
set(e[i],px[i],py[i],pz[i]);
131
132 }
133
135
136 more=((channel!=-1)&&(channel!=p->
getChannel()));
137
138
139 }while(more);
140
142
143 for(i=0;i<ndaugjs;i++){
145 }
146 return ;
147}
sprintf(cut,"kal_costheta0_em>-0.93&&kal_costheta0_em<0.93&&kal_pxy0_em>=0.05+%d*0.1&&kal_pxy0_em<0.15+%d*0.1&&NGch>=2", j, j)
void lugive_(const char *cnfgstr, int length)
void continuum_(double *, int *, int *, int *, double *, double *, double *, double *)
************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)
void set(int i, double d)