#include <G4GEMProbability.hh>

Inherits G4VEmissionProbability.

Inherited by G4AlphaGEMProbability, G4B10GEMProbability, G4B11GEMProbability, G4B12GEMProbability, G4B13GEMProbability, G4B8GEMProbability, G4Be10GEMProbability, G4Be11GEMProbability, G4Be12GEMProbability, G4Be7GEMProbability, G4Be9GEMProbability, G4C10GEMProbability, G4C11GEMProbability, G4C12GEMProbability, G4C13GEMProbability, G4C14GEMProbability, G4C15GEMProbability, G4C16GEMProbability, G4DeuteronGEMProbability, G4F17GEMProbability, G4F18GEMProbability, G4F19GEMProbability, G4F20GEMProbability, G4F21GEMProbability, G4He3GEMProbability, G4He6GEMProbability, G4He8GEMProbability, G4Li6GEMProbability, G4Li7GEMProbability, G4Li8GEMProbability, G4Li9GEMProbability, G4Mg22GEMProbability, G4Mg23GEMProbability, G4Mg24GEMProbability, G4Mg25GEMProbability, G4Mg26GEMProbability, G4Mg27GEMProbability, G4Mg28GEMProbability, G4N12GEMProbability, G4N13GEMProbability, G4N14GEMProbability, G4N15GEMProbability, G4N16GEMProbability, G4N17GEMProbability, G4Na21GEMProbability, G4Na22GEMProbability, G4Na23GEMProbability, G4Na24GEMProbability, G4Na25GEMProbability, G4Ne18GEMProbability, G4Ne19GEMProbability, G4Ne20GEMProbability, G4Ne21GEMProbability, G4Ne22GEMProbability, G4Ne23GEMProbability, G4Ne24GEMProbability, G4NeutronGEMProbability, G4O14GEMProbability, G4O15GEMProbability, G4O16GEMProbability, G4O17GEMProbability, G4O18GEMProbability, G4O19GEMProbability, G4O20GEMProbability, G4ProtonGEMProbability, and G4TritonGEMProbability.

Collaboration diagram for G4GEMProbability:

[legend]

Public Member Functions
	G4GEMProbability (G4int anA, G4int aZ, G4double aSpin)

virtual	~G4GEMProbability ()

G4double	EmissionProbability (const G4Fragment &fragment, G4double anEnergy)

void	Dump () const

G4int	GetZ_asInt (void) const

G4int	GetA_asInt (void) const

G4double	GetZ (void) const

G4double	GetA (void) const

G4double	GetSpin (void) const

void	SetCoulomBarrier (const G4VCoulombBarrier *aCoulombBarrierStrategy)

G4double	GetCoulombBarrier (const G4Fragment &fragment) const

G4double	CalcAlphaParam (const G4Fragment &) const

G4double	CalcBetaParam (const G4Fragment &) const

Public Member Functions inherited from G4VEmissionProbability
	G4VEmissionProbability ()

virtual	~G4VEmissionProbability ()

void	Initialise ()

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

Protected Attributes
G4double	fPlanck

std::vector< G4double >	ExcitEnergies

std::vector< G4double >	ExcitSpins

std::vector< G4double >	ExcitLifetimes

Protected Attributes inherited from G4VEmissionProbability
G4int	OPTxs

G4bool	useSICB

G4double	LevelDensity

G4Pow *	fG4pow

G4PairingCorrection *	fPairCorr

Detailed Description

Definition at line 54 of file G4GEMProbability.hh.

Constructor & Destructor Documentation

G4GEMProbability::G4GEMProbability	(	G4int	anA,
		G4int	aZ,
		G4double	aSpin
	)

Definition at line 59 of file G4GEMProbability.cc.

                                                                        : 
   theA(anA), theZ(aZ), Spin(aSpin), theCoulombBarrierPtr(0)
 {
   theEvapLDPptr = new G4EvaporationLevelDensityParameter;
   fG4pow = G4Pow::GetInstance(); 
   fPlanck= CLHEP::hbar_Planck*fG4pow->logZ(2);
   fPairCorr = G4PairingCorrection::GetInstance();
 }

Here is the call graph for this function:

G4GEMProbability::~G4GEMProbability ( )

virtual

Definition at line 68 of file G4GEMProbability.cc.

 {
   delete theEvapLDPptr;
 }

Member Function Documentation

G4double G4GEMProbability::CalcAlphaParam ( const G4Fragment & fragment ) const

inline

Definition at line 202 of file G4GEMProbability.hh.

 {
   //JMQ 190709 values according to Furihata's paper (based on notes added 
   //on proof in Dostrovskii's paper)
   G4double res;
   if(GetZ_asInt() == 0) {
     res = 0.76+1.93/fG4pow->Z13(fragment.GetA_asInt()-GetA_asInt());
   } else {
     res = 1.0 + CCoeficient(fragment.GetZ_asInt()-GetZ_asInt());
   }
   return res;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4GEMProbability::CalcBetaParam ( const G4Fragment & fragment ) const

inline

Definition at line 216 of file G4GEMProbability.hh.

 {
   //JMQ 190709 values according to Furihata's paper (based on notes added 
   //on proof in Dostrovskii's paper)
   G4double res;
   if(GetZ_asInt() == 0) {
     res = (1.66/fG4pow->Z23(fragment.GetA_asInt()-GetA_asInt())-0.05)*CLHEP::MeV/
        CalcAlphaParam(fragment);
   } else {
     res = -GetCoulombBarrier(fragment);
   }
   return res;
 }

Here is the call graph for this function:

void G4GEMProbability::Dump ( ) const

Definition at line 261 of file G4GEMProbability.cc.

 {
   G4double mass = G4NucleiProperties::GetNuclearMass(theA, theZ);
   G4double efermi = 0.0;
   if(theA > 1) { 
     efermi = G4NucleiProperties::GetNuclearMass(theA-1, theZ)
       + neutron_mass_c2 - mass;
   }
   G4int nlev = ExcitEnergies.size();
   G4cout << "GEM: List of Excited States for Isotope Z= " 
      << theZ << " A= " << theA << " Nlevels= " << nlev 
      << " Efermi(MeV)= " << efermi
      << G4endl;
   for(G4int i=0; i< nlev; ++i) {
     G4cout << "Z= " << theZ << " A= " << theA 
        << " Mass(GeV)= " << mass/GeV
        << " Eexc(MeV)= " << ExcitEnergies[i]
        << " Time(ns)= " <<  ExcitLifetimes[i]/ns 
        << G4endl;
   }
   G4cout << G4endl;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Implements G4VEmissionProbability.

Definition at line 73 of file G4GEMProbability.cc.

 {
   G4double probability = 0.0;
     
   if (MaximalKineticEnergy > 0.0 && fragment.GetExcitationEnergy() > 0.0) {
     G4double CoulombBarrier = GetCoulombBarrier(fragment);
       
     probability = 
       CalcProbability(fragment,MaximalKineticEnergy,CoulombBarrier);
 
     // Next there is a loop over excited states for this channel 
     // summing probabilities
     size_t nn = ExcitEnergies.size();
     if (0 < nn) {
       G4double SavedSpin = Spin;
       for (size_t i = 0; i <nn; ++i) {
     Spin = ExcitSpins[i];
     // substract excitation energies
     G4double Tmax = MaximalKineticEnergy - ExcitEnergies[i];
     if (Tmax > 0.0) {
       G4double width = CalcProbability(fragment,Tmax,CoulombBarrier);
       //JMQ April 2010 added condition to prevent reported crash
       // update probability
       if (width > 0. && fPlanck < width*ExcitLifetimes[i]) {
         probability += width;
       }
     }
       }
       // Restore Spin
       Spin = SavedSpin;
     }
   }
   //  Normalization = probability;
   return probability;
 }

Here is the call graph for this function:

G4double G4GEMProbability::GetA ( void ) const

inline

Definition at line 151 of file G4GEMProbability.hh.

 { 
   return theA;
 }

G4int G4GEMProbability::GetA_asInt ( void ) const

inline

Definition at line 141 of file G4GEMProbability.hh.

 { 
   return theA;
 }

Here is the caller graph for this function:

G4double G4GEMProbability::GetCoulombBarrier ( const G4Fragment & fragment ) const

inline

Definition at line 173 of file G4GEMProbability.hh.

 {
   G4double res = 0.0;
   if (theCoulombBarrierPtr) {
     G4int Acomp = fragment.GetA_asInt();
     G4int Zcomp = fragment.GetZ_asInt();
     res = theCoulombBarrierPtr->GetCoulombBarrier(Acomp-theA, Zcomp-theZ,
                 fragment.GetExcitationEnergy() -
                 fPairCorr->GetPairingCorrection(Acomp,Zcomp));
   }
   return res;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4GEMProbability::GetSpin ( void ) const

inline

Definition at line 156 of file G4GEMProbability.hh.

 { 
   return Spin; 
 }

G4double G4GEMProbability::GetZ ( void ) const

inline

Definition at line 146 of file G4GEMProbability.hh.

 { 
   return theZ; 
 }

G4int G4GEMProbability::GetZ_asInt ( void ) const

inline

Definition at line 136 of file G4GEMProbability.hh.

 { 
   return theZ; 
 }

Here is the caller graph for this function:

void G4GEMProbability::SetCoulomBarrier ( const G4VCoulombBarrier * aCoulombBarrierStrategy )

inline

Definition at line 167 of file G4GEMProbability.hh.

 {
   theCoulombBarrierPtr = aCoulombBarrierStrategy;
 }

Here is the caller graph for this function:

Member Data Documentation

std::vector<G4double> G4GEMProbability::ExcitEnergies

protected

Definition at line 126 of file G4GEMProbability.hh.

std::vector<G4double> G4GEMProbability::ExcitLifetimes

protected

Definition at line 132 of file G4GEMProbability.hh.

std::vector<G4double> G4GEMProbability::ExcitSpins

protected

Definition at line 129 of file G4GEMProbability.hh.

G4double G4GEMProbability::fPlanck

protected

Definition at line 123 of file G4GEMProbability.hh.

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

geant4.10.03.p01/source/processes/hadronic/models/de_excitation/gem_evaporation/include/G4GEMProbability.hh
geant4.10.03.p01/source/processes/hadronic/models/de_excitation/gem_evaporation/src/G4GEMProbability.cc

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation