88 if(edep <= 0.0) {
return 0.0; }
98 evis /= (1.0 + bfactor*edep/
109 if(0.0 == p->
GetPDGCharge() || eloss < 0.0 || length <= 0.0) {
115 eloss /= (1.0 + bfactor*eloss/length);
129 nloss /= (1.0 + bfactor*nloss/
range);
131 evis = eloss + nloss;
165 G4cout <<
"### G4EmSaturation::FindG4BirksCoefficient for "
183 G4Exception(
"G4EmSaturation::InitialiseBirksCoefficient",
"em0001",
194 if(curBirks == 0.0) {
204 if(curBirks == 0.0) {
return; }
213 for (
size_t i=0; i<nelm; ++i) {
214 const G4Element* elm = (*theElementVector)[i];
216 G4double w = Z*Z*theAtomNumDensityVector[i];
234 G4cout <<
"### Birks coefficients used in run time" <<
G4endl;
241 << br*
MeV/
mm <<
" mm/MeV" <<
" "
254 G4cout <<
"### Birks coefficients for Geant4 materials" <<
G4endl;
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
G4EmSaturation(G4int verb)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
void DumpG4BirksCoefficients()
static size_t GetNumberOfMaterials()
static G4ParticleTable * GetParticleTable()
static constexpr double MeV
static G4MaterialTable * GetMaterialTable()
static constexpr double mm
const G4ParticleDefinition * electron
const G4ParticleDefinition * proton
const G4ElementVector * GetElementVector() const
G4double GetPDGCharge() const
const G4double * GetVecNbOfAtomsPerVolume() const
void SetBirksConstant(G4double value)
static std::vector< G4double > effCharges
virtual G4double VisibleEnergyDeposition(const G4ParticleDefinition *, const G4MaterialCutsCouple *, G4double length, G4double edepTotal, G4double edepNIEL=0.0) const
G4double GetAtomicMassAmu(const G4String &symb) const
static constexpr double proton_mass_c2
void InitialiseBirksCoefficient(const G4Material *)
const G4String & GetName() const
static constexpr double g
static std::vector< G4double > massFactors
static constexpr double cm2
const G4ParticleDefinition const G4Material *G4double range
G4double GetBirksConstant() const
G4double GetRange(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
void DumpBirksCoefficients()
std::vector< G4Material * > G4MaterialTable
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
static std::vector< G4String > g4MatNames
std::vector< G4Element * > G4ElementVector
G4IonisParamMat * GetIonisation() const
static std::vector< G4double > g4MatData
virtual ~G4EmSaturation()
G4int GetPDGEncoding() const
G4GLOB_DLL std::ostream G4cout
static G4LossTableManager * Instance()
void InitialiseG4Saturation()
const G4Material * GetMaterial() const
size_t GetNumberOfElements() const
void InitialiseG4materials()
static G4NistManager * Instance()
G4double GetDensity() const
G4double FindG4BirksCoefficient(const G4Material *)