G4EvaporationProbability Class Reference

#include <G4EvaporationProbability.hh>

Inheritance diagram for G4EvaporationProbability:

Public Member Functions
	G4EvaporationProbability (G4int anA, G4int aZ, G4double aGamma, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4EvaporationProbability ()
G4int	GetZ (void) const
G4int	GetA (void) const
G4double	ProbabilityDistributionFunction (const G4Fragment &aFragment, G4double K)
G4double	EmissionProbability (const G4Fragment &fragment, G4double anEnergy)
Public Member Functions inherited from G4VEmissionProbability
	G4VEmissionProbability ()
virtual	~G4VEmissionProbability ()
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Protected Member Functions
	G4EvaporationProbability ()
virtual G4double	CrossSection (const G4Fragment &fragment, G4double K)=0
virtual G4double	CalcAlphaParam (const G4Fragment &fragment)=0
virtual G4double	CalcBetaParam (const G4Fragment &fragment)=0

Additional Inherited Members
Protected Attributes inherited from G4VEmissionProbability
G4int	OPTxs
G4bool	useSICB
G4Pow *	fG4pow
G4PairingCorrection *	fPairCorr
G4EvaporationLevelDensityParameter *	theEvapLDPptr

Detailed Description

Definition at line 44 of file G4EvaporationProbability.hh.

Constructor & Destructor Documentation

G4EvaporationProbability::G4EvaporationProbability	(	G4int	anA,
		G4int	aZ,
		G4double	aGamma,
		G4VCoulombBarrier *	aCoulombBarrier
	)

Definition at line 49 of file G4EvaporationProbability.cc.

  : theA(anA),
    theZ(aZ),
    Gamma(aGamma),
    theCoulombBarrierptr(aCoulombBarrier) 
{}

G4EvaporationProbability::~G4EvaporationProbability ( )

virtual

Definition at line 65 of file G4EvaporationProbability.cc.

{}

G4EvaporationProbability::G4EvaporationProbability ( )

protected

Definition at line 58 of file G4EvaporationProbability.cc.

  : theA(0),
    theZ(0),
    Gamma(0.0),
    theCoulombBarrierptr(0) 
{}

Member Function Documentation

virtual G4double G4EvaporationProbability::CalcAlphaParam ( const G4Fragment &  fragment )

protectedpure virtual

virtual G4double G4EvaporationProbability::CalcBetaParam ( const G4Fragment &  fragment )

protectedpure virtual

virtual G4double G4EvaporationProbability::CrossSection ( const G4Fragment &  fragment, G4double  K  )

protectedpure virtual

Referenced by ProbabilityDistributionFunction().

G4double G4EvaporationProbability::EmissionProbability ( const G4Fragment &  fragment, G4double  anEnergy  )

virtual

Implements G4VEmissionProbability.

Definition at line 69 of file G4EvaporationProbability.cc.

References G4Fragment::GetExcitationEnergy().

{
  G4double probability = 0.0;

  if (anEnergy > 0.0 && fragment.GetExcitationEnergy() > 0.0) {
    probability = CalculateProbability(fragment, anEnergy);

  }
  return probability;
}

G4int G4EvaporationProbability::GetA ( void  ) const

inline

Definition at line 55 of file G4EvaporationProbability.hh.

{ return theA;} 

G4int G4EvaporationProbability::GetZ ( void  ) const

inline

Definition at line 53 of file G4EvaporationProbability.hh.

{ return theZ; }

G4double G4EvaporationProbability::ProbabilityDistributionFunction	(	const G4Fragment &	aFragment,
		G4double	K
	)

Definition at line 210 of file G4EvaporationProbability.cc.

References G4Fragment::ComputeGroundStateMass(), CrossSection(), G4VEmissionProbability::fPairCorr, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4PairingCorrection::GetPairingCorrection(), G4Fragment::GetZ_asInt(), python.hepunit::hbarc, G4EvaporationLevelDensityParameter::LevelDensityParameter(), python.hepunit::millibarn, python.hepunit::pi, and G4VEmissionProbability::theEvapLDPptr.

{ 
  G4int ResidualA = fragment.GetA_asInt() - theA;
  G4int ResidualZ = fragment.GetZ_asInt() - theZ;  
  G4double U = fragment.GetExcitationEnergy();
  //G4cout << "### G4EvaporationProbability::ProbabilityDistributionFunction" << G4endl;
  //G4cout << "FragZ= " << fragment.GetZ_asInt() << " FragA= " << fragment.GetA_asInt()
  //     << " Z= " << theZ << "  A= " << theA << G4endl;
  //G4cout << "PC " << fPairCorr << "   DP " << theEvapLDPptr << G4endl;

  // if(K <= theCoulombBarrierptr->GetCoulombBarrier(ResidualA,ResidualZ,U)) return 0.0;

  G4double delta1 = fPairCorr->GetPairingCorrection(ResidualA,ResidualZ);
 
  G4double delta0 = fPairCorr->GetPairingCorrection(fragment.GetA_asInt(),
                            fragment.GetZ_asInt());

  
  G4double ParticleMass = fragment.ComputeGroundStateMass(theZ,theA);
  G4double ResidualMass = fragment.ComputeGroundStateMass(ResidualZ,ResidualA);

  G4double theSeparationEnergy = ParticleMass + ResidualMass 
    - fragment.GetGroundStateMass();

  G4double a0 = theEvapLDPptr->LevelDensityParameter(fragment.GetA_asInt(),
                             fragment.GetZ_asInt(),
                             U - delta0);

  G4double a1 = theEvapLDPptr->LevelDensityParameter(ResidualA, ResidualZ,
                             U - theSeparationEnergy - delta1);
  
  
  G4double E0 = U - delta0;

  G4double E1 = U - theSeparationEnergy - delta1 - K;

  if (E1<0.) { return 0.; }

  //JMQ 14/02/09 BUG fixed: hbarc should be in the denominator instead of hbar_Planck 
  //Without 1/hbar_Panck remains as a width

  //G4double Prob=Gamma*ParticleMass/((pi*hbarc)*(pi*hbarc)*std::exp(2*std::sqrt(a0*E0)))
  //  *K*CrossSection(fragment,K)*std::exp(2*std::sqrt(a1*E1))*millibarn;

  static const G4double pcoeff = millibarn/((pi*hbarc)*(pi*hbarc)); 

  // Fixed numerical problem
  G4double Prob = pcoeff*Gamma*ParticleMass*std::exp(2*(std::sqrt(a1*E1) - std::sqrt(a0*E0)))
    *K*CrossSection(fragment,K);

  return Prob;
}

---

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

• G4EvaporationProbability.hh
• G4EvaporationProbability.cc