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	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.

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

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 34 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::csNN

(

G4double

e

)

Definition at line 83 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; 
}

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::csPN

(

G4double

e

)

Definition at line 95 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; 
}

Here is the caller graph for this function:

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

Definition at line 109 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::G4cbrt

(

G4double

x

)

Definition at line 117 of file G4InuclSpecialFunctions.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::G4cbrt

(

G4int

n

)

Definition at line 121 of file G4InuclSpecialFunctions.cc.

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

Here is the call graph for this function:

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

Definition at line 151 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 178 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::getAL

(

G4int

A

)

Definition at line 79 of file G4InuclSpecialFunctions.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::inuclRndm

(

)

Definition at line 125 of file G4InuclSpecialFunctions.cc.

                                            { 
  return G4UniformRand(); 
} 

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::nucleiLevelDensity

(

G4int

A

)

Definition at line 33 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;
}

Here is the caller graph for this function:

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

(

)

Definition at line 144 of file G4InuclSpecialFunctions.cc.

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

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::randomGauss

(

G4double

sigma

)

Definition at line 129 of file G4InuclSpecialFunctions.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 55 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));
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4InuclSpecialFunctions::randomPHI

(

)

Definition at line 140 of file G4InuclSpecialFunctions.cc.

                                            { 
  return twopi * inuclRndm();
} 

Here is the call graph for this function:

Here is the caller graph for this function: