Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ParticleHPThermalBoost Class Reference

#include <G4ParticleHPThermalBoost.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 42 of file G4ParticleHPThermalBoost.hh.

Member Function Documentation

G4double G4ParticleHPThermalBoost::GetThermalEnergy ( const G4HadProjectile aP,
const G4Element anE,
G4double  aT 
)
inline

Definition at line 45 of file G4ParticleHPThermalBoost.hh.

48  {
49  G4double theA = anE->GetN();
50  G4double theZ = anE->GetZ();
51  return GetThermalEnergy(aP, theA ,theZ, aT);
52  }
G4double GetN() const
Definition: G4Element.hh:135
G4double GetZ() const
Definition: G4Element.hh:131
G4double GetThermalEnergy(const G4HadProjectile &aP, const G4Element *anE, G4double aT)
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 54 of file G4ParticleHPThermalBoost.hh.

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

Here is the call graph for this function:


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