G4PreCompoundIon Class Reference

#include <G4PreCompoundIon.hh>

Inheritance diagram for G4PreCompoundIon:

Public Member Functions
	G4PreCompoundIon (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4PreCompoundIon ()
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)
virtual G4double	CalcEmissionProbability (const G4Fragment &aFragment)=0
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	ProbabilityDistributionFunction (G4double eKin, const G4Fragment &aFragment)
virtual G4double	CrossSection (G4double ekin)=0
virtual G4double	GetRj (G4int NumberParticles, G4int NumberCharged)=0
virtual G4double	FactorialFactor (G4int N, G4int P)=0
virtual G4double	CoalescenceFactor (G4int A)=0
virtual G4double	GetAlpha ()=0
G4double	GetBeta ()
G4double	GetOpt0 (G4double ekin)
virtual G4double	ProbabilityDistributionFunction (G4double K, const G4Fragment &aFragment)=0
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 42 of file G4PreCompoundIon.hh.

Constructor & Destructor Documentation

G4PreCompoundIon()

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

Definition at line 46 of file G4PreCompoundIon.cc.

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

~G4PreCompoundIon()

G4PreCompoundIon::~G4PreCompoundIon ( )

virtual

Definition at line 55 of file G4PreCompoundIon.cc.

{}

Member Function Documentation

CoalescenceFactor()

virtual G4double G4PreCompoundIon::CoalescenceFactor ( G4int  A )

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

CrossSection()

virtual G4double G4PreCompoundIon::CrossSection ( G4double  ekin )

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

FactorialFactor()

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

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

GetAlpha()

virtual G4double G4PreCompoundIon::GetAlpha ( )

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by GetOpt0().

GetBeta()

G4double G4PreCompoundIon::GetBeta ( )

inlineprotected

Definition at line 86 of file G4PreCompoundIon.hh.

{
  return -GetCoulombBarrier();
}

Referenced by GetOpt0().

GetOpt0()

G4double G4PreCompoundIon::GetOpt0 ( G4double  ekin )

inlineprotected

Definition at line 92 of file G4PreCompoundIon.hh.

{
  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);
}

Referenced by G4PreCompoundAlpha::CrossSection(), G4PreCompoundDeuteron::CrossSection(), G4PreCompoundHe3::CrossSection(), and G4PreCompoundTriton::CrossSection().

GetRj()

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

protectedpure virtual

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

ProbabilityDistributionFunction()

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

protectedvirtual

Implements G4PreCompoundFragment.

Definition at line 58 of file G4PreCompoundIon.cc.

{
  if ( !IsItPossible(aFragment) ) { return 0.0; }
  G4double efinal = eKin + GetBindingEnergy();
  //G4cout << "Efinal= " << efinal << " Ekin= " << eKin << G4endl;
  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;
 
  G4double g0 = (6.0/pi2)*aFragment.GetA_asInt()*theParameters->GetLevelDensity();
  G4double g1 = (6.0/pi2)*GetRestA()*theParameters->GetLevelDensity();
 
  //JMQ 06/02/209  This is  THE BUG that was killing cluster emission
  // G4double gj = (6.0/pi2)*GetA() *
  //   G4PreCompoundParameters::GetAddress()->GetLevelDensity();
 
  G4double gj = g1;
 
  G4double A0 = G4double(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;
  //G4cout << "E0= " << E0 << G4endl;
  if (E0 <= 0.0) { return 0.0; }
 
  G4double E1 = (std::max(0.0,GetMaximalKineticEnergy() - 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);
 
  //G4cout << "rj= " << rj << " xs= " << xs << G4endl;
 
  // JMQ 10/02/09 reshaping of the formula (unnecessary std::pow elimitated)
  /*
  G4double r0 = theParameters->Getr0();
  G4double pA = (3.0/4.0) * std::sqrt(std::max(0.0, 2.0/(GetReducedMass()*
        (eKin+GetBindingEnergy()))))/(pi * r0 * r0 *r0* GetRestA())* 
                eKin*CrossSection(eKin)*millibarn* 
                CoalescenceFactor(aFragment.GetA_asInt()) * FactorialFactor(N,P)* 
       GetRj(aFragment.GetNumberOfParticles(), aFragment.GetNumberOfCharged());
 
  G4double pB = std::pow((g1*E1)/(g0*E0),N-GetA()-1.0)*(g1/g0);
  G4double pC = std::pow((gj*Ej)/(g0*E0),GetA()-1.0)*(gj/g0)/E0; 
  pA *= pB * pC;
  */
  
  G4double pA = fact*eKin*xs*rj 
    * CoalescenceFactor(aFragment.GetA_asInt()) * FactorialFactor(N,P)
    * std::sqrt(2.0/(GetReducedMass()*efinal)) 
    * g4pow->powN(g1*E1/(g0*E0), N-A-1)
    * g4pow->powN(gj*Ej/(g0*E0), A-1)*gj*g1/(g0*g0*E0*GetRestA()); 
   
  return pA;
}

---

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

• G4PreCompoundIon.hh
• G4PreCompoundIon.cc