59 {
60
63
64 if ( first )
65 {
67
68 float val = 0.6;
70 char vak[20];
71 sprintf( vak,
"PARJ(13)=%f", val );
72 std::string temp( vak );
73 lugive_( temp.c_str(), strlen( temp.c_str() ) );
74 }
75 EvtVector4R p4[100];
76
78
79 int flavor;
80
81 int i, more;
82 int ndaugjs;
83 int kf[100];
84 EvtId id[100];
86
88
89 double px[100], py[100], pz[100], e[100];
90
91 if ( p->
getNDaug() != 0 ) {
return; }
92 do {
93
95
96 for ( i = 0; i < ndaugjs; i++ )
97 {
98
100
102
103
104
105
106
107 if ( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] >= e[i] * e[i] )
108 { e[i] = sqrt( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] ) + 0.0000000000001; }
109
110 p4[i].
set( e[i], px[i], py[i], pz[i] );
111 }
112
114
115 more = ( ( channel != -1 ) && ( channel != p->
getChannel() ) );
116
117 } while ( more );
118
120
121 for ( i = 0; i < ndaugjs; i++ ) { p->
getDaug( i )->
init(
id[i], p4[i] ); }
122 return;
123}
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)