void parentMomentum(double w,double y,double &E,double &px,double &py,double&pz);
double pp(double E);
- void twoBodyDecay(starlightConstants::particleTypeEnum &ipid,double E,double W,double px0,double py0,double pz0,double &px1,double &py1,double&pz1,double &px2,double &py2,/*double &py2,*/double &pz2,int &iFbadevent);
+ void twoBodyDecay(starlightConstants::particleTypeEnum &ipid,double E,double W,
+ double px0,double py0,double pz0,
+ double &px1,double &py1,double &pz1, double &pt1,
+ double &px2,double &py2,double &pz2, double &pt2,
+ double &mdec,
+ int &iFbadevent);
// void transform(double betax,double betay,double betaz,double &E,double &px,double &py,double &pz,int &iFbadevent);
void thephi(double W,double px,double py,double pz,double E,double &theta,double &phi);