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

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

(

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;
    
    RegisterPhysics(new WLSExtraPhysics());
    RegisterPhysics(fOpticalPhysics = new WLSOpticalPhysics(fAbsorptionOn));

    RegisterPhysics(new G4RadioactiveDecayPhysics());

    fStepMaxProcess = new WLSStepMax();
}

Here is the call graph for this function:

WLSPhysicsList::~WLSPhysicsList ( )

virtual

Definition at line 112 of file WLSPhysicsList.cc.

{
    delete fMessenger;

    delete fStepMaxProcess;
}

Member Function Documentation

void WLSPhysicsList::AddStepMax

(

)

Definition at line 312 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:

void WLSPhysicsList::ClearPhysics

(

)

Make sure that the physics list is empty.

Definition at line 121 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:

void WLSPhysicsList::ConstructParticle ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 132 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:

void WLSPhysicsList::ConstructProcess ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 154 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:

WLSStepMax * WLSPhysicsList::GetStepMaxProcess

(

)

Definition at line 305 of file WLSPhysicsList.cc.

{
  return fStepMaxProcess;
}

void WLSPhysicsList::RemoveFromPhysicsList

(

const G4String &

name

)

Remove specific physics from physics list.

Definition at line 222 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:

void WLSPhysicsList::SetAbsorption

(

G4bool

toggle

)

Definition at line 244 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:

void WLSPhysicsList::SetCutForElectron

(

G4double

cut

)

Definition at line 282 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:

void WLSPhysicsList::SetCutForGamma

(

G4double

cut

)

Definition at line 274 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:

void WLSPhysicsList::SetCutForPositron

(

G4double

cut

)

Definition at line 290 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:

void WLSPhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 255 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:

void WLSPhysicsList::SetNbOfPhotonsCerenkov

(

G4int

maxNumber

)

Definition at line 331 of file WLSPhysicsList.cc.

{
   fOpticalPhysics->SetNbOfPhotonsCerenkov(maxNumber);
}

Here is the call graph for this function:

void WLSPhysicsList::SetStepMax

(

G4double

step

)

Definition at line 298 of file WLSPhysicsList.cc.

{
  fStepMaxProcess->SetStepMax(step);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void WLSPhysicsList::SetVerbose

(

G4int

verbose

)

Definition at line 338 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:

---

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

• source/geant4.10.03.p03/examples/extended/optical/wls/include/WLSPhysicsList.hh
• source/geant4.10.03.p03/examples/extended/optical/wls/src/WLSPhysicsList.cc