Geant4  10.02.p01
G4DNAElastic.cc
Go to the documentation of this file.
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 // add elastic scattering processes of proton, hydrogen, helium, alpha+, alpha++
27 
28 #include "G4DNAElastic.hh"
29 #include "G4SystemOfUnits.hh"
30 #include "G4Positron.hh"
31 
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 
34 using namespace std;
35 
37  G4VEmProcess(processName, type), isInitialised(false)
38 {
40 }
41 
42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 
45 {
46 }
47 
48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
49 
51 {
54 
55  return (&p == G4Electron::Electron() || &p == G4Positron::Positron()
56  || &p == G4Proton::Proton() || &p == instance->GetIon("hydrogen")
57  || &p == instance->GetIon("alpha++")
58  || &p == instance->GetIon("alpha+")
59  || &p == instance->GetIon("helium"));
60 
61 }
62 
63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
64 
66 {
67  if(!isInitialised)
68  {
69  isInitialised = true;
70  SetBuildTableFlag(false);
71 
73 
74  if(name == "e-")
75  {
76  if(!EmModel())
77  {
81  }
82  AddEmModel(1, EmModel());
83  }
84 
85  else if(name == "proton" || name == "hydrogen")
86  {
87  if(!EmModel())
88  {
92  }
93  AddEmModel(1, EmModel());
94  }
95 
96  // "alpha" must be explicitely used, not alpha++
97  else if(name == "helium" || name == "alpha" || name == "alpha+")
98  {
99  if(!EmModel())
100  {
102  EmModel()->SetLowEnergyLimit(0 * eV);
104  }
105  AddEmModel(1, EmModel());
106  }
107  }
108 }
109 
110 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
111 
113 {
114  G4cout << " Total cross sections computed from " << EmModel()->GetName()
115  << " model" << G4endl;
116 }
117 
118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static const double MeV
Definition: G4SIunits.hh:211
void SetBuildTableFlag(G4bool val)
G4VEmModel * EmModel(G4int index=1) const
G4String name
Definition: TRTMaterials.hh:40
virtual G4bool IsApplicable(const G4ParticleDefinition &)
Definition: G4DNAElastic.cc:50
G4DNAElastic(const G4String &processName="DNAElastic", G4ProcessType type=fElectromagnetic)
Definition: G4DNAElastic.cc:36
G4bool isInitialised
Definition: G4DNAElastic.hh:64
virtual void InitialiseProcess(const G4ParticleDefinition *)
Definition: G4DNAElastic.cc:65
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:725
virtual ~G4DNAElastic()
Definition: G4DNAElastic.cc:44
void SetEmModel(G4VEmModel *, G4int index=1)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
virtual void PrintInfo()
static G4Proton * Proton()
Definition: G4Proton.cc:93
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
static G4DNAGenericIonsManager * Instance(void)
static const double eV
Definition: G4SIunits.hh:212
static G4Positron * Positron()
Definition: G4Positron.cc:94
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=0)
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
const G4String & GetName() const
Definition: G4VEmModel.hh:795
static MCTruthManager * instance
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:732
G4ParticleDefinition * GetIon(const G4String &name)
G4ProcessType