G4HETCFragment Class Reference

#include <G4HETCFragment.hh>

Inheritance diagram for G4HETCFragment:

Collaboration diagram for G4HETCFragment:

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

Protected Member Functions
virtual G4double	K (const G4Fragment &aFragment)=0
virtual G4double	GetSpinFactor () const =0
virtual G4double	GetAlpha () const =0
virtual G4double	GetBeta () const =0
G4double	BetaRand (G4int N, G4int L) const

Private Member Functions
G4double	IntegrateEmissionProbability (G4double &Low, G4double &Up, const G4Fragment &aFragment)
	G4HETCFragment ()
	G4HETCFragment (const G4HETCFragment &right)
const G4HETCFragment &	operator= (const G4HETCFragment &right)
G4int	operator== (const G4HETCFragment &right) const
G4int	operator!= (const G4HETCFragment &right) const

Private Attributes
G4double	r2norm

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
G4PreCompoundParameters *	theParameters
G4Pow *	g4pow
G4int	theA
G4int	theZ
G4int	theResA
G4int	theResZ
G4int	theFragA
G4int	theFragZ
G4double	theResA13
G4double	theBindingEnergy
G4double	theMaxKinEnergy
G4double	theResMass
G4double	theReducedMass
G4double	theMass
G4double	theEmissionProbability
G4double	theCoulombBarrier
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 40 of file G4HETCFragment.hh.

Constructor & Destructor Documentation

G4HETCFragment() [1/3]

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

Definition at line 39 of file G4HETCFragment.cc.

  : G4VPreCompoundFragment(part, aCoulombBarrier)
{
  G4double r0 = theParameters->Getr0();
  r2norm = r0*r0/(CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc*CLHEP::hbarc);
}

Here is the call graph for this function:

~G4HETCFragment()

G4HETCFragment::~G4HETCFragment ( )

virtual

Definition at line 47 of file G4HETCFragment.cc.

{}

Here is the call graph for this function:

G4HETCFragment() [2/3]

G4HETCFragment::G4HETCFragment ( )

private

G4HETCFragment() [3/3]

G4HETCFragment::G4HETCFragment ( const G4HETCFragment &  right )

private

Member Function Documentation

BetaRand()

G4double G4HETCFragment::BetaRand ( G4int  N, G4int  L  ) const

inlineprotected

Definition at line 79 of file G4HETCFragment.hh.

{
  G4double Y1 = G4RandGamma::shoot(N,1);
  G4double Y2 = G4RandGamma::shoot(L,1);
  
  return Y1/(Y1+Y2);
}

Here is the call graph for this function:

Here is the caller graph for this function:

CalcEmissionProbability()

G4double G4HETCFragment::CalcEmissionProbability ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 51 of file G4HETCFragment.cc.

{
  if (GetEnergyThreshold() <= 0.0) 
    {
      theEmissionProbability = 0.0;
      return 0.0;
    }    
  // Coulomb barrier is the lower limit 
  // of integration over kinetic energy
  theEmissionProbability = 
    IntegrateEmissionProbability(theCoulombBarrier,theMaxKinEnergy,aFragment);
    
  return theEmissionProbability;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetAlpha()

virtual G4double G4HETCFragment::GetAlpha ( ) const

protectedpure virtual

Implements G4VPreCompoundFragment.

Implemented in G4HETCNeutron, G4HETCDeuteron, G4HETCAlpha, G4HETCProton, and G4HETCTriton.

Here is the caller graph for this function:

GetBeta()

virtual G4double G4HETCFragment::GetBeta ( ) const

protectedpure virtual

Implements G4VPreCompoundFragment.

Implemented in G4HETCNeutron, G4HETCDeuteron, G4HETCAlpha, G4HETCProton, and G4HETCTriton.

Here is the caller graph for this function:

GetSpinFactor()

virtual G4double G4HETCFragment::GetSpinFactor ( ) const

protectedpure virtual

Implemented in G4HETCNeutron, G4HETCDeuteron, G4HETCAlpha, G4HETCProton, and G4HETCTriton.

Here is the caller graph for this function:

IntegrateEmissionProbability()

G4double G4HETCFragment::IntegrateEmissionProbability ( G4double &  Low, G4double &  Up, const G4Fragment &  aFragment  )

private

Definition at line 67 of file G4HETCFragment.cc.

{    
  G4double U = aFragment.GetExcitationEnergy();

  G4int P  = aFragment.GetNumberOfParticles();
  G4int H  = aFragment.GetNumberOfHoles();
  G4int N  = P + H;
  G4int Pb = P - theA;
  G4int Nb = Pb + H;
  if (Nb <= 0.0) { return 0.0; }
  G4double ga = (6.0/pi2)*theFragA*theParameters->GetLevelDensity();
  G4double gb = (6.0/pi2)*theResA*theParameters->GetLevelDensity();

  G4double A  = G4double(P*P+H*H+P-3*H)/(4.0*ga);
  G4double Ab = G4double(Pb*Pb+H*H+Pb-3*H)/(4.0*gb);
  U = std::max(U-A,0.0);
  if (U <= 0.0) { return 0.0; }

  G4int Pf = P;
  G4int Hf = H;
  G4int Nf = N-1;
  for (G4int i = 1; i < theA; ++i)
    {
      Pf *= (P-i);
      Hf *= (H-i);
      Nf *= (N-1-i);
    }

  G4double X = std::max(Up - Ab + GetBeta(),0.0);
  G4double Y = std::max(Up - Ab - Low, 0.0);

  G4double Probability = r2norm*GetSpinFactor()*theReducedMass*GetAlpha() 
    *g4pow->Z23(theResA)*Pf*Hf*Nf*K(aFragment)*(X/Nb - Y/(Nb+1))
    *U*g4pow->powN(gb*Y,Nb)/g4pow->powN(ga*U,N);

  return Probability;
}

Here is the call graph for this function:

Here is the caller graph for this function:

K()

virtual G4double G4HETCFragment::K ( const G4Fragment &  aFragment )

protectedpure virtual

Implemented in G4HETCNeutron, G4HETCDeuteron, G4HETCAlpha, G4HETCProton, and G4HETCTriton.

Here is the caller graph for this function:

operator!=()

G4int G4HETCFragment::operator!= ( const G4HETCFragment &  right ) const

private

operator=()

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

private

operator==()

G4int G4HETCFragment::operator== ( const G4HETCFragment &  right ) const

private

Member Data Documentation

r2norm

G4double G4HETCFragment::r2norm

private

Definition at line 76 of file G4HETCFragment.hh.

---

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

• G4HETCFragment.hh
• G4HETCFragment.cc