G4EmDNAPhysics_option7 Class Reference

#include <G4EmDNAPhysics_option7.hh>

Inheritance diagram for G4EmDNAPhysics_option7:

Public Member Functions
	G4EmDNAPhysics_option7 (G4int ver=1, const G4String &name="")
virtual	~G4EmDNAPhysics_option7 ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")
	G4VPhysicsConstructor (const G4String &name, G4int physics_type)
virtual	~G4VPhysicsConstructor ()
virtual void	ConstructParticle ()=0
virtual void	ConstructProcess ()=0
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const
virtual void	TerminateWorker ()

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Types inherited from G4VPhysicsConstructor
using	PhysicsBuilder_V = G4VPCData::PhysicsBuilders_V
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
PhysicsBuilder_V	GetBuilders () const
void	AddBuilder (G4PhysicsBuilderInterface *bld)
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel
G4String	namePhysics
G4int	typePhysics
G4ParticleTable *	theParticleTable
G4int	g4vpcInstanceID
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 37 of file G4EmDNAPhysics_option7.hh.

Constructor & Destructor Documentation

G4EmDNAPhysics_option7()

G4EmDNAPhysics_option7::G4EmDNAPhysics_option7 ( G4int  ver = 1, const G4String &  name = ""  )

explicit

Definition at line 94 of file G4EmDNAPhysics_option7.cc.

                                                                         :
    G4VPhysicsConstructor("G4EmDNAPhysics_option7"), verbose(ver)
{
  G4EmParameters* param = G4EmParameters::Instance();
  param->SetDefaults();
  param->SetFluo(true);  
  param->SetAuger(true);  
  param->SetAugerCascade(true);  
  param->SetDeexcitationIgnoreCut(true);
  param->ActivateDNA();
 
  SetPhysicsType(bElectromagnetic);
}

~G4EmDNAPhysics_option7()

G4EmDNAPhysics_option7::~G4EmDNAPhysics_option7 ( )

virtual

Definition at line 110 of file G4EmDNAPhysics_option7.cc.

{
}

Member Function Documentation

ConstructParticle()

void G4EmDNAPhysics_option7::ConstructParticle ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 116 of file G4EmDNAPhysics_option7.cc.

{
// bosons
  G4Gamma::Gamma();
 
// leptons
  G4Electron::Electron();
  G4Positron::Positron();
 
// baryons
  G4Proton::Proton();
 
  G4GenericIon::GenericIonDefinition();
 
  G4DNAGenericIonsManager * genericIonsManager;
  genericIonsManager = G4DNAGenericIonsManager::Instance();
  genericIonsManager->GetIon("alpha++");
  genericIonsManager->GetIon("alpha+");
  genericIonsManager->GetIon("helium");
  genericIonsManager->GetIon("hydrogen");
 
}

ConstructProcess()

void G4EmDNAPhysics_option7::ConstructProcess ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 141 of file G4EmDNAPhysics_option7.cc.

{
  if(verbose > 1)
  {
    G4cout << "### " << GetPhysicsName() << " Construct Processes " << G4endl;
  }
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
  auto myParticleIterator=GetParticleIterator();
  myParticleIterator->reset();
  while( (*myParticleIterator)() )
  {
    G4ParticleDefinition* particle = myParticleIterator->value();
    G4String particleName = particle->GetParticleName();
 
    if (particleName == "e-")
    {
      //  *** Electron solvation ***
      // Either kills the electron track or transform the electron to
      // a solvated electron, depending on whether the chemistry is
      // activated
      
      G4DNAElectronSolvation* solvation =
        new G4DNAElectronSolvation("e-_G4DNAElectronSolvation");
      auto therm = G4DNASolvationModelFactory::GetMacroDefinedModel();
      therm->SetHighEnergyLimit(10.*eV); // limit of the Uehara's model
      solvation->SetEmModel(therm);
      ph->RegisterProcess(solvation, particle);
      
      //  *** Elastic scattering ***
      
      G4DNAElastic* theDNAElasticProcess =
        new G4DNAElastic("e-_G4DNAElastic");
      G4DNAUeharaScreenedRutherfordElasticModel* emElast =
        new G4DNAUeharaScreenedRutherfordElasticModel();
      emElast->SetHighEnergyLimit(1*MeV);
      theDNAElasticProcess->SetEmModel(emElast);
      ph->RegisterProcess(theDNAElasticProcess, particle);
 
      //  *** Excitation ***
       
      G4DNAExcitation* theDNAExcitationProcess =
          new G4DNAExcitation("e-_G4DNAExcitation");
 
      // 1-st model should be Set
      G4DNAEmfietzoglouExcitationModel* emExc = new G4DNAEmfietzoglouExcitationModel();
      emExc->SetActivationLowEnergyLimit(8*eV);
      emExc->SetActivationHighEnergyLimit(10.*keV);
      theDNAExcitationProcess->SetEmModel(emExc);
 
      // 2-nd model should be Add
      G4DNABornExcitationModel* bornExc = new G4DNABornExcitationModel();
      bornExc->SetActivationLowEnergyLimit(10*keV);
      bornExc->SetActivationHighEnergyLimit(1.*MeV);
      theDNAExcitationProcess->AddEmModel(2, bornExc);
      
      ph->RegisterProcess(theDNAExcitationProcess, particle);
 
      //  *** Ionisation ***
       
      G4DNAIonisation* theDNAIonisationProcess =
          new G4DNAIonisation("e-_G4DNAIonisation");
 
      // 1-st model should be Set
      G4DNAEmfietzoglouIonisationModel* emIonModel = new G4DNAEmfietzoglouIonisationModel();
      emIonModel->SetActivationLowEnergyLimit(10.*eV);
      emIonModel->SetActivationHighEnergyLimit(10.*keV);
      theDNAIonisationProcess->SetEmModel(emIonModel);
 
      // 2-nd model should be Add
      G4DNABornIonisationModel* bornIon = new G4DNABornIonisationModel();
      bornIon->SetActivationLowEnergyLimit(10*keV);
      bornIon->SetActivationHighEnergyLimit(1.*MeV);
      theDNAIonisationProcess->AddEmModel(2,bornIon);
      
      ph->RegisterProcess(theDNAIonisationProcess, particle);
      
      //  *** Vibrational excitation ***
      ph->RegisterProcess(new G4DNAVibExcitation("e-_G4DNAVibExcitation"), particle);
      
      //  *** Attachment ***
      ph->RegisterProcess(new G4DNAAttachment("e-_G4DNAAttachment"), particle);
    
    } else if ( particleName == "proton" ) {
      ph->RegisterProcess(new G4DNAElastic("proton_G4DNAElastic"), particle);
      ph->RegisterProcess(new G4DNAExcitation("proton_G4DNAExcitation"), particle);
      ph->RegisterProcess(new G4DNAIonisation("proton_G4DNAIonisation"), particle);
      ph->RegisterProcess(new G4DNAChargeDecrease("proton_G4DNAChargeDecrease"), particle);
 
    } else if ( particleName == "hydrogen" ) {
      ph->RegisterProcess(new G4DNAElastic("hydrogen_G4DNAElastic"), particle);
      ph->RegisterProcess(new G4DNAExcitation("hydrogen_G4DNAExcitation"), particle);
      ph->RegisterProcess(new G4DNAIonisation("hydrogen_G4DNAIonisation"), particle);
      ph->RegisterProcess(new G4DNAChargeIncrease("hydrogen_G4DNAChargeIncrease"), particle);
 
    } else if ( particleName == "alpha" ) {
      ph->RegisterProcess(new G4DNAElastic("alpha_G4DNAElastic"), particle);
      ph->RegisterProcess(new G4DNAExcitation("alpha_G4DNAExcitation"), particle);
      ph->RegisterProcess(new G4DNAIonisation("alpha_G4DNAIonisation"), particle);
      ph->RegisterProcess(new G4DNAChargeDecrease("alpha_G4DNAChargeDecrease"), particle);
 
    } else if ( particleName == "alpha+" ) {
      ph->RegisterProcess(new G4DNAElastic("alpha+_G4DNAElastic"), particle);
      ph->RegisterProcess(new G4DNAExcitation("alpha+_G4DNAExcitation"), particle);
      ph->RegisterProcess(new G4DNAIonisation("alpha+_G4DNAIonisation"), particle);
      ph->RegisterProcess(new G4DNAChargeDecrease("alpha+_G4DNAChargeDecrease"), particle);
      ph->RegisterProcess(new G4DNAChargeIncrease("alpha+_G4DNAChargeIncrease"), particle);
 
    } else if ( particleName == "helium" ) {
      ph->RegisterProcess(new G4DNAElastic("helium_G4DNAElastic"), particle);
      ph->RegisterProcess(new G4DNAExcitation("helium_G4DNAExcitation"), particle);
      ph->RegisterProcess(new G4DNAIonisation("helium_G4DNAIonisation"), particle);
      ph->RegisterProcess(new G4DNAChargeIncrease("helium_G4DNAChargeIncrease"), particle);
    }
    // Extension to HZE proposed by Z. Francis
    else if ( particleName == "GenericIon" ) {
      ph->RegisterProcess(new G4DNAIonisation("GenericIon_G4DNAIonisation"), particle);
    }
    
    // Warning : the following particles and processes are needed by EM Physics builders
    // They are taken from the default Livermore Physics list
    // These particles are currently not handled by Geant4-DNA
    
      // e+
      
    else if (particleName == "e+") {
 
      // Identical to G4EmStandardPhysics_option3
      
      G4eMultipleScattering* msc = new G4eMultipleScattering();
      msc->SetStepLimitType(fUseDistanceToBoundary);
      G4eIonisation* eIoni = new G4eIonisation();
      eIoni->SetStepFunction(0.2, 100*um);      
 
      ph->RegisterProcess(msc, particle);
      ph->RegisterProcess(eIoni, particle);
      ph->RegisterProcess(new G4eBremsstrahlung(), particle);
      ph->RegisterProcess(new G4eplusAnnihilation(), particle);
 
    } else if (particleName == "gamma") {
    
      // photoelectric effect - Livermore model only
      G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();
      thePhotoElectricEffect->SetEmModel(new G4LivermorePhotoElectricModel());
      ph->RegisterProcess(thePhotoElectricEffect, particle);
 
      // Compton scattering - Livermore model only
      G4ComptonScattering* theComptonScattering = new G4ComptonScattering();
      theComptonScattering->SetEmModel(new G4LivermoreComptonModel());
      ph->RegisterProcess(theComptonScattering, particle);
 
      // gamma conversion - Livermore model below 80 GeV
      G4GammaConversion* theGammaConversion = new G4GammaConversion();
      theGammaConversion->SetEmModel(new G4LivermoreGammaConversionModel());
      ph->RegisterProcess(theGammaConversion, particle);
 
      // default Rayleigh scattering is Livermore
      G4RayleighScattering* theRayleigh = new G4RayleighScattering();
      ph->RegisterProcess(theRayleigh, particle);
    }
    
   // Warning : end of particles and processes are needed by EM Physics builders 
    
  }
 
  // Deexcitation
  //
  G4VAtomDeexcitation* de = new G4UAtomicDeexcitation();
  G4LossTableManager::Instance()->SetAtomDeexcitation(de);
}

---

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

• G4EmDNAPhysics_option7.hh
• G4EmDNAPhysics_option7.cc