ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
PhysListEmStandardSSM.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PhysListEmStandardSSM.cc
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 //
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 
33 #include "PhysListEmStandardSSM.hh"
34 
35 #include "G4ParticleDefinition.hh"
36 #include "G4ProcessManager.hh"
37 
38 #include "G4ComptonScattering.hh"
39 #include "G4GammaConversion.hh"
40 #include "G4PhotoElectricEffect.hh"
41 
42 #include "G4CoulombScattering.hh"
46 
47 #include "G4eIonisation.hh"
48 #include "G4eBremsstrahlung.hh"
49 #include "G4eplusAnnihilation.hh"
50 
51 #include "G4MuIonisation.hh"
52 #include "G4MuBremsstrahlung.hh"
53 #include "G4MuPairProduction.hh"
54 
55 #include "G4hIonisation.hh"
56 #include "G4ionIonisation.hh"
57 
58 #include "G4SystemOfUnits.hh"
59 
60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61 
63  : G4VPhysicsConstructor(name)
64 {}
65 
66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
67 
69 {}
70 
71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72 
74 {
75  // Add standard EM Processes
76 
78  particleIterator->reset();
79  while( (*particleIterator)() ){
81  G4ProcessManager* pmanager = particle->GetProcessManager();
82  G4String particleName = particle->GetParticleName();
83 
84  if (particleName == "gamma") {
85  // gamma
89 
90  } else if (particleName == "e-") {
91  //electron
92  pmanager->AddProcess(new G4eIonisation, -1, 1, 1);
93  pmanager->AddProcess(new G4eBremsstrahlung, -1, 2, 2);
94 
98  //model->SetLowEnergyThreshold(10*eV);
99  model->SetPolarAngleLimit(0.0);
100  cs->AddEmModel(0, model);
101  pmanager->AddDiscreteProcess(cs);
102 
103  } else if (particleName == "e+") {
104  //positron
105 
106  pmanager->AddProcess(new G4eIonisation, -1, 1, 1);
107  pmanager->AddProcess(new G4eBremsstrahlung, -1, 2, 2);
108  pmanager->AddProcess(new G4eplusAnnihilation, 0,-1, 3);
109 
113  model->SetPolarAngleLimit(0.0);
114  cs->AddEmModel(0, model);
115  pmanager->AddDiscreteProcess(cs);
116 
117  } else if (particleName == "mu+" ||
118  particleName == "mu-" ) {
119  //muon
120  pmanager->AddProcess(new G4MuIonisation, -1, 1, 1);
121  pmanager->AddProcess(new G4MuBremsstrahlung, -1, 2, 2);
122  pmanager->AddProcess(new G4MuPairProduction, -1, 3, 3);
125  model->SetPolarAngleLimit(0.0);
126  cs->AddEmModel(0, model);
127  pmanager->AddDiscreteProcess(cs);
128 
129 
130  } else if (particleName == "alpha" || particleName == "He3") {
131  pmanager->AddProcess(new G4ionIonisation, -1, 1, 1);
134  cs->SetBuildTableFlag(false);
135  pmanager->AddDiscreteProcess(cs);
136 
137  } else if (particleName == "GenericIon" ) {
138  pmanager->AddProcess(new G4ionIonisation, -1, 1, 1);
141  cs->SetBuildTableFlag(false);
142  pmanager->AddDiscreteProcess(cs);
143 
144  } else if ((!particle->IsShortLived()) &&
145  (particle->GetPDGCharge() != 0.0) &&
146  (particle->GetParticleName() != "chargedgeantino")) {
147  //all others charged particles except geantino
148  pmanager->AddProcess(new G4hIonisation, -1, 1, 1);
149  pmanager->AddDiscreteProcess(new G4CoulombScattering);
150  }
151  }
152 }
153 
154 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
155