89 if(edep <= 0.0) {
return 0.0; }
107 if(nloss < 0.0) { nloss = 0.0; }
111 if(0.0 == p->
GetPDGCharge() || eloss < 0.0 || length <= 0.0) {
117 if(eloss > 0.0) { eloss /= (1.0 + bfactor*eloss/length); }
120 if(nloss > 0.0 &&
proton) {
129 nloss /= (1.0 + bfactor*nloss/range);
132 evis = eloss + nloss;
149 G4cout <<
"### G4EmSaturation::FindG4BirksCoefficient for "
203 ed <<
"Birks constants are not defined for material " << name
204 <<
" ! \n Define Birks constants for the material"
205 <<
" or not apply saturation.";
206 G4Exception(
"G4EmSaturation::InitialiseBirksCoefficient",
"em0088",
219 for (
size_t i=0; i<nelm; ++i) {
220 const G4Element* elm = (*theElementVector)[i];
222 G4double w = Z*Z*theAtomNumDensityVector[i];
237 G4cout <<
"### G4EmSaturation::FindBirksCoefficient Birks coefficient for "
248 G4cout <<
"### Birks coeffitients used in run time" <<
G4endl;
252 << br*
MeV/
mm <<
" mm/MeV" <<
" "
265 G4cout <<
"### Birks coeffitients for Geant4 materials" <<
G4endl;
std::vector< const G4Material * > matPointers
const G4ParticleDefinition * electron
G4IonisParamMat * GetIonisation() const
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
static G4LossTableManager * Instance()
std::vector< G4double > g4MatData
const G4ParticleDefinition * proton
std::vector< G4Element * > G4ElementVector
std::ostringstream G4ExceptionDescription
G4double FindBirksCoefficient(const G4Material *)
void InitialiseG4materials()
const G4String & GetName() const
void InitialiseBirksCoefficient(const G4Material *)
virtual G4double VisibleEnergyDeposition(const G4ParticleDefinition *, const G4MaterialCutsCouple *, G4double length, G4double edepTotal, G4double edepNIEL=0.0)
void SetBirksConstant(G4double value)
virtual ~G4EmSaturation()
const G4Material * curMaterial
std::vector< G4double > effCharges
const G4double w[NPOINTSGL]
std::vector< G4String > g4MatNames
G4int GetPDGEncoding() const
void DumpBirksCoefficients()
const G4ElementVector * GetElementVector() const
static G4NistManager * Instance()
G4double GetBirksConstant() const
const G4double * GetVecNbOfAtomsPerVolume() const
std::vector< G4String > matNames
G4GLOB_DLL std::ostream G4cout
G4double FindG4BirksCoefficient(const G4Material *)
std::vector< G4double > massFactors
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
void DumpG4BirksCoefficients()
static G4ParticleTable * GetParticleTable()
G4LossTableManager * manager
G4double GetAtomicMassAmu(const G4String &symb) const
size_t GetNumberOfElements() const
G4EmSaturation(G4int verb)
G4double GetPDGCharge() const
G4double GetRange(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
const G4Material * GetMaterial() const