G4PreCompoundIon Class Reference

#include <G4PreCompoundIon.hh>

Inheritance diagram for G4PreCompoundIon:

Public Member Functions
	G4PreCompoundIon (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
	~G4PreCompoundIon () override=default
	G4PreCompoundIon (const G4PreCompoundIon &right)=delete
const G4PreCompoundIon &	operator= (const G4PreCompoundIon &right)=delete
G4bool	operator== (const G4PreCompoundIon &right) const =delete
G4bool	operator!= (const G4PreCompoundIon &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 G4PreCompoundFragment &	operator= (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
G4ReactionProduct *	GetReactionProduct () 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 G4LorentzVector &	GetMomentum () const
void	SetMomentum (const G4LorentzVector &lv)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)
	G4VPreCompoundFragment (const G4VPreCompoundFragment &right)=delete
const G4VPreCompoundFragment &	operator= (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
virtual G4double	FactorialFactor (G4int N, G4int P) const =0
virtual G4double	CoalescenceFactor (G4int A) 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
G4NuclearLevelData *	fNucData
G4DeexPrecoParameters *	theParameters
G4Pow *	g4calc
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 41 of file G4PreCompoundIon.hh.

Constructor & Destructor Documentation

G4PreCompoundIon() [1/2]

G4PreCompoundIon::G4PreCompoundIon	(	const G4ParticleDefinition *	part,
		G4VCoulombBarrier *	aCoulombBarrier )

Definition at line 47 of file G4PreCompoundIon.cc.

  : G4PreCompoundFragment(part,aCoulombBarrier)
{
  G4double r0 = theParameters->GetR0();
  fact = 0.75*CLHEP::millibarn/(CLHEP::pi*r0*r0*r0);
}

~G4PreCompoundIon()

G4PreCompoundIon::~G4PreCompoundIon ( )

overridedefault

G4PreCompoundIon() [2/2]

G4PreCompoundIon::G4PreCompoundIon ( const G4PreCompoundIon & right )

delete

Member Function Documentation

CoalescenceFactor()

virtual G4double G4PreCompoundIon::CoalescenceFactor ( G4int A ) const

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

FactorialFactor()

virtual G4double G4PreCompoundIon::FactorialFactor ( G4int N, G4int P ) const

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

GetRj()

virtual G4double G4PreCompoundIon::GetRj ( G4int NumberParticles, G4int NumberCharged ) const

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

operator!=()

G4bool G4PreCompoundIon::operator!= ( const G4PreCompoundIon & right ) const

delete

operator=()

const G4PreCompoundIon & G4PreCompoundIon::operator= ( const G4PreCompoundIon & right )

delete

operator==()

G4bool G4PreCompoundIon::operator== ( const G4PreCompoundIon & right ) const

delete

ProbabilityDistributionFunction()

G4double G4PreCompoundIon::ProbabilityDistributionFunction ( G4double eKin, const G4Fragment & aFragment )

overrideprotectedvirtual

Implements G4PreCompoundFragment.

Definition at line 56 of file G4PreCompoundIon.cc.

{
  G4double efinal = eKin + theBindingEnergy;
  if(efinal <= 0.0 ) { return 0.0; } 
 
  G4double U = aFragment.GetExcitationEnergy();
  G4int P = aFragment.GetNumberOfParticles();
  G4int H = aFragment.GetNumberOfHoles();
  G4int A = GetA();
  G4int N = P + H;
 
  static const G4double sixoverpi2 = 6.0/CLHEP::pi2;
  G4double g0 = sixoverpi2*fNucData->GetLevelDensity(theFragZ, theFragA, U);
  G4double g1 = sixoverpi2*fNucData->GetLevelDensity(theResZ, theResA, 0.0);
 
  G4double gj = g1;
 
  G4double A0 = (P*P+H*H+P-3*H)/(4.0*g0);
  G4double A1 = std::max(0.0,(A0*g0 + A*(A-2*P-1)*0.25)/g1); 
 
  G4double E0 = U - A0;
  if (E0 <= 0.0) { return 0.0; }
 
  G4double E1 = std::max(0.0,theMaxKinEnergy - eKin - A1); 
 
  G4double Aj = A*(A+1)/(4.0*gj); 
  G4double Ej = std::max(0.0,efinal - Aj); 
 
  G4double rj = GetRj(P, aFragment.GetNumberOfCharged());
  G4double xs = CrossSection(eKin);
  
  G4double pA = fact*eKin*xs*rj 
    * CoalescenceFactor(theFragA) * FactorialFactor(N,P)
    * std::sqrt(2.0/(theReducedMass*efinal)) 
    * g4calc->powN(g1*E1/(g0*E0), N-A-1)
    * g4calc->powN(gj*Ej/(g0*E0), A-1)*gj*g1/(g0*g0*E0*theResA); 
   
  return pA;
}

---

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

• G4PreCompoundIon.hh
• G4PreCompoundIon.cc