Geant4  10.02
G4DNAExcitation.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 // $Id: G4DNAExcitation.cc 91992 2015-08-13 07:20:24Z gcosmo $
27 
28 #include "G4DNAExcitation.hh"
30 #include "G4SystemOfUnits.hh"
31 #include "G4Positron.hh"
32 
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 
35 using namespace std;
36 
38  G4ProcessType type) :
39  G4VEmProcess(processName, type), isInitialised(false)
40 {
42 }
43 
44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45 
47 {
48 }
49 
50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
51 
53 {
54 
57 
58  return (&p == G4Electron::Electron() || &p == G4Positron::Positron()
60  || &p == instance->GetIon("hydrogen")
61  || &p == instance->GetIon("alpha++")
62  || &p == instance->GetIon("alpha+")
63  || &p == instance->GetIon("helium"));
64 }
65 
66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
67 
69 {
70  if(!isInitialised)
71  {
72  isInitialised = true;
73  SetBuildTableFlag(false);
74 
76 
77  if(name == "e-")
78  {
79 
80  // Emfietzoglou model
81  /*
82  if(!EmModel()) SetEmModel(new G4DNAEmfietzoglouExcitationModel);
83  EmModel()->SetLowEnergyLimit(8.23*eV);
84  EmModel()->SetHighEnergyLimit(10*MeV);
85  */
86  // Born model
87  if(!EmModel())
88  {
90  SetEmModel(born);
91  born->SetLowEnergyLimit(9 * eV);
92  born->SetHighEnergyLimit(1 * MeV);
93  }
94  AddEmModel(1, EmModel());
95  }
96 
97  else if(name == "e+")
98  {
99  if(!EmModel())
100  {
102  SetEmModel(lepts);
103  lepts->SetLowEnergyLimit(1 * eV);
104  lepts->SetHighEnergyLimit(1 * MeV);
105  }
106  AddEmModel(1, EmModel());
107  }
108 
109  else if(name == "proton")
110  {
111  if(!EmModel(1)) // MK: Is this a correct test ?
112  {
115  SetEmModel(miller, 1);
116  miller->SetLowEnergyLimit(10 * eV);
117  miller->SetHighEnergyLimit(500 * keV);
118 
120  SetEmModel(born, 2);
121  born->SetLowEnergyLimit(500 * keV);
122  born->SetHighEnergyLimit(100 * MeV);
123  }
124 
125  AddEmModel(1, EmModel(1));
126  if(EmModel(2)) AddEmModel(2, EmModel(2));
127  }
128 
129  else if(name == "hydrogen")
130  {
131  if(!EmModel())
132  {
135  SetEmModel(miller);
136  miller->SetLowEnergyLimit(10 * eV);
137  miller->SetHighEnergyLimit(500 * keV);
138  }
139  AddEmModel(1, EmModel());
140  }
141 
142  else if(name == "alpha" || name == "alpha+" || name == "helium")
143  {
144  if(!EmModel())
145  {
148  SetEmModel(miller);
149  miller->SetLowEnergyLimit(1 * keV);
150  miller->SetHighEnergyLimit(400 * MeV);
151  }
152  AddEmModel(1, EmModel());
153  }
154  }
155 }
156 
157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
158 
160 {
161  if(EmModel(2))
162  {
163  G4cout << " Total cross sections computed from " << EmModel(1)->GetName()
164  << " and " << EmModel(2)->GetName() << " models" << G4endl;
165  }
166  else
167  {
168  G4cout << " Total cross sections computed from "
169  << EmModel()->GetName()
170  << G4endl;
171  }
172 }
173 
174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static const double MeV
Definition: G4SIunits.hh:211
void SetBuildTableFlag(G4bool val)
virtual void InitialiseProcess(const G4ParticleDefinition *)
G4DNAExcitation(const G4String &processName="DNAExcitation", G4ProcessType type=fElectromagnetic)
G4VEmModel * EmModel(G4int index=1) const
G4String name
Definition: TRTMaterials.hh:40
virtual void PrintInfo()
static G4Proton * ProtonDefinition()
Definition: G4Proton.cc:88
const G4String & GetParticleName() const
G4DNABornExcitationModel1 G4DNABornExcitationModel
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:725
virtual ~G4DNAExcitation()
void SetEmModel(G4VEmModel *, G4int index=1)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
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
static const double keV
Definition: G4SIunits.hh:213
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:732
G4ParticleDefinition * GetIon(const G4String &name)
G4ProcessType
virtual G4bool IsApplicable(const G4ParticleDefinition &)