#include <G4HETCFragment.hh>

Inheritance 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	GetKineticEnergy (const G4Fragment &aFragment)=0

G4ReactionProduct *	GetReactionProduct () const

G4int	GetA () const

G4int	GetZ () const

G4int	GetRestA () const

G4int	GetRestZ () const

G4double	ResidualA13 () const

G4double	GetCoulombBarrier () const

G4double	GetBindingEnergy () const

G4double	GetMaximalKineticEnergy () const

G4double	GetEnergyThreshold () const

G4double	GetEmissionProbability () const

G4double	GetNuclearMass () const

G4double	GetRestNuclearMass () const

G4double	GetReducedMass () const

const G4LorentzVector &	GetMomentum () const

void	SetMomentum (const G4LorentzVector &value)

const G4String	GetName () const

void	SetOPTxs (G4int)

void	UseSICB (G4bool)

Protected Member Functions
virtual G4double	K (const G4Fragment &aFragment)=0

virtual G4double	GetSpinFactor ()=0

virtual G4double	GetAlpha ()=0

virtual G4double	GetBeta ()=0

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

Protected Member Functions inherited from G4VPreCompoundFragment
G4bool	IsItPossible (const G4Fragment &aFragment) const

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
G4PreCompoundParameters *	theParameters

G4Pow *	g4pow

G4double	theEmissionProbability

G4double	theCoulombBarrier

G4int	OPTxs

G4bool	useSICB

Detailed Description

Definition at line 40 of file G4HETCFragment.hh.

Constructor & Destructor Documentation

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

Definition at line 39 of file G4HETCFragment.cc.

References G4PreCompoundParameters::Getr0(), and G4VPreCompoundFragment::theParameters.

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

G4HETCFragment::~G4HETCFragment ( )

virtual

Definition at line 47 of file G4HETCFragment.cc.

48 {}

Member Function Documentation

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

inlineprotected

Definition at line 80 of file G4HETCFragment.hh.

References G4INCL::DeJongSpin::shoot().

Referenced by G4HETCChargedFragment::GetKineticEnergy(), and G4HETCNeutron::GetKineticEnergy().

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

G4double G4HETCFragment::CalcEmissionProbability ( const G4Fragment & aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 51 of file G4HETCFragment.cc.

References G4VPreCompoundFragment::GetEnergyThreshold(), G4VPreCompoundFragment::GetMaximalKineticEnergy(), G4VPreCompoundFragment::theCoulombBarrier, and G4VPreCompoundFragment::theEmissionProbability.

 {
   if (GetEnergyThreshold() <= 0.0) 
     {
       theEmissionProbability = 0.0;
       return 0.0;
     }    
   // Coulomb barrier is the lower limit 
   // of integration over kinetic energy
   G4double LowerLimit = theCoulombBarrier;
   
   // Excitation energy of nucleus after fragment emission is the upper limit
   // of integration over kinetic energy
   G4double UpperLimit = GetMaximalKineticEnergy();
   
   theEmissionProbability = 
     IntegrateEmissionProbability(LowerLimit,UpperLimit,aFragment);
     
   return theEmissionProbability;
 }