G4EmSaturation Class Reference

#include <G4EmSaturation.hh>

Public Member Functions
	G4EmSaturation ()
virtual	~G4EmSaturation ()
G4double	VisibleEnergyDeposition (const G4ParticleDefinition *, const G4MaterialCutsCouple *, G4double length, G4double edepTotal, G4double edepNIEL=0.0)
G4double	VisibleEnergyDeposition (const G4Step *)
G4double	FindG4BirksCoefficient (const G4Material *)
void	SetVerbose (G4int)
void	DumpBirksCoefficients ()
void	DumpG4BirksCoefficients ()

---

Detailed Description

Definition at line 69 of file G4EmSaturation.hh.

---

Constructor & Destructor Documentation

G4EmSaturation::G4EmSaturation

(

)

Definition at line 58 of file G4EmSaturation.cc.

References G4NistManager::Instance().

{
  verbose = 1;
  manager = 0;

  curMaterial = 0;
  curBirks    = 0.0;
  curRatio    = 1.0;
  curChargeSq = 1.0;
  nMaterials  = 0;

  electron = 0;
  proton   = 0;
  nist     = G4NistManager::Instance();

  Initialise(); 
}

G4EmSaturation::~G4EmSaturation

(

)

[virtual]

Definition at line 78 of file G4EmSaturation.cc.

{}

---

Member Function Documentation

void G4EmSaturation::DumpBirksCoefficients

(

)

Definition at line 233 of file G4EmSaturation.cc.

References G4cout, and G4endl.

{
  if(nMaterials > 0) {
    G4cout << "### Birks coeffitients used in run time" << G4endl;
    for(G4int i=0; i<nMaterials; ++i) {
      G4double br = matPointers[i]->GetIonisation()->GetBirksConstant();
      G4cout << "   " << matNames[i] << "     " 
             << br*MeV/mm << " mm/MeV" << "     "
             << br*matPointers[i]->GetDensity()*MeV*cm2/g 
             << " g/cm^2/MeV" 
             << G4endl;
    }
  }
}

void G4EmSaturation::DumpG4BirksCoefficients

(

)

Definition at line 250 of file G4EmSaturation.cc.

References G4cout, and G4endl.

{
  if(nG4Birks > 0) {
    G4cout << "### Birks coeffitients for Geant4 materials" << G4endl;
    for(G4int i=0; i<nG4Birks; ++i) {
      G4cout << "   " << g4MatNames[i] << "   " 
             << g4MatData[i]*MeV/mm << " mm/MeV" << G4endl;
    }
  }
}

G4double G4EmSaturation::FindG4BirksCoefficient

(

const G4Material *

)

Definition at line 136 of file G4EmSaturation.cc.

References G4cout, G4endl, and G4Material::GetName().

{
  G4String name = mat->GetName();
  // is this material in the vector?
  
  for(G4int j=0; j<nG4Birks; ++j) {
    if(name == g4MatNames[j]) {
      if(verbose > 0) 
        G4cout << "### G4EmSaturation::FindG4BirksCoefficient for "
               << name << " is " << g4MatData[j]*MeV/mm << " mm/MeV "
               << G4endl;
      return g4MatData[j];
    }
  }
  return FindBirksCoefficient(mat);
}

void G4EmSaturation::SetVerbose

(

G4int

)

[inline]

Definition at line 132 of file G4EmSaturation.hh.

Referenced by G4LossTableManager::SetVerbose().

{
  verbose = val;
}

G4double G4EmSaturation::VisibleEnergyDeposition

(

const G4Step *

)

[inline]

Definition at line 137 of file G4EmSaturation.hh.

References G4Track::GetMaterialCutsCouple(), G4Step::GetNonIonizingEnergyDeposit(), G4Track::GetParticleDefinition(), G4Step::GetStepLength(), G4Step::GetTotalEnergyDeposit(), G4Step::GetTrack(), and VisibleEnergyDeposition().

{
  G4Track* track = step->GetTrack();
  return VisibleEnergyDeposition(track->GetParticleDefinition(),
                                 track->GetMaterialCutsCouple(),
                                 step->GetStepLength(),
                                 step->GetTotalEnergyDeposit(),
                                 step->GetNonIonizingEnergyDeposit());
}

G4double G4EmSaturation::VisibleEnergyDeposition	(	const G4ParticleDefinition *	,
		const G4MaterialCutsCouple *	,
		G4double	length,
		G4double	edepTotal,
		G4double	edepNIEL = 0.0
	)

Definition at line 83 of file G4EmSaturation.cc.

References G4MaterialCutsCouple::GetMaterial(), G4ParticleDefinition::GetPDGEncoding(), G4LossTableManager::GetRange(), and G4Proton::Proton().

Referenced by G4Scintillation::PostStepDoIt(), and VisibleEnergyDeposition().

{
  if(edep <= 0.0) { return 0.0; }

  G4double evis = edep;
  G4double bfactor = FindBirksCoefficient(couple->GetMaterial());

  if(bfactor > 0.0) { 

    G4int pdgCode = p->GetPDGEncoding();
    // atomic relaxations for gamma incident
    if(22 == pdgCode) {
      evis /= (1.0 + bfactor*edep/manager->GetRange(electron,edep,couple));

      // energy loss
    } else {

      // protections
      G4double nloss = niel;
      if(nloss < 0.0) nloss = 0.0;
      G4double eloss = edep - nloss;

      // neutrons
      if(2112 == pdgCode || eloss < 0.0 || length <= 0.0) {
        nloss = edep;
        eloss = 0.0;
      }

      // continues energy loss
      if(eloss > 0.0) { eloss /= (1.0 + bfactor*eloss/length); }
 
      // non-ionizing energy loss
      if(nloss > 0.0) {
        if(!proton) { proton = G4Proton::Proton(); }
        G4double escaled = nloss*curRatio;
        G4double range = manager->GetRange(proton,escaled,couple)/curChargeSq; 
        nloss /= (1.0 + bfactor*nloss/range);
      }

      evis = eloss + nloss;
    }
  }
  
  return evis;
}

---

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

• G4EmSaturation.hh
• G4EmSaturation.cc

---