#include <G4PreCompoundNucleon.hh>

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

Public Member Functions inherited from G4VPreCompoundFragment
	G4VPreCompoundFragment (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4VPreCompoundFragment ()

void	Initialize (const G4Fragment &aFragment)

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	ProbabilityDistributionFunction (G4double eKin, const G4Fragment &aFragment)

virtual G4double	CrossSection (G4double ekin)=0

virtual G4double	GetRj (G4int NumberParticles, G4int NumberCharged)=0

virtual G4double	GetAlpha ()=0

virtual G4double	GetBeta ()=0

G4double	GetOpt0 (G4double ekin)

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 41 of file G4PreCompoundNucleon.hh.

Constructor & Destructor Documentation

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

Definition at line 48 of file G4PreCompoundNucleon.cc.

   : G4PreCompoundFragment(part,aCoulombBarrier)
 {
   fact = 2*CLHEP::millibarn/(CLHEP::pi2*CLHEP::hbarc*CLHEP::hbarc*CLHEP::hbarc);
 }

G4PreCompoundNucleon::~G4PreCompoundNucleon ( )

virtual

Definition at line 55 of file G4PreCompoundNucleon.cc.

56 {}

Member Function Documentation

virtual G4double G4PreCompoundNucleon::CrossSection ( G4double ekin )

protectedpure virtual

Implemented in G4PreCompoundNeutron, and G4PreCompoundProton.

Referenced by ProbabilityDistributionFunction().

virtual G4double G4PreCompoundNucleon::GetAlpha ( )

protectedpure virtual

Implemented in G4PreCompoundNeutron, and G4PreCompoundProton.

Referenced by GetOpt0().

virtual G4double G4PreCompoundNucleon::GetBeta ( )

protectedpure virtual

Implemented in G4PreCompoundNeutron, and G4PreCompoundProton.

Referenced by GetOpt0().

G4double G4PreCompoundNucleon::GetOpt0 ( G4double ekin )

inlineprotected

Definition at line 82 of file G4PreCompoundNucleon.hh.

References GetAlpha(), GetBeta(), G4PreCompoundParameters::Getr0(), G4VPreCompoundFragment::ResidualA13(), and G4VPreCompoundFragment::theParameters.

Referenced by G4PreCompoundNeutron::CrossSection(), and G4PreCompoundProton::CrossSection().

 {
   G4double r0 = theParameters->Getr0()*ResidualA13();
   // cross section is now given in mb (r0 is in mm) for the sake of consistency
   //with the rest of the options
   return 1.e+25*CLHEP::pi*r0*r0*ResidualA13()*GetAlpha()*(1.+GetBeta()/K);
 }

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

protectedpure virtual

Implemented in G4PreCompoundNeutron, and G4PreCompoundProton.

Referenced by ProbabilityDistributionFunction().

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

protectedvirtual

Implements G4PreCompoundFragment.

Definition at line 59 of file G4PreCompoundNucleon.cc.

 {
   if ( !IsItPossible(aFragment) ) { return 0.0; }
 
   G4double U = aFragment.GetExcitationEnergy();
   G4int P = aFragment.GetNumberOfParticles();
   G4int H = aFragment.GetNumberOfHoles();
   G4int N = P + H;
 
   G4double g0 = (6.0/pi2)*aFragment.GetA_asInt()*theParameters->GetLevelDensity();
   G4double g1 = (6.0/pi2)*GetRestA()*theParameters->GetLevelDensity();
   
   G4double A0 = G4double(P*P+H*H+P-3*H)/(4.0*g0);
   G4double A1 = (A0 - 0.5*P)/g1;  
 
   G4double E0 = U - A0;
   if (E0 <= 0.0) { return 0.0; }
 
   G4double E1 = U - eKin - GetBindingEnergy() - A1;
   if (E1 <= 0.0) { return 0.0; }
 
   G4double rj = GetRj(P, aFragment.GetNumberOfCharged());
   G4double xs = CrossSection(eKin);
 
   if (rj <0.0 || xs < 0.0) {  
     std::ostringstream errOs;
     G4cout<<"WARNING:  NEGATIVE VALUES "<<G4endl;     
     errOs << "Rj=" << rj <<"  xsec("
       <<eKin/MeV<<" MeV)= "<< xs <<"  A= "<<GetA()<<"  Z= "<<GetZ()
       <<G4endl;
     throw G4HadronicException(__FILE__, __LINE__, errOs.str());
   }
 
   G4double Probability = fact * GetReducedMass() * rj * xs * eKin * P * (N-1) 
     * g4pow->powN(g1*E1/(g0*E0),N-2) * g1 / (E0*g0*g0);
   
   /*
   G4double Probability = 1.0/pi2*2.0/(hbarc*hbarc*hbarc) * GetReducedMass()
     * GetRj(aFragment.GetNumberOfParticles(), aFragment.GetNumberOfCharged())  
     * eKin*CrossSection(eKin)*millibarn * P*(N-1.0) * 
    std::pow(g1*E1/(g0*E0),N-2.0)/E0 * g1/(g0*g0);
   */
 
   return Probability;
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation