G4InuclSpecialFunctions Namespace Reference

Classes
class	paraMaker

Functions
G4double	bindingEnergy (G4int A, G4int Z)
G4double	bindingEnergyAsymptotic (G4int A, G4int Z)
G4double	FermiEnergy (G4int A, G4int Z, G4int ntype)
G4double	getAL (G4int A)
G4double	csNN (G4double e)
G4double	csPN (G4double e)
G4double	G4cbrt (G4double x)
G4double	G4cbrt (G4int n)
G4double	randomInuclPowers (G4double ekin, const G4double(&coeff)[4][4])
G4double	randomGauss (G4double sigma)
G4double	randomPHI ()
std::pair< G4double, G4double >	randomCOS_SIN ()
G4double	nucleiLevelDensity (G4int A)
G4LorentzVector	generateWithFixedTheta (G4double ct, G4double p, G4double mass=0.)
G4LorentzVector	generateWithRandomAngles (G4double p, G4double mass=0.)

Function Documentation

bindingEnergy()

G4double G4InuclSpecialFunctions::bindingEnergy	(	G4int	A,
		G4int	Z )

Definition at line 36 of file bindingEnergy.cc.

                                                                {
  // NOTE:  Test condition copied from G4NucleiProperties.cc; not encapsulated
  if (A < 1 || Z < 0 || Z > A) return 0.;
 
  return G4NucleiProperties::GetBindingEnergy(A, Z);
}

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), G4Fissioner::deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4CascadeDeexciteBase::explosion(), G4NucleiModel::fillBindingEnergies(), G4CascadeRecoilMaker::goodNucleus(), and G4NucleiModel::initializeCascad().

bindingEnergyAsymptotic()

G4double G4InuclSpecialFunctions::bindingEnergyAsymptotic	(	G4int	A,
		G4int	Z )

Definition at line 33 of file bindingEnergyAsymptotic.cc.

                                                                          {
  // calculates the nuclei binding energy 
  // using smooth liquid high energy formula  
  G4double X = (1.0 - 2.0*Z/A); X *= X;
  G4double X1 = G4cbrt(A);
  G4double X2 = X1 * X1;
  G4double X3 = 1.0 / X1;
  G4double X4 = 1.0 / X2;
  G4double X5 = (1.0 - 0.62025 * X4); X5 *= X5;
  G4double Z13 = G4cbrt(Z);
 
  G4double DM = 17.035 * (1.0 - 1.846 * X) * A -
    25.8357 * (1.0 - 1.712 * X) * X2 * X5 -
    0.779 * Z * (Z - 1) * X3 * 
    (1.0 - 1.5849 * X4 + 1.2273 / A + 1.5772 * X4 * X4) +
    0.4328 * Z13*Z13*Z13*Z13 * X3 * 
    (1.0 - 0.57811 * X3 - 0.14518 * X4 + 0.496 / A);
 
  return DM;
}

Referenced by G4Fissioner::deExcite().

csNN()

G4double G4InuclSpecialFunctions::csNN

(

G4double

e

)

Definition at line 82 of file G4InuclSpecialFunctions.cc.

                                                 {
  G4double snn;
 
  if (e < 40.0) {
    snn = -1174.8 / (e * e) + 3088.5 / e + 5.3107;
  } else {
    snn = 93074.0 / (e * e) - 11.148 / e + 22.429;
  }
 
  return snn; 
}

Referenced by G4NonEquilibriumEvaporator::deExcite().

csPN()

G4double G4InuclSpecialFunctions::csPN

(

G4double

e

)

Definition at line 94 of file G4InuclSpecialFunctions.cc.

                                                 {
  G4double spn;
 
  if (e < 40.0) {
    spn = -5057.4 / (e * e) + 9069.2 / e + 6.9466;
  } else {
    spn = 239380.0 / (e * e) + 1802.0 / e + 27.147;
  }
 
  return spn; 
}

Referenced by G4NonEquilibriumEvaporator::deExcite().

FermiEnergy()

G4double G4InuclSpecialFunctions::FermiEnergy	(	G4int	A,
		G4int	Z,
		G4int	ntype )

Definition at line 108 of file G4InuclSpecialFunctions.cc.

                                                                           {
  G4Pow* g4pow = G4Pow::GetInstance();
  const G4double C = 55.4 / g4pow->Z23(A);
  G4double arg = (ntype==0) ? g4pow->Z23(A-Z) : g4pow->Z23(Z);
 
  return C * arg;
}

Referenced by G4NonEquilibriumEvaporator::deExcite().

G4cbrt() [1/2]

G4double G4InuclSpecialFunctions::G4cbrt

(

G4double

x

)

Definition at line 116 of file G4InuclSpecialFunctions.cc.

                                                   {
  return x==0 ? 0. : (x<0?-1.:1.)*G4Exp(G4Log(std::fabs(x))/3.);
}

Referenced by bindingEnergyAsymptotic(), G4EquilibriumEvaporator::deExcite(), G4Fissioner::deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4NucleiModel::fillPotentials(), G4NucleiModel::generateModel(), G4NucleiModel::generateNucleonMomentum(), getAL(), G4CascadeColliderBase::inelasticInteractionPossible(), and G4IntraNucleiCascader::initialize().

G4cbrt() [2/2]

G4double G4InuclSpecialFunctions::G4cbrt

(

G4int

n

)

Definition at line 120 of file G4InuclSpecialFunctions.cc.

                                                {
  return n==0 ? 0. : (n<0?-1.:1.)*G4Pow::GetInstance()->Z13(std::abs(n));
}

generateWithFixedTheta()

G4LorentzVector G4InuclSpecialFunctions::generateWithFixedTheta	(	G4double	ct,
		G4double	p,
		G4double	mass = 0. )

Definition at line 147 of file G4InuclSpecialFunctions.cc.

                                       {
  G4double phi = randomPHI();
  G4double pt = p * std::sqrt(std::fabs(1.0 - ct * ct));
 
  // Buffers to avoid memory thrashing
  static G4ThreadLocal G4ThreeVector *pvec_G4MT_TLS_ = 0;
  if (!pvec_G4MT_TLS_) {
    pvec_G4MT_TLS_ = new G4ThreeVector;
    G4AutoDelete::Register(pvec_G4MT_TLS_);
  }
  G4ThreeVector &pvec = *pvec_G4MT_TLS_;
 
  static G4ThreadLocal G4LorentzVector *momr_G4MT_TLS_ = 0;
  if (!momr_G4MT_TLS_) {
    momr_G4MT_TLS_ = new G4LorentzVector;
    G4AutoDelete::Register(momr_G4MT_TLS_);
  }
  G4LorentzVector &momr = *momr_G4MT_TLS_;
 
  pvec.set(pt*std::cos(phi), pt*std::sin(phi), p*ct);
  momr.setVectM(pvec, mass);
 
  return momr;
}

Referenced by G4CascadeFinalStateAlgorithm::FillDirManyBody(), G4CascadeFinalStateAlgorithm::FillDirThreeBody(), and G4NucleiModel::initializeCascad().

generateWithRandomAngles()

G4LorentzVector G4InuclSpecialFunctions::generateWithRandomAngles	(	G4double	p,
		G4double	mass = 0. )

Definition at line 174 of file G4InuclSpecialFunctions.cc.

                                                                           {
  std::pair<G4double, G4double> COS_SIN = randomCOS_SIN();
  G4double phi = randomPHI();
  G4double pt = p * COS_SIN.second;
 
  // Buffers to avoid memory thrashing  
  static G4ThreadLocal G4ThreeVector *pvec_G4MT_TLS_ = 0;
  if (!pvec_G4MT_TLS_) {
    pvec_G4MT_TLS_ = new G4ThreeVector;
    G4AutoDelete::Register(pvec_G4MT_TLS_);
  }
  G4ThreeVector &pvec = *pvec_G4MT_TLS_;
 
  static G4ThreadLocal G4LorentzVector *momr_G4MT_TLS_ = 0;
  if (!momr_G4MT_TLS_) {
    momr_G4MT_TLS_ = new G4LorentzVector;
    G4AutoDelete::Register(momr_G4MT_TLS_);
  }
  G4LorentzVector &momr = *momr_G4MT_TLS_;
 
  pvec.set(pt*std::cos(phi), pt*std::sin(phi), p*COS_SIN.first);
  momr.setVectM(pvec, mass);
 
  return momr;
}

Referenced by G4EquilibriumEvaporator::deExcite(), G4Fissioner::deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4NucleiModel::generateNucleonMomentum(), and G4NucleiModel::initializeCascad().

getAL()

G4double G4InuclSpecialFunctions::getAL

(

G4int

A

)

Definition at line 78 of file G4InuclSpecialFunctions.cc.

                                               {
  return 0.76 + 2.2 / G4cbrt(A);
}

Referenced by G4EquilibriumEvaporator::deExcite(), and G4NonEquilibriumEvaporator::deExcite().

nucleiLevelDensity()

G4double G4InuclSpecialFunctions::nucleiLevelDensity

(

G4int

A

)

Definition at line 32 of file nucleiLevelDensity.cc.

                                                            {
  const G4double NLD[226] = {
    // 20 - 29
    3.94, 3.84, 3.74, 3.64, 3.55, 4.35, 4.26, 4.09, 3.96, 4.18,
    // 30 - 39
    4.39, 4.61, 4.82, 4.44, 4.44, 4.43, 4.42, 5.04, 5.66, 5.8,
    // 40 - 49
    5.95, 5.49, 6.18, 7.11, 6.96, 7.2, 7.73, 6.41, 6.85, 6.77,
    // 50 - 59
    6.91, 7.3, 7.2, 6.86, 8.06, 7.8, 7.82, 8.41, 8.13, 7.19,
    // 60 - 69
    8.35, 8.13, 8.02, 8.93, 8.9, 9.7, 9.65, 10.55, 9.38, 9.72,
    // 70 - 79
    10.66, 11.98, 12.76, 12.1, 12.86, 13.0, 12.81, 12.8, 12.65, 12.0,
    // 80 - 89
    12.69, 14.05, 13.33, 13.28, 13.22, 13.17, 8.66, 11.03, 10.4, 13.47,
    // 90 - 99
    10.17, 12.22, 11.62, 12.95, 13.15, 13.57, 12.87, 16.2, 14.71, 15.69,
    // 100 - 109
    14.1, 18.56, 16.22, 16.7, 17.13, 17.0, 16.86, 16.2, 15.61, 16.8,
    // 110 - 119
    17.93, 17.5, 16.97, 17.3, 17.6, 15.78, 16.8, 17.49, 16.03, 15.08,
    // 120 - 129
    16.74, 17.74, 17.43, 18.1, 17.1, 19.01, 17.02, 17.0, 17.02, 18.51,
    // 130 - 139
    17.2, 16.8, 16.97, 16.14, 16.91, 17.69, 15.5, 14.56, 14.35, 16.5,
    // 140 - 149
    18.29, 17.8, 17.05, 21.31, 19.15, 19.5, 19.78, 20.3, 20.9, 21.9,
    // 150 - 159
    22.89, 25.68, 24.6, 24.91, 23.24, 22.9, 22.46, 21.98, 21.64, 21.8,
    // 160 - 169
    21.85, 21.7, 21.69, 23.7, 21.35, 23.03, 20.66, 21.81, 20.77, 22.2,
    // 170 - 179
    22.58, 22.55, 21.45, 21.16, 21.02, 20.87, 22.09, 22.0, 21.28, 23.05,
    // 180 - 189
    21.7, 21.18, 22.28, 23.0, 22.11, 23.56, 22.83, 24.88, 22.6, 23.5,
    // 190 - 199
    23.89, 23.9, 23.94, 21.16, 22.3, 21.7, 21.19, 20.7, 20.29, 21.32,
    // 200 - 209
    19.0, 17.93, 17.85, 15.7, 13.54, 11.9, 10.02, 10.48, 10.28, 11.72,
    // 210 - 219
    13.81, 14.7, 15.5, 16.3, 17.2, 18.0, 18.9, 19.7, 20.6, 21.4,
    // 220 - 229
    22.3, 23.1, 24.0, 24.8, 25.6, 26.5, 27.3, 28.2, 29.0, 29.9,
    // 230 - 239
    30.71, 30.53, 31.45, 29.6, 30.2, 30.65, 30.27, 29.52, 30.08, 29.8,
    // 240 - 245
    29.87, 30.25, 30.5, 29.8, 29.17, 28.67};
 
  return (a>=20) ? NLD[a-20] : 0.1*a;
}

Referenced by G4Fissioner::deExcite().

randomCOS_SIN()

std::pair< G4double, G4double > G4InuclSpecialFunctions::randomCOS_SIN

(

)

Definition at line 140 of file G4InuclSpecialFunctions.cc.

                                                                   {
  G4double CT = 1.0 - 2.0*G4UniformRand();
 
  return std::pair<G4double, G4double>(CT, std::sqrt(1.0 - CT*CT));
}

Referenced by generateWithRandomAngles().

randomGauss()

G4double G4InuclSpecialFunctions::randomGauss

(

G4double

sigma

)

Definition at line 125 of file G4InuclSpecialFunctions.cc.

                                                            {
  const G4double eps = 1.0e-6;
  G4double r1 = G4UniformRand();
  r1 = r1 > eps ? r1 : eps;
  G4double r2 = G4UniformRand();
  r2 = r2 > eps ? r2 : eps;
  r2 = r2 < 1.0 - eps ? r2 : 1.0 - eps; 
 
  return sigma * std::sin(twopi * r1) * std::sqrt(-2.0 * G4Log(r2)); 
} 

Referenced by G4Fissioner::deExcite().

randomInuclPowers()

G4double G4InuclSpecialFunctions::randomInuclPowers	(	G4double	ekin,
		const G4double(&)	coeff[4][4] )

Definition at line 54 of file G4InuclSpecialFunctions.cc.

                                                      {
  G4Pow* theG4Pow = G4Pow::GetInstance();
 
  G4double S = G4UniformRand();       // Random fraction for expansion
 
  G4double C, V;
  G4double PQ=0., PR=0.;
  for (G4int i=0; i<4; i++) {
    V = 0.0;
    for (G4int k=0; k<4; k++) {
      C = coeff[i][k];
      V += C * theG4Pow->powN(ekin, k);
    }
 
    PQ += V;
    PR += V * theG4Pow->powN(S, i);
  }
 
  return std::sqrt(S) * (PR + (1-PQ)*(S*S*S*S));
}

Referenced by G4InuclParamAngDst::GetCosTheta(), and G4InuclParamMomDst::GetMomentum().

randomPHI()

G4double G4InuclSpecialFunctions::randomPHI

(

)

Definition at line 136 of file G4InuclSpecialFunctions.cc.

                                            { 
  return twopi*G4UniformRand();
} 

Referenced by generateWithFixedTheta(), generateWithRandomAngles(), and G4NucleiModel::initializeCascad().