104 G4double Theta = std::acos(cosTheta);
107 G4ThreeVector pFinal1(std::sin(Theta)*std::cos(phi), std::sin(Theta)*std::sin(phi), cosTheta);
109 G4double pInCM = std::sqrt((S-(m10+m20)*(m10+m20))*(S-(m10-m20)*(m10-m20))/(4.*S));
110 pFinal1 = pFinal1 * pInCM;
120 p4Final1 *= toLabFrame;
121 p4Final2 *= toLabFrame;
130 finalTracks->push_back(final1);
131 finalTracks->push_back(final2);