133 for (
size_t j=0; j<msc; ++j) {
137 for (
size_t k=0; k<emp; ++k) {
142 for (
size_t a=0;
a<mod; ++
a) {
144 for (
size_t b=0; b<fmod; ++b) {
152 for (
size_t b=0; b<fmod; ++b) {
233 G4cout <<
"G4LossTableManager::Register G4VEnergyLossProcess : "
288 for (
G4int i=0; i<
n; ++i) {
292 G4cout <<
"G4LossTableManager::Register G4VMultipleScattering : "
304 for (
size_t i=0; i<msc; ++i) {
315 for (
G4int i=0; i<
n; ++i) {
319 G4cout <<
"G4LossTableManager::Register G4VEmProcess : "
331 for (
size_t i=0; i<emp; ++i) {
342 G4cout <<
"G4LossTableManager::Register G4VEmModel : "
352 for (
size_t i=0; i<
n; ++i) {
363 G4cout <<
"G4LossTableManager::Register G4VEmFluctuationModel : "
373 for (
size_t i=0; i<
n; ++i) {
384 if(!p || !part) {
return; }
389 G4cout <<
"G4LossTableManager::RegisterExtraParticle "
412 G4cout <<
"G4LossTableManager::PreparePhysicsTable for "
423 G4cout <<
"====== G4LossTableManager::PreparePhysicsTable start ====="
452 G4cout <<
"G4LossTableManager::PreparePhysicsTable for "
461 G4cout <<
"====== G4LossTableManager::PreparePhysicsTable start ====="
481 G4cout <<
"G4LossTableManager::PreparePhysicsTable for "
491 G4cout <<
"====== G4LossTableManager::PreparePhysicsTable start ====="
520 G4cout <<
"### G4LossTableManager::LocalPhysicsTable() for "
535 G4cout <<
"===== G4LossTableManager::LocalPhysicsTable() for run "
587 G4cout <<
"### G4LossTableManager::LocalPhysicsTable end"
592 G4cout <<
"%%%%% All dEdx and Range tables for worker are ready for run "
605 G4cout <<
"### G4LossTableManager::BuildPhysicsTable() for "
617 G4cout <<
"===== G4LossTableManager::BuildPhysicsTable() for run "
629 G4cout <<
"### G4LossTableManager start initilisation for first particle "
655 G4cout <<
" base particle "
697 G4cout <<
"### G4LossTableManager::BuildPhysicsTable end: "
703 G4cout <<
"%%%%% All dEdx and Range tables are built for master run= "
740 <<
" tables are assigned"
757 G4cout <<
"G4LossTableManager::BuildTables() for "
761 std::vector<G4PhysicsTable*> t_list;
762 std::vector<G4VEnergyLossProcess*> loss_list;
763 std::vector<G4bool> build_flags;
774 for (i=0; i<
n_loss; ++i) {
782 for(
G4int j=0; j<nvec; ++j) {
784 if(ptr == (*pvec)[j]) {
809 loss_list.push_back(p);
810 build_flags.push_back(val);
815 G4int n_dedx = t_list.size();
816 if (0 == n_dedx || !em) {
817 G4cout <<
"G4LossTableManager WARNING: no DEDX processes for "
824 G4cout <<
"G4LossTableManager::BuildTables() start to build range tables"
825 <<
" and the sum of " << n_dedx <<
" processes"
828 <<
" nSubRegions= " << nSubRegions;
877 std::vector<G4PhysicsTable*> listSub;
878 std::vector<G4PhysicsTable*> listCSDA;
880 for (i=0; i<n_dedx; ++i) {
886 if (0 < nSubRegions) {
889 listSub.push_back(dedx);
898 listCSDA.push_back(dedx);
902 if (0 < nSubRegions) {
904 if (1 < listSub.size()) {
927 G4cout <<
"G4LossTableManager::BuildTables: Tables are built for "
949 ed <<
"Energy loss process not found for " << aParticle->
GetParticleName()
951 G4Exception(
"G4LossTableManager::ParticleHaveNoLoss",
"em0001",
981 for (
size_t k=0; k<emp; ++k) {
1056 if(0.0 < v1 && 0.0 < v2 && v2 < 1.e+50) {
1063 }
else if(v1 <= 0.0) {
1097 const std::vector<G4VEnergyLossProcess*>&
1112 const std::vector<G4VMultipleScattering*>&
1168 G4String ss =
"G4LossTableManager::" + tit;
G4EmConfigurator * EmConfigurator()
void SetRandomStep(G4bool val)
void BuildRangeTable(const G4PhysicsTable *dedxTable, G4PhysicsTable *rangeTable, G4bool isIonisation=false)
void SetDEDXTable(G4PhysicsTable *p, G4EmTableType tType)
void SetIntegral(G4bool val)
G4int WorkerVerbose() const
static G4LossTableManager * Instance()
void InitialiseAtomicDeexcitation()
G4VAtomDeexcitation * atomDeexcitation
std::ostringstream G4ExceptionDescription
void SetVerbose(G4int verb)
void SetIonisation(G4bool val)
G4VSubCutProducer * subcutProducer
G4double maxRangeVariation
std::vector< G4VEmModel * > mod_vector
G4PhysicsTable * SubLambdaTable() const
G4ElectronIonPair * emElectronIonPair
void DeRegister(G4VEnergyLossProcess *p)
G4PhysicsTable * RangeTableForLoss() const
void SetLambdaBinning(G4int val)
G4EmSaturation * emSaturation
void push_back(G4PhysicsVector *)
static G4ThreadLocal G4LossTableManager * instance
void AddCollaborativeProcess(G4VEnergyLossProcess *)
void SetStepFunction(G4double v1, G4double v2)
G4PhysicsTable * CSDARangeTable() const
const G4String & GetName() const
void SetNumberOfBins(G4int val)
G4EnergyLossMessenger * GetMessenger()
void CopyTables(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *)
void SetMaxEnergyForCSDARange(G4double val)
std::vector< G4PhysicsTable * > range_vector
G4PhysicsTable * IonisationTableForSubsec() const
G4PhysicsTable * IonisationTable() const
G4ProcessManager * GetProcessManager() const
const std::vector< G4VEmProcess * > & GetEmProcessVector()
G4PhysicsTable * BuildLambdaTable(G4EmTableType tType=fRestricted)
void SetMaxEnergy(G4double val)
static G4PhysicsTable * PreparePhysicsTable(G4PhysicsTable *physTable)
const G4String & GetParticleName() const
std::map< PD, G4VEnergyLossProcess *, std::less< PD > > loss_map
G4EmConfigurator * emConfigurator
void SetVerboseLevel(G4int)
void SetBuildCSDARange(G4bool val)
void SetInitialisationFlag(G4bool flag)
G4VEnergyLossProcess * BuildTables(const G4ParticleDefinition *aParticle)
G4bool BuildCSDARange() const
G4PhysicsTable * LambdaTable() const
void ParticleHaveNoLoss(const G4ParticleDefinition *aParticle)
void SetInverseRangeTable(G4PhysicsTable *p)
const G4ParticleDefinition * SecondaryParticle() const
void SetSubCutProducer(G4VSubCutProducer *)
void SetParameters(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *)
G4GLOB_DLL std::ostream G4cout
G4PhysicsTable * DEDXTable() const
std::vector< G4VEmFluctuationModel * > fmod_vector
G4VEnergyLossProcess * currentLoss
G4ElectronIonPair * ElectronIonPair()
G4EmSaturation * EmSaturation()
void SetMaxEnergyForCSDARange(G4double val)
G4bool startInitialisation
G4int NumberOfSubCutoffRegions() const
void Register(G4VEnergyLossProcess *p)
const G4ParticleDefinition * BaseParticle() const
void SetDEDXBinningForCSDARange(G4int val)
void BuildInverseRangeTable(const G4PhysicsTable *rangeTable, G4PhysicsTable *invRangeTable, G4bool isIonisation=false)
G4PhysicsTable * DEDXTableForSubsec() const
const G4String & GetProcessName() const
const std::vector< G4VEnergyLossProcess * > & GetEnergyLossProcessVector()
void SetSubCutoff(G4bool val, const G4Region *r=0)
void RegisterExtraParticle(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *p)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
std::vector< G4PhysicsTable * > inv_range_vector
void BuildDEDXTable(G4PhysicsTable *dedxTable, const std::vector< G4PhysicsTable * > &)
void SetVerbose(G4int value)
void SetMinSubRange(G4double val)
G4EmParameters * theParameters
G4EmCorrections * emCorrections
void SetMinEnergy(G4double val)
void SetLambdaTable(G4PhysicsTable *p)
G4bool all_tables_are_built
void SetMaxEnergy(G4double val)
const G4String & GetName() const
void SetStepFunction(G4double v1, G4double v2)
G4PhysicsTable * InverseRangeTable() const
void DumpG4BirksCoefficients()
void BuildPhysicsTable(const G4ParticleDefinition *aParticle)
G4LossTableBuilder * tableBuilder
const G4ParticleDefinition * Particle() const
void SetLinearLossLimit(G4double val)
void SetMinEnergy(G4double val)
std::vector< G4PhysicsTable * > dedx_vector
void PrepareModels(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *p)
static G4EmParameters * Instance()
std::vector< G4VEnergyLossProcess * > loss_vector
void PreparePhysicsTable(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *p, G4bool theMaster)
void SetCSDARangeTable(G4PhysicsTable *pRange)
void PrintEWarning(G4String, G4double)
void SetVerbose(G4int val)
void LocalPhysicsTables(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *p)
void SetDEDXBinning(G4int val)
void SetMaxEnergyForMuons(G4double val)
static G4Electron * Electron()
std::vector< PD > part_vector
void SetIntegral(G4bool val)
void SetLossFluctuations(G4bool val)
const G4String & GetName() const
std::vector< G4bool > isActive
void SetSecondaryRangeTable(G4PhysicsTable *p)
std::vector< G4VEmProcess * > emp_vector
G4PhysicsTable * DEDXunRestrictedTable() const
void SetSubLambdaTable(G4PhysicsTable *p)
std::vector< PD > base_part_vector
void ActivateSubCutoff(G4bool val, const G4Region *region=0)
void SetRangeTableForLoss(G4PhysicsTable *p)
std::vector< G4VMultipleScattering * > msc_vector
void SetAtomDeexcitation(G4VAtomDeexcitation *)
G4EnergyLossMessenger * theMessenger
G4PhysicsTable * BuildDEDXTable(G4EmTableType tType=fRestricted)
const std::vector< G4VMultipleScattering * > & GetMultipleScatteringVector()
G4ProcessVector * GetProcessList() const
G4bool stepFunctionActive
std::vector< G4bool > tables_are_built
void SetSplineFlag(G4bool flag)
G4bool IsIonisationProcess() const