Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4NeutronHPThermalBoost Class Reference

#include <G4NeutronHPThermalBoost.hh>

Public Member Functions

G4double GetThermalEnergy (const G4HadProjectile &aP, const G4Element *anE, G4double aT)
 
G4double GetThermalEnergy (const G4HadProjectile &aP, G4double theA, G4double theZ, G4double aT)
 

Detailed Description

Definition at line 39 of file G4NeutronHPThermalBoost.hh.

Member Function Documentation

◆ GetThermalEnergy() [1/2]

◆ GetThermalEnergy() [2/2]

G4double G4NeutronHPThermalBoost::GetThermalEnergy ( const G4HadProjectile aP,
G4double  theA,
G4double  theZ,
G4double  aT 
)
inline

Definition at line 51 of file G4NeutronHPThermalBoost.hh.

54 {
55 // prepare neutron
56 G4double eKinetic = aP.GetKineticEnergy();
57 G4ReactionProduct theNeutron( const_cast<G4ParticleDefinition *>(aP.GetDefinition()) );
58 theNeutron.SetMomentum( aP.Get4Momentum().vect() );
59 theNeutron.SetKineticEnergy( eKinetic );
60 G4ThreeVector neuVelo = (1./aP.GetDefinition()->GetPDGMass())*theNeutron.GetMomentum();
61
62 // prepare properly biased thermal nucleus
63 G4Nucleus aNuc;
64 G4double eps = 0.0001;
65 G4double eleMass;
66 eleMass = ( G4NucleiProperties::GetNuclearMass( static_cast<G4int>(theA+eps) , static_cast<G4int>(theZ+eps) ) ) / G4Neutron::Neutron()->GetPDGMass();
67
68 G4ReactionProduct aThermalNuc = aNuc.GetBiasedThermalNucleus(eleMass, neuVelo, aT);
69
70 // boost to rest system and return
71 G4ReactionProduct boosted;
72 boosted.Lorentz(theNeutron, aThermalNuc);
73 return boosted.GetKineticEnergy();
74 }
int G4int
Definition: G4Types.hh:66
Hep3Vector vect() const
const G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
const G4LorentzVector & Get4Momentum() const
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4ReactionProduct GetBiasedThermalNucleus(G4double aMass, G4ThreeVector aVelocity, G4double temp=-1) const
Definition: G4Nucleus.cc:108
G4double GetKineticEnergy() const
void Lorentz(const G4ReactionProduct &p1, const G4ReactionProduct &p2)

The documentation for this class was generated from the following file: