68 theNeutron.SetKineticEnergy( ke );
89 G4double nEnergy = theNeutron.GetTotalEnergy();
96 G4double cmsMom = std::sqrt(the3CMS*the3CMS);
97 G4double sqrts = std::sqrt((totE-cmsMom)*(totE+cmsMom));
101 theNeutron.Lorentz(theNeutron, theCMS);
102 theTarget.
Lorentz(theTarget, theCMS);
103 G4double en = theNeutron.GetTotalMomentum();
108 tempVector.
setX( std::cos(theta)*std::sin(cms_theta)*std::cos(cms_phi)
109 +std::sin(theta)*std::cos(phi)*std::cos(cms_theta)*std::cos(cms_phi)
110 -std::sin(theta)*std::sin(phi)*std::sin(cms_phi) );
111 tempVector.
setY( std::cos(theta)*std::sin(cms_theta)*std::sin(cms_phi)
112 +std::sin(theta)*std::cos(phi)*std::cos(cms_theta)*std::sin(cms_phi)
113 +std::sin(theta)*std::sin(phi)*std::cos(cms_phi) );
114 tempVector.
setZ( std::cos(theta)*std::cos(cms_theta)
115 -std::sin(theta)*std::cos(phi)*std::sin(cms_theta) );
117 theNeutron.SetMomentum(tempVector);
124 theNeutron.Lorentz(theNeutron, -1.*theCMS);
125 theTarget.
Lorentz(theTarget, -1.*theCMS);
128 if ( theNeutron.GetKineticEnergy() <= 0 )