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 )
static constexpr double MeV
G4double GetExcitationEnergyOfExcitedIsomer(G4int, G4int, G4int)
static G4IonTable * GetIonTable()
void SetMomentumChange(const G4ThreeVector &aV)
void SetMomentum(const G4double x, const G4double y, const G4double z)
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
G4double GetPDGMass() const
void SetTotalEnergy(const G4double en)
G4LENDManager * lend_manager
void SetMass(const G4double mas)
void SetEnergyChange(G4double anEnergy)
G4double GetTotalEnergy() const
G4int GetNucleusEncoding(G4int iZ, G4int iA, G4int iM)
G4double GetKineticEnergy() const
static G4Pow * GetInstance()
const G4Material * GetMaterial() const
G4double powA(G4double A, G4double y) const
const G4Isotope * GetIsotope()
G4HadFinalState * returnUnchanged(const G4HadProjectile &aTrack, G4HadFinalState *theResult)
double getElasticFinalState(double e_in, double temperature, double(*rng)(void *), void *rngState)
static constexpr double twopi
G4ErrorTarget * theTarget
G4double GetTotalMomentum() const
G4ThreeVector GetMomentum() const
ThreeVector shoot(const G4int Ap, const G4int Af)
const G4LorentzVector & Get4Momentum() const
G4double GetKineticEnergy() const
const G4ParticleDefinition * GetDefinition() const
static constexpr double k_Boltzmann
void SetMomentum(const G4ThreeVector &momentum)
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static G4LENDManager * GetInstance()
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4GIDI_target * get_target_from_map(G4int nuclear_code)
G4HadFinalState theParticleChange
G4double GetTemperature() const