F04PhysicsList Class Reference

#include <F04PhysicsList.hh>

Inheritance diagram for F04PhysicsList:

Collaboration diagram for F04PhysicsList:

Public Member Functions
	F04PhysicsList (G4String)
virtual	~F04PhysicsList ()
virtual void	SetCuts ()
void	SetCutForGamma (G4double)
void	SetCutForElectron (G4double)
void	SetCutForPositron (G4double)
void	SetStepMax (G4double)
F04StepMax *	GetStepMaxProcess ()
void	AddStepMax ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
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
G4double	fMaxChargedStep
F04PhysicsListMessenger *	fMessenger

Static Private Attributes
static G4ThreadLocal F04StepMax *	fStepMaxProcess = 0

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 43 of file F04PhysicsList.hh.

Constructor & Destructor Documentation

F04PhysicsList()

F04PhysicsList::F04PhysicsList

(

G4String

physName

)

Definition at line 70 of file F04PhysicsList.cc.

                                                : G4VModularPhysicsList()
{
    G4LossTableManager::Instance();

    defaultCutValue  = 1.*mm;
    fCutForGamma     = defaultCutValue;
    fCutForElectron  = defaultCutValue;
    fCutForPositron  = defaultCutValue;

    fMessenger = new F04PhysicsListMessenger(this);

    SetVerboseLevel(1);

//    G4PhysListFactory factory;
    G4VModularPhysicsList* phys = NULL;
    if (physName == "QGSP_BERT") {
       phys = new QGSP_BERT;
    } else {
       phys = new FTFP_BERT;
    }

//    if (factory.IsReferencePhysList(physName))
//       phys =factory.GetReferencePhysList(physName);

    // Physics List is defined via environment variable PHYSLIST
//    if (!phys) phys = factory.ReferencePhysList();

    if (!phys) G4Exception("WLSPhysicsList::WLSPhysicsList","InvalidSetup",
                              FatalException,"PhysicsList does not exist");

    for (G4int i = 0; ; ++i) {
       G4VPhysicsConstructor* elem =
                  const_cast<G4VPhysicsConstructor*> (phys->GetPhysics(i));
       if (elem == NULL) break;
       G4cout << "RegisterPhysics: " << elem->GetPhysicsName() << G4endl;
       RegisterPhysics(elem);
    }

    RegisterPhysics(new F04ExtraPhysics());
    RegisterPhysics(new G4OpticalPhysics());

    fMaxChargedStep = DBL_MAX;
}

Here is the call graph for this function:

~F04PhysicsList()

F04PhysicsList::~F04PhysicsList ( )

virtual

Definition at line 116 of file F04PhysicsList.cc.

{
    delete fMessenger;

    //delete fStepMaxProcess;
}

Member Function Documentation

AddStepMax()

void F04PhysicsList::AddStepMax

(

)

Definition at line 308 of file F04PhysicsList.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:

ConstructParticle()

void F04PhysicsList::ConstructParticle ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 125 of file F04PhysicsList.cc.

{
    G4VModularPhysicsList::ConstructParticle();

    G4GenericIon::GenericIonDefinition();

    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 F04PhysicsList::ConstructProcess ( void  )

virtual

Reimplemented from G4VModularPhysicsList.

Definition at line 148 of file F04PhysicsList.cc.

{
    G4VModularPhysicsList::ConstructProcess();

    fStepMaxProcess = new F04StepMax();
    G4AutoDelete::Register(fStepMaxProcess);
    
    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()

F04StepMax * F04PhysicsList::GetStepMaxProcess

(

)

Definition at line 301 of file F04PhysicsList.cc.

{
  return fStepMaxProcess;
}

SetCutForElectron()

void F04PhysicsList::SetCutForElectron

(

G4double

cut

)

Definition at line 277 of file F04PhysicsList.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 F04PhysicsList::SetCutForGamma

(

G4double

cut

)

Definition at line 269 of file F04PhysicsList.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 F04PhysicsList::SetCutForPositron

(

G4double

cut

)

Definition at line 285 of file F04PhysicsList.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 F04PhysicsList::SetCuts ( )

virtual

Reimplemented from G4VUserPhysicsList.

Definition at line 250 of file F04PhysicsList.cc.

{
    if (verboseLevel >0) {
        G4cout << "F04PhysicsList::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:

SetStepMax()

void F04PhysicsList::SetStepMax

(

G4double

step

)

Definition at line 293 of file F04PhysicsList.cc.

{
  fMaxChargedStep = step ;
  fStepMaxProcess->SetStepMax(fMaxChargedStep);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

fCutForElectron

G4double F04PhysicsList::fCutForElectron

private

Definition at line 71 of file F04PhysicsList.hh.

fCutForGamma

G4double F04PhysicsList::fCutForGamma

private

Definition at line 70 of file F04PhysicsList.hh.

fCutForPositron

G4double F04PhysicsList::fCutForPositron

private

Definition at line 72 of file F04PhysicsList.hh.

fMaxChargedStep

G4double F04PhysicsList::fMaxChargedStep

private

Definition at line 74 of file F04PhysicsList.hh.

fMessenger

F04PhysicsListMessenger* F04PhysicsList::fMessenger

private

Definition at line 77 of file F04PhysicsList.hh.

fStepMaxProcess

G4ThreadLocal F04StepMax * F04PhysicsList::fStepMaxProcess = 0

staticprivate

Definition at line 75 of file F04PhysicsList.hh.

---

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

• F04PhysicsList.hh
• F04PhysicsList.cc