80 for (
size_t k=0 ; k < nMatElements ; k++ )
82 if(atomDensities[k] > maxdens) {
83 maxdens=atomDensities[k];
96 G4double el=30.724*z1*z2*std::sqrt(zpow)*asum/a2;
97 G4double fl=0.0793*z23[z1]*std::sqrt(z2*asum*asum*asum/(a1*a1*a1*a2))
98 /std::pow(zpow, 0.75);
102 1.0/(1+fl*(3.4008*std::pow(eps, 0.16667)+0.40244*std::pow(eps, 0.75)+
eps));
const G4double * GetVecNbOfAtomsPerVolume() const
static constexpr double Avogadro
static constexpr double eV
G4LindhardRobinsonPartition()
virtual G4double PartitionNIEL(G4int z1, G4double a1, const G4Material *material, G4double energy) const
static constexpr double amu
static const G4double eps
const G4Element * GetElement(G4int iel) const
size_t GetNumberOfElements() const