F04PhysicsList Class Reference

#include <Doxymodules_field.h>

Inheritance diagram for F04PhysicsList:

Public Member Functions
	F04PhysicsList (G4String)
	~F04PhysicsList () override
void	SetStepMax (G4double)
F04StepMax *	GetStepMaxProcess ()
void	AddStepMax ()
void	ConstructParticle () override
void	ConstructProcess () override

Private Attributes
G4double	fMaxChargedStep = DBL_MAX
F04PhysicsListMessenger *	fMessenger = nullptr

Static Private Attributes
static G4ThreadLocal F04StepMax *	fStepMaxProcess = nullptr

Detailed Description

Definition at line 110 of file Doxymodules_field.h.

Constructor & Destructor Documentation

F04PhysicsList()

F04PhysicsList::F04PhysicsList

(

G4String

physName

)

Definition at line 72 of file F04PhysicsList.cc.

                                                : G4VModularPhysicsList()
{
    G4LossTableManager::Instance();
 
    defaultCutValue  = 1.*mm;
 
    fMessenger = new F04PhysicsListMessenger(this);
 
    SetVerboseLevel(1);
 
//    G4PhysListFactory factory;
    G4VModularPhysicsList* phys = nullptr;
    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("F04PhysicsList::F04PhysicsList","InvalidSetup",
                              FatalException,"PhysicsList does not exist");
 
    for (G4int i = 0; ; ++i) {
       auto  elem =
                  const_cast<G4VPhysicsConstructor*> (phys->GetPhysics(i));
       if (elem == nullptr) break;
       G4cout << "RegisterPhysics: " << elem->GetPhysicsName() << G4endl;
       RegisterPhysics(elem);
    }
 
    RegisterPhysics(new G4StepLimiterPhysics());
    RegisterPhysics(new G4OpticalPhysics());
}

~F04PhysicsList()

F04PhysicsList::~F04PhysicsList ( )

override

Definition at line 113 of file F04PhysicsList.cc.

{
    delete fMessenger;
 
    //delete fStepMaxProcess;
}

Member Function Documentation

SetStepMax()

void F04PhysicsList::SetStepMax

(

G4double

step

)

Definition at line 256 of file F04PhysicsList.cc.

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

GetStepMaxProcess()

F04StepMax * F04PhysicsList::GetStepMaxProcess

(

)

Definition at line 264 of file F04PhysicsList.cc.

{
  return fStepMaxProcess;
}

AddStepMax()

void F04PhysicsList::AddStepMax

(

)

Definition at line 271 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);
      }
  }
}

ConstructParticle()

void F04PhysicsList::ConstructParticle ( )

override

Definition at line 122 of file F04PhysicsList.cc.

{
    G4VModularPhysicsList::ConstructParticle();
 
    G4GenericIon::GenericIonDefinition();
 
    auto  muonPlusDecayTable = new G4DecayTable();
    muonPlusDecayTable -> Insert(new
                           G4MuonDecayChannelWithSpin("mu+",0.986));
    muonPlusDecayTable -> Insert(new
                           G4MuonRadiativeDecayChannelWithSpin("mu+",0.014));
    G4MuonPlus::MuonPlusDefinition() -> SetDecayTable(muonPlusDecayTable);
 
    auto  muonMinusDecayTable = new G4DecayTable();
    muonMinusDecayTable -> Insert(new
                            G4MuonDecayChannelWithSpin("mu-",0.986));
    muonMinusDecayTable -> Insert(new
                            G4MuonRadiativeDecayChannelWithSpin("mu-",0.014));
    G4MuonMinus::MuonMinusDefinition() -> SetDecayTable(muonMinusDecayTable);
}

ConstructProcess()

void F04PhysicsList::ConstructProcess ( )

override

Definition at line 145 of file F04PhysicsList.cc.

{
    G4VModularPhysicsList::ConstructProcess();
 
    fStepMaxProcess = new F04StepMax();
    G4AutoDelete::Register(fStepMaxProcess);
 
    auto  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);
    }
 
    G4VProcess* process = processTable->
         FindProcess("muMinusCaptureAtRest",G4MuonMinus::MuonMinus());
 
    if (pmanager) {
       if (process) pmanager->RemoveProcess(process);
       process = new G4MuonMinusCapture(new G4MuMinusCapturePrecompound());
       pmanager->AddRestProcess(process);
    }
 
    auto  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();
}

Member Data Documentation

fMaxChargedStep

G4double F04PhysicsList::fMaxChargedStep = DBL_MAX

private

Definition at line 64 of file F04PhysicsList.hh.

fStepMaxProcess

G4ThreadLocal F04StepMax * F04PhysicsList::fStepMaxProcess = nullptr

staticprivate

Definition at line 65 of file F04PhysicsList.hh.

fMessenger

F04PhysicsListMessenger* F04PhysicsList::fMessenger = nullptr

private

Definition at line 67 of file F04PhysicsList.hh.

---

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

• .doxygen/Doxymodules_field.h
• extended/field/field04/include/F04PhysicsList.hh
• extended/field/field04/src/F04PhysicsList.cc