IORTPhysicsList Class Reference

#include <IORTPhysicsList.hh>

Inheritance diagram for IORTPhysicsList:

Public Member Functions
	IORTPhysicsList ()
virtual	~IORTPhysicsList ()
void	ConstructParticle ()
void	SetCuts ()
void	SetCutForGamma (G4double)
void	SetCutForElectron (G4double)
void	SetCutForPositron (G4double)
void	SetDetectorCut (G4double cut)
void	AddPhysicsList (const G4String &name)
void	ConstructProcess ()
void	AddStepMax ()
IORTStepMax *	GetStepMaxProcess ()
void	AddPackage (const G4String &name)
Public Member Functions inherited from G4VModularPhysicsList
	G4VModularPhysicsList ()
virtual	~G4VModularPhysicsList ()
void	RegisterPhysics (G4VPhysicsConstructor *)
const G4VPhysicsConstructor *	GetPhysics (G4int index) const
const G4VPhysicsConstructor *	GetPhysics (const G4String &name) const
const G4VPhysicsConstructor *	GetPhysicsWithType (G4int physics_type) const
void	ReplacePhysics (G4VPhysicsConstructor *)
void	RemovePhysics (G4VPhysicsConstructor *)
void	RemovePhysics (G4int type)
void	RemovePhysics (const G4String &name)
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const
Public Member Functions inherited from G4VUserPhysicsList
	G4VUserPhysicsList ()
virtual	~G4VUserPhysicsList ()
	G4VUserPhysicsList (const G4VUserPhysicsList &)
G4VUserPhysicsList &	operator= (const G4VUserPhysicsList &)
void	Construct ()
void	UseCoupledTransportation (G4bool vl=true)
void	SetDefaultCutValue (G4double newCutValue)
G4double	GetDefaultCutValue () const
void	BuildPhysicsTable ()
void	PreparePhysicsTable (G4ParticleDefinition *)
void	BuildPhysicsTable (G4ParticleDefinition *)
G4bool	StorePhysicsTable (const G4String &directory=".")
G4bool	IsPhysicsTableRetrieved () const
G4bool	IsStoredInAscii () const
const G4String &	GetPhysicsTableDirectory () const
void	SetPhysicsTableRetrieved (const G4String &directory="")
void	SetStoredInAscii ()
void	ResetPhysicsTableRetrieved ()
void	ResetStoredInAscii ()
void	DumpList () const
void	DumpCutValuesTable (G4int flag=1)
void	DumpCutValuesTableIfRequested ()
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
void	SetCutsWithDefault ()
void	SetCutValue (G4double aCut, const G4String &pname)
G4double	GetCutValue (const G4String &pname) const
void	SetCutValue (G4double aCut, const G4String &pname, const G4String &rname)
void	SetParticleCuts (G4double cut, G4ParticleDefinition *particle, G4Region *region=0)
void	SetParticleCuts (G4double cut, const G4String &particleName, G4Region *region=0)
void	SetCutsForRegion (G4double aCut, const G4String &rname)
void	ResetCuts ()
	obsolete methods More...
void	SetApplyCuts (G4bool value, const G4String &name)
G4bool	GetApplyCuts (const G4String &name) const
void	RemoveProcessManager ()
void	AddProcessManager (G4ParticleDefinition *newParticle, G4ProcessManager *newManager=0)
void	CheckParticleList ()
void	DisableCheckParticleList ()
G4int	GetInstanceID () const
void	InitializeWorker ()

Additional Inherited Members
Static Public Member Functions inherited from G4VModularPhysicsList
static const G4VMPLManager &	GetSubInstanceManager ()
Static Public Member Functions inherited from G4VUserPhysicsList
static const G4VUPLManager &	GetSubInstanceManager ()
Protected Types inherited from G4VModularPhysicsList
typedef G4VMPLData::G4PhysConstVectorData	G4PhysConstVector
Protected Member Functions inherited from G4VModularPhysicsList
	G4VModularPhysicsList (const G4VModularPhysicsList &)
G4VModularPhysicsList &	operator= (const G4VModularPhysicsList &)
Protected Member Functions inherited from G4VUserPhysicsList
void	AddTransportation ()
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
void	BuildIntegralPhysicsTable (G4VProcess *, G4ParticleDefinition *)
virtual void	RetrievePhysicsTable (G4ParticleDefinition *, const G4String &directory, G4bool ascii=false)
void	InitializeProcessManager ()
Protected Attributes inherited from G4VModularPhysicsList
G4int	verboseLevel
G4int	g4vmplInstanceID
Protected Attributes inherited from G4VUserPhysicsList
G4ParticleTable *	theParticleTable
G4int	verboseLevel
G4double	defaultCutValue
G4bool	isSetDefaultCutValue
G4ProductionCutsTable *	fCutsTable
G4bool	fRetrievePhysicsTable
G4bool	fStoredInAscii
G4bool	fIsCheckedForRetrievePhysicsTable
G4bool	fIsRestoredCutValues
G4String	directoryPhysicsTable
G4bool	fDisableCheckParticleList
G4int	g4vuplInstanceID
Static Protected Attributes inherited from G4VModularPhysicsList
static G4RUN_DLL G4VMPLManager	G4VMPLsubInstanceManager
Static Protected Attributes inherited from G4VUserPhysicsList
static G4RUN_DLL G4VUPLManager	subInstanceManager

Detailed Description

Definition at line 51 of file IORTPhysicsList.hh.

Constructor & Destructor Documentation

IORTPhysicsList::IORTPhysicsList

(

)

Definition at line 110 of file IORTPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, G4LossTableManager::Instance(), python.hepunit::mm, and G4VModularPhysicsList::SetVerboseLevel().

                                 : G4VModularPhysicsList()
{
  G4LossTableManager::Instance();
  defaultCutValue = 0.01 *mm; //1.*mm;
  cutForGamma     = defaultCutValue;
  cutForElectron  = defaultCutValue;
  cutForPositron  = defaultCutValue;

  helIsRegistered  = false;
  bicIsRegistered  = false;
  biciIsRegistered = false;
  locIonIonInelasticIsRegistered = false;
  radioactiveDecayIsRegistered = false;

  stepMaxProcess  = 0;

  pMessenger = new IORTPhysicsListMessenger(this);

  SetVerboseLevel(1);

  // EM physics
  emPhysicsList = new G4EmStandardPhysics_option3(1);
  emName = G4String("emstandard_opt3");

  // Decay physics and all particles
  decPhysicsList = new G4DecayPhysics();
}

IORTPhysicsList::~IORTPhysicsList ( )

virtual

Definition at line 139 of file IORTPhysicsList.cc.

{
  delete pMessenger;
  delete emPhysicsList;
  delete decPhysicsList;
  for(size_t i=0; i<hadronPhys.size(); i++) {delete hadronPhys[i];}
}

Member Function Documentation

void IORTPhysicsList::AddPackage

(

const G4String &

name

)

void IORTPhysicsList::AddPhysicsList

(

const G4String &

name

)

Definition at line 177 of file IORTPhysicsList.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), and G4VModularPhysicsList::verboseLevel.

Referenced by IORTPhysicsListMessenger::SetNewValue().

{

    if (verboseLevel>1) {
    G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
    }
    if (name == emName) return;

    /////////////////////////////////////////////////////////////////////////////
    //   ELECTROMAGNETIC MODELS
    /////////////////////////////////////////////////////////////////////////////
    if (name == "standard_opt3") {
    emName = name;
    delete emPhysicsList;
    emPhysicsList = new G4EmStandardPhysics_option3();
    G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
    G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmStandardPhysics_option3" << G4endl;


    } else if (name == "LowE_Livermore") {   
    emName = name;
    delete emPhysicsList;
    emPhysicsList = new G4EmLivermorePhysics();
    G4RunManager::GetRunManager()-> PhysicsHasBeenModified();
    G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmLivermorePhysics" << G4endl;

    } else if (name == "LowE_Penelope") {    
    emName = name;
    delete emPhysicsList;
    emPhysicsList = new G4EmPenelopePhysics();
    G4RunManager::GetRunManager()-> PhysicsHasBeenModified();
    G4cout << "THE FOLLOWING ELECTROMAGNETIC PHYSICS LIST HAS BEEN ACTIVATED: G4EmPenelopePhysics" << G4endl;

    /////////////////////////////////////////////////////////////////////////////
    //   HADRONIC MODELS
    /////////////////////////////////////////////////////////////////////////////
    } else if (name == "Elastic")
    {
    if(!helIsRegistered) 
    {
        G4cout << "THE FOLLOWING HADRONIC ELASTIC PHYSICS LIST HAS BEEN ACTIVATED: G4HadronElasticPhysics()" << G4endl;
        hadronPhys.push_back( new G4HadronElasticPhysics());
        helIsRegistered = true;
    }
    else  G4cout << "AN ELASTIC PHYSICS HAS BEEN ALREADY ACTIVATED!" << G4endl;
    }
    else if (name == "DElastic")
    {
    if(!helIsRegistered) 
    {
        hadronPhys.push_back( new G4HadronDElasticPhysics());
        helIsRegistered = true;
    }
    else  G4cout << "AN ELASTIC PHYSICS HAS BEEN ALREADY ACTIVATED!" << G4endl;

    }
    else if (name == "HElastic")
    {
    if(!helIsRegistered) 
    {
        hadronPhys.push_back( new G4HadronHElasticPhysics());
        helIsRegistered = true;
    }
    else  G4cout << "AN ELASTIC PHYSICS HAS BEEN ALREADY ACTIVATED!" << G4endl;

    }
    else if (name == "Em_extra_physics")
    {
        hadronPhys.push_back( new G4EmExtraPhysics());
    }
    else if (name == "Stopping_physics")
    {
        hadronPhys.push_back( new G4StoppingPhysics());
    }
    else if (name == "Neutron_tracking_cut")
    {
        hadronPhys.push_back( new G4NeutronTrackingCut());
    }
    else if (name == "Hadron_QGSP_BIC")
    {
        hadronPhys.push_back( new G4HadronPhysicsQGSP_BIC());
       // helIsRegistered = true;
    }
    else if (name == "Hadron_QBBC") 
    {
        hadronPhys.push_back(new G4HadronInelasticQBBC());
        //bicIsRegistered = true;
        G4cout << "THE FOLLOWING HADRONIC INELASTIC PHYSICS LIST HAS BEEN ACTIVATED: G4HadronInelasticQBBC()" << G4endl;
    } 

    else if (name == "binary") 
    {
    hadronPhys.push_back(new G4HadronInelasticQBBC());
    //bicIsRegisted = true;
    G4cout << "THE FOLLOWING HADRONIC INELASTIC PHYSICS LIST HAS BEEN ACTIVATED: G4HadronInelasticQBBC()" << G4endl;
    }
    
    else if (name == "binary_ion")
    {
        hadronPhys.push_back(new G4IonBinaryCascadePhysics());
        //biciIsRegistered = true;
    }
/*
    else if (name == "local_ion_ion_inelastic")
    {
        hadronPhys.push_back(new LocalIonIonInelasticPhysic());
        locIonIonInelasticIsRegistered = true;
    }
    else if (name == "local_incl_ion_ion_inelastic")
    {
        hadronPhys.push_back(new LocalINCLIonIonInelasticPhysic());
        locIonIonInelasticIsRegistered = true;
    } 
*/
    else if (name == "decay")
    {
    hadronPhys.push_back(new G4DecayPhysics());
    //radioactiveDecayIsRegistered = true;
    }
    else if (name == "radioactive_decay" && !radioactiveDecayIsRegistered )
    {
    hadronPhys.push_back(new G4RadioactiveDecayPhysics());
    radioactiveDecayIsRegistered = true;

    // The following is the construction of the QGSP_BIC_EMY Reference physics list 
    // reconstructed here like a builder: it should be identical to the
    // one contained inside the $G4INSTALL/physics_lists/lists folder
    }
    else if (name == "QGSP_BIC_EMY")
    {
    AddPhysicsList("emstandard_opt3");
    hadronPhys.push_back( new G4EmExtraPhysics());
    hadronPhys.push_back( new G4HadronElasticPhysics());
    hadronPhys.push_back( new G4StoppingPhysics());
    hadronPhys.push_back( new G4IonBinaryCascadePhysics());
    hadronPhys.push_back( new G4NeutronTrackingCut());
    hadronPhys.push_back( new G4HadronPhysicsQGSP_BIC());
    hadronPhys.push_back( new G4DecayPhysics());

    }
    else {

    G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
        << " is not defined"
        << G4endl;
    }
}

void IORTPhysicsList::AddStepMax

(

)

Definition at line 326 of file IORTPhysicsList.cc.

References G4ProcessManager::AddDiscreteProcess(), G4ParticleDefinition::GetProcessManager(), IORTStepMax::IsApplicable(), and theParticleIterator.

Referenced by ConstructProcess().

{
  // Step limitation seen as a process
  stepMaxProcess = new IORTStepMax();

  theParticleIterator->reset();
  while ((*theParticleIterator)()){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();

    if (stepMaxProcess->IsApplicable(*particle) && pmanager)
      {
    pmanager ->AddDiscreteProcess(stepMaxProcess);
      }
  }
}

void IORTPhysicsList::ConstructParticle ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 148 of file IORTPhysicsList.cc.

References G4VPhysicsConstructor::ConstructParticle().

{
  decPhysicsList->ConstructParticle();
}

void IORTPhysicsList::ConstructProcess ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 154 of file IORTPhysicsList.cc.

References G4EmConfigurator::AddModels(), AddStepMax(), G4VUserPhysicsList::AddTransportation(), and G4VPhysicsConstructor::ConstructProcess().

{
  // transportation
  AddTransportation();

  // electromagnetic physics list
  emPhysicsList->ConstructProcess();
  em_config.AddModels();

  // decay physics list
  decPhysicsList->ConstructProcess();

  // hadronic physics lists
  for(size_t i=0; i<hadronPhys.size(); i++) {
    hadronPhys[i] -> ConstructProcess();
  }

  // step limitation (as a full process)
  //
  AddStepMax();
}

IORTStepMax* IORTPhysicsList::GetStepMaxProcess ( )

inline

Definition at line 69 of file IORTPhysicsList.hh.

{return stepMaxProcess;};

void IORTPhysicsList::SetCutForElectron

(

G4double

cut

)

Definition at line 371 of file IORTPhysicsList.cc.

References G4Electron::Electron(), and G4VUserPhysicsList::SetParticleCuts().

Referenced by IORTPhysicsListMessenger::SetNewValue().

{
  cutForElectron = cut;
  SetParticleCuts(cutForElectron, G4Electron::Electron());
}

void IORTPhysicsList::SetCutForGamma

(

G4double

cut

)

Definition at line 364 of file IORTPhysicsList.cc.

References G4Gamma::Gamma(), and G4VUserPhysicsList::SetParticleCuts().

Referenced by IORTPhysicsListMessenger::SetNewValue().

{
  cutForGamma = cut;
  SetParticleCuts(cutForGamma, G4Gamma::Gamma());
}

void IORTPhysicsList::SetCutForPositron

(

G4double

cut

)

Definition at line 378 of file IORTPhysicsList.cc.

References G4Positron::Positron(), and G4VUserPhysicsList::SetParticleCuts().

Referenced by IORTPhysicsListMessenger::SetNewValue().

{
  cutForPositron = cut;
  SetParticleCuts(cutForPositron, G4Positron::Positron());
}

void IORTPhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 344 of file IORTPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, G4VUserPhysicsList::DumpCutValuesTable(), G4BestUnit, G4cout, G4endl, G4VUserPhysicsList::SetCutValue(), SetDetectorCut(), and G4VModularPhysicsList::verboseLevel.

{

  if (verboseLevel >0){
    G4cout << "PhysicsList::SetCuts:";
    G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
  }

  // set cut values for gamma at first and for e- second and next for e+,
  // because some processes for e+/e- need cut values for gamma
  SetCutValue(cutForGamma, "gamma");
  SetCutValue(cutForElectron, "e-");
  SetCutValue(cutForPositron, "e+");

  // Set cuts for detector
  SetDetectorCut(defaultCutValue); 
  if (verboseLevel>0) DumpCutValuesTable();
}

void IORTPhysicsList::SetDetectorCut

(

G4double

cut

)

Definition at line 384 of file IORTPhysicsList.cc.

References G4ProductionCuts::GetIndex(), G4RegionStore::GetInstance(), and G4RegionStore::GetRegion().

Referenced by SetCuts().

{
  G4String regionName = "DetectorLog";
  G4Region* region = G4RegionStore::GetInstance()->GetRegion(regionName);

  G4ProductionCuts* cuts = new G4ProductionCuts ;
  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("gamma"));
  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("e-"));
  cuts -> SetProductionCut(cut,G4ProductionCuts::GetIndex("e+"));
  region -> SetProductionCuts(cuts);
}

---

The documentation for this class was generated from the following files:

• IORTPhysicsList.hh
• IORTPhysicsList.cc