G4InuclSpecialFunctions Namespace Reference

Data Structures
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	inuclRndm ()
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

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

Definition at line 37 of file bindingEnergy.cc.

References G4NucleiProperties::GetBindingEnergy().

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), G4eIonisationSpectrum::AverageEnergy(), G4RDeIonisationSpectrum::AverageEnergy(), G4QMDNucleus::CalEnergyAndAngularMomentumInCM(), G4NonEquilibriumEvaporator::deExcite(), G4BigBanger::deExcite(), G4Fissioner::deExcite(), G4EquilibriumEvaporator::deExcite(), G4CascadeDeexciteBase::explosion(), G4NucleiModel::fillBindingEnergies(), G4AtomicTransitionManager::G4AtomicTransitionManager(), G4RDAtomicTransitionManager::G4RDAtomicTransitionManager(), G4ProtonField::GetBarrier(), G4SigmaPlusField::GetField(), G4KaonZeroField::GetField(), G4SigmaMinusField::GetField(), G4SigmaZeroField::GetField(), G4KaonMinusField::GetField(), G4KaonPlusField::GetField(), G4PionZeroField::GetField(), G4PionPlusField::GetField(), G4AntiProtonField::GetField(), G4PionMinusField::GetField(), G4CascadeRecoilMaker::goodNucleus(), G4NucleiModel::initializeCascad(), G4LowEnergyPhotoElectric::PostStepDoIt(), G4LowEnergyIonisation::PostStepDoIt(), G4hImpactIonisation::PostStepDoIt(), G4eIonisationSpectrum::Probability(), G4RDeIonisationSpectrum::Probability(), G4DeltaAngle::SampleDirection(), G4eIonisationSpectrum::SampleEnergy(), G4RDeIonisationSpectrum::SampleEnergy(), G4LivermorePolarizedPhotoElectricModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), G4DNARuddIonisationModel::SampleSecondaries(), G4DNABornIonisationModel::SampleSecondaries(), G4LivermorePhotoElectricModel::SampleSecondaries(), G4LivermoreIonisationModel::SampleSecondaries(), G4PenelopePhotoElectricModel::SampleSecondaries(), G4KleinNishinaModel::SampleSecondaries(), G4MuElecInelasticModel::SampleSecondaries(), G4MicroElecInelasticModel::SampleSecondaries(), G4PEEffectFluoModel::SampleSecondaries(), G4PenelopeComptonModel::SampleSecondaries(), G4PenelopeIonisationModel::SampleSecondaries(), and G4KM_OpticalEqRhs::SetFactor().

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

  return G4NucleiProperties::GetBindingEnergy(A, Z);
}

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

Definition at line 33 of file bindingEnergyAsymptotic.cc.

References G4cbrt().

Referenced by G4Fissioner::deExcite().

                                                                          {
  // 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 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 * G4cbrt(Z*Z*Z*Z) * X3 * 
    (1.0 - 0.57811 * X3 - 0.14518 * X4 + 0.496 / A);

  return DM;
}

G4double G4InuclSpecialFunctions::csNN

(

G4double

e

)

Definition at line 79 of file G4InuclSpecialFunctions.cc.

Referenced by G4NonEquilibriumEvaporator::deExcite().

                                                 {
  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; 
}

G4double G4InuclSpecialFunctions::csPN

(

G4double

e

)

Definition at line 91 of file G4InuclSpecialFunctions.cc.

Referenced by G4NonEquilibriumEvaporator::deExcite().

                                                 {
  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; 
}

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

Definition at line 105 of file G4InuclSpecialFunctions.cc.

References G4cbrt().

Referenced by G4NonEquilibriumEvaporator::deExcite().

                                                                           {
  const G4double C = 55.4;
  G4double arg = (ntype==0) ? G4double(A-Z)/A : G4double(Z)/A;

  return C * G4cbrt(arg*arg);   // 2/3 power
}

G4double G4InuclSpecialFunctions::G4cbrt

(

G4double

x

)

Definition at line 112 of file G4InuclSpecialFunctions.cc.

References G4Exp(), and G4Log().

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

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

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

Definition at line 142 of file G4InuclSpecialFunctions.cc.

References G4ThreadLocal, randomPHI(), CLHEP::Hep3Vector::set(), and CLHEP::HepLorentzVector::setVectM().

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

                                       {
  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;
  G4ThreeVector &pvec = *pvec_G4MT_TLS_;

  static G4ThreadLocal G4LorentzVector *momr_G4MT_TLS_ = 0;
  if (!momr_G4MT_TLS_) momr_G4MT_TLS_ = new G4LorentzVector;
  G4LorentzVector &momr = *momr_G4MT_TLS_;

  pvec.set(pt*std::cos(phi), pt*std::sin(phi), p*ct);
  momr.setVectM(pvec, mass);

  return momr;
}

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

Definition at line 163 of file G4InuclSpecialFunctions.cc.

References G4ThreadLocal, randomCOS_SIN(), randomPHI(), CLHEP::Hep3Vector::set(), and CLHEP::HepLorentzVector::setVectM().

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

                                                                           {
  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;
  G4ThreeVector &pvec = *pvec_G4MT_TLS_;

  static G4ThreadLocal G4LorentzVector *momr_G4MT_TLS_ = 0;
  if (!momr_G4MT_TLS_) momr_G4MT_TLS_ = new G4LorentzVector;
  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;
}

G4double G4InuclSpecialFunctions::getAL

(

G4int

A

)

Definition at line 75 of file G4InuclSpecialFunctions.cc.

References G4cbrt().

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

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

G4double G4InuclSpecialFunctions::inuclRndm

(

)

Definition at line 116 of file G4InuclSpecialFunctions.cc.

References G4UniformRand.

Referenced by G4NonEquilibriumEvaporator::deExcite(), G4Fissioner::deExcite(), G4EquilibriumEvaporator::deExcite(), G4CascadeFinalStateAlgorithm::GenerateCosTheta(), G4NucleiModel::generateInteractionLength(), G4NucleiModel::generateInteractionPartners(), G4NucleiModel::generateNucleonMomentum(), G4InuclParamAngDst::GetCosTheta(), G4NucleiModel::initializeCascad(), randomCOS_SIN(), randomGauss(), randomPHI(), and G4IntraNucleiCascader::setupCascade().

                                            { 
  return G4UniformRand(); 
} 

G4double G4InuclSpecialFunctions::nucleiLevelDensity

(

G4int

A

)

Definition at line 33 of file nucleiLevelDensity.cc.

References test::a.

Referenced by G4Fissioner::deExcite().

                                                            {
  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;
}

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

(

)

Definition at line 135 of file G4InuclSpecialFunctions.cc.

References inuclRndm().

Referenced by generateWithRandomAngles().

                                                                   {
  G4double CT = 1.0 - 2.0 * inuclRndm();

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

G4double G4InuclSpecialFunctions::randomGauss

(

G4double

sigma

)

Definition at line 120 of file G4InuclSpecialFunctions.cc.

References G4Log(), inuclRndm(), and python.hepunit::twopi.

Referenced by G4Fissioner::deExcite().

                                                            {
  const G4double eps = 1.0e-6;
  G4double r1 = inuclRndm();
  r1 = r1 > eps ? r1 : eps;
  G4double r2 = inuclRndm();
  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)); 
} 

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

Definition at line 51 of file G4InuclSpecialFunctions.cc.

References G4UniformRand, G4Pow::GetInstance(), and G4Pow::powN().

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

                                                      {
  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));
}

G4double G4InuclSpecialFunctions::randomPHI

(

)

Definition at line 131 of file G4InuclSpecialFunctions.cc.

References inuclRndm(), and python.hepunit::twopi.

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

                                            { 
  return twopi * inuclRndm();
}