#include <G4EmSaturation.hh>

Public Member Functions
	G4EmSaturation (G4int verb)

virtual	~G4EmSaturation ()

virtual G4double	VisibleEnergyDeposition (const G4ParticleDefinition , const G4MaterialCutsCouple , G4double length, G4double edepTotal, G4double edepNIEL=0.0) const

void	InitialiseG4Saturation ()

G4double	FindG4BirksCoefficient (const G4Material *)

void	DumpBirksCoefficients ()

void	DumpG4BirksCoefficients ()

G4double	VisibleEnergyDepositionAtAStep (const G4Step *) const

void	SetVerbose (G4int)

Detailed Description

Definition at line 70 of file G4EmSaturation.hh.

Constructor & Destructor Documentation

◆ G4EmSaturation()

G4EmSaturation::G4EmSaturation ( G4int verb )

explicit

Definition at line 62 of file G4EmSaturation.cc.

{
  verbose = verb;
 
  nWarnings = nG4Birks = 0;
 
  electron = nullptr;
  proton   = nullptr;
  nist     = G4NistManager::Instance();
}

◆ ~G4EmSaturation()

G4EmSaturation::~G4EmSaturation ( )

virtual

Definition at line 75 of file G4EmSaturation.cc.

76{}

Member Function Documentation

◆ DumpBirksCoefficients()

void G4EmSaturation::DumpBirksCoefficients ( )

Definition at line 236 of file G4EmSaturation.cc.

{
  G4cout << "### Birks coefficients used in run time" << G4endl;
  const G4MaterialTable* mtable = G4Material::GetMaterialTable();
  for(G4int i=0; i<nMaterials; ++i) {
    const G4Material* mat = (*mtable)[i];
    G4double br = mat->GetIonisation()->GetBirksConstant();
    if(br > 0.0) {
      G4cout << "   " << mat->GetName() << "     " 
         << br*MeV/mm << " mm/MeV" << "     "
         << br*mat->GetDensity()*MeV*cm2/g 
         << " g/cm^2/MeV  massFactor=  " << massFactors[i]  
         << " effCharge= " << effCharges[i] << G4endl;
    }
  }
}

Referenced by InitialiseG4Saturation().

◆ DumpG4BirksCoefficients()

void G4EmSaturation::DumpG4BirksCoefficients ( )

Definition at line 255 of file G4EmSaturation.cc.

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

◆ FindG4BirksCoefficient()

G4double G4EmSaturation::FindG4BirksCoefficient ( const G4Material * mat )

Definition at line 159 of file G4EmSaturation.cc.

{
  if(0 == nG4Birks) {  InitialiseG4materials(); }
 
  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 0.0;
}

◆ InitialiseG4Saturation()

void G4EmSaturation::InitialiseG4Saturation ( )

Definition at line 143 of file G4EmSaturation.cc.

{
  nMaterials = G4Material::GetNumberOfMaterials();
  massFactors.resize(nMaterials, 1.0);
  effCharges.resize(nMaterials, 1.0);
 
  if(0 == nG4Birks) {  InitialiseG4materials(); }
 
  for(G4int i=0; i<nMaterials; ++i) {
    InitialiseBirksCoefficient((*G4Material::GetMaterialTable())[i]);
  }
  if(verbose > 0) { DumpBirksCoefficients(); }
}

Referenced by G4EmParameters::SetBirksActive().

◆ SetVerbose()

void G4EmSaturation::SetVerbose ( G4int val )

inline

Definition at line 133 of file G4EmSaturation.hh.

{
  verbose = val;
}

◆ VisibleEnergyDeposition()

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

virtual

Definition at line 80 of file G4EmSaturation.cc.

{
  // no energy deposition
  if(edep <= 0.0) { return 0.0; }
 
  // zero step length may happens only if step limiter process
  // is applied, in that case saturation should not be applied
  if(length <= 0.0) { return edep; }
 
  G4double evis = edep;
  G4double bfactor = couple->GetMaterial()->GetIonisation()->GetBirksConstant();
 
  if(bfactor > 0.0) { 
 
    // atomic relaxations for gamma incident
    if(22 ==  p->GetPDGEncoding()) {
      //G4cout << "%% gamma edep= " << edep/keV << " keV " <<manager << G4endl; 
      evis /= (1.0 + bfactor*edep/
        G4LossTableManager::Instance()->GetRange(electron,edep,couple));
 
      // energy loss
    } else {
 
      // protections
      G4double nloss = std::max(niel, 0.0);
      G4double eloss = edep - nloss;
 
      // neutrons and neutral hadrons
      if(0.0 == p->GetPDGCharge() || eloss < 0.0) {
        nloss = edep;
        eloss = 0.0;
      } else {
 
    // continues energy loss
    eloss /= (1.0 + bfactor*eloss/length); 
      }
      // non-ionizing energy loss
      if(nloss > 0.0) {
        G4int idx = couple->GetMaterial()->GetIndex();
        G4double escaled = nloss*massFactors[idx];
        /*
        G4cout << "%% p edep= " << nloss/keV << " keV  Escaled= " 
               << escaled << " MeV  in " << couple->GetMaterial()->GetName()
               << "  " << p->GetParticleName()
               << G4endl; 
        */
        G4double range = G4LossTableManager::Instance()
          ->GetRange(proton,escaled,couple)/effCharges[idx]; 
        nloss /= (1.0 + bfactor*nloss/range);
      }
      evis = eloss + nloss;
    }
  }
  return evis;
}

Referenced by VisibleEnergyDepositionAtAStep().

◆ VisibleEnergyDepositionAtAStep()

G4double G4EmSaturation::VisibleEnergyDepositionAtAStep ( const G4Step * step ) const