G4GEMProbability Class Reference

#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:

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	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
G4Pow *	fG4pow
G4PairingCorrection *	fPairCorr
G4EvaporationLevelDensityParameter *	theEvapLDPptr

Private Member Functions
G4double	CalcProbability (const G4Fragment &fragment, G4double MaximalKineticEnergy, G4double V)
G4double	CCoeficient (G4int) const
G4double	I0 (G4double t)
G4double	I1 (G4double t, G4double tx)
G4double	I2 (G4double s0, G4double sx)
G4double	I3 (G4double s0, G4double sx)
	G4GEMProbability ()
	G4GEMProbability (const G4GEMProbability &right)
const G4GEMProbability &	operator= (const G4GEMProbability &right)
G4bool	operator== (const G4GEMProbability &right) const
G4bool	operator!= (const G4GEMProbability &right) const

Private Attributes
G4Pow *	fG4pow
G4PairingCorrection *	fPairCorr
G4VLevelDensityParameter *	theEvapLDPptr
G4int	theA
G4int	theZ
G4double	Spin
const G4VCoulombBarrier *	theCoulombBarrierPtr

Detailed Description

Definition at line 54 of file G4GEMProbability.hh.

Constructor & Destructor Documentation

G4GEMProbability() [1/3]

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::~G4GEMProbability ( )

virtual

Definition at line 68 of file G4GEMProbability.cc.

{
  delete theEvapLDPptr;
}

G4GEMProbability() [2/3]

G4GEMProbability::G4GEMProbability ( )

private

G4GEMProbability() [3/3]

G4GEMProbability::G4GEMProbability ( const G4GEMProbability &  right )

private

Member Function Documentation

CalcAlphaParam()

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:

CalcBetaParam()

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:

Here is the caller graph for this function:

CalcProbability()

G4double G4GEMProbability::CalcProbability ( const G4Fragment &  fragment, G4double  MaximalKineticEnergy, G4double  V  )

private

Definition at line 110 of file G4GEMProbability.cc.

{
  G4int A = fragment.GetA_asInt();
  G4int Z = fragment.GetZ_asInt();

  G4int ResidualA = A - theA;
  G4int ResidualZ = Z - theZ;
  G4double U = fragment.GetExcitationEnergy();
  
  G4double NuclearMass = fragment.ComputeGroundStateMass(theZ, theA);
  
  G4double Alpha = CalcAlphaParam(fragment);
  G4double Beta = CalcBetaParam(fragment);
    
  //                       ***RESIDUAL***
  //JMQ (September 2009) the following quantities refer to the RESIDUAL:
  
  G4double delta0 = fPairCorr->GetPairingCorrection(ResidualA, ResidualZ);  
  
  G4double a = theEvapLDPptr->
    LevelDensityParameter(ResidualA,ResidualZ,MaximalKineticEnergy+V-delta0);
  G4double Ux = (2.5 + 150.0/G4double(ResidualA))*MeV;
  G4double Ex = Ux + delta0;
  G4double T  = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
  //JMQ fixed bug in units
  G4double E0 = Ex - T*(G4Log(T/MeV) - G4Log(a*MeV)/4.0 
    - 1.25*G4Log(Ux/MeV) + 2.0*std::sqrt(a*Ux));
  //                      ***end RESIDUAL ***
  //                       ***PARENT***
  //JMQ (September 2009) the following quantities refer to the PARENT:
     
  G4double deltaCN = fPairCorr->GetPairingCorrection(A, Z); 
  G4double aCN     = theEvapLDPptr->LevelDensityParameter(A, Z, U-deltaCN);
  G4double UxCN    = (2.5 + 150.0/G4double(A))*MeV;
  G4double ExCN    = UxCN + deltaCN;
  G4double TCN     = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
  //                     ***end PARENT***

  G4double Width;
  G4double InitialLevelDensity;
  G4double t = MaximalKineticEnergy/T;
  if ( MaximalKineticEnergy < Ex ) {
    //JMQ 190709 bug in I1 fixed (T was  missing)
    Width = (I1(t,t)*T + (Beta+V)*I0(t))/G4Exp(E0/T);
    //JMQ 160909 fix:  InitialLevelDensity has been taken away 
    //(different conditions for initial CN..) 
  } else {

    //VI minor speedup
    G4double expE0T = G4Exp(E0/T);
    static const G4double sqrt2 = std::sqrt(2.0);

    G4double tx = Ex/T;
    G4double s0 = 2.0*std::sqrt(a*(MaximalKineticEnergy-delta0));
    G4double sx = 2.0*std::sqrt(a*(Ex-delta0));
    // VI: protection against FPE exception
    if(s0 > 350.) { s0 = 350.; }
    Width = I1(t,tx)*T/expE0T + I3(s0,sx)*G4Exp(s0)/(sqrt2*a);

    // VI this cannot happen!
    // For charged particles (Beta+V) = 0 beacuse Beta = -V
    //if (theZ == 0) {
    //  Width += (Beta+V)*(I0(tx)/expE0T + 2.0*sqrt2*I2(s0,sx)*G4Exp(s0));
    //}
  }
  
  //JMQ 14/07/2009 BIG BUG : NuclearMass is in MeV => hbarc instead of 
  //                         hbar_planck must be used
  //    G4double g = (2.0*Spin+1.0)*NuclearMass/(pi2* hbar_Planck*hbar_Planck);
  G4double gg = (2.0*Spin+1.0)*NuclearMass/(pi2* hbarc*hbarc);
  
  //JMQ 190709 fix on Rb and  geometrical cross sections according to 
  //           Furihata's paper (JAERI-Data/Code 2001-105, p6)
  G4double Rb = 0.0;
  if (theA > 4) 
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      G4double Aj = fG4pow->Z13(theA);
      Rb = 1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85;
      Rb *= fermi;
    }
  else if (theA>1)
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      G4double Aj = fG4pow->Z13(theA);
      Rb=1.5*(Aj+Ad)*fermi;
    }
  else
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      Rb = 1.5*Ad*fermi;
    }
  G4double GeometricalXS = pi*Rb*Rb; 
  //end of JMQ fix on Rb by 190709
  
  //JMQ 160909 fix: initial level density must be calculated according to the 
  // conditions at the initial compound nucleus 
  // (it has been removed from previous "if" for the residual) 

  if ( U < ExCN ) 
    {
      //JMQ fixed bug in units
      //VI moved the computation here
      G4double E0CN = ExCN - TCN*(G4Log(TCN/MeV) - 0.25*G4Log(aCN*MeV) 
                  - 1.25*G4Log(UxCN/MeV) 
                  + 2.0*std::sqrt(aCN*UxCN));

      InitialLevelDensity = (pi/12.0)*G4Exp((U-E0CN)/TCN)/TCN;
    } 
  else 
    {
      //VI speedup
      G4double x  = U-deltaCN;
      G4double x1 = std::sqrt(aCN*x);

      InitialLevelDensity = (pi/12.0)*G4Exp(2*x1)/(x*std::sqrt(x1));
    }

  //JMQ 190709 BUG : pi instead of sqrt(pi) must be here according 
  // to Furihata's report:
  Width *= pi*gg*GeometricalXS*Alpha/(12.0*InitialLevelDensity); 
   
  return Width;
}

Here is the call graph for this function:

Here is the caller graph for this function:

CCoeficient()

G4double G4GEMProbability::CCoeficient ( G4int  aZ ) const

inlineprivate

Definition at line 186 of file G4GEMProbability.hh.

{
  //JMQ 190709 C's values from Furihata's paper 
  //(notes added on proof in Dostrovskii's paper) 
  //data = {{20, 0.}, {30, -0.06}, {40, -0.10}, {50, -0.10}};
  G4double C = 0.0;
  if (aZ >= 50){
    C=-0.10/G4double(theA);
  } else if (aZ > 20) {
    C=(0.123482-0.00534691*aZ-0.0000610624*aZ*aZ+5.93719*1e-7*aZ*aZ*aZ+
       1.95687*1e-8*aZ*aZ*aZ*aZ)/G4double(theA);
  }
  return C;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Dump()

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:

EmissionProbability()

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:

GetA()

G4double G4GEMProbability::GetA ( void  ) const

inline

Definition at line 151 of file G4GEMProbability.hh.

{ 
  return theA;
}

GetA_asInt()

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:

GetCoulombBarrier()

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:

GetSpin()

G4double G4GEMProbability::GetSpin ( void  ) const

inline

Definition at line 156 of file G4GEMProbability.hh.

{ 
  return Spin; 
}

Here is the caller graph for this function:

GetZ()

G4double G4GEMProbability::GetZ ( void  ) const

inline

Definition at line 146 of file G4GEMProbability.hh.

{ 
  return theZ; 
}

GetZ_asInt()

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:

I0()

G4double G4GEMProbability::I0 ( G4double  t )

inlineprivate

Definition at line 230 of file G4GEMProbability.hh.

{
  return G4Exp(t) - 1.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

I1()

G4double G4GEMProbability::I1 ( G4double  t, G4double  tx  )

inlineprivate

Definition at line 235 of file G4GEMProbability.hh.

{
  return (t - tx + 1.0)*G4Exp(tx) - t - 1.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

I2()

G4double G4GEMProbability::I2 ( G4double  s0, G4double  sx  )

inlineprivate

Definition at line 241 of file G4GEMProbability.hh.

{
  G4double S = 1.0/std::sqrt(s0);
  G4double Sx = 1.0/std::sqrt(sx);
  
  G4double p1 = S*S*S*( 1.0 + S*S*( 1.5 + 3.75*S*S) );
  G4double p2 = Sx*Sx*Sx*( 1.0 + Sx*Sx*( 1.5 + 3.75*Sx*Sx) )*G4Exp(sx-s0);
  
  return p1-p2;
}

Here is the call graph for this function:

I3()

G4double G4GEMProbability::I3 ( G4double  s0, G4double  sx  )

private

Definition at line 239 of file G4GEMProbability.cc.

{
  G4double s2 = s0*s0;
  G4double sx2 = sx*sx;
  G4double S = 1.0/std::sqrt(s0);
  G4double S2 = S*S;
  G4double Sx = 1.0/std::sqrt(sx);
  G4double Sx2 = Sx*Sx;
  
  G4double p1 = S *(2.0 + S2 *( 4.0 + S2 *( 13.5 + S2 *( 60.0 + S2 * 325.125 ))));
  G4double p2 = Sx*Sx2 *(
             (s2-sx2) + Sx2 *(
                      (1.5*s2+0.5*sx2) + Sx2 *(
                                   (3.75*s2+0.25*sx2) + Sx2 *(
                                                  (12.875*s2+0.625*sx2) + Sx2 *(
                                                                (59.0625*s2+0.9375*sx2) + Sx2 *(324.8*s2+3.28*sx2))))));
  
  p2 *= G4Exp(sx-s0);
  
  return p1-p2; 
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator!=()

G4bool G4GEMProbability::operator!= ( const G4GEMProbability &  right ) const

private

operator=()

const G4GEMProbability& G4GEMProbability::operator= ( const G4GEMProbability &  right )

private

operator==()

G4bool G4GEMProbability::operator== ( const G4GEMProbability &  right ) const

private

SetCoulomBarrier()

void G4GEMProbability::SetCoulomBarrier ( const G4VCoulombBarrier *  aCoulombBarrierStrategy )

inline

Definition at line 167 of file G4GEMProbability.hh.

{
  theCoulombBarrierPtr = aCoulombBarrierStrategy;
}

Member Data Documentation

ExcitEnergies

std::vector<G4double> G4GEMProbability::ExcitEnergies

protected

Definition at line 126 of file G4GEMProbability.hh.

ExcitLifetimes

std::vector<G4double> G4GEMProbability::ExcitLifetimes

protected

Definition at line 132 of file G4GEMProbability.hh.

ExcitSpins

std::vector<G4double> G4GEMProbability::ExcitSpins

protected

Definition at line 129 of file G4GEMProbability.hh.

fG4pow

G4Pow* G4GEMProbability::fG4pow

private

Definition at line 107 of file G4GEMProbability.hh.

fPairCorr

G4PairingCorrection* G4GEMProbability::fPairCorr

private

Definition at line 108 of file G4GEMProbability.hh.

fPlanck

G4double G4GEMProbability::fPlanck

protected

Definition at line 123 of file G4GEMProbability.hh.

Spin

G4double G4GEMProbability::Spin

private

Definition at line 116 of file G4GEMProbability.hh.

theA

G4int G4GEMProbability::theA

private

Definition at line 112 of file G4GEMProbability.hh.

theCoulombBarrierPtr

const G4VCoulombBarrier* G4GEMProbability::theCoulombBarrierPtr

private

Definition at line 119 of file G4GEMProbability.hh.

theEvapLDPptr

G4VLevelDensityParameter* G4GEMProbability::theEvapLDPptr

private

Definition at line 110 of file G4GEMProbability.hh.

theZ

G4int G4GEMProbability::theZ

private

Definition at line 113 of file G4GEMProbability.hh.

---

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

• G4GEMProbability.hh
• G4GEMProbability.cc