Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4PreCompoundNucleon Class Referenceabstract

#include <G4PreCompoundNucleon.hh>

+ Inheritance diagram for G4PreCompoundNucleon:

Public Member Functions

 G4PreCompoundNucleon (const G4ParticleDefinition *, G4VCoulombBarrier *p=nullptr)
 
 ~G4PreCompoundNucleon () override=default
 
 G4PreCompoundNucleon (const G4PreCompoundNucleon &right)=delete
 
const G4PreCompoundNucleonoperator= (const G4PreCompoundNucleon &right)=delete
 
G4bool operator== (const G4PreCompoundNucleon &right) const =delete
 
G4bool operator!= (const G4PreCompoundNucleon &right) const =delete
 
- Public Member Functions inherited from G4PreCompoundFragment
 G4PreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
 
 ~G4PreCompoundFragment () override=default
 
G4double CalcEmissionProbability (const G4Fragment &aFragment) override
 
G4double SampleKineticEnergy (const G4Fragment &aFragment) override
 
 G4PreCompoundFragment (const G4PreCompoundFragment &right)=delete
 
const G4PreCompoundFragmentoperator= (const G4PreCompoundFragment &right)=delete
 
G4bool operator== (const G4PreCompoundFragment &right) const =delete
 
G4bool operator!= (const G4PreCompoundFragment &right) const =delete
 
- Public Member Functions inherited from G4VPreCompoundFragment
 G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *)
 
virtual ~G4VPreCompoundFragment ()
 
void Initialize (const G4Fragment &aFragment)
 
G4bool IsItPossible (const G4Fragment &aFragment) const
 
G4ReactionProductGetReactionProduct () const
 
G4int GetA () const
 
G4int GetZ () const
 
G4int GetRestA () const
 
G4int GetRestZ () const
 
G4double GetBindingEnergy () const
 
G4double GetEnergyThreshold () const
 
G4double GetEmissionProbability () const
 
G4double GetNuclearMass () const
 
G4double GetRestNuclearMass () const
 
const G4LorentzVectorGetMomentum () const
 
void SetMomentum (const G4LorentzVector &lv)
 
void SetOPTxs (G4int opt)
 
void UseSICB (G4bool use)
 
 G4VPreCompoundFragment (const G4VPreCompoundFragment &right)=delete
 
const G4VPreCompoundFragmentoperator= (const G4VPreCompoundFragment &right)=delete
 
G4bool operator== (const G4VPreCompoundFragment &right) const =delete
 
G4bool operator!= (const G4VPreCompoundFragment &right) const =delete
 

Protected Member Functions

G4double ProbabilityDistributionFunction (G4double eKin, const G4Fragment &) override
 
virtual G4double GetRj (G4int NumberParticles, G4int NumberCharged) const =0
 
- Protected Member Functions inherited from G4PreCompoundFragment
G4double CrossSection (G4double ekin) const
 
- Protected Member Functions inherited from G4VPreCompoundFragment
virtual G4double GetAlpha () const =0
 
virtual G4double GetBeta () const
 

Additional Inherited Members

- Protected Attributes inherited from G4VPreCompoundFragment
G4NuclearLevelDatafNucData
 
G4DeexPrecoParameterstheParameters
 
G4Powg4calc
 
G4int theA
 
G4int theZ
 
G4int theResA {0}
 
G4int theResZ {0}
 
G4int theFragA {0}
 
G4int theFragZ {0}
 
G4double theResA13 {0.0}
 
G4double theBindingEnergy {0.0}
 
G4double theMinKinEnergy {0.0}
 
G4double theMaxKinEnergy {0.0}
 
G4double theResMass {0.0}
 
G4double theReducedMass {0.0}
 
G4double theMass
 
G4double theEmissionProbability {0.0}
 
G4double theCoulombBarrier {0.0}
 
G4int OPTxs {3}
 
G4bool useSICB {true}
 

Detailed Description

Definition at line 40 of file G4PreCompoundNucleon.hh.

Constructor & Destructor Documentation

◆ G4PreCompoundNucleon() [1/2]

G4PreCompoundNucleon::G4PreCompoundNucleon ( const G4ParticleDefinition * part,
G4VCoulombBarrier * p = nullptr )

Definition at line 47 of file G4PreCompoundNucleon.cc.

49 : G4PreCompoundFragment(part,aCoulombBarrier)
50{}
G4PreCompoundFragment(const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)

◆ ~G4PreCompoundNucleon()

G4PreCompoundNucleon::~G4PreCompoundNucleon ( )
overridedefault

◆ G4PreCompoundNucleon() [2/2]

G4PreCompoundNucleon::G4PreCompoundNucleon ( const G4PreCompoundNucleon & right)
delete

Member Function Documentation

◆ GetRj()

virtual G4double G4PreCompoundNucleon::GetRj ( G4int NumberParticles,
G4int NumberCharged ) const
protectedpure virtual

◆ operator!=()

G4bool G4PreCompoundNucleon::operator!= ( const G4PreCompoundNucleon & right) const
delete

◆ operator=()

const G4PreCompoundNucleon & G4PreCompoundNucleon::operator= ( const G4PreCompoundNucleon & right)
delete

◆ operator==()

G4bool G4PreCompoundNucleon::operator== ( const G4PreCompoundNucleon & right) const
delete

◆ ProbabilityDistributionFunction()

G4double G4PreCompoundNucleon::ProbabilityDistributionFunction ( G4double eKin,
const G4Fragment & aFragment )
overrideprotectedvirtual

Implements G4PreCompoundFragment.

Definition at line 52 of file G4PreCompoundNucleon.cc.

55{
56 G4double U = aFragment.GetExcitationEnergy();
57 G4int P = aFragment.GetNumberOfParticles();
58 G4int H = aFragment.GetNumberOfHoles();
59 G4int N = P + H;
60
61 static const G4double sixoverpi2 = 6.0/CLHEP::pi2;
63 G4double g1 = sixoverpi2*fNucData->GetLevelDensity(theResZ, theResA, 0.0);
64
65 G4double A0 = (P*P+H*H+P-3*H)/(4.0*g0);
66 G4double A1 = (A0 - 0.5*P)/g1;
67
68 G4double E0 = U - A0;
69 if (E0 <= 0.0) { return 0.0; }
70
71 G4double E1 = U - eKin - theBindingEnergy - A1;
72 if (E1 <= 0.0) { return 0.0; }
73
74 G4double rj = GetRj(P, aFragment.GetNumberOfCharged());
75 G4double xs = CrossSection(eKin);
76
77 if (rj < 0.0 || xs < 0.0) { return 0.0; }
78
79 static const G4double fact = 2*CLHEP::millibarn
80 /(CLHEP::pi2*CLHEP::hbarc*CLHEP::hbarc*CLHEP::hbarc);
81 G4double Probability = fact * theReducedMass * rj * xs * eKin * P * (N-1)
82 * g4calc->powN(g1*E1/(g0*E0),N-2) * g1 / (E0*g0*g0);
83
84 //G4cout << "N=" << N << " g0=" << g0 << " g1=" << g1 << " E0=" << E0
85 // << " E1=" << E1 << " prob=" << Probability << G4endl;
86
87 return Probability;
88}
double G4double
Definition G4Types.hh:83
int G4int
Definition G4Types.hh:85
G4int GetNumberOfParticles() const
G4int GetNumberOfHoles() const
G4double GetExcitationEnergy() const
G4int GetNumberOfCharged() const
G4double GetLevelDensity(G4int Z, G4int A, G4double U)
G4double powN(G4double x, G4int n) const
Definition G4Pow.cc:162
G4double CrossSection(G4double ekin) const
virtual G4double GetRj(G4int NumberParticles, G4int NumberCharged) const =0
#define N
Definition crc32.c:57
G4double Probability(const G4int A, const G4FermiFragment *f1, const G4FermiFragment *f2, const G4double mass, const G4double exc)

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