Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4PreCompoundNucleon Class Referenceabstract

#include <G4PreCompoundNucleon.hh>

Inheritance diagram for G4PreCompoundNucleon:
Collaboration diagram for G4PreCompoundNucleon:

Public Member Functions

 G4PreCompoundNucleon (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
 
virtual ~G4PreCompoundNucleon ()
 
- Public Member Functions inherited from G4PreCompoundFragment
 G4PreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
 
virtual ~G4PreCompoundFragment ()
 
G4double CalcEmissionProbability (const G4Fragment &aFragment)
 
G4double SampleKineticEnergy (const G4Fragment &aFragment)
 
- Public Member Functions inherited from G4VPreCompoundFragment
 G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
 
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 &value)
 
void SetOPTxs (G4int)
 
void UseSICB (G4bool)
 

Protected Member Functions

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

Additional Inherited Members

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

Detailed Description

Definition at line 41 of file G4PreCompoundNucleon.hh.

Constructor & Destructor Documentation

G4PreCompoundNucleon::G4PreCompoundNucleon ( const G4ParticleDefinition part,
G4VCoulombBarrier aCoulombBarrier 
)

Definition at line 47 of file G4PreCompoundNucleon.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4PreCompoundNucleon::~G4PreCompoundNucleon ( )
virtual

Definition at line 52 of file G4PreCompoundNucleon.cc.

53 {}

Member Function Documentation

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

Implemented in G4PreCompoundNeutron, and G4PreCompoundProton.

Here is the caller graph for this function:

G4double G4PreCompoundNucleon::ProbabilityDistributionFunction ( G4double  eKin,
const G4Fragment aFragment 
)
protectedvirtual

Implements G4PreCompoundFragment.

Definition at line 56 of file G4PreCompoundNucleon.cc.

58 {
59  G4double U = aFragment.GetExcitationEnergy();
60  G4int P = aFragment.GetNumberOfParticles();
61  G4int H = aFragment.GetNumberOfHoles();
62  G4int N = P + H;
63 
64  static const G4double sixoverpi2 = 6.0/CLHEP::pi2;
66  G4double g1 = sixoverpi2*theResA*theParameters->GetLevelDensity();
67 
68  G4double A0 = G4double(P*P+H*H+P-3*H)/(4.0*g0);
69  G4double A1 = (A0 - 0.5*P)/g1;
70 
71  G4double E0 = U - A0;
72  if (E0 <= 0.0) { return 0.0; }
73 
74  G4double E1 = U - eKin - theBindingEnergy - A1;
75  if (E1 <= 0.0) { return 0.0; }
76 
77  G4double rj = GetRj(P, aFragment.GetNumberOfCharged());
78  G4double xs = CrossSection(eKin);
79 
80  if (rj <0.0 || xs < 0.0) { return 0.0; }
81 
82  static const G4double fact = 2*CLHEP::millibarn
84  G4double Probability = fact * theReducedMass * rj * xs * eKin * P * (N-1)
85  * g4calc->powN(g1*E1/(g0*E0),N-2) * g1 / (E0*g0*g0);
86 
87  return Probability;
88 }
const int N
Definition: mixmax.h:43
G4double powN(G4double x, G4int n) const
Definition: G4Pow.cc:128
static constexpr double pi2
Definition: SystemOfUnits.h:57
static constexpr double millibarn
Definition: SystemOfUnits.h:86
static constexpr double hbarc
G4int GetNumberOfParticles() const
Definition: G4Fragment.hh:345
G4double GetLevelDensity() const
int G4int
Definition: G4Types.hh:78
static double P[]
G4DeexPrecoParameters * theParameters
G4int GetNumberOfHoles() const
Definition: G4Fragment.hh:365
G4double CrossSection(G4double ekin) const
virtual G4double GetRj(G4int NumberParticles, G4int NumberCharged) const =0
double G4double
Definition: G4Types.hh:76
G4int GetNumberOfCharged() const
Definition: G4Fragment.hh:350
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:283

Here is the call graph for this function:


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