WLSPhysicsList Class Reference

#include <WLSPhysicsList.hh>

Inheritance diagram for WLSPhysicsList:

Collaboration 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 ()

Private Attributes
G4double	fCutForGamma
G4double	fCutForElectron
G4double	fCutForPositron
WLSStepMax *	fStepMaxProcess
WLSOpticalPhysics *	fOpticalPhysics
WLSPhysicsListMessenger *	fMessenger
G4bool	fAbsorptionOn
G4VMPLData::G4PhysConstVectorData *	fPhysicsVector

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 ()
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
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::WLSPhysicsList

(

G4String

physName

)

Definition at line 69 of file WLSPhysicsList.cc.

                                                : 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();
}

Here is the call graph for this function:

~WLSPhysicsList()

WLSPhysicsList::~WLSPhysicsList ( )

virtual

Definition at line 121 of file WLSPhysicsList.cc.

{
    delete fMessenger;

    delete fStepMaxProcess;
}

Member Function Documentation

AddStepMax()

void WLSPhysicsList::AddStepMax

(

)

Definition at line 321 of file WLSPhysicsList.cc.

{
  // Step limitation seen as a process

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

      if (fStepMaxProcess->IsApplicable(*particle) && !particle->IsShortLived())
      {
         if (pmanager) pmanager ->AddDiscreteProcess(fStepMaxProcess);
      }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ClearPhysics()

void WLSPhysicsList::ClearPhysics

(

)

Make sure that the physics list is empty.

Definition at line 130 of file WLSPhysicsList.cc.

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

Here is the caller graph for this function:

ConstructParticle()

void WLSPhysicsList::ConstructParticle ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 141 of file WLSPhysicsList.cc.

{
    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);

}

Here is the call graph for this function:

ConstructProcess()

void WLSPhysicsList::ConstructProcess ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 163 of file WLSPhysicsList.cc.

{
    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();

}

Here is the call graph for this function:

GetStepMaxProcess()

WLSStepMax * WLSPhysicsList::GetStepMaxProcess

(

)

Definition at line 314 of file WLSPhysicsList.cc.

{
  return fStepMaxProcess;
}

RemoveFromPhysicsList()

void WLSPhysicsList::RemoveFromPhysicsList

(

const G4String &

name

)

Remove specific physics from physics list.

Definition at line 231 of file WLSPhysicsList.cc.

{
    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);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetAbsorption()

void WLSPhysicsList::SetAbsorption

(

G4bool

toggle

)

Definition at line 253 of file WLSPhysicsList.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

SetCutForElectron()

void WLSPhysicsList::SetCutForElectron

(

G4double

cut

)

Definition at line 291 of file WLSPhysicsList.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

SetCutForGamma()

void WLSPhysicsList::SetCutForGamma

(

G4double

cut

)

Definition at line 283 of file WLSPhysicsList.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

SetCutForPositron()

void WLSPhysicsList::SetCutForPositron

(

G4double

cut

)

Definition at line 299 of file WLSPhysicsList.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

SetCuts()

void WLSPhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 264 of file WLSPhysicsList.cc.

{
    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();
}

Here is the call graph for this function:

SetNbOfPhotonsCerenkov()

void WLSPhysicsList::SetNbOfPhotonsCerenkov

(

G4int

maxNumber

)

Definition at line 340 of file WLSPhysicsList.cc.

{
   fOpticalPhysics->SetNbOfPhotonsCerenkov(maxNumber);
}

Here is the call graph for this function:

SetStepMax()

void WLSPhysicsList::SetStepMax

(

G4double

step

)

Definition at line 307 of file WLSPhysicsList.cc.

{
  fStepMaxProcess->SetStepMax(step);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetVerbose()

void WLSPhysicsList::SetVerbose

(

G4int

verbose

)

Definition at line 347 of file WLSPhysicsList.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

fAbsorptionOn

G4bool WLSPhysicsList::fAbsorptionOn

private

Definition at line 91 of file WLSPhysicsList.hh.

fCutForElectron

G4double WLSPhysicsList::fCutForElectron

private

Definition at line 82 of file WLSPhysicsList.hh.

fCutForGamma

G4double WLSPhysicsList::fCutForGamma

private

Definition at line 81 of file WLSPhysicsList.hh.

fCutForPositron

G4double WLSPhysicsList::fCutForPositron

private

Definition at line 83 of file WLSPhysicsList.hh.

fMessenger

WLSPhysicsListMessenger* WLSPhysicsList::fMessenger

private

Definition at line 89 of file WLSPhysicsList.hh.

fOpticalPhysics

WLSOpticalPhysics* WLSPhysicsList::fOpticalPhysics

private

Definition at line 87 of file WLSPhysicsList.hh.

fPhysicsVector

G4VMPLData::G4PhysConstVectorData* WLSPhysicsList::fPhysicsVector

private

Definition at line 93 of file WLSPhysicsList.hh.

fStepMaxProcess

WLSStepMax* WLSPhysicsList::fStepMaxProcess

private

Definition at line 85 of file WLSPhysicsList.hh.

---

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

• WLSPhysicsList.hh
• WLSPhysicsList.cc