Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/src/PhysListEmStandard.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 //
 27 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 28 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 29 
 30 #include "PhysListEmStandard.hh"
 31 
 32 #include "G4BuilderType.hh"
 33 #include "G4ComptonScattering.hh"
 34 #include "G4GammaConversion.hh"
 35 #include "G4IonParametrisedLossModel.hh"
 36 #include "G4KleinNishinaModel.hh"
 37 #include "G4LivermorePhotoElectricModel.hh"
 38 #include "G4LossTableManager.hh"
 39 #include "G4MscStepLimitType.hh"
 40 #include "G4MuBremsstrahlung.hh"
 41 #include "G4MuIonisation.hh"
 42 #include "G4MuMultipleScattering.hh"
 43 #include "G4MuPairProduction.hh"
 44 #include "G4NuclearStopping.hh"
 45 #include "G4ParticleDefinition.hh"
 46 #include "G4PhotoElectricEffect.hh"
 47 #include "G4PhysicsListHelper.hh"
 48 #include "G4ProcessManager.hh"
 49 #include "G4RayleighScattering.hh"
 50 #include "G4SystemOfUnits.hh"
 51 #include "G4UAtomicDeexcitation.hh"
 52 #include "G4eBremsstrahlung.hh"
 53 #include "G4eIonisation.hh"
 54 #include "G4eMultipleScattering.hh"
 55 #include "G4eplusAnnihilation.hh"
 56 #include "G4hBremsstrahlung.hh"
 57 #include "G4hIonisation.hh"
 58 #include "G4hMultipleScattering.hh"
 59 #include "G4hPairProduction.hh"
 60 #include "G4ionIonisation.hh"
 61 
 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 63 
 64 PhysListEmStandard::PhysListEmStandard(const G4String& name) : G4VPhysicsConstructor(name)
 65 {
 66   G4EmParameters* param = G4EmParameters::Instance();
 67   param->SetDefaults();
 68   param->SetVerbose(0);
 69   param->SetMinEnergy(10 * eV);
 70   param->SetMaxEnergy(10 * TeV);
 71   param->SetNumberOfBinsPerDecade(10);
 72   param->SetStepFunction(0.1, 100 * um);
 73   param->SetStepFunctionMuHad(0.1, 50 * um);
 74   param->SetStepFunctionLightIons(0.1, 20 * um);
 75   param->SetStepFunctionIons(0.1, 1 * um);
 76   param->SetMscStepLimitType(fUseDistanceToBoundary);
 77   param->SetDeexcitationIgnoreCut(true);
 78   SetPhysicsType(bElectromagnetic);
 79 }
 80 
 81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 82 
 83 // PhysListEmStandard::~PhysListEmStandard()
 84 //{}
 85 
 86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 87 
 88 void PhysListEmStandard::ConstructProcess()
 89 {
 90   G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 91 
 92   // Add standard EM Processes
 93   //
 94   auto particleIterator = GetParticleIterator();
 95   particleIterator->reset();
 96   while ((*particleIterator)()) {
 97     G4ParticleDefinition* particle = particleIterator->value();
 98     G4String particleName = particle->GetParticleName();
 99 
100     if (particleName == "gamma") {
101       ph->RegisterProcess(new G4RayleighScattering, particle);
102       ph->RegisterProcess(new G4PhotoElectricEffect, particle);
103       G4ComptonScattering* cs = new G4ComptonScattering;
104       cs->SetEmModel(new G4KleinNishinaModel());
105       ph->RegisterProcess(cs, particle);
106       ph->RegisterProcess(new G4GammaConversion, particle);
107     }
108     else if (particleName == "e-") {
109       ph->RegisterProcess(new G4eMultipleScattering(), particle);
110       ph->RegisterProcess(new G4eIonisation(), particle);
111       ph->RegisterProcess(new G4eBremsstrahlung(), particle);
112     }
113     else if (particleName == "e+") {
114       ph->RegisterProcess(new G4eMultipleScattering(), particle);
115       ph->RegisterProcess(new G4eIonisation(), particle);
116       ph->RegisterProcess(new G4eBremsstrahlung(), particle);
117       ph->RegisterProcess(new G4eplusAnnihilation(), particle);
118     }
119     else if (particleName == "mu+" || particleName == "mu-") {
120       ph->RegisterProcess(new G4MuMultipleScattering(), particle);
121       ph->RegisterProcess(new G4MuIonisation(), particle);
122       ph->RegisterProcess(new G4MuBremsstrahlung(), particle);
123       ph->RegisterProcess(new G4MuPairProduction(), particle);
124     }
125     else if (particleName == "proton" || particleName == "pi-" || particleName == "pi+") {
126       ph->RegisterProcess(new G4hMultipleScattering(), particle);
127       ph->RegisterProcess(new G4hIonisation(), particle);
128       ph->RegisterProcess(new G4hBremsstrahlung(), particle);
129       ph->RegisterProcess(new G4hPairProduction(), particle);
130     }
131     else if (particleName == "alpha" || particleName == "He3") {
132       ph->RegisterProcess(new G4hMultipleScattering(), particle);
133       ph->RegisterProcess(new G4ionIonisation(), particle);
134       ph->RegisterProcess(new G4NuclearStopping(), particle);
135     }
136     else if (particleName == "GenericIon") {
137       ph->RegisterProcess(new G4hMultipleScattering(), particle);
138       G4ionIonisation* ionIoni = new G4ionIonisation();
139       ionIoni->SetEmModel(new G4IonParametrisedLossModel());
140       ph->RegisterProcess(ionIoni, particle);
141       ph->RegisterProcess(new G4NuclearStopping(), particle);
142     }
143     else if ((!particle->IsShortLived()) && (particle->GetPDGCharge() != 0.0)
144              && (particle->GetParticleName() != "chargedgeantino"))
145     {
146       // all others charged particles except geantino
147       ph->RegisterProcess(new G4hMultipleScattering(), particle);
148       ph->RegisterProcess(new G4hIonisation(), particle);
149     }
150   }
151 
152   // Deexcitation
153   //
154   G4VAtomDeexcitation* de = new G4UAtomicDeexcitation();
155   G4LossTableManager::Instance()->SetAtomDeexcitation(de);
156 }
157 
158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
159