164 G4cout <<
"G4EmCalculator::GetDEDX: E(MeV)= " << kinEnergy/
MeV
165 <<
" DEDX(MeV/mm)= " << res*
mm/
MeV
169 <<
" isIon= " <<
isIon
188 G4cout <<
" G4EmCalculator::GetRangeFromRestrictedDEDX: E(MeV)= "
190 <<
" range(mm)= " << res/
mm
209 ed <<
"G4EmCalculator::GetCSDARange: CSDA table is not built; "
210 <<
" use UI command: /process/eLoss/CSDARange true";
211 G4Exception(
"G4EmCalculator::GetCSDARange",
"em0077",
220 G4cout <<
" G4EmCalculator::GetCSDARange: E(MeV)= " << kinEnergy/
MeV
221 <<
" range(mm)= " << res/
mm
258 G4cout <<
"G4EmCalculator::GetKinEnergy: Range(mm)= " << range/
mm
259 <<
" KinE(MeV)= " << res/
MeV
291 G4cout <<
"G4EmCalculator::GetXSPerVolume: E(MeV)= " << kinEnergy/
MeV
292 <<
" cross(cm-1)= " << res*
cm
296 G4cout <<
" idx= " << idx <<
" Escaled((MeV)= "
332 if(x > 0.0) { res = 1.0/
x; }
334 G4cout <<
"G4EmCalculator::GetMeanFreePath: E(MeV)= " << kinEnergy/
MeV
335 <<
" MFP(mm)= " << res/
mm
366 G4cout <<
"### G4EmCalculator: Inverse Range Table for "
384 <<
" e(MeV)= " << kinEnergy/
MeV <<
" cut(MeV)= " << cut/
MeV
394 G4cout <<
" ICRU73 ion E(MeV)= " << kinEnergy <<
" ";
396 <<
" DEDX(MeV*cm^2/g)= "
408 <<
" Escaled(MeV)= " << escaled;
412 if(
verbose > 1) {
G4cout <<
" no basePart E(MeV)= " << kinEnergy <<
" "; }
416 <<
" DEDX(MeV*cm^2/g)= "
437 G4cout <<
"At boundary energy(MeV)= " << eth/
MeV
438 <<
" DEDX(MeV/mm)= " << res1*
mm/
MeV
446 if(res1 > 0.0 && escaled > 0.0) {
447 res *= (1.0 + (res0/res1 - 1.0)*eth/escaled);
464 G4cout <<
"After Corrections: DEDX(MeV/mm)= " << res*
mm/
MeV
472 G4cout <<
"Sum: E(MeV)= " << kinEnergy/
MeV
473 <<
" DEDX(MeV/mm)= " << res*
mm/
MeV
475 <<
" cut(MeV)= " << cut/
MeV
479 <<
" isIon=" <<
isIon
498 const std::vector<G4VEnergyLossProcess*> vel =
505 for(
G4int i=0; i<
n; ++i) {
511 dedx +=
ComputeDEDX(kinEnergy,part,(vel[i])->GetProcessName(),mat,cut);
531 const std::vector<G4VEnergyLossProcess*> vel =
545 for(
G4int i=0; i<
n; ++i) {
556 dedx +=
ComputeDEDX(kinEnergy,part,(vel[i])->GetProcessName(),
588 <<
" NuclearDEDX(MeV/mm)= " << res*
mm/
MeV
589 <<
" NuclearDEDX(MeV*cm^2/g)= "
618 G4cout <<
"G4EmCalculator::ComputeXSPerVolume: E(MeV)= " << kinEnergy/
MeV
619 <<
" cross(cm-1)= " << res*
cm
620 <<
" cut(keV)= " << cut/
keV
655 <<
" cross(barn)= " << res/
barn
657 <<
" Z= " << Z <<
" A= " << A/(
g/
mole) <<
" g/mole"
689 <<
" cross(barn)= " << res/
barn
691 <<
" Z= " << Z <<
" shellIdx= " << shellIdx <<
G4endl;
710 if(res > 0.0) { res = 1.0/res; }
743 if(x > 0.0) { mfp = 1.0/
x; }
746 <<
" MFP(mm)= " << mfp/
mm
762 ConvertRangeToEnergy(part, mat, range);
808 G4cout <<
"\n G4EmCalculator::UpdateParticle: isIon 1 "
811 <<
" e= " << kinEnergy <<
G4endl;
841 G4cout <<
"### WARNING: G4EmCalculator::FindParticle fails to find "
865 G4cout <<
"### WARNING: G4EmCalculator::FindMaterial fails to find "
877 if(reg !=
"" && reg !=
"world") {
903 size_t nr = store->size();
905 for(
size_t i=0; i<nr; ++i) {
907 material, ((*store)[i])->GetProductionCuts());
908 if(couple) {
break; }
915 ed <<
"G4EmCalculator::FindCouple: fail for material <"
917 if(region) { ed <<
" and region " << region->
GetName(); }
918 G4Exception(
"G4EmCalculator::FindCouple",
"em0078",
1021 G4cout <<
"G4EmCalculator::FindEmModel WARNING: no particle"
1022 <<
" or materail defined; particle: " << p <<
G4endl;
1031 G4cout <<
"## G4EmCalculator::FindEmModel for " << partname
1033 <<
") and " << processName <<
" at E(MeV)= " << scaledEnergy
1035 if(p != part) {
G4cout <<
" GenericIon is the base particle" <<
G4endl; }
1146 const std::vector<G4VEnergyLossProcess*> v =
1149 for(
G4int i=0; i<
n; ++i) {
1150 if((v[i])->GetProcessName() == processName) {
1168 const std::vector<G4VEmProcess*> v =
1171 for(
G4int i=0; i<
n; ++i) {
1172 if((v[i])->GetProcessName() == processName) {
1190 const std::vector<G4VMultipleScattering*> v =
1193 for(
G4int i=0; i<
n; ++i) {
1194 if((v[i])->GetProcessName() == processName) {
1214 for(
G4int i=0; i<nproc; ++i) {
1215 if(processName == (*pv)[i]->GetProcessName()) {
1233 for(
G4int i=0; i<
n; ++i) {
1234 if((*pv)[i] == proc) {
1269 for(
size_t i=0; i<
nn; ++i) {
G4PhysicsTable * LambdaTable() const
G4VEmProcess * FindDiscreteProcess(const G4ParticleDefinition *, const G4String &processName)
G4ProductionCuts * GetProductionCuts() const
const G4PhysicsTable * currentLambda
G4double LowEnergyLimit() const
void FindLambdaTable(const G4ParticleDefinition *, const G4String &processName, G4double kinEnergy)
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4Material * FindOrBuildMaterial(const G4String &name, G4bool isotopes=true, G4bool warning=false)
const G4ParticleDefinition * theGenericIon
G4VEmModel * SelectModel(G4double kinEnergy, size_t idx)
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
virtual void InitialiseForElement(const G4ParticleDefinition *, G4int Z)
const G4String & GetName() const
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
static G4LossTableManager * Instance()
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
void SetDynamicMassCharge(G4double massratio, G4double charge2ratio)
std::ostringstream G4ExceptionDescription
G4double GetDEDX(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
const G4Material * FindMaterial(const G4String &)
G4Material * FindOrBuildSimpleMaterial(G4int Z, G4bool warning=false)
const G4Material * cutMaterial
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
std::vector< const G4MaterialCutsCouple * > localCouples
const G4ParticleDefinition * currentParticle
G4PhysicsTable * RangeTableForLoss() const
G4double EffectiveChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4bool BuildCSDARange() const
const G4String & GetName() const
void PrintRangeTable(const G4ParticleDefinition *)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
virtual G4double GetShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=0)=0
virtual void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4bool UpdateCouple(const G4Material *, G4double cut)
G4double GetDensity() const
G4double ComputeTotalDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, G4double cut=DBL_MAX)
G4double ComputeElectronicDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *mat, G4double cut=DBL_MAX)
G4PhysicsTable * GetCrossSectionTable()
G4double EffectiveChargeCorrection(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4ionEffectiveCharge * ionEffCharge
const G4MaterialCutsCouple * currentCouple
G4double GetCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=0)
G4double GetMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, const G4Region *r=0)
void SetupMaterial(const G4Material *)
void PrintInverseRangeTable(const G4ParticleDefinition *)
G4double GetRange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
const G4Element * GetElement(G4int iel) const
G4ProcessManager * GetProcessManager() const
const std::vector< G4VEmProcess * > & GetEmProcessVector()
G4DynamicParticle dynParticle
static G4NistManager * Instance()
const G4String & GetParticleName() const
G4double ComputeShellIonisationCrossSectionPerAtom(const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy, const G4Material *mat=0)
G4double GetCSDARange(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
G4double GetDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
static const G4double reg
static G4RegionStore * GetInstance()
G4double GetCSDARange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
G4double GetShellIonisationCrossSectionPerAtom(const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy)
G4PhysicsTable * LambdaTable() const
G4VEnergyLossProcess * FindEnergyLossProcess(const G4ParticleDefinition *)
G4IonTable * GetIonTable() const
virtual void InitialiseForMaterial(const G4ParticleDefinition *, const G4Material *)
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
G4VEnergyLossProcess * currentProcess
const G4MaterialCutsCouple * FindCouple(const G4Material *, const G4Region *r=0)
G4PhysicsTable * DEDXTable() const
G4double ComputeDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=DBL_MAX)
G4double NuclearDEDX(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy, G4bool fluct=true)
const G4ParticleDefinition * baseParticle
G4EmCorrections * EmCorrections()
G4bool FindEmModel(const G4ParticleDefinition *, const G4String &processName, G4double kinEnergy)
G4double ComputeNuclearDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *)
G4double GetRangeFromRestricteDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
const G4String & GetParticleType() const
const G4ParticleDefinition * BaseParticle() const
G4double GetEnergy(const G4ParticleDefinition *aParticle, G4double range, const G4MaterialCutsCouple *couple)
const G4String & GetProcessName() const
void CheckMaterial(G4int Z)
G4bool ActiveForParticle(const G4ParticleDefinition *part, G4VProcess *proc)
const std::vector< G4VEnergyLossProcess * > & GetEnergyLossProcessVector()
G4VMultipleScattering * FindMscProcess(const G4ParticleDefinition *, const G4String &processName)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
void SetKineticEnergy(G4double aEnergy)
virtual G4double ComputeShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=0)=0
G4double GetRangeFromRestricteDEDX(const G4ParticleDefinition *aParticle, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
std::vector< const G4Material * > localMaterials
G4VEnergyLossProcess * FindEnLossProcess(const G4ParticleDefinition *, const G4String &processName)
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
G4VEnergyLossProcess * GetEnergyLossProcess(const G4ParticleDefinition *)
static G4ProductionCutsTable * GetProductionCutsTable()
static G4Positron * Positron()
static G4GenericIon * GenericIon()
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idxRegion) const
const G4ParticleDefinition * FindParticle(const G4String &)
G4PhysicsTable * InverseRangeTable() const
G4double GetPDGMass() const
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
static G4ParticleTable * GetParticleTable()
G4double ComputeMeanFreePath(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
G4double ComputeCrossSectionPerShell(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4int Z, G4int shellIdx, G4double cut=0.0)
G4double GetKinEnergy(G4double range, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
const G4Region * FindRegion(const G4String &)
G4String currentMaterialName
const G4ParticleDefinition * FindIon(G4int Z, G4int A)
G4String currentProcessName
G4String currentParticleName
const G4double x[NPOINTSGL]
G4double ComputeCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
virtual G4double ComputeCrossSectionPerShell(const G4ParticleDefinition *, G4int Z, G4int shellIdx, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4double ComputeGammaAttenuationLength(G4double kinEnergy, const G4Material *)
G4bool GetProcessActivation(G4VProcess *aProcess) const
static G4Electron * Electron()
G4double ComputeCrossSectionPerAtom(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4double Z, G4double A, G4double cut=0.0)
const G4String & GetName() const
G4VProcess * FindProcess(const G4ParticleDefinition *part, const G4String &processName)
size_t GetNumberOfElements() const
void SetVerbose(G4int val)
void PrintDEDXTable(const G4ParticleDefinition *)
G4VAtomDeexcitation * AtomDeexcitation()
const G4ParticleDefinition * lambdaParticle
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
G4double ComputeEnergyCutFromRangeCut(G4double range, const G4ParticleDefinition *, const G4Material *)
const G4Material * currentMaterial
G4VEmModel * currentModel
G4double GetPDGCharge() const
G4bool UpdateParticle(const G4ParticleDefinition *, G4double kinEnergy)
G4LossTableManager * manager
G4double ComputeDEDXForCutInRange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *mat, G4double rangecut=DBL_MAX)
const std::vector< G4VMultipleScattering * > & GetMultipleScatteringVector()
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idx) const
G4ProcessVector * GetProcessList() const