PhysicsList Class Reference

#include <PhysicsList.hh>

Inheritance diagram for PhysicsList:

[legend]

Collaboration diagram for PhysicsList:

[legend]

List of all members.

Public Member Functions
	PhysicsList ()
	~PhysicsList ()
Protected Member Functions
void	ConstructParticle ()
void	ConstructProcess ()
void	SetCuts ()
void	ConstructEM ()
void	ConstructHad ()
void	ConstructDecay ()
void	AddStepMax ()

---

Detailed Description

Definition at line 16 of file PhysicsList.hh.

---

Constructor & Destructor Documentation

PhysicsList::PhysicsList

(

)

Definition at line 16 of file PhysicsList.cc.

                        :  G4VUserPhysicsList()
{
  defaultCutValue = 1.0*mm;
  SetVerboseLevel(1);
}

PhysicsList::~PhysicsList

(

)

Definition at line 24 of file PhysicsList.cc.

{}

---

Member Function Documentation

void PhysicsList::ConstructParticle

(

)

[protected]

Definition at line 36 of file PhysicsList.cc.

{
  // pseudo-particles
  G4Geantino::GeantinoDefinition();

  // gamma
  G4Gamma::GammaDefinition();
  
  // e+/-
  G4Electron::ElectronDefinition();
  G4Positron::PositronDefinition();
  
  // muon
  G4MuonPlus::MuonPlusDefinition();
  G4MuonMinus::MuonMinusDefinition();
  // neutrino
  G4NeutrinoE::NeutrinoEDefinition();
  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
  G4NeutrinoMu::NeutrinoMuDefinition();
  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
  
  G4Proton::ProtonDefinition();
  G4Neutron::NeutronDefinition();

  G4LeptonConstructor pLeptonConstructor;
  pLeptonConstructor.ConstructParticle();
  
  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();

  G4ShortLivedConstructor pShortLivedConstructor;
  pShortLivedConstructor.ConstructParticle();  

  G4BosonConstructor  pBosonConstructor;
  pBosonConstructor.ConstructParticle();

  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  

}

void PhysicsList::ConstructProcess

(

)

[protected]

Definition at line 82 of file PhysicsList.cc.

{
  AddTransportation();
  ConstructEM();
  ConstructDecay();
  ConstructHad();

  // step limitation (as a full process)
  //  
  AddStepMax();      
}

void PhysicsList::SetCuts

(

)

[protected]

Definition at line 319 of file PhysicsList.cc.

{
  //G4VUserPhysicsList::SetCutsWithDefault method sets 
  //the default cut value for all particle types 
  //
  SetCutsWithDefault();
     
  if (verboseLevel>0) DumpCutValuesTable();
}

void PhysicsList::ConstructEM

(

)

[protected]

Definition at line 114 of file PhysicsList.cc.

Referenced by ConstructProcess().

{
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();
     
    if (particleName == "gamma") {
      // gamma         
      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
      pmanager->AddDiscreteProcess(new G4ComptonScattering);
      pmanager->AddDiscreteProcess(new G4GammaConversion);
      
    } else if (particleName == "e-") {
      //electron
      pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
      pmanager->AddProcess(new G4eIonisation,        -1, 2,2);
      pmanager->AddProcess(new G4eBremsstrahlung,    -1, 3,3);      

    } else if (particleName == "e+") {
      //positron
      pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
      pmanager->AddProcess(new G4eIonisation,        -1, 2,2);
      pmanager->AddProcess(new G4eBremsstrahlung,    -1, 3,3);
      pmanager->AddProcess(new G4eplusAnnihilation,   0,-1,4);

    } else if( particleName == "mu+" || 
               particleName == "mu-"    ) {
      //muon  
      pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
      pmanager->AddProcess(new G4MuIonisation,       -1, 2,2);
      pmanager->AddProcess(new G4MuBremsstrahlung,   -1, 3,3);
      pmanager->AddProcess(new G4MuPairProduction,   -1, 4,4);
                  
    } else if (particle->GetPDGCharge() != 0.0 && 
               ! particle->IsShortLived()      &&
               particle->GetParticleName() != "chargedgeantino" ) {
      pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
      pmanager->AddProcess(new G4hIonisation,        -1, 2,2);
    }
  }
}

void PhysicsList::ConstructHad

(

)

[protected]

Definition at line 198 of file PhysicsList.cc.

Referenced by ConstructProcess().

{
// add hadronic processes for Proton
//
  G4ParticleDefinition* proton = G4Proton::Proton();
  G4ProcessManager* proton_manager = proton->GetProcessManager();

  // process: elastic scattering
  //
  G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess();
  
  // model:  use LElastic, i.e. parameterized model
  //
  G4LElastic* theElasticModel = new G4LElastic();
  theElasticProcess->RegisterMe(theElasticModel);
  proton_manager->AddDiscreteProcess(theElasticProcess);

  // process: inelastic interaction
  //   
  G4ProtonInelasticProcess* protonInelasticprocess 
                                    =  new G4ProtonInelasticProcess();
                                    
  // model:  use Binary cascade as model (good up to several GeV)
  //
  G4BinaryCascade* theModel = new G4BinaryCascade();  
  protonInelasticprocess->RegisterMe(theModel);
  
  // cross Section
  //
  protonInelasticprocess->AddDataSet(new G4NucleonNuclearCrossSection());
  
  // add process
  //
  proton_manager->AddDiscreteProcess(protonInelasticprocess);

// add hadronic processes for Neutron 
//  
  G4ParticleDefinition* neutron = G4Neutron::Neutron();
  G4ProcessManager* neutron_manager = neutron->GetProcessManager();

  // process:  elastic scattering
  //
  neutron_manager->AddDiscreteProcess(theElasticProcess);

  // process: inelastic interaction
  //     
  G4NeutronInelasticProcess* neutronInelasticprocess 
                                  =  new G4NeutronInelasticProcess();
  
  // model
  //                              
  neutronInelasticprocess->RegisterMe(theModel);
  
  // cross Section
  //
  neutronInelasticprocess->AddDataSet(new G4NucleonNuclearCrossSection());
  
  // add process
  //  
  neutron_manager->AddDiscreteProcess(neutronInelasticprocess);

// process : Capture
//  
  G4HadronCaptureProcess * theNeutronCapture = new G4HadronCaptureProcess();
  G4LCapture * theNeutronCaptureModel = new G4LCapture();
  theNeutronCapture->RegisterMe(theNeutronCaptureModel);
  neutron_manager->AddDiscreteProcess(theNeutronCapture);

// process : Fission
//    
  G4HadronFissionProcess * theNeutronFission = new G4HadronFissionProcess();
  G4LFission * theNeutronFissionModel = new G4LFission();
  theNeutronFission->RegisterMe(theNeutronFissionModel);
  neutron_manager->AddDiscreteProcess(theNeutronFission);

//  code for higher energies
/* 

   G4TheoFSGenerator * theHEModel = new G4TheoFSGenerator();
   G4FTFModel * theStringModel = new G4FTFModel();
   G4ExcitedStringDecay * theStringFrag = 
                    new G4ExcitedStringDecay(new G4LundStringFragmentation());
   theStringModel->SetFragmentationModel(theStringFrag);
        
   theHEModel->SetTransport(new G4BinaryCascade());
   theHEModel->SetMinEnergy(4.*GeV);
   theHEModel->SetMaxEnergy(100*TeV);
   theHEModel->SetHighEnergyGenerator(theStringModel);
   
   theModel->SetMaxEnergy(5*GeV);  // reduce use of casacde to below 5 GeV

   protonInelasticprocess->RegisterMe(theHEModel);
  
   neutronInelasticprocess->RegisterMe(theHEModel);
   
*/    
}

void PhysicsList::ConstructDecay

(

)

[protected]

Definition at line 162 of file PhysicsList.cc.

Referenced by ConstructProcess().

{
  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    if (theDecayProcess->IsApplicable(*particle)) {
      pmanager ->AddProcess(theDecayProcess);
      // set ordering for PostStepDoIt and AtRestDoIt
      pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
      pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
    }
  }  
}

void PhysicsList::AddStepMax

(

)

[protected]

Definition at line 301 of file PhysicsList.cc.

Referenced by ConstructProcess().

{
  // Step limitation seen as a process
  G4StepLimiter* stepMaxProcess = new G4StepLimiter();

  theParticleIterator->reset();
  while ((*theParticleIterator)()){
      G4ParticleDefinition* particle = theParticleIterator->value();
      G4ProcessManager* pmanager = particle->GetProcessManager();
      if (particle->GetPDGCharge() != 0.0)
        {
          pmanager ->AddDiscreteProcess(stepMaxProcess);
        }
  }
}

---

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

• include/PhysicsList.hh
• src/PhysicsList.cc