G4PEEffectModel Class Reference

#include <G4PEEffectModel.hh>

Inheritance diagram for G4PEEffectModel:

Public Member Functions
	G4PEEffectModel (const G4ParticleDefinition *p=0, const G4String &nam="PhotoElectric")
virtual	~G4PEEffectModel ()
virtual void	Initialise (const G4ParticleDefinition *, const G4DataVector &)
virtual G4double	ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double, G4double)
virtual G4double	CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)

---

Detailed Description

Definition at line 61 of file G4PEEffectModel.hh.

---

Constructor & Destructor Documentation

G4PEEffectModel::G4PEEffectModel	(	const G4ParticleDefinition *	p = 0,
		const G4String &	nam = "PhotoElectric"
	)

Definition at line 66 of file G4PEEffectModel.cc.

References G4Electron::Electron(), G4cout, G4endl, G4Gamma::Gamma(), and G4VEmModel::SetAngularDistribution().

  : G4VEmModel(nam)
{
  G4cout << "### G4PEEffectModel is obsolete "
         << "and will be removed for the next release." << G4endl;

  theGamma    = G4Gamma::Gamma();
  theElectron = G4Electron::Electron();
  fminimalEnergy = 1.0*eV;
  fParticleChange = 0;

  // default generator
  SetAngularDistribution(new G4SauterGavrilaAngularDistribution());
}

G4PEEffectModel::~G4PEEffectModel

(

)

[virtual]

Definition at line 84 of file G4PEEffectModel.cc.

{}

---

Member Function Documentation

G4double G4PEEffectModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	,
		G4double	kinEnergy,
		G4double	Z,
		G4double	A,
		G4double	,
		G4double
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 99 of file G4PEEffectModel.cc.

References G4SandiaTable::GetSandiaCofPerAtom().

Referenced by G4AdjointPhotoElectricModel::AdjointCrossSectionPerAtom().

{
  G4double* SandiaCof = G4SandiaTable::GetSandiaCofPerAtom((G4int)Z, energy);

  G4double energy2 = energy*energy;
  G4double energy3 = energy*energy2;
  G4double energy4 = energy2*energy2;

  return SandiaCof[0]/energy  + SandiaCof[1]/energy2 +
    SandiaCof[2]/energy3 + SandiaCof[3]/energy4;
}

G4double G4PEEffectModel::CrossSectionPerVolume	(	const G4Material *	,
		const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	cutEnergy,
		G4double	maxEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 116 of file G4PEEffectModel.cc.

References G4SandiaTable::GetSandiaCofForMaterial(), and G4Material::GetSandiaTable().

{
  G4double* SandiaCof = 
    material->GetSandiaTable()->GetSandiaCofForMaterial(energy);
                                
  G4double energy2 = energy*energy;
  G4double energy3 = energy*energy2;
  G4double energy4 = energy2*energy2;
          
  return SandiaCof[0]/energy  + SandiaCof[1]/energy2 +
    SandiaCof[2]/energy3 + SandiaCof[3]/energy4; 
}

void G4PEEffectModel::Initialise	(	const G4ParticleDefinition *	,
		const G4DataVector &
	)			[virtual]

Implements G4VEmModel.

Definition at line 89 of file G4PEEffectModel.cc.

References G4VEmModel::GetParticleChangeForGamma(), and G4VEmModel::SetDeexcitationFlag().

{
  // always false before the run
  SetDeexcitationFlag(false);
  if(!fParticleChange) { fParticleChange = GetParticleChangeForGamma(); }
}

void G4PEEffectModel::SampleSecondaries	(	std::vector< G4DynamicParticle * > *	,
		const G4MaterialCutsCouple *	,
		const G4DynamicParticle *	,
		G4double	tmin,
		G4double	maxEnergy
	)			[virtual]

Implements G4VEmModel.

Definition at line 134 of file G4PEEffectModel.cc.

References G4InuclSpecialFunctions::bindingEnergy(), fStopAndKill, G4VEmModel::GetAngularDistribution(), G4Element::GetAtomicShell(), G4DynamicParticle::GetKineticEnergy(), G4MaterialCutsCouple::GetMaterial(), G4DynamicParticle::GetMomentumDirection(), G4Element::GetNbOfAtomicShells(), G4VParticleChange::ProposeLocalEnergyDeposit(), G4VParticleChange::ProposeTrackStatus(), G4VEmAngularDistribution::SampleDirection(), G4VEmModel::SelectRandomAtom(), and G4ParticleChangeForGamma::SetProposedKineticEnergy().

{
  const G4Material* aMaterial = couple->GetMaterial();

  G4double energy = aDynamicPhoton->GetKineticEnergy();
  G4ParticleMomentum PhotonDirection = aDynamicPhoton->GetMomentumDirection();

  // select randomly one element constituing the material.
  const G4Element* anElement = SelectRandomAtom(aMaterial,theGamma,energy);
  
  //
  // Photo electron
  //

  // Select atomic shell
  G4int nShells = anElement->GetNbOfAtomicShells();
  G4int i  = 0;  
  while ((i<nShells) && (energy<anElement->GetAtomicShell(i))) { ++i; }

  G4double edep = energy;

  // no shell available
  if (i < nShells) {
  
    G4double bindingEnergy = anElement->GetAtomicShell(i);
    G4double elecKineEnergy = energy - bindingEnergy;

    if (elecKineEnergy > fminimalEnergy) {

      edep = bindingEnergy;
      G4ThreeVector elecDirection =
        GetAngularDistribution()->SampleDirection(aDynamicPhoton, 
                                                  elecKineEnergy,
                                                  i, 
                                                  couple->GetMaterial());

      G4DynamicParticle* aParticle = 
        new G4DynamicParticle(theElectron, elecDirection, elecKineEnergy);
      fvect->push_back(aParticle);

    } 
  }
    
  fParticleChange->SetProposedKineticEnergy(0.);
  fParticleChange->ProposeTrackStatus(fStopAndKill);
  if(edep > 0.0) {
    fParticleChange->ProposeLocalEnergyDeposit(edep);
  }
}

---

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

• G4PEEffectModel.hh
• G4PEEffectModel.cc

---