G4ElectronIonPair Class Reference

#include <G4ElectronIonPair.hh>

Public Member Functions
	G4ElectronIonPair ()
virtual	~G4ElectronIonPair ()
G4double	MeanNumberOfIonsAlongStep (const G4ParticleDefinition *, const G4Material *, G4double edepTotal, G4double edepNIEL=0.0)
G4double	MeanNumberOfIonsAlongStep (const G4Step *)
G4int	SampleNumberOfIonsAlongStep (const G4Step *)
std::vector< G4ThreeVector > *	SampleIonsAlongStep (const G4Step *)
G4int	ResidualeChargePostStep (const G4ParticleDefinition *, const G4TrackVector *secondary=0, G4int processSubType=-1)
G4int	ResidualeChargePostStep (const G4Step *)
G4double	FindG4MeanEnergyPerIonPair (const G4Material *)
void	DumpMeanEnergyPerIonPair ()
void	DumpG4MeanEnergyPerIonPair ()
void	SetVerbose (G4int)

Detailed Description

Definition at line 74 of file G4ElectronIonPair.hh.

Constructor & Destructor Documentation

G4ElectronIonPair::G4ElectronIonPair

(

)

Definition at line 58 of file G4ElectronIonPair.cc.

{
  verbose = 1;
  curMaterial = 0;
  curMeanEnergy = 0.0;
  nMaterials = 0;
  invFanoFactor = 1.0/0.2;
  Initialise();
}

G4ElectronIonPair::~G4ElectronIonPair ( )

virtual

Definition at line 70 of file G4ElectronIonPair.cc.

{}

Member Function Documentation

void G4ElectronIonPair::DumpG4MeanEnergyPerIonPair

(

)

Definition at line 187 of file G4ElectronIonPair.cc.

References python.hepunit::eV, G4cout, and G4endl.

{
  if(nMaterials > 0) {
    G4cout << "### G4ElectronIonPair: mean energy per ion pair "
       << " for Geant4 materials" << G4endl;
    for(G4int i=0; i<nMaterials; ++i) {
      G4cout << "   " << g4MatNames[i] << "    Epair= " 
         << g4MatData[i]/eV << " eV" << G4endl;
    }
  }
}

void G4ElectronIonPair::DumpMeanEnergyPerIonPair

(

)

Definition at line 168 of file G4ElectronIonPair.cc.

References python.hepunit::eV, G4cout, G4endl, G4Material::GetIonisation(), G4Material::GetMaterialTable(), G4IonisParamMat::GetMeanEnergyPerIonPair(), G4Material::GetName(), G4Material::GetNumberOfMaterials(), and test::x.

{
  G4int nmat = G4Material::GetNumberOfMaterials();
  const G4MaterialTable* mtable = G4Material::GetMaterialTable();
  if(nmat > 0) {
    G4cout << "### G4ElectronIonPair: mean energy per ion pair avalable:" << G4endl;
    for(G4int i=0; i<nmat; ++i) {
      const G4Material* mat = (*mtable)[i];
      G4double x = mat->GetIonisation()->GetMeanEnergyPerIonPair();
      if(x > 0.0) {
    G4cout << "   " << mat->GetName() << "   Epair=  " 
           << x/eV << " eV" << G4endl;
      }
    }
  }
}

G4double G4ElectronIonPair::FindG4MeanEnergyPerIonPair

(

const G4Material *

mat

)

Definition at line 145 of file G4ElectronIonPair.cc.

References python.hepunit::eV, G4cout, G4endl, G4Material::GetIonisation(), G4Material::GetName(), and G4IonisParamMat::SetMeanEnergyPerIonPair().

Referenced by MeanNumberOfIonsAlongStep().

{
  G4String name = mat->GetName();
  G4double res  = 0.0;

  // is this material in the vector?  
  for(G4int j=0; j<nMaterials; j++) {
    if(name == g4MatNames[j]) {
      res = g4MatData[j];
      mat->GetIonisation()->SetMeanEnergyPerIonPair(res);
      if(verbose > 0) {
    G4cout << "### G4ElectronIonPair::FindG4MeanEnergyPerIonPair for "
           << name << " Epair= " << res/eV << " eV is set"
           << G4endl;
      }
      break;
    }
  }
  return res;
}

G4double G4ElectronIonPair::MeanNumberOfIonsAlongStep	(	const G4ParticleDefinition *	part,
		const G4Material *	material,
		G4double	edepTotal,
		G4double	edepNIEL = 0.0
	)

Definition at line 75 of file G4ElectronIonPair.cc.

References FindG4MeanEnergyPerIonPair(), G4Material::GetIonisation(), G4IonisParamMat::GetMeanEnergyPerIonPair(), G4ParticleDefinition::GetPDGCharge(), and eplot::material.

Referenced by HistoManager::AddEnergy(), MeanNumberOfIonsAlongStep(), and SampleNumberOfIonsAlongStep().

{
  G4double nion = 0.0;

  // NIEL does not provide ionisation clusters
  if(edep > niel) {

    // neutral particles do not produce ionisation along step
    if(part->GetPDGCharge() != 0.0) {

      // select material
      if(material != curMaterial) {
    curMaterial = material;
    curMeanEnergy = material->GetIonisation()->GetMeanEnergyPerIonPair();

    // if mean energy is not defined then look into G4 DB
    if(0.0 == curMeanEnergy) {
      curMeanEnergy = FindG4MeanEnergyPerIonPair(material);
    } 
      }
      if(curMeanEnergy > 0.0) { nion = (edep - niel)/curMeanEnergy; }
    }
  }
  return nion;
}

G4double G4ElectronIonPair::MeanNumberOfIonsAlongStep ( const G4Step *  step )

inline

Definition at line 139 of file G4ElectronIonPair.hh.

References G4StepPoint::GetMaterial(), G4Step::GetNonIonizingEnergyDeposit(), G4Track::GetParticleDefinition(), G4Step::GetPreStepPoint(), G4Step::GetTotalEnergyDeposit(), G4Step::GetTrack(), and MeanNumberOfIonsAlongStep().

{
  return MeanNumberOfIonsAlongStep(step->GetTrack()->GetParticleDefinition(),
                   step->GetPreStepPoint()->GetMaterial(),
                   step->GetTotalEnergyDeposit(),
                   step->GetNonIonizingEnergyDeposit());
}

G4int G4ElectronIonPair::ResidualeChargePostStep	(	const G4ParticleDefinition *	,
		const G4TrackVector *	secondary = 0,
		G4int	processSubType = -1
	)

Definition at line 133 of file G4ElectronIonPair.cc.

Referenced by ResidualeChargePostStep().

{
  G4int nholes = 0;

  if(2 == subType || 12 == subType || 13 == subType) { nholes = 1; }
  return nholes;
}

G4int G4ElectronIonPair::ResidualeChargePostStep ( const G4Step *  step )

inline

Definition at line 157 of file G4ElectronIonPair.hh.

References G4Track::GetParticleDefinition(), G4Step::GetPostStepPoint(), G4StepPoint::GetProcessDefinedStep(), G4VProcess::GetProcessSubType(), G4Step::GetSecondary(), G4Step::GetTrack(), and ResidualeChargePostStep().

{
  G4int subtype = -1;
  const G4VProcess* proc = step->GetPostStepPoint()->GetProcessDefinedStep();
  if(proc) { subtype = proc->GetProcessSubType(); }
  return ResidualeChargePostStep(step->GetTrack()->GetParticleDefinition(),
                 step->GetSecondary(),
                 subtype);
}

std::vector< G4ThreeVector > * G4ElectronIonPair::SampleIonsAlongStep

(

const G4Step *

step

)

Definition at line 108 of file G4ElectronIonPair.cc.

References G4cout, G4endl, G4UniformRand, G4StepPoint::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), SampleNumberOfIonsAlongStep(), and test::v.

{
  std::vector<G4ThreeVector>* v = 0;

  G4int nion = SampleNumberOfIonsAlongStep(step);

  // sample ionisation along step
  if(nion > 0) {

    v = new std::vector<G4ThreeVector>;
    G4ThreeVector prePos = step->GetPreStepPoint()->GetPosition();
    G4ThreeVector deltaPos = step->GetPostStepPoint()->GetPosition() - prePos;  
    for(G4int i=0; i<nion; ++i) {
      v->push_back( prePos + deltaPos*G4UniformRand() );
    }
    if(verbose > 1 ) { 
      G4cout << "### G4ElectronIonPair::SampleIonisationPoints: "
         << v->size() << "  ion pairs are added" << G4endl;
    }
  }
  return v;
}

G4int G4ElectronIonPair::SampleNumberOfIonsAlongStep ( const G4Step *  step )

inline

Definition at line 148 of file G4ElectronIonPair.hh.

References G4lrint(), MeanNumberOfIonsAlongStep(), and G4INCL::DeJongSpin::shoot().

Referenced by HistoManager::AddEnergy(), and SampleIonsAlongStep().

{
  // use gamma distribution with mean value n=meanion and 
  // dispersion D=meanion/invFanoFactor
  G4double meanion = MeanNumberOfIonsAlongStep(step);
  return G4lrint(G4RandGamma::shoot(meanion*invFanoFactor,invFanoFactor));
} 

void G4ElectronIonPair::SetVerbose ( G4int  val )

inline

Definition at line 167 of file G4ElectronIonPair.hh.

Referenced by G4EmManager::SetVerbose(), and G4LossTableManager::SetVerbose().

{
  verbose = val;
}

---

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

• G4ElectronIonPair.hh
• G4ElectronIonPair.cc