#include <WLSPhysicsList.hh>

Inheritance diagram for WLSPhysicsList:

Public Member Functions
	WLSPhysicsList (G4String)

virtual	~WLSPhysicsList ()

void	SetCuts ()

void	SetCutForGamma (G4double)

void	SetCutForElectron (G4double)

void	SetCutForPositron (G4double)

void	SetStepMax (G4double)

WLSStepMax *	GetStepMaxProcess ()

void	AddStepMax ()

void	RemoveFromPhysicsList (const G4String &)
	Remove specific physics from physics list. More...

void	ClearPhysics ()
	Make sure that the physics list is empty. More...

virtual void	ConstructParticle ()

virtual void	ConstructProcess ()

void	SetAbsorption (G4bool)

void	SetNbOfPhotonsCerenkov (G4int)

void	SetVerbose (G4int)

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 47 of file WLSPhysicsList.hh.

Constructor & Destructor Documentation

WLSPhysicsList::WLSPhysicsList ( G4String physName )

Definition at line 69 of file WLSPhysicsList.cc.

References G4VUserPhysicsList::defaultCutValue, G4cout, G4endl, G4VModularPhysicsList::GetInstanceID(), G4VModularPhysicsList::GetPhysics(), G4VPhysicsConstructor::GetPhysicsName(), G4VModularPhysicsList::GetSubInstanceManager(), G4LossTableManager::Instance(), python.hepunit::mm, G4VUPLSplitter< T >::offset, and G4VModularPhysicsList::RegisterPhysics().

                                                 : G4VModularPhysicsList()
 {
     G4LossTableManager::Instance();
 
     defaultCutValue  = 1.*mm;
     fCutForGamma     = defaultCutValue;
     fCutForElectron  = defaultCutValue;
     fCutForPositron  = defaultCutValue;
 
 //    G4PhysListFactory factory;
     G4VModularPhysicsList* phys = NULL;
     if (physName == "QGSP_BERT_HP") {
        phys = new QGSP_BERT_HP;
     } else {
        phys = new FTFP_BERT;
     }
 //    if (factory.IsReferencePhysList(physName)) {
 //       phys = factory.GetReferencePhysList(physName);
 //       if(!phys)G4Exception("WLSPhysicsList::WLSPhysicsList","InvalidSetup",
 //                            FatalException,"PhysicsList does not exist");
        fMessenger = new WLSPhysicsListMessenger(this);
 //    }
 
     for (G4int i = 0; ; ++i) {
        G4VPhysicsConstructor* elem =
                   const_cast<G4VPhysicsConstructor*> (phys->GetPhysics(i));
        if (elem == NULL) break;
        G4cout << "RegisterPhysics: " << elem->GetPhysicsName() << G4endl;
        RegisterPhysics(elem);
     }
 
     fAbsorptionOn = true;
     
     //This looks complex, but it is not:
     //Get from base-class the pointer of the phsyicsVector
     //to be used. Remember: G4VModularPhysicsList is now a split class.
     //Why G4VModularPhysicsList::RegisterPhysics method is not used instead?
     //If possible we can remove this...
     fPhysicsVector =
                 GetSubInstanceManager().offset[GetInstanceID()].physicsVector;
     
     fPhysicsVector->push_back(new WLSExtraPhysics());
     fPhysicsVector->push_back(fOpticalPhysics =
                                         new WLSOpticalPhysics(fAbsorptionOn));
 
     fPhysicsVector->push_back(new G4RadioactiveDecayPhysics());
 
     fStepMaxProcess = new WLSStepMax();
 }

WLSPhysicsList::~WLSPhysicsList ( )

virtual

Definition at line 121 of file WLSPhysicsList.cc.

 {
     delete fMessenger;
 
     delete fStepMaxProcess;
 }

Member Function Documentation

void WLSPhysicsList::AddStepMax ( )

Definition at line 321 of file WLSPhysicsList.cc.

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

Referenced by ConstructProcess().

 {
   // Step limitation seen as a process
 
   theParticleIterator->reset();
   while ((*theParticleIterator)()){
       G4ParticleDefinition* particle = theParticleIterator->value();
       G4ProcessManager* pmanager = particle->GetProcessManager();
 
       if (fStepMaxProcess->IsApplicable(*particle) && !particle->IsShortLived())
       {
          if (pmanager) pmanager ->AddDiscreteProcess(fStepMaxProcess);
       }
   }
 }

void WLSPhysicsList::ClearPhysics ( )

Make sure that the physics list is empty.

Definition at line 130 of file WLSPhysicsList.cc.

Referenced by WLSPhysicsListMessenger::SetNewValue().

 {
     for (G4PhysConstVector::iterator p  = fPhysicsVector->begin();
                                      p != fPhysicsVector->end(); ++p) {
         delete (*p);
     }
     fPhysicsVector->clear();
 }

void WLSPhysicsList::ConstructParticle ( void )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 141 of file WLSPhysicsList.cc.

References G4VModularPhysicsList::ConstructParticle(), JA::Insert(), G4MuonMinus::MuonMinusDefinition(), and G4MuonPlus::MuonPlusDefinition().

 {
     G4VModularPhysicsList::ConstructParticle();
 
     G4DecayTable* MuonPlusDecayTable = new G4DecayTable();
     MuonPlusDecayTable -> Insert(new
                            G4MuonDecayChannelWithSpin("mu+",0.986));
     MuonPlusDecayTable -> Insert(new
                            G4MuonRadiativeDecayChannelWithSpin("mu+",0.014));
     G4MuonPlus::MuonPlusDefinition() -> SetDecayTable(MuonPlusDecayTable);
 
     G4DecayTable* MuonMinusDecayTable = new G4DecayTable();
     MuonMinusDecayTable -> Insert(new
                             G4MuonDecayChannelWithSpin("mu-",0.986));
     MuonMinusDecayTable -> Insert(new
                             G4MuonRadiativeDecayChannelWithSpin("mu-",0.014));
     G4MuonMinus::MuonMinusDefinition() -> SetDecayTable(MuonMinusDecayTable);
 
 }

void WLSPhysicsList::ConstructProcess ( void )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 163 of file WLSPhysicsList.cc.

References G4ProcessManager::AddProcess(), AddStepMax(), G4VModularPhysicsList::ConstructProcess(), G4INCL::ClusterDecay::decay(), G4ProcessTable::FindProcess(), G4ParticleDefinition::GetProcessManager(), G4ProcessTable::GetProcessTable(), idxAtRest, idxPostStep, G4MuonMinus::MuonMinus(), G4MuonPlus::MuonPlus(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), G4ProcessManager::RemoveProcess(), G4ProcessManager::SetProcessOrdering(), and SetVerbose().

 {
     G4VModularPhysicsList::ConstructProcess();
 
     SetVerbose(0);
 
     G4DecayWithSpin* decayWithSpin = new G4DecayWithSpin();
 
     G4ProcessTable* processTable = G4ProcessTable::GetProcessTable();
 
     G4VProcess* decay;
     decay = processTable->FindProcess("Decay",G4MuonPlus::MuonPlus());
 
     G4ProcessManager* pManager;
     pManager = G4MuonPlus::MuonPlus()->GetProcessManager();
 
     if (pManager) {
       if (decay) pManager->RemoveProcess(decay);
       pManager->AddProcess(decayWithSpin);
       // set ordering for PostStepDoIt and AtRestDoIt
       pManager ->SetProcessOrdering(decayWithSpin, idxPostStep);
       pManager ->SetProcessOrdering(decayWithSpin, idxAtRest);
     }
 
     decay = processTable->FindProcess("Decay",G4MuonMinus::MuonMinus());
 
     pManager = G4MuonMinus::MuonMinus()->GetProcessManager();
 
     if (pManager) {
       if (decay) pManager->RemoveProcess(decay);
       pManager->AddProcess(decayWithSpin);
       // set ordering for PostStepDoIt and AtRestDoIt
       pManager ->SetProcessOrdering(decayWithSpin, idxPostStep);
       pManager ->SetProcessOrdering(decayWithSpin, idxAtRest);
     }
 
     G4PionDecayMakeSpin* poldecay = new G4PionDecayMakeSpin();
 
     decay = processTable->FindProcess("Decay",G4PionPlus::PionPlus());
 
     pManager = G4PionPlus::PionPlus()->GetProcessManager();
 
     if (pManager) {
       if (decay) pManager->RemoveProcess(decay);
       pManager->AddProcess(poldecay);
       // set ordering for PostStepDoIt and AtRestDoIt
       pManager ->SetProcessOrdering(poldecay, idxPostStep);
       pManager ->SetProcessOrdering(poldecay, idxAtRest);
     }
 
     decay = processTable->FindProcess("Decay",G4PionMinus::PionMinus());
 
     pManager = G4PionMinus::PionMinus()->GetProcessManager();
 
     if (pManager) {
       if (decay) pManager->RemoveProcess(decay);
       pManager->AddProcess(poldecay);
       // set ordering for PostStepDoIt and AtRestDoIt
       pManager ->SetProcessOrdering(poldecay, idxPostStep);
       pManager ->SetProcessOrdering(poldecay, idxAtRest);
     }
 
     AddStepMax();
 
 }

WLSStepMax * WLSPhysicsList::GetStepMaxProcess ( )

Definition at line 314 of file WLSPhysicsList.cc.

 {
   return fStepMaxProcess;
 }

void WLSPhysicsList::RemoveFromPhysicsList ( const G4String & name )

Remove specific physics from physics list.

Definition at line 231 of file WLSPhysicsList.cc.

References FatalException, G4Exception(), and G4VPhysicsConstructor::GetPhysicsName().

Referenced by SetAbsorption(), and WLSPhysicsListMessenger::SetNewValue().

 {
     G4bool success = false;
     for (G4PhysConstVector::iterator p  = fPhysicsVector->begin();
                                      p != fPhysicsVector->end(); ++p) {
         G4VPhysicsConstructor* e = (*p);
         if (e->GetPhysicsName() == name) {
            fPhysicsVector->erase(p);
            success = true;
            break;
         }
     }
     if (!success) {
        G4ExceptionDescription message;
        message << "PhysicsList::RemoveFromEMPhysicsList "<< name << "not found";
        G4Exception("example WLSPhysicsList::RemoveFromPhysicsList()",
        "ExamWLSPhysicsList01",FatalException,message);
     }
 }

void WLSPhysicsList::SetAbsorption ( G4bool toggle )

Definition at line 253 of file WLSPhysicsList.cc.

References WLSOpticalPhysics::ConstructProcess(), and RemoveFromPhysicsList().

Referenced by WLSPhysicsListMessenger::SetNewValue().

 {
        fAbsorptionOn = toggle;
        RemoveFromPhysicsList("Optical");
        fPhysicsVector->
                     push_back(fOpticalPhysics = new WLSOpticalPhysics(toggle));
        fOpticalPhysics->ConstructProcess();
 }

void WLSPhysicsList::SetCutForElectron ( G4double cut )

Definition at line 291 of file WLSPhysicsList.cc.

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

Referenced by WLSPhysicsListMessenger::SetNewValue().

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

void WLSPhysicsList::SetCutForGamma ( G4double cut )

Definition at line 283 of file WLSPhysicsList.cc.

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

Referenced by WLSPhysicsListMessenger::SetNewValue().

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

void WLSPhysicsList::SetCutForPositron ( G4double cut )

Definition at line 299 of file WLSPhysicsList.cc.

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

Referenced by WLSPhysicsListMessenger::SetNewValue().

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

void WLSPhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 264 of file WLSPhysicsList.cc.

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

 {
     if (verboseLevel >0) {
         G4cout << "WLSPhysicsList::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(fCutForGamma, "gamma");
     SetCutValue(fCutForElectron, "e-");
     SetCutValue(fCutForPositron, "e+");
 
     if (verboseLevel>0) DumpCutValuesTable();
 }

void WLSPhysicsList::SetNbOfPhotonsCerenkov ( G4int maxNumber )

Definition at line 339 of file WLSPhysicsList.cc.

References WLSOpticalPhysics::SetNbOfPhotonsCerenkov().

 {
    fOpticalPhysics->SetNbOfPhotonsCerenkov(maxNumber);
 }

void WLSPhysicsList::SetStepMax ( G4double step )

Definition at line 307 of file WLSPhysicsList.cc.

References WLSStepMax::SetStepMax().

Referenced by WLSPhysicsListMessenger::SetNewValue().

 {
   fStepMaxProcess->SetStepMax(step);
 }

void WLSPhysicsList::SetVerbose ( G4int verbose )

Definition at line 346 of file WLSPhysicsList.cc.

References WLSOpticalPhysics::GetAbsorptionProcess(), WLSOpticalPhysics::GetBoundaryProcess(), WLSOpticalPhysics::GetCerenkovProcess(), WLSOpticalPhysics::GetMieHGScatteringProcess(), WLSOpticalPhysics::GetRayleighScatteringProcess(), WLSOpticalPhysics::GetScintillationProcess(), and G4VProcess::SetVerboseLevel().

Referenced by ConstructProcess(), and WLSPhysicsListMessenger::SetNewValue().

 {
    fOpticalPhysics->GetCerenkovProcess()->SetVerboseLevel(verbose);
    fOpticalPhysics->GetScintillationProcess()->SetVerboseLevel(verbose);
    fOpticalPhysics->GetAbsorptionProcess()->SetVerboseLevel(verbose);
    fOpticalPhysics->GetRayleighScatteringProcess()->SetVerboseLevel(verbose);
    fOpticalPhysics->GetMieHGScatteringProcess()->SetVerboseLevel(verbose);
    fOpticalPhysics->GetBoundaryProcess()->SetVerboseLevel(verbose);
 }

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation