ExtPhysicsList.cxx

Go to the documentation of this file.

//
// File: ExtPhysicsList.cc
// Date: 2005.3.9
// Author: L.L.Wang
//
 
#include "ExtPhysicsList.h"
 
#include "ExtMuIonisation.h"
#include "ExteIonisation.h"
#include "ExthIonisation.h"
// #include "ExtMultipleScattering.h"
 
#include "G4ParticleTypes.hh"
 
// #include "G4MultipleScattering.hh"
#include "G4MuBremsstrahlung.hh"
#include "G4MuIonisation.hh"
#include "G4eBremsstrahlung.hh"
#include "G4eIonisation.hh"
#include "G4hIonisation.hh"
 
#include "G4ProcessManager.hh"
 
ExtPhysicsList::ExtPhysicsList() { defaultCutValue = 0.5 * mm; }
 
ExtPhysicsList::~ExtPhysicsList() { ; }
 
void ExtPhysicsList::ConstructParticle() {
  // In this method, static member functions should be called
  // for all particles which you want to use.
  // This ensures that objects of these particle types will be
  // created in the program.
  G4Gamma::GammaDefinition();
  G4Electron::ElectronDefinition();
  G4Positron::PositronDefinition();
  G4MuonPlus::MuonPlusDefinition();
  G4MuonMinus::MuonMinusDefinition();
  G4PionPlus::PionPlusDefinition();
  G4PionMinus::PionMinusDefinition();
  G4KaonPlus::KaonPlusDefinition();
  G4KaonMinus::KaonMinusDefinition();
  G4Proton::ProtonDefinition();
  G4AntiProton::AntiProtonDefinition();
}
 
void ExtPhysicsList::ConstructProcess() {
  // Define transportation process (must be invoked)
  AddTransportation();
 
  auto theParticleIterator = GetParticleIterator();
  // Add Energy Loss Process
  theParticleIterator->reset();
  while ( ( *theParticleIterator )() )
  {
    G4ParticleDefinition* particle     = theParticleIterator->value();
    G4ProcessManager*     pmanager     = particle->GetProcessManager();
    G4String              particleName = particle->GetParticleName();
    if ( particleName == "e-" || particleName == "e+" )
    {
      //         pmanager->AddProcess( new G4MultipleScattering,-1,1,1);
      //         pmanager->AddProcess( new G4eIonisation,-1,2,2);
      pmanager->AddProcess( new ExteIonisation, -1, 1, 1 );
      //         pmanager->AddProcess( new G4eBremsstrahlung,-1,2,2);
    }
    else if ( particleName == "mu-" || particleName == "mu+" )
    {
      //         pmanager->AddProcess( new G4MultipleScattering,-1,1,1);
      //         pmanager->AddProcess( new G4MuIonisation,-1,2,2);
      pmanager->AddProcess( new ExtMuIonisation, -1, 1, 1 );
      //         pmanager->AddProcess( new G4MuBremsstrahlung,-1,2,2);
    }
    else if ( ( !particle->IsShortLived() ) && ( particle->GetPDGCharge() != 0.0 ) &&
              ( particleName != "chargedgeantino" ) )
    {
      //         pmanager->AddProcess( new G4MultipleScattering,-1,1,1);
      //         pmanager->AddProcess( new G4hIonisation,-1,2,2);
      //         pmanager->AddProcess( new ExtMultipleScattering,-1,1,1);
      pmanager->AddProcess( new ExthIonisation, -1, 1, 1 );
    }
  }
}
 
void ExtPhysicsList::SetCuts() {
  // " G4VUserPhysicsList::SetCutsWithDefault" method sets
  //   the default cut value for all particle types
  SetCutsWithDefault();
}