G4PreCompoundIon Class Reference

#include <G4PreCompoundIon.hh>

Inheritance diagram for G4PreCompoundIon:

Public Member Functions
	G4PreCompoundIon (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4PreCompoundIon ()
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)

---

Detailed Description

Definition at line 42 of file G4PreCompoundIon.hh.

---

Constructor & Destructor Documentation

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

Definition at line 47 of file G4PreCompoundIon.cc.

References G4PreCompoundParameters::Getr0(), G4INCL::Math::pi, and G4VPreCompoundFragment::theParameters.

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

G4PreCompoundIon::~G4PreCompoundIon

(

)

[virtual]

Definition at line 55 of file G4PreCompoundIon.cc.

{}

---

Member Function Documentation

virtual G4double G4PreCompoundIon::CoalescenceFactor

(

G4int

A

)

[protected, pure virtual]

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

virtual G4double G4PreCompoundIon::CrossSection

(

G4double

ekin

)

[protected, pure virtual]

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

virtual G4double G4PreCompoundIon::FactorialFactor	(	G4int	N,
		G4int	P
	)			[protected, pure virtual]

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

virtual G4double G4PreCompoundIon::GetAlpha

(

)

[protected, pure virtual]

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by GetOpt0().

G4double G4PreCompoundIon::GetBeta

(

)

[inline, protected]

Definition at line 86 of file G4PreCompoundIon.hh.

References G4VPreCompoundFragment::GetCoulombBarrier().

Referenced by GetOpt0().

{
  return -GetCoulombBarrier();
}

G4double G4PreCompoundIon::GetOpt0

(

G4double

ekin

)

[inline, protected]

Definition at line 92 of file G4PreCompoundIon.hh.

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

Referenced by G4PreCompoundTriton::CrossSection(), G4PreCompoundHe3::CrossSection(), G4PreCompoundDeuteron::CrossSection(), and G4PreCompoundAlpha::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 G4PreCompoundIon::GetRj	(	G4int	NumberParticles,
		G4int	NumberCharged
	)			[protected, pure virtual]

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, and G4PreCompoundTriton.

Referenced by ProbabilityDistributionFunction().

G4double G4PreCompoundIon::ProbabilityDistributionFunction	(	G4double	eKin,
		const G4Fragment &	aFragment
	)			[protected, virtual]

Implements G4PreCompoundFragment.

Definition at line 59 of file G4PreCompoundIon.cc.

References CoalescenceFactor(), CrossSection(), FactorialFactor(), G4VPreCompoundFragment::g4pow, G4VPreCompoundFragment::GetA(), G4Fragment::GetA_asInt(), G4VPreCompoundFragment::GetBindingEnergy(), G4Fragment::GetExcitationEnergy(), G4PreCompoundParameters::GetLevelDensity(), G4VPreCompoundFragment::GetMaximalKineticEnergy(), G4Fragment::GetNumberOfCharged(), G4Fragment::GetNumberOfHoles(), G4Fragment::GetNumberOfParticles(), G4VPreCompoundFragment::GetReducedMass(), G4VPreCompoundFragment::GetRestA(), GetRj(), G4VPreCompoundFragment::IsItPossible(), G4Pow::powN(), and G4VPreCompoundFragment::theParameters.

{
  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

---