G4LindhardRobinsonPartition Class Reference

#include <G4LindhardPartition.hh>

Inheritance diagram for G4LindhardRobinsonPartition:

Collaboration diagram for G4LindhardRobinsonPartition:

Public Member Functions
	G4LindhardRobinsonPartition ()
virtual	~G4LindhardRobinsonPartition ()
virtual G4double	PartitionNIEL (G4int z1, G4double a1, const G4Material *material, G4double energy) const
Public Member Functions inherited from G4VNIELPartition
	G4VNIELPartition ()
virtual	~G4VNIELPartition ()

Public Attributes
G4double	z23 [120]
size_t	max_z

Detailed Description

Definition at line 56 of file G4LindhardPartition.hh.

Constructor & Destructor Documentation

G4LindhardRobinsonPartition::G4LindhardRobinsonPartition

(

)

Definition at line 66 of file G4LindhardPartition.cc.

{
  max_z = 120;
  for(size_t i=1; i<max_z; i++) {z23[i]=std::pow((G4double)i, 2./3.);}
}

virtual G4LindhardRobinsonPartition::~G4LindhardRobinsonPartition ( )

inlinevirtual

Definition at line 60 of file G4LindhardPartition.hh.

{ }

Member Function Documentation

G4double G4LindhardRobinsonPartition::PartitionNIEL ( G4int  z1, G4double  a1, const G4Material *  material, G4double  energy  ) const

virtual

Implements G4VNIELPartition.

Definition at line 72 of file G4LindhardPartition.cc.

{
  size_t nMatElements = material->GetNumberOfElements();
        
  const G4double *atomDensities=material->GetVecNbOfAtomsPerVolume();
  G4double maxdens=0.0;
  size_t maxindex=0;
  for (size_t k=0 ; k < nMatElements ; k++ )
    {
      if(atomDensities[k] > maxdens) {
        maxdens=atomDensities[k];
        maxindex=k;
      }
    }
  const G4Element *element=material->GetElement(maxindex);

  G4int z2=G4int(element->GetZ());
        
  G4double a2=element->GetA()/(Avogadro*amu);
        
  G4double zpow=z23[z1]+z23[z2];
  G4double asum=a1+a2;
        
  G4double el=30.724*z1*z2*std::sqrt(zpow)*asum/a2;
  G4double fl=0.0793*z23[z1]*std::sqrt(z2*asum*asum*asum/(a1*a1*a1*a2))
              /std::pow(zpow, 0.75);
  G4double eps=(energy/eV)*(1.0/el);

  return 
    1.0/(1+fl*(3.4008*std::pow(eps, 0.16667)+0.40244*std::pow(eps, 0.75)+eps));
}

Here is the call graph for this function:

Member Data Documentation

size_t G4LindhardRobinsonPartition::max_z

Definition at line 67 of file G4LindhardPartition.hh.

G4double G4LindhardRobinsonPartition::z23[120]

Definition at line 66 of file G4LindhardPartition.hh.

---

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

• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm7/include/G4LindhardPartition.hh
• source/geant4.10.03.p03/examples/extended/electromagnetic/TestEm7/src/G4LindhardPartition.cc