PhysicsList Class Reference

#include <PhysicsList.hh>

Inheritance diagram for PhysicsList:

Public Member Functions
	PhysicsList (DetectorConstruction *)
	~PhysicsList () override
void	ConstructParticle () override
void	ConstructProcess () override
void	AddPhysicsList (const G4String &name)
void	AddDecay ()
void	AddRadioactiveDecay ()
void	AddStepMax ()
void	GetRange (G4double)
G4VProcess *	GetProcess (const G4String &) const

Private Attributes
G4VPhysicsConstructor *	fEmPhysicsList = nullptr
G4String	fEmName = " "
DetectorConstruction *	fDet = nullptr
PhysicsListMessenger *	fMessenger = nullptr

Detailed Description

Definition at line 50 of file PhysicsList.hh.

Constructor & Destructor Documentation

PhysicsList()

PhysicsList::PhysicsList

(

DetectorConstruction *

det

)

Definition at line 84 of file PhysicsList.cc.

  : fDet(det)
{
  fMessenger = new PhysicsListMessenger(this);
  SetVerboseLevel(1);
 
  // EM physics
  AddPhysicsList("emstandard_opt3");    
  
  // fix lower limit for cut
  G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(10*eV, 1*GeV);
  SetDefaultCutValue(1*mm);
}

~PhysicsList()

PhysicsList::~PhysicsList ( )

override

Definition at line 100 of file PhysicsList.cc.

{
  delete fMessenger;
}

Member Function Documentation

ConstructParticle()

void PhysicsList::ConstructParticle ( )

override

Definition at line 107 of file PhysicsList.cc.

{
    G4BosonConstructor  pBosonConstructor;
    pBosonConstructor.ConstructParticle();
 
    G4LeptonConstructor pLeptonConstructor;
    pLeptonConstructor.ConstructParticle();
 
    G4MesonConstructor pMesonConstructor;
    pMesonConstructor.ConstructParticle();
 
    G4BaryonConstructor pBaryonConstructor;
    pBaryonConstructor.ConstructParticle();
 
    G4IonConstructor pIonConstructor;
    pIonConstructor.ConstructParticle();
 
    G4ShortLivedConstructor pShortLivedConstructor;
    pShortLivedConstructor.ConstructParticle();  
}

ConstructProcess()

void PhysicsList::ConstructProcess ( )

override

Definition at line 130 of file PhysicsList.cc.

{
  // Transportation
  //
  AddTransportation();
 
  // Electromagnetic physics list
  //
  fEmPhysicsList->ConstructProcess();
    
  // Decay Process
  //
  AddDecay();
    
  // Decay Process
  //
  AddRadioactiveDecay();  
 
  // step limitation (as a full process)
  //  
  AddStepMax();
  
  // example of Get process
  auto process = GetProcess("RadioactiveDecay");
  if (process != nullptr) {
    G4cout << "\n  GetProcess : " << process->GetProcessName() << G4endl;
  }
}

AddPhysicsList()

void PhysicsList::AddPhysicsList

(

const G4String &

name

)

Definition at line 161 of file PhysicsList.cc.

{
  if (verboseLevel>0) {
    G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
  }
  
  if (name == fEmName) return;
 
  if (name == "local") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new PhysListEmStandard(name);
    
  } else if (name == "emstandard_opt0") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysics();
 
  } else if (name == "emstandard_opt1") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysics_option1();
 
  } else if (name == "emstandard_opt2") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysics_option2();
 
  } else if (name == "emstandard_opt3") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysics_option3();
    
  } else if (name == "emstandard_opt4") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysics_option4();
        
  } else if (name == "emstandardSS") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysicsSS();
    
  } else if (name == "emstandardGS") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysicsGS();
    
  } else if (name == "emstandardWVI") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmStandardPhysicsWVI();
    
  } else if (name == "emlivermore") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmLivermorePhysics();
    
  } else if (name == "empenelope") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmPenelopePhysics();
            
  } else if (name == "emlowenergy") {
    fEmName = name;
    delete fEmPhysicsList;
    fEmPhysicsList = new G4EmLowEPPhysics();
            
  } else {
 
    G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
           << " is not defined"
           << G4endl;
  }
  
  // Em options
  //
  G4EmParameters::Instance()->SetBuildCSDARange(true);        
  G4EmParameters::Instance()->SetGeneralProcessActive(false);
}

AddDecay()

void PhysicsList::AddDecay

(

)

Definition at line 244 of file PhysicsList.cc.

{
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
    
  // Decay Process
  //
  G4Decay* fDecayProcess = new G4Decay();
 
  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while( (*particleIterator)() ){
    G4ParticleDefinition* particle = particleIterator->value();
    if (fDecayProcess->IsApplicable(*particle) && !particle->IsShortLived()) 
      ph->RegisterProcess(fDecayProcess, particle);    
  }
}

AddRadioactiveDecay()

void PhysicsList::AddRadioactiveDecay

(

)

Definition at line 263 of file PhysicsList.cc.

{  
  G4RadioactiveDecay* radioactiveDecay = new G4RadioactiveDecay();
  
  G4bool armFlag = false;
  radioactiveDecay->SetARM(armFlag);                //Atomic Rearangement
 
  // atomic de-excitation module
  if (armFlag) {
    G4EmParameters::Instance()->SetAuger(true);
    G4EmParameters::Instance()->SetDeexcitationIgnoreCut(true);
  }
  
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();  
  ph->RegisterProcess(radioactiveDecay, G4GenericIon::GenericIon());
  
  // mandatory for G4NuclideTable
  //
  const G4double meanLife = 1*picosecond, halfLife = meanLife*std::log(2);
  G4NuclideTable::GetInstance()->SetThresholdOfHalfLife(halfLife);
}

AddStepMax()

void PhysicsList::AddStepMax

(

)

Definition at line 287 of file PhysicsList.cc.

{
  // Step limitation seen as a process
  StepMax* stepMaxProcess = new StepMax();
 
  auto particleIterator=GetParticleIterator();
  particleIterator->reset();
  while ((*particleIterator)()){
    G4ParticleDefinition* particle = particleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
 
    if (stepMaxProcess->IsApplicable(*particle) && !particle->IsShortLived())
      pmanager->AddDiscreteProcess(stepMaxProcess);
  }
}

GetRange()

void PhysicsList::GetRange

(

G4double

val

)

Definition at line 305 of file PhysicsList.cc.

{
  G4LogicalVolume* lBox = fDet->GetWorld()->GetLogicalVolume();
  const G4MaterialCutsCouple* couple = lBox->GetMaterialCutsCouple();
  const G4Material* currMat = lBox->GetMaterial();
 
  G4ParticleDefinition* part;
  G4double cut;
  part = G4Electron::Electron();
  cut = G4LossTableManager::Instance()->GetRange(part,val,couple);
  G4cout << "material : " << currMat->GetName()       << G4endl;
  G4cout << "particle : " << part->GetParticleName()  << G4endl;
  G4cout << "energy   : " << G4BestUnit(val,"Energy") << G4endl;
  G4cout << "range    : " << G4BestUnit(cut,"Length") << G4endl;
}

GetProcess()

G4VProcess * PhysicsList::GetProcess

(

const G4String &

processName

)

const

Definition at line 323 of file PhysicsList.cc.

{
  G4ParticleDefinition* particle = G4GenericIon::GenericIon();
  G4ProcessVector* procList = particle->GetProcessManager()->GetProcessList();
  G4int nbProc = particle->GetProcessManager()->GetProcessListLength();
  for (G4int k=0; k<nbProc; k++) {
    G4VProcess* process = (*procList)[k];
    if (process->GetProcessName() == processName) return process;
  }
  return nullptr;
}

Member Data Documentation

fEmPhysicsList

G4VPhysicsConstructor* PhysicsList::fEmPhysicsList = nullptr

private

Definition at line 68 of file PhysicsList.hh.

fEmName

G4String PhysicsList::fEmName = " "

private

Definition at line 69 of file PhysicsList.hh.

fDet

DetectorConstruction* PhysicsList::fDet = nullptr

private

Definition at line 70 of file PhysicsList.hh.

fMessenger

PhysicsListMessenger* PhysicsList::fMessenger = nullptr

private

Definition at line 71 of file PhysicsList.hh.

---

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

• extended/electromagnetic/TestEm1/include/PhysicsList.hh
• extended/electromagnetic/TestEm1/src/PhysicsList.cc