109 defaultCutValue(1.0 *
mm),
110 isSetDefaultCutValue(false),
111 fRetrievePhysicsTable(false),
112 fStoredInAscii(true),
113 fIsCheckedForRetrievePhysicsTable(false),
114 fIsRestoredCutValues(false),
115 directoryPhysicsTable(
"."),
118 fDisableCheckParticleList(false)
172 :verboseLevel(right.verboseLevel),
173 defaultCutValue(right.defaultCutValue),
174 isSetDefaultCutValue(right.isSetDefaultCutValue),
175 fRetrievePhysicsTable(right.fRetrievePhysicsTable),
176 fStoredInAscii(right.fStoredInAscii),
177 fIsCheckedForRetrievePhysicsTable(right.fIsCheckedForRetrievePhysicsTable),
178 fIsRestoredCutValues(right.fIsRestoredCutValues),
179 directoryPhysicsTable(right.directoryPhysicsTable),
182 fDisableCheckParticleList(right.fDisableCheckParticleList)
209 if (
this != &right) {
231 if (newParticle == 0)
return;
232 G4Exception(
"G4VUserPhysicsList::AddProcessManager",
234 "This method is obsolete");
243 #ifdef G4MULTITHREADED
245 G4ParticleTable::lockCount++;
262 G4cout <<
"G4VUserPhysicsList::InitializeProcessManager: creating ProcessManager to "
281 G4cout <<
"G4VUserPhysicsList::InitializeProcessManager: copying ProcessManager to "
290 #ifdef G4MULTITHREADED
302 #ifdef G4MULTITHREADED
304 G4ParticleTable::lockCount++;
317 if (pmanager!=0)
delete pmanager;
320 G4cout <<
"G4VUserPhysicsList::RemoveProcessManager: ";
321 G4cout <<
"remove ProcessManager from ";
331 #ifdef G4MULTITHREADED
372 G4cout <<
"G4VUserPhysicsList::SetDefaultCutValue: negative cut values"
373 <<
" :" << value/
mm <<
"[mm]" <<
G4endl;
390 G4cout <<
"G4VUserPhysicsList::SetDefaultCutValue:"
391 <<
"default cut value is changed to :"
405 G4cout <<
"G4VUserPhysicsList::GetCutValue "
406 <<
" : No Default Region " <<
G4endl;
411 "No Default Region");
435 G4cout <<
"G4VUserPhysicsList::SetCutValue "
436 <<
" : No Region of " << rname <<
G4endl;
474 G4cout <<
"G4VUserPhysicsList::SetParticleCuts: negative cut values"
475 <<
" :" << cut/
mm <<
"[mm]"
476 <<
" for "<< particleName <<
G4endl;
487 G4cout <<
"G4VUserPhysicsList::SetParticleCuts "
488 <<
" : No Default Region " <<
G4endl;
491 G4Exception(
"G4VUserPhysicsList::SetParticleCuts ",
493 "No Default Region");
507 G4cout <<
"G4VUserPhysicsList::SetParticleCuts: "
508 <<
" :" << cut/
mm <<
"[mm]"
509 <<
" for "<< particleName <<
G4endl;
531 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable"
532 <<
" Retrieve Cut Table failed !!" <<
G4endl;
535 G4Exception(
"G4VUserPhysicsList::BuildPhysicsTable",
537 "Fail to retrieve Production Cut Table");
541 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable"
542 <<
" Retrieve Cut Table successfully " <<
G4endl;
549 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable"
550 <<
" does not retrieve Cut Table but calculate " <<
G4endl;
570 if( particle!=GammaP &&
587 G4cout <<
"#### G4VUserPhysicsList::BuildPhysicsTable() - BuildPhysicsTable("
596 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable "
597 <<
"Physics table can not be retreived and will be calculated "
606 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable "
607 <<
" Retrieve Physics Table for "
618 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable "
619 <<
"Calculate Physics Table for "
630 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable "
631 <<
" : No Process Manager for "
634 <<
" should be created in your PhysicsList" <<
G4endl;
637 G4Exception(
"G4VUserPhysicsList::BuildPhysicsTable",
639 "No process manager");
650 G4cout <<
"G4VUserPhysicsList::BuildPhysicsTable "
651 <<
" : No Process Vector for "
655 G4Exception(
"G4VUserPhysicsList::BuildPhysicsTable",
657 "No process Vector");
663 G4cout <<
" ProcessManager : " << pManager <<
" ProcessManagerShadow : " << pManagerShadow <<
G4endl;
664 for(
G4int iv1=0;iv1<pVector->
size();iv1++)
665 {
G4cout <<
" " << iv1 <<
" - " << (*pVector)[iv1]->GetProcessName() <<
G4endl; }
666 G4cout <<
"--------------------------------------------------------------" <<
G4endl;
669 for(
G4int iv2=0;iv2<pVectorShadow->
size();iv2++)
670 {
G4cout <<
" " << iv2 <<
" - " << (*pVectorShadow)[iv2]->GetProcessName() <<
G4endl; }
673 for (
G4int j=0; j < pVector->
size(); ++j) {
678 if ( pManagerShadow == pManager )
680 (*pVector)[j]->BuildPhysicsTable(*particle);
684 (*pVector)[j]->BuildWorkerPhysicsTable(*particle);
706 G4cout<<
"G4VUserPhysicsList::PreparePhysicsTable "
707 <<
": No Process Manager for "
710 <<
" should be created in your PhysicsList" <<
G4endl;
713 G4Exception(
"G4VUserPhysicsList::PreparePhysicsTable",
715 "No process manager");
728 G4cout <<
"G4VUserPhysicsList::PreparePhysicsTable "
729 <<
": No Process Vector for "
733 G4Exception(
"G4VUserPhysicsList::PreparePhysicsTable",
735 "No process Vector");
738 for (
G4int j=0; j < pVector->
size(); ++j) {
744 if ( pManagerShadow == pManager )
746 (*pVector)[j]->PreparePhysicsTable(*particle);
750 (*pVector)[j]->PrepareWorkerPhysicsTable(*particle);
776 G4cout <<
"G4VUserPhysicsList::BuildIntegralPhysicsTable "
777 <<
" BuildPhysicsTable is invoked for "
794 if ((idx++ % 4) == 3) {
829 G4Exception(
"G4VUserPhysicsList::StorePhysicsTable",
831 "Fail to store Cut Table");
836 G4cout <<
"G4VUserPhysicsList::StorePhysicsTable "
837 <<
" Store material and cut values successfully" <<
G4endl;
850 for ( j=0; j < pVector->
size(); ++j) {
851 if (!(*pVector)[j]->StorePhysicsTable(particle,dir,ascii)){
852 G4String comment =
"Fail to store physics table for ";
853 comment += (*pVector)[j]->GetProcessName();
855 G4Exception(
"G4VUserPhysicsList::StorePhysicsTable",
889 for ( j=0; j < pVector->
size(); ++j) {
891 (*pVector)[j]->RetrievePhysicsTable(particle,directory,ascii);
896 G4cout <<
"G4VUserPhysicsList::RetrievePhysicsTable "
897 <<
" Fail to retrieve Physics Table for "
898 << (*pVector)[j]->GetProcessName() <<
G4endl;
899 G4cout <<
"Calculate Physics Table for "
903 (*pVector)[j]->BuildPhysicsTable(*particle);
906 for ( j=0; j < pVector->
size(); ++j) {
918 G4cout <<
"G4VUserPhysicsList::SetApplyCuts for " << name <<
G4endl;
976 G4cout <<
"G4VUserPhysicsList::SetVerboseLevel :"
991 G4cout <<
"G4VUserPhysicsList::ResetCuts() is obsolete."
992 <<
" This method gives no effect and you can remove it. "<<
G4endl;
G4double GetCutValue(const G4String &pname) const
void SetDefaultCutValue(G4double newCutValue)
#define G4MT_theMessenger
void BuildIntegralPhysicsTable(G4VProcess *, G4ParticleDefinition *)
G4ProductionCuts * GetProductionCuts() const
void SetApplyCuts(G4bool value, const G4String &name)
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
void DumpCutValuesTableIfRequested()
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
void RemoveProcessManager()
G4PhysicsListHelper * _thePLHelper
void SetProcessManager(G4ProcessManager *aProcessManager)
#define fDisplayThreshold
void SetCutValue(G4double aCut, const G4String &pname)
G4bool fIsRestoredCutValues
void SetEnergyRange(G4double lowedge, G4double highedge)
void PreparePhysicsTable(G4ParticleDefinition *)
G4bool RegisterProcess(G4VProcess *process, G4ParticleDefinition *particle)
G4ProductionCutsTable * fCutsTable
G4double GetProductionCut(G4int index) const
G4bool fRetrievePhysicsTable
#define fIsPhysicsTableBuilt
G4bool fIsCheckedForRetrievePhysicsTable
void SetProductionCut(G4double cut, G4int index=-1)
G4bool GetApplyCuts(const G4String &name) const
void SetCutsForRegion(G4double aCut, const G4String &rname)
G4ParticleTable * theParticleTable
const G4String & GetParticleSubType() const
G4ParticleDefinition * GetGenericIon() const
G4UserPhysicsListMessenger * _theMessenger
G4ProcessManager * GetProcessManager() const
G4int GetInstanceID() const
void InitializeProcessManager()
const G4String & GetParticleName() const
G4bool IsGeneralIon() const
static G4RegionStore * GetInstance()
G4bool _fIsPhysicsTableBuilt
void SetPhysicsTableRetrieved(const G4String &directory="")
void SetParticleCuts(G4double cut, G4ParticleDefinition *particle, G4Region *region=0)
void DumpCutValuesTable(G4int flag=1)
void AddProcessManager(G4ParticleDefinition *newParticle, G4ProcessManager *newManager=0)
G4GLOB_DLL std::ostream G4cout
G4ProcessManager * GetMasterProcessManager() const
static G4PART_DLL G4ThreadLocal G4int slavetotalspace
void SetVerboseLevel(G4int value)
G4bool RetrieveCutsTable(const G4String &directory, G4bool ascii=false)
G4int GetInstanceID() const
void SetVerboseLevel(G4int value)
G4ParticleTable::G4PTblDicIterator * _theParticleIterator
G4VUserPhysicsList & operator=(const G4VUserPhysicsList &)
const G4String & GetProcessName() const
G4bool StorePhysicsTable(const G4String &directory=".")
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4bool isSetDefaultCutValue
G4bool IsShortLived() const
static G4ProductionCutsTable * GetProductionCutsTable()
void UseCoupledTransportation(G4bool vl=true)
void SetApplyCutsFlag(G4bool)
void SetCutsWithDefault()
void DeleteAllParticles()
static G4ParticleTable * GetParticleTable()
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
static G4RUN_DLL G4ThreadLocal T * offset
static G4PhysicsListHelper * GetPhysicsListHelper()
G4bool StoreCutsTable(const G4String &directory, G4bool ascii=false)
static G4RUN_DLL G4VUPLManager subInstanceManager
static const G4VUPLManager & GetSubInstanceManager()
const XML_Char int const XML_Char * value
void ResetCuts()
obsolete methods
G4bool GetApplyCutsFlag() const
G4int CreateSubInstance()
G4bool fDisableCheckParticleList
G4String directoryPhysicsTable
virtual void RetrievePhysicsTable(G4ParticleDefinition *, const G4String &directory, G4bool ascii=false)
virtual ~G4VUserPhysicsList()
G4PTblDicIterator * GetIterator() const
#define theParticleIterator
void SetMasterProcessManager(G4ProcessManager *aNewPM)
G4ProcessVector * GetProcessList() const