169 G4cout <<
"G4EmCalculator::GetDEDX: E(MeV)= " << kinEnergy/
MeV
170 <<
" DEDX(MeV/mm)= " << res*
mm/
MeV
174 <<
" isIon= " <<
isIon
193 G4cout <<
" G4EmCalculator::GetRangeFromRestrictedDEDX: E(MeV)= "
195 <<
" range(mm)= " << res/
mm
214 ed <<
"G4EmCalculator::GetCSDARange: CSDA table is not built; "
215 <<
" use UI command: /process/eLoss/CSDARange true";
216 G4Exception(
"G4EmCalculator::GetCSDARange",
"em0077",
225 G4cout <<
" G4EmCalculator::GetCSDARange: E(MeV)= " << kinEnergy/
MeV
226 <<
" range(mm)= " << res/
mm
263 G4cout <<
"G4EmCalculator::GetKinEnergy: Range(mm)= " << range/
mm
264 <<
" KinE(MeV)= " << res/
MeV
310 G4cout <<
"G4EmCalculator::GetXSPerVolume: E(MeV)= " << kinEnergy/
MeV
311 <<
" cross(cm-1)= " << res*
cm
315 G4cout <<
" idx= " << idx <<
" Escaled((MeV)= "
352 if(x > 0.0) { res = 1.0/
x; }
354 G4cout <<
"G4EmCalculator::GetMeanFreePath: E(MeV)= " << kinEnergy/
MeV
355 <<
" MFP(mm)= " << res/
mm
386 G4cout <<
"### G4EmCalculator: Inverse Range Table for "
404 <<
" e(MeV)= " << kinEnergy/
MeV <<
" cut(MeV)= " << cut/
MeV
414 G4cout <<
" ICRU73 ion E(MeV)= " << kinEnergy <<
" ";
416 <<
" DEDX(MeV*cm^2/g)= "
428 <<
" Escaled(MeV)= " << escaled;
432 if(
verbose > 1) {
G4cout <<
" no basePart E(MeV)= " << kinEnergy <<
" "; }
436 <<
" DEDX(MeV*cm^2/g)= "
457 G4cout <<
"At boundary energy(MeV)= " << eth/
MeV
458 <<
" DEDX(MeV/mm)= " << res1*
mm/
MeV
466 if(res1 > 0.0 && escaled > 0.0) {
467 res *= (1.0 + (res0/res1 - 1.0)*eth/escaled);
484 G4cout <<
"After Corrections: DEDX(MeV/mm)= " << res*
mm/
MeV
492 G4cout <<
"Sum: E(MeV)= " << kinEnergy/
MeV
493 <<
" DEDX(MeV/mm)= " << res*
mm/
MeV
495 <<
" cut(MeV)= " << cut/
MeV
499 <<
" isIon=" <<
isIon
518 const std::vector<G4VEnergyLossProcess*> vel =
525 for(
G4int i=0; i<
n; ++i) {
531 dedx +=
ComputeDEDX(kinEnergy,part,(vel[i])->GetProcessName(),mat,cut);
551 const std::vector<G4VEnergyLossProcess*> vel =
565 for(
G4int i=0; i<
n; ++i) {
576 dedx +=
ComputeDEDX(kinEnergy,part,(vel[i])->GetProcessName(),
615 <<
" NuclearDEDX(MeV/mm)= " << res*
mm/
MeV
616 <<
" NuclearDEDX(MeV*cm^2/g)= "
646 G4cout <<
"G4EmCalculator::ComputeXSPerVolume: E(MeV)= " << kinEnergy/
MeV
647 <<
" cross(cm-1)= " << res*
cm
648 <<
" cut(keV)= " << aCut/
keV
685 <<
" cross(barn)= " << res/
barn
687 <<
" Z= " << Z <<
" A= " << A/(
g/
mole) <<
" g/mole"
688 <<
" cut(keV)= " << aCut/
keV
721 <<
" cross(barn)= " << res/
barn
723 <<
" Z= " << Z <<
" shellIdx= " << shellIdx
724 <<
" cut(keV)= " << aCut/
keV
744 if(res > 0.0) { res = 1.0/res; }
777 if(x > 0.0) { mfp = 1.0/
x; }
780 <<
" MFP(mm)= " << mfp/
mm
796 ConvertRangeToEnergy(part, mat, range);
842 G4cout <<
"\n G4EmCalculator::UpdateParticle: isIon 1 "
845 <<
" e= " << kinEnergy <<
G4endl;
875 G4cout <<
"### WARNING: G4EmCalculator::FindParticle fails to find "
899 G4cout <<
"### WARNING: G4EmCalculator::FindMaterial fails to find "
911 if(reg !=
"" && reg !=
"world") {
937 size_t nr = store->size();
939 for(
size_t i=0; i<nr; ++i) {
941 material, ((*store)[i])->GetProductionCuts());
942 if(couple) {
break; }
949 ed <<
"G4EmCalculator::FindCouple: fail for material <"
951 if(region) { ed <<
" and region " << region->
GetName(); }
952 G4Exception(
"G4EmCalculator::FindCouple",
"em0078",
1061 G4cout <<
"G4EmCalculator::FindEmModel WARNING: no particle"
1062 <<
" or materail defined; particle: " << p <<
G4endl;
1071 G4cout <<
"## G4EmCalculator::FindEmModel for " << partname
1073 <<
") and " << processName <<
" at E(MeV)= " << scaledEnergy
1075 if(p != part) {
G4cout <<
" GenericIon is the base particle" <<
G4endl; }
1186 const std::vector<G4VEnergyLossProcess*> v =
1189 for(
G4int i=0; i<
n; ++i) {
1190 if((v[i])->GetProcessName() == processName) {
1208 const std::vector<G4VEmProcess*> v =
1211 for(
G4int i=0; i<
n; ++i) {
1212 if((v[i])->GetProcessName() == processName) {
1230 const std::vector<G4VMultipleScattering*> v =
1233 for(
G4int i=0; i<
n; ++i) {
1234 if((v[i])->GetProcessName() == processName) {
1254 for(
G4int i=0; i<nproc; ++i) {
1255 if(processName == (*pv)[i]->GetProcessName()) {
1273 for(
G4int i=0; i<
n; ++i) {
1274 if((*pv)[i] == proc) {
1309 for(
size_t i=0; i<
nn; ++i) {
G4double ComputeNuclearDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *)
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
G4bool GetProcessActivation(G4VProcess *aProcess) const
G4double GetRangeFromRestricteDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
static G4GenericIon * GenericIon()
void PrintInverseRangeTable(const G4ParticleDefinition *)
const G4PhysicsTable * currentLambda
G4double GetMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=nullptr)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
void FindLambdaTable(const G4ParticleDefinition *, const G4String &processName, G4double kinEnergy, G4int &proctype)
G4PhysicsTable * DEDXTable() const
const G4Material * currentMaterial
std::ostringstream G4ExceptionDescription
const G4Region * FindRegion(const G4String &)
void SetupMaterial(const G4Material *)
static G4ParticleTable * GetParticleTable()
G4VAtomDeexcitation * AtomDeexcitation()
static constexpr double MeV
G4double EffectiveChargeCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
const G4ParticleDefinition * lambdaParticle
void SetKineticEnergy(G4double aEnergy)
static constexpr double keV
static constexpr double mm
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
const G4String & GetName() const
G4VEnergyLossProcess * GetEnergyLossProcess(const G4ParticleDefinition *)
const G4ParticleDefinition * theGenericIon
G4double GetDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
G4bool BuildCSDARange() const
void CheckMaterial(G4int Z)
G4PhysicsTable * RangeTableForLoss() const
G4double ComputeDEDXForCutInRange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *mat, G4double rangecut=DBL_MAX)
const G4String & GetParticleName() const
const G4String & GetParticleType() const
G4double ComputeElectronicDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *mat, G4double cut=DBL_MAX)
G4IonTable * GetIonTable() const
G4double GetPDGCharge() const
const G4ParticleDefinition * BaseParticle() const
G4DynamicParticle dynParticle
const G4ParticleDefinition * FindIon(G4int Z, G4int A)
void SetDynamicMassCharge(G4double massratio, G4double charge2ratio)
G4ionEffectiveCharge * ionEffCharge
G4ProcessVector * GetProcessList() const
void SetCurrentCouple(const G4MaterialCutsCouple *)
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
G4double EffectiveChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4double GetPDGMass() const
G4VProcess * FindProcess(const G4ParticleDefinition *part, const G4String &processName)
const std::vector< G4VEnergyLossProcess * > & GetEnergyLossProcessVector()
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
G4double LowEnergyLimit() const
virtual G4double ComputeShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=nullptr)=0
const G4String & GetName() const
G4VMultipleScattering * FindMscProcess(const G4ParticleDefinition *, const G4String &processName)
static constexpr double g
const G4Material * FindMaterial(const G4String &)
G4EmCorrections * EmCorrections()
G4PhysicsTable * GetCrossSectionTable()
G4String currentProcessName
static constexpr double cm2
const std::vector< G4VEmProcess * > & GetEmProcessVector()
virtual void InitialiseForElement(const G4ParticleDefinition *, G4int Z)
G4double GetShellIonisationCrossSectionPerAtom(const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy)
G4PhysicsTable * InverseRangeTable() const
const G4ParticleDefinition const G4Material *G4double range
void PrintDEDXTable(const G4ParticleDefinition *)
G4PhysicsTable * LambdaTable() const
G4double GetKinEnergy(G4double range, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
G4double ComputeMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
const G4ParticleDefinition * FindParticle(const G4String &)
G4Material * FindOrBuildMaterial(const G4String &name, G4bool isotopes=true, G4bool warning=false)
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
G4VEmModel * SelectModel(G4double kinEnergy, size_t idx)
G4VEnergyLossProcess * currentProcess
double A(double temperature)
G4double ComputeTotalDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, G4double cut=DBL_MAX)
G4Material * FindOrBuildSimpleMaterial(G4int Z, G4bool warning=false)
const G4String & GetProcessName() const
static G4Positron * Positron()
static constexpr double eV
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
virtual G4double ComputeCrossSectionPerShell(const G4ParticleDefinition *, G4int Z, G4int shellIdx, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4VEmModel * currentModel
G4double ComputeGammaAttenuationLength(G4double kinEnergy, const G4Material *)
G4bool UpdateCouple(const G4Material *, G4double cut)
G4double ComputeCrossSectionPerAtom(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4double Z, G4double A, G4double cut=0.0)
static G4Electron * Electron()
void PrintRangeTable(const G4ParticleDefinition *)
std::vector< const G4Material * > localMaterials
virtual void InitialiseForMaterial(const G4ParticleDefinition *, const G4Material *)
G4ProductionCuts * GetProductionCuts() const
G4VEmModel * GetModelByIndex(G4int idx=0, G4bool ver=false) const
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idxRegion) const
G4VEnergyLossProcess * FindEnLossProcess(const G4ParticleDefinition *, const G4String &processName)
const std::vector< G4VMultipleScattering * > & GetMultipleScatteringVector()
G4double ComputeCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
static G4ProductionCutsTable * GetProductionCutsTable()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
G4double GetTransportMeanFreePath(const G4ParticleDefinition *part, G4double kinEnergy)
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
static G4RegionStore * GetInstance()
static constexpr double nm
const G4String & GetName() const
G4ProcessManager * GetProcessManager() const
G4double ComputeCrossSectionPerShell(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4int Z, G4int shellIdx, G4double cut=0.0)
G4double GetCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=nullptr)
void SetFluctuationFlag(G4bool val)
G4PhysicsTable * LambdaTable() const
static constexpr double barn
G4double GetCSDARange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
G4double GetEnergy(const G4ParticleDefinition *aParticle, G4double range, const G4MaterialCutsCouple *couple)
const G4MaterialCutsCouple * currentCouple
G4String currentMaterialName
G4String currentParticleName
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
static constexpr double cm
G4GLOB_DLL std::ostream G4cout
G4LossTableManager * manager
static G4LossTableManager * Instance()
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
G4double GetRangeFromRestricteDEDX(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
G4bool FindEmModel(const G4ParticleDefinition *, const G4String &processName, G4double kinEnergy)
G4VEmProcess * FindDiscreteProcess(const G4ParticleDefinition *, const G4String &processName)
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idx) const
G4double GetRange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=nullptr)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static const G4double reg
static constexpr double mole
std::vector< const G4MaterialCutsCouple * > localCouples
const G4MaterialCutsCouple * FindCouple(const G4Material *, const G4Region *r=nullptr)
G4EmParameters * theParameters
void SetVerbose(G4int val)
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
const G4ParticleDefinition * baseParticle
G4bool UpdateParticle(const G4ParticleDefinition *, G4double kinEnergy)
const G4Element * GetElement(G4int iel) const
G4VEnergyLossProcess * FindEnergyLossProcess(const G4ParticleDefinition *)
const G4ParticleDefinition * currentParticle
virtual G4double GetShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=nullptr)=0
static constexpr double GeV
G4double LowestElectronEnergy() const
G4double ComputeDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=DBL_MAX)
G4double ComputeShellIonisationCrossSectionPerAtom(const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy, const G4Material *mat=nullptr)
const G4Material * cutMaterial
G4double GetCSDARange(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
static G4EmParameters * Instance()
G4double CrossSectionPerVolume(G4double kineticEnergy, const G4MaterialCutsCouple *couple)
size_t GetNumberOfElements() const
G4double ComputeEnergyCutFromRangeCut(G4double range, const G4ParticleDefinition *, const G4Material *)
G4bool ActiveForParticle(const G4ParticleDefinition *part, G4VProcess *proc)
G4double GetDEDX(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
static constexpr double gram
static G4NistManager * Instance()
G4double GetDensity() const