G4hQAOModel Class Reference

#include <G4hQAOModel.hh>

Inheritance diagram for G4hQAOModel:

Public Member Functions
	G4hQAOModel ()
	~G4hQAOModel ()
G4bool	HasMaterial (const G4Material *)
G4double	StoppingPower (const G4Material *material, G4double kineticEnergy)
G4double	ElectronicStoppingPower (G4double z, G4double kineticEnergy) const
Public Member Functions inherited from G4VhElectronicStoppingPower
	G4VhElectronicStoppingPower ()
virtual	~G4VhElectronicStoppingPower ()
virtual G4double	StoppingPower (const G4Material *material, G4double kineticEnergy)=0
virtual G4bool	HasMaterial (const G4Material *material)=0
virtual G4double	ElectronicStoppingPower (G4double z, G4double kineticEnergy) const =0

Additional Inherited Members
Protected Member Functions inherited from G4VhElectronicStoppingPower
G4double	HeEffChargeSquare (const G4double z, const G4double kineticEnergyHe) const
G4double	GetHeMassAMU () const

Detailed Description

Definition at line 57 of file G4hQAOModel.hh.

Constructor & Destructor Documentation

G4hQAOModel()

G4hQAOModel::G4hQAOModel

(

)

Definition at line 59 of file G4hQAOModel.cc.

                        :G4VhElectronicStoppingPower(),
  numberOfMaterials(6),
  sizeL0(67),
  sizeL1(22),
  sizeL2(14),
  currentMaterial(0),
  currentElement(0)
{
  theZieglerFactor = eV*cm2*1.0e-15;
  thePlasmonFactor = 28.816 * 28.816  * 1e-6;
}

~G4hQAOModel()

G4hQAOModel::~G4hQAOModel

(

)

Definition at line 73 of file G4hQAOModel.cc.

{}

Member Function Documentation

ElectronicStoppingPower()

G4double G4hQAOModel::ElectronicStoppingPower ( G4double  z, G4double  kineticEnergy  ) const

virtual

Implements G4VhElectronicStoppingPower.

Definition at line 100 of file G4hQAOModel.cc.

{
  G4int Z = (G4int)z;
  G4int nbOfShell = GetNumberOfShell(Z);
  if(nbOfShell < 1) nbOfShell = 1;
 
  G4double dedx=0.0;
 
  G4double v = c_light * std::sqrt( 2.0 * kineticEnergy / proton_mass_c2 );
  G4double coeff = twopi*proton_mass_c2*z / (electron_mass_c2*theZieglerFactor) ;
  G4double fBetheVelocity = fine_structure_const * c_light / v;
  coeff *= fine_structure_const * fine_structure_const * hbarc_squared /
           kineticEnergy ;
 
  G4double l0Term = 0, l1Term = 0, l2Term = 0;
 
  for (G4int nos = 0 ; nos < nbOfShell ; nos++){
 
    G4double l0 = 0, l1 = 0, l2 = 0;
    G4double NormalizedEnergy = ( 2.0 * electron_mass_c2 * v * v  ) /
                          ( c_squared * GetShellEnergy(Z, nos) );
    G4double shStrength = GetShellStrength(Z, nos);
 
    l0 = GetL0(NormalizedEnergy);
    l0Term += shStrength  * l0;
 
    l1 = GetL1(NormalizedEnergy);
    l1Term += shStrength * l1;
 
    l2 = GetL2(NormalizedEnergy);
    l2Term += shStrength * l2;
 
/*
    //if(Z == 6){
    G4cout << nos << ". "
           << " E(MeV)= " << kineticEnergy/MeV
           << " normE= "  << NormalizedEnergy
           << " sh en= "  << GetShellEnergy(Z,nos)
           << " str= "  << shStrength
           << " v0/v= " << fBetheVelocity
           << " l0= " << l0Term
           << " l1= " << l1Term
           << " l2= " << l2Term
           << G4endl;
    //   }
*/
  }
 
  dedx = coeff * (l0Term - fBetheVelocity*l1Term + fBetheVelocity*fBetheVelocity*l2Term);
 
  //G4cout << " E(MeV)= " << kineticEnergy/MeV
  //       << " dedx(Mev/mm)= " << dedx*mm/MeV << G4endl;
 
  if(dedx < 0.0) dedx = 0.0;
  return dedx;
}

Referenced by StoppingPower().

HasMaterial()

G4bool G4hQAOModel::HasMaterial ( const G4Material *  )

inlinevirtual

Implements G4VhElectronicStoppingPower.

Definition at line 66 of file G4hQAOModel.hh.

{return true;};

StoppingPower()

G4double G4hQAOModel::StoppingPower ( const G4Material *  material, G4double  kineticEnergy  )

virtual

Implements G4VhElectronicStoppingPower.

Definition at line 78 of file G4hQAOModel.cc.

{
  G4double eloss = 0.0 ;
  G4int numberOfElements = material->GetNumberOfElements() ;
  const G4ElementVector* theElementVector   = material->GetElementVector();
  const G4double* theAtomicNumDensityVector = material->GetAtomicNumDensityVector() ;
  currentMaterial = material;
 
  //  loop for the elements in the material
  for (G4int i=0; i<numberOfElements; i++)
    {
      currentElement = (*theElementVector)[i] ;
      G4double z = currentElement->GetZ() ;
      eloss += ElectronicStoppingPower(z,kineticEnergy)*theAtomicNumDensityVector[i];
    }
 
  return eloss;
}

---

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

• G4hQAOModel.hh
• G4hQAOModel.cc