G4EmSaturation Class Reference

#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 71 of file G4EmSaturation.hh.

Constructor & Destructor Documentation

G4EmSaturation::G4EmSaturation ( G4int  verb )

explicit

Definition at line 63 of file G4EmSaturation.cc.

{
  verbose = verb;

  nWarnings = nG4Birks = 0;

  electron = nullptr;
  proton   = nullptr;
  nist     = G4NistManager::Instance();
}

Here is the call graph for this function:

G4EmSaturation::~G4EmSaturation ( )

virtual

Definition at line 76 of file G4EmSaturation.cc.

{}

Member Function Documentation

void G4EmSaturation::DumpBirksCoefficients

(

)

Definition at line 232 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;
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4EmSaturation::DumpG4BirksCoefficients

(

)

Definition at line 251 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;
    }
  }
}

G4double G4EmSaturation::FindG4BirksCoefficient

(

const G4Material *

mat

)

Definition at line 155 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;
}

Here is the call graph for this function:

void G4EmSaturation::InitialiseG4Saturation

(

)

Definition at line 139 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(); }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4EmSaturation::SetVerbose ( G4int  val )

inline

Definition at line 134 of file G4EmSaturation.hh.

{
  verbose = val;
}

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

virtual

Definition at line 81 of file G4EmSaturation.cc.

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

  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 || length <= 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

inline

Definition at line 141 of file G4EmSaturation.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/electromagnetic/utils/include/G4EmSaturation.hh
• source/geant4.10.03.p03/source/processes/electromagnetic/utils/src/G4EmSaturation.cc